Publications

2022

• The Complex-Scaled Half-Space Matching Method
  
  • Bonnet-Ben Dhia Anne-Sophie
  • Chandler-Wilde Simon N.
  • Fliss Sonia
  • Hazard Christophe
  • Perfekt Karl-Mikael
  • Tjandrawidjaja Yohanes
  SIAM Journal on Mathematical Analysis, Society for Industrial and Applied Mathematics, 2022, 54 (1), pp.512-557. The Half-Space Matching (HSM) method has recently been developed as a new method for the solution of 2D scattering problems with complex backgrounds, providing an alternative to Perfectly Matched Layers (PML) or other artificial boundary conditions. Based on half-plane representations for the solution, the scattering problem is rewritten as a system of integral equations in which the unknowns are restrictions of the solution to the boundaries of a finite number of overlapping half-planes contained in the domain: this integral equation system is coupled to a standard finite element discretisation localised around the scatterer. While satisfactory numerical results have been obtained for real wavenumbers, wellposedness and equivalence to the original scattering problem have been established only for complex wavenumbers. In the present paper, by combining the HSM framework with a complex-scaling technique, we provide a new formulation for real wavenumbers which is provably well-posed and has the attraction for computation that the complex-scaled solutions of the integral equation system decay exponentially at infinity. The analysis requires the study of double-layer potential integral operators on intersecting infinite lines, and their analytic continuations. The effectiveness of the method is validated by preliminary numerical results. (10.1137/20M1387122)
  DOI : 10.1137/20M1387122
• Local transparent boundary conditions for wave propagation in fractal trees (ii): error and complexity analysis
  
  • Joly Patrick
  • Kachanovska Maryna
  SIAM Journal on Numerical Analysis, Society for Industrial and Applied Mathematics, 2022, 60 (2). This work is dedicated to a refined error analysis of the high-order transparent boundary conditions introduced in the companion work [8] for the weighted wave equation on a fractal tree. The construction of such boundary conditions relies on truncating the meromorphic series that represents the symbol of the Dirichlet-to-Neumann operator. The error induced by the truncation depends on the behaviour of the eigenvalues and the eigenfunctions of the weighted Laplacian on a self-similar metric tree. In this work we quantify this error by computing asymptotics of the eigenvalues and bounds for Neumann traces of the eigenfunctions. We prove the sharpness of the obtained bounds for a class of self-similar trees. (10.1137/20M1357524)
  DOI : 10.1137/20M1357524
• The Morozov's principle applied to data assimilation problems
  
  • Bourgeois Laurent
  • Dardé Jérémi
  ESAIM: Mathematical Modelling and Numerical Analysis, Société de Mathématiques Appliquées et Industrielles (SMAI) / EDP, 2022. This paper is focused on the Morozov's principle applied to an abstract data assimilation framework, with particular attention to three simple examples: the data assimilation problem for the Laplace equation, the Cauchy problem for the Laplace equation and the data assimilation problem for the heat equation. Those ill-posed problems are regularized with the help of a mixed type formulation which is proved to be equivalent to a Tikhonov regularization applied to a well-chosen operator. The main issue is that such operator may not have a dense range, which makes it necessary to extend well-known results related to the Morozov's choice of the regularization parameter to that unusual situation. The solution which satisfies the Morozov's principle is computed with the help of the duality in optimization, possibly by forcing the solution to satisfy given a priori constraints. Some numerical results in two dimensions are proposed in the case of the data assimilation problem for the Laplace equation. (10.1051/m2an/2022061)
  DOI : 10.1051/m2an/2022061
• On the approximation of electromagnetic fields by edge finite elements. Part 4: analysis of the model with one sign-changing coefficient
  
  • Ciarlet Patrick
  Numerische Mathematik, Springer Verlag, 2022, 152, pp.223-257. In electromagnetism, in the presence of a negative material surrounded by a classical material, the electric permittivity, and possibly the magnetic permeability, can exhibit a sign-change at the interface. In this setting, the study of electromagnetic phenomena is a challenging topic. We focus on the time-harmonic Maxwell equations in a bounded set $\Omega$ of ${\mathbb R}^3$, and more precisely on the numerical approximation of the electromagnetic fields by edge finite elements. Special attention is paid to low-regularity solutions, in terms of the Sobolev scale $({\boldsymbol{H}}^{\mathtt{s}}(\Omega))_{\mathtt{s}>0}$. With the help of T-coercivity, we address the case of one sign-changing coefficient, both for the model itself, and for its discrete version. Optimal a priori error estimates are derived. (10.1007/s00211-022-01315-x)
  DOI : 10.1007/s00211-022-01315-x