Publications

2026

• Poisson-type problems with transmission conditions at boundaries of infinite metric trees
  
  • Kachanovska Maryna
  • Naderi Kiyan
  • Pankrashkin Konstantin
  Journal of Mathematical Analysis and Applications, Elsevier, 2026, 557 (1), pp.130261. The paper introduces a Poisson-type problem on a mixed-dimensional structure combining a Euclidean domain and a lower-dimensional self-similar component touching along a compact surface (interface). The lower-dimensional piece is a so-called infinite metric tree (one-dimensional branching structure), and the key ingredient of the study is a rigorous definition of the gluing conditions between the two components. These constructions are based on the recent concept of embedded trace maps and some abstract machineries derived from a suitable Green-type formula. The problem is then reduced to the study of Fredholm properties of a linear combination of Dirichlet-to-Neumann maps for the tree and the Euclidean domain, which yields desired existence and uniqueness results. One also shows that large finite sections of the tree can be used for an efficient approximation of solutions (10.1016/j.jmaa.2025.130261)
  DOI : 10.1016/j.jmaa.2025.130261
• Stability of time stepping methods for discontinuous Galerkin discretizations of Friedrichs' systems
  
  • Imperiale Sébastien
  • Joly Patrick
  • Rodríguez Jerónimo
  , 2025. In this work we study new various energy-based theoretical results on the stability of s-stages, s-th order explicit Runge-Kutta integrators as well as a modified leap-frog scheme applied to discontinuous Galerkin discretizations of transient linear symmetric hyperbolic Friedrichs' systems. We restrict the present study to conservative systems and Cauchy problems.
• Slip optimization on arbitrary 3D microswimmers: a reduced-dimension and boundary-integral framework
  
  • Bonnet Marc
  • Das Kausik
  • Veerapaneni Shravan
  • Zhu Hai
  , 2026. This article presents a computational framework for determining the optimal slip velocity of a microswimmer with arbitrary three-dimensional geometry suspended in a viscous fluid. The objective is to minimize the hydrodynamic power dissipation required to maintain unit speed along the net swimming direction. By exploiting the linearity of the Stokes equations and the Lorentz reciprocal theorem, we derive an explicit linear operator that maps the tangential surface slip velocity to the resulting rigid-body translational and rotational velocities, effectively decoupling the hydrodynamic boundary value problem from the optimization loop. The a priori infinite-dimensional search space for the slip optimization is reduced to the finite dimension $r$ of rigid-body motions by finding an appropriate subspace of the operator's domain. This reduces the PDE-constrained optimization to a low-dimensional programming problem that can be solved at negligible computational cost once the system matrices are assembled. The optimization algorithm requires 2$r$ auxiliary flow problems that are solved numerically using a high-order boundary integral method. We validate the accuracy of the proposed method and present optimal slip profiles and swimming trajectories for a variety of microswimmer shapes. We investigate the effect of some common geometrical symmetries of the swimmer shape on the resulting optimal motion, and in particular present a modified version of the slip optimization algorithm for axisymmetric shapes, where tangential rigid-body velocities may occur
• A posteriori error estimates for mixed finite element discretization of the multigroup Neutron Simplified Transport equations with Robin boundary condition
  
  • Ciarlet Patrick
  • Do Minh-Hieu
  • Gervais Mario
  • Madiot François
  , 2026. We analyse a posteriori error estimates for the discretization with mixed finite elements on simplicial or Cartesian meshes of the multigroup neutron simplified transport (SPN ) equations, in the case where a Robin (or Fourier type) boundary condition is imposed on the boundary. This boundary condition is of particular importance in neutronics, since it corresponds to the well-known vacuum boundary condition. We provide guaranteed and locally efficient estimators. In particular, a specific estimator is designed to handle the Robin boundary condition. We also develop the theory in the case of mixed imposed boundary conditions, of Dirichlet, Neumann or Fourier type. The approach is further extended to a Domain Decomposition Method, the so-called DD+L 2 jumps method. In this framework, the adaptive mesh refinement strategy is implemented for a discretization using Cartesian meshes on each subdomain. Numerical experiments illustrate the theory.
