Sine-Gordon方程的混合有限体积元方法及数值模拟

研究了在Dirichlet边界条件和Neumann边界条件下一维sine-Gordon方程的混合有限体积元方法.通过引入将试探函数空间映射到检验函数空间的迁移算子γh,结合混合有限元方法和有限体积元方法,构造了半离散格式,时间显式和隐式全离散混合有限体积元格式.给出了显格式离散解的稳定性分析,并得到了三种格式的最优阶误差估计.最后,给出数值算例来验证理论分析结果和数值格式的有效性.     

四阶抛物偏微分方程的H1-Galerkin混合元方法及数值模拟

到目前为止, H¹-Galerkin 混合有限元方法研究的问题仅局限于二阶发展方程. 然而对于高阶发展方程, 特别是重要的四阶发展方程问题的研究却没有出现. 本文首次提出四阶发展方程的H¹-Galerkin 混合有限元方法, 为了给出理论分析的需要, 我们考虑四阶抛物型发展方程. 通过引进三个适当的中间辅助变量, 形成四个一阶方程组成的方程组系统, 提出四阶抛物型方程的H¹-Galerkin 混合有限元方法. 得到了一维情形下的半离散和全离散格式的最优收敛阶误差估计和多维情形的半离散格式误差估计, 并采用迭代方法证明了全离散格式的稳定性. 最后, 通过数值例子验证了提出算法的可行性. 在一维情况下我们能够同时得到未知纯量函数、一阶导数、负二阶导数和负三阶导数的最优逼近解, 这一点是以往混合元方法所不能得到的.     

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

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

分数阶Burgers方程的有限元计算(英文)

建立了一维和二维分数阶Burgers方程的有限元格式.时间分数阶导数使用L1方法离散,空间方向使用有限元方法离散.通过选择合适的基函数,将离散后的方程转化成一个非线性代数方程组,并应用牛顿迭代方法求解.数值实验显示出了方法的有效性.     

热传导型半导体器件的二次元特征有限体积元方法及分析：H^1模估计

该文构造了热传导型半导体器件的全离散特征有限体积元格式,将特征线方法与有限体积元方法相结合,采用Lagrange型分片二次多项式空间和分片常数函数空间分别作为试探函数和检验函数空间,并进行误差分析,得到了最优阶 H1模误差估计结果.     

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

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

应用广义Laguerre函数的Navier-Stokes方程外部问题混合谱方法

研究应用广义Laguerre函数的四阶非线性偏微分方程外部问题混合谱方法.构造了圆外Navier-Stokes方程流函数形式的混合谱方法,数值结果显示了该方法在空间方向的谱精度.     

二维热传导型半导体器件的特征有限体积元方法和H1模分析

将特征线方法与有限体积元方法相结合,采用分片线性函数和分片常数函数分别作为有限体积元方法的试探函数和检验函数空间,构造了热传导型半导体器件的全离散特征有限体积元格式.并进行收敛性分析,在一般的条件下得到了最优阶H1模误差估计结果.     

空间分数阶电报方程的格子Boltzmann方法

应用格子Boltzmann方法（LBM）对Riemann Liouville空间分数阶电报方程进行了数值模拟研究.首先,将分数阶算子中的积分项进行离散化处理,并进行了收敛阶分析.然后,构建了带修正函数项的一维三速度(D1Q3)的LBM演化模型.利用Chapman Enskog多尺度技术和Taylor展开技术,推导出各平衡态分布函数和修正函数的具体表达式,准确地从所建的演化模型恢复出宏观方程.最后,数值计算结果表明该模型是稳定、有效的.     

空间分数阶偏微分方程非标准有限差分方法的稳定性和收敛性（英文）

本文研究空间分数阶偏微分方程非标准有限差分方法数值解的相关问题.采用Grünwald-Letnikov公式和平移Grünwald-Letnikov公式分别对两个空间分数阶导数进行离散.再运用带有时间和空间步长的分母函数构造非标准有限差分方法.进而利用von Neumann分析方法对差分格式的稳定性和收敛性进行研究,获得了一些新的结果.数值例子验证了非标准有限差分方法用于求解空间分数阶偏微分方程的有效性.     

An expanded mixed covolume method for sobolev equation with convection term on triangular grids

The expanded mixed covolume method for the two‐dimensional Sobolev equation with convection term is developed and studied. This method uses the lowest‐order Raviart‐Thomas mixed finite element space as the trial function space. By introducing a transfer operator γ_h which maps the trial function space into the test function space and combining expanded mixed finite element with mixed covolume method, the continuous‐in‐time, discrete‐in‐time expanded mixed covolume schemes are constructed, and optimal error estimates for these schemes are obtained. Numerical results are given to examine the validity and effectiveness of the proposed schemes.© 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A characteristics-mixed covolume method for a convection-dominated transport problem

