奇异摄动对流扩散方程组的参数一致的二阶有限差分格式

针对一类奇异摄动对流扩散方程组,在■网格上构造了经典的迎风有限差分格式,并利用闸函数方法证明了数值方法为一阶收敛.在此基础上,设计了一个Richardson外推格式,并严格证明了外推方法的精度为二阶一致收敛.数值实验验证了本文的理论结果.     

一类奇异摄动对流扩散方程组的自适应网格方法的收敛性分析

基于经典的迎风有限差分方法,本文讨论一类奇异摄动对流扩散方程组的自适应网格算法.首先,利用Taylor级数展开,给出离散格式的局部截断误差估计.然后,利用等分布原理和极大值原理,证明基于弧长控制函数的自适应网格算法是一阶收敛的.最后的数值实验验证了本文的理论结果.     

基于有限体积法的非结构网格大涡模拟离散方法研究

非结构网格下的大涡模拟是解决复杂几何体高Reynolds(雷诺)数流动的有效途径.首先,基于有限体积法,研究了对流项和扩散项非结构网格下的离散方法.研究结果表明:基于TVD(total variation diminishing)限制器的限制中心差分格式保证了对流项的二阶精度并抑制了非物理振荡,同时,线性迎风格式虽然稳定,但数值耗散过大,且不能保证有界,中心差分格式引起了周期性非物理振荡;扩散项的超松弛非正交修正减小了网格非正交带来的离散误差,但修正系数须根据网格非正交的程度进行合理选取.为验证所述离散方法对大涡模拟的适用性,数值计算了Re=1.14×10~6下的非定常三维小球绕流,计算方法包括:计算网格用基于Delaunay三角剖分和Netgen前沿推进算法的四面体非结构网格;湍流模型用改进的延迟分离涡大涡模型;在离散格式的选取上,对流项用限制中心差分,扩散项加入非正交修正,插值格式用最小二乘法,时间项用二阶后向差分.计算结果表明,所用离散方法稳定收敛并且与实验数据基本吻合.     

非饱和水流问题的迎风差分法及其数值模拟

本文考察了非饱和水流问题模型方程的守恒型迎风差分法.我们基于有限体积方法建立的非饱和流动的守恒形式,分别提出了一阶和二阶迎风差分格式,并对差分格式进行了误差估计,给出了收敛性定理.最后,数值模拟验证了计算格式的有效性.     

求解奇异摄动问题的一个耦合差分格式

本文考察奇异摄动问题(1.1).在一特殊的非均匀网格上,将不稳定、二阶精度的中心差格式和稳定、一阶精度的Abrahamsson-Keller-Kreiss箱子格式相耦合,得到了一个二阶一致收敛的差分格式.最后给出了数值结果.     

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

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

二维奇异摄动对流扩散方程的自适应移动网格算法

针对二维奇异摄动对流扩散方程,在任意网格下给出了经典的迎风有限差分格式.利用二元多项式插值技术,推导出一阶最大范数的后验误差估计,并以此设计了一个自适应网格生成算法.数值实验表明本文构造的自适应移动网格算法是有效的.     

平面驱动方腔流的混沌研究

应用二维涡量-流函数形式的不可压N-S方程组的一致四阶精度的紧致格式,对高Re下平面驱动方腔问题数值模拟.利用混沌时间序列分析的手段,定性、定量的研究高Re下平面驱动方腔内流动系统,从规则状态到混沌状态的转变,并详细地给出了其混沌特征.     

迎风紧致格式与驱动方腔流动问题的直接数值模拟

本文给出了一种求解不可压缩流动问题的高精度差分格式,即迎风紧致格式.出发方程采用二维非定常原始变量Naiver-Stokes方程组.在差分方程中,对流项采用三阶精度的迎风紧致差分,其余空间导数项采用四阶紧致差分.本文利用该差分格式在等距网格上数值模拟了驱动方腔流动中的分离涡运动.在257×257的细网格上,Re数最高计算到10000.Re≤5000时的计算结果与前人结果符合得很好.当Re≥7500时发现流动不存在定常层流解而为非定常周期性解,并首次给出了非定常解的结果。     

两点混合边值问题的紧有限体积格式

