Systems of parabolic differential equations are studied in the paper. Two a posteriori error estimates for the approximate solution obtained by the finite element method of lines are presented. A statement on the rate of convergence of the approximation of error by estimator to the error is proved.

This paper is devoted to the study of a posteriori error estimates for the scalar nonlinear convection-diffusion-reaction equation $c_t+\nabla ·\left(𝐮f\left(c\right)\right)-\nabla ·\left(D\nabla c\right)+\lambda c=0$. The estimates for the error between the exact solution and an upwind finite volume approximation to the solution are derived in the $L^1$-norm, independent of the diffusion parameter $D$. The resulting a posteriori error estimate is used to define an grid adaptive solution algorithm for the finite volume scheme. Finally numerical experiments underline the applicability...

This paper is devoted to the study of a posteriori error estimates for the scalar nonlinear convection-diffusion-reaction equation $c_t+\nabla ·\left(𝐮f\left(c\right)\right)-\nabla ·\left(D\nabla c\right)+\lambda c=0$. The estimates for the error between the exact solution and an upwind finite volume approximation to the solution are derived in the L1-norm, independent of the diffusion parameter D. The resulting a posteriori error estimate is used to define an grid adaptive solution algorithm for the finite volume scheme. Finally numerical experiments underline the applicability...

We deal with a posteriori error estimates of the discontinuous Galerkin method applied to the nonstationary heat conduction equation. The problem is discretized in time by the backward Euler scheme and a posteriori error analysis is based on the Helmholtz decomposition.

We will investigate the possibility to use superconvergence results for the mixed finite element discretizations of some time-dependent partial differential equations in the construction of a posteriori error estimators. Since essentially the same approach can be followed in two space dimensions, we will, for simplicity, consider a model problem in one space dimension.

We consider the efficient and reliable solution of linear-quadratic optimal control problems governed by parametrized parabolic partial differential equations. To this end, we employ the reduced basis method as a low-dimensional surrogate model to solve the optimal control problem and develop a posteriori error estimation procedures that provide rigorous bounds for the error in the optimal control and the associated cost functional. We show that our approach can be applied to problems involving...

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...

The goal of this paper is to construct a first-order upwind scheme for solving the system of partial differential equations governing the one-dimensional flow of two superposed immiscible layers of shallow water fluids. This is done by generalizing a numerical scheme presented by Bermúdez and Vázquez-Cendón [3, 26, 27] for solving one-layer shallow water equations, consisting in a $Q$-scheme with a suitable treatment of the source terms. The difficulty in the two layer system comes from the coupling...

The goal of this paper is to construct a first-order upwind scheme for solving the system of partial differential equations governing the one-dimensional flow of two superposed immiscible layers of shallow water fluids. This is done by generalizing a numerical scheme presented by Bermúdez and Vázquez-Cendón [3, 6, 27] for solving one-layer shallow water equations, consisting in a Q-scheme with a suitable treatment of the source terms. The difficulty in the two layer system comes from the coupling...

We propose a quasi-Newton algorithm for solving fluid-structure interaction problems. The basic idea of the method is to build an approximate tangent operator which is cost effective and which takes into account the so-called added mass effect. Various test cases show that the method allows a significant reduction of the computational effort compared to relaxed fixed point algorithms. We present 2D and 3D fluid-structure simulations performed either with a simple 1D structure model or with shells...

We propose a quasi-Newton algorithm for solving fluid-structure interaction problems. The basic idea of the method is to build an approximate tangent operator which is cost effective and which takes into account the so-called added mass effect. Various test cases show that the method allows a significant reduction of the computational effort compared to relaxed fixed point algorithms. We present 2D and 3D fluid-structure simulations performed either with a simple 1D structure model or with...

Existence of a solution to the quasi-variational inequality problem arising in a model for sand surface evolution has been an open problem for a long time. Another long-standing open problem concerns determining the dual variable, the flux of sand pouring down the evolving sand surface, which is also of practical interest in a variety of applications of this model. Previously, these problems were solved for the special case in which the inequality is simply variational. Here, we introduce a regularized...

It is shown that the approximating equations whose existence is required in the author's previous work on partially regular weak solutions can be constructed without any additional assumption about the equation itself. This leads to a variation of a Galerkin method.