固液两相流中的K—ε双方程湍流模式

建立了固液两相流中更一般的Ｋ－ε双方程湍流模式。模化了固相和液相的连续方程、动量方程及Ｋ方程和ε方程，该湍流模型考虑了固液两相间速度的滑移，颗粒间的作用及相间作用。使用本文所建立的湍流模型，数值预测了一管湍流中的沙水混合流动，其预测结果与实验结果比较一致。     

固液两相流中的K－ε双方程湍流模式

建立了固液两相流中更一般的Ｋ－ε双方程湍流模式．模化了固相和液相的连续方程、动量方程及Ｋ方程和ε方程．该湍流模型考虑了固液两相间速度的滑移，颗粒间的作用及相间作用．使用本文所建立的湍流模型，数值预测了一管湍流中的沙水混合流动，其预测结果与实验结果比较一致．     

可压缩均匀各向同性湍流的直接数值模拟

采用8阶精度的中心差分格式及7阶精度的迎风偏斜格式对Re_l = 72~153, M_t= 0.2~0.7的均匀各向同性湍流进行了直接数值模拟, 建立了湍流数据库. 与他人的计算结果吻合十分理想, 说明方法的有效性. 数值结果表明, 采用适当的迎风型差分格式可以克服起动问题(start-up problem)对湍流Mach数的限制, 提高可计算的湍流Mach数, 是可压湍流直接数值模拟的有效方法. 分析了压缩性效应对湍流统计量的影响, 发现压缩性使得湍动能的衰减加快. 探讨了可压湍流中微激波产生的机理, 对流场进行了标度律分析. 发现在本文的Reynolds数和湍流Mach数条件下, 流场中扩展自相似性仍然成立, 同时发现压缩性对标度指数影响不大.     

通过求解输运方程计算壁面距离

壁面距离在当代湍流模化中仍然扮演着关键角色,然而苦于遍历计算壁面距离的高昂代价,该文考虑了求解偏微分方程的途径.基于Eikonal方程构造出类Euler形式的输运方程,这样,可以直接利用求解Euler和Navier-Stokes方程的CFD程序使用的高效数值格式和部分代码.基于北航的MI-CFD(CFD for missles)数值平台,详尽地介绍了该输运方程在直角坐标下的求解过程.使用隐式LIJSGS时间推进和迎风空间离散,发现该方程具有鲁棒快速的收敛特性.为了保证精度,网格度量系数必须也迎风插值计算.讨论了初始条件和边界条件的特殊处理.成功应用该壁面距离求解方法计算了几个含1-1对应网格和重叠网格的复杂外形.     

在两相粘性跨音速喷管流动中薄层方程的一种隐式求解方法

本略去沿流动方向的粘性，将任意曲线坐标系中无量纲化的Ｎ－Ｓ方程简化为薄层方程。采用隐式近似因子分解法求解气相控制方程，采用特征线法跟踪颗粒，然后获得两相跨音速湍流充分耦合的数值方法。其中，颗粒尺寸是分级的，用参考平面中的拟特征线法处理喷管的粘性亚音速进口边界条件，湍流采用代数模型。该计算方法应用于火箭喷管两相粘流计算，并预估了固体火箭发动机的推力和比冲，计算与试验结果吻合很好。中还讨论了不同颗     

颗粒湍能输运方程的两相湍流模型和平面闭式两相射流的数值模拟

在现有的两相湍流数值模拟中,对颗粒湍流普遍采用以局部追随概念为基础的代数模型,其预报结果在很多情况下与实验不符。本文提出了以颗粒湍能输运方程为基础的κ-ε-kκ两相湍流模型,并以平面闭式两相射流为例进行了数值模拟,预报结果与实验符合良好,表明此模型明显优于k-ε-A.P.的颗粒湍流代数模型。     

三维可压缩Boltzmann模型及其在低Mach数湍流中的应用

改进了有限差分格子Boltzmann方法(FDLBM),以直接数值模拟气动噪声.基于LB求解器特性,采用动力学方程中的恒定对流速度以实施高阶迎风差分,提高了声波和湍流的分辨率.通过建立一个新的三维粒子模型,计算得到了任意比热容的三维可压缩Navier-Stokes系统.此外,利用Bhatnagar-Gross-Krook (BGK)碰撞算子,通过引入热流量修正,实现了Prandtl数的可变性.在激波管内弱声波以及伴随有温度梯度的Taylor-Couette层流的验证计算中,提出的新方法结果良好.此外也对NACA0012翼型绕流进行了三维模拟.其中,Reynolds数、Mach数和攻角分别取2× 105,8.75×10-2以及9°.计算发现,在机翼前缘附近的分离气流位置,以及表面压力波动强度的Maeh数依赖性方面,数值计算结果与实验结果相吻合.     

三维可压缩Boltzmann模型及其在低Mach数湍流中的应用（英文）

