线性对流占优扩散方程的后验误差估计

对于线性对流占优扩散方程,采用特征线有限元方法离散时间导数项和对流项,用分片线性有限元离散空间扩散项,并给出了一致的后验误差估计,其中估计常数不依赖与扩散项系数。     

对流占优扩散方程的分裂特征混合有限元方法

利用修正的特征线方法,构建一类求解对流占优扩散方程的分裂特征混合有限元算法.在新的算法中,混合系统的系数矩阵对称正定,且原未知函数u与流函数σ=-ε▽u可分离求解.推导了加权能量模意义下的最优阶误差估计,并给出数值算例验证理论上的分析结果.     

对流扩散方程特征线三角元法的一致估计

利用三角形线性元的积分恒等式,给出了二维非定常对流占优扩散方程的特征线有限元解和真解的一致最优估计,并利用插值后处理算子,得到了有限元解梯度的一致超收敛估计,即只与初值和右端项有关,而与ε无关.     

对流扩散方程的三层ENO-MMOCAA差分方法

本文把多步修正特征线法[1],MMOCAA差分方法[2]及ENO插值[3]相结合,提出了求解对流扩散方程的多步ENO-MMOCAA差分方法.该方法关于时间及空间都具有二阶以上的精度且可避免在解的大梯度附近产生振荡.本文给出了格式的误差估计及数值算例.     

线性对流占优扩散方程的扩展特征混合有限元法

本文针对线性对流占优扩散方程提出了一种新型数值模拟方法一扩展特征混合有限元法,即对对流部分沿特征线方向离散,而对扩散部分采用扩展混合有限元方法,同时高精度逼近未知函数,未知函数的梯度及伴随向量函数,通过严格的数值分析,得到其最优L^2模误差估计。     

对流扩散方程的三层WENO-MMOCAA差分方法

本文把三层修正特征线法,MMOCAA 差分方法及WENO 插值相结合,提出了求解对流扩散方程的三层WENO-MMOCAA 差分格式.此格式关于时间具有二阶精度,关于空间具有二阶以上精度且可避免基于二次以上Lagrange 插值的三层MMOCAA 差分方法在解的大梯度附近所产生的振荡.本文使用新的分析方法,给出了格式的误差估计.本文的数值算例表明新格式可消除振荡.     

对流扩散方程一类改进的特征线修正有限元方法

１引言在地下水污染,地下渗流驱动,核污染,半导体等问题的数值模拟中,均涉及抛物型对流扩散方程（或方程组）的数值求解问题．这些对流扩散型偏微分方程（或方程组）具有共同的特点：对流的影响远大于扩散的影响,即对流占优性,对流占优性给问题的数值求解带来许多困难,因此对流占优问题的有效数值解法一直是计算数学中重要的研究内容．用通常的差分法或有限元法进行数值求解将出现数值振荡．为了克服数值振荡,提出各种迎风方法和修正的特征方法并在这些问题上得到成功的实际应用、８０年代,Ｄｏｕｇｌａｓ和Ｒｕｓｓｅｌｌ［２］等…     

二维非均匀水沙模型的多步特征-混合元数值模拟

通过对方程的对流部分采用沿着特征线方向向后两步差分格式进行离散,而对扩散部分采用混合有限元格式进行离散,从而利用多步特征-混合有限元方法对平面非均匀水沙模型进行了数值模拟,给出了相应的误差分析及数值算例.     

一类对流扩散方程的特征-修正差分格式及最大模估计

1 引 言许多物理问题的数学模型都可归结为对流扩散型方程（组）,这类模型具有较强的双曲性质.为了在数值计算中更好地反映这种性质,Douglas,Russell[1]提出了特征法的思想.特征法的主要特点是将模型方程的对流项隐含在解沿特征方向的导数之中,以沿特征方向的向后差商逼近取代了通常方法中沿时间方向的向后差商逼近.由于解在特征方向上的变     

对流扩散方程特征线双线性元的一致估计

利用双线性元的积分恒等式,给出了二维非定常对流占优扩散方程的特征线有限元解和真解的一致误差估计,并利用插值后处理算子给出了有限元解梯度的一致超收敛估计,即上述误差与ε无关,而仅与右端f和初值u_0有关.     

带广义边界条件的四阶抛物型方程的混合间断时空有限元法

构造具有广义边界条件的四阶线性抛物型方程的混合间断时空有限元格式,利用混合有限元方法将高阶方程降阶,利用空间连续而时间允许间断的时空有限元方法离散方程,证明了离散解的存在唯一性,稳定性和收敛性,并给出数值算例验证了方法的有效性.     

一类四阶抛物型积分-微分方程的混合间断时空有限元法

构造四阶抛物型积分-微分方程的混合间断时空有限元格式,利用混合有限元方法将高阶方程降阶,利用空间连续而时间允许间断的时空有限元方法离散方程,证明离散解的稳定性,存在唯一性和收敛性.     

Discontinuous Galerkin finite element heterogeneous multiscale method for advection–diffusion problems with multiple scales

