Publications

2026

• Early-Reverberation Imaging Functions for Bounded Elastic Domains
  
  • Ducasse Eric
  • Rodriguez Samuel
  • Bonnet Marc
  Acta Acustica, EDP Sciences, 2026, 10, pp.2. For the ultrasonic inspection of bounded elastic structures, finite-duration imaging functions are derived in the Fourier-Laplace domain.The signals involved are exponentially windowed, so that early reflections are taken into account more strongly than later ones in the imaging methodology.Applying classical approaches to the general case of anisotropic elasticity, we express the Fréchet derivatives of the relevant data-misfit functional with respect to arbitrary perturbations of the mass density and stiffnesses in terms of forward and adjoint solutions.Their definitions incorporate the exponentially decaying weighting. The proposed finite-duration imaging functions are then defined on that basis.As some areas of the structure are less insonified than others, it is necessary to define normalized imaging functions to compensate for these variations.Our approach in particular aims to overcome the difficulty of dealing with bounded domains containing defects not located in direct line of sight from the transducers and measured signals of long duration.For this initiation work, we demonstate the potential of the proposed method on a two-dimensional test case featuring the imaging of mass and elastic stiffness variations in a region of a bounded isotropic medium that is not directly visible from the transducers. (10.1051/aacus/2025069)
  DOI : 10.1051/aacus/2025069
• Hybrid FEM/IPDG semi-implicit schemes for time domain electromagnetic wave propagation in non cylindrical coaxial cables
  
  • Beni Hamad Akram
  • Imperiale Sébastien
  • Joly Patrick
  ESAIM: Mathematical Modelling and Numerical Analysis, Société de Mathématiques Appliquées et Industrielles (SMAI) / EDP, 2026. In this work, we develop an efficient numerical method for solving 3D Maxwell's equations in non-cylindrical coaxial cables. The main challenge arises from the elongated geometry of the computational domain, which induces strong anisotropy between the longitudinal direction (along the cable) and the transverse directions (within the cross-sections). This leads to the use of highly anisotropic meshes, where the longitudinal mesh size is much larger than the transverse one.<p>Our objective is to design a numerical scheme that is explicit in the longitudinal direction, with a CFL stability condition depending only on the longitudinal mesh size. In a previous work, we achieved this for cylindrical cables by employing prismatic edge elements, 1D quadrature for longitudinal mass lumping, and a hybrid explicit/implicit time discretization. The present paper extends this approach to non-cylindrical cables, addressing several new difficulties with the following key ingredients: (1) representing the cable as a deformation of a reference cylindrical cable and employing mapping techniques between the physical and reference domains; (2) using an anisotropic space discretization that combines an interior penalty discontinuous Galerkin (IPDG) method in the transverse directions with a conforming finite element method in the longitudinal direction; (3) utilizing prismatic edge elements on a prismatic mesh of the reference cable; and (4) adapting the construction of the hybrid explicit-implicit time discretization to the new structure of the semidiscrete problem. From a theoretical perspective, the main difficulty lies in the stability analysis, which requires extending and adapting standard techniques for DG methods in space and energy methods in time.</p>
• Multiscale modeling for a class of high-contrast heterogeneous sign-changing problems
  
  • Chung Eric T.
  • Ciarlet Patrick
  • Jin Xingguang
  • Ye Changqing
  Journal of Computational Mathematics -International Edition-, Global Science Press, 2026. The mathematical formulation of sign-changing problems involves a linear second-order partial differential equation in the divergence form, where the coefficient can assume positive and negative values in different subdomains. These problems find their physical background in negative-index metamaterials, either as inclusions embedded into common materials as the matrix or vice versa. In this paper, we propose a numerical method based on the constraint energy minimizing generalized multiscale finite element method (CEMGMsFEM) specifically designed for sign-changing problems. The construction of auxiliary spaces in the original CEM-GMsFEM is tailored to accommodate the sign-changing setting. The numerical results demonstrate the effectiveness of the proposed method in handling sophisticated coefficient profiles and the robustness of coefficient contrast ratios. Under several technical assumptions and by applying the T-coercivity theory, we establish the inf-sup stability and provide an a priori error estimate for the proposed method.
