采用分段湍流模式研究绕水翼的空化流动

结合空化流动特点,建立了一个包含空间尺度信息的分段湍流模式。计算中,应用基于质量传输空化模型,分别采用三种湍流模型计算了绕Clark-y型水翼云状空化流动,得到了随时间变化的空泡形态以及升、阻力等流场和动力特性。通过与实验结果的对比,表明这三种湍流模型均能捕捉云状空化区域的空泡形态和空泡脱落的非定常细节。分段湍流模式能够更好地调整流场内的湍流黏性,更精确地预测空穴长度和空穴尾部水汽分布,与实验结果吻合较好。     

二维空化、超空化流动的数值模拟

通过引入混合密度函数发展了一种新的空化模型,运用Navier-Stokes方程加k-ε两方程湍流模型对Ep-pler E817和NACA 0015水翼的空化流动进行了数值模拟,得到了从空化到超空化的一系列变化过程的数值模拟结果。通过与实验数据进行对比,发现二者吻合较好。     

空化、超空化流动的数值模拟方法研究

基于结构化网格,运用可压缩流N—S方程及k-ε湍流模型对流场进行求解,在低压区域引入一种基于混合密度函数的空化模型对轴对称体的空化、超空化流动进行了数值模拟．通过将半球圆柱的计算结果与实验数据和前人的计算结果进行对比,验证了所发展的数值方法的可靠性．最后,采用非定常的数值方法,研究了钝头体射弹的空化、超空化流动特性,并模拟了其超空泡的发展过程．     

考虑当地涡旋运动修正的湍流模型在非定常空化湍流流场计算中的应用

为了更精确地模拟当地涡旋运动对复杂非定常空化湍流流动的影响,考虑当地涡旋运动特性对标准k-ε湍流模型进行了修正;并采用该修正模型对二维Clark-Y水翼周围非定常空化湍流场进行了计算;通过与实验结果进行对比,分析并验证了该修正模型计算结果的有效性.研究表明:该修正湍流模型根据当地涡旋效应对湍动能输运方程中的生成项进行修正,可以有效解决标准k-ε模型对湍动能的过度预测问题,得到的空化核心区域时均湍动能和湍流黏性的预测结果降低了约30％,与实际值更为接近;标准k-ε模型计算所得的时均空泡长度不大于0.6倍弦长,而修正模型计算所得的时均空穴可以覆盖水翼吸力面;修正模型可以准确预测水翼吸力面表面空泡的非定常脱落过程,包括附着空穴的拉伸和断裂,以及脱落空泡的逐渐消失.     

PANS模型在空化湍流数值计算中的应用

采用一种基于标准k-ε模型改进的局部时均化模型（Partially-Averaged Navier-Stokes Model,PANS）,并应用于空化流动计算。控制不同的模型参数,分别对绕平头轴对称回转体和Clark-Y型水翼的空化流动进行模拟,并与实验结果进行对比。结果表明：PANS模型中未分解湍动能比率fk的取值对预测空化流动的数值计算精度有重要影响,改变fk的取值可实现对不同滤波尺度范围内的求解;随着fk值的减小PANS的预测精度逐步提高,能在相对较大范围内求解较小尺度的湍流运动过程中,预测到湍流运动中强烈的非定常特性;同时可以比较准确地预测空化流场结构和动力特性。     

一种修正的低温流体空化流动计算模型

为了更准确地预测低温流体的空化流动特性, 基于Kubota空化模型, 对蒸发和凝结源项进行修正, 建立了一种考虑热力学效应的空化模型. 分别采用原始和修正的Kubota空化模型, 计算了绕对称回转体液氮的空化流动, 通过与实验结果的比较对修正的空化模型进行了评价. 结果表明, 与原Kubota空化模型比较, 修正的空化模型由于考虑了热力学效应, 计算获得的蒸发量减小, 凝结量增大, 空穴长度减小, 空穴界面形态呈模糊状态.计算结果与实验结果更加一致, 说明修正的空化模型能准确的描述低温流体空化过程的质量传输过程, 能够更准确模拟低温流体中的空化流动特性.      

Turbulence and cavitation models for time-dependent turbulent cavitating flows

Cavitation typically occurs when the fluid pressure is lower than the vapor pressure at a local thermodynamic state,and the flow is frequently unsteady and turbulent.To assess the state-of-the-art of computational capabilities for unsteady cavitating flows,different cavitation and turbulence model combinations are conducted.The selected cavitation models include several widely-used models including one based on phenomenological argument and the other utilizing interface dynamics.The kε turbulence model with additional implementation of the filter function and density correction function are considered to reduce the eddy viscosity according to the computed turbulence length scale and local fluid density respectively.We have also blended these alternative cavitation and turbulence treatments,to illustrate that the eddy viscosity near the closure region can significantly influence the capture of detached cavity.From the experimental validations regarding the force analysis,frequency,and the cavity visualization,no single model combination performs best in all aspects.Furthermore,the implications of parameters contained in different cavitation models are investigated.The phase change process is more pronounced around the detached cavity,which is better illus-trated by the interfacial dynamics model.Our study provides insight to aid further modeling development.     

