Search

In this paper, we propose a SAV method for the Cahn–Hilliard type phase field model with new dynamic boundary conditions. By using the continuous finite element method in space and the backward difference method in time, the fully discrete numerical schemes preserving the energy law are constructed. Numerical examples show that the proposed scheme can simulate the phase field model well even in a relatively rough grid. [ABSTRACT FROM AUTHOR]

In this paper, we study a weak Galerkin (WG) finite element method for semiconductor device simulations. We consider the one-dimensional drift–diffusion (DD) and high-field (HF) models, which involves not only first derivative convection terms but also second derivative diffusion terms, as well as a coupled Poisson potential equation. The main difficulties in the analysis include the treatment of the nonlinearity and coupling of the models. The weak Galerkin finite element method adopts piecewise polynomials of degree k for the approximations of electron concentration and electric potential in the interior of elements, and piecewise polynomials of degree k + 1 for the discrete weak derivative space. The optimal order error estimates in a discrete H 1 norm and the standard L 2 norm are derived. Numerical experiments are presented to illustrate our theoretical analysis. Moreover, numerical schemes also work out for the discontinuous diffusion coefficient problems. [ABSTRACT FROM AUTHOR]