??????????????????У????

为避免用四阶紧致格式求解泊松方程所具有编程复杂和难以实现的困难,对 传统的五点二阶中心差分格式进行改进;通过增加对残差的校正计算,提出了一种新型具有 四阶精度的两步预估校正格式. 新格式虽需要增加一定的计算量,但它的格式精度高,边界 条件处理极简单,易于编程实现. 通过数值实验,结果证明上述格式的确具有易于编程和计 算精度高的优点. 预估校正格式很容易推广到其他复杂情形.     

二维泊松方程的两步预估校正格式

为避免用四阶紧致格式求解泊松方程所具有编程复杂和难以实现的困难,对传统的五点二阶中心差分格式进行改进;通过增加对残差的校正计算,提出了一种新型具有四阶精度的两步预估校正格式.新格式虽需要增加一定的计算量,但它的格式精度高,边界条件处理极简单,易于编程实现.通过数值实验,结果证明上述格式的确具有易于编程和计算精度高的优点.预估校正格式很容易推广到其他复杂情形.     

求解非线性结构动力方程的预估校正-辛时间子域法

提出了求解非线性结构动力方程的预估校正-辛时间子域法。首先,将结构非线性动力方程转换为状态空间方程,在任一时间子域内利用改进的欧拉法对各离散时刻的状态变量值进行预估和校正。然后,将离散的非线性项用Lagrange插值多项式展开并视为外荷载,结合辛时间子域法即可求解非线性动力系统的响应。这种方法不必对状态矩阵求逆,无需计算高阶导数,计算简单,格式统一,易于编程。算例结果表明,本文方法具有较高的计算精度、效率和稳定性,是一种求解非线性结构动力方程的有效方法。     

对流扩散方程的迎风变换及相应有限差分方法

本文提出所谓迎风变换,将对流扩散方程分解为对流迎风函数和扩散方程,并构造相应的有限差分格式。对流迎风函数以简明的指数解析形式反映对流扩散现象的迎风效应,原则上消除了源于不对称对流算子的困难,能够便利对流扩散方程的数值求解。有限差分格式具有二阶精度和无条件稳定性,算例表明其准确性、收敛速度及对边界层效应的适应能力均明显优于中心差分格式和迎风差分格式。     

三维对流扩散方程的高精度全隐式多重网格方法

提出了数值求解三维非定常变系数对流扩散方程的一种高精度全隐紧致差分格式,该格式在空间上具有四阶精度,时间具有二阶精度。为了克服传统迭代法在每一个时间步上迭代求解隐格式时收敛速度慢的缺点,采用多重网格加速技术,建立了适用于本文高精度全隐紧致格式的多重网格算法,从而大大加快了迭代收敛速度。数值实验结果验证了本文方法的精确性、稳定性和对高网格雷诺数问题的强适应性。     

燃烧室两相流场亚网格燃烧模型的研究

在三维任意曲线坐标系下采用不同的亚网格燃烧模型对环形燃烧室火焰筒气液两相湍流瞬态反应流进行大涡模拟．计算中所采用的数学度模型有：k方程亚网格尺度模型估算亚网格湍流黏性;热通量辐射模型估算辐射换热,分别采用亚网格EBU燃烧模型（E-A model）、亚网格二阶矩输运方程模型（SOM）和亚网格二阶矩代数模型（SOM-A）估算化学反应速率．并在非交错网格系统下气相采用SIMPLE算法和混合差分格式求解,液相采用Lagrange处理,并用PSIC算法对其进行求解．通过实验结果和计算结果的比较,表明在三维任意曲线坐标系下对燃烧室火焰简两相湍流油雾燃烧流场进行大涡模拟,3种不同的亚网格燃烧模型都能真实反映两相湍流化学反应流流动及实际燃烧过程,而采用亚网格二阶矩输运方程模型稍优于其他两种亚网格燃烧模型．     

高超声速三维化学非平衡流动简化NS方程数值解

