三维对流扩散方程的高精度全隐式多重网格方法

提出了数值求解三维非定常变系数对流扩散方程的一种高精度全隐紧致差分格式,该格式在空间上具有四阶精度,时间具有二阶精度。为了克服传统迭代法在每一个时间步上迭代求解隐格式时收敛速度慢的缺点,采用多重网格加速技术,建立了适用于本文高精度全隐紧致格式的多重网格算法,从而大大加快了迭代收敛速度。数值实验结果验证了本文方法的精确性、稳定性和对高网格雷诺数问题的强适应性。     

全机绕流Euler方程多重网格分区计算方法

全机三维复杂形状绕流数值求解只能采用分区求解的方法,本文采用可压缩Euler方程有限体积方法以及多重网格分区方法对流场进行分区计算。数值方法采用改进的van Leer迎风型矢通量分裂格式和MUSCL方法,基于有限体积方法和迎风型矢通量分裂方法,建立一套处理子区域内分界面的耦合条件。各个子区域之间采用显式耦合条件,区域内部采用隐式格式和局部时间步长等,以加快收敛速度。计算结果飞机表面压力分布等气动力特性与实验值进行了比较,二者基本吻合。计算结果表明采用分析“V”型多重网格方法,能提高计算效率,加快收敛速度达到接近一个量级。根据全机数值计算结果和可视化结果讨论了流场背风区域旋涡的形成过程。     

一种新的求解对流扩散边值问题的无网格方法

基于配点法和楔形基函数,提出了一种新的求解对流扩散边值问题的无网格方法。通过一维和二维的问题验证了该数值方法的可行性;并根据数值算例和分析,可以看到该数值方法能达到满意的收敛效果。该数值方法的隐格式形式能够有效地消除对流占优问题的数值振荡现象,是一种真正的无网格方法。     

一种基于块雅可比迭代的高阶FR格式隐式方法

最近, 基于非结构网格的高阶通量重构格式(flux reconstruction, FR)因其构造简单且通用性强而受到越来越多人的关注. 但将FR格式应用于大规模复杂流动的模拟时仍面临计算开销大、求解时间长等问题. 因此, 亟需发展与之相适应的高效隐式求解方法和并行计算技术. 本文提出了一种基于块Jacobi迭代的高阶FR格式求解定常二维欧拉方程的单GPU隐式时间推进方法. 由于直接求解FR格式空间和隐式时间离散后的全局线性方程组效率低下并且内存占用很大. 而通过块雅可比迭代的方式, 能够改变全局线性方程组左端矩阵的特征, 克服影响求解并行性的相邻单元依赖问题, 使得只需要存储和计算对角块矩阵. 最终将求解全局线性方程组转化为求解一系列局部单元线性方程组, 进而又可利用LU分解法在GPU上并行求解这些小型局部线性方程组. 通过二维无黏Bump流动和NACA0012无黏绕流两个数值实验表明, 该隐式方法计算收敛所用的迭代步数和计算时间均远小于使用多重网格加速的显式Runge-Kutta格式, 且在计算效率方面至少有一个量级的提升.      

基于高阶和高效格式的悬停旋翼可压缩无粘绕流的计算

发展了一种基于高精度和高效格式计算悬停旋翼复杂绕流的隐式有限体积方法。控制方程为Euler方程,其中对流项通量的左右状态量采用五阶加权基本无振荡(WENO)格式重构,时间推进应用隐式LU-SGS算法,为进一步加速收敛,采用三层V循环多重网格松弛方法。考虑到多重网格方法的思想以及五阶WENO格式涉及6个网格单元,建议仅在最细网格上使用WENO格式。计算结果表明本文方法能有效捕捉激波,对尾迹也有较高分辨率,基于隐式LU-SGS算法的多重网格方法能有效提高计算效率。     

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

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

非结构网格上p型多重网格法流场数值模拟

在二维、三维非结构网榕上,针对间断Galerkin方法计算量大、收敛慢的缺点将p型多重网格方法应用于该方法求解跨音速Euler方程,提高计算效率。p型多重网格方法是通过对不同阶次多项式近似解进行递归迭代求解,来达到加速收敛。文中对高阶近似(p＞0)使用显式格式,最低阶近似(p=0)采用隐式格式。NACA0012翼型和O...     

应用于计算气动声学的优化有限紧致格式

