On Preservation of Positivity in Some Finite Element Methods for the Heat Equation

2015 ◽  
Vol 15 (4) ◽  
pp. 417-437 ◽  
Author(s):  
Panagiotis Chatzipantelidis ◽  
Zoltan Horváth ◽  
Vidar Thomée

AbstractWe consider the initial boundary value problem for the homogeneous heat equation, with homogeneous Dirichlet boundary conditions. By the maximum principle the solution is nonnegative for positive time if the initial data are nonnegative. We complement in a number of ways earlier studies of the possible extension of this fact to spatially semidiscrete and fully discrete piecewise linear finite element discretizations, based on the standard Galerkin method, the lumped mass method, and the finite volume element method. We also provide numerical examples that illustrate our findings.

2019 ◽  
Vol 40 (4) ◽  
pp. 2415-2449
Author(s):  
D C Antonopoulos ◽  
V A Dougalis ◽  
G Kounadis

Abstract We consider a simple initial-boundary-value problem for the shallow water equations in one space dimension. We discretize the problem in space by the standard Galerkin finite element method on a quasiuniform mesh and in time by the classical four-stage, fourth order, explicit Runge–Kutta scheme. Assuming smoothness of solutions, a Courant number restriction and certain hypotheses on the finite element spaces, we prove $L^{2}$ error estimates that are of fourth-order accuracy in the temporal variable and of the usual, due to the nonuniform mesh, suboptimal order in space. We also make a computational study of the numerical spatial and temporal orders of convergence, and of the validity of a hypothesis made on the finite element spaces.


Author(s):  
C. M. Chen ◽  
V. Thomée

AbstractFor the heat equation in two space dimensions we consider semidiscrete and totally discrete variants of the lumped mass modification of the standard Galerkin method, using piecewise linear approximating functions, and demonstrate error estimates of optimal order in L2 and of almost optimal order in L∞.


2003 ◽  
Vol 3 (1) ◽  
pp. 45-58 ◽  
Author(s):  
Dejan Bojović

Abstract In this paper we consider the first initial boundary-value problem for the heat equation with variable coefficients in a domain (0; 1)x(0; 1)x(0; T]. We assume that the solution of the problem and the coefficients of the equation belong to the corresponding anisotropic Sobolev spaces. Convergence rate estimate which is consistent with the smoothness of the data is obtained.


1995 ◽  
Vol 05 (03) ◽  
pp. 351-365 ◽  
Author(s):  
V. SHUTYAEV ◽  
O. TRUFANOV

This paper is concerned with the numerical analysis of the mathematical model for a semiconductor device with the use of the Boltzmann equation. A mixed initial-boundary value problem for nonstationary Boltzmann-Poisson system in the case of one spatial variable is considered. A numerical algorithm for solving this problem is constructed and justified. The algorithm is based on an iterative process and the finite element method. A numerical example is presented.


2013 ◽  
Vol 5 (04) ◽  
pp. 548-568 ◽  
Author(s):  
Tao Lin ◽  
Yanping Lin ◽  
Xu Zhang

AbstractThis article extends the finite element method of lines to a parabolic initial boundary value problem whose diffusion coefficient is discontinuous across an interface that changes with respect to time. The method presented here uses immersed finite element (IFE) functions for the discretization in spatial variables that can be carried out over a fixed mesh (such as a Cartesian mesh if desired), and this feature makes it possible to reduce the parabolic equation to a system of ordinary differential equations (ODE) through the usual semi-discretization procedure. Therefore, with a suitable choice of the ODE solver, this method can reliably and efficiently solve a parabolic moving interface problem over a fixed structured (Cartesian) mesh. Numerical examples are presented to demonstrate features of this new method.


2018 ◽  
Vol 28 (06) ◽  
pp. 1199-1231
Author(s):  
Gerardo Huaroto ◽  
Wladimir Neves

In this paper, we study a fractional type degenerate heat equation posed in bounded domains. We show the existence of solutions for measurable and bounded non-negative initial data, and homogeneous Dirichlet boundary condition. The nonlocal diffusion effect relies on an inverse of the [Formula: see text]-fractional Laplacian operator, and the solvability is proved for any [Formula: see text].


Sign in / Sign up

Export Citation Format

Share Document