Publications

2012

• Operator Factorization for Multiple-Scattering Problems and an Application to Periodic Media
  
  • Coatléven Julien
  • Joly Patrick
  Communications in Computational Physics, Global Science Press, 2012, 11 (2), pp.303-318. This work concerns multiple-scattering problems for time-harmonic equations in a reference generic media. We consider scatterers that can be sources, obstacles or compact perturbations of the reference media. Our aim is to restrict the computational domain to small compact domains containing the scatterers. We use Robin-to-Robin (RtR) operators (in the most general case) to express boundary conditions for the interior problem. We show that one can always factorize the RtR map using only operators defined using single-scatterer problems. This factorization is based on a decomposition of the diffracted field, on the whole domain where it is defined. Assuming that there exists a good method for solving single-scatterer problems, it then gives a convenient way to compute RtR maps for a random number of scatterers.
• Propagation of guided waves through weak penetrable scatterers
  
  • Maurel Agnes
  • Mercier Jean-François
  Journal of the Acoustical Society of America, Acoustical Society of America, 2012, 131 (3), pp.1874-1889. The scattering of a scalar wave propagating in a waveguide containing weak penetrable scatterers is inspected in the Born approximation. The scatterers are of arbitrary shape and present a contrast both in density and in wavespeed (or bulk modulus), a situation that can be translated in the context of SH waves, water waves, or transverse electric/transverse magnetic polarized electromagnetic waves. For small size inclusions compared to the waveguide height, analytical expressions of the transmission and reflection coefficients are derived, and compared to results of direct numerical simulations. The cases of periodically and randomly distributed inclusions are considered in more detail, and compared with unbounded propagation through inclusions. Comparisons with previous results valid in the low frequency regime are proposed. © 2012 Acoustical Society of America. (10.1121/1.3682037)
  DOI : 10.1121/1.3682037
• On the use of sampling methods to identify cracks in acoustic waveguides
  
  • Bourgeois Laurent
  • Lunéville Éric
  Inverse Problems, IOP Publishing, 2012, 28 (10), pp.105011.1-105011.18. We consider the identification of cracks in an acoustic 2D/3D waveguide with the help of sampling methods such as the linear sampling method or the factorization method. A modal version of these sampling methods is used. Our paper emphasizes the fact that if one a priori knows the type of boundary condition which actually applies on the crack, then we shall adapt the formulation of our sampling method to such boundary conditions in order to improve the efficiency of the method. The need for such adaptation is proved theoretically and illustrated numerically with the help of 2D examples. We also show by using our modal formulation that the factorization method is applicable in a waveguide with the same data as the linear sampling method. © 2012 IOP Publishing Ltd. (10.1088/0266-5611/28/10/105011)
  DOI : 10.1088/0266-5611/28/10/105011
• An elementary introduction to the construction and the analysis of Perfectly Matched Layers for time domain wave propagation
  
  • Joly Patrick
  SeMA Journal: Boletin de la Sociedad Española de Matemática Aplicada, Springer, 2012, 57, pp.5-48.
• T-coercivity for scalar interface problems between dielectrics and metamaterials
  
  • Bonnet-Ben Dhia Anne-Sophie
  • Chesnel Lucas
  • Ciarlet Patrick
  ESAIM: Mathematical Modelling and Numerical Analysis, Société de Mathématiques Appliquées et Industrielles (SMAI) / EDP, 2012, 46, pp.1363-1387. Some electromagnetic materials have, in a given frequency range, an effective dielectric permittivity and/or a magnetic permeability which are real-valued negative coefficients when dissipation is neglected. They are usually called metamaterials. We study a scalar transmission problem between a classical dielectric material and a metamaterial, set in an open, bounded subset of R^d, with d=2,3. Our aim is to characterize occurences where the problem is well-posed within the Fredholm (or coercive + compact) framework. For that, we build some criteria, based on the geometry of the interface between the dielectric and the metamaterial. The proofs combine simple geometrical arguments with the approach of T-coercivity, introduced by the first and third authors and co-worker. Furthermore, the use of localization techniques allows us to derive well-posedness under conditions that involve the knowledge of the coefficients only near the interface. When the coefficients are piecewise constant, we establish the optimality of the criteria. (10.1051/m2an/2012006)
  DOI : 10.1051/m2an/2012006