对于含间断的计算气动声学问题,数值计算的格式不仅要求低耗散低色散的设计,对短波具有较高的分辨率,还要求能捕捉激波.中心紧致格式具有高精度,具有无耗散和低色散特征,但不能捕捉间断和激波;WENO格式处理间断较为成功,而耗散和色散误差相对较大.有限紧致格式可以将紧致格式与WENO格式相结合构造成混合格式,利用光滑因子之间的关系对激波区域进行自动判断,将传统的全域求解的紧致格式划分为有限的局部紧致求解,间断点上的激波捕捉铜梁自动作为局部紧致求解的边界通量,在在光滑区域具有紧致格式的高精度低耗散性质,在激波附近不产生非物理振荡.本文利用有限紧致格式思想,构造了新的适合于气动声学问题的优化有限紧致格式,将其应用于计算气动声学一维标准测试问题,对相关格式的模拟性能进行了评估,显示该格式在宽频声波传播和含有间断的声波传播模拟方面具有优势.     

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

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

二维振荡叶栅非定常粘性流动数值模拟

采用显式四步Runge-Kutta格式，结合Baldwon-Lomax紊流模型求解Navier-Stokes方程，借助运动网格技术，完成了对二维振荡叶栅非定常粘性流动的数值模拟。为了加速求解过程，引入了变系数隐式残差光顺方法，取得了较好效果。数值结果与已公布的数据有很好的一致性。     

Alternating segment explicit-implicit scheme for nonlinear third-order KdV equation

A group of asymmetric difference schemes to approach the Korteweg-de Vries(KdV)equation is given here.According to such schemes,the full explicit difference scheme and the fun implicit one,an alternating segment explicit-implicit difference scheme for solving the KdV equation is constructed.The scheme is linear unconditionally stable by the analysis of linearization procedure,and is used directly on the parallel computer. The numerical experiments show that the method has high accuracy.     

Improvements on the accuracy of derivative estimation from DPIV velocity measurements

Accuracy of out-of-plane vorticity estimation from in-plane experimental velocity measurements is investigated with particular application to digital particle image velocimetry (DPIV). Simulations of known flow fields are used to quantify errors associated with amplification of the velocity measurement noise and method bias error due to spatial sampling resolution. A novel, adaptable, hybrid estimation scheme combining implicit compact finite difference and Richardson extrapolation schemes is proposed for improved vorticity estimation. The scheme delivers higher-order truncation error with less noise amplification than an explicit second order finite difference scheme. Finally, a complete framework for predicting, a priori, the random, bias, and total error of the vorticity estimation on the basis of the error of the resolved velocities and the choice of differentiation scheme is developed and presented. A portion of this work was presented at ASME IMECE 2003 conference An erratum to this article is available at .     

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.     

High-order numerical methods of fractional-order Stokes’ first problem for heated generalized second grade fluid

The high-order implicit finite difference schemes for solving the fractionalorder Stokes’ first problem for a heated generalized second grade fluid with the Dirichlet boundary condition and the initial condition are given. The stability, solvability, and convergence of the numerical scheme are discussed via the Fourier analysis and the matrix analysis methods. An improved implicit scheme is also obtained. Finally, two numerical examples are given to demonstrate the effectiveness of the mentioned schemes.     

Development of an Artificial Compressibility Methodology with Implicit LU-SGS Method

A numerical method has been developed to solve the steady and unsteady incompressible Navier-Stokes equations in a two-dimensional, curvilinear coordinate system. The solution procedure is based on the method of artificial compressibility and uses a third-order flux-difference splitting upwind differencing scheme for convective terms and second-order center difference for viscous terms. A time-accurate scheme for unsteady incompressible flows is achieved by using an implicit real time discretization and a dual-time approach, which introduces pseudo-unsteady terms into both the mass conservation equation and momentum equations. An efficient fully implicit algorithm LU-SGS, which was originally derived for the compressible Eulur and Navier-Stokes equations by Jameson and Toon [1], is developed for the pseudo-compressibility formulation of the two dimensional incompressible Navier-Stokes equations for both steady and unsteady flows. A variety of computed results are presented to validate the present scheme. Numerical solutions for steady flow in a square lid-driven cavity and over a backward facing step and for unsteady flow in a square driven cavity with an oscillating lid and in a circular tube with a smooth expansion are respectively presented and compared with experimental data or other numerical results.     

基于全隐式无分裂算法求解三维N-S方程

