二维对流反应方程的高精度多重网格方法

利用一阶偏导数项的四阶紧致差分算子,直接推导出了数值求解二维对流反应方程的一种新的高精度紧致差分格式。为了提高差分方程的求解效率,采用多重网格加速技术,建立了与之相适应的多重网格V循环算法。数值实验结果验证了本文方法的精确性和可靠性。     

求解对流扩散方程的紧致修正方法

提出了求解对流扩散方程的紧致修正方法,该方法是在低阶离散格式的源项中,引入紧致修正项,从而构造高阶紧致修正格式,并进行求解.采用紧致修正方法对典型的对流扩散方程进行计算.结果表明,紧致修正方法虽然与二阶经典差分方法建立在相同的结点数上,但紧致修正方法的精度与紧致方法的精度相同,均具有四阶精度.所以紧致修正方法可以在少网...     

非结构网格下非定常流场双时间步长的加权ENO-强紧致杂交高分辨率格式

针对三维非定常、可压缩流场的Navier-Stokes方程组,本文提出一种新的双时间步长高精度快速迭代格式。该格式在时间上具有二阶精度,在空间离散上不低于三阶。在对流项与粘性项的处理上,本格式分别采用了加权ENO-强紧致格式与紧致四阶精度格式的思想。几个典型算例的实践表明:计算结果与相关实验数据比较吻合,初步表明了该算法可以在非结构网格下具有高效率与高分辨率的特征。     

三维定常问题的高阶紧致差分格式向非定常问题的推广

将已经建立的求解三维定常对流扩散方程的高阶紧致差分格式直接推广到三维非定常对流扩散方程的数值求解,时间导数项利用二阶向后欧拉差分公式,所得到的高阶隐式紧致差分格式时间为二阶精度,空间为四阶精度,并且是无条件稳定的.数值实验结果验证了本文方法的精确性和稳健性.     

边界层流动中湍斑的直接数值模拟

用Navier-Stokes方程直接数值模拟平板边界层流动中湍斑的形成和演化过程.发展了模拟湍斑的高精度、高分辨率的高效计算方法,包括推出四阶时间分裂法以提高精度;提出三维耦合差分方法,用于关于压力的泊松方程和关于速度的亥姆霍兹方程的空间离散,建立其四阶三维耦合中心差分格式;并采用四阶紧致迎风差分格式,避免了一般四阶中心差分格式不适用于边界邻域的困难和提高了分辨率;精心地处理各种边界条件,以保持精度和稳定.该方法适用于包含边界邻域的整个区域内的湍斑模拟.通过模拟平板边界层流动中湍斑的复杂演化过程,显示了湍斑的基本特征.     

对流扩散方程的指数型摄动差分法

改进了作者所提出的对流扩散方程四阶指数型摄动差分格式,并阐明其在高Reynolds数适应性和节省计算量方面的显著优点。指数型摄动差分法经改进后具有较为简便的形式,克服了其他紧致高阶格式不能使用于高Reynolds数问题的致命弱点。文中针对计算流体力学的基本困难,作一至三维流动模型方程和自然对流传热问题的精细计算,且以双精制算法检验格式的四阶精度,表明摄动差分法能在较粗的网格下给出相当准确的结果,十分显著地节省计算机时,并对"激波"和"边界层"等高Reynolds数效应有极高的分辨能力。     

非定常对流扩散方程的高精度多重网格方法

由已有的求解定常对流扩散方程的高阶紧致差分格式出发,直接推导出了数值求解非定常对流扩散方程的一种高阶隐式紧致差分格式,其时间为二阶精度,空间为四阶精度,并且是无条件稳定的。为了加快传统迭代法在求解隐格式时在每一个时间步上的迭代收敛速度,采用了多重网格加速技术。数值实验结果验证了本文方法的高阶精度、高效性及高稳定性。     

不同离散格式下的离心泵内部流动数值模拟研究

