基于同位网格下求解N-S方程的快速算法

在有限容积法基础上建立了基于同位网格的SIMPLEM算法。此算法使初始压力场与速度场耦合,让压力场和速度场同时更好地满足动量方程和连续性方程,且兼顾考虑扩散对流项对计算节点速度修正值的影响及源项与速度场之间的同步性,详细给出了算法的推导过程且对方腔顶盖驱动流进行了数值模拟。计算节点的布置采用同位网格技术,界面流速通过动量插值确定,在不同条件下讨论了迭代次数与残差的关系和不同算法的收敛性,同时验证了算法及程序是准确和可信的。     

流固耦合分析的一种改进CBS有限元算法

针对流固耦合问题, 发展了一种基于任意拉格朗日-欧拉(ALE)描述有限元法的弱耦合分区算法. 运用半隐式特征线分裂算法求解Navier-Stokes方程, 在压力Poisson 方程中引入质量源项以满足几何守恒律; 运用子块移动技术更新动态网格, 并配以光滑处理防止网格质量下降; 采用Newmark-β 法求解结构运动方程. 为保持流体-结构界面处速度和动量守恒, 利用修正结合界面边界条件方法求解界面处速度通量和动量通量. 运用本方法分别模拟了不同雷诺数下单圆柱横向和两向流致振动、串列双圆柱两向流致振动. 计算表明, 本文方法计算效率高, 计算结果与已有实验和数值计算数据吻合.     

不同网格扭曲率下压力修正全隐算法——IDEAL求解性能研究

2008年,本文作者和陶文铨等提出了一种用于速度和压力耦合求解的高效稳定压力修正全隐算法IDEAL,该算法通过在每个迭代层次上对压力方程进行两次内迭代计算,完全克服了SIMPLE算法的两个假设,充分满足了速度和压力之间的耦合,从而大大提高了计算的收敛性和健壮性。为了进一步实现IDEAL算法的推广应用,本文基于三维倾斜方腔顶盖驱动流动,研究了IDEAL算法在不同网格扭曲率下的求解特性。研究发现,在不同网格扭曲率下,IDEAL算法的健壮性和收敛性均优于SIMPLE算法,特别在高网格扭曲率情况下,IDEAL算法求解性能更加优于SIMPLE算法。在不同网格扭曲率下,IDEAL算法健壮性保持不变,几乎可以在任意速度亚松弛因子下获得收敛的解,同时IDEAL算法最短计算耗时较SIMPLE算法减少了56%~89%,验证了IDEAL算法的优越性。     

不同网格扭曲率下压力修正全隐算法——IDEAL求解性能研究

2008年,本文作者和陶文铨等提出了一种用于速度和压力耦合求解的高效稳定压力修正全隐算法IDEAL,该算法通过在每个迭代层次上对压力方程进行两次内迭代计算,完全克服了SIMPLE算法的两个假设,充分满足了速度和压力之间的耦合,从而大大提高了计算的收敛性和健壮性.为了进一步实现IDEAL算法的推广应用,本文基于三维倾斜方腔顶盖驱动流动,研究了IDEAL算法在不同网格扭曲率下的求解特性.研究发现,在不同网格扭曲率下,IDEAL算法的健壮性和收敛性均优于SIMPLE算法,特别在高网格扭曲率情况下,IDEAL算法求解性能更加优于SIMPLE算法.在不同网格扭曲率下,IDEAL算法健壮性保持不变,几乎可以在任意速度亚松弛因子下获得收敛的解,同时IDEAL算法最短计算耗时较SIMPLE算法减少了56%~89%,验证了IDEAL算法的优越性.     

航行体水下发射流固耦合效应分析

对于水下发射过程来说,掌握水动力载荷形成机理与结构响应特征是一个亟待解决的问题.研究该问题需要考虑含相变的复杂多相流动,变约束的结构运动以及这二者之间的耦合效应.本文采用松耦合的方法,以流体求解器为主体,将自编的固体结构程序接入流体求解器中,在每个时间步长内分别对流体动力学方程和固体结构动力学方程进行求解,通过流固界面之间的数据交换实现耦合计算.其中,流体求解器基于雷诺平均纳维斯托克斯方程,采用单流体模型处理多相流问题,引入空化模型描述空化相变,采用修正的湍流模型模拟混合物的湍流效应,并采用动网格技术处理移动边界问题.航行体的刚体运动和结构振动分开求解.结构求解器采用等效梁模型描述结构的振动,通过坐标变换给出了随体坐标系下的结构振动方程,求解方法采用时域积分法.所建立的流固耦合方法不仅能够捕捉到自然空化的演化情况,还可获得航行体所受水动力、结构振动响应以及截面的弯矩,获得了实验的验证.基于该方法研究了结构刚度、发射速度对空泡溃灭与结构振动耦合效应的影响规律.结果表明,同步溃灭是影响结构载荷的主要因素,包括溃灭压力幅值,溃灭压力作用位置,以及溃灭压力与结构振动的相位关系.     

