Special Session 144: PDEs and Irregular Interfaces: New Frontiers for Industrial Applications

A fast integral equation solver for self-similar and prefractal domains with corners

Ludvig af Klinteberg
Malardalen University (MDU)
Sweden
Co-Author(s):    Ludvig af Klinteberg
Abstract:
Boundary integral equation methods are highly efficient tools for solving elliptic boundary value problems on domains with complex geometries. Nonetheless, prefractal and self-similar domains with corners, such as the Koch curve, remain challenging because of corner singularities in the solution. In this talk, we show how state-of-the-art numerical methods can overcome these difficulties, yielding a framework for the efficient solution of discretized boundary value problems that lead to dense linear systems with millions of degrees of freedom. This makes it possible to solve the problem on very fine prefractal approximations and to study the convergence of the resulting numerical solutions in the fractal limit.

On the traces of harmonic functions $H^{1/2}$ and $H^{3/2}$ in Lipschitz domains

Cherif Amrouche
Universite de Pau et des Pays de l`Adour
France
Co-Author(s):    Mohand Moussaoui
Abstract:
In this work, we revisit the following estimate due to Dahlberg \cite{Dahl}. Let $\textit{\textbf x}_0$ a fixed point in a bounded Lipschitz domain $\Omega$. Then there exists a constant $C > 0$ such that if $u$ is a harmonic function in $\Omega$ and vanishes at $\textit{\textbf x}_0$, then \begin{equation*} C^{-1} \Vert u \Vert_{L^2(\Gamma)} \leq \Big(\int_\Omega \varrho\vert \nabla u \vert^2\Big)^{1/2} \leq C \Vert u \Vert_{L^2(\Gamma)}, \end{equation*} where $\varrho$ is the distance to the boundary of $\Omega$. Using Grisvard`s work and interpolation theory for subspaces, we complete the solvability of the inhomogeneous Dirichlet problem: $$ (\mathscr{L}_D^0)\ \ \ \ -\Delta u = f\quad \ \mbox{in}\ \Omega \quad \mbox{and } \quad u = 0 \ \ \mbox{on }\Gamma, $$ in a framework of fractional Sobolev spaces $H^s(\Omega)$, when $\Omega$ is a polygon or a polyhedron domain and $1/2 \leq s \leq 2$. Thanks to these regularity results and an explicit function given by Ne$\mathrm{\check{c}}$as, we show that the above inequalities cannot be valid in their current form. On the other hand, we identify a functional space which satisfies the embeddings $H^{1/2}_{00}(\Omega)\hookrightarrow E(\nabla;\, \Omega) \hookrightarrow H^{1/2}(\Omega)$ and the trace operator $\gamma_0$ from $E(\nabla;\, \Omega)$ into $L^2(\Gamma)$ is well-defined and continuous. This leads to an alternative to the functions $H^{1/2}(\Omega)$, non necessarily harmonic, for having a trace in $L^2(\Gamma)$ and also to a new characterization of $H^{1/2}_{00}(\Omega)$ as the kernel of this operator. However, we show that if the domain $\Omega$ is of class $\mathscr{C}^{1, 1}$, then the above inequalities are valid.

