An approach of dynamic mesh adaptation for simulating 3‐dimensional unsteady moving‐immersed‐boundary flows

In this paper, we present an approach of dynamic mesh adaptation for simulating complex 3‐dimensional incompressible moving‐boundary flows by immersed boundary methods. Tetrahedral meshes are adapted by a hierarchical refining/coarsening algorithm. Regular refinement is accomplished by dividing 1 tetrahedron into 8 subcells, and irregular refinement is only for eliminating the hanging points. Merging the 8 subcells obtained by regular refinement, the mesh is coarsened. With hierarchical refining/coarsening, mesh adaptivity can be achieved by adjusting the mesh only 1 time for each adaptation period. The level difference between 2 neighboring cells never exceeds 1, and the geometrical quality of mesh does not degrade as the level of adaptive mesh increases. A predictor‐corrector scheme is introduced to eliminate the phase lag between adapted mesh and unsteady solution. The error caused by each solution transferring from the old mesh to the new adapted one is small because most of the nodes on the 2 meshes are coincident. An immersed boundary method named local domain‐free discretization is employed to solve the flow equations. Several numerical experiments have been conducted for 3‐dimensional incompressible moving‐boundary flows. By using the present approach, the number of mesh nodes is reduced greatly while the accuracy of solution can be preserved.     

机翼大攻角下失速颤振的气动弹性研究

本文采用流固耦合的数值方法研究了机翼在0°～50°攻角下的颤振。计算结果表明,随着来流攻角α0的增大,机翼的固有频率对颤振的影响越来越大,颤振由线性的强迫振动逐渐发展成为非线性的自激振动,而当α0增加到一定程度以后,大尺度分离流交替地从机翼头部和尾部产生,机翼和流场会发生共振,引起机翼的失速颤振。     

1+1/2对转涡轮叶排轴向间距对性能影响的研究

木文对1+1/2对转涡轮中轴向间距影响进行了初步研究。对不同轴向间距叶栅流场进行了大量时间精确模拟。研究表明,1+1/2对转涡轮高低压转叶间轴向间距成为追求高效率的制约因素,轴向间距为高压转叶喉部宽度是关键点。因此有关小轴向间距下强激波/叶排干扰的研究应该加强。     

叶型优化在非定常流动条件下效果的数值分析

本文对叶型优化在非定常流动条件下的效果进行数值研究。叶型优化基于比较成熟的在定常流动条件下的反问题解法，非定常流动则考虑到动静叶栅之间的相互作用，非定常流动的求解在Ｎ－Ｓ方程的基础上采用分区计算的方法来完成。数值计算的结果表明，定常流动条件下优化得到的叶型在非定常流动条件下同样具有较好的气动性能。     

A time‐accurate pseudo‐compressibility approach based on an unstructured hybrid finite volume technique applied to unsteady turbulent premixed flame propagation

A preconditioning approach based on the artificial compressibility formulation is extended to solve the governing equations for unsteady turbulent reactive flows with heat release, at low Mach numbers, on an unstructured hybrid grid context. Premixed reactants are considered and a flamelet approach for combustion modelling is adopted using a continuous quenched mean reaction rate. An overlapped cell‐vertex finite volume method is adopted as a discretisation scheme. Artificial dissipation terms for hybrid grids are explicitly added to ensure a stable, discretised set of equations. A second‐order, explicit, hybrid Runge–Kutta scheme is applied for the time marching in pseudo‐time. A time derivative of the dependent variable is added to recover the time accuracy of the preconditioned set of equations. This derivative is discretised by an implicit, second‐order scheme. The resulting scheme is applied to the calculation of an infinite planar (one‐dimensional) turbulent premixed flame propagating freely in reactants whose turbulence is supposed to be frozen, homogeneous and isotropic. The accuracy of the results obtained with the proposed method proves to be excellent when compared to the data available in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

进口畸变条件下改善压气机性能的数值研究

依据"非定常两代流型"理论,验证了二维叶栅在均匀进气条件下,施加激励对其性能所起的改善作用,尔后重点研究了进口具有总压的周向畸变时,施加激励对其气动性能的影响。计算结果表明,在有无进口畸变时,施加一定的激励都会改善流场结构,气动参数有积极变化,压升和效率增加,特别是在流量比较小时,变化尤为明显.对于进口畸变和施加激励影响压气机性能的内在原因,从流场结构的变化方面进行了初步分析。     