• Analysis of a two-level domain decomposition preconditioner for the time-harmonic Maxwell equations in anisotropic media
  
  • Bonazzoli Marcella
  • Ciarlet Patrick
  • Modave Axel
  • Rappaport Ari
  , 2026. We analyze a domain decomposition preconditioner, namely a two-level additive Schwarz method, for the time-harmonic Maxwell equations in anisotropic media. The material law is described by a tensor-valued electric permittivity ε, magnetic permeability µ and conductivity σ which are assumed to be uniformly symmetric positive definite in the physical domain. Convergence estimates for the preconditioned GMRES solver are obtained through bounds on the norm and the field-of-values (FOV) of the preconditioned operator. Our purpose is to extend the convergence analysis available for scalar and constant coefficients established in Bonazzoli et al. [5] to this tensorial setting. While the overall argument follows the additive Schwarz framework therein, the anisotropic case requires substantial new ingredients. Among these are a coefficient-weighted discrete Helmholtz decomposition, regularity estimates adapted to the anisotropic setting, and a stronger "high frequency regime" assumption. The latter allows control of unsigned terms that vanish via orthogonality in the scalar case. These tools are crucial for the main technical result: bounding the FOV away from the origin through estimates explicit in the frequency and anisotropy parameters, under suitable resolution assumptions.
• Waves within a network of slowly time-modulated interfaces: time-dependent effective properties, reciprocity and high-order dispersion
  
  • Darche Michaël
  • Assier Raphaël
  • Guenneau Sebastien
  • Lombard Bruno
  • Touboul Marie
  , 2026. We consider wave propagation through a 1D periodic network of slowly time-modulated interfaces. Each interface is modelled by time-dependent spring-mass jump conditions, where mass and rigidity interface parameters are modulated in time. Low-frequency homogenisation yields a leading-order model described by an effective time-dependent wave equation, i.e. a wave equation with effective mass density and Young's modulus which are homogeneous in space but depend on time. This means that time-dependent bulk effective properties can be created by an array where only interfaces are modulated in time. The occurrence of k-gaps in case of a periodic modulation is also analysed. Second-order homogenisation is then performed and leads to an effective model which is reciprocal but encapsulates higher-order dispersive effects. These findings and the limitations of the models are illustrated through time-domain simulations.
• Accelerating the Method of Reflections with Domain Decomposition techniques for Boundary Integral Equations in Multiple Scattering
  
  • Chaillat Stéphanie
  • Darbas Marion
  • Gander Martin J
  • Halpern Laurence
  , 2026. The Method of Reflections was historically introduced to obtain approximate solu-tions as series expansions for the motion of particles in suspension. It can however equally well be used for solving multiple scattering problems numerically. We show for Helmholtz multiple scattering problems that the Method of Reflections, whether applied in its alternating or parallel version, suffers from convergence problems when scatterers are close. We use boundary integral equations to formulate the methods, and then identify them as algebraic Schwarz methods, thereby interpreting them as boundary domain decomposition techniques. This connection allows us to introduce remedies such as overlap (which can be partial, covering only the illuminating region of the obstacles) and coarse spaces from domain decomposition into the Method of Reflections. This leads to substantially accelerated variants, and also naturally makes them suitable preconditioners for GMRES. These new approaches are particularly efficient for closeby obstacles. Moreover, numerical experiments show that the number of iterations remains robust with respect to the wavenumber.
• Fluid-structure Green's functions via BEM/BEM coupling for flow induced noise in arbitrary elastic geometries
  
  • Pacaut Louise
  • Chaillat Stéphanie
  • Mercier Jean-François
  • Serre Gilles
  , 2026. We address the challenge of efficiently simulating the noise generated by the interaction of a turbulent flow noise with complex elastic structures, a coupled fluid/structure interaction (FSI) problem. Current approaches typically separate vibro-acoustic and hydro-acoustic contributions, limiting the accuracy of hydrodynamic noise predictions. To overcome this limitation, we develop a numerical method for computing a Green's function tailored to the coupled FSI problem, enabling a monolithic prediction of the radiated noise without separating the two components. This approach not only improves the accuracy of hydrodynamic noise simulations but also significantly reduces computational costs. The Green's function is constructed using a novel integral formulation and solved numerically via a coupled fast BEM/ BEM solver.