A Functional Framework for Maxwell`s Equations on $H^1$-Extension Domains

Romain CERVERA
Universite Paris Saclay
France
Co-Author(s):    Anna Rozanova-Pierrat, Patrick Ciarlet, Alexander Teplyaev
Abstract:
The mathematical analysis of boundary value problems for Maxwell's equations traditionally requires at least a Lipschitz-continuous boundary to properly define trace operators. Building upon the work of Creo, Lancia, Vernole, Hinz, and Teplyaev (2018) regarding fractal boundaries, we address the challenge of formulating electromagnetic problems on highly irregular geometries. Our focus is on the large class of $H^1$-extension domains, where classical tools relying on surface measures and standard boundary regularity break down. To overcome this, we propose a boundary measure-free approach that relies entirely on the intrinsic Hilbert structure of the trace space. Within this framework, we construct the notion of normal and tangential trace operators for the spaces $\mathbf{H}(\operatorname{div}, \Omega)$ and $\mathbf{H}(\operatorname{\mathbf{curl}}, \Omega)$, providing a rigorous basis for solving electromagnetic boundary value problems in geometries where classical theory is not applicable.

Acoustic impedance scattering on extension domains

Gabriel CLARET
Paris-Saclay University
France
Co-Author(s):    Anna Rozanova-Pierrat
Abstract:
From sonars to medical imaging or even robot vacuum cleaners, acoustic scattering is at the core of numerous applications. In an acoustic scattering problem, an incident acoustic field is scattered by an obstacle in such a way that the total field -- superposition of the incident and scattered ones -- satisfies a given boundary condition. Modelling that problem, a crucial aspect is the geometry of the obstacle, and it is known that irregular shapes such as fractals are most relevant to represent the behaviour of real-life objects. In this talk, we discuss the acoustic scattering problem by an obstacle described by a Sobolev extension domain. Those domains can be smooth, Lipschitz, but also fractal, multi-fractal... That obstacle is endowed with an impedance-type boundary condition, understood with but also without a boundary measure. Relying on tools from potential theory, we study the well-posedness of the impedance scattering problem in both cases, which we restate in terms of boundary equations. As an application of our analysis, we consider the inverse scattering problem consisting in identifying the obstacle and the boundary condition from many scattered plane waves.

Obstacle problems for nonlinear fractional operators in irregular domains

Simone Creo
Sapienza Universita' di Roma
Italy
Co-Author(s):    Salvatore Fragapane
Abstract:
In this talk we study obstacle problems for the regional fractional $p$-Laplacian in a domain $\Omega\subset\mathbb{R}^2$ having fractal boundary of Koch type. We first prove well-posedness results for the solution of the obstacle problem, as well as two equivalent formulations. Then, we study corresponding approximating obstacle problems in a sequence of domains $\Omega_n\subset\mathbb{R}^2$ having pre-fractal boundary, for $n\in\mathbb{N}$. After proving the well-posedness of the approximating obstacle problems, we perform the asymptotic analysis for both $n\to+\infty$ and $p\to+\infty$.

Integral equation methods for ice floes and other surface wave scattering problems

Tristan Goodwill
University of Chicago
USA
Co-Author(s):    
Abstract:
Ocean waves cause large deformations floating ice sheets and ice floes, which can cause icequakes and more generally impact the stability of the cryosphere. These waves are also an example of bulk-surface problems involving complicated boundary conditions, nonlocal effects, and infinite or extended interfaces. In this talk, I present novel integral equation methods and fast solvers for a broad class of surface wave problems. The integral equation is Fredholm second kind and only posed on a compact region of interest. I demonstrate the effectiveness of these methods on a variety of glaciological applications.

A discontinuous Galerkin method on a fractal domain

David P Hewett
University College London
England
Co-Author(s):    Andrea Moiola, Sergio Gomez
Abstract:
When using classical FEM to solve a PDE problem in a domain with a complicated or non-smooth boundary, one typically requires a mesh with a large number of elements in order to accurately capture the boundary geometry. In this talk we present a discontinuous Galerkin (dG) FEM for the solution of PDE problems in a domain with a fractal boundary (the Dirichlet Poisson problem in the Koch snowflake) in which the boundary geometry is captured exactly, by using a mesh comprising elements which themselves have fractal boundary. On such meshes the classical dG normal derivative interface terms cannot be defined in the usual way because inter-element interfaces may not have a well-defined normal vector. We show how this can be addressed by introducing proxies for the normal derivative terms, based on integrals over subsets of the mesh elements. We prove well-posedness of the resulting dG FEM and provide a partial error analysis. We also discuss practical implementation of the method and provide numerical results demonstrating its effectiveness. This is joint work with Sergio Gomez and Andrea Moiola.

On the Quadrature Error Estimates with PINN-Based Densities in Boundary Integral Methods

Ismail Labaali
Sapienza university
Italy
Co-Author(s):    
Abstract:
We derive quadrature error estimates for the evaluation of boundary integral potentials with logarithmic singularities, considering standard discretization techniques such as Gauss--Legendre and trapezoidal rules on general geometries. The analysis provides explicit asymptotic expressions for the quadrature error, highlighting its dependence on the distance of the evaluation point from the boundary and on the number of quadrature points. These estimates offer a practical tool for predicting accuracy in nearly singular regimes and for guiding adaptive strategies. As a complementary study, we investigate the applicability of these error estimates when the density is not known analytically but is approximated via physics-informed neural networks (PINNs). The learned density is used to reconstruct the potential through standard quadrature, and the resulting error is compared with the theoretical predictions. We examine whether the derived estimates remain valid in this data-driven setting and identify conditions under which they accurately capture the observed behavior. In particular, we show that the consistency of the estimates depends on the smoothness and local accuracy of the PINN density, especially near the boundary

Robin Green Function Estimates and a Model of Mammalian Lungs

Marco Michetti
Gothenburg University
Sweden
Co-Author(s):    Guy David, Stefano Decio, Max Engelstein, Marcel Filoche and Svitlana Mayboroda
Abstract:
In this talk we present delicate properties of the Green function with Robin boundary conditions, in particular, elucidating the nature of the passage between the Dirichlet-like and Neumann-like behavior. This yields sharp quantifiable bounds on the corresponding harmonic measure and proves the phase transition in the behavior of the total flow earlier conjectured in physics literature in concert with the efficacy of mammalian lungs.

High order cubature for iterated function system

Zo{\\"\\i}s Moitier
ENSTA
France
Co-Author(s):    Patrick Joly, Maryna Kachanovska
Abstract:
The talk is motivated by application to fractal antenna engineering, where antennas with self-similar shapes operate across multiple frequencies. Recent works by Chandler-Wilde, Gibbs, Hewett, Moiola and their co-workers propose boundary integral formulations for solving Helmholtz scattering problems on fractal screens with Dirichlet boundary conditions. These boundary integral equations utilize the Hausdorff measure on fractals instead of the standard Lebesgue`s measure. This makes the evaluation of the boundary integral difficult from a numerical viewpoint. To tackle this point, we propose new interpolatory high-order tensor cubature formula on fractals, based on Chebyshev points on an interval. These formulas allow computing integrals of restrictions of regular functions to fractals with a high accuracy. We will discuss the construction of such cubature (in particular computation of the weights) and their properties (asymptotic behavior of weights).

