An Improved Linearity-Preserving Cell-Centered Scheme for Nonlinear Diffusion Problems on General Meshes

In this paper, we suggest a new vertex interpolation algorithm to improve an existing cell-centered finite volume scheme for nonlinear diffusion problems on general meshes. The new vertex interpolation algorithm is derived by applying a special limit procedure to the well-known MPFA-O method. Since the MPFA-O method for 3D cases has been addressed in some studies, the new vertex interpolation algorithm can be extended to 3D cases naturally. More interesting is that the solvability of the corresponding local system is proved under some assumptions. Additionally, we modify the edge flux approximation by an edge-based discretization of diffusion coefficient, and thus the improved scheme is free of the so-called numerical heat-barrier issue suffered by many existing cell-centered or hybrid schemes. The final scheme allows arbitrary continuous or discontinuous diffusion coefficients and can be applicable to arbitrary star-shaped polygonal meshes. A second-order convergence rate for the approximate solution and a first-order accuracy for the flux are observed in numerical experiments. In the comparative experiments with some existing vertex interpolation algorithms, the new algorithm shows obvious improvement on highly distorted meshes.     

A Finite Volume Scheme for Savage-Hutter Equations on Unstructured Grids

A Godunov-type finite volume scheme on unstructured grids is proposed to numerically solve the Savage-Hutter equations in curvilinear coordinate. We show the direct observation that the model isn't a Galilean invariant system. At the cell boundary, the modified Harten-Lax-van Leer (HLL) approximate Riemann solver is adopted to calculate the numerical flux. The modified HLL flux is not troubled by the lack of Galilean invariance of the model and it is helpful to handle discontinuities at free interface. Rigidly the system is not always a hyperbolic system due to the dependence of flux on the velocity gradient. Even so, our numerical results still show quite good agreements with reference solutions. The simulations for granular avalanche flows with shock waves indicate that the scheme is applicable.     

一维定常型对流占优扩散方程的一类迎风有限体积格式

本文针对一维定常型对流占优扩散方程提出了一类迎风有限体积格式 .该格式对对流项具有二阶精度 ,对扩散项保持一阶精度 ,符合对流占优扩散问题强对流、弱扩散的特点 .     

Meshfree Finite Volume Element Method for Constrained Optimal Control Problem Governed by Random Convection Diffusion Equations

In this paper, we investigate a stochastic meshfree finite volume element method for an optimal control problem governed by the convection diffusion equations with random coefficients. There are two contributions of this paper. Firstly, we establish a scheme to approximate the optimality system by using the finite volume element method in the physical space and the meshfree method in the probability space, which is competitive for high-dimensional random inputs. Secondly, the a priori error estimates are derived for the state,the co-state and the control variables. Some numerical tests are carried out to confirm the theoretical results and demonstrate the efficiency of the proposed method.     

A High-Order Finite Difference Scheme for 3D Unsteady Convection Diffusion Reaction Equations北大核心CSCD

针对三维非稳态对流扩散反应方程,构造了一种高精度紧致有限差分格式,对空间的离散采用四阶紧致差分方法,对时间的离散采用Taylor级数展开和余项修正技术,所提格式在时间上的精度为二阶、在空间上的精度为四阶。利用Fourier稳定性分析法证明了该格式是无条件稳定的。最后给出数值算例验证了理论结果。     

时间分数阶扩散方程的一个新的高阶有限差分/谱逼近

研究时间分数阶扩散方程,结合时间方向的有限差分格式和空间方向的Legendre Collocation谱方法,构造了一个高阶稳定数值格式.数值算例表明该格式是无条件稳定和长时间稳定的,其收敛阶为O(△t3-α+N-m),其中△t,N和m分别是时间步长,空间多项式阶数以及精确解的正则度.     

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

讨论了拟线性对流占优扩散问题的数值模拟.对对流部分采用特征线格式进行离散,以消除流动锋线前沿的数值弥散现象,保证格式的稳定性;而对扩散部分采用扩展混合有限元方法,同时逼近未知函数,未知函数的梯度及伴随向量函数.理论分析和数值算例表明,此方法是稳定的,具有最优L2逼近精度.     

半导体瞬态问题的修正迎风有限体积格式

半导体器件的瞬时状态由包含三个拟线性偏微分方程所组成的方程组的初边值问题来描述.其中电子位势方程足椭圆型的,电子和空穴浓度方程是对流扩散型的.对电子位势方程采用一次元有限体积法米逼近,对电子浓度和空穴浓度方程采用修正的迎风有限体积方法来逼近,并进行详细的理论分析,关于位势得到O(h Δt)阶的H1模误差估计结果,关于浓度得到O(h2 Δt)阶的L2模误差估计结果.最后,给出数值例子.     

Stability of Semi-implicit Finite Volume Scheme for Level Set Like Equation

We study numerical methods for level set like equations arising in image processing and curve evolution problems. Semi-implicit finite volume-element type schemes are constructed for the general level set like equation(image selective smoothing model) given by Alvarez et al.(Alvarez L, Lions P L, Morel J M. Image selective smoothing and edge detection by nonlinear diffusion II. SIAM J. Numer.Anal., 1992, 29: 845–866). Through the reasonable semi-implicit discretization in time and co-volume method for space approximation, we give finite volume schemes,unconditionally stable in L∞ and W1,2(W1,1) sense in isotropic(anisotropic) diffusion domain.     