改进了有限差分格子Boltzmann方法(FDLBM),以直接数值模拟气动噪声.基于LB求解器特性,采用动力学方程中的恒定对流速度以实施高阶迎风差分,提高了声波和湍流的分辨率.通过建立一个新的三维粒子模型,计算得到了任意比热容的三维可压缩Navier-Stokes系统.此外,利用Bhatnagar-Gross-Krook(BGK)碰撞算子,通过引入热流量修正,实现了Prandtl数的可变性.在激波管内弱声波以及伴随有温度梯度的Taylor-Couette层流的验证计算中,提出的新方法结果良好.此外也对NACA0012翼型绕流进行了三维模拟.其中,Reynolds数、Mach数和攻角分别取2×105,8.75×10-2以及9°.计算发现,在机翼前缘附近的分离气流位置,以及表面压力波动强度的Mach数依赖性方面,数值计算结果与实验结果相吻合.     

突然扩张方管中三维湍流流动的数值模拟

本文运用ＳＩＭＰＬＥＣ算法计算了突然扩张方管中的三维湍流流动，湍流模型采用ｋ－ε模型．计算结果详细反映了突然扩张方管中三维湍流流场．从本文结果可以看出，由于突然扩张方管几何形状非轴对称，且尺寸有限，边壁对流场的作用是不可忽略的．以往文献中常见的二维突然扩张湍流的数值模拟结果与三维情况有较大差别，在靠近边壁的区域差别很大，因此对于突然扩张方管中湍流流动的数值模拟应用三维模拟．本文计算所得突然扩张截面后主回流区长度与实验结果接近．本文方法可为数值模拟突然扩张方管中湍流流场及各物理参数的分布提供有效工具．     

超声速流场中6自由度物体运动的模拟研究

物体在流场中自由运动的模拟有很广泛的应用,文章描述计算6自由度(6DOF)刚体在超声速流场中自由运动的一种方法.流体部分求解LES方程,亚网格模型为拉伸涡模型.激波和刚体边界周围区域采用迎风型WENO格式,湍流区域采用低数值耗散的TCD格式.时间推进采用三阶的SSP R-K法.刚体采用6自由度模型,刚体姿态用四元数来表示,控制方程为常微分方程,采用四阶Runge-Kutta法求解.文章给出若干算例来验证程序的有效性,结果理想.     

Turbulent flow around single concentric long capsule in a pipe

A numerical solution was developed for the equations governing the turbulent flow around single concentric long capsule in a pipe. First, a turbulence model was established for the concentric annulus between the capsule and the pipe to simulate the flow as axi-symmetric, two dimensional, steady flow without edge effect. Second, the same case was considered taking into account the edge effect. Finally, turbulence modelling was established to simulate the case as a three dimensional steady flow, with a view of investigating the validity of axi-symmetric flow assumption. Three different turbulence models were used: an algebraic model (Baldwin–Lomax model) and two types of two-equation models (k–ε and k–ω). Obtained results of pressure gradient along the capsule were compared with available experimental data to verify the used models. In addition, experimental data of the velocity profiles of other investigators were also used in this concern. The results predicted by the three different turbulence models were shown to agree well with the experimental data, though precision differed from one to another.     

Symbolic computation of high order compact difference schemes for three dimensional linear elliptic partial differential equations with variable coefficients

We present a symbolic computation procedure for deriving various high order compact difference approximation schemes for certain three dimensional linear elliptic partial differential equations with variable coefficients. Based on the Maple software package, we approximate the leading terms in the truncation error of the Taylor series expansion of the governing equation and obtain a 19 point fourth order compact difference scheme for a general linear elliptic partial differential equation. A test problem is solved numerically to validate the derived fourth order compact difference scheme. This symbolic derivation method is simple and can be easily used to derive high order difference approximation schemes for other similar linear elliptic partial differential equations.     

动脉血管流动计算的伽辽金有限元法研究

得到大动脉三维模型的过二重分叉的二维截定常流的NS方程有限元解,采用了物理坐标系统换到曲线边界贴休坐标系的数学技巧,以支流至主动脉流率为参数,计算了雷诺数为1000的壁面切应力,所得结果与前人的工作（包括实验数据）进行了比较,发现与他们的结果非常接近,改进了Sharma和Kapoor(1995)的工作,相比之下,所用的数值方法上更经济,适用的雷诺数更大。     

Laguerre polynomial approach for solving linear delay difference equations

This paper presents a numerical method for the approximate solution of mth-order linear delay difference equations with variable coefficients under the mixed conditions in terms of Laguerre polynomials. The aim of this article is to present an efficient numerical procedure for solving mth-order linear delay difference equations with variable coefficients. Our method depends mainly on a Laguerre series expansion approach. This method transforms linear delay difference equations and the given conditions into matrix equation which corresponds to a system of linear algebraic equation. The reliability and efficiency of the proposed scheme are demonstrated by some numerical experiments and performed on the computer algebraic system Maple.     