CFD模拟方法的发展成就与展望

本文综述了计算流体力学(computational fluid dynamics, CFD),尤其是计算空气动力学的发展概况.从计算方法、网格技术、湍流模型、大涡模拟等方面分别总结了CFD所取得的成就,分析了当前存在的问题、困惑,展望了其发展趋势.在CFD计算方法中,主要介绍了中心格式、迎风格式、TVD格式、WENO格式、紧致格式以及间断Galerkin有限元方法,对不同方法的原理和特性进行了系统阐述.网格技术包括结构网格、非结构网格、混合网格以及重叠网格,重点讨论了重叠网格的若干关键技术.在湍流模型中,对目前的模型进行分类介绍,包括线性涡黏性模型、二阶矩模型、非线性模型等,还介绍了转捩模型、DES方法以及SAS方法等.在大涡模拟方法中,就其中若干相关的研究方向进行了探讨,包括滤波方法、亚格子模型、收敛标准、数值格式等. 文中还包含了作者在相关领域的若干研究成果.     

低温介质质量传输空化模型的代理分析及优化

利用基于代理模型的整体敏感度分析方法,对考虑模型参数和物质属性的低温介质质量传输空化模型进行了评价,从而获得影响预测精确度的主要因素,并对空化模型参数进行了校准。结果表明：凝结系数对压力和温度预测精度的影响较大,整体敏感度分别为44%和29%;蒸发系数的影响较小,整体敏感度均小于5%。通过代理模型优化分析获得了Merkle空化模型对压力和温度的预测精度随凝结系数的变化趋势：当凝结系数小于20时,随凝结系数的增大,压力预测精度提高,温度预测精度降低;当凝结系数大于20时,压力和温度的预测精度均随凝结系数的增大而降低。优化后的蒸发系数和凝结系数分别为3.8和19.43,满足Pareto最优解,模型预测能力得到提高。     

Modeling the indentation, penetration and cavitation of elastic membranes

An energy argument is used to predict the penetration of an elastic membrane by a rigid cylindrical indenter. We study a sequence of axisymmetric equilibrium deformations beginning with the indentation of a solid disc, followed by penetration of the membrane at a critical value of the indenter displacement, and ending in a state of cavitation in which the membrane is detached from the indenter and a central traction-free hole is maintained by boundary data alone.     

Modeling droplet dispersion and interphase turbulent kinetic energy transfer using a new dual-timescale Langevin model

Dispersion of spray droplets and the modulation of turbulence in the ambient gas by the dispersing droplets are two coupled phenomena that are closely linked to the evolution of global spray characteristics, such as the spreading rate of the spray and the spray cone angle. Direct numerical simulations (DNS) of turbulent gas flows laden with sub-Kolmogorov size particles, in the absence of gravity, report that dispersion statistics and turbulent kinetic energy (TKE) evolve on different timescales. Furthermore, each timescale behaves differently with Stokes number, a non-dimensional flow parameter (defined in this context as the ratio of the particle response time to the Kolmogorov timescale of turbulence) that characterizes how quickly a particle responds to turbulent fluctuations in the carrier or gas phase. A new dual-timescale Langevin model (DLM) composed of two coupled Langevin equations for the fluctuating velocities, one for each phase, is proposed. This model possesses a unique feature that the implied TKE and velocity autocorrelation in each phase evolve on different timescales. Consequently, this model has the capability of simultaneously predicting the disparate Stokes number trends in the evolution of dispersion statistics, such as velocity autocorrelations, and TKE in each phase. Predictions of dispersion statistics and TKE from the new model show good agreement with published DNS of non-evaporating and evaporating droplet-laden turbulent flow.     

Experimental and numerical investigation of turbulent cross-flow in a staggered tube bundle

This paper presents the results of measurements and numerical predictions of turbulent cross-flow in a staggered tube bundle. The bundle consists of transverse and longitudinal pitch-to-diameter ratios of 3.8 and 2.1, respectively. The experiments were conducted using a particle image velocimetry technique, in a flow of water in a channel at a Reynolds number of 9300 based on the inlet velocity and the tube diameter. A commercial CFD code, ANSYS CFX V10.0, is used to predict the turbulent flow in the bundle. The steady and isothermal Reynolds–Averaged Navier–Stokes (RANS) equations were used to predict the turbulent flow using each of the following four turbulence models: a k-epsilon, a standard k-omega, a k-omega-based shear stress transport, and an epsilon-based second moment closure. The epsilon-based models used a scalable wall function and the omega-based models used a wall treatment that switches automatically between low-Reynolds and standard wall function formulations.