本文针对常系数和变系数两点混合边值问题提出一种紧有限体积格式,该格式形成的线性代数方程组具有三对角性质,可以使用追赶法求解.证明格式按照H1半范数具有四阶收敛精度.利用节点计算值,给出单元中点值和一阶导数值的高精度后处理计算公式,这两个公式同样具有四阶精度.数值算例验证了理论分析的正确性,并说明了格式的有效性.     

An upwind finite volume scheme and its maximum‐principle‐preserving ADI splitting for unsteady‐state advection‐diffusion equations

We develop an upwind finite volume (UFV) scheme for unsteady‐state advection‐diffusion partial differential equations (PDEs) in multiple space dimensions. We apply an alternating direction implicit (ADI) splitting technique to accelerate the solution process of the numerical scheme. We investigate and analyze the reason why the conventional ADI splitting does not satisfy maximum principle in the context of advection‐diffusion PDEs. Based on the analysis, we propose a new ADI splitting of the upwind finite volume scheme, the alternating‐direction implicit, upwind finite volume (ADFV) scheme. We prove that both UFV and ADFV schemes satisfy maximum principle and are unconditionally stable. We also derive their error estimates. Numerical results are presented to observe the performance of these schemes. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 211–226, 2003     

High order upwind scheme for modelling turbulent shallow water flow in hydraulic structures

An unstructured finite volume model for quasi-2D free surface flow with wet-dry fronts and turbulence modelling is presented. The convective flux is discretised with either a an hybrid second-order/first-order scheme, or a fully second order scheme, both of them upwind Godunov's schemes based on Roe's average. The hybrid scheme uses a second order discretisation for the two unit discharge components, whilst keeping a first order discretisation for the water depth [2]. In such a way the numerical diffusion is much reduced, without a significant reduction on the numerical stability of the scheme, obtaining in such a way accurate and stable results. It is important to keep the numerical diffusion to a minimum level without loss of numerical stability, specially when modelling turbulent flows, because the numerical diffusion may interfere with the real turbulent diffusion. In order to avoid spurious oscillations of the free surface when the bathymetry is irregular, an upwind discretisation of the bed slope source term [4] with second order corrections is used [2]. In this way a fully second order scheme which gives an exact balance between convective flux and bed slope in the hydrostatic case is obtained. The k – ε equations are solved with either an hybrid or a second order scheme. In all the numerical simulations the importance of using a second order upwind spatial discretisation has been checked [1]. A first order scheme may give rather good predictions for the water depth, but it introduces too much numerical diffusion and therefore, it excessively smooths the velocity profiles. This is specially important when comparing different turbulence models, since the numerical diffusion introduced by a first order upwind scheme may be of the same order of magnitude as the turbulent diffusion. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Second order modified method of characteristics mixed defect-correction finite element method for time dependent Navier–Stokes problems

In this paper, a second order modified method of characteristics defect-correction (SOMMOCDC) mixed finite element method for the time dependent Navier–Stokes problems is presented. In this method, the hyperbolic part (the temporal and advection term) are treated by a second order characteristics tracking scheme, and the non-linear term is linearized at the same time. Then, we solve the equations with an added artificial viscosity term and correct this solution by using the defect-correction technique. The error analysis shows that this method has a good convergence property. In order to show the efficiency of the SOMMOCDC mixed finite element method, we first present some numerical results of an analytical solution problem, which agrees very well with our theoretical results. Then, we give some numerical results of lid-driven cavity flow with the Reynolds number Re = 5,000, 7,500 and 10,000. From these numerical results, we can see that the schemes can result in good accuracy, which shows that this method is highly efficient.     

任意精度的三点紧致显格式及其在CFD中的应用

通过在泰勒级数展开中运用逐阶迭代的方法,推导出了空间任意精度的三点紧致显格式的表达式,又由Fourier分析法得到了格式的数值弥散和耗散特性．与以往的高精度紧致差分格式不同,提出的格式不用隐式求解代数方程组并且可以达到任意精度．通过方波问题和顶盖方腔流的算例表明,格式在稀疏网格下可以得到很高的精度,不仅能节省计算量,而且易于编程,有很高的计算效率．     