• A general-purpose global regularization method for 3D volume integral operators
  
  • Anderson Thomas G.
  • Bonnet Marc
  • Faria Luiz M.
  • Pérez-Arancibia Carlos
  , 2026. Singular volume integral operators associated with constant-coefficient partial differential operators extend the applicability of potential theory to inhomogeneous problems, for example arising from nonlinearities or variable coefficients. Typically the PDE kernels in these operators give rise to singularities at all O(1/h^3 ) volume discretization/evaluation points in a mesh of characteristic size h, while the slowly-decaying nature of such kernels give rise to long-range interactions that require coupling to fast summation algorithms. The presented method uses Green's identities to regularize a wide variety of both scalar-valued and vector-valued volume integral operators by use of a certain regularizing volume density interpolant. The analysis shows how the regularizing effect of the interpolant is global in the sense that the interpolation quality increases in an exactly compensatory fashion as the distance to the Green's function singularity decreases. High-order convergence estimates with tabulated simplex quadratures are established, including with exact representation of curved domains.
• Trapped modes in electromagnetic waveguides
  
  • Bonnet-Ben Dhia Anne-Sophie
  • Chesnel Lucas
  • Fliss Sonia
  Integral Equations and Operator Theory, Springer Verlag, 2026. We consider the Maxwell's equations with perfect electric conductor boundary conditions in three-dimensional unbounded domains which are the union of a bounded resonator and one or several semi-infinite waveguides. We are interested in the existence of electromagnetic trapped modes, i.e. $L^2$ solutions of the problem without source term. These trapped modes are associated to eigenvalues of the Maxwell's operator, that can be either below the essential spectrum or embedded in it. First for homogeneous waveguides, we present different families of geometries for which we can prove the existence of eigenvalues. Then we exhibit certain non homogeneous waveguides with local perturbations of the dielectric constants that support trapped modes. Let us mention that some of the mechanisms we propose are very specific to Maxwell's equations and have no equivalent for the scalar Dirichlet or Neumann Laplacians.
• On some coupled local and nonlocal diffusion models
  
  • Borthagaray Juan Pablo
  • Ciarlet Patrick
  Mathematical Models and Methods in Applied Sciences, World Scientific Publishing, 2026. We study problems in which a local model is coupled with a nonlocal one. We propose two energies: both of them are based on the same classical weighted $H^1$-semi norm to model the local part, while two different weighted $H^s$-semi norms, with $s \in (0, 1)$, are used to model the nonlocal part. The corresponding strong formulations are derived. In doing so, one needs to develop some technical tools, such as suitable integration by parts formulas for operators with variable diffusivity, and one also needs to study the mapping properties of the Neumann operators that arise. In contrast to problems coupling purely local models, in which one requires transmission conditions on the interface between the subdomains, the presence of a nonlocal operator may give rise to nonlocal fluxes. These nonlocal fluxes may enter the problem as a source term, thereby changing its structure. Finally, we focus on a specific problem, that we consider most relevant, and study regularity of solutions and finite element discretizations. We provide numerical experiments to illustrate the most salient features of the models. (10.1142/S0218202526500442)
  DOI : 10.1142/S0218202526500442
• Predicting topologically protected interface state with high-frequency homogenization
  
  • Touboul Marie
  • Lombard Bruno
  • Coutant Antonin
  Comptes-Rendus-de-l'Academie-des-Sciences, 2026, 354, pp.269-291. When two semi-infinite periodic media are joined together, a localized interface mode may exist, whose frequency belongs to their common band gap. Moreover, if certain spatial symmetries are satisfied, this mode is topologically protected and thus is robust to defects. A method has recently been proposed to identify the existence and the frequency of this mode, based on the computation of surface impedances at all the frequencies in the gap. In this work, we approximate the surface impedances thanks to highfrequency effective models, and therefore get a prediction of topologically protected interface states while only computing the solution of an eigenvalue problem at the edges of the bandgaps. We also show that the nearby eigenvalues high-frequency effective models give rise to a better approximation of the surface impedance.