A discontinuous Galerkin finite element heterogeneous multiscale method is proposed for advection–diffusion problems with highly oscillatory coefficients. The method is based on a coupling of a discontinuous Galerkin discretization for an effective advection–diffusion problem on a macroscopic mesh, whose a priori unknown data are recovered from micro finite element calculations on sampling domains within each macro element. The computational work involved is independent of the high oscillations in the problem at the smallest scale. The stability of our method (depending on both macro and micro mesh sizes) is established for both diffusion dominated and advection dominated regimes without any assumptions about the type of heterogeneities in the data. Fully discrete a priori error bounds are derived for locally periodic data. Numerical experiments confirm the theoretical error estimates.     

A Modified Weak Galerkin Finite Element Method for Sobolev Equation

For Sobolev equation, we present a new numerical scheme based on a modified weak Galerkin finite element method, in which differential operators are approximated by weak forms through the usual integration by parts. In particular, the numerical method allows the use of discontinuous finite element functions and arbitrary shape of element. Optimal order error estimates in discrete $H^1$ and $L^2$ norms are established for the corresponding modified weak Galerkin finite element solutions. Finally, some numerical results are given to verify theoretical results.     

An approximation of incompressible miscible displacement in porous media by mixed finite elements and symmetric finite volume element method of characteristics

A nonlinear system of two coupled partial differential equations models miscible displacement of one incompressible fluid by another in a porous medium. A sequential implicit time‐stepping procedure is defined, in which the pressure and Darcy velocity of the mixture are approximated by a mixed finite element method and the concentration is approximated by a combination of a modified symmetric finite volume element method and the method of characteristics. Optimal order convergence in H¹ and in L² are proved for full discrete schemes. Finally, some numerical experiments are presented. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A characteristic mixed element method for displacement problems of compressible flow in porous media

A new characteristic mixed element scheme is formulated to solve numerically displacement problems of compressible fluids in porous media. A new mixed finite element method is introduced to solve the pressure equation of parabolic type, in which the mixed element system is symmetric positive definite and the pressure equation is separated from the flux equation. The modified method of characteristics is used to treat convection-dominated diffusion equations of the concentrations. The convergence with optimal accuracy is proved under the general condition. Project supported in part by China State Major Key Project for Basic Researches, Doctoral Station Foundation and TCTPF of China State Education Commission.     

多孔介质二相驱动问题一种修正特征有限元方法的迭代解法及其收敛性分析

0 引言 多孔介质二相驱动问题的数学模型是由压力方程与浓度方程组成的偏微分方程组的初边值问题.关于该问题的数值解问题,已有大量的文献.为了得到最优的L~2-模误差估计,好多方法用混合元方法解压力方程.我们知道,混合元法得到的方程组系数矩阵是非正定的,从而解混合元比解标准元要困难得多,虽然许多人研究了混合元方法的求解问题,但到目前为止,还没有看到令人满意的好的算法.为了避开对混合元的求解,著名学者T.F.Russell考虑了用标准有限元方法解压力方程,用特征有限元方法解浓度方程的求解方法及其迭代解法,对只有分子扩散的二相驱动问题得到了最优的L~2模误差估计,对有机械弥散的一般二相驱动问题得不到最优的L~2模误差估计,同时在收敛性证明中要求压力有限元空间的指数至少是二.     

The modified method of characteristics with mixed finite element domain decomposition procedures for the transient behavior of a semiconductor device

For the transient behavior of a semiconductor device, the modified method of characteristics with mixed finite element domain decomposition procedures applicable to parallel arithmetic is put forward. The electric potential equation is described by the mixed finite element method, and the electric, hole concentration and heat conduction equations are treated by the modified method of characteristics finite element domain decomposition methods. Some techniques, such as calculus of variations, domain decomposition, characteristic method, energy method, negative norm estimate and prior estimates and techniques are employed. Optimal order estimates in L² norm are derived for the error in the approximation solution. Thus the well‐known theoretical problem has been thoroughly and completely solved.© 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 353–368 2012     

A uniform estimate for the MMOC for two‐dimensional advection‐diffusion equations

We prove an optimal‐order error estimate in a weighted energy norm for the modified method of characteristics (MMOC) and the modified method of characteristics with adjusted advection (MMOCAA) for two‐dimensional time‐dependent advection‐diffusion equations, in the sense that the generic constants in the estimates depend on certain Sobolev norms of the true solution but not on the scaling diffusion parameter ε. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

求解Brinkman方程的修正弱有限元方法

本文针对Brinkman方程引入了一种修正弱Galerkin(MWG)有限元方法.我们通过具有两个离散弱梯度算子的变分形式来逼近模型. 在MWG方法中, 分别用次数为$k$和$k-1$的不连续分段多项式来近似速度函数$u$和压力函数$p$. MWG方法的主要思想是用内部函数的平均值代替边界函数. 因此, 与WG方法相比, MWG方法在不降低准确性的同时, 具有更少的自由度, 对于任意次数不超过$k-1$ 的多项式,MWG方法均可以满足稳定性条件. MWG 方法具有高度的灵活性, 它允许在具有一定形状正则性的任意多边形或多面体上使用不连续函数. 针对$H^1$和$L^22$范数下的速度和压力近似解, 建立了最优阶误差估计. 数值算例表明了该方法的准确性, 收敛性和稳定性.