The upwind finite difference fractional steps method for combinatorial system of dynamics of fluids in porous media and its application

For combinatorial system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set. Some techniques, such as implicit-explicit difference scheme, calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates, are adopted. Optimal order estimates in L ² norm are derived to determine the error in the second order approximate solution. This method has already been applied to the numerical simulation of migration-accumulation of oil resources. Keywords: combinatorial system, multilayer dynamics of fluids in porous media, two-class upwind finite difference fractional steps method, convergence, numerical simulation of energy sources.     

An adaptive refinement method for convection‐diffusion problems

An adaptive finite difference method for singularly perturbed convection‐diffusion problems is presented. The method is introduced using a first‐order upwind scheme and a suitable error estimator based on the first derivatives. To obtain the grid structure needed for the cross stencil a special refinement strategy is considered. To avoid the slave points we change the stencil at the interface points from a cross to a skew one. After the convergence of the refinement algorithm we use a combination of a first order upwind and a second order central schemes to achieve higher order of convergence. Several numerical examples show the efficiency of our treatment. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Central ADER schemes for hyperbolic conservation laws

In this paper we first briefly review the very high order ADER methods for solving hyperbolic conservation laws. ADER methods use high order polynomial reconstruction of the solution and upwind fluxes as the building block. They use a first order upwind Godunov and the upwind second order weighted average (WAF) fluxes. As well known the upwind methods are more accurate than central schemes. However, the superior accuracy of the ADER upwind schemes comes at a cost, one must solve exactly or approximately the Riemann problems (RP). Conventional Riemann solvers are usually complex and are not available for many hyperbolic problems of practical interest. In this paper we propose to use two central fluxes, instead of upwind fluxes, as the building block in ADER scheme. These are the monotone first order Lax-Friedrich (LXF) and the third order TVD flux. The resulting schemes are called central ADER schemes. Accuracy of the new schemes is established. Numerical implementations of the new schemes are carried out on the scalar conservation laws with a linear flux, nonlinear convex flux and non-convex flux. The results demonstrate that the proposed scheme, with LXF flux, is comparable to those using first and second order upwind fluxes while the scheme, with third order TVD flux, is superior to those using upwind fluxes. When compared with the state of art ADER schemes, our central ADER schemes are faster, more accurate, Riemann solver free, very simple to implement and need less computer memory. A way to extend these schemes to general systems of nonlinear hyperbolic conservation laws in one and two dimensions is presented.     

THE UPWIND FINITE ELEMENT SCHEME AND MAXIMUM PRINCIPLE FOR NONLINEAR CONVECTION-DIFFUSION PROBLEM

In this paper, a kind of partial upwind finite element scheme is studied for twodimensional nonlinear convection-diffusion problem. Nonlinear convection term approximated by partial upwind finite element method considered over a mesh dual to the triangular grid, whereas the nonlinear diffusion term approximated by Galerkin method. A linearized partial upwind finite element scheme and a higher order accuracy scheme are constructed respectively. It is shown that the numerical solutions of these schemes preserve discrete maximum principle. The convergence and error estimate are also given for both schemes under some assumptions. The numerical results show that these partial upwind finite element scheme are feasible and accurate.     

A second order characteristics finite element scheme for natural convection problems

In this paper a second order characteristics finite element scheme is applied to the numerical solution of natural convection problems. Firstly, after recalling the mathematical model, a second order time discretization of the material time derivative is introduced. Next, fully discretized schemes are proposed by using finite element methods. Numerical results for the two-dimensional problem of buoyancy-driven flow in a square cavity with differentially heated side walls are given and compared with a reference solution.     

A new combined finite element‐upwind finite volume method for convection‐dominated diffusion problems

In this article, we develop a combined finite element‐weighted upwind finite volume method for convection‐dominated diffusion problems in two dimensions, which discretizes the diffusion term with the standard finite element scheme, and the convection and source terms with the weighted upwind finite volume scheme. The developed method leads to a totally new scheme for convection‐dominated problems, which overcomes numerical oscillation, avoids numerical dispersion, and has high‐order accuracy. Stability analyses of the scheme are given for the problems with constant coefficients. Numerical experiments are presented to illustrate the stability and optimal convergence of our proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 799–818, 2016