In this paper, we propose a characteristics-mixed covolume method for approximating the solution to a convection dominated transport problem. The method is a combination of characteristic approximation to handle the convection term in time and mixed covolume method spatial approximation to deal with the diffusion term. The velocity and press are approximated by the lowest order Raviart-Thomas mixed finite element space on rectangles. The projection of a mixed covolume element is introduced. We prove its first order optimal rate of convergence for the approximate velocities in the L² norm as well as for the approximate pressures in the L² norm.     

线性抛物型积分微分方程的扩展混合体积元方法

1 引言 考虑线性抛物型积分微分方程初边值问题: {pt(x,t)-▽.{A(x,t)▽p(x,t) +∫t0 B(x,t,τ)▽p(x,τ)dτ}=f(x,t),(x,t)∈Ω×(0,T],(1.1) p(x,0):p0(x), x∈Ω, p(x,t)=0, (x,t)∈(a)Ω×(0,T]. 这里x=(x,y),Ω=(a,b)×(c,d),(e)Ω是区域Ω的边界,p为未知函数,A=(aij)2×2为已知的对称正定矩阵,B=(bij)2×2为已知矩阵,而且aij,bij,(aij)t(i,j=1,2)光滑有界,f∈L2(Ω).     

Superconvergence of mixed covolume method for elliptic problems on triangular grids

In this paper, we consider the superconvergence of a mixed covolume method on the quasi-uniform triangular grids for the variable coefficient-matrix Poisson equations. The superconvergence estimates between the solution of the mixed covolume method and that of the mixed finite element method have been obtained. With these superconvergence estimates, we establish the superconvergence estimates and the L^∞$L^{\infty }$-error estimates for the mixed covolume method for the elliptic problems. Based on the superconvergence of the mixed covolume method, under the condition that the triangulation is uniform, we construct a post-processing method for the approximate velocity which improves the order of approximation of the approximate velocity.     

A general mixed covolume framework for constructing conservative schemes for elliptic problems

We present a general framework for the finite volume or covolume schemes developed for second order elliptic problems in mixed form, i.e., written as first order systems. We connect these schemes to standard mixed finite element methods via a one-to-one transfer operator between trial and test spaces. In the nonsymmetric case (convection-diffusion equation) we show one-half order convergence rate for the flux variable which is approximated either by the lowest order Raviart-Thomas space or by its image in the space of discontinuous piecewise constants. In the symmetric case (diffusion equation) a first order convergence rate is obtained for both the state variable (e.g., concentration) and its flux. Numerical experiments are included.

     

Mixed covolume method for parabolic problems on triangular grids

In this paper, we present a mixed covolume method for parabolic equations on triangular grids. This method use the lowest order Raviart–Thomas (R–T) mixed finite element space as the trial space. We prove the optimal order of convergence for the approximate pressure and velocity in L²-norm. Furthermore, we obtain the quasi-optimal error estimates for the approximate pressure in L^∞-norm.     

Symmetric mixed covolume methods for parabolic problems

We present a mixed covolume method for a system of first order partial differential equations resulting from the mixed formulation of the general self‐adjoint parabolic problem with a variable nondiagonal diffusion tensor. The lowest order Raviart‐Thomas mixed element space on rectangles is used. We prove the first order optimal rate of convergence for approximate pressure as well as for approximate velocity. We also prove the second order superconvergence both for approximate velocity and pressure in certain discrete norms. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 561–583, 2002     

An expanded mixed covolume method for elliptic problems

We consider the mixed covolume method combining with the expanded mixed element for a system of first‐order partial differential equations resulting from the mixed formulation of a general self‐adjoint elliptic problem with a full diffusion tensor. The system can be used to model the transport of a contaminant carried by a flow in porous media. We use the lowest order Raviart‐Thomas mixed element space. We show the first‐order error estimate for the approximate solution in L² norm. We show the superconvergence both for pressure and velocity in certain discrete norms. We also get a finite difference scheme by using proper approximate integration formulas. Finally we give some numerical examples. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A discontinuous mixed covolume method for elliptic problems

We develop a discontinuous mixed covolume method for elliptic problems on triangular meshes. An optimal error estimate for the approximation of velocity is obtained in a mesh-dependent norm. First-order L²-error estimates are derived for the approximations of both velocity and pressure.     

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

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