We derive a posteriori estimates for a discretization in space of the standard Cahn–Hilliard equation with a double obstacle free energy. The derived estimates are robust and efficient, and in practice are combined with a heuristic time step adaptation. We present numerical experiments in two and three space dimensions and compare our method with an existing heuristic spatial mesh adaptation algorithm.
We deal with the numerical solution of the nonstationary heat conduction equation with mixed Dirichlet/Neumann boundary conditions. The backward Euler method is employed for the time discretization and the interior penalty discontinuous Galerkin method for the space discretization. Assuming shape regularity, local quasi-uniformity, and transition conditions, we derive both a posteriori upper and lower error bounds. The analysis is based on the Helmholtz decomposition, the averaging interpolation...
We deal with a nonstationary semilinear singularly perturbed convection–diffusion problem. We discretize this problem by discontinuous Galerkin method in space and by midpoint rule in time. We present diffusion–uniform error estimates with sketches of proofs.
In this paper we consider a model shape optimization problem. The state variable solves an elliptic equation on a domain with one part of the boundary described as the graph of a control function. We prove higher regularity of the control and develop a priori error analysis for the finite element discretization of the shape optimization problem under consideration. The derived a priori error estimates are illustrated on two numerical examples.
Mathematical models for option pricing often result in partial differential equations. Recent enhancements are models driven by Lévy processes, which lead to a partial differential equation with an additional integral term. In the context of model calibration, these partial integro differential equations need to be solved quite frequently. To reduce the computational cost the implementation of a reduced order model has shown to be very successful numerically. In this paper we give a priori error...
In this paper we solve the time-dependent incompressible Navier-Stokes equations by splitting the non-linearity and incompressibility, and using discontinuous or continuous finite element methods in space. We prove optimal error estimates for the velocity and suboptimal estimates for the pressure. We present some numerical experiments.
In this paper we solve the time-dependent incompressible Navier-Stokes equations by splitting the non-linearity and incompressibility, and using discontinuous or continuous finite element methods in space. We prove optimal error estimates for the velocity and suboptimal estimates for the pressure. We present some numerical experiments.
The paper investigates the Galerkin method for an initial boundary value problem for heat convection equations. New error estimates for the approximate solutions and their derivatives in strong norm are obtained.
We address the numerical approximation of the two-phase Stefan problem and discuss an adaptive finite element method based on rigorous a posteriori error estimation and refinement/coarsening. We also investigate how to restrict coarsening for the resulting method to be stable and convergent. We review implementation issues associated with bisection and conclude with simulations of a persistent corner singularity, for which adaptivity is an essential tool.
This work represents a first step towards the simulation of the motion of water in a complex hydrodynamic configuration, such as a channel network or a river delta, by means of a suitable “combination” of different mathematical models. In this framework a wide spectrum of space and time scales is involved due to the presence of physical phenomena of different nature. Ideally, moving from a hierarchy of hydrodynamic models, one should solve throughout the whole domain the most complex model (with...
In this paper, a dynamic viscoelastic problem is numerically studied. The variational problem is written in terms of the velocity field and it leads to a parabolic linear variational equation. A fully discrete scheme is introduced by using the finite element method to approximate the spatial variable and an Euler scheme to discretize time derivatives. An a priori error estimates result is recalled, from which the linear convergence is derived under suitable regularity conditions. Then, an a posteriori...
In this paper, a dynamic viscoelastic problem is numerically studied. The variational problem is written in terms of the velocity field and it leads to a parabolic linear variational equation. A fully discrete scheme is introduced by using the finite element method to approximate the spatial variable and an Euler scheme to discretize time derivatives. An a priori error estimates result is recalled, from which the linear convergence is derived under suitable regularity conditions. Then, an a posteriori error...
In this article we introduce an adaptive multi-level method in space and time for convection diffusion problems. The scheme is based on a multi-level spatial splitting and the use of different time-steps. The temporal discretization relies on the characteristics method. We derive an a posteriori error estimate and design a corresponding adaptive algorithm. The efficiency of the multi-level method is illustrated by numerical experiments, in particular for a convection-dominated problem.
We present the current Reduced Basis framework for the efficient numerical approximation of parametrized steady Navier–Stokes equations. We have extended the existing setting developed in the last decade (see e.g. [S. Deparis, SIAM J. Numer. Anal. 46 (2008) 2039–2067; A. Quarteroni and G. Rozza, Numer. Methods Partial Differ. Equ. 23 (2007) 923–948; K. Veroy and A.T. Patera, Int. J. Numer. Methods Fluids 47 (2005) 773–788]) to more general affine and nonaffine parametrizations (such as volume-based...
We study error estimates and their convergence rates for approximate solutions of spectral Galerkin type for the equations for the motion of a viscous chemical active fluid in a bounded domain. We find error estimates that are uniform in time and also optimal in the L2-norm and H1-norm. New estimates in the H(-1)-norm are given.
We describe a numerical method for the equation in with Dirichlet boundary and initial conditions which is a combination of the method of characteristics and the finite-difference method. We prove both an a priori local error-estimate of a high order and stability. Example 3.3 indicates that our approximate solutions are disturbed only by a minimal amount of the artificial diffusion.