对流扩散方程的迎风变换及相应有限差分方法

本文提出所谓迎风变换,将对流扩散方程分解为对流迎风函数和扩散方程,并构造相应的有限差分格式。对流迎风函数以简明的指数解析形式反映对流扩散现象的迎风效应,原则上消除了源于不对称对流算子的困难,能够便利对流扩散方程的数值求解。有限差分格式具有二阶精度和无条件稳定性,算例表明其准确性、收敛速度及对边界层效应的适应能力均明显优于中心差分格式和迎风差分格式。     

一种基于Associated Hermite正交函数求解对流扩散方程的算法

提出了一种基于AH(Associated Hermite)正交基函数求解对流扩散方程的无条件稳定算法。该算法将方程的时间项通过Hermite多项式作为正交基函数进行展开,利用Galerkin方法消除时间变量项,从而导出有限维AH域隐式差分方程,突破了传统显式差分格式稳定性条件的限制,最后通过对AH域展开系数的求解得到该对流扩散方程的数值解。在数值算例中,将该算法与传统显示差分法和交替方向隐式差分法进行对比分析,数值计算结果表明,算法无条件稳定且其计算精度与时间步长无关,对于具有精细结构的对流换热问题,该算法具有明显的效率优势,且保持了较高的精度。     

数值求解Navier—Stokes方程的迎风紧致差分格式

从迎风紧致逼近^[1]出发,提出数值求解可压Navier-Stokes方程的一种高精度的数值方法。利用Steger-Warming的通量分裂技术^[2]将守恒型方程中的流通向量分裂成两部分,在此基础上据风向构造逼近于无粘项的三阶迎风紧致有限差分格式。对方程中的粘性部分采用通常的二阶差分逼近。所建立的差分格式被用来数值求解了三维粘性绕流问题。     

高精度迎风紧致差分格式与热流计算研究的注记

利用高精度差分格式求解了可压缩 N-S方程球头热流问题。分析了不同差分格式在对球头粘性绕流热流计算中存在的问题 ,并分析了相应的网格雷诺数。在利用高精度迎风紧致 [1 ] 格式求解粘性绕流热流问题时 ,采用 Steger-Warming[2 ]的通量分裂技术将守恒型方程中的流通向量分裂成两部分 ,在此基础上据风向构造逼近于无粘项的高精度迎风格式。对方程中的粘性部分采用中心差分格式。数值结果表明 :高精度差分格式能在较大的网格雷诺数下较好地计算球头驻点热流     

求解浅水波方程的半离散中心迎风方法

将Jin's的界面方法应用到求解双曲守恒型方程的半离散中心迎风方法中,给出了一种新的求解浅水波方程的半离散中心迎风差分方法。对于源项,不是采用传统的单元均值而是采用单元界面处的值来近似,使所得格式对稳定态的求解是均衡的。且已证明所给的二阶精度的求解格式保持水深的非负性,这一特性使其能够较好的处理干河床问题。使用该方法产生的数值粘性(与O(Δ2r-1)同阶)要比交错的中心格式小(与O(Δx2r/Δt)同阶),而且由于数值粘性与时间步长无关,从而时间步长可根据稳定性需要尽可能的小,因此适用于稳定态的求解。     

径向辛体系一种新的差分格式

通过对Hellinger-Reissner变分原理进行坐标变换,将径向模拟为时间,导向辛体系,得到Hamilton对偶方程组.将微分形式的有限差分法引入弹性力学极坐标系下径向辛体系,把对偶方程组中的微分方程直接改用差分方程代替,推导出极坐标系下问题的辛差分方程,从而得到一种全新的径向辛体系差分格式.求解方程组,可直接得到位移和应力.编程计算曲梁等算例,结果表明该辛差分格式是有效的,丰富了弹性力学辛体系差分法的内容.     

方柱绕流的数值模拟

采用有限差分法,对雷诺数为2.2×10~4的方柱绕流进行了大涡模拟(简称LES)。运用时间分裂控制(Split-Operator)法,将N-S方程分为对流步、扩散步和传播步。对Smagorinsky假设在近壁区的发散问题用两层模型进行处理。对流项用迎风—中心差分格式模拟,压力方程用SOR法迭代求解。计算得到的沿对称线的时均顺流向速度与文献上的实验结果进行了比较,结果吻合较好,同时还对绕方柱流的流场结构进行了分析研究。     