二阶矩两相湍流模型中的改进算法

建立了融合SIMPLEC算法在内的考虑了浓度修正影响的颗粒压力修正方程。提出了二阶矩湍流模型下考虑浓度修正值影响的两相湍流流动的算法,并将它和没有考虑浓度修正值影响的二阶矩湍流模型进行了对照。结果表明在二阶矩模型中是否考虑浓度修正影响会不同程度地影响流场的速度以及浓度等参数分布,考虑浓度修正影响的二阶矩湍流模型更能有效地预测稠密两相湍流流动。     

多介质可压缩流体动力学界面捕捉方法

研究多介质流界面捕捉方法的主要目的是消除多介质流体在界面处压力、速度可能出现的非物理振荡现象 ,并通过流体动力学方程和界面捕捉方程的耦合 ,将多介质流体动力学计算形式上转化为单介质流体计算 ,从而可以采用对计算单介质有效的高精度计算方法来处理多介质流动问题。推广了Shyue界面捕捉和其等效方程的推导方法 ,给出的结果可以适用于具有状态方程 p =( ,e,a1 , ,an) +( ,e,b1 , ,bn)e的介质 ,并通过了数值试验验证。     

基于Boltzmann模型方程的气体运动论统一算法研究

模型方程出发,研究确立含流态控制参数可描述不同流域气体流动特征的气体分子速度分布函数方程; 研究发展气体运动论离散速度坐标法, 借助非定常时间分裂数值计算方法和NND差分格式, 结合DSMC方法关于分子运动与碰撞去耦技术, 发展直接求解速度分布函数的气体运动论耦合迭代数值格式; 研制可用于物理空间各点宏观流动取矩的离散速度数值积分方法, 由此提出一套能有效模拟稀薄流到连续流不同流域气体流动问题统一算法. 通过对不同Knudsen数下一维激波内流动、二维圆柱、三维球体绕流数值计算表明, 计算结果与有关实验数据及其它途径研究结果(如DSMC模拟值、N-S数值解)吻合较好, 证实气体运动论统一算法求解各流域气体流动问题的可行性. 尝试将统一算法进行HPF并行化程序设计, 基于对球体绕流及类``神舟'返回舱外形绕流问题进行HPF初步并行试算, 显示出统一算法具有很好的并行可扩展性, 可望建立起新型的能有效模拟各流域飞行器绕流HPF并行算法研究方向. 通过将气体运动论统一算法推广应用于微槽道流动计算研究, 已初步发展起可靠模拟二维短微槽道流动数值算法; 通过对Couette流、Poiseuille流、压力驱动的二维短槽道流数值模拟, 证实该算法对微槽道气体流动问题具有较强的模拟能力, 可望发展起基于Boltzmann模型方程能可靠模拟MEMS微流动问题气体运动论数值计算方法研究途径.      

基于贴体网格的VOF方法数模流场研究

提出了一种基于VOF方法的模拟具有复杂边界形状结构物附近流场的新算法,BFC—SIMPLE—VOF算法。采用坐标变换方法实现了任意复杂区域的结构化网格划分,在贴体网格下对二维不可压缩粘性流体的控制方程进行了离散。提出了基于交错网格的修正SIMPLE算法来迭代求解压力一速度场,修正了贴体坐标下的界面跟踪方法（VOF方法）...     

多介质流体非守恒律欧拉方程组的数值计算方法

对多介质流体在界面处满足的Euler方程进行了探讨 ,方程组中增加了描述材料参数间断性质的对流形式非守恒律方程组。以波传播算法为基础 ,通过Roe方程近似求解Riemann问题 ,同时采用相同的数值差分格式求解流体动力学Euler方程组和界面方程组。该方法可以有效消除多介质流体在界面处压力、速度可能出现的非物理振荡。给出了部分典型一维和二维数值计算结果。     