为研究离散格式对离心泵性能预测精度的影响,本文以自吸式离心泵为计算模型,采用Realizableκ-ε湍流模式进行三维内流场的数值模拟研究,分析了从零流量到最大工作流量下的内部流动和水力性能。建立了考虑内部间隙影响的自吸式离心泵全三维计算模型,分析了动量方程对流项采用一阶差分和二阶差分格式对计算精度的影响,同时分析了压力项的Standard和PRESTO离散格式对计算精度的影响。结果表明,在小流量工况下,采用二阶迎风格式具有较高的计算精度,而在大流量工况下采用一阶迎风格式更为合适。该结果可为准确预测离心泵全工况外特性提供参考依据。     

基于中心差分的对流扩散方程四阶紧凑格式

在经典中心差分格式的基础上,提出对流扩散方程的四阶紧凑差分格式。具体方法是,先就一维情形,将中心差分格式改造为不受网格Reynolds数限制的恒稳二阶格式,再在不增加相关网格点的前提下,通过格式中对流系数和源项的摄动处理,使稳格式的精度提高至四阶。本文并作一、二、三维流动模型方程及高Rayleigh数自然对流传热问题的数值求解,例示本文格式的优良性态。     

热声波数值模拟的虚假振荡研究

采用可压缩流动的SIMPLE算法对一维封闭空腔内由边界突然加热所引起的非稳态热声波进行了数值模拟,对流-扩散项采用了中心差分、一阶迎风差分、QUICK、及MIJSCL等不同格式。计算表明各种格式均存在不同程度的虚假振荡现象,其大小与热声波的强度及离散格式的形式等多种因素有关。这些结果对热声波的进一步研究及高效可靠的对流差分格式的开发具有重要意义。     

二维半线性扩散反应方程的高精度全隐格式及其多重网格方法

针对二维非定常半线性扩散反应方程，空间导数项采用四阶紧致差分公式离散，时间导数项采用四阶向后Euler公式进行离散，提出一种无条件稳定的高精度五层全隐格式.格式截断误差为O（τ⁴+τ²h²+h⁴），即时间和空间均具有四阶精度.对于第一、二、三时间层采用Crank-Nicolson方法进行离散，并采用Richardson外推公式将启动层时间精度外推到四阶.建立适用于该格式的多重网格方法，加快在每个时间层上迭代求解代数方程组的收敛速度，提高计算效率.最后通过数值实验验证格式的精确性和稳定性以及多重网格方法的高效性.     

基于多步高阶差分格式求解时域Maxwell方程

提出基于无穷维哈密尔顿系统及分裂算子理论的多步高阶差分格式,求解时域Maxwell方程.在时间方向上,针对Maxwell方程采用不同阶数的辛算法进行差分离散;在空间方向上,采用四阶差分格式进行差分离散.探讨多步高阶差分格式的稳定性及数值色散性,最后给出数值计算结果.结果表明,五级四阶格式为最有效的多步高阶差分格式,具有高精度、占用较少的计算机资源等优点,适用于长时间的数值模拟.     

The Influence of Inlet Disturbances on the Post-Stall Behaviors in Compression System

The object of this paper is to provide a reliable tool to carry out the parametrical studies of post-stall behaviors in multistage axial compression systems. An adapted version of the 1.5D Euler equations with additional source terms is discretized with a finite volume method and are solved in time by a fourth-order Runge–Kutta scheme. The equations are discretized at mid-span both inside the blade rows and the non-bladed regions. The source terms express the blade-flow interactions and are estimated by calculating the velocity triangles for each blade row. Additional source terms are introduced to represent the effects of inlet disturbances on post-stall behaviors and the physical analysis is therefore proposed to explain the phenomenon.     

三维对流扩散方程非等距网格上的四阶紧致格式及其多重网格方法

提出了数值求解三维变系数对流扩散方程非等距网格上的四阶精度19点紧致差分格式,为了提高求解效率,采用多重网格方法求解高精度格式所形成的大型代数方程组。数值实验结果表明本文方法对于不同的网格雷诺数问题,在精确性、稳定性和减少计算工作量方面均明显优于7点中心差分格式。     

Richtmyer-Meshkov不稳定性强化混合参变机理

在不同参数条件下, 计算分析了H₂O和N₂等混合物界面上激波诱导Richtmyer-Meshkov(R-M)不稳定性过程.采用有限差分方法数值求解了二维可压缩Navier-Stokes方程, 对流项以5阶特征紧致-WENO混合格式离散, 输运项以6阶对称紧致格式离散, 时间方向以3阶显式Runge-Kutta方法推进.研究表明, 界面振幅和激波强度增大, 均可增强界面附近涡量场, 强化混合.      

