三维半导体问题的迎风有限体积格式

半导体器件的瞬时状态由包含三个拟线性偏微分方程所组成的方程组的初边值问题来描述.其中电子位势方程是椭圆型的,电子和空穴浓度方程是对流扩散型的.作者对三维半导体模型问题采用四面体网格上的有限体积元方法进行逼近,具体地,对电子位势方程采用一次元有限体积法来逼近,对电子浓度和空穴浓度方程采用迎风有限体积方法来逼近,并进行了详细的理论分析,得到了O(h+\Delta t)阶的L^2模误差估计结果.     

半导体瞬态问题的混合迎风有限体积格式

半导体器件的瞬时状态由包含3个拟线性偏微分方程所组成的方程组的初边值问题来描述.在三角剖分的基础上,对椭圆型的电子位势方程采用混合有限体积元法来逼近,对对流扩散型的电子浓度和空穴浓度方程采用迎风有限体积元方法来逼近,并进行了详细的理论分析,得到了最优阶的误差估计结果.最后,针对混合有限体积元法和迎风有限体积元法分别单独使用以及两种方法结合使用的情形给出了不同的数值算例.     

热传导型半导体瞬态问题的迎风有限体积元方法

半导体瞬态问题的数学模型是由四个方程组成的非线性偏微分方程组的初边值问题所决定．其中电子浓度和空穴浓度方程往往是对流占优扩散问题,普通的方法已不适用,为此本文用迎风格式处理对流项部分,提出一种全离散迎风有限体积元方法,并进行收敛性分析,在最一般的情况下得到了一阶精度L2模误差估计结果．     

三维半导体器件问题在时空局部加密复合网格上的有限差分格式

半导体器件的瞬时状态由三个方程组成的非线性偏微分方程组的初边值问题决定,电子位势方程是椭圆型的,电子和空穴浓度方程是抛物型的．依据实际数值模拟的需要,提出了一类三维半导体问题在时间和空间上进行局部加密的复合网格上的有限差分形式,并给出了电子和空穴浓度的最大模误差估计,最后给出了数值算例．     

三维热传导型半导体问题的特征有限元方法和分析

本文研究三维热传导型半导体瞬态问题的特征有限元方法及其理论分析,其数学模型是一类非线性偏微分方程的初边值问题,对电子位势方程提出Ｇａｌｅｒｋｉｎ逼近;对电子,空穴浓度方程采用特征有限元逼近;对热传导方程采用对时间向后差分的Ｇａｌｅｒｋｉｎ逼近．应用微分方程先验估计理论和技巧得到了最优阶Ｌ＾２误差估计。     

热传导型半导体问题的配置方法及误差分析

热传导型半导体瞬态问题的数学模型是一类非线性偏微分方程的初边值问题．电子位势方程是椭圆型的，电子、空穴浓度方程及热传导方程是抛物型的．该文给出求解的配置方法，得到次优犔２模误差估计，并将配置法和Ｇａｌｅｒｋｉｎ有限元方法进行数值结果比较．     

二维热传导型半导体问题的一个二阶格式

热传导型半导体器件瞬态问题的数学模型由四个拟线性偏微分方程所组成的方程组的初边值问题来描述。其中电子位势方程是椭圆型的,电子和空穴浓度方程是对流扩散型的,温度方程为热传导型的。本文对二维热传导型半导体的一类混合初边值问题利用降阶法给出了一个二阶差分格式,并对其进行了详细的理论分析,得到了离散的犾２ 误差估计结果。     

半导体问题的交替方向有限元方法

该文用交替方向有限元方法求解半导体问题的Energy Trans port (ET)模型。对模型中椭圆型的电子位势方程采用交替方向迭代法，对流占优扩散的电子浓度和空穴浓度方程采用特征交替方向有限元方法，热传导方程利用Patch逼近采用交替方向有限元方法求解。利用微分方程的先验估计理论和技巧，分别得到了椭圆型方程和抛物型方程的最优H+1和L+2误差估计。     

半正定两相驱动问题的多步有限体积方法及其理论分析

考虑多维半正定两相驱动方程的初边值问题,在非结构网格上构造多步的迎风有限体积格式,利用微分方程先验估计理论证明了格式的离散模形式的误差估计为D(△t~2 h),其中△t和h分别表示时空步长．数值算例进一步验证了格式的有效性．     

三维热传导型半导体问题的交替方向特征有限元方法

提出交替方向特征有限元方法,对电场位势方程采用混合元格式,对电子,空穴浓度方程采用交替方向特征有限元格式,对温度方程提出交替方向格式.应用向量积计算及先验估计理论和技巧,得到最佳的L2误差估计.     

带非局部边界条件的非线性抛物方程的隐式差分解法

在准静态弹性力学中常遇到求解带有非局部边界条件的抛物方程初边值问题.本文构造了一个数值求解带有非局部边界条件的非线性抛物方程的隐式差分格式,利用离散泛函分析的知识和不动点定理证明了差分解是存在的,且在离散最大模意义下关于时间步长一阶收敛,关于空间步长二阶收敛,并给出了数值算例.     

A superconvergent $L^{\infty}$-error estimate of RT1 mixed methods for elliptic control problems with an integral constraint

In this paper, we investigate the superconvergence property of mixed finite element methods for a linear elliptic control problem with an integral constraint. The state and co-state are approximated by the order $k=1$ Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions. A superconvergent approximation of the control variable $u$ will be constructed by a projection of the discrete adjoint state. It is proved that this approximation have convergence order $h^{2}$ in $L^{\infty}$-norm. Finally, a numerical example is given to demonstrate the theoretical results.     

A global superconvergent $L^{\infty}$-error estimate of mixed finite element methods for semilinear elliptic optimal control problems

