Publications

2021

• Local transparent boundary conditions for wave propagation in fractal trees (I). Method and numerical implementation
  
  • Joly Patrick
  • Kachanovska Maryna
  SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2021. This work is dedicated to the construction and analysis of high-order transparentboundary conditions for the weighted wave equation on a fractal tree, which models sound propaga-tion inside human lungs. This article follows the works [9, 6], aimed at the analysis and numerical treatment of the model, as well as the construction of low-order and exact discrete boundary conditions. The method suggested in the present work is based on the truncation of the meromorphicseries that represents the symbol of the Dirichlet-to-Neumann operator, in the spirit of the absorbingboundary conditions of B. Engquist and A. Majda. We analyze its stability and convergence, as wellas present computational aspects of the method. Numerical results confirm theoretical findings (10.1137/20M1362334)
  DOI : 10.1137/20M1362334
• An automatic PML for acoustic finite element simulations in convex domains of general shape
  
  • Bériot Hadrien
  • Modave Axel
  International Journal for Numerical Methods in Engineering, Wiley, 2021, 122 (5), pp.1239-1261. This article addresses the efficient finite element solution of exterior acoustic problems with truncated computational domains surrounded by perfectly matched layers (PMLs). The PML is a popular nonreflecting technique that combines accuracy, computational efficiency, and geometric flexibility. Unfortunately, the effective implementation of the PML for convex domains of general shape is tricky because of the geometric parameters that are required to define the PML medium. In this work, a comprehensive implementation strategy is proposed. This approach, which we call the automatically matched layer (AML) implementation, is versatile and fully automatic for the end‐user. With the AML approach, the mesh of the layer is extruded, the required geometric parameters are automatically obtained during the extrusion step, and the practical implementation relies on a simple modification of the Jacobian matrix in the elementwise integrals. The AML implementation is validated and compared with other implementation strategies using numerical benchmarks in two and three dimensions, considering computational domains with regular and nonregular boundaries. A three‐dimensional application with a generally shaped domain generated using a convex hull is proposed to illustrate the interest of the AML approach for realistic industrial cases. (10.1002/nme.6560)
  DOI : 10.1002/nme.6560
• Stability and Convergence Analysis of Time-domain Perfectly Matched Layers for The Wave Equation in Waveguides
  
  • Bécache Eliane
  • Kachanovska Maryna
  SIAM Journal on Numerical Analysis, Society for Industrial and Applied Mathematics, 2021. This work is dedicated to the proof of stability and convergence of the Bérenger's perfectly matched layers in the waveguides for an arbitrary L ∞ damping function. The proof relies on the Laplace domain techniques and an explicit representation of the solution to the PML problem in the waveguide. A bound for the PML error that depends on the absorption parameter and the length of the PML is presented. Numerical experiments confirm the theoretical findings. (10.1137/20M1330543)
  DOI : 10.1137/20M1330543
• Homogenization of Maxwell's equations and related scalar problems with sign-changing coefficients
  
  • Bunoiu Renata
  • Chesnel Lucas
  • Ramdani Karim
  • Rihani Mahran
  Annales de la Faculté des Sciences de Toulouse. Mathématiques., Université Paul Sabatier _ Cellule Mathdoc, 2021, 30 (5), pp.1075-1119. In this work, we are interested in the homogenization of time-harmonic Maxwell's equations in a composite medium with periodically distributed small inclusions of a negative material. Here a negative material is a material modelled by negative permittivity and permeability. Due to the sign-changing coefficients in the equations, it is not straightforward to obtain uniform energy estimates to apply the usual homogenization techniques. The goal of this article is to explain how to proceed in this context. The analysis of Maxwell's equations is based on a precise study of two associated scalar problems: one involving the sign-changing permittivity with Dirichlet boundary conditions, another involving the sign-changing permeability with Neumann boundary conditions. For both problems, we obtain a criterion on the physical parameters ensuring uniform invertibility of the corresponding operators as the size of the inclusions tends to zero. In the process, we explain the link existing with the so-called Neumann-Poincaré operator, complementing the existing literature on this topic. Then we use the results obtained for the scalar problems to derive uniform energy estimates for Maxwell's system. At this stage, an additional difficulty comes from the fact that Maxwell's equations are also sign-indefinite due to the term involving the frequency. To cope with it, we establish some sort of uniform compactness result. (10.5802/afst.1694)
  DOI : 10.5802/afst.1694
