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 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.
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 -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.
We consider highly accurate schemes for nonlinear time fractional Schrödinger equations (NTFSEs). While an strategy is employed for approximating the Caputo fractional derivative in the temporal direction, compact CCD finite difference approaches are incorporated in the space. A highly effective hybrid -CCD method is implemented successfully. The accuracy of this linearized scheme is order six in space, and order in time, where is the order of the Caputo fractional derivative involved. It...
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.
We present in this paper a stability study concerning finite volume schemes applied to the two-dimensional Maxwell system, using rectangular or triangular meshes. A stability condition is proved for the first-order upwind scheme on a rectangular mesh. Stability comparisons between the Yee scheme and the finite volume formulation are proposed. We also compare the stability domains obtained when considering the Maxwell system and the convection equation.