A numerical study of an unsteady laminar flow in a doubly constricted 3D vessel

Unsteady flow dynamics in doubly constricted 3D vessels have been investigated under pulsatile flow conditions for a full cycle of period T. The coupled non‐linear partial differential equations governing the mass and momentum of a viscous incompressible fluid has been numerically analyzed by a time accurate Finite Volume Scheme in an implicit Euler time marching setting. Roe's flux difference splitting of non‐linear terms and the pseudo‐compressibility technique employed in the current numerical scheme makes it robust both in space and time. Computational experiments are carried out to assess the influence of Reynolds' number and the spacing between two mild constrictions on the pressure drop across the constrictions. The study reveals that the pressure drop across a series of mild constrictions can get physiologically critical and is also found to be sensitive both to the spacing between the constrictions and the oscillatory nature of the inflow profile. The flow separation zone on the downstream constriction is seen to detach from the diverging wall of the constriction leading to vortex shedding with 3D features earlier than that on the wall in the spacing between the two constrictions. Copyright © 2002 John Wiley & Sons, Ltd.     

沉浮,俯仰及控制面振动引起的跨音速非定常气动载荷计算

本文将文献[9]提出改进的通量分裂方法,应用于随时间变化的贴体网格中,建立了可用于求解非定常Euler方程的通量分裂方法.该方法是以连续的特征值分离为基础,它具有方法简单,便于推广使用的特点.同时克服了Steger-Warming通量分裂方法存在的问题.对通量分裂后的Euler方程.利用MUSCL型迎风差分建立了具有二阶精度的有限体积方程.文中以NACA64A—10翼型为例,对其在跨音速流场中进行沉浮、俯仰及带有振动控制面引起的非定常气动载荷进行了计算.部分计算结果与相应的实验结果进行了比较,吻合良好     

A projection method for solving incompressible viscous flows on domains with moving boundaries

In this paper, a projection method is presented for solving the flow problems in domains with moving boundaries. In order to track the movement of the domain boundaries, arbitrary‐Lagrangian–Eulerian (ALE) co‐ordinates are used. The unsteady incompressible Navier–Stokes equations on the ALE co‐ordinates are solved by using a projection method developed in this paper. This projection method is based on the Bell's Godunov‐projection method. However, substantial changes are made so that this algorithm is capable of solving the ALE form of incompressible Navier–Stokes equations. Multi‐block structured grids are used to discretize the flow domains. The grid velocity is not explicitly computed; instead the volume change is used to account for the effect of grid movement. A new method is also proposed to compute the freestream capturing metrics so that the geometric conservation law (GCL) can be satisfied exactly in this algorithm. This projection method is also parallelized so that the state of the art high performance computers can be used to match the computation cost associated with the moving grid calculations. Several test cases are solved to verify the performance of this moving‐grid projection method. Copyright © 2004 John Wiley Sons, Ltd.     

Interfacial dynamics‐based modelling of turbulent cavitating flows,Part‐2: Time‐dependent computations

The interfacial dynamics‐based cavitation model, developed in Part‐1, is further employed for unsteady flow computations. The pressure‐based operator‐splitting algorithm (PISO) is extended to handle the time‐dependent cavitating flows with particular focus on the coupling of the cavitation and turbulence models, and the large density ratio associated with cavitation. Furthermore, the compressibility effect is important for unsteady cavitating flows because in a water–vapour mixture, depending on the composition, the speed of sound inside the cavity can vary by an order of magnitude. The implications of the issue of the speed of the sound are assessed with alternative modelling approaches. Depending on the geometric confinement of the nozzle, compressibility model and cavitation numbers, either auto‐oscillation or quasi‐steady behaviour is observed. The adverse pressure gradient in the closure region is stronger at the maximum cavity size. One can also observe that the mass transfer process contributes to the cavitation dynamics. Compared to the steady flow computations, the velocity and vapour volume fraction distributions within the cavity are noticeably improved with time‐dependent computations. Copyright © 2004 John Wiley & Sons, Ltd.