In this paper, we discuss the superconvergence of mixed finite element methods for a semilinear elliptic control problem with an integral constraint. The state and costate are approximated by the order $k=1$ Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions. Approximation of the optimal control of the continuous optimal control problem will be constructed by a projection of the discrete adjoint state. It is proved that this approximation has convergence order $h^{2}$ in $L^{\infty}$-norm. Finally, a numerical example is given to demonstrate the theoretical results.     

CONVERGENCE OF AN ALTERNATING A-φ SCHEME FOR QUASI-MAGNETOSTATIC EDDY CURRENT PROBLEM

We propose in this paper an alternating A-$\phi$ method for the quasi-magnetostatic eddy current problem by means of finite element approximations. Bounds for continuous and discrete error in finite time are given. And it is verified that provided the time step $\tau$ is sufficiently small, the proposed algorithm yields for finite time $T$ an error of $O(h+\tau^{1/2})$ in the $L^2$-norm for the magnetic field $H(= \mu^{-1} \nabla \times A)$, where $h$ is the mesh size, $\mu$ the magnetic permeability.     

An Approximation of Three-Dimensional Semiconductor Devices by Mixed Finite Element Method and Characteristics-Mixed Finite Element Method

The mathematical model for semiconductor devices in three space dimensions are numerically discretized. The system consists of three quasi-linear partial differential equations about three physical variables: the electrostatic potential, the electron concentration and the hole concentration. We use standard mixed finite element method to approximate the elliptic electrostatic potential equation. For the two convection-dominated concentration equations, a characteristics-mixed finite element method is presented. The scheme is locally conservative. The optimal $L^2$-norm error estimates are derived by the aid of a post-processing step. Finally, numerical experiments are presented to validate the theoretical analysis.     

A Two-Grid Finite Element Approximation for Nonlinear Time Fractional Two-Term Mixed Sub-Diffusion and Diffusion Wave Equations

In this paper, we develop a two-grid method (TGM) based on the FEM for 2D nonlinear time fractional two-term mixed sub-diffusion and diffusion wave equations. A two-grid algorithm is proposed for solving the nonlinear system, which consists of two steps: a nonlinear FE system is solved on a coarse grid, then the linearized FE system is solved on the fine grid by Newton iteration based on the coarse solution. The fully discrete numerical approximation is analyzed, where the Galerkin finite element method for the space derivatives and the finite difference scheme for the time Caputo derivative with order $\alpha\in(1,2)$ and $\alpha_{1}\in(0,1)$. Numerical stability and optimal error estimate $O(h^{r+1}+H^{2r+2}+\tau^{\min\{3-\alpha,2-\alpha_{1}\}})$ in $L^{2}$-norm are presented for two-grid scheme, where $t,$ $H$ and $h$ are the time step size, coarse grid mesh size and fine grid mesh size, respectively. Finally, numerical experiments are provided to confirm our theoretical results and effectiveness of the proposed algorithm.     

对流扩散方程的解

关注如下的对流扩散方程 $$ u_{t}=\text{div}(|\nabla u^{m}|^{p-2}\nabla u^{m})+\sum_{i=1}^{N}\frac{\partial b_{i}(u^{m})}{\partial x_{i}} $$ 的初边值问题. 若 $p>1+\frac{1}{m}$, 通过考虑正则化问题的解 $u_{k}$, 利用 Moser 迭代技巧, 得到了$u_{k}$ 的 $L^{\infty}$ 模与 梯度 $\nabla u_{k}$ 的 $L^{p}$ 模的局部有界性. 利用紧致性定理, 得到了对流扩散方程本身解的存在性. 若 $p<1+\frac{1}{m},\ p>2$ 或者 $p=1+\frac{1}{m}$, 利用类似的方法可以得到解的存在性. 证明了解的唯一性, 同时讨论了正性和熄灭性等解的性质.     

具有调和曲率与有限$L^p$曲率的完备流形

Let(M~n, g)(n ≥ 3) be an n-dimensional complete Riemannian manifold with harmonic curvature and positive Yamabe constant. Denote by R and R?m the scalar curvature and the trace-free Riemannian curvature tensor of M, respectively. The main result of this paper states that R?m goes to zero uniformly at infinity if for p ≥ n, the L~p-norm of R?m is finite.As applications, we prove that(M~n, g) is compact if the L~p-norm of R?m is finite and R is positive, and(M~n, g) is scalar flat if(M~n, g) is a complete noncompact manifold with nonnegative scalar curvature and finite L~p-norm of R?m. We prove that(M~n, g) is isometric to a spherical space form if for p ≥n/2, the L~p-norm of R?m is sufficiently small and R is positive.In particular, we prove that(M~n, g) is isometric to a spherical space form if for p ≥ n, R is positive and the L~p-norm of R?m is pinched in [0, C), where C is an explicit positive constant depending only on n, p, R and the Yamabe constant.     

Compact finite difference schemes of the time fractional Black-Scholes model

In this paper, three compact difference schemes for the time-fractional Black-Scholes model governing European option pricing are presented. Firstly, in order to obtain the fourth-order accuracy in space by applying the Pad\''{e} approximation, we eliminate the convection term of the B-S equation by an exponential transformation. Then the time fractional derivative is approximated by $L1$ formula, $L2 - 1_\sigma$ formula and $L1 - 2$ formula respectively, and three compact difference schemes with oders $O(\Delta t^{2-\alpha}+h ^4)$, $O(\Delta t^{2}+h ^4)$ and $O(\Delta t^{3-\alpha}+h ^4)$ are constructed. Finally, numerical example is carried out to verify the accuracy and effectiveness of proposed methods, and the comparisons of various schemes are given. The paper also provides numerical studies including the effect of fractional orders and the effect of different parameters on option price in time-fractional B-S model.