本文研究用时间稳定法求解高超声速三维化学非平衡流动简化ＮＳ方程，把整个非平衡流体动力学问题分裂为流体力学问题和化学问题。流体力学与化学问题之间的耦合通过迭代方式加以处理。考虑高温空气分解与电离，七个化学组元，并认为化学反应是以有限速率进行的．差分格式采用无自由参数的混合反扩散格式，激波及壁面处理应用特征线理论．对来流攻角为１０°，马赫数为２６，飞行高度为８５ｋｍ的钝锥体，给出了流场参数和物面加热率的数值计算结果．     

二维对流扩散方程的高精度全隐式多重网格方法

提出了数值求解二维非定常变系数对流扩散方程的一种时间二阶、空间四阶精度的三层全隐紧致差分格式。为了加快迭代求解隐格式时在每一个时间步上的收敛速度，采用多重网格加速技术，建立了适用于本文高精度金隐紧致格式的多重网格算法。数值实验结果验证了本文方法的精确性、稳定性和对高网格雷诺数问题的强适应性。     

低浓度固液两相流的颗粒相动理学模型

用广义Fokker-Planck扩散模型描述液相湍动对颗粒的挟带作用，用修正的BGK模型描述粒间碰撞效应，建立了封闭的颗粒相PDF输运方程．运用Chapman-Enskog迭代法求得方程的二阶近似解，获得颗粒相脉动速度二阶矩和三阶矩闭合关系．模型与颗粒流模型相容，与液相湍流闭合模型是否相容依赖于扩散模型的具体形式，并据此比较了不同的涡一颗粒作用模型．模型与二维明渠流轻质沙和天然沙试验资料符合很好．表明细小粒径颗粒能够充分跟随水流运动；大粒径颗粒的相间平均速度差和壁面滑移速度明显，近壁区内的颗粒沿流向和垂向脉动强度都可能大于水流，并存在一定程度的颗粒碰撞效应．     

靠近排列的两个交错方柱三维绕流的数值模拟

通过用时间分裂算法求解Navier－Stokes方程,对中等Reynolds数下靠近排列的两个交错方柱三维绕流进行了数值模拟,其中,中间速度场用四阶Adams格式计算,压力场通过结合近似因子分解方法AF1与稳定的双共轭梯度方法Bi-CGSTAB进行迭代求解．数值模拟发现当两个方柱靠得较近时,有互相吸引趋势,而且上游方柱的Strouhal数较大．方柱的交错排列方式对绕流影响明显．计算结果与实验定性吻合,而且比用MAC－AF1方法计算的结果好．     

NUMERICAL STUDIES ON THE MIXING OF CH_4 AND KEROSENE INJECTED INTO A SUPERSONIC FLOW WITH H_2 PILOT INJECTION

Ｎｏｍｅｎｃｌａｔｕｒｅｃｉ 　Ｍａｓｓｆｒａｃｔｉｏｎ ;ｃｉｘ,ｃｉｙ,ｃｉｚ 　Ｄｅｒｉｖａｔｉｖｅｓｏｆｇｒａｄｉｅｎｔｏｆｃｉｃｐ 　ＣｏｎｓｔａｎｔｐｒｅｓｓｕｒｅｓｐｅｃｉｆｉｃｈｅａｔｃＤ 　Ｄｒａｇｃｏｅｆｆｉｃｉｅｎｔｄ　ＤｒｏｐｌｅｔｄｉａｍｅｔｅｒＤ　ＴｏｔａｌｄｉｆｆｕｓｉｏｎｃｏｅｆｆｉｃｉｅｎｔＥ　ＴｏｔａｌｅｎｅｒｇｙｐｅｒｕｎｉｔｖｏｌｕｍｅＦ　Ｄｒａｇｆｏｒｃｅｏｆｄｒｏｐｌｅｔｓｆｘ,ｆｙ,ｆｚ 　ＣｏｍｐｏｎｅｎｔｓｏｆＦ^Ｈｇ,^Ｈｌ 　Ｓｏｕｒｃｅｔｅｒｍｓｏｆｇａ…     