Simulation of planar converging flow of a Leonov viscoelastic fluid

An efficient finite element algorithm is presented to simulate the planar converging flow for the viscoelastic fluid of the Leonov model. The governing equation set, composed of the continuity, momentum and constitutive equations for the Leonov fluid flow, is conveniently decoupled and a two-stage cyclic iteration technique is employed to solve the velocity and elastic strain fields separately. Artificial viscosity terms are imposed on the momentum equations to relax the elastic force and data smoothing is performed on the iterative calculations for velocities to further stabilize the numerical computations. The calculated stresses agree qualitatively with the experimental measurements and other numerically simulated results available in the literature. Computations were successful to moderately high values of Deborah number of about 27·5.     

A robust algorithm for computing fluid flows on highly non‐smooth staggered grids

A new method for computing the fluid flow in complex geometries using highly non‐smooth and non‐orthogonal staggered grid is presented. In a context of the SIMPLE algorithm, pressure and physical tangential velocity components are used as dependent variables in momentum equations. To reduce the sensitivity of the curvature terms in response to coordinate line orientation change, these terms are exclusively computed using Cartesian velocity components in momentum equations. The method is then used to solve some fairly complicated 2‐D and 3‐D flow field using highly non‐smooth grids. The accuracy of results on rough grids (with sharp grid line orientation change and non‐uniformity) was found to be high and the agreement with previous experimental and numerical results was quite good. Copyright © 2008 John Wiley & Sons, Ltd.     

A promising boundary element formulation for three‐dimensional viscous flow

In this paper, a new set of boundary‐domain integral equations is derived from the continuity and momentum equations for three‐dimensional viscous flows. The primary variables involved in these integral equations are velocity, traction, and pressure. The final system of equations entering the iteration procedure only involves velocities and tractions as unknowns. In the use of the continuity equation, a complex‐variable technique is used to compute the divergence of velocity for internal points, while the traction‐recovery method is adopted for boundary points. Although the derived equations are valid for steady, unsteady, compressible, and incompressible problems, the numerical implementation is only focused on steady incompressible flows. Two commonly cited numerical examples and one practical pipe flow problem are presented to validate the derived equations. Copyright © 2004 John Wiley & Sons, Ltd.     

Justification of the drift-flux model for two-phase flow in a circular pipe

A flow of a gas-liquid dispersed mixture in a circular pipe with a variable inclination to the horizon, as applied to oil and gas flows in wells, is considered. Within the framework of a multi- fluid approach, the equations of an asymptotic drift-flux model, which contains an algebraic relation between the phase velocities and one momentum equation for the volume-averaged velocity of the mixture, are derived. It is shown that the drift-flux model in this formulation strictly follows from the balance laws under assumption of inertialess velocity slip of the phases in case of validity of one of the following conditions: (i) the dispersed-phase volume fraction is small; (ii) the phase velocity slip may be neglected; or (iii) the flow regime is inertialess and the acceleration of the mixture can be neglected. A numerical algorithm based on the SIMPLE method is implemented for solving the obtained equations of the drift-flux model. The possibility of modeling the gravitational segregation and the pressure buildup in a shut-in well and transient slug flows is demonstrated.     

Finite volume multigrid method of the planar contraction flow of a viscoelastic fluid

This paper reports on a numerical algorithm for the steady flow of viscoelastic fluid. The conservative and constitutive equations are solved using the finite volume method (FVM) with a hybrid scheme for the velocities and first‐order upwind approximation for the viscoelastic stress. A non‐uniform staggered grid system is used. The iterative SIMPLE algorithm is employed to relax the coupled momentum and continuity equations. The non‐linear algebraic equations over the flow domain are solved iteratively by the symmetrical coupled Gauss–Seidel (SCGS) method. In both, the full approximation storage (FAS) multigrid algorithm is used. An Oldroyd‐B fluid model was selected for the calculation. Results are reported for planar 4:1 abrupt contraction at various Weissenberg numbers. The solutions are found to be stable and smooth. The solutions show that at high Weissenberg number the domain must be long enough. The convergence of the method has been verified with grid refinement. All the calculations have been performed on a PC equipped with a Pentium III processor at 550 MHz. Copyright © 2001 John Wiley & Sons, Ltd.     