• Crouzeix-Raviart elements on simplicial meshes in $d$ dimensions
  
  • Bohne Nis-Erik
  • Ciarlet Patrick
  • Sauter Stefan
  Foundations of Computational Mathematics, Springer Verlag, 2026. In this paper we introduce Crouzeix-Raviart elements of general polynomial order $k$ and spatial dimension $d\geq2$ for simplicial finite element meshes. We give explicit representations of the non-conforming basis functions and prove that the conforming companion space, i.e., the conforming finite element space of polynomial order $k$ is contained in the Crouzeix-Raviart space. We prove a direct sum decomposition of the Crouzeix-Raviart space into (a subspace of) the conforming companion space and the span of the non-conforming basis functions. Degrees of freedom are introduced which are bidual to the basis functions and give rise to the definition of a local approximation/interpolation operator. In two dimensions or for $k=1$, these freedoms can be split into simplex and $(d-1)$ dimensional facet integrals in such a way that, in a basis representation of Crouzeix-Raviart functions, all coefficients which belong to basis functions related to lower-dimensional faces in the mesh are determined by these facet integrals. It will also be shown that such a set of degrees of freedom does not exist in higher space dimension and $k&gt;1$.
• A hybridizable discontinuous Galerkin method with transmission variables for time-harmonic electromagnetic problems
  
  • Rappaport Ari
  • Chaumont-Frelet Théophile
  • Modave Axel
  SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2026, 48 (4), pp.B629-B650. The CHDG method is a hybridizable discontinuous Galerkin (HDG) finite element method suitable for the iterative solution of time-harmonic wave propagation problems. Hybrid unknowns corresponding to transmission variables are introduced at the element interfaces and the physical unknowns inside the elements are eliminated, resulting in a hybridized system with favorable properties for fast iterative solution. In this paper, we extend the CHDG method, initially studied for the Helmholtz equation, to the time-harmonic Maxwell equations. We prove that the local problems stemming from hybridization are well-posed and that the fixed-point iteration naturally associated to the hybridized system is contractive. We propose a 3D implementation with a discrete scheme based on nodal basis functions. The resulting solver and different iterative strategies are studied with several numerical examples using a high-performance parallel C++ code. (10.1137/25M1756818)
  DOI : 10.1137/25M1756818
• Solving numerically the two-dimensional time harmonic Maxwell problem with sign-changing coefficients
  
  • Chaaban Farah
  • Ciarlet Patrick
  • Rihani Mahran
  ESAIM: Mathematical Modelling and Numerical Analysis, Société de Mathématiques Appliquées et Industrielles (SMAI) / EDP, 2026. We are investigating the numerical solution to the 2D time-harmonic Maxwell equations in the presence of a classical medium and a metamaterial, that is with sign-changing coefficients. As soon as the problem has a (unique) solution, we are able to build a converging numerical approximation based on the finite element method, for which there is no constraint on the meshes related to the sign-changing behavior. To that aim, we use Lagrange finite elements to approximate the scalar potentials appearing in the Helmholtz decomposition of the vector-valued electromagnetic fields. Convergence in strong norm is proven for the fields. Numerical examples illustrate the theory.
• Multiscale methods for wave propagation in materials with sign-changing coefficients
  
  • Chung Eric T.
  • Ciarlet Patrick
  • Jin Xingguang
  • Ye Changqing
  Multiscale Modeling and Simulation: A SIAM Interdisciplinary Journal, Society for Industrial and Applied Mathematics, 2026. From a mathematical perspective, the extraordinary properties of metamaterials are often reflected in the coefficients of the governing partial differential equations (PDEs). These coefficients may fall outside the assumptions of classical theory, particularly when the effective dielectric permittivity and/or magnetic permeability are negative. This situation can transform a coercive operator into a non-coercive one, potentially leading to ill-posedness. In this paper, we utilize the Constraint Energy Minimizing Generalized Multiscale Finite Element Method (CEM-GMsFEM), specifically designed for time-harmonic electromagnetic wave problems, where the construction of auxiliary spaces in the original CEM-GMsFEM is tailored to accommodate the sign-changing setting. Based on the framework of T-coercivity theory and resolution conditions, we establish the inf-sup stability and provide an a priori error estimate for the proposed method. The numerical results demonstrate the effectiveness and robustness of our approach in handling such sophisticated coefficient profiles.
