展向喷嘴湍流横向射流大涡模拟及其统计分析

采用大涡模拟方法数值模拟了展向椭圆喷嘴的湍流横向射流,对其大尺度结构的时空演化和湍流脉动速度场的时间序列分析、频谱分析、PDF分析以及时、空截面上的统计平均特性进行分析.结果表明,在射流出口附近的下游核心区中速度脉动剧烈,显现出明显的湍流特征.除了三维涡环脱落、扭曲、变形、摆动所对应频率之外,还存在很宽的湍流基频,它与在喷嘴出口附近产生的三维涡环的时空演化过程密切相关.由于展向椭圆喷嘴的湍流横向射流中的三维涡环快速脱落和强相互作用导致射流尾迹中的强湍流脉动,展向椭圆喷嘴湍流横向射流的PDF空间演化特征结构复杂.在射流核心区的湍流偏应力变化平缓,其统计平均值分布接近左右对称.展向椭圆喷嘴的湍流横向射流脉动速度场具有极为复杂的统计行为,与流向椭圆喷嘴相比具有更好的掺混能力.     

高压燃气射流在整装液体中扩展过程的实验和数值模拟

为研究高压燃气射流在整装液体中的扩展规律,探索提高射流稳定性的控制方法,设计了4级渐扩型圆柱观察室.采用高速录像系统,观察了燃气射流在液体模拟工质中扩展与掺混的过程,获得了射流扩展形态的序列照片,测量了Taylor空腔的扩展速度,对比了观察室内壁形状对射流扩展过程的影响.建立了二维轴对称气液两相湍流模型,模拟射流扩展过...     

零重力条件下低温射流抑制大尺寸液氢储罐热分层的数值研究

液氢是一种常用的沸点低、易蒸发的空间低温推进剂. 空间微重力环境中浮力对流被极大减弱,当推进剂储罐壁面存在局部漏热时,储罐内部气液两相流体系会出现环绕漏热源的热分层现象,引起局部过热沸腾,导致储罐内部压力急剧增大,危害系统结构安全. 利用低温射流抑制储罐热分层现象是一种有效手段. 低温流体通过设置在储罐内部的射流喷嘴与储罐内部的流体混合,消减局部高温,实现温度的均匀化. 采用全充满的二维大尺寸储罐模型,对零重力条件下液氢储罐内局部漏热引起的热分层现象开展了数值模拟,主要分析了位于靠近储罐底部的漏热带以及出口衔接段漏热带漏热形成的局部热分层现象的抑制和消除,并研究了不同低温射流条件对于消除零重力条件下液氢储罐内部热分层效果的影响. 研究结果表明对于大尺寸储罐,当采用圆形射流喷嘴且低温射流条件相同时,射流喷嘴的位置对罐体内部热分层消除效果影响不是很明显. 当射流喷嘴位于储罐内部同一相对位置且入射流量相同时,圆形射流喷嘴因出流方向更集中,罐内流场演变更快,消除热分层比半球形射流喷嘴更有效.      

入口速度比对射流拟序结构的影响

为了深入了解湍流流动机理以及湍流拟序结构发现过程的影响因素，本文采用大涡模拟方法对不同入口射流伴流速度比的平面湍射流流动进行了数值模拟。采用分步投影法求解动量方程，亚格子项采用标准Smagorinsky亚格子模式模拟，压力泊松方程采用修正的循环消去法快速求解，空间方程采用二阶精度的差分格式，在时间方向上采用二阶精度的显式差分格式。模拟结果给出了平面射流中湍流拟序结构的瞬态发展演变过程，分析了入口速度比对射流拟序结构发展演化过程及宏观流场形态的影响。为进一步研究射流拟序结构及其在湍流流动中的作用提供了基础。     

基于ALE,CEL和SPH方法的球形破片高速冲击充液结构对比研究

为研究ALE,CEL和SPH方法在高速冲击流固耦合动力学数值分析中的差异性,开展球形破片高速冲击充液结构数值模拟研究。建立经文献资料验证的ALE,CEL和SPH三种动力学模型,研究了流体压力变化、形成的空腔尺寸、破片速度衰减变化和充液结构变形等模拟精度,并分析相应的计算成本。结果表明,ALE,CEL和SPH三种方法均能有效模拟破片高速冲击充液结构的流固耦合动力学过程;ALE方法预测的空腔尺寸精度较高;CEL方法预测的流体压力、破片速度衰减和充液结构变形精度较高;SPH方法预测的空腔尺寸、破片速度衰减精度较高;当网格尺寸一致时,SPH方法计算时长约为ALE和CEL方法的两倍,但SPH方法前后处理更加简便。     

大速差射流预燃室内三维回流两相湍流的数值模拟

本文由多流体两相流模型、气相湍流κ-ε模型和颗粒湍流代数模型出发,成功地模拟了真实形状大速差射流预燃室中三维湍流回流两相流动,得到了这类复杂的气固两相流中不同纵横截面上气相速度场、颗粒速度场及浓度场和两相湍流度场的分布,并且获得了与实验定性一致的合理结果,揭示了预燃室中气固两相流动与混合的主要物理特征,探讨了大速差射流技术稳焰和强化燃烧的两相流动机理。     

出口条件对圆形湍流射流远场区自保持性的影响

报道出口条件对圆形湍流射流自保持性影响的实验研究结果. 对来自渐缩和长管两种不同结构喷嘴的射流,在相同雷诺数条件下,沿轴线进行了速度测量; 研究的统计量包括平均速度、湍流强度、高阶矩、能谱和积分尺度. 实验结果表明,渐缩喷嘴射流比长管射流发展得更快、更容易达到自保持状态. 通过对比发现,在两射流的速度(温度)场中,平均速度(温度)、湍流强度、偏斜因子和平坦度因子都存在明显的异同. 同时发现两射流的积分尺度随轴向距离的增加都成线性增长,且在渐缩喷嘴射流中增长得更快. 通过对比两射流的边界层厚度、径向与轴向湍流强度的比值、湍动能能谱图并结合前人的研究结果,对两射流湍流场所表现出的不同的统计学行为给出了合理的解释.      

大速差射流预燃室内三维回流两相湍流的数值模拟

本文由多流体两相流模型、气相湍流κ-ε模型和颗粒湍流代数模型出发,成功地模拟了真实形状大速差射流预燃室中三维湍流回流两相流动,得到了这类复杂的气固两相流中不同纵横截面上气相速度场、颗粒速度场及浓度场和两相湍流度场的分布,并且获得了与实验定性一致的合理结果,揭示了预燃室中气固两相流动与混合的主要物理特征,探讨了大速差射流技术稳焰和强化燃烧的两相流动机理。     

围压下自振空化射流冲蚀性能实验研究

自振空化射流是利用瞬态流和水声学原理调制而成的一种新型射流,为研究围压下自激振荡空化射流的冲蚀破碎规律,利用高压釜装置测量了1.0mm出口直径的风琴管自振空化喷嘴在各种射流参数情况下冲蚀铝试样的冲蚀质量,并与同等条件下锥形喷嘴冲蚀效率进行了对比。测量结果表明,冲蚀质量基本与射流压力成正比;存在最优喷距和围压,使得冲蚀效果最佳,在本实验条件,分别为喷嘴出口直径的5～7倍和2MPa左右;相同条件下,自振空化喷嘴冲蚀质量约为同等条件下锥形喷嘴冲蚀的1～2倍,这为自振空化射流提高钻井速度等实际应用提供了实验依据。     

喷嘴内部流道型线对射流流场的影响

为研究喷嘴内部流道型线对射流流场的气体动特性参数的影响,根据可压缩流体轴对称N－S方程,利用RNGκ－ε湍流模式和有限体积法,并采用四边形非结构网格,对不同内部流道型线的喷嘴自由射流进行数值模拟。对于轴对称等直径圆管喷嘴,进口处的流道型线对射流流道参数的分布有较大的影响;轴对称收缩喷嘴的收缩角大小主要影响射流出口附近的流动,对流动具有不同的阻滞效果。计算结果与实验吻合较好,若要获得较佳的外部流场参数,优化喷嘴内部流道设计十分重要。     

Forced convective film condensation inside vertical tubes

A theoretical study of forced convective film condensation inside vertical tubes is presented. We propose a unified procedure for predicting the pressure gradient and condensation heat transfer coefficient of a vapor flowing turbulently in the core and associated with laminar or turbulent film on the tube wall. The analysis for the vapor flows is performed under the condition that the velocity profiles are locally self-similar. The laminar and turbulent film models equate the gravity, pressure and viscous forces, and consider the effect of interfacial shear. The transition from laminar to turbulent film depends not only on the liquid Reynolds number but also on the interfacial shear stress. In this work we also proposed a new eddy viscosity model which is divided into three regions: the inner region in liquid condensate near the wall; the interface region including both liquid and vapor; and the outer region for the vapor core. Comparisons of the theory with some published experimental data showed good agreement.     

Computations for a jet impinging obliquely on a flat surface

A SIMPLE-C algorithm and Jones-Launder k-ε two-equation turbulence model are used to simulate a two-dimensional jet impinging obliquely on a flat surface. Both the continuity and momentum equations for the unsteady state are cast into suitable finite difference equations. The pressure, velocity, turbulent kinetic energy and turbulent energy dissipation rate distributions are solved and show good agreement with various experimental data. The calculations show that the flow field structure of the jet impinging obliquely on a flat surface is strongly affected by the oblique impingement angle. The maximum pressure zone of the obliquely impinging jet flow field moves towards the left as the oblique impingement angle is decreased.     

Turbulent flow structure and bubble distribution in an axisymmetric nonisothermal impinging gas-liquid jet

The flow structure of a bubbly impinging jet in the presence of heat transfer between the two-phase flow and the surface is numerically investigated on the basis of the Eulerian approach. The model uses the system of Reynolds-averaged Navier–Stokes equations in the axisymmetric approximation written with account for the inverse effect of the bubbles on the average and fluctuating flow parameters. The influence of the gas volumetric flow rate ratio and the dimensions of the bubbles on the flow structure in a gas-liquid impinging jet is studied, In the presence of gas bubbles the liquid velocity is higher than the corresponding value in the single-phase flow. A considerable, more than twofold, anisotropy between the axial and radial turbulent fluctuations in the gas-liquid impinging jet is shown to exist. An addition of air bubbles leads to a considerable growth in the liquid velocity fluctuations in the two-phase flow (up to 50% compared with the single-fluid liquid impinging jet). An increase in the disperse phase dimensions leads to intensification of turbulence of the liquid.     

Experimental study of an impinging jet with different swirl rates

A stereo PIV technique using advanced pre- and post-processing algorithms is implemented for the experimental study of the local structure of turbulent swirling impinging jets. The main emphasis of the present work is the analysis of the influence of swirl rate on the flow structure. During measurements, the Reynolds number was 8900, the nozzle-to-plate distance was equal to three nozzle diameters and the swirl rate was varied from 0 to 1.0. For the studied flows, spatial distributions of the mean velocity and statistical moments (including triple moments) of turbulent pulsations were measured.

The influence of the PIV finite spatial resolution on the measured dissipation rate and velocity moments was analyzed and compared with theoretical predictions. For this purpose, a special series of 2D PIV measurements was carried out with vector spacing up to several Kolmogorov lengthscales.

All terms of the axial mean momentum and the turbulent kinetic energy budget equations were obtained for the cross-section located one nozzle diameter from the impinging plate. For the TKE budget, the dissipation term was directly calculated from the instantaneous velocity fields, thereby allowing the pressure diffusion term to be found as a residual one. It was found that the magnitude of pressure diffusion decreased with the growth of the swirl rate. In general, the studied swirling impinging jets had a greater spread rate and a more rapid decay in absolute velocity when compared to the non-swirling jet.     

Measurements of turbulent flow in axially finned rod bundles

The hydraulic characteristics in subchannels of axially finned rod bundles installed in the Korea Multipurpose Research Reactor (KMRR) were measured using one-component laser-Doppler velocimetry (LDV). Pressure drops for each component, time mean axial velocity, and axial turbulent intensity were measured. Then the friction factors in rod bundles were estimated from the measured pressure drops. The turbulent crossflow mixing rate between neighboring subchannels was evaluated from the measured data. The results show that (1) the friction factors for axially finned rod bundles are less than those given by Moody's correlation for smooth tubes, and (2) as the flow develops, the turbulent crossflow mixing rate between neighboring subchannels decreases and, in the developed region, the values level out.     

压力梯度对湍流逆梯度输运的影响

研究Birkhoff系统的一般Lie对称性导致的非Noether守恒量. 得到非Noether守恒 量的存在定理，举例说明结果的应用.     

On the Capability of PIV-Based Wall Pressure Estimation for an Impinging Jet Flow

A PIV-based pressure estimation methodology is used to compute the wall pressure from the velocity field of a turbulent impinging jet flow. A simplified formulation (2D-2C) is applied to velocity fields issued from PIV data. The ability of the method to qualitatively estimate the wall pressure signature of a 3D unsteady impinging jet flow using only two velocity components in a plane is demonstrated. Nevertheless, the 2D flow assumption used in the context of planar measurements involves an underestimation of the wall pressure values all along the radial direction. The formulation based on the full integral formalism (3D-3C), computed from DNS data without any assumption on the flow, provides a reference solution. The contributions of the surface and volume integrals to the pressure coefficient are assessed. It is shown that the most important contribution to the wall pressure comes from the volume integral. Then the underestimation observed for the simplified formulation is mostly linked with the assumptions considered for the source term computation. The effect of each assumption is quantitatively analysed with the help of the DNS data and some ways to improve the simplified methodology are finally proposed.     

Investigation of the turbulent model for pressure fluctuations with spectral theory

The pressure fluctuations in turbulent shear flows are investigated with the theory of spectral analysis.An expression for pressure spectra is analytically derived in terms of velocity spectra.This derivation is based on a formal solution of the Navier-Stokes equation and quasi-normal assumption to express the third and fourth order velocity correlations in terms of double velocity correlation.Then,a turbulent model for the computation of pressure fluctuation intensity with Renolds stress and mean flow velocity gradients is established.The turbulent constants in the model are calculated from the assumptions about the general behaviour of velocity spectra in high Renolds number flows.Comparison with direct simulation of turbulent boundary layer is made.It is found that the turbulent-turbulent,cross correlation,and turbulent-shear source terms for mean square value of pressure fluctuation are about the same magnitude.     

Measurements in the vicinity of a stagnation point

This paper presents measurements of a plane jet impinging onto a normal flat plate placed up to five jet widths from the jet outlet. The small spacing ensured that the stagnation streamline remained in the potential core of the jet. The plate shear stress distribution compared well to that from an analytical solution for the laminar development of the plate boundary layer whose external velocity was determined from the measured pressure. By comparing the shear stress measured under the present low level of free stream turbulence (0.35%) at the jet exit with that of Tu and Wood [Exp. Thermal Fluid Sci. 13 (1996) 364–373] made at about 4%, it is concluded that the turbulence level at the nozzle exit has only a second-order influence on the surface shear stress around the stagnation point. Some spanwise non-uniformity was observed in the plate shear stress, but this was confined largely to the transition region. The mean velocity, Reynolds stresses, and fluctuating pressure were measured along the stagnation streamline using a fast-response pressure probe. A significant increase in the streamwise normal stress and the mean square of the pressure fluctuations occurred before they were eventually attenuated by the plate. This increase occurred in the region where the streamwise velocity was decreasing close to the plate causing extra energy production through the normal stresses. Spectra of the velocity and pressure fluctuations showed that the increase in level was mainly due to the low frequency motion, whereas the subsequent decrease occurred at higher frequencies.     

A parametric study of adverse pressure gradient turbulent boundary layers

There are many open questions regarding the behaviour of turbulent boundary layers subjected to pressure gradients and this is confounded by the large parameter space that may affect these flows. While there have been many valuable investigations conducted within this parameter space, there are still insufficient data to attempt to reduce this parameter space. Here, we consider a parametric study of adverse pressure gradient turbulent boundary layers where we restrict our attention to the pressure gradient parameter, β, the Reynolds number and the acceleration parameter, K. The statistics analyzed are limited to the streamwise fluctuating velocity. The data show that the mean velocity profile in strong pressure gradient boundary layers does not conform to the classical logarithmic law. Moreover, there appears to be no measurable logarithmic region in these cases. It is also found that the large-scale motions scaling with outer variables are energised by the pressure gradient. These increasingly strong large-scale motions are found to be the dominant contributor to the increase in turbulence intensity (scaled with friction velocity) with increasing pressure gradient across the boundary layer.