Potential flow solver for ship hydrodynamics problems with fully nonlinear free surface boundary conditions

Andrea Mola
IMT School for Advanced Studies Lucca
Italy
Co-Author(s):    Andrea Mola
Abstract:
This contribution will present a mathematical model for the simulation of the flow past a ship sailing in calm or wavy water. The fluid dynamic model is based on the potential flow theory and is governed by the Laplace equation, here complemented by fully nonlinear boundary conditions on the water free surface and non penetration boundary conditions on the ship hull. The spatial discretization of the Laplace equation is carried out by means of the Boundary Element Method (BEM), while the discretization of the time dependent fully nonlinear free surface boundary conditions is based on the Finite Element Method (FEM). Such a combined FEM-BEM spatial discretization strategy results in a Differential Algebraic Equations (DAE) system, that is time integrated via implicit Backward Difference Formula (BDF) method with variable degree and time step. The numerical solver is implemented in a stand alone C++ program which efficiently leverages on several open source numerical libraries. Validation test cases based on comparison with experimental measurements obtained on industrial benchmark hulls will be presented. In addition, some results on the application on competition rowing boat hulls will be shown and discussed.

Identification of constant coefficients in a model of linear anisotropic subdiffusion

Gianluca Mola
Sorbonne University Abu Dhabi
Italy
Co-Author(s):    Simone Creo, Maria Rosaria Lancia, Andrea Mola and Silvia Romanelli
Abstract:
Let $\left(H, \langle \cdot, \cdot \rangle \right)$ be a Hilbert space and $A_{i}:D(A_i) \to H$ ($i = 1,\cdots,n$) be a family of nonnegative and self-adjoint operators mutually commuting. We study the inverse problem consisting in the identification of the function $u:[0,\infty) \to H$ and $n$ constants $\lambda_{1},\cdots,\lambda_{n} > 0$ that fulfill the initial-value problem $$ \partial_{t}^{\alpha}u(t) + \lambda_{1} A_{1}u(t) + \cdots + \lambda_{n} A_{n}u(t) = 0, \quad t > 0, \quad u(0) = x, $$ and the additional energy measurements at $t=T$ $$ \left\langle A_{1} u(T),u(T)\right\rangle = \varphi_{1}, \quad \cdots \quad, \left\langle A_{n} u(T),u(T)\right\rangle = \varphi_{n}. $$ Here, $\partial_{t}^{\alpha}$ denotes the \emph{Caputo fractional derivative} of order $0 < \alpha < 1$, which accounts for models of subdiffusion phenomena. Under suitable assumptions on the initial datum $x \in H$, we shall provide a uniqueness result. Existence is also achieved, assuming that the overdeterminating conditions $\varphi_{1},\cdots,\varphi_{n} > 0$ belong to a proper set, which is described in concrete cases by means of numerical simulations for different values of $\alpha$. Finally, we prove the convergence of the solution $(u,\lambda_{1},\cdots,\lambda_{n})$ as $\alpha \to 1^{-}$, thus linking the fractional derivative model with the classic derivative one. Applications to heat conduction and plane elasticity are considered.

