对流扩散方程的绝对稳定高阶中心差分格式

将作者提出的数值摄动算法改进为区分离散单元内上游和下游并分别对通量进行高精度重构的双重数值摄动算法,与原(单重)摄动算法相比,双重摄动算法既提高了格式精度又明显扩大了格式的稳定域范围.利用双重摄动算法,即分别利用上游和下游基点变量的摄动重构将高阶流体力学关系及迎风机制耦合进二阶中心格式之中,由此构建了对流扩散方程的对网格Reynolds数的任意值均稳定(绝对稳定)高精度(四阶和八阶精度)三基点中心TVD差分格式,通过解析分析以及3个算例计算证实了构建格式的优良性能;3个算例包括一维线性、非线性(Burgers方程)和二维变系数对流扩散方程.数值计算表明:构建的格式在粗网格下不振荡,构建格式在粗网格时的最大误差L_∞和均方误差L_2与二阶中心格式在细网格时的相应误差一致,对线性方程,构建格式在细网格下可达到L_2精度阶.     

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

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

用物理黏性构建高阶不振荡对流扩散差分格式

利用数值摄动算法, 通过扩散格式数值摄动重构把对流扩散方程的2阶中心差分格式(2-CDS)重构为高精度高分辨率格式, 解析分析和模型方程计算证实了新格式的高精度不振荡性质. 新格式是把物理黏性使流动光滑化的扩散运动规律引入2-CDS 中的结果. 该法显然与构建高级离散格式的常见方法不同. 证实: 数值摄动重构中引入扩散运动规律的结果格式与引入对流运动规律(下游不影响上游的规律)的结果格式一致, 说明对离散方程的数值摄动运算, 在维持原格式结构形式不动的条件下, 不仅能提高格式精度和稳健性, 且可揭示对流离散运动规律与扩散离散运动规律之间的内在关联;同时证实, 文中提出和使用的上、下游分裂方法是构建高精度不振荡离散格式的一个有效方法.     

对流扩散方程QUICK格式的数值摄动高精度重构格式

利用高智提出的数值摄动算法, 把对流扩散方程的常用QUICK格式(黏性和对流项分别用二阶中心和QUICK格式离散)进行了高精度重构, 包括利用离散单元内所有结点的全域重构和分别利用离散单元内上下游结点的上下游重构, 得到两类新的更高阶精度的数值摄动重构格式, 称为高的QUICK格式(G-QUICK格式). G-QUICK格式与QUICK格式相比简单性相当, 但精度更高; 全域重构G-QUICK格式和QUICK格式均为条件稳定, 上下游重构得到一些绝对稳定的G-QUICK格式. 解析分析和数值算例均证实了G-QUICK格式的优良性能, 上下游重构的G-QUICK格式为在对流扩散方程的QUICK格式中避免使用人工黏性提供了新途径.      

同位网格摄动有限体积格式求解浮力驱动方腔流

利用对流扩散方程的摄动有限体积格式,在Rayleigh数从10$^{3}$ 到10$^{8}$的范围内对浮力驱动方腔流动问题作了数值模拟. 对流扩散方程的摄动 有限体积格式具有一阶迎风格式的简洁形式,使用相同的基点,重构近似精度高,特别是两 相邻控制体中心到公共界面的距离相等或不相等,PFV格式公式相同等优点. 在数值模拟中, 无论均匀网格还是非均匀网格均获得与DSC方法、自适应有限元法、多重网格法等Benchmark 解相符较好的数值结果,证明UPFV格式对高Rayleigh数对流传热问题的适用性和有效性.     

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

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

数值耗散对压气机流动分离涡模拟的影响研究

分离涡模拟DES是压气机流动模拟中常用的高保真湍流模式。为了使DES准确解析湍流,数值耗散必须限制在合理范围内。然而,当前的压气机流动DES类研究工作中仍然普遍采用高耗散的迎风格式。本文首先基于DES类方法计算的各向同性衰减湍流结果,定量比较了多种不同数值格式的耗散,证实了高耗散迎风格式严重低估中高波数湍流能量。高阶重构格式可以一定程度上改善该问题,但能量耗散仍然过高。本文在高阶重构的基础上,进一步引入自适应耗散函数修改Riemann求解器,构造了自适应耗散格式。该格式在全波数范围都能准确地预测湍流能谱。将该格式配合DES类方法模拟跨声速离心压气机流动,其预测的压比相比于三阶迎风格式,更加接近实验结果。此外,自适应耗散格式显著提高了中小尺度流动结构的分辨率。分析表明,在使用DES类方法模拟压气机流动时,有必要采用数值耗散较低的离散格式,以准确预测压气机总体性能和流动结构。本文构造的自适应耗散格式是一种良好选择。     

