Publications

2024

• Exploring low-rank structure for an inverse scattering problem with far-field data
  
  • Zhou Yuyuan
  • Audibert Lorenzo
  • Meng Shixu
  • Zhang Bo
  SIAM Journal on Applied Mathematics, Society for Industrial and Applied Mathematics, 2024. The inverse scattering problem exhibits an inherent low-rank structure due to its ill-posed nature; however developing low-rank structures for the inverse scattering problem remains challenging. In this work, we introduce a novel low-rank structure tailored for solving the inverse scattering problem. The particular low-rank structure is given by the generalized prolate spheroidal wave functions, computed stably and accurately via a Sturm-Liouville problem. We first process the far-field data to obtain a post-processed data set within a disk domain. Subsequently, the post-processed data are projected onto a low-rank space given by the low-rank structure. The unknown is approximately solved in this low-rank space, by dropping higher-order terms. The low-rank structure leads to a H\"{o}lder-logarithmic type stability estimate for arbitrary unknown functions, and a Lipschitz stability estimate for unknowns belonging to a finite dimensional low-rank space. Various numerical experiments are conducted to validate its performance, encompassing assessments of resolution capability, robustness against randomly added noise and modeling errors, and demonstration of increasing stability.
• Time-harmonic wave propagation in junctions of two periodic half-spaces
  
  • Amenoagbadji Pierre
  • Fliss Sonia
  • Joly Patrick
  , 2024. We are interested in the Helmholtz equation in a junction of two periodic half-spaces. When the overall medium is periodic in the direction of the interface, Fliss and Joly (2019) proposed a method which consists in applying a partial Floquet-Bloch transform along the interface, to obtain a family of waveguide problems parameterized by the Floquet variable. In this paper, we consider two model configurations where the medium is no longer periodic in the direction of the interface. Inspired by the works of Gérard-Varet and Masmoudi (2011, 2012), and Blanc, Le Bris, and Lions (2015), we use the fact that the overall medium has a so-called quasiperiodic structure, in the sense that it is the restriction of a higher dimensional periodic medium. Accordingly, the Helmholtz equation is lifted onto a higher dimensional problem with coefficients that are periodic along the interface. This periodicity property allows us to adapt the tools previously developed for periodic media. However, the augmented PDE is elliptically degenerate (in the sense of the principal part of its differential operator) and thus more delicate to analyse.
• About entropy penalized stochastic control problems
  
  • Bourdais Thibaut
  , 2024. This thesis focuses on stochastic optimal control problems with constraints on the marginal laws of the state process. These problems have constituted a very active research area in the past few years since they offer many practical and theoretical perspectives. In particular, this work is motivated by demand side management applications in power systems. More specifically, we are interested in controlling the overall electrical consumption of agents on a one-day period while imposing that the terminal marginal law of the state of the agents is equal to the initial law. The interest in this formulation is to define, in a seasonal environment, a periodic control procedure which can be reconducted on each period.This thesis explores an original way of solving stochastic optimal control problems with constraints in law based on a reformulation as constrained optimization problems on the space of probability measures. The decision variable is then splitted into two decision variables, each one considering only a specific part of the constraints while the deviation between these two probability measures is penalized by adding a relative entropy term. We take advantage of this penalized version to propose an alternating minimization procedure to approximate a solution to the original problem. This procedure involves solving sequentially two simple subproblems, each one consisting in minimizing the objective function over one variable while the other is fixed. The first subproblem amounts to a pointwise minimization of a running cost function. The solution of the second subproblem can be expressed as the so called exponential twist of a reference (Markovian) probability measure whose study constitutes a substantial part of this thesis. This alternating procedure is proved to converge to an approximate solution of the original control problem under various assumptions allowing for instance non convex running costs with respect to the control variable, as well as jump diffusions dynamics.