• Well-posed homogenized strain-gradient models for linear elastodynamics and elastostatics in arbitrary periodic media
  
  • Cornaggia Rémi
  • Bonnet Marc
  • Rosi Giuseppe
  • El Ouafa Saad
  • Auffray Nicolas
  , 2026. This work develops well-posed homogenized strain-gradient models for linear elastostatics and elastodynamics in periodic media, with a primary focus on elastic wave propagation. Using the classical two-scale asymptotic expansion method, we carry out second-order periodic homogenization for media in $\Rbb^d$ ($d = 2, 3$), with no restriction on the periodicity cell geometry or material distribution. Reciprocity identities applied to suitably chosen pairs of cell solutions provide alternative expressions for the effective stiffness and inertia tensors arising at the leading, first and second orders, substantially reducing the number of cell problems that must actually be solved. Since direct two-scale homogenization beyond leading order generically yields ill-posed effective operators, a Boussinesq-trick procedure is introduced, involving a tunable scalar weight, to recast the resulting fourth-order partial differential equation as a valid strain-gradient elasticity (SGE) model possessing the requisite symmetry, sign-definiteness and coercivity properties. These properties are then used, via the Hille-Yosida theorem, to establish the well-posedness of the corresponding transient initial-value and forced-response problems. Several practically relevant special cases are examined, including centrosymmetric cells, homogeneous mass density and homogeneous elasticity, each yielding simplified model structures. Numerical illustrations on three two-dimensional periodicity cells (square, hexagonal and a non-centrosymmetric chiral lattice) compare the resulting dispersion relations against reference Floquet-Bloch computations and assess transient wave propagation, demonstrating the model's capacity to capture anisotropic and dispersive effects beyond classical elasticity while preserving mathematical well-posedness.
• Squirmers with arbitrary shape and slip: modeling, simulation, and optimization
  
  • Das Kausik
  • Zhu Hai
  • Bonnet Marc
  • Veerapaneni Shravan
  , 2026. We consider arbitrary-shaped microswimmers of spherical topology and propose a framework for expressing their slip velocity in terms of tangential basis functions defined on the boundary of the swimmer using the Helmholtz decomposition. Given a time-independent slip velocity profile, we show that the trajectory followed by the microswimmer is a circular helix. We derive analytical expressions for the translational and rotational velocities of a prolate spheroid swimmer in terms of its Helmholtz decomposition modes and explore the effect of aspect ratio on these rigid body velocities. Then, for a given arbitrary swimmer shape of spherical topology, we investigate which slip profile minimizes the total power loss. A partial minimization is performed in which the direction of net motion of the swimmer is prescribed, followed by a global optimization procedure in which the best net motion direction is determined. The optimization results suggest that the competition between linear and rotational optimal motion is linked to symmetries in the shape of the microswimmer.
• A HDG method with transmission variables for time-harmonic wave propagation problems with constant coefficients
  
  • Pescuma Simone
  • Gabard Gwénaël
  • Chaumont-Frelet Théophile
  • Modave Axel
  , 2026. Iterative finite element solvers for time-harmonic wave problems are notoriously slow to converge, owing to fundamental properties of these problems. We present a variant of the hybridizable discontinuous Galerkin (HDG) method that is better suited to fast iterative solution. Unlike the standard hybridization strategy, which eliminates physical unknowns by introducing an auxiliary numerical flux on element faces, our approach instead introduces a transmission variable on those faces. For Helmholtz problems, this reformulation, known as CHDG, has been shown to significantly accelerate the convergence of iterative schemes relative to standard HDG. The present work extends CHDG to a general framework covering wave propagation problems with constant coefficients, capable of handling diverse wave types in a unified manner. We prove that the resulting hybridized system is well-posed and amenable to fixed-point iteration. As a practical application, we apply the method to the time-harmonic linearized Euler equations with a uniform subsonic mean flow. The method is illustrated through two-dimensional numerical benchmarks involving both sound and vorticity waves, with a systematic comparison of the convergence behaviour of several iterative schemes across a range of configurations.