Numerical scheme for solving system of fractional partial differential equations with Volterra‐type integral term through two‐dimensional block‐pulse functions

In the current study, an approximate scheme is established for solving the fractional partial differential equations (FPDEs) with Volterra integral terms via two‐dimensional block‐pulse functions (2D‐BPFs). According to the definitions and properties of 2D‐BPFs, the original problem is transformed into a system of linear algebra equations. By dispersing the unknown variables for these algebraic equations, the numerical solutions can be obtained. Besides, the proof of the convergence of this system is given. Finally, several numerical experiments are presented to test the feasibility and effectiveness of the proposed method.     

A high order finite difference method with Richardsonextrapolation for 3D convection diffusion equation

In this paper, we extend the Sun and Zhang’s [24] work on high order finite difference method, which is based on the Richardson extrapolation technique and an operator interpolation scheme for the one and two dimensional steady convection diffusion equations to the three dimensional case. Firstly, we employ a fourth order compact difference scheme to get the fourth order accurate solution on the fine and the coarse grids. Then, we use the Richardson extrapolation technique by combining the two approximate solutions to get a sixth order accurate solution on coarse grid. Finally, we apply an operator interpolation scheme to achieve the sixth order accurate solution on the fine grid. During this process, we use alternating direction implicit (ADI) method to solve the resulting linear systems. Numerical experiments are conducted to verify the accuracy and effectiveness of the present method.     

Singly diagonally implicit runge‐kutta method for time‐dependent reaction‐diffusion equation

In this article, an efficient fourth‐order accurate numerical method based on Padé approximation in space and singly diagonally implicit Runge‐Kutta method in time is proposed to solve the time‐dependent one‐dimensional reaction‐diffusion equation. In this scheme, we first approximate the spatial derivative using the second‐order central finite difference then improve it to fourth‐order by applying Padé approximation. A three stage fourth‐order singly diagonally implicit Runge‐Kutta method is then used to solve the resulting system of ordinary differential equations. It is also shown that the scheme is unconditionally stable, and is suitable for stiff problems. Several numerical examples are solved by the scheme and the efficiency and accuracy of the new scheme are compared with two widely used high‐order compact finite difference methods. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1423–1441, 2011     

Numerical simulation of axisymmetric afterbody flows with jet exhaust

A fifth-order accurate compact difference scheme was used to compute the flow over an axisymmetric afterbody with jet exhaust. The solution was based on the mass-averaged Navier-Stokes equations combined with a two-parameter differential model of turbulence. The computations were performed on a specially generated mesh such that the flow in the exterior and interior of the nozzle could be described simultaneously. Numerical results are presented for various external flow conditions and various chamber pressures.     

An extension ofA-stability to alternating direction implicit methods

Completely implicit, noniterative, finite-difference schemes have recently been developed by several authors for nonlinear, multidimensional systems of hyperbolic and mixed hyperbolic-parabolic partial differential equations. The method of Douglas and Gunn or the method of approximate factorization can be used to reduce the computational problem to a sequence of one-dimensional or alternating direction implicit (ADI) steps. Since the eigenvalues of partial differential equations (for example, the equations of compressible fluid dynamics) are often widely distributed with large imaginary parts,A-stable integration formulas provide ideal time-differencing approximations. In this paper it is shown that if anA-stable linear multistep method is used to integrate a model two-dimensional hyperbolic-parabolic partial differential equation, then one can always construct an ADI scheme by the method of approximate factorization which is alsoA-stable, i.e., unconditionally stable. A more restrictive result is given for three spatial dimensions. Since necessary and sufficient conditions forA-stability can easily be determined by using the theory of positive real functions, the stability analysis of the factored partial difference equations is reduced to a simple algebraic test.The main results of this paper were presented at the SIAM National Meeting, Madison, Wis., May 24 to 26, 1978, and section 9 was part of a presentation at the 751st Meeting of the American Mathematical Society, San Luis Obispo, California, Nov. 11 to 12, 1977.     

Numerical treatment for solving fractional Riccati differential equation

This paper presents an accurate numerical method for solving fractional Riccati differential equation (FRDE). The proposed method so called fractional Chebyshev finite difference method (FCheb-FDM). In this technique, we approximate FRDE with a finite dimensional problem. The method is based on the combination of the useful properties of Chebyshev polynomials approximation and finite difference method. The Caputo fractional derivative is replaced by a difference quotient and the integral by a finite sum. By this method the given problem is reduced to a problem for solving a system of algebraic equations, and by solving this system, we obtain the solution of FRDE. Special attention is given to study the convergence analysis and estimate an error upper bound of the obtained approximate formula. Illustrative examples are included to demonstrate the validity and applicability of the proposed technique.