Publications

2021

• General-purpose kernel regularization of boundary integral equations via density interpolation
  
  • Maltez Faria Luiz
  • Pérez-Arancibia Carlos
  • Bonnet Marc
  Computer Methods in Applied Mechanics and Engineering, Elsevier, 2021, 378, pp.113703. This paper presents a general high-order kernel regularization technique applicable to all four integral operators of Calder\'on calculus associated with linear elliptic PDEs in two and three spatial dimensions. Like previous density interpolation methods, the proposed technique relies on interpolating the density function around the kernel singularity in terms of solutions of the underlying homogeneous PDE, so as to recast singular and nearly singular integrals in terms of bounded (or more regular) integrands. We present here a simple interpolation strategy which, unlike previous approaches, does not entail explicit computation of high-order derivatives of the density function along the surface. Furthermore, the proposed approach is kernel- and dimension-independent in the sense that the sought density interpolant is constructed as a linear combination of point-source fields, given by the same {Green's function} used in the integral equation formulation, thus making the procedure applicable, in principle, to any PDE with known {Green's function}. For the sake of definiteness, we focus here on Nystr\"om methods for the (scalar) Laplace and Helmholtz equations and the (vector) elastostatic and time-harmonic elastodynamic equations. The method's accuracy, flexibility, efficiency, and compatibility with fast solvers are demonstrated by means of a variety of large-scale three-dimensional numerical examples. (10.1016/j.cma.2021.113703)
  DOI : 10.1016/j.cma.2021.113703
• Effective wave motion in periodic discontinua near spectral singularities at finite frequencies and wavenumbers
  
  • Guzina Bojan B
  • Bonnet Marc
  Wave Motion, Elsevier, 2021, 103, pp.102729. We consider the effective wave motion, at spectral singularities such as corners of the Brillouin zone and Dirac points, in periodic continua intercepted by compliant interfaces that pertain to e.g. masonry and fractured materials. We assume the Bloch-wave form of the scalar wave equation (describing anti-plane shear waves) as a point of departure, and we seek an asymptotic expansion about a reference point in the wavenumber-frequency space-deploying wavenumber separation as the perturbation parameter. Using the concept of broken Sobolev spaces to cater for the presence of kinematic discontinuities, we next define the "mean" wave motion via inner product between the Bloch wave and an eigenfunction (at specified wavenumber and frequency) for the unit cell of periodicity. With such projection-expansion approach, we obtain an effective field equation, for an arbitrary dispersion branch, near apexes of "wavenumber quadrants" featured by the first Brillouin zone. For completeness, we investigate asymptotic configurations featuring both (a) isolated, (b) repeated, and (c) nearby eigenvalues. In the case of repeated eigenvalues, we find that the "mean" wave motion is governed by a system of wave equations and Dirac equations, whose size is given by the eigenvalue multiplicity, and whose structure is determined by the participating eigenfunctions, the affiliated cell functions, and the direction of wavenumber perturbation. One of these structures is shown to describe the so-called Dirac points-apexes of locally conical dispersion surfaces-that are relevant to the generation of topologically protected waves. In situations featuring clusters of tightly spaced eigenvalues, the effective model is found to entail a Diraclike system of equations that generates "blunted" conical dispersion surfaces. We illustrate the analysis by numerical simulations for two periodic configurations in R 2 that showcase the asymptotic developments in terms of (i) wave dispersion, (ii) forced wave motion, and (iii) frequency-and wavenumber-dependent phonon behavior. (10.1016/j.wavemoti.2021.102729)
  DOI : 10.1016/j.wavemoti.2021.102729
• On a surprising instability result of Perfectly Matched Layers for Maxwell's equations in 3D media with diagonal anisotropy
  
  • Bécache Eliane
  • Fliss Sonia
  • Kachanovska Maryna
  • Kazakova Maria
  Comptes Rendus. Mathématique, Académie des sciences (Paris), 2021. The analysis of Cartesian Perfectly Matched Layers (PMLs) in the context of time-domain electromagnetic wave propagation in a 3D unbounded anisotropic homogeneous medium modelled by a diagonal dielectric tensor is presented. Contrary to the 3D scalar wave equation or 2D Maxwell's equations some diagonal anisotropies lead to the existence of backward waves giving rise to instabilities of the PMLs. Numerical experiments confirm the presented result. (10.5802/crmath.165)
  DOI : 10.5802/crmath.165
• 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