煤油液滴直径对两相旋转爆轰发动机流场的影响

为探究煤油液滴不同初始直径对气液两相旋转爆轰发动机流场的影响,假设初始注入的煤油液滴具有均匀直径,考虑雾化破碎、蒸发等过程,建立了非定常两相爆轰的Eulerian-Lagrangian模型,进行了液态煤油/高温空气爆轰的非预混二维数值模拟。结果表明：在初始液滴直径为1～70μm的工况范围,燃烧室内均形成了单个稳定传播的旋转爆轰波;全局当量比为1时,爆轰波前的空气区域大于液滴煤油的蒸气区域,导致波前燃料空气混合不均匀,波前均存在富油区和贫油区,两相速度差导致分离出的空气形成低温条带;当煤油液滴的初始直径较小时,波前的反应物混合过程主要受蒸发的影响,爆轰波可稳定传播;当直径减小至1μm时,煤油液滴在入口处即蒸发,旋转爆轰波表现为气相传播的特性,爆轰波结构平整;当煤油液滴的初始直径较大时,波前的反应物混合过程主要受液滴破碎的影响;对于相同的燃料质量流量,在不同初始煤油液滴直径工况下,煤油液滴最大的停留时间均占爆轰波传播时间尺度的80%以上;爆轰波前燃料预蒸发为气相的占比越高,爆轰波的传播速度越高;初始液滴直径为10～70μm的工况范围内,爆轰波的速度随初始直径的增大先升高后降低。     

Numerical Studies on the Mixing of CH4 and Kerosene Injected into a Supersonic Flow with H2 Pilot Injection

Two-fluid model and divisional computation techniques were used. The multi-species gas fully N-S equations were solved by upwind TVD scheme. Liquid phase equations were solved by NND scheme. The phases-interaction ODE equations were solved by 2nd Runge-Kutta approach. The favorable agreement is obtained between computational results and PLIF experimental results of iodized air injected into a supersonic flow. Then, the numerical studies were carried out on the mixing of CH ₄ and kerosene injected into a supersonic flow with H ₂ pilot injection. The results indicate that the penetration of kerosene approaches maximum when it is injected from the second injector. But the kerosene is less diffused compared with the gas fuels. The free droplet region appears in the flow field. The mixing mechanism of CH ₄ with H ₂ pilot injection is different from that of kerosene. In the staged duct, H ₂ can be entrained into both recirculation zones produced by the step and injectors. But CH ₄ can only be carried into the recirculation between the injectors. Therefore, initiations of H ₂ and CH ₄ can occur in those regions. The staged duct is better in enhancing mixing and initiation with H ₂ pilot flame.     

Simulation of the primary breakup of a high-viscosity liquid jet by a coaxial annular gas flow

The conversion of low-grade fossil and biogenic energy resources (petcoke, biomass) to a synthesis gas in a high pressure entrained flow gasification process opens a wide spectrum for high efficient energy conversion processes. The synthesis gas can be used for production of methane (SNG), liquid fuels (BtL, CtL) or as fuel for operation of a gas turbine in a combined cycle power plant (IGCC). The production of a tar free high quality syngas is a challenging objective especially due to the fact that typical liquid or suspension fuels for entrained flow gasifiers feature viscosities up to 1000 mPas. Fuel droplet conversion at typical entrained flow gasification conditions is characterized by heat up, evaporation and subsequent degradation of the vapour phase. To guarantee a high fuel conversion rate in the gasifier an efficient atomization of the fuel is required. Mainly twin-fluid burner nozzles are used for atomization of those typically high viscous fuels. The present study is focused on the assessment of the accuracy of CFD computations for the primary breakup of high-viscosity liquids using an external mixing twin fluid nozzle. In a first step experiments were performed with a Newtonian glycerol-water-mixture featuring a liquid viscosity of 400 mPas. Jet breakup was investigated using a high speed camera as well as PIV and LDA-System for a detailed investigation of the flow field. In a second step the experimental results serve as reference data to assess the accuracy of CFD computations. Compressible large eddy simulations (LES) were performed to capture the morphology of the primary breakup as well as the important flow field characteristics. A Volume of Fluid (VOF) approach was used to track the unsteady evolution and breakup of the liquid jet. Comparison of experimental and numerical results showed good agreement with respect to breakup frequency, velocity fields and morphology.     

Numerical investigation of turbulent-jet primary breakup using one-dimensional turbulence

Primary breakup to form droplets at liquid surfaces is an important fundamental process to study as it determines the initial properties of the dispersed phase, which affect mixing rates, secondary breakup, droplet collisions, and flow separation within the dispersed flow region. Primary breakup can be regarded as one of the least developed model components for simulating and predicting liquid jet breakup. However, it is of paramount importance in many technical applications, e.g. fuel injection in engines and spray painting. This paper presents a numerical investigation of primary breakup of a turbulent liquid jet in still air at standard conditions using the one-dimensional turbulence (ODT) modeling framework. ODT is a stochastic model that simulates turbulent flow evolution along a notional 1D line of sight by applying instantaneous maps to represent the effect of individual turbulent eddies on property profiles. An important feature of ODT is the resolution of all relevant scales, both temporal and spatial. The restriction to one spatial dimension in ODT permits affordable high resolution of interfacial and single-phase property gradients, which is key to capturing the local behavior of the breakup process and allows simulations at high Reynolds and Weber numbers that are currently not accessible to direct numerical simulations (DNS).This paper summarizes our extensions of the ODT model to simulate geometrically simple jet breakup problems, including representations of Rayleigh wave breakup, turbulent breakup, and shear-driven breakup. Each jet breakup simulation consists of a short temporal channel section to initialize a turbulent velocity profile at the nozzle exit followed by an adjacent jet section. The simulations are carried out for jet exit Reynolds number of 11,500, 23,000, 46,000 and 92,000 while the Weber number is varied within the range 10²–10⁷. We present results on breakup statistics including spatial locations of droplet release, droplet sizes and liquid core length. The results on primary breakup are compared to experimental results and models.     

Evaporation of bi-component droplets in a heated, highly turbulent flow

This work aims to understand the phenomena that occur in a combustion chamber where multi-component fuel droplets are injected. Many evaporation models exist but the influence of turbulence on spray vaporization is not yet well understood. This study gives a useful database to improve these models. The objective of the work is to measure the dispersion and the evaporation of bi-component (octane/3-pentanone) droplets and the resulting vapor mixing in a well-known, heated, highly turbulent channel flow. The carrier flow shows high turbulence levels, flat profiles for the mean velocity and the velocity fluctuations. The injected droplets have a large variety of behaviors due to the large polydispersion and to the turbulence. The evolution of 3-pentanone liquid concentration, mass flux, and droplet clusters are described. Mean concentration, fluctuations of concentration, and mixing of the vapor phase are characterized.     

Shock wave interactions with particles and liquid fuel droplets

Shock waves traveling through a multiphase flow environment are studied numerically using the Flux Corrected Transport (FCT) algorithm. Both solid particles and liquid droplets are used as the dispersed phase with their trajectories being computed using a Lagrangian tracking scheme. The phases are coupled by including source terms which account for mass transfer, momentum, and energy exchange from the dispersed phase in the governing equations of motion for the gas phase. For solid particles, droplet size effects are examined at constant mass loading. Deceleration of the shock wave is observed with effects increasing with decreasing particle size. The equilibrium velocity attained is found to agree with analytical results for an equivalent dense gas with a modified specific heat ratio. For liquid droplets, a droplet breakup model is introduced and the results show a faster attenuation rate than with the solid particle model. The inclusion of vaporization to the breakup model is seen to increase the attenuation rate but does not alter the final equilibrium velocity. When an energy release model is used in the simulations, behavior resembling a detonation is observed under certain conditions, with energy release coupling with and accelerating the shock front. Received 17 July 2000 / Accepted 20 August 2002 / Published online 4 December 2002 Correspondence to: Dr. K. Kailasanath (e-mail: kailas@lcp.nrl.navy.mil)     

Flow characteristics of spray impingement in PFI injection systems

The present paper addresses an experimental study of the dynamic exchanges between the impact of an intermittent spray and the liquid film formed over the target, based on detailed phase Doppler anemometry (PDA) measurements of droplet size, velocity and volume flux in the vicinity of the impact. The flow configuration is that of a pulsed injector spraying gasoline onto a flat disc to simulate the port fuel injection (PFI) of an internal combustion engine operating at cold-start conditions. The measurements evidence that the outcome of impact cannot be accurately predicted based on the characteristics of the free spray, but requires precise knowledge of the flow structure, induced by the target. The implications for spray–wall interaction modelling are then discussed based on the application of conservation equations to the mass, momentum and energy exchanged between the impinging droplets and the liquid film. The results show that the liquid film starts to form in the vicinity of the stagnation region at early stages of injection and a non-negligible proportion of droplets impinging at outer regions splash after interaction with the film. Film disruption is mainly driven by the intermittent axial momentum of impinging droplets, which enhances the vertical oscillations. The radial momentum imparted to the liquid film at the stagnation region is fed back onto secondary droplets emerging later during the injection cycle at outwards locations, where momentum of impacting droplets is much smaller. As a consequence, although the number of splashed droplets is enhanced by normal momentum, their size and ejection velocity depends more on the radial spread induced onto the liquid film and, hence, on the radial momentum at impact. The analysis further shows that existing spray–wall interaction models can be improved if the dynamic exchanges between the impacting spray and the liquid film are accounted.     

Shock–Free Breakup of Droplets. Temporal Characteristics

In the present paper, we consider the shock–free breakup of droplets in their encounter with a layer (sheet) of a moving gas in the absence of pressure perturbations when the droplets are affected by a short U–shaped pulse of aerodynamic forces. Under a high pressure of the ambient gas medium p₀ = 20—80 bar, the droplets (ethanol or liquid oxygen) have a chance to break up after stay in a thing (2—5 mm thick) gas layer (jet) moving with a velocity of 1—10 m/sec. A distinctive feature of the process is that the characteristic time of droplet deformation and the period of natural oscillations coincide with the residence time for the droplets in the region of their interaction with the gas stream. Empirical formulas are proposed for determination of the total breakup time and the duration of the droplet disintegration stage in shock–free breakup.     

A STUDY ON DROPLET VELOCITY OF ETHANOL DURING ELECTROSPRAYING PROCESS AT CONE-JET MODE

研究液滴在静电喷雾下的速度特性是理解喷雾形态的形成及演化的关键.结合锥射流模式下乙醇静电喷雾实验数据,建立了静电喷雾二维轴对称模型.基于离散相液滴运动方程、连续相空气运动方程、电场方程以及用户自定义函数,进行了数值求解,获得了锥射流模式下的乙醇静电喷雾形态、空间电场分布以及液滴速度场分布.考虑了不同空气入口流速的影响,得到了乙醇/空气同轴射流静电喷雾形态的变化规律.结果表明,喷雾外围液滴与空气流场有较强的相互作用,导致喷雾中轴线附近的液滴速度分布变化较小,而在喷雾外围处的液滴速度分布沿径向剧烈变化;随着空气入口速度的增大,乙醇/空气同轴射流静电喷雾形态先趋于发散,当空气入口速度大于喷雾外围液滴轴向速度时,喷雾形态则趋于聚拢.因此,除改变施加电压、液体流量和电极结构外,通过控制空气入口速度来影响喷雾液滴速度场,也可获得不同的静电喷雾效果.