• Viscosity solutions of centralized control problems in measure spaces
  
  • Aussedat Averil
  • Jerhaoui Othmane
  • Zidani Hasnaa
  ESAIM: Control, Optimisation and Calculus of Variations, EDP Sciences, 2024, 30, pp.91. This work focuses on a control problem in the Wasserstein space of probability measures over Rd. Our aim is to link this control problem to a suitable Hamilton-Jacobi-Bellman (HJB) equation. We explore a notion of viscosity solution using test functions that are locally Lipschitz and locally semiconvex or semiconcave functions. This regularity allows to define a notion of viscosity and a Hamiltonian function relying on directional derivatives. Using a generalization of Ekeland's principle, we show that the corresponding HJB equation admits a comparison principle, and deduce that the value function is the unique solution in this viscosity sense. The PDE tools are developed in the general framework of Measure Differential Equations. (10.1051/cocv/2024081)
  DOI : 10.1051/cocv/2024081
• Construction and analysis of spectral signatures for defects in complex media
  
  • Pourre Fabien
  , 2024. Concrete is widely used in construction, particularly in the reactor building of nuclear power plants. Monitoring its evolution and identifying any defects that could compromise its proper functioning can be achieved through non-destructive testing. Concrete is composed of aggregates, but their high concentration and proximity pose challenges for classical methods, such as the Linear Sampling Method, which fail to produce quantitative results or exploitable images.The objective of this thesis is to build an imaging algorithm that can estimatethe density of the aggregates and recover the local distribution of those smallinhomogeneities.To address this challenge, inspiration was drawn from a monotonicity property of the Transmissioneigenvalues. They correspond to the frequencies for which an incident wave exists such that the scattered field is trivial outside the scatterers.Instead of comparing the scattered field to the vacuum, this thesis introduces a new approach: comparing it, at a fixed wavenumber, to a numerical scattering problem, referred to as the background. This led to the introduction of a new class of eigenvalues known as the f-averaged Steklov eigenvalues. Each of these eigenvalues is associatedwith an artificial background problem that contains a resonator and is solution to a simple spectral problem inside the resonator. In addition, each f-averaged Steklov eigenvalue is monotonically increasing with respect to the number of inhomogeneities inside that resonator. These spectral signatures can be recovered from the data. By computing theseeigenvalues for various positions of the resonator, we can estimate the variation ofthe local density of the inhomogeneities in an unknown medium.
• A rounding and clustering-based exact algorithm for the p-center problem
  
  • Alès Zacharie
  • Duran-Mateluna Cristian
  • Elloumi Sourour
  Computers and Operations Research, Elsevier, 2024, 183, pp.107185. The p-center problem consists of selecting p facilities from a set of possible sites and allocating a set of clients to them in such a way that the maximum distance between a client and the facility to which it is allocated is minimized. This paper proposes a new scalable exact solution algorithm based on client clustering and an iterative distance rounding procedure. The client clustering enables to initialize and update a subset of clients for which the p-center problem is iteratively solved. The rounding drastically reduces the number of distinct distances considered at each iteration. Our algorithm is tested on 396 benchmark instances with up to 1.9 million clients and facilities. Our results show that our approach outperforms existing methods run on the same computer except when p is smaller than 5. In this case, however, we optimally solve all instances in less than 2 minutes on average. (10.1016/j.cor.2025.107185)
  DOI : 10.1016/j.cor.2025.107185
