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.     

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.     

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.     

Calculation of impulsively started incompressible viscous flows

The paper's focus is the calculation of unsteady incompressible 2D flows past airfoils. In the framework of the primitive variable Navier–Stokes equations, the initial and boundary conditions must be assigned so as to be compatible, to assure the correct prediction of the flow evolution. This requirement, typical of all incompressible flows, viscous or inviscid, is often violated when modelling the flow past immersed bodies impulsively started from rest. Its fulfillment can however be restored by means of a procedure enforcing compatibility, consisting in a pre‐processing of the initial velocity field, here described in detail. Numerical solutions for an impulsively started multiple airfoil have been obtained using a finite element incremental projection method. The spatial discretization chosen for the velocity and pressure are of different order to satisfy the inf–sup condition and obtain a smooth pressure field. Results are provided to illustrate the effect of employing or not the compatibility procedure, and are found in good agreement with those obtained with a non‐primitive variable solver. In addition, we introduce a post‐processing procedure to evaluate an alternative pressure field which is found to be more accurate than the one resulting from the projection method. This is achieved by considering an appropriate ‘unsplit’ version of the momentum equation, where the velocity solution of the projection method is substituted. Copyright © 2004 John Wiley & Sons, Ltd.     

Use of intermittent jets to enhance flux in crossflow filtration

This paper deals with the influence of a new flow unsteadiness on the permeate fluxes in crossflow filtration. A pneumatically controlled valve generates intermittent jets from the main flow, causing the formation of large vortices moving downstream along the tubular membrane. The main results of the numerical calculation of such flows are given. The experimental study was carried out by filtering a bentonite suspension through an ultrafiltration mineral membrane. Time evolutions of flux were achieved in steady and unsteady operating conditions. Results concerning the influence and limits of the nozzle to tube diameter ratio and the jet velocities are discussed. The applicability of such an unsteady flow is examined with a view to effects on energy consumption and possible viscosity effects.     

关于河道一维非恒定流水温预测模型的研究

作者在全面地分析影响河道水温变化各种因素的基础上，根据大渡河流域实测资料所推荐的有关计算纵向离散系数的经验公式．确定了一种适宜于河道一维非恒定流水温的预测模型．作者运用此模型预测了水库下泄低温水、火电厂温排水对下游水环境的影响及其应用价值，其结果可为制定水库生态用水调度方案和流域开发的环境影响评价提供科学依据．     

动态油膜-转子系统的支座松动故障研究

采用张文等建立的短轴承非稳态非线性油膜力模型,利用数值模拟等方法,分析了带有一端支座松动故障的转子 轴承系统的复杂运动现象,讨论了转速比变化时此类系统中具有周期、拟周期、混沌等多种形式的运动·研究结果表明:当转子 轴承系统出现松动故障时,即便是在较高转速区也会出现周期3分频运动,而且,松动端轴颈的运动轨迹以及支座的相轨迹都呈现出特有的形状,这为有效识别转子 轴承系统的支座松动故障提供了理论参考·