• Metamaterials and Fluid Flows
  
  • Avallone Francesco
  • Bosia Federico
  • Chen Yi
  • Colombo Giada
  • Craster Richard
  • de Ponti Jacopo Maria
  • Fabbiane Nicolò
  • Haberman Michael
  • Hussein Mahmoud
  • Hwang Wontae
  • Iemma Umberto
  • Juhl Abigail
  • Kadic Muamer
  • Kotsonis Marios
  • Laude Vincent
  • Marquet Olivier
  • Mery Fabien
  • Michelis Theodoros
  • Nouh Mostafa
  • Ragni Daniele
  • Touboul Marie
  • Wegener Martin
  • Krushynska Anastasiia
  Nature Communications, Nature Publishing Group, 2026. (10.1038/s41467-026-70163-2)
  DOI : 10.1038/s41467-026-70163-2
• Discretization in multilayered media with high contrasts: is it all about the boundaries?
  
  • Carvalho Camille
  • Chaillat Stéphanie
  • Tsogka Chrysoula
  • Cortes Elsie A
  , 2026. Wave propagation in multilayered media with high material contrasts poses significant numerical challenges, as large variations in wavenumbers lead to strong reflections and complex transmission of the incoming wave field. To address these difficulties, we employ a boundary integral formulation thereby avoiding volumetric discretization. In this framework, the accuracy of the numerical solution depends strongly on how the material interfaces are discretized. In this work, we demonstrate that standard meshing strategies based on resolving the maximum wavenumber across the domain become computationally inefficient in multilayered configurations, where high wavenumbers are confined to localized subdomains. Through a systematic study of multilayer transmission problems, we show that no simple discretization rule based on the maximum wavenumber or material contrasts emerges. Instead, the wavenumber of the background (exterior) medium plays a dominant role in determining the optimal boundary resolution. Building on these insights, we propose an adaptive approach that achieves uniform accuracy and efficient computation across multiple layers. Numerical experiments for a range of multilayer configurations demonstrate the scalability and robustness of the proposed approach.
• Htool-DDM: A C++ library for parallel solvers and compressed linear systems.
  
  • Marchand Pierre
  • Tournier Pierre-Henri
  • Jolivet Pierre
  Journal of Open Source Software, Open Journals, 2026, 11 (118), pp.9279. (10.21105/joss.09279)
  DOI : 10.21105/joss.09279
• Automated far-field sound field estimation combining robotized acoustic measurements and the boundary elements method
  
  • Pascal Caroline
  • Marchand Pierre
  • Chapoutot Alexandre
  • Doaré Olivier
  Acta Acustica, EDP Sciences, 2026. The identification and reconstruction of acoustic fields radiated by unknown structures is usually performed using either Sound Field Estimation or Near-field Acoustic Holography techniques. The latter turns out to be especially useful when data is only available close to the source, but information throughout the whole space is needed. Yet, the lack of amendable and efficient implementations of state-of-the-art solutions, as well as the laborious and often lengthy deployment of acoustic measurements continue to be significant obstacles to the practical application of such methods. The purpose of this work is to address both problems. First, a completely automated metrology setup is proposed, in which a robotic arm is used to gather extensive and accurately positioned acoustic data without any human intervention. The impact of the robot on acoustic pressure measurements is cautiously evaluated, and proved to remain limited below 1 kHz. The Sound Field Estimation is then tackled using the Boundary Element Method, and implemented using the FreeFEM software. Numerically simulated measurements have allowed us to assess the method accuracy, which matches theoretically expected results and proves to remain robust against positioning inaccuracies, provided that the robot is carefully calibrated. The overall solution has been successfully tested using actual robotized measurements of an unknown loudspeaker, with a reconstruction error of less than 30 %. (10.1051/aacus/2026017)
  DOI : 10.1051/aacus/2026017