基于非结构化同位网格的SIMPLE算法

通过基于非结构化网格的有限体积法对二维稳态Navier—Stokes方程进行了数值求解。其中对流项采用延迟修正的二阶格式进行离散；扩散项的离散采用二阶中心差分格式；对于压力-速度耦合利用SIMPLE算法进行处理；计算节点的布置采用同位网格技术，界面流速通过动量插值确定。本文对方腔驱动流、倾斜腔驱动流和圆柱外部绕流问题进行了计算，讨论了非结构化同位网格有限体积法在实现SIMPLE算法时，迭代次数与欠松弛系数的关系、不同网格情况的收敛性、同结构化网格的对比以及流场尾迹结构。通过和以往结果比较可知，本文的方法是准确和可信的。     

Finite element solution of multi-dimensional two-phase flow through porous media with arbitrary heating conditions

Mechanistic models for flow regime transitions and drag forces proposed in an earlier work are employed to predict two-phase flow characteristics in multi-dimensional porous layers. The numerical scheme calls for elimination of velocities in favor of pressure and void fraction. The momentum equations for vapor and liquid then can be reduced to a system of two partial differential equations (PDEs) which must be solved simultaneously for pressure and void fraction.

Solutions are obtained both in two-dimensional cartesian and in axi-symmetric coordinate systems. The porous layers in both cases are composed of regions with different permeabilities. The finite element method is employed by casting the PDEs in their equivalent variational forms. Two classes of boundary conditions (specified pressure and specified fluid fluxes) can be incorporated in the solution. Volumetric heating can be included as a source term. The numerical procedure is thus suitable for a wide variety of geometry and heating conditions. Numerical solutions are also compared with available experimental data.     

Evolution of the balance equations in saturated thermoelastic porous media following abrupt simultaneous changes in pressure and temperature

A mathematical model is developed for saturated flow of a Newtonian fluid in a thermoelastic, homogeneous, isotropic porous medium domain under nonisothermal conditions. The model contains mass, momentum and energy balance equations. Both the momentum and energy balance equations have been developed to include a Forchheimer term which represents the interaction at the solid-fluid interface at high Reynolds numbers. The evolution of these equations, following an abrupt change in both fluid pressure and temperature, is presented. Using a dimensional analysis, four evolution periods are distinguished. At the very first instant, pressure, effective stress, and matrix temperature are found to be disturbed with no attenuation. During this stage, the temporal rate of pressure change is linearly proportional to that of the fluid temperature. In the second time period, nonlinear waves are formed in terms of solid deformation, fluid density, and velocities of phases. The equation describing heat transfer becomes parabolic. During the third evolution stage, the inertial and the dissipative terms are of equal order of magnitude. However, during the fourth time period, the fluid's inertial terms subside, reducing the fluid's momentum balance equation to the form of Darcy's law. During this period, we note that the body and surface forces on the solid phase are balanced, while mechanical work and heat conduction of the phases are reduced.     

Numerical study of concentrated fluid–particle suspension flow in a wavy channel

The particle migration effects and fluid–particle interactions occurring in the flow of highly concentrated fluid–particle suspension in a spatially modulated channel have been investigated numerically using a finite volume method. The mathematical model is based on the momentum and continuity equations for the suspension flow and a constitutive equation accounting for the effects of shear‐induced particle migration in concentrated suspensions. The model couples a Newtonian stress/shear rate relationship with a shear‐induced migration model of the suspended particles in which the local effective viscosity is dependent on the local volume fraction of solids. The numerical procedure employs finite volume method and the formulation is based on diffuse‐flux model. Semi‐implicit method for pressure linked equations has been used to solve the resulting governing equations along with appropriate boundary conditions. The numerical results are validated with the analytical expressions for concentrated suspension flow in a plane channel. The results demonstrate strong particle migration towards the centre of the channel and an increasing blunting of velocity profiles with increase in initial particle concentration. In the case of a stenosed channel, the particle concentration is lowest at the site of maximum constriction, whereas a strong accumulation of particles is observed in the recirculation zone downstream of the stenosis. The numerical procedure applied to investigate the effects of concentrated suspension flow in a wavy passage shows that the solid particles migrate from regions of high shear rate to low shear rate with low velocities and this phenomenon is strongly influenced by Reynolds numbers and initial particle concentration. Copyright © 2008 John Wiley & Sons, Ltd.