基于多块结构网格,本文研究和发展了三维N-S方程的全隐式无分裂算法.对流项的离散运用Roe格式,粘性项的离散利用中心型格式.在每一次隐式时间迭代中,运用GMRES方法直接求解隐式离散引起的大型稀疏线性方组.为了降低内存需求以及矩阵与向量之间的运算操作数,Jacobian矩阵的一种逼近方法被应用在本文的算法之中.计算结果与实验结果基本吻合,表明本文的全隐式无分裂方法是有效和可行的.     

Operator-splitting methods for the incompressible Navier-Stokes equations on non-staggered grids. Part 1: First-order schemes

New implicit finite difference schemes for solving the time-dependent incompressible Navier-Stokes equations using primitive variables and non-staggered grids are presented in this paper. A priori estimates for the discrete solution of the methods are obtained. Employing the operator approach, some requirements on the difference operators of the scheme are formulated in order to derive a scheme which is essentially consistent with the initial differential equations. The operators of the scheme inherit the fundamental properties of the corresponding differential operators and this allows a priori estimates for the discrete solution to be obtained. The estimate is similar to the corresponding one for the solution of the differential problem and guarantees boundedness of the solution. To derive the consistent scheme, special approximations for convective terms and div and grad operators are employed. Two variants of time discretization by the operator-splitting technique are considered and compared. It is shown that the derived scheme has a very weak restriction on the time step size. A lid-driven cavity flow has been predicted to examine the stability and accuracy of the schemes for Reynolds number up to 3200 on the sequence of grids with 21 × 21, 41 × 41, 81 × 81 and 161 × 161 grid points.     

A dual-porosity model for gas reservoir flow incorporating adsorption behaviour—Part II. Numerical algorithm and example applications

In the present study, an algorithm is presented for the dual-porosity model formulated in Part I of this series. The resultant flow equation with the dual-porosity formulation is of an integro-(partial) differential equation involving differential terms for the Darcy flow in large fractures and integrals in time for diffusion within matrix blocks. The algorithm developed here to solve this equation involves a step-by-step finite difference procedure combined with a quadrature scheme. The quadrature scheme, used for the integral terms, is based on the trapezoidal method which is of second-order precision. This order of accuracy is consistent with the first- and second-order finite difference approximations used here to solve the differential terms in the derived flow equation. In an approach consistent with many petroleum reservoir and groundwater numerical flow models, the example formulation presented uses a first-order implicit algorithm. A two-dimensional example is also demonstrated, with the proposed model and numerical scheme being directly incorporated into the commercial gas reservoir simulator SIMED II that is based on a fully implicit finite difference approach. The solution procedure is applied to several problems to demonstrate its performance. Results from the derived dual-porosity formulation are also compared to the classic Warren–Root model. Whilst some of this work confirmed previous findings regarding Warren–Root inaccuracies at early times, it was also found that inaccuracy can re-enter the Warren–Root results whenever there are changes in boundary conditions leading to transient variation within the domain.     

Implementing a high‐order accurate implicit operator scheme for solving steady incompressible viscous flows using artificial compressibility method

This paper uses a fourth‐order compact finite‐difference scheme for solving steady incompressible flows. The high‐order compact method applied is an alternating direction implicit operator scheme, which has been used by Ekaterinaris for computing two‐dimensional compressible flows. Herein, this numerical scheme is efficiently implemented to solve the incompressible Navier–Stokes equations in the primitive variables formulation using the artificial compressibility method. For space discretizing the convective fluxes, fourth‐order centered spatial accuracy of the implicit operators is efficiently obtained by performing compact space differentiation in which the method uses block‐tridiagonal matrix inversions. To stabilize the numerical solution, numerical dissipation terms and/or filters are used. In this study, the high‐order compact implicit operator scheme is also extended for computing three‐dimensional incompressible flows. The accuracy and efficiency of this high‐order compact method are demonstrated for different incompressible flow problems. A sensitivity study is also conducted to evaluate the effects of grid resolution and pseudocompressibility parameter on accuracy and convergence rate of the solution. The effects of filtering and numerical dissipation on the solution are also investigated. Test cases considered herein for validating the results are incompressible flows in a 2‐D backward facing step, a 2‐D cavity and a 3‐D cavity at different flow conditions. Results obtained for these cases are in good agreement with the available numerical and experimental results. The study shows that the scheme is robust, efficient and accurate for solving incompressible flow problems. Copyright © 2010 John Wiley & Sons, Ltd.