On Rumples and Rolls: Efficient Representations for Elastodynamics in Ice Sheets

Peter Nekrasov
University of Chicago
USA
Co-Author(s):    Peter Nekrasov, Jeremy G. Hoskins
Abstract:
Elastic waves are involved in a number of important phenomena, including the bending and flexing of floating ice sheets driven by the ocean. These elastic effects are frequently modeled using a time-harmonic wave equation with zero traction boundary conditions, which correspond to a free or unconstrained surface. In this talk, we introduce a novel approach for handling coupled hydro-elastic waves that reduce the problem to a Fredholm second kind integral equation on the boundary. These techniques are applied to ongoing fieldwork around the rumple zones of the McMurdo Ice Shelf in order to better understand wave localization and propagation in the polar regions.

The Drude-Born-Fedorov system on anisotropic fractal porous media

Eric Stachura
Kennesaw State University
USA
Co-Author(s):    Ioannis Stratis
Abstract:
In this talk, I will discuss a mathematical framework for wave propagation in anisotropic fractal porous media. These materials are fractal structures that could have different dimensions in each direction and appear, for instance, in the modeling of composite structures with fractal type microstructures. I will discuss the vector framework to study electromagnetic wave propagation and apply it to the special case of the Drude-Born-Fedorov system--a particular model of an electromagnetically bi-anisotropic medium which incorporates chirality. This is joint work with I. G. Stratis (National and Kapodistrian University of Athens).

The radiation condition for Helmholtz equations above (locally perturbed) periodic surfaces

Ruming Zhang
Technical University of Berlin
Germany
Co-Author(s):    Ruming Zhang
Abstract:
The radiation condition is the key question in the mathematical modelling for scattering problems in unbounded domains. Mathematically, it plays the role as the boundary condition at the infinity, which guarantees the well-posedness of the mathematical problem; physically, it describes the far-field asymptotic behaviour of the physical waves. In this paper, we focus on the radiation conditions for scattering problems above (locally perturbed) periodic surfaces. According to Hu et al. (2021), the radiating solution satisfies the Sommerfeld radiation condition: \[ \frac{\partial u}{\partial r}-i k u=o\left(r^{-1/2}\right) \] Although there are literature which have studied this problem, there is no specific method for dealing with periodic structures. Due to this reason, the important properties for the periodic structures may be ignored. Moreover, the existing method is not extendable to bi-periodic structures in three dimensional spaces. In this talk, we study the radiation condition for the time-harmonic scattering problem with periodic surfaces, which is modelled by the Helmholtz equation. We introduce a novel method based on the Floquet-Bloch transform, which, to the best of the author`s knowledge, is the first method that works particularly for periodic media. With this method, we improve the Sommerfeld radiation condition for the scattered field from periodic media to: \[ \frac{\partial u}{\partial r}-i k u=O\left(r^{-3/2}\right) \] More importantly, the prospect of extending this method to 3D cases is optimistic.