• Helmholtz equation in periodic media with a line defect
  
  • Coatléven Julien
  Journal of Computational Physics, Elsevier, 2012, 231 (4), pp.1675-1704. We consider the Helmholtz equation in an unbounded periodic media perturbed by an unbounded defect whose structure is compatible with the periodicity of the underlying media. We exhibit a method coupling Dirichlet-to-Neumann maps with the Lippmann-Schwinger equation approach to solve this problem, where the Floquet-Bloch transform in the direction of the defect plays a central role. We establish full convergence estimates that makes the link between the rate of decay of a function and the good behavior of a quadrature rule to approximate the inverse Floquet-Bloch transform. Finally we exhibit a few numerical results to illustrate the efficiency of the method. © 2011 Elsevier Inc. (10.1016/j.jcp.2011.10.022)
  DOI : 10.1016/j.jcp.2011.10.022
• Evaluation of 3-D Singular and Nearly Singular Integrals in Galerkin BEM for Thin Layers
  
  • Lenoir Marc
  • Salles Nicolas
  SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2012, 36, pp.3057-3078. An explicit method for the evaluation of singular and near-singular integrals arising in three-dimensional Galerkin BEM is presented. It is based on a recursive reduction of the dimension of the integration domain leading to a linear combination of one-dimensional regular integrals, which can be exactly evaluated. This method has appealing properties in terms of reliability, precision, and flexibility. The results we present here are devoted to the case of thin layers for the Helmholtz equation, a situation where the panels are close and parallel, known to be difficult in terms of accuracy. Nevertheless, the method applies as well to two-dimensional BEM, secant planes, or even volume integral equations. A MATLAB implementation of the formulas presented here is available online. (10.1137/120866567)
  DOI : 10.1137/120866567
• The Variational Theory of Complex Rays for three-dimensional Helmholtz problems
  
  • Kovalevsky Louis
  • Ladevèze Pierre
  • Riou Hervé
  • Bonnet Marc
  Journal of Computational Acoustics, World Scientific Publishing, 2012, 20, pp.125021 (25 pages). This article proposes an extension of the Variational Theory of Complex Rays (VTCR) to three-dimensional linear acoustics, The VTCR is a Trefftz-type approach designed for mid-frequency range problems and has been previously investigated for structural dynamics and 2D acoustics. The proposed 3D formulation is based on a discretization of the amplitude portrait using spherical harmonics expansions. This choice of discretization allows to substantially reduce the numerical integration work by taking advantage of well-known analytical properties of the spherical harmonics. It also permits (like with the previous 2D Fourier version) an effective \emph{a priori} selection method for the discretization parameter in each sub-region, and allows to estimate the directivity of the pressure field by means of a natural definition of rescaled amplitude portraits. The accuracy and performance of the proposed formulation are demonstrated on a set of numerical examples that include results on an actual case study from the automotive industry. (10.1142/S0218396X1250021X)
  DOI : 10.1142/S0218396X1250021X
• An adaptive algorithm for cohesive zone model and arbitrary crack propagation
  
  • Chiaruttini Vincent
  • Geoffroy Dominique
  • Riolo Vincent
  • Bonnet Marc
  Revue Européenne de Mécanique Numérique/European Journal of Computational Mechanics, Hermès / Paris : Lavoisier, 2012, 21, pp.208-218. This paper presents an approach to the numerical simulation of crack propagation with cohesive models for the case of structures subjected to mixed mode loadings. The evolution of the crack path is followed by using an adaptive method: with the help of a macroscopic branching criterion based on the calculation of an energetic integral, the evolving crack path is remeshed as the crack evolves in the simulation. Special attention is paid to the unknown fields transfer approach that is crucial for the success of the computational treatment. This approach has been implemented in the finite element code Z-Set (jointly developed by Onera and Ecole des Mines) and is tested on two examples, one featuring a straight crack path and the other involving a complex crack propagation under critical monotonous loading monotonous. (10.1080/17797179.2012.744544)
  DOI : 10.1080/17797179.2012.744544