The experimental results revealed extremely high levels of turbulence production by the normal stresses, as well as regions of negative turbulence production. The convective transport by mean flow and turbulent diffusion were observed to be significantly higher than in classical turbulent boundary layers. As a result, turbulence production is generally not in equilibrium with its dissipation rate. In spite of these characteristics, it was observed that the Reynolds normal stresses approximated from the k-based two-equation models were in a closer agreement with experiments than values obtained from the second moment closure. The results show that none of the turbulence models was able to consistently reproduce the mean and turbulent quantities reasonably well. The omega-based models predicted the mean velocities better in the developing region while the epsilon-based models gave better results in the region where the flow is becoming spatially periodic.     

Numerical study of shear band instability and effect of cavitation on the response of a specimen under undrained biaxial loading

This paper presents an extension of the local second gradient model to multiphasic materials (solids particles, air, water) and including the cavitation phenomenon. This new development was made in order to model the response of saturated dilatant materials under deviatoric stress and undrained conditions and possibly, in future, the behavior of unsaturated soils. Some experiments have showed the significance of cavitation for the hydromechanical response of materials. However, to date and as far as we are aware, no attempt was made to implement the cavitation as a phase change mechanism with a control of pore pressure. The first part of the results section explores the effects of permeability, dilation angle and loading rate on the stability of shear bands during a localization event. The reasons underlying the band instability are discussed in detail, which helps defining the conditions required to maintain stability and investigating the effects of cavitation without parasite effect of materials parameters or loading rate. The model showed that, if a uniform response is obtained, cavitation triggers localization. However, in case of a localized solution, cavitation follows the formation of the shear band, with the two events being quite distinct.     

Level set method for numerical simulation of a cavitation bubble,its growth,collapse and rebound near a rigid wall

A level set method of non-uniform grids is used to simulate the whole evolution of a cavitation bubble, including its growth, collapse and rebound near a rigid wall. Single-phase Navier–Stokes equation in the liquid region is solved by MAC projection algorithm combined with second-order ENO scheme for the advection terms. The moving interface is captured by the level set function, and the interface velocity is resolved by “one-side” velocity extension from the liquid region to the bubble region, complementing the second-order weighted least squares method across the interface and projection inside bubble. The use of non-uniform grid overcomes the difficulty caused by the large computational domain and very small bubble size. The computation is very stable without suffering from large flow-field gradients, and the results are in good agreements with other studies. The bubble interface kinematics, dynamics and its effect on the wall are highlighted, which shows that the code can effectively capture the “shock wave”-like pressure and velocity at jet impact, toroidal bubble, and complicated pressure structure with peak, plateau and valley in the later stage of bubble oscillating. The project supported by the National Natural Science Foundation of China (10272032 and 10672043). The English text was polished by Keren Wang.     

Modeling of dynamics, heat transfer, and combustion in two-phase turbulent flows: 1. Isothermal flows

This paper presents a review of authors' collective works in the field of two-phase flow modeling done in the past few decades. The paper is aimed at the construction of mathematical models for simulation of particle-laden turbulent flows. A kinetic equation was obtained for the probability density function (PDF) of the particle velocity distribution in turbulent flows. The proposed kinetic equation describes both the interaction of particles with turbulent eddies of the carrier phase and particle-particle collisions. This PDF equation is used for the derivation of different schemes describing turbulent momentum transfer in the dispersed particle phase. The turbulent characteristics of the gaseous phase are calculated on the basis of the k - turbulence model with a modulation effect of particles on the turbulence.

The constructed models have been applied to the calculation of various two-phase gas-particle turbulent flows in jets and channels as well as particle deposition in tubes and separators. For validating the theoretical and numerical results, a wide range of comparisons with experimental data from Russian and foreign sources has been done.     

Creation, transportation, and coalescence of liquid drops by means of a light beam

A novel principle of manipulation of discrete drops using concentration-capillary forces controlled by the thermal action of a light beam is proposed. The drops are created by the light beam in a thin layer of absorbing solution and in a film of that solution beneath an air bubble in the cell. The possibility of transporting both a single drop and a drop in an air bubble by means of a light beam is demonstrated. For the first time two drops are made to coalesce on a solid substrate by bringing them into contact by means of a light beam.     

Modeling of dynamics, heat transfer, and combustion in two-phase turbulent flows: 2. Flows with heat transfer and combustion

The objective of this part of the paper is to summarize the information concerning the authors' works in the field of simulation of two-phase gas-particle turbulent flows with heat transfer and combustion. A kinetic equation had been derived for the probability density function (PDF) of the particle velocity, temperature, and mass distributions in turbulent flows. This PDF equation is used for the construction of the governing conservation equations of mass, momentum, and heat transfer in the dispersed particle phase.The numerical scheme incorporates two-phase fluid dynamics, convective and radiative heat transfer, and combustion. The proposed models have been applied to the calculation of various particle-laden turbulent flows in jets, combustion and gasification chambers, and furnaces.