有厚度平板尾缘可压缩剪切层中的涡结构数值模拟

采用一个新型Fu Ma高精度UCD5 SCD6紧致差分算法,通过直接求解二维Navier Stokes方程,成功实现了有厚度平板尾缘可压缩剪切层中涡结构的数值模拟,并考查了平板厚度对其的影响.计算对流马赫数M_c=0.3,平板厚度分别为1,2,3,4个参考长度.结果表明,增加平板厚度可促使平板尾缘可压缩剪切层中的涡提前卷起,有利于两股气流混合.     

Fully HOC Scheme for Mixed Convection Flow in a Lid-Driven Cavity Filled with a Nanofluid

A fully higher-order compact (HOC) finite difference scheme on the 9-point two-dimensional (2D) stencil is formulated for solving the steady-state laminar mixed convection flow in a lid-driven inclined square enclosure filled with water-$Al_2O_3$ nanofluid. Two cases are considered depending on the direction of temperature gradient imposed (Case I, top and bottom; Case II, left and right). The developed equations are given in terms of the stream function-vorticity formulation and are non-dimensionalized and then solved numerically by a fourth-order accurate compact finite difference method. Unlike other compact solution procedure in literature for this physical configuration, the present method is fully compact and fully higher-order accurate. The fluid flow, heat transfer and heat transport characteristics were illustrated by streamlines, isotherms and averaged Nusselt number. Comparisons with previously published work are performed and found to be in excellent agreement. A parametric study is conducted and a set of graphical results is presented and discussed to elucidate that significant heat transfer enhancement can be obtained due to the presence of nanoparticles and that this is accentuated by inclination of the enclosure at moderate and large Richardson numbers.     

Compact fourth-order finite volume method for numerical solutions of Navier–Stokes equations on staggered grids

The development of a compact fourth-order finite volume method for solutions of the Navier–Stokes equations on staggered grids is presented. A special attention is given to the conservation laws on momentum control volumes. A higher-order divergence-free interpolation for convective velocities is developed which ensures a perfect conservation of mass and momentum on momentum control volumes. Three forms of the nonlinear correction for staggered grids are proposed and studied. The accuracy of each approximation is assessed comparatively in Fourier space. The importance of higher-order approximations of pressure is discussed and numerically demonstrated. Fourth-order accuracy of the complete scheme is illustrated by the doubly-periodic shear layer and the instability of plane-channel flow. The efficiency of the scheme is demonstrated by a grid dependency study of turbulent channel flows by means of direct numerical simulations. The proposed scheme is highly accurate and efficient. At the same level of accuracy, the fourth-order scheme can be ten times faster than the second-order counterpart. This gain in efficiency can be spent on a higher resolution for more accurate solutions at a lower cost.     

A compact finite difference scheme for the fractional sub-diffusion equations

In this paper, a compact finite difference scheme for the fractional sub-diffusion equations is derived. After a transformation of the original problem, the L1 discretization is applied for the time-fractional part and fourth-order accuracy compact approximation for the second-order space derivative. The unique solvability of the difference solution is discussed. The stability and convergence of the finite difference scheme in maximum norm are proved using the energy method, where a new inner product is introduced for the theoretical analysis. The technique is quite novel and different from previous analytical methods. Finally, a numerical example is provided to show the effectiveness and accuracy of the method.     

Multigrid method and fourth-order compact difference discretization scheme with unequal meshsizes for 3D poisson equation

A fourth-order compact difference discretization scheme with unequal meshsizes in different coordinate directions is employed to solve a three-dimensional (3D) Poisson equation on a cubic domain. Two multgrid methods are developed to solve the resulting sparse linear systems. One is to use the full-coarsening multigrid method with plane Gauss–Seidel relaxation, which uses line Gauss–Seidel relaxation to compute each planewise solution. The other is to construct a partial semi-coarsening multigrid method with the traditional point or plane Gauss–Seidel relaxations. Numerical experiments are conducted to test the computed accuracy of the fourth-order compact difference scheme and the computational efficiency of the multigrid methods with the fourth-order compact difference scheme.