基于MUSTA格式的全非线性Boussinesq波浪传播数值模型

建立了求解二维全非线性布氏（Boussinesq）水波方程的有限差分/有限体积混合数值格式. 针对守恒形式的控制方程,采用有限体积方法并结合 MUSTA格式计算数值通量, 剩余项则采用有限差分方法求解, 采用具有总变差减小（totalvariation diminishing, TVD）性质的三阶龙格-库塔法进行时间积分.该格式具备间断捕捉、程序实现简单、数值稳定性强、海岸动边界以及波浪破碎处理方便和可调参数少等优点.利用典型算例对数值模型进行了验证,计算结果与实验数据吻合较好.      

A FULLY NONLINEAR BOUSSINESQ MODELFOR WAVE PROPAGATION BASED ON MUSTA SCHEME

建立了求解二维全非线性布氏（Boussinesq）水波方程的有限差分/有限体积混合数值格式. 针对守恒形式的控制方程,采用有限体积方法并结合 MUSTA格式计算数值通量, 剩余项则采用有限差分方法求解, 采用具有总变差减小（totalvariation diminishing, TVD）性质的三阶龙格-库塔法进行时间积分.该格式具备间断捕捉、程序实现简单、数值稳定性强、海岸动边界以及波浪破碎处理方便和可调参数少等优点.利用典型算例对数值模型进行了验证,计算结果与实验数据吻合较好.     

分块隐式有限差分法计算弯管紊流

本文利用贴体坐标系中分块隐式有限差分法计算矩形截面90°弯管中不可压恒定紊流.在计算中,雷诺方程的数值离散采用混合差分格式,局部联立求解雷诺方程和连续方程而得到速度压力解.在全流场的迭求解过程中采用对称联立Gauss—Seidel法.利用标准K-ε紊流模型模拟紊流.计算结果与有关试验进行了对比.     

A study of upstream-weighted high-order differencing for approximation to flow convection

This paper is concerned with a number of upstream-weighted second- and third-order difference schemes. Also considered are the conventional upwind and central difference schemes for comparison. It commences with a general difference equation which unifies all the given first-, second- and third-order schemes. The various schemes are evaluated through the use of the general equation. The unboundedness and accuracy of the solutions by the difference schemes are assessed via various analyses: examination of the coefficients of the difference equation, Taylor series truncation error analysis, study of the upstream connection to numerical diffusion, single-cell analysis. Finally, the difference schemes are tested on one- and two-dimensional model problems. It is shown that the high-order schemes suffer less from the problem of numerical diffusion than the first-order upwind difference scheme. However, unboundedness cannot be avoided in the solutions by these schemes. Among them the linear upwind difference scheme presents the best compromise between numerical diffusion and solution unboundedness.     

对流项占优问题的MLPG/SUPG方法数值模拟

在计算对流项占优问题时易产生假扩散,本文把流线型迎风格式应用于MLPG方法中可以减少对流项的影响,通过两个典型例子（旋转流场问题和Brezzi问题）验证该格式的精度与有效性,并与文献中的迎风格式的计算结果进行比较,计算结果表明,该方法能有效地克服假扩散现象,有较好的稳定性和较高的计算精度。     

对流扩散方程的摄动有限体积(PFV)方法及讨论

在有限体积(FV)方法的重构近似中,引入数值摄动处理,即把界面数值通量摄动展开成网格间距的幂级数,并利用积分方程自身的性质求出幂级数的系数,同时获得高精度迎风和中心型摄动有限体积(PFV)格式.对标量输运方程给出积分近似为二阶、重构近似为二、三和四阶迎风和中心型PFV格式,这些PFV格式的结构形式及使用基点数与一阶迎风格式完全一致,迎风PFV格式满足对流有界准则;二阶和四阶中心PFV格式对网格Peclet数的任意值均为正型格式,比常用的二阶中心格式优越.用一维标量输运和方腔流动算例说明PFV格式的优良性能,并把PFV方法与性质相近的摄动有限差分(PFD)方法及相关的高精度方法作了对比分析.     

PERTURBATIONAL FINITE VOLUME METHOD FOR THE SOLUTION OF 2-D NAVIER-STOKES EQUATIONS ON UNSTRUCTURED AND STRUCTURED COLOCATED MESHES