求解多维双曲守恒律方程组的四阶半离散格式

提出了求解多维双曲守恒律方程组的四阶半离散格式。该方法以中心加权基本无振荡(CWENO)重构为基础,同时考虑到在R iemann扇内波传播的局部速度,从而回避了计算过程中的网格交错,建立了数值耗散较小的介于迎风格式和中心格式之间的半离散格式。本文的四阶半离散格式是Kurganov等人的三阶半离散格式的高阶推广。大量的数值算例充分说明了本文方法的高分辨率和稳定性。     

摄动有限差分方法研究进展

振动有限差分（ＰＦＤ）方法,既离散徽商项也离散非微商项（包括微商系数）,在微商用直接差分近似的前提下提高差分格式的精度和分辨率．ＰＦＤ方法包括局部线化微分方程的摄动精确数值解（ＰＥＮＳ）方法和摄动数值解（ＰＮＳ）方法以及考虑非线性近似的摄动高精度差分（ＰＨＤ）方法。论述了这些方法的基本思想、具体技巧、若干方程（对流扩散方程、对流扩散反应方程、双曲方程、抛物方程和ＫｄＶ方程）的ＰＥＮＳ、ＰＮＳ和ＰＨＤ格式,它们的性质及数值实验．并与有关的数值方法作了必要的比较．最后提出值得进一步研究的一些课题．      

求解一维理想磁流体方程的移动网格熵稳定格式

磁流体方程的数值求解在等离子体物理学、天体物理研究以及流动控制等领域具有重要意义,本文构造了用于求解理想磁流体动力学方程的基于移动网格的熵稳定格式,此方法将Roe型熵稳定格式与自适应移动网格算法结合,空间方向采用熵稳定格式对磁流体动力学方程进行离散,利用变分法构造网格演化方程并通过Gauss-Seidel迭代法对其迭代求解实现网格的自适应分布,在此基础上采用守恒型插值公式实现新旧节点上的量值传递,利用三阶强稳定Runge-Kutta方法将数值解推进到下一时间层。数值实验表明,该算法能有效捕捉解的结构(特别是激波和稀疏波),分辨率高,通用性好,具有强鲁棒性。     

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.     

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…     

Developing a physical influence upwind scheme for pressure-based cell-centered finite volume methods

In the framework of a cell-centered finite volume method (FVM), the advection scheme plays the most important role in developing FVMs to solve complicated fluid flow problems for a wide range of Reynolds numbers. Advection schemes have been widely developed for FVMs employing pressure-velocity coupling methodology in the incompressible flow limit. In this regard, the physical influence upwind scheme (PIS) is developed for a cell-centered finite volume coupled solver (FVCS) using a pressure-weighted interpolation method for linking the pressure and velocity fields. The well-known exponential differencing scheme and skew upwind differencing scheme are also deployed in the current FVCS and their numerical results are presented. The accuracy and convergence of the present PIS are evaluated solving flow in a lid-driven square cavity, a lid-driven skewed cavity, and over a backward-facing step (BFS). The flow within the lid-driven square cavity is numerically solved at Reynolds numbers from 400 to 10 000 on a relatively coarse mesh with respect to other reported solutions. The lid-driven skewed cavity test case at Reynolds number of 1000 demonstrates the numerical performance of the present PIS on nonorthogonal grids. The flow over a BFS at Reynolds number of 800 is numerically solved to examine capabilities of current FVCS employing the current PIS in inlet-outlet flow computations. The numerical results obtained by the current PIS are in excellent agreement with those of benchmark solutions of corresponding test cases. Incorporating implicit role of pressure terms in a pressure-weighted interpolation method and development of PIS provides satisfactory solution convergence alongside the numerical accuracy for the current FVCS. A particular numerical verification is performed for the V velocity calculation within the BFS flow field, which confirms the reliability of present PIS.     

