Andrea Moiola

Andrea Moiola, Nicola Galante, Emile Parolin

Classical Trefftz methods approximate Helmholtz solutions using propagative plane waves and are subject to strong numerical instabilities. Evanescent plane wave bases can substantially mitigate this phenomenon. We propose a simple recipe to select such basis functions. We show that the numerical results obtained by the Ultraweak Variational Formulation (UWVF) greatly improve thanks to this choice. More details and examples will soon be available in [Galante, Moiola, Parolin 2026].

Sergio Gómez, David Hewett, Andrea Moiola

We formulate, analyse, and implement a discontinuous Galerkin finite element method (DG-FEM) for the approximation of the solution of an elliptic boundary value problem in a domain with fractal boundary. We consider the case of the Poisson equation in the Koch snowflake domain with zero Dirichlet boundary conditions, but our methodology can be generalised to other cases. Rather than first approximating the snowflake domain by a polygonal "prefractal'' and then applying a standard DG-FEM on the prefractal, we define a DG-FEM on the snowflake itself, using a geometry-conforming mesh (a fractal tiling) consisting of fractal elements, each similar to the original snowflake. Fluxes across inter-element boundaries, which are fractal curves, are represented in a weak way by integrals over element subdomains. We show how, for local polynomial basis functions, these integrals can be evaluated exactly using the similarity of the elements. We prove well-posedness and quasi-optimality of the method, and provide a partial convergence analysis. We present numerical results for piecewise linear and piecewise quadratic basis functions, which demonstrate the effectiveness of the method.

Chiara Perinati, Lise-Marie Imbert-Gérard, Andrea Moiola et al.

We present and analyze a discontinuous Galerkin method for the numerical solution of a class of second-order linear mixed-type partial differential equations, i.e. equations that change their nature from elliptic to hyperbolic through the computational domain. Well-posedness of the discrete problem is established via coercivity in an energy norm, achieved through the Morawetz multiplier technique. We derive $hp$-a priori error estimates in the energy norm, which we use to prove convergence rates for standard and quasi-Trefftz polynomial spaces. Numerical experiments validate the theoretical results.