IntroductionThefinitevolume (FV)methodusestheintegralformoftheconservationequationasitsstartingpointandcanutilizeconvenientlydiversifiedgrids(structuredandunstructuredgrids)andissuitableforverycomplexgeometry ,whicharewhyitispopularwithengineeringandhasbeenwidelyusedinagreatvarietyofcommercialsoftwareofcomputationalfluiddynamics.Relativetothefiniteelement (FE)methodandthefinitedifferential (FD)method ,thedisadvantageofFVmethodisthatitisnothigheraccuracy .FVmethodisofsecondlevelapproximatio…     

极坐标系下的Legendre谱元方法求解Poisson-型方程

r=0处的坐标奇异性是求解极坐标下Poisson-型方程的关键。本文提出一种极坐标系下基于Galerkin变分的Legendre谱元方法用于求解圆形区域内的Poisson-型方程,物理区域的径向和周向划分若干单元,计算单元均采用Legendre多项式展开;圆心所在单元的径向使用LGR(Legendre Gauss Radau)积分点,其他单元径向使用LGL(Legendre Gauss Lobatto)积分点,从而避免了极点处1/r坐标奇异性,周向单元均采用LGL积分点。利用区域分解技术,可以避免节点在极点附近聚集;最后求解了多个Dirichlet或Neumann边界条件下的Poisson-型方程算例。数值结果表明,谱元方法具有很高的精度。     

Internal three-dimensional viscous flow solutions using the vorticity-potential method

A numerical solution procedure for internal three-dimensional viscous flow is proposed in this paper. The formulation is based on the non-primitive variables, the vorticity and potentials, on a curvilinear grid. A new upwind difference scheme is introduced to overcome the convective instabilities arising in the central difference scheme for the vorticity transport equations, while keeping false diffusion to a minimum level. Developing flows in both straight and curved square ducts are simulated to validate the procedure. The results are compared with both experimental measurements and analytical solutions.     

A vector upstream differencing scheme for problems in fluid flow involving significant source terms in steady-state linear systems

This paper considers a finite difference scheme for modelling the convection/diffusion equation in strongly convective flow regimes including circumstances in which significant source terms are present. The main objective is to provide an alternative approach to central and/or upwind difference methods which for various reasons are unsatisfactory. To illustrate the main features of the scheme, an assessment of its accuracy is made by means of a Taylor expansion analysis and a study of its performance in two model problems. As a demonstration of its generality for use in large-scale practical problems, some numerical results are presented for the prediction of the temperature distribution in a flow through a partially blocked heated rod bundle. The main conclusions are that in almost all practical circumstances results obtained using the scheme are not susceptible to false diffusion or spatial oscillations, which are, respectively, the inherent weaknesses in many upwind and central difference scheme formulations, and in general its use results in improved overall accuracy.     

Effect of spatial discretization schemes on numerical solutions of viscoelastic fluid flows

This study examines the effect of discretization schemes for the convection term in the constitutive equation on numerical solutions of viscoelastic fluid flows. For this purpose, a temporally evolving mixing layer, a two-dimensional vortex pair interacting with a wall, and a fully developed turbulent channel flow are selected as test cases, and eight different discretization schemes are considered. Among them, the first-order upwind difference scheme (UD) and artificial diffusion scheme (AD), which are commonly used in the literature, show most stable and smooth solutions even for highly extensional flows. However, the stress fields are smeared too much by these schemes and the corresponding flow fields are quite different from those obtained by higher-order upwind difference schemes. Among higher-order upwind difference schemes investigated in this study, a third-order compact upwind difference scheme (CUD3) with locally added AD shows stable and most accurate solutions for highly extensional flows even at relatively high Weissenberg numbers.     

A finite element method to solve the Navier-Stokes equations using the method of characteristics

We present a simple and efficient finite element method to solve the Navier-Stokes equations in primitive variables V, p. It uses (a) an explicit advection step, by upwind differencing. Improvement with regard to the classical upwind differencing scheme of the first order is realized by accurate calculation of the characteristic curve across several elements, and higher order interpolation; (b) an implicit diffusion step, avoiding any theoretical limitation on the time increment, and (c) determination of the pressure field by solving the Poisson equation. Two laminar flow calculations are presented and compared to available numerical and experimental results.