• Optimal slip velocities of micro-swimmers with arbitrary axisymmetric shapes
  
  • Guo Hanliang
  • Zhu Hai
  • Liu Ruowen
  • Bonnet Marc
  • Veerapaneni Shravan
  Journal of Fluid Mechanics, Cambridge University Press (CUP), 2021, 910, pp.A26. This article presents a computational approach for determining the optimal slip velocities on any given shape of an axisymmetric micro-swimmer suspended in a viscous fluid. The objective is to minimize the power loss to maintain a target swimming speed, or equivalently to maximize the efficiency of the micro-swimmer. Owing to the linearity of the Stokes equations governing the fluid motion, we show that this PDE-constrained optimization problem reduces to a simpler quadratic optimization problem, whose solution is found using a high-order accurate boundary integral method. We consider various families of shapes parameterized by the reduced volume and compute their swimming efficiency. {Among those, prolate spheroids were found to be the most efficient micro-swimmer shapes for a given reduced volume. We propose a simple shape-based scalar metric that can determine whether the optimal slip on a given shape makes it a pusher, a puller, or a neutral swimmer.} (10.1017/jfm.2020.969)
  DOI : 10.1017/jfm.2020.969
• Variational Methods for Acoustic Radiation in a Duct with a Shear Flow and an Absorbing Boundary
  
  • Mercier Jean-François
  SIAM Journal on Applied Mathematics, Society for Industrial and Applied Mathematics, 2021, 81 (6), pp.2658-2683. The well-posedness of the acoustic radiation in a 2D duct in presence of both a shear flow and an absorbing wall described by the Myers boundary condition is studied thanks to variational methods. Without flow the problem is found well-posed for any impedance value. The presence of a flow complicates the results. With a uniform flow the problem is proven to be always of the Fredholm type but is found well-posed only when considering a dissipative radiation problem. With a general shear flow, the Fredholm property is recovered for a weak enough shear and the dissipative radiation problem requires to introduce extra conditions to be well-posed: enough dissipation, a large enough frequency and non-intuitive conditions on the impedance value. (10.1137/20M1384026)
  DOI : 10.1137/20M1384026
• Optimal Ciliary Locomotion of Axisymmetric Microswimmers
  
  • Guo Hanliang
  • Zhu Hai
  • Liu Ruowen
  • Bonnet Marc
  • Veerapaneni Shravan
  Journal of Fluid Mechanics, Cambridge University Press (CUP), 2021, 927, pp.A22. Many biological microswimmers locomote by periodically beating the densely-packed cilia on their cell surface in a wave-like fashion. While the swimming mechanisms of ciliated microswimmers have been extensively studied both from the analytical and the numerical point of view, the optimization of the ciliary motion of microswimmers has received limited attention, especially for non-spherical shapes. In this paper, using an envelope model for the microswimmer, we numerically optimize the ciliary motion of a ciliate with an arbitrary axisymmetric shape. The forward solutions are found using a fast boundary integral method, and the efficiency sensitivities are derived using an adjoint-based method. Our results show that a prolate microswimmer with a 2:1 aspect ratio shares similar optimal ciliary motion as the spherical microswimmer, yet the swimming efficiency can increase two-fold. More interestingly, the optimal ciliary motion of a concave microswimmer can be qualitatively different from that of the spherical microswimmer, and adding a constraint to the ciliary length is found to improve, on average, the efficiency for such swimmers. (10.1017/jfm.2021.744)
  DOI : 10.1017/jfm.2021.744