Augmented upwind numerical schemes for the groundwater transport advection‐dispersion equation with local operators

When solute transport is advection‐dominated, the advection‐dispersion equation approximates to a hyperbolic‐type partial differential equation, and finite difference and finite element numerical approximation methods become prone to artificial oscillations. The upwind scheme serves to correct these responses to produce a more realistic solution. The upwind scheme is reviewed and then applied to the advection‐dispersion equation with local operators for the first‐order upwinding numerical approximation scheme. The traditional explicit and implicit schemes, as well as the Crank‐Nicolson scheme, are developed and analyzed for numerical stability to form a comparison base. Two new numerical approximation schemes are then proposed, namely, upwind–Crank‐Nicolson scheme, where only for the advection term is applied, and weighted upwind‐downwind scheme. These newly developed schemes are analyzed for numerical stability and compared to the traditional schemes. It was found that an upwind–Crank‐Nicolson scheme is appropriate if the Crank‐Nicolson scheme is only applied to the advection term of the advection‐dispersion equation. Furthermore, the proposed explicit weighted upwind‐downwind finite difference numerical scheme is an improvement on the traditional explicit first‐order upwind scheme, whereas the implicit weighted first‐order upwind‐downwind finite difference numerical scheme is stable under all assumptions when the appropriate weighting factor (θ) is assigned.     

New numerical schemes based on a criterion for constructing essentially stable and accurate numerical schemes for convection-dominated equations

In order to obtain stable and accurate numerical solutions for the convection-dominated steady transport equations, we propose a criterion for constructing numerical schemes for the convection term that the roots of the characteristic equation of the resulting difference equation have poles. By imposing this criterion on the difference coefficients of the convection term, we construct two numerical schemes for the convection-dominated equations. One is based on polynomial differencing and the other on locally exact differencing. The former scheme coincides with the QUICK scheme when the mesh Reynolds number (Rm) is $\mathop \[{\textstyle{{\rm 8} \over {\rm 3}}}\] $, which is the critical value for its stability, while it approaches the second-order upwind scheme as Rm goes to infinity. Hence the former scheme interpolates a stable scheme between the QUICK scheme at Rm = $\mathop \[{\textstyle{{\rm 8} \over {\rm 3}}}\] $ and the second-order upwind scheme at Rm = ∞. Numerical solutions with the present new schemes for the one-dimensional, linear, steady convection-diffusion equations showed good results.     

不可压N-S方程高效算法及二维槽道湍流分析

构造了基于非等距网格的迎风紧致格式,并将其与三阶精度的Adams半隐方法相结合,构造了求解不可压N－S方程高效算法。该算法利用基于交错网格的离散形式的压力Poisson方程求解压力项,解决了边界处的残余散度问题;同时还利用快速Fourier变换将方程的隐式部分解耦,离散后的代数方程组利用追赶法求解,大大减少了计算量。通过对二维槽道流动的数值模拟,证实了所构造的数值方法具有精度高,稳定性好,能抑制混淆误差等优点,同时具有很高的计算效率,是进行壁湍流直接数值模拟的有效方法。在数值模拟的基础上对二维槽道流动进行了分析,得到了Reynolds数从6000到15000的二维流动饱和态解（所谓“二维槽道湍流”）;定性及定量结果均与他人的数值计算结果吻合十分理想。对流场进行了分析,指出了“二维湍流”与三维湍流统计特性的区别。     

群速度控制格式及二维Riemann解

强激波和强接触间断的数值模拟一直是计算流体力学里一个富有挑战性的课题,它们是很多实际流动的基础。三阶迎风紧致格式是一种具有较高分辨率的高精度方法,但是在计算激波时仍有数值振荡产生。本文根据数值解的群速度特性,在三阶迎风紧致格式的基础上提出了一种群速度控制格式,使得能够正确模拟含有强激波和强接触间断的复杂流动。在此基础上构造了求解包含大压力比和密度比的二维界面问题的数值方法。计算结果表明,方法对激波和接触间断的分辨效果是令人满意的。