• Étude analytique et numérique de problèmes inverses en diffraction acoustique pour la conception de microphones spatiaux
  
  • Lerévérend Dorian
  , 2024. CONTEXTE : Ces travaux de thèse sont motivés par un besoin concret identifié par l'association Mon Cartable Connecté. Dans le cadre de l'amélioration de son dispositif de téléprésence pour aider à la scolarisation à distance d'enfants hospitalisés, l'association souhaite implémenter un son spatialisé et immersif. En raison de l'impact prouvé d'un tel dispositif sur les interactions sociales entre élèves, il faudra que la solution proposée n'évoque pas une forme humaine. Enfin, elle devra respecter des contraintes d'encombrement et de coût. Parmi les technologies existantes, nous pouvons citer les microphones ambisoniques (Zoom H3-VR, Zylia Pro, ...) et les microphones en formes de têtes humaines (Kemar, Neumann KU-100, ...). Les premiers produisent un signal spatialisé (avec de nombreux canaux) mais délivrent un son neutre, ce qui limite l'immersion. Les seconds reproduisent l'écoute naturelle humaine sans aucun traitement numérique mais génèrent seulement un son binaural (avec deux canaux). Nous proposons donc de réunir ces deux technologies dans un objet dont la forme serait différente d'une tête humaine. ENJEUX : Une première difficulté est le fait que la technologie ambisonique repose sur les fonctions harmoniques surfaciques de l'appareil d'enregistrement utilisé. En effet, ces dernières sont connues analytiquement pour des sphères et des ellipses uniquement. Nous décidons de laisser cette remarque de côté et nous nous concentrons sur la seconde difficulté : la question de l'existence d'obstacles ayant des géométries différentes et présentant la même réponse impulsionnelle est un problème ouvert. Nous pouvons tout de même définir un modèle simplifié afin de proposer des méthodes de résolution. Nous considérons trois grandeurs principales pour caractériser l'acoustique d'un objet : sa géométrie, ses propriétés d'absorption et la position de ses points d'écoute (tympans ou microphones). Il existe de nombreuses bases de données donnant libre accès à des géométries de têtes et aux réponses fréquentielles associées. Nous pouvons alors définir plusieurs problèmes d'identification des valeurs d'un ou plusieurs de ces paramètres. Nous les résolvons en minimisant une fonction de coût par descente de gradient. STRUCTURE DU MANUSCRIT : Nous commençons par supposer que la tête est un obstacle impénétrable. Nous développons un solveur pour le problème direct de diffraction acoustique avec Matlab. Cet algorithme performant est basé sur l'implémentation des matrices hiérarchiques dans Gypsilab (développé en Matlab par Matthieu Aussal). Il nous permet de calculer rapidement des réponses en fréquence (HRTFs) réalistes avec différentes conditions d'impédance jusqu'à 20kHz. Voyant que l'erreur entre les résultats numériques et les mesures expérimentales peut être réduite par un choix judicieux de condition de bord, nous identifions l'impédance optimale de têtes humaines jusqu'à 15kHz à partir de données expérimentales sans phase. Se pose alors la question des formes équivalentes. Nous déterminons des triplets (forme, impédance, tympans) donnant des HRTFs assez proches de celles de têtes humaines jusqu'à 4kHz. Enfin, nous supposons que la tête est un obstacle pénétrable et changeons notre modèle. Nous transposons les développements réalisés précédemment à ce nouveau cadre. Ainsi, nous identifions les indices de réfraction de têtes humaines jusqu'à 4kHz. Nous terminons avec la détermination numérique d'obstacles équivalents définis par des triplets (forme, indice de réfraction, tympans) dont les réponses fréquentielles approchent des HRTFs mesurées expérimentalement jusqu'à 4kHz.
• Construction of transparent conditions for electromagnetic waveguides
  
  • Bonnet-Ben Dhia Anne-Sophie
  • Chesnel Lucas
  • Fliss Sonia
  • Parigaux Aurélien
  , 2024. We are interested in the numerical resolution of diffraction problems in closed electromagnetic waveguides by means of finite elements methods. To proceed, we need to truncate the domain and design adapted transparent conditions on the artificial boundary to avoid spurious reflections. When the guide is homogeneous in the transverse section, this can be done by writing an Electric-to-Magnetic condition based on a modal decomposition of the field. The latter takes a rather simple form thanks to the orthogonality of transverse modes. For guides that are heterogeneous in the transverse section, the transverse modes are no longer orthogonal but satisfy bi-orthogonality relations linked to the Poynting energy flux. Modal decompositions are more delicate to derive and it may happen that certain modes have phase and group velocities of different sign, which prevents the use of Perfectly Matched Layers. Adapting techniques already developed in elasticity, we derive a new transparent condition based on a Poynting-to-Magnetic operator with overlap. To illustrate the method, we present numerical results obtained with Nédélec finite elements using the XLiFE++ library.
• Notes de cours sur les équations de Maxwell et leur approximation
  
  • Ciarlet Patrick
  , 2024, pp.151.