• On Simultaneous Identification of the Shape and Generalized Impedance Boundary Condition in Obstacle Scattering
  
  • Bourgeois Laurent
  • Chaulet Nicolas
  • Haddar Houssem
  SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2012, 34 (3), pp.A1824-A1848. We consider the inverse obstacle scattering problem of determining both the shape and the "equiva- lent impedance" from far field measurements at a fixed frequency. In this work, the surface impedance is represented by a second order surface differential operator (refer to as generalized impedance boundary condition) as opposed to a scalar function. The generalized impedance boundary condition can be seen as a more accurate model for effective impedances and is widely used in the scattering problem for thin coatings. Our approach is based on a least square optimization technique. A major part of our analysis is to characterize the derivative of the cost function with respect to the boundary and this complex surface impedance configuration. In particular, we provide an extension of the notion of shape derivative to the case where the involved impedance parameters do not need to be surface traces of given functions, which leads (in general) to a non-vanishing tangential boundary perturbation. The efficiency of considering this type of derivative is illustrated by several 2D numerical experiments based on a (classical) steepest descent method. The feasibility of retrieving both the shape and the impedance parameters is also discussed in our numerical experiments. (10.1137/110850347)
  DOI : 10.1137/110850347
• Usual Anderson localization restored in bilayered left- and right-handed structures
  
  • Maurel Agnès
  • Ourir Abdelwaheb
  • Mercier Jean-François
  • Pagneux Vincent
  Physical Review B: Condensed Matter and Materials Physics (1998-2015), American Physical Society, 2012, 85 (20). We present a study of the attenuation length in a one-dimensional array of alternating left- and right-handed materials in which both the permittivities and the permeabilities are disordered. This type of structure has been shown to present an anomaly in the attenuation length when only permeabilities are disordered. We derive a simple analytical expression of the attenuation length, when the disorder in the refraction index is due to perturbations in both the permeability and the permittivity. Our expression is able to explain the transition to the anomalous behavior when perturbation only in the permeability or only in the permittivity is considered. Besides, we show that the anomaly is dramatically affected when considering perturbations in permeability and permittivity. The coupling effects are able to restore the ordinary localization length. © 2012 American Physical Society. (10.1103/physrevb.85.205138)
  DOI : 10.1103/physrevb.85.205138
• Interior transmission eigenvalue problem for Maxwell's equations: The T-coercivity as an alternative approach
  
  • Chesnel Lucas
  Inverse Problems, IOP Publishing, 2012, 28 (6). In this paper, we examine the interior transmission problem for Maxwells equations in the case where both and , the physical parameters of the scattering medium, differ from 0 and 0 modelling the background medium. Using the T-coercivity method, we propose an alternative approach to the classical techniques to prove that this problem is of Fredholm type and that the so-called transmission eigenvalues form at most a discrete set. The T-coercivity approach allows us to deal with cases where 0 and 0 can change sign. We also provide results of localization and FaberKrahn-type inequalities for the transmission eigenvalues. © 2012 IOP Publishing Ltd. (10.1088/0266-5611/28/6/065005)
  DOI : 10.1088/0266-5611/28/6/065005