• Analysis of time-harmonic electromagnetic problems with elliptic material coefficients
  
  • Ciarlet Patrick
  • Modave Axel
  Mathematical Methods in the Applied Sciences, Wiley, 2026, 49, pp.3797-3815. We consider time-harmonic electromagnetic problems with material coefficients represented by elliptic fields, covering a wide range of complex and anisotropic material media. The properties of elliptic fields are analyzed, with particular emphasis on scalar fields and normal tensor fields. Time-harmonic electromagnetic problems with general elliptic material fields are then studied. Well-posedness results for classical variational formulations with different boundary conditions are reviewed, and hypotheses for the coercivity of the corresponding sesquilinear forms are investigated. Finally, the proposed framework is applied to examples of media used in the literature: isotropic lossy media, geometric media, and gyrotropic media. (10.1002/mma.70318)
  DOI : 10.1002/mma.70318
• Comprehensive dataset of continuously monitored hydrodynamic and physico-chemical parameters in various karst systems, France
  
  • Cinkus Guillaume
  • Bondu Raphaël
  • Aliouache Mohammed
  • Jourde Hervé
  • Mazzilli Naomi
  • Arfib Bruno
  • Bailly-Comte Vincent
  • Batiot-Guilhe Christelle
  • Celle Hélène
  • Delbart Célestine
  • Fournier Matthieu
  • Labat David
  • Peyraube Nicolas
  • Steinmann Marc
  • Valdès-Lao Danièle
  • Baudement Cécile
  • Béranger Sandra
  • Binet Stéphane
  • Boetsch Anne
  • Brunet Pascal
  • Carrière Simon Damien
  • Chalikakis Konstantinos
  • Charlier Jean-Baptiste
  • Cousquer Yohann
  • Défarge Christian
  • Ducasse Joshua
  • Emblanch Christophe
  • Fabre Juliette
  • Fleury Perrine
  • Garin Thibaut
  • Huon Julien
  • Jardani Abderrahim
  • Johannet Anne
  • Jouves Johan
  • Jozja Nevila
  • Justy Lucile
  • Ladouche Bernard
  • Lamarque Thierry
  • Lastennet Roland
  • Léonardi Véronique
  • Lobry Olivier
  • Lorette Guillaume
  • Loup Christophe
  • Marchand Pierre
  • Manière Louis
  • Maréchal Jean-Christophe
  • Martin Lucie
  • Muller Rémi
  • Massei Nicolas
  • Naessens Fabien
  • Ollivier Chloé
  • Probst Anne
  • Seidel Jean-Luc
  • Serene Leïla
  • Sivelle Vianney
  • Zappelli Alexandre
  Scientific Data, Nature Publishing Group, 2026. Karst aquifers are a crucial source of water, supplying approximately 10% of the global population and often serving as the sole water resource in certain regions. These aquifers are characterized by highly heterogeneous flow dynamics and exhibit significant temporal variability in both hydrodynamic and physico-chemical conditions. Continuous monitoring of these parameters is essential for advancing our understanding of karst aquifer functioning; however, comprehensive, high-frequency datasets remain limited. We present a comprehensive dataset covering 13 karst springs monitored across nine observatories of the French Karst National Observatory Service (SNO KARST), spanning various hydroclimatic regions (oceanic, mountainous, Mediterranean). The SNO KARST aims to strengthen knowledge-sharing and to promote cross-disciplinary research on karst systems at the national scale. The dataset includes: (1) hydrodynamic data (water level, discharge), and (2) physico-chemical data (water temperature, electric conductivity, pH, dissolved oxygen, turbidity, Total Organic Carbon (TOC), Dissolved Organic Carbon (DOC), nitrate, and organic matter fluorescence). Spanning over a decade of continuous monitoring, such a dataset is required for the analysis of the hydrological and physico-chemical dynamics of karst aquifers, the assessment of their vulnerability to pollution and climate change, and the modeling of hydrodynamic and hydrochemical variables, ultimately aiming to improve the management and preservation of these critical water resources in contrasted contexts. (10.1038/s41597-026-07561-0)
  DOI : 10.1038/s41597-026-07561-0