• SpinDoctor-IVIM: A virtual imaging framework for intravoxel incoherent motion MRI
  
  • Lashgari Mojtaba
  • Yang Zheyi
  • Bernabeu Miguel O
  • Li Jing-Rebecca
  • Frangi Alejandro F
  Medical Image Analysis, Elsevier, 2024, 99, pp.103369. <div><p>Intravoxel incoherent motion (IVIM) imaging is increasingly recognised as an important tool in clinical MRI, where tissue perfusion and diffusion information can aid disease diagnosis, monitoring of patient recovery, and treatment outcome assessment. Currently, the discovery of biomarkers based on IVIM imaging, similar to other medical imaging modalities, is dependent on long preclinical and clinical validation pathways to link observable markers derived from images with the underlying pathophysiological mechanisms. To speed up this process, virtual IVIM imaging is proposed. This approach provides an efficient virtual imaging tool to design, evaluate, and optimise novel approaches for IVIM imaging. In this work, virtual IVIM imaging is developed through a new finite element solver, SpinDoctor-IVIM, which extends SpinDoctor, a diffusion MRI simulation toolbox. SpinDoctor-IVIM simulates IVIM imaging signals by solving the generalised Bloch-Torrey partial differential equation. The input velocity to SpinDoctor-IVIM is computed using HemeLB, an established Lattice Boltzmann blood flow simulator. Contrary to previous approaches, SpinDoctor-IVIM accounts for volumetric microvasculature during blood flow simulations, incorporates diffusion phenomena in the intravascular space, and accounts for the permeability between the intravascular and extravascular spaces. The above-mentioned features of the proposed framework are illustrated with simulations on a realistic microvasculature model.</p></div> (10.1016/j.media.2024.103369)
  DOI : 10.1016/j.media.2024.103369
• Coupled Boundary Element and Finite Volume Methods for Modeling Fluid-Induced Seismicity in Fault Networks within Low-Permeability Rocks
  
  • Romanet Pierre
  • Scuderi Marco Maria
  • Ampuero Jean-Paul
  • Chaillat Stéphanie
  • Cappa Frederic
  Geophysical Journal International, Oxford University Press (OUP), 2024, 243 (3). To better understand the mechanics of injection-induced seismicity, we developed a two-dimensional numerical code to simulate both seismic and aseismic slip on non-planar faults and fault networks driven by fluid diffusion along permeable faults. Our approach integrates a boundary element method to model fault slip governed by rate-and-state friction with a finite volume method for simulating fluid diffusion along fault networks. We demonstrate the method's capabilities with two illustrative examples: (1) fluid injection inducing slow slip on a primary rough, rate-strengthening fault, which subsequently triggers microseismicity on secondary, smaller faults, and (2) fluid injection on a single fault in a network of intersecting faults, leading to fluid diffusion and reactivation of slip throughout the network. In both cases, the simulated slow slip migrates more rapidly than the fluid pressure diffusion front. The observed migration patterns of microseismicity in the first example and slow slip in the second example resemble diffusion processes but involve diffusivity values that differ significantly from the fault hydraulic diffusivity. These results support the conclusion that the microseismicity front is not a direct proxy for the fluid diffusion front and cannot be used to directly infer hydraulic diffusivity, consistently with some decametric scale in-situ experiments of fault activation under controlled conditions. This work highlights the importance of distinguishing between mechanical and hydrological processes in the analysis of induced seismicity, providing a powerful tool for improving our understanding of fault behavior in response to fluid injection, in particular when a network of faults is involved. (10.1093/gji/ggaf377)
  DOI : 10.1093/gji/ggaf377
• A posteriori error estimates for the DD+$L^2$ jumps method on the Neutron Diffusion equations
  
  • Ciarlet Patrick
  • Do Minh-Hieu
  • Gervais Mario
  • Madiot François
  , 2024. We analyse a posteriori error estimates for the discretization of the neutron diffusion equations with a Domain Decomposition Method, the so-called DD+$L^2$ jumps method. We provide guaranteed and locally efficient estimators on a base block equation, the one-group neutron diffusion equation. Classically, one introduces a Lagrange multiplier to account for the jumps on the interface. This Lagrange multiplier is used for the reconstruction of the physical variables. Remarkably, no reconstruction of the Lagrange multiplier is needed to achieve the optimal a posteriori estimates.
• High-order numerical integration on self-affine sets
  
  • Joly Patrick
  • Kachanovska Maryna
  • Moitier Zoïs
  , 2024. We construct an interpolatory high-order cubature rule to compute integrals of smooth functions over self-affine sets with respect to an invariant measure. The main difficulty is the computation of the cubature weights, which we characterize algebraically, by exploiting a self-similarity property of the integral. We propose an \( h \)-version and a \( p \)-version of the cubature, present an error analysis and conduct numerical experiments.