• Asymptotic analysis at any order of Helmholtz's problem in a corner with a thin layer: an algebraic approach
  
  • Baudet Cédric
  Asymptotic Analysis, IOS Press, 2026. We consider the Helmholtz equation in an angular sector partially covered by a homogeneous layer of small thickness, denoted ε. We propose in this work an asymptotic expansion of the solution with respect to ε at any order. This is done using matched asymptotic expansion, which consists here in introducing different asymptotic expansions of the solution in three subdomains: the vicinity of the corner, the layer and the rest of the domain. These expansions are linked through matching conditions. The presence of the corner makes these matching conditions delicate to derive because the fields have singular behaviors. Our approach is to reformulate these matching conditions purely algebraically by writing all asymptotic expansions as formal series. By using algebraic calculus we reduce the matching conditions to scalar relations linking the singular behaviors of the fields. These relations have a convolutive structure and involve some coefficients that can be computed analytically. Our asymptotic expansion is justified rigorously with error estimates. (10.1177/09217134251389983)
  DOI : 10.1177/09217134251389983
• Wave propagation in the frequency regime in one-dimensional quasiperiodic media -Limiting absorption principle
  
  • Amenoagbadji Pierre
  • Fliss Sonia
  • Joly Patrick
  , 2026. <div><p>We study the one-dimensional Helmholtz equation with (possibly perturbed) quasiperiodic coefficients. Quasiperiodic functions are the restriction of higher dimensional periodic functions along a certain (irrational) direction. In classical settings, for real-valued frequencies, this equation is generally not well-posed: existence of solutions in L 2 is not guaranteed and uniqueness in L ∞ may fail. This is a well-known difficulty of Helmholtz equations, but it has never been addressed in the quasiperiodic case. We tackle this issue by using the limiting absorption principle, which consists in adding some imaginary part (also called absorption) to the frequency in order to make the equation well-posed in L 2 , and then defining the physically relevant solution by making the absorption tend to zero. In previous work, we introduced a definition of the solution of the equation with absorption based on Dirichlet-to-Neumann (DtN) boundary conditions. This approach offers two key advantages: it facilitates the limiting process and has a direct numerical counterpart. In this work, we first explain why the DtN boundary conditions have to be replaced by Robin-to-Robin boundary conditions to make the absorption go to zero. We then prove, under technical assumptions on the frequency, that the limiting absorption principle holds and we propose a numerical method to compute the physical solution.</p></div>
• Discrete FEM-BEM coupling with the Generalized Optimized Schwarz Method
  
  • Boisneault Antonin
  • Bonazzoli Marcella
  • Claeys Xavier
  • Marchand Pierre
  , 2026. The present contribution aims at developing a non-overlapping Domain Decomposition (DD) approach to the solution of acoustic wave propagation boundary value problems based on the Helmholtz equation, on both bounded and unbounded domains. This DD solver, called Generalized Optimized Schwarz Method (GOSM), is a substructuring method, that is, the unknowns of an iteration are associated with the subdomains interfaces. We extend the analysis presented in a previous paper of one of the author to a fully discrete setting. We do not consider only a specific set of boundary conditions, but a whole class including, e.g., Dirichlet, Neumann, and Robin conditions. Our analysis will also cover interface conditions corresponding to a Finite Element Method - Boundary Element Method (FEM-BEM) coupling. In particular, we shall focus on three classical FEM-BEM couplings, namely the Costabel, Johnson-Nédélec and Bielak-MacCamy couplings. As a remarkable outcome, the present contribution yields well-posed substructured formulations of these classical FEM-BEM couplings for wavenumbers different from classical spurious resonances. We also establish an explicit relation between the dimensions of the kernels of the initial variational formulation, the local problems and the substructured formulation. That relation especially holds for any wavenumber for the substructured formulation of Costabel FEM-BEM coupling, which allows us to prove that the latter formulation is well-posed even at spurious resonances. Besides, we introduce a systematically geometrically convergent iterative method for the Costabel FEM-BEM coupling, with estimates on the convergence speed.
• 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
• 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. 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.
• 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.
• 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
• 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.