• Transparent boundary conditions for evolution equations in infinite periodic strips
  
  • Coatléven Julien
  SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2012, 34 (3), pp.1563-1583. We consider the solution of a generic equation $\gamma\rho(\mathbf{x})\partial^p_tu(\mathbf{x},t)-\Delta u(\mathbf{x},t) +V(\mathbf{x})u(\mathbf{x},t) = f(\mathbf{x},t)$, $\mathbf{x} = (x,y)$, for $t>0$, $p=1,2$ in a domain $\Omega$ which is infinite in $x$ and bounded in $y$. We assume that $f(\cdot,t)$ is supported for all $t>0$ in $\Omega_0 = \{\mathbf{x} \in \Omega \; | \; -a_- < x < a_+\}$ and that $\rho(\mathbf{x})$ and $V(\mathbf{x})$ are x-periodic in $\Omega \setminus \Omega_0$. We consider the associated $\theta$-scheme in time, to obtain a semidiscretized problem. We then show how to obtain for each time step exact boundary conditions on the vertical segments, $\Gamma_0^- = \{\mathbf{x}\in \Omega\; | \; x=-a_-\}$ and $\Gamma_0^+ = \{\mathbf{x}\in \Omega \;| \; x=a_+\}$, that will enable us to find the solution on $\Omega_0 \cup \Gamma_0^+ \cup \Gamma_0^-$. Then the solution can be extended in $\Omega$ in a straightforward manner from the values on $\Gamma_0^-$ and $\Gamma_0^+$. The method is based on the solution of local problems on a single periodicity cell, solved during an initialization step. The exact boundary conditions as well as the extension operators can be obtained for each time step through elementary computations using the solution of these local cell problems. (10.1137/110838030)
  DOI : 10.1137/110838030
• Uniform controllability of scalar conservation laws in the vanishing viscosity limit
  
  • Léautaud Matthieu
  SIAM Journal on Control and Optimization, Society for Industrial and Applied Mathematics, 2012, 50 (3), pp.1661-1699. We deal with viscous perturbations of scalar conservation laws on a bounded interval with a general flux function f and a small dissipation coefficient Ɛ. Acting on this system on both endpoints of the interval, we prove global exact controllability to constant states with nonzero speed. More precisely, we construct boundary controls so that the solution is driven to the targeted constant state, and we moreover require these controls to be uniformly bounded as Ɛ → 0+ in an appropriate space. For general (nonconvex) flux functions this can be done for sufficiently large time, and for convex fluxes f, we have a precise estimate on the minimal time needed to control. © 2012 Society for Industrial and Applied Mathematics. (10.1137/100803043)
  DOI : 10.1137/100803043
• A low frequency model for acoustic propagation in a 2D flow duct: numerical computation
  
  • Joubert Lauris
  • Joly Patrick
  Communications in Computational Physics, Global Science Press, 2012, 11 (2), pp.508-524. In this paper we study a low frequency model for acoustic propagation in a 2D flow duct. For some Mach profile flow, we are able to give a well-posedness theorem. Its proof relies on a quasi-explicit expression of the solution which provides us an efficient numerical method. We give and comment numerical results for particular linear, tangent and quadratic profiles. Finally, we give a numerical validation of our asymptotic model.
• Giens 2011
  
  • Bonnet Marc
  • Cornuault Christian
  • Pagano Stéphane
  , 2012.
• Giens 2011
  
  • Bonnet Marc
  • Cornuault Christian
  • Pagano Stéphane
  , 2012.
• Application of the multi-level time-harmonic fast multipole BEM to 3-D visco-elastodynamics
  
  • Grasso Eva
  • Chaillat Stéphanie
  • Bonnet Marc
  • Semblat Jean-François
  Engineering Analysis with Boundary Elements, Elsevier, 2012, 36, pp.744-758. This article extends previous work by the authors on the single- and multi-domain time-harmonic elastodynamic multi-level fast multipole BEM formulations to the case of weakly dissipative viscoelastic media. The underlying boundary integral equation and fast multipole formulations are formally identical to that of elastodynamics, except that the wavenumbers are complex-valued due to attenuation. Attention is focused on evaluating the multipole decomposition of the viscoelastodynamic fundamental solution. A damping-dependent modification of the selection rule for the multipole truncation parameter, required by the presence of complex wavenumbers, is proposed. It is empirically adjusted so as to maintain a constant accuracy over the damping range of interest in the approximation of the fundamental solution, and validated on numerical tests focusing on the evaluation of the latter. The proposed modification is then assessed on 3D single-region and multi-region visco-elastodynamic examples for which exact solutions are known. Finally, the multi-region formulation is applied to the problem of a wave propagating in a semi-infinite medium with a lossy semi-spherical inclusion (seismic wave in alluvial basin). These examples involve problem sizes of up to about $3\,10^{5}$ boundary unknowns. (10.1016/j.enganabound.2011.11.015)
  DOI : 10.1016/j.enganabound.2011.11.015