• Fading regularization method for an inverse boundary value problem associated with the biharmonic equation
  
  • Boukraa Mohamed Aziz
  • Caillé Laëtitia
  • Delvare Franck
  Journal of Computational and Applied Mathematics, Elsevier, 2024, 457, pp.116285. In this paper, we propose a numerical algorithm that combines the fading regularization method with the method of fundamental solutions (MFS) to solve a Cauchy problem associated with the biharmonic equation. We introduce a new stopping criterion for the iterative process and compare its performance with previous criteria. Numerical simulations using MFS validate the accuracy of this stopping criterion for both compatible and noisy data and demonstrate the convergence, stability, and efficiency of the proposed algorithm, as well as its ability to deblur noisy data. (10.1016/j.cam.2024.116285)
  DOI : 10.1016/j.cam.2024.116285
• $C^{ 0,1}$ -Itô chain rules and generalized solutions of parabolic PDEs
  
  • Ciccarella Carlo
  • Russo Francesco
  , 2024. In this paper we first establish an It\^o formula for a finite quadratic variation process $X$ expanding $f(t,X_t),$ when $f$ is of class $C^2$ in space and is absolutely continuous in time. Second, via a Fukushima-Dirichlet decomposition we obtain an explicit chain rule for $f(t,X_t)$, when $X$ is a continuous semimartingale and $f$ is a ``quasi-strong solution'' (in the sense of approximation of classical solutions) of a parabolic PDE.
• Propagation of ultrasounds in random multi-scale media and effecitve speed of sound estimation
  
  • Goepfert Quentin
  , 2024. Ultrasounds are widely used in medical imaging modalities. Originally, the ultrasound devices were built to image the internal structure of the tissues. In recent years, a change of paradigm operated and the goal is now also to assess physical parameters that can be used for medical diagnosis.The speed of acoustic waves inside soft tissues can be used for diagnosis of breast cancers or hepatic steatosis. Moreover, it determines the quality of the tomographic reconstruction of the tissues. Indeed, the images are usually computed by backpropagating the measured echoes at the speed of sound in water. However, the discrepancy between the speed of sound in water and the actual speed of sound inside the tissues results in nonphysical artifacts on the image.In order to establish a quantitative estimator of the propagation speed of sound inside the soft tissues, it is necessary to deeply understand the scattering of the medium. It is commonly admitted that the backscattered echoes are produced by numerous unresolved scatterers inside the medium (cell nuclei, mitochondria...). The scattering is then often modeled by the Born approximation. However, this model does not capture the variation of the effective speed of sound inside the tissue due to the unresolved scatterers. The goal of this thesis is thus to establish a propagation model that takes into account the variations of the effective speed of sound inside the tissues. Then, we will theoretically study the estimators previously introduced by Alexandre Aubry in his work.The tissue is here modeled as a bounded homogeneous mediumin which lie unresolved scatterers. As their distribution is unknown and inaccessible, their number and position is modeled as a random process. To obtain a simple form of the backscattered field, the techniques and tools developed for the quantitative stochastic homogenization theory will be used and a high-order asymptotic expansion will be proven.An asymptotic analysis of the imaging functional is carried out by using the high-order asymptotic expansion. Furthermore, the theoretical study of the estimators introduced by Alexandre Aubry and his team confirms and justifies some of the experimental results. In particular, it is possible to recover the effective speed of sound by a local spatial average of the imaging function.Numerical simulation supports each and every major result proven in this thesis.
