This work concerns the numerical approximation with an explicit first-order finite volume method of inviscid, nonequilibrium, high-temperature flows in multiple space dimensions. It is devoted to the analysis of the numerical scheme for the approximation of the hyperbolic system in homogeneous form. We derive a general framework for the design of numerical schemes for this model from numerical schemes for the monocomponent compressible Euler equations for a polytropic gas. Under a very simple condition on the adiabatic exponent of the polytropic gas, the scheme for the multicomponent system enjoys the same properties as the one for the monocomponent system: discrete entropy inequality, positivity of the partial densities and internal energies, discrete maximum principle on the mass fractions, and discrete minimum principle on the entropy. Our approach extends the relaxation of energy (Coquel and Perthame in SIAM J. Numer. Anal. 35:2223–2249, 1998) to the multicomponent Euler system. In the limit of instantaneous relaxation we show that the solution formally converges to a unique and stable equilibrium solution to the multicomponent Euler equations. We then use this framework to design numerical schemes from three schemes for the polytropic Euler system: the Godunov exact Riemann solver, and the HLL and pressure relaxation based approximate Riemann solvers. Numerical experiments in one and two space dimensions on flows with discontinuous solutions support the conclusions of our analysis and highlight stability, robustness and convergence of the schemes.

The scattering of electromagnetic waves from obstacles with wave-material interaction in thin layers on the surface is described by generalized impedance boundary conditions, which provide effective approximate models. In particular, this includes a thin coating around a perfect conductor and the skin effect of a highly conducting material. The approach taken in this work is to derive, analyse and discretize a system of time-dependent boundary integral equations that determines the tangential traces of the scattered electric and magnetic fields. In a familiar second step, the fields are evaluated in the exterior domain by a representation formula, which uses the time-dependent potential operators of Maxwell’s equations. The time-dependent boundary integral equation is discretized with Runge–Kutta based convolution quadrature in time and Raviart–Thomas boundary elements in space. Using the frequency-explicit bounds from the well-posedness analysis given here together with known approximation properties of the numerical methods, the full discretization is proved to be stable and convergent, with explicitly given rates in the case of sufficient regularity. Taking the same Runge–Kutta based convolution quadrature for discretizing the time-dependent representation formulas, the optimal order of convergence is obtained away from the scattering boundary, whereas an order reduction occurs close to the boundary. The theoretical results are illustrated by numerical experiments.

We consider the numerical approximation of the ill-posed data assimilation problem for stationary convection–diffusion equations and extend our previous analysis in Burman et al. (Numer. Math. 144:451–477, 2020) to the convection-dominated regime. Slightly adjusting the stabilized finite element method proposed for dominant diffusion, we draw upon a local error analysis to obtain quasi-optimal convergence along the characteristics of the convective field through the data set. The weight function multiplying the discrete solution is taken to be Lipschitz continuous and a corresponding super approximation result (discrete commutator property) is proven. The effect of data perturbations is included in the analysis and we conclude the paper with some numerical experiments.

We study the performance of some preconditioning techniques for a class of block three-by-three linear systems of equations arising from finite element discretizations of the coupled Stokes–Darcy flow problem. In particular, we investigate preconditioning techniques including block preconditioners, constraint preconditioners, and augmented Lagrangian-based ones. Spectral and field-of-value analyses are established for the exact versions of these preconditioners. The result of numerical experiments are reported to illustrate the performance of inexact variants of the various preconditioners used with flexible GMRES in the solution of a 3D test problem with large jumps in the permeability.

The present work proposes a second-order time splitting scheme for a linear dispersive equation with a variable advection coefficient subject to transparent boundary conditions. For its spatial discretization, a dual Petrov–Galerkin method is considered which gives spectral accuracy. The main difficulty in constructing a second-order splitting scheme in such a situation lies in the compatibility condition at the boundaries of the sub-problems. In particular, the presence of an inflow boundary condition in the advection part results in order reduction. To overcome this issue a modified Strang splitting scheme is introduced that retains second-order accuracy. For this numerical scheme a stability analysis is conducted. In addition, numerical results are shown to support the theoretical derivations.

The present paper is devoted to the convergence analysis of a class of asymptotic preserving particle schemes (Filbet and Rodrigues in SIAM J. Numer. Anal. 54(2):1120–1146, 2016) for the Vlasov equation with a strong external magnetic field. In this regime, classical Particle-in-Cell methods are subject to quite restrictive stability constraints on the time and space steps, due to the small Larmor radius and plasma frequency. The asymptotic preserving discretization that we are going to study removes such a constraint while capturing the large-scale dynamics, even when the discretization (in time and space) is too coarse to capture fastest scales. Our error bounds are explicit regarding the discretization, stiffness parameter, initial data and time.

We present variational approximations of boundary value problems for curvature flow (curve shortening flow) and elastic flow (curve straightening flow) in two-dimensional Riemannian manifolds that are conformally flat. For the evolving open curves we propose natural boundary conditions that respect the appropriate gradient flow structure. Based on suitable weak formulations we introduce finite element approximations using piecewise linear elements. For some of the schemes a stability result can be shown. The derived schemes can be employed in very different contexts. For example, we apply the schemes to the Angenent metric in order to numerically compute rotationally symmetric self-shrinkers for the mean curvature flow. Furthermore, we utilise the schemes to compute geodesics that are relevant for optimal interface profiles in multi-component phase field models.

In this paper we consider a finite-volume approximation for the Cahn–Hilliard equation with dynamic boundary conditions. The convergence of the scheme is proved in Nabet (IMA J Numer Anal 36(4): 1898–1942, 2016), we prove here an error estimate for the fully-discrete scheme. We also give numerical simulations which validate the theoretical result.

We design and analyze a C^0 interior penalty method for the approximation of classical solutions of the Dirichlet boundary value problem of the Monge–Ampère equation on convex polygonal domains. The method is based on an enhanced cubic Lagrange finite element that enables the enforcement of the convexity of the approximate solutions. Numerical results that corroborate the a priori and a posteriori error estimates are presented. It is also observed from numerical experiments that this method can capture certain weak solutions.