同位非结构和结构网格摄动有限体积法（PFV）求解二维Navier-Stokes方程组

根据NS方程组的一阶迎风和二阶中心有限体积（UFV和CFV）格式,导出NS方程组迎风和中心摄动有限体积(UPFV和CPFV)格式．该格式通过把控制体界面质量通量摄动展开成网格间距的幂级数,并由守恒方程本身求得幂级数系数而获得．迎风和中心摄动有限体积格式使用了与一阶迎风和二阶中心格式相同的基点数和相同的表达形式,宜于计算机编程．顶盖驱动方腔流和驻点流标量输运的数值实验证明,迎风PFV格式比一阶UFV、二阶CFV格式有更高的精度,更高的分辨率．尤其是良好的鲁棒特性．对顶盖驱动方腔流,在Re数从102到104范围内,亚松弛系数可在0.3～0.8任取,收敛性能良好．     

A Second Order Nonconforming Triangular Mixed Finite Element Scheme for the Stationary Navier-Stokes Equations

In this paper,a nonconforming triangular mixed finite element scheme with second order convergence behavior is proposed for the stationary Navier-Stokes equations.The new nonconforming triangular element is taken as approximation space for the velocity and the linear element for the pressure.The convergence analysis is presented and optimal error estimates of both broken H1-norm and L2-norm for velocity as well as the L2-norm for the pressure are derived.     

A Finite Volume Implicit Euler Scheme for the Linearized Shallow Water Equations: Stability and Convergence

The linearized shallow water equations are discretized in space by a finite volume method and in time by an implicit Euler scheme. Stability and convergence of the scheme are proved.     

Sobolev方程的一类各向异性非协调有限元逼近

在各向异性网格下,分别讨论了Sobolev方程在半离散和全离散格式下的一类非协调有限元逼近,得到了与传统有限元方法相同的误差估计和一些超逼近性质.同时在半离散格式下,通过构造具有各向异性特征的插值后处理算子得到了整体超收敛结果.     

Convergence of Finite Volume Schemes for Hamilton-Jacobi Equations with Dirichlet Boundary Conditions

We study numerical methods for time-dependent Hamilton-Jacobi equations with weak Dirichlet boundary conditions. We first propose a new class of abstract monotone approximation schemes and get a convergence rate of 1/2 . Then, according to the abstract convergence results, by newly constructing monotone finite volume approximations on interior and boundary points, we obtain convergent finite volume schemes for time-dependent Hamilton-Jacobi equations with weak Dirichlet boundary conditions. Finally give some numerical results.     

A New Exponential Compact Scheme for the Two-Dimensional Unsteady Nonlinear Burgers' and Navier-Stokes Equations in Polar Cylindrical Coordinates

In this article, a new compact difference scheme is proposed in exponential form to solve two-dimensional unsteady nonlinear Burgers' and Navier-Stokes equations of motion in polar cylindrical coordinates by using half-step discretization. At each time level by using only nine grid points in space, the proposed scheme gives accuracy of order four in space and two in time. The method is directly applicable to the equations having singularities at boundary points. Stability analysis is explained in detail and many benchmark problems like Burgers', Navier-Stokes and Taylor-vortex problems in polar cylindrical coordinates are solved to verify the accuracy and efficiency of the scheme.     

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

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

A Continuous-Stage Modified Leap-Frog Scheme for High-Dimensional Semi-Linear Hamiltonian Wave Equations

Among the typical time integrations for PDEs, Leap-frog scheme is the well-known method which can easily be used. A most welcome feature of the Leap-frog scheme is that it has very simple scheme and is easy to be implemented. The main purpose of this paper is to propose and analyze an improved Leap-frog scheme, the so-called continuous-stage modified Leap-frog scheme for high-dimensional semi-linear Hamiltonian wave equations. To this end, under the assumption of periodic boundary conditions, we begin with the formulation of the nonlinear Hamiltonian equation as an abstract second-order ordinary differential equation (ODE) and its operator-variation-of-constants formula (the Duhamel Principle). Then the continuous-stage modified Leap-frog scheme is formulated. Accordingly, the convergence, energy preservation, symplecticity conservation and long-time behaviour of explicit schemes are rigorously analysed. Numerical results demonstrate the remarkable advantage and efficiency of the improved Leap-frog scheme compared with the existing mostly used numerical schemes in the literature.     

无结构网格有限体积方法及在热对流中的应用研究

首次将无结构三角网格的有限体积方法和压强连接半隐式算法相结合 ,用于求解非平行壁管道中的热对流问题 .并由此分析了化学汽相淀积薄膜生长的均匀性问题 .计算结果对于分析一类管道中热和动量输运现象均有普遍指导意义     

Superconvergence and Asymptotic Expansions for Bilinear Finite Volume Element Approximations

Aiming at the isoparametric bilinear finite volume element scheme, we initially derive an asymptotic expansion and a high accuracy combination formula of the derivatives in the sense of pointwise by employing the energy-embedded method on uniform grids. Furthermore, we prove that the approximate derivatives are convergent of order two. Finally, numerical examples verify the theoretical results.