• Guided modes in a hexagonal periodic graph like domain
  
  • Delourme Bérangère
  • Fliss Sonia
  Multiscale Modeling and Simulation: A SIAM Interdisciplinary Journal, Society for Industrial and Applied Mathematics, 2024, 22 (3), pp.1196-1245. This paper deals with the existence of guided waves and edge states in particular two-dimensional media obtained by perturbing a reference periodic medium with honeycomb symmetry. This reference medium is a thin periodic domain (the thickness is denoted δ &gt; 0) with an hexagonal structure, which is close to an honeycomb quantum graph. In a first step, we show the existence of Dirac points (conical crossings) at arbitrarily large frequencies if δ is chosen small enough. We then perturbe the domain by cutting the perfectly periodic medium along the so-called zigzag direction, and we consider either Dirichlet or Neumann boundary conditions on the cut edge. In the two cases, we prove the existence of edges modes as well as their robustness with respect to some perturbations, namely the location of the cut and the thickness of the perturbed edge. In particular, we show that different locations of the cut lead to almost-non dispersive edge states, the number of locations increasing with the frequency. All the results are obtained via asymptotic analysis and semi-explicit computations done on the limit quantum graph. Numerical simulations illustrate the theoretical results. (10.1137/23M1600177)
  DOI : 10.1137/23M1600177
• The Sphericity Paradox and the Role of Hoop Stresses in Free Subduction on a Sphere
  
  • Chaillat Stéphanie
  • Gerardi Gianluca
  • Li Yida
  • Chamolly Alexander
  • Li Zhong‐hai
  • Ribe Neil M.
  Journal of Geophysical Research : Solid Earth, American Geophysical Union, 2024, 129 (9), pp.e2024JB029500. Oceanic plates are doubly curved spherical shells, which influences how they respond to loading during subduction. Here we study a viscous fluid model for gravity‐driven subduction of a shell comprising a spherical plate and an attached slab. The shell is 100–1,000 times more viscous than the upper mantle. We use the boundary‐element method to solve for the flow. Solutions of an axisymmetric model show that the effect of sphericity on the flexure of shells is greater for smaller shells that are more nearly flat (the “sphericity paradox”). Both axisymmetric and three‐dimensional models predict that the deviatoric membrane stress in the slab should be dominated by the longitudinal normal stress (hoop stress), which is typically about twice as large as the downdip stress and of opposite sign. Our models also predict that concave‐landward slabs can exhibit both compressive and tensile hoop stress depending on the depth, whereas the hoop stress in convex slabs is always compressive. We test these two predictions against slab shape and earthquake focal mechanism data from the Mariana subduction zone, assuming that the deviatoric stress in our flow models corresponds to that implied by centroid moment tensors. The magnitude of the hoop stress exceeds that of the downdip stress for about half the earthquakes surveyed, partially verifying our first prediction. Our second prediction is supported by the near‐absence of earthquakes under tensile hoop stress in the portion of the slab having convex geometry. (10.1029/2024JB029500)
  DOI : 10.1029/2024JB029500
• Notes de cours sur les méthodes variationnelles pour l'analyse et la résolution de problèmes non coercifs
  
  • Bonnet-Ben Dhia Anne-Sophie
  • Ciarlet Patrick
  , 2024.
• The T-coercivity approach for mixed problems
  
  • Barré Mathieu
  • Ciarlet Patrick
  Comptes Rendus. Mathématique, Académie des sciences (Paris), 2024. Classically, the well-posedness of variational formulations of mixed linear problems is achieved through the inf-sup condition on the constraint. In this note, we propose an alternative framework to study such problems by using the T-coercivity approach to derive a global inf-sup condition. Generally speaking, this is a constructive approach that, in addition, drives the design of suitable approximations. As a matter of fact, the derivation of the uniform discrete inf-sup condition for the approximate problems follows easily from the study of the original problem. To support our view, we solve a series of classical mixed problems with the T-coercivity approach. Among others, the celebrated Fortin Lemma appears naturally in the numerical analysis of the approximate problems.
• Accelerated iterative DG finite element solvers for large-scale time-harmonic acoustic problems
  
  • Modave Axel
  , 2024. Finite element methods are widely used to solve time-harmonic wave propagation problems, but solving large cases can be extremely difficult even with the computational power of parallel computers. In this work, the linear system resulting from the finite element discretization is solved with iterative solution methods, which are efficient in parallel but can require a large number of iterations. In standard discontinuous Galerkin (DG) methods, the numerical solution is discontinuous at the interfaces between the elements. In hybridizable DG methods, additional unknowns are introduced at the interfaces between the finite elements, and the physical unknowns are eliminated from the global system, resulting in a hybridized system. We have recently proposed a new strategy, called CHDG, where the additional unknowns correspond to transmission variables, whereas in the standard approach they are numerical fluxes. This strategy improves the properties of the hybridized system for faster iterative solution procedures. In this talk, we present and study a 3D CHDG implementation with nodal finite element basis functions. The resulting scheme has properties amenable to efficient parallel computing. Numerical results are presented to validate the method, and preliminary 3D computational results are proposed. (10.3397/IN_2024_2877)
  DOI : 10.3397/IN_2024_2877
• Far-field sound field estimation using robotized measurements and the boundary elements method
  
  • Pascal Caroline
  • Marchand Pierre
  • Chapoutot Alexandre
  • Doaré Olivier
  , 2024, 270 (11), pp.816-827. Sound Field Estimation (SFE) is a numerical technique widely used to identify and reconstruct the acoustic fields radiated by unknown structures. In particular, SFE proves to be useful when data is only available close to the source, but information in the whole space is required. However, the practical implementation of this method is still hindered by two major drawbacks: the lack of efficient implementation of existing numerical methodologies, and the time-consuming and tedious roll-out of acoustic measurements. This paper aims to provide a solution to both issues. First, the measurements step is fully automated by using a robotic arm, able to accurately gather geometric and acoustic data without any human assistance. In this matter, a particular attention has been paid to the impact of the robot on the acoustic pressure measurements. The sound field prediction is then tackled using the Boundary Element Method (BEM), and implemented using the FreeFEM++ BEM library. Numerically simulated measurements have allowed us to assess the method accuracy, and the overall solution has been successfully tested using actual robotized measurements of an unknown loudspeaker (10.3397/IN_2024_2661)
  DOI : 10.3397/IN_2024_2661
• Computation of Green's functions for the acoustic scattering by an elastic structure excited by a turbulent flow in water
  
  • Pacaut Louise
  • Serre Gilles
  • Mercier Jean-François
  • Chaillat Stéphanie
  • Cotté Benjamin
  , 2024, 270 (5), pp.5995-6006. To model the hydrodynamic noise produced by an elastic ship hull or propeller excited by a turbulent boundary layer, we need an efficient method to compute the acoustic scattering by an elastic body surrounded by a fluid. In 3D, Boundary Element Methods (BEM) are used to reduce the computational costs, for both the fluid and the elastic body. A natural way to compute the boundary integral representation (BIR) of the sound pressure is to use formulations based on the free space acoustic and elastic Green's functions. However, since the turbulent flow along the elastic body is known only statistically, the use of these Green's functions would be too expensive. A remedy is to compute a Green's function adapted to the physical problem, thus satisfying the transmission conditions of the fluid-structure problem. This so-called "tailored Green's function" is determined by solving a coupled acoustic-elastic problem with the BEM, and leads to a simplified BIR of the sound pressure compatible with a stochastic source term. We first validate the computation of the tailored Green's function over a classic spherical geometry. Then we compare the scattering of multiple quadrupoles by elastic or rigid NACA0012 profiles. (10.3397/IN_2024_3671)
  DOI : 10.3397/IN_2024_3671
• Shape and parameter identification by the linear sampling method for a restricted Fourier integral operator
  
  • Audibert Lorenzo
  • Meng Shixu
  Inverse Problems, IOP Publishing, 2024, 40 (9), pp.095007. In this paper we provide a new linear sampling method based on the same data but a different definition of the data operator for two inverse problems: the multi-frequency inverse source problem for a fixed observation direction and the Born inverse scattering problems. We show that the associated regularized linear sampling indicator converges to the average of the unknown in a small neighborhood as the regularization parameter approaches to zero. We develop both a shape identification theory and a parameter identification theory which are stimulated, analyzed, and implemented with the help of the prolate spheroidal wave functions and their generalizations. We further propose a prolate-based implementation of the linear sampling method and provide numerical experiments to demonstrate how this linear sampling method is capable of reconstructing both the shape and the parameter. (10.1088/1361-6420/ad5e18)
  DOI : 10.1088/1361-6420/ad5e18
• Multiscale modeling for a class of high-contrast heterogeneous sign-changing problems
  
  • Ye Changqing
  • Jin Xingguang
  • Ciarlet Patrick
  • Chung Eric T.
  , 2024. 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 (CEM-GMsFEM) 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.