Mixed Finite Element-Second Order Upwind Fractional Step Difference Scheme of Darcy–Forchheimer Miscible Displacement and Its Numerical Analysis

2021 ◽  
Vol 86 (2) ◽  
Author(s):  
Yirang Yuan ◽  
Changfeng Li ◽  
Qing Yang
2017 ◽  
Vol 10 (3) ◽  
pp. 541-561
Author(s):  
Yirang Yuan ◽  
Qing Yang ◽  
Changfeng Li ◽  
Tongjun Sun

AbstractA coupled mathematical system of four quasi-linear partial differential equations and the initial-boundary value conditions is presented to interpret transient behavior of three dimensional semiconductor device with heat conduction. The electric potential is defined by an elliptic equation, the electron and hole concentrations are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite element approximation is used to get the electric field potential and one order of computational accuracy is improved. Two concentration equations and the heat conduction equation are solved by a fractional step scheme modified by a second-order upwind difference method, which can overcome numerical oscillation, dispersion and computational complexity. This changes the computation of a three dimensional problem into three successive computations of one-dimensional problem where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by prior estimate theory and other special techniques of partial differential equations. This type of parallel method is important in numerical analysis and is most valuable in numerical application of semiconductor device and it can successfully solve this international famous problem.


Mathematics ◽  
2020 ◽  
Vol 8 (2) ◽  
pp. 187
Author(s):  
Yaxin Hou ◽  
Cao Wen ◽  
Hong Li ◽  
Yang Liu ◽  
Zhichao Fang ◽  
...  

In this article, some high-order time discrete schemes with an H 1 -Galerkin mixed finite element (MFE) method are studied to numerically solve a nonlinear distributed-order sub-diffusion model. Among the considered techniques, the interpolation approximation combined with second-order σ schemes in time is used to approximate the distributed order derivative. The stability and convergence of the scheme are discussed. Some numerical examples are provided to indicate the feasibility and efficiency of our schemes.


2006 ◽  
Vol 6 (2) ◽  
pp. 154-177 ◽  
Author(s):  
E. Emmrich ◽  
R.D. Grigorieff

AbstractIn this paper, we study the convergence of the finite difference discretization of a second order elliptic equation with variable coefficients subject to general boundary conditions. We prove that the scheme exhibits the phenomenon of supraconvergence on nonuniform grids, i.e., although the truncation error is in general of the first order alone, one has second order convergence. All error estimates are strictly local. Another result of the paper is a close relationship between finite difference scheme and linear finite element methods combined with a special kind of quadrature. As a consequence, the results of the paper can be viewed as the introduction of a fully discrete finite element method for which the gradient is superclose. A numerical example is given.


Sign in / Sign up

Export Citation Format

Share Document