Related papers: Towards modular Hierarchical Poincar\'{e}-Steklov …
We consider the design of structure-preserving discretization methods for the solution of systems of boundary controlled Partial Differential Equations (PDEs) thanks to the port-Hamiltonian formalism. We first provide a novel general…
We present a first order system least squares (FOSLS) method for the Helmholtz equation at high wave number k, which always deduces Hermitian positive definite algebraic system. By utilizing a non-trivial solution decomposition to the dual…
Thanks to modern manufacturing technologies, heterogeneous materials with complex inner structures (e.g., foams) can be easily produced. However, their utilization is not straightforward, as the classical constitutive laws are not…
In the presence of the supersymmetric CP phases, e^+ e^- \to Z^{*} \to h A scattering is analyzed, with special emphasis on Z^{*} h A vertex, taking into account radiatve corrections due to dominant top quark and top squark loops in the…
We present a new framework for expressing finite element methods on multiple intersecting meshes: multimesh finite element methods. The framework enables the use of separate meshes to discretize parts of a computational domain that are…
As the use of spectral/$hp$ element methods, and high-order finite element methods in general, continues to spread, community efforts to create efficient, optimized algorithms associated with fundamental high-order operations have grown.…
This paper presents a convergence analysis of the evolving surface finite element method (ESFEM) applied to the original Eyles-King-Styles model of tumour growth. The model consists of a Poisson equation in the bulk, a forced mean curvature…
Cluster perturbation theory in combination with the Lanczos method is used to compute the one-electron spectral function of the Holstein polaron in one and two dimensions. It is shown that the method allows reliable calculations using…
In this paper, we present a finite element method (FEM) framework enhanced by an operator-adapted wavelet decomposition algorithm designed for the efficient analysis of multiscale electromagnetic problems. Usual adaptive FEM approaches,…
One approach for the simulation of metamaterials is to extend an associated continuum theory concerning its kinematic equations, and the relaxed micromorphic continuum represents such a model. It incorporates the Curl of the nonsymmetric…
We discuss a system of third order PDEs for strictly convex smooth functions on domains of Euclidean space. We argue that it may be understood as a closure of sorts of the first order prolongation of a family of second order PDEs. We…
The Poisson-Boltzmann equation is widely used to model molecular electrostatics; however, it is usually solved in linearised form because the sinh nonlinearity is challenging, limiting its applicability in highly charged systems such as…
We compute small rational models for configuration spaces of points on oriented surfaces, as right modules over the framed little disks operad. We do this by splitting these surfaces in unions of several handles. We first describe rational…
The $h$-version of the finite-element method ($h$-FEM) applied to the high-frequency Helmholtz equation has been a classic topic in numerical analysis since the 1990s. It is now rigorously understood that (using piecewise polynomials of…
In the first part of the paper we introduce some geometric tools needed to describe slow-fast Hamiltonian systems on smooth manifolds. We start with a smooth Poisson bundle $p: M\to B$ of a regular (i.e. of constant rank) Poisson manifold…
This paper extends the Method of Particular Solutions (MPS) to the computation of eigenfrequencies and eigenmodes of plates. Specific approximation schemes are developed, with plane waves (MPS-PW) or Fourier-Bessel functions (MPS-FB). This…
We combine the parameterization method for invariant manifolds with the finite element method for elliptic PDEs,to obtain a new computational framework for high order approximation of invariant manifolds attached to unstable equilibrium…
Using the point fusion procedure we obtain the new integrable systems from the Elliptic Schlesinger system (ESS). These new systems have the pole orders higher than one in the matrix of the Lax operator. Quadratic Poisson algebras on the…
Numerical integration of the stiffness matrix in higher order finite element (FE) methods is recognized as one of the heaviest computational tasks in a FE solver. The problem becomes even more relevant when computing the Gram matrix in the…
The nonlocal Cahn-Hilliard equation provides a natural extension of the classical model for phase separation by incorporating long-range interactions through a singular convolution kernel. While this formulation admits a rich existence and…