气泡-液体闭式多股射流的大涡模拟研究

采用气泡-液体两相流动的欧拉-拉氏大涡模拟,研究了矩形通道内多股射流形成的气泡-液体两相湍流流动,得到了气液两相湍流瞬态结构,产生和发展过程。研究结果发现气泡和液体都有瞬态大涡结构,气泡脉动比液体的强。大涡模拟统计结果给出了有无气泡两种情况下的液体湍流脉动速度均方根值分布。瞬态和统计结果都表明,气泡增大了液体湍流,液体湍流来源于其自身的剪切产生和气泡的作用。这与二阶矩模型模拟结果定性一致。     

A pressure-based algorithm for multi-phase flow at all speeds

A new finite volume-based numerical algorithm for predicting incompressible and compressible multi-phase flow phenomena is presented. The technique is equally applicable in the subsonic, transonic, and supersonic regimes. The method is formulated on a non-orthogonal coordinate system in collocated primitive variables. Pressure is selected as a dependent variable in preference to density because changes in pressure are significant at all speeds as opposed to variations in density, which become very small at low Mach numbers. The pressure equation is derived from overall mass conservation. The performance of the new method is assessed by solving the following two-dimensional two-phase flow problems: (i) incompressible turbulent bubbly flow in a pipe, (ii) incompressible turbulent air–particle flow in a pipe, (iii) compressible dilute gas–solid flow over a flat plate, and (iv) compressible dusty flow in a converging diverging nozzle. Predictions are shown to be in excellent agreement with published numerical and/or experimental data.     

Thermal and efficiency analysis of the CO2 laser cutting process

In CO₂ laser gas-assisted cutting process, modeling of the interaction mechanism is important. Consequently, the present study treats the complete problem of the interaction of the melting surface with the boundary layer and describes the behavior of the melting layer. In the analysis, gas–liquid interface parameters are developed and relationships between the parameters influencing the cutting action are identified theoretically. To achieve this, effects of momentum and gas–liquid interface shear stress, due to the assisting gas jet, are considered. The approximate magnitude of the heat absorbed is estimated and melting layer thickness is predicted. An experiment is carried out and the theoretical predictions are compared with the experimental findings. First and second law efficiencies of the cutting process are predicted, which may, then, be used to improve the process. It is found that the assisting jet velocity increases the first and second law efficiencies of the CO₂ laser cutting process.     

Characterization of chaotic dynamic behavior in the gas–liquid slug flow using directed weighted complex network analysis

Understanding the nonlinear and complex dynamics underlying the gas–liquid slug flow is a significant but challenging problem. We systematically carried out gas–liquid two-phase flow experiments for measuring the time series of flow signals, which is studied in terms of the mapping from time series to complex networks. In particular, we construct directed weighted complex networks (DWCN) from time series and then associate different aspects of chaotic dynamics with the topological indices of the DWCN and further demonstrate that the DWCN can be exploited to detect unstable periodic orbits of low periods. Examples using time series from classical chaotic systems are provided to demonstrate the effectiveness of our approach. We construct and analyze numbers of DWCNs for different gas–liquid flow patterns and find that our approach can quantitatively distinguish different experimental gas–liquid flow patterns. Furthermore, the DWCN analysis indicates that slug flow shows obvious chaotic behavior and its unstable periodic orbits reflect the intermittent quasi-periodic oscillation behavior between liquid slug and large gas slug. These interesting and significant findings suggest that the directed weighted complex network can potentially be a powerful tool for uncovering the underlying dynamics leading to the formation of the gas–liquid slug flow.     

垂直上升管泡状流压力波动的多尺度分析

从理想的单气泡运动物理模型和离散正交小波与统计相结合的信号处理技术研究了泡状流压力波动的机理与多尺度时频规律,得到以下结论:泡状流的压力波动主要来源于气泡的运动和诱导湍流及液相湍流,其中在低流速下前者产生的压力波动幅度大于后两者;存在一个临界频率,低于该频率压力信号均方根随频率的增加而增加,高于该频率压力信号均方根随频率的增加而减小;发生流型转变时压力信号的时频特征发生明显改变,最大均方根出现在更低的频率区。     

气泡-液体多股射流二阶矩两相湍流模型的模拟

本文采用二阶矩两相湍流模型模拟了气泡－液体闭式多股射流。研究结果显示出预报的两相速度、雷诺应力和气泡体积分数分布的合理性，以及在速度较高和剪切较强的情况下，气泡湍流脉动仍然大于液体湍流脉动，气泡增强液体湍流，气液两相湍流脉动均存在各向异性，而且气泡脉动的各向异性比液体脉动的大等规律。     

Simulation of turbulent non-isothermal polydisperse bubbly flow behind a sudden tube expansion

The results of numerical simulation of the structure of non-isothermal polydisperse bubbly turbulent flow and heat transfer behind a sudden tube expansion are presented. The study was carried out at a change in the initial diameter of the air bubbles within d _m1 = 1–5 mm and their volumetric void fraction β = 0–10 %. Small bubbles are available in almost the entire cross section of the tube, while the large bubbles pass mainly through the flow core. An increase in the size of dispersed phase causes the growth of turbulence in the liquid phase due to flow turbulization, when there is a separated flow of liquid past the large bubbles. Adding the air bubbles causes a significant reduction in the length of the separation zone and heat transfer enhancement, and these effects increase with increasing bubble size and their gas volumetric flow rate ratio.     

Stochastic modeling of atomizing spray in a complex swirl injector using large eddy simulation

Large-eddy simulation of an atomizing spray issuing from a gas-turbine injector is performed. The filtered Navier–Stokes equations with dynamic subgrid scale model are solved on unstructured grids to compute the swirling turbulent flow through complex passages of the injector. The collocated grid, incompressible flow algorithm on arbitrary shaped unstructured grids developed by Mahesh et al. (J. Comp. Phys. 197 (2004) 215–240) is used in this work. A Lagrangian point-particle formulation with a stochastic model for droplet breakup is used for the liquid phase. Following Kolmogorov’s concept of viewing solid particle-breakup as a discrete random process, the droplet breakup is considered in the framework of uncorrelated breakup events, independent of the initial droplet size. The size and number density of the newly produced droplets is governed by the Fokker–Planck equation for the evolution of the pdf of droplet radii. The parameters of the model are obtained dynamically by relating them to the local Weber number and resolved scale turbulence properties. A hybrid particle-parcel is used to represent the large number of spray droplets. The predictive capability of the LES together with Lagrangian droplet dynamics models to capture the droplet dispersion characteristics, size distributions, and the spray evolution is examined in detail by comparing it with the spray patternation study for the gas-turbine injector. The present approach is computationally efficient and captures the global features of the fragmentary process of liquid atomization in complex configurations.     

On the generation and maintenance of waves and turbulence in simulations of free-surface turbulence

One of the challenges in numerical simulation of wave–turbulence interaction is the precise setup and maintenance of wave and turbulence fields. In this paper, we investigate techniques for the generation and suppression of specific surface wave modes, the generation of turbulence in an inhomogeneous physical domain with a wavy boundary-fitted grid, and the generation and maintenance of waves and turbulence during the complex wave–turbulence interaction process. We apply surface pressure to generate and suppress waves. Based on the solution of linearized Cauchy–Poisson problem, we derive three pressure expressions, which lead to a δ-function method, a time-segment method, and a gradual method. Numerical experiments show that these methods generate waves as specified and eliminate spurious waves effectively. The nonlinear wave effect is accounted for with a time-relaxation method. For turbulence generation, we extend the linear forcing method to an inhomogeneous physical domain with a curvilinear computational grid. Effects of force distribution and computational grid distortion are examined. For wave–turbulence interaction, we develop an algorithm to instantaneously identify specific progressive and standing waves. To precisely control the wave amplitude in a complex turbulent flow field, we further develop an energy controlling method. Finally, a simulation example of wave–turbulence interaction is presented. Results show that turbulence has unique features in the presence of waves. Velocity fluctuations are found to be strongly dependent on the wave phase; variations of these fluctuations are explained by the pressure–strain correlation associated with the wave-induced strain field.     

数值模拟两相汽蚀流动的新模型和算法

提出了数值模拟两相汽蚀流动现象的新汽蚀模型和算法.提出的汽蚀模型和算法耦合了考虑紊流粘性效应的Reynolds-Averaged Navier-Stokes方程求解方法,可以自动模拟空泡起始点、空泡长度和汽蚀空泡形状.在流场计算和界面修正的迭代计算过程中,跟踪并得到液相/气相界面.数值模拟了圆锥形圆柱体和半球形圆柱体在不同汽蚀系数下的汽蚀流动现象,得到了与实验值完全吻合的数值模拟结果,并且与已发表的数值研究结果进行了比较.计算结果表明提出的汽蚀模型和算法能够有效地模拟汽蚀流动中的气泡界面和气泡长度,汽蚀模型和算法的正确性和实用性得到了相应的验证.     

An experimental investigation of the propagation mechanism of critical deflagration waves that lead to the onset of detonation

The reflection of a CJ detonation from a perforated plate is used to generate high speed deflagrations downstream in order to investigate the critical conditions that lead to the onset of detonation. Different perforated plates were used to control the turbulence in the downstream deflagration waves. Streak Schlieren photography, ionization probes and pressure transducers are used to monitor the flow field and the transition to detonation. Stoichiometric mixtures of acetylene–oxygen and propane–oxygen were tested at low initial pressures. In some cases, acetylene–oxygen was diluted with 80% argon in order to render the mixture more “stable” (i.e., more regular detonation cell structure). The results show that prior to successful detonation initiation, a deflagration is formed that propagates at about half the CJ detonation velocity of the mixture. This “critical” deflagration (which propagates at a relatively constant velocity for a certain duration prior to the onset of detonation) is comprised of a leading shock wave followed by an extended turbulent reaction zone. The critical deflagration speed is not dependent on the turbulence characteristics of the perforated plate but rather on the energetics of the mixture like a CJ detonation (i.e., the deflagration front is driven by the expansion of the combustion products). Hence, the critical deflagration is identified as a CJ deflagration. The high intensity turbulence that is required to sustain its propagation is maintained via chemical instabilities in the reaction zone due to the coupling of pressure fluctuations with the energy release. Therefore, in “unstable” mixtures, critical deflagrations can be supported for long durations, whereas in “stable” mixtures, deflagrations decay as the initial plate generated turbulence decays. The eventual onset of detonation is postulated to be a result of the amplification of pressure waves (i.e., turbulence) that leads to the formation of local explosion centers via the SWACER mechanism during the pre-detonation period.     

Study on bubble formation in non-Newtonian fluids by laser image technique

A laser image system for studying bubble formation at the orifice submerged in liquid was established. The process of bubble formation can be directly visualized and real-time recorded through computer by means of He–Ne laser as light source using the beam expanding and light amplification technique. The behaviors of bubble formation in polyacrylamide (PAM) aqueous solutions were investigated at temperature 293.15 K and orifice diameters 1, 1.5 and 2 mm, respectively, the chamber volume was 90 cm³ and the gas flowrate from 0.1 to 0.6 cm³/s. The influences of mass concentration of solution and orifice diameter on bubble detachment volume were investigated experimentally and the results show that bubble detachment volume increases with solution concentration and orifice diameter in the experimental range of this work.     

Speckle photography applied to statistical analysis of turbulence

A speckle photographic technique is used for visualizing the planar distribution of the refractive deflection angles of light transmitted through the compressible turbulent flow. Both double and single (prolonged) exposure speckle photography are applied for statistical analysis of such flows. Using single (prolonged) exposure speckle photography (SPESP), instantaneous quantitative measurement of 2-D distribution of turbulence intensity in a flame is performed. Anisotropy of turbulence is visualized by a diffraction halo form and quantitatively evaluated by measuring the diffraction halo diameters. Using double exposure speckle photography (DESP), quantitative visualization of the planar distribution of the refractive deflection angles of the light transmitted through the compressible turbulent flow is done. Turbulent structures are visible in the patterns of the deflection angles isolines. The 2-D correlation functions of these deflection angles are constructed and analyzed. The 3-D density correlation functions are evaluated using the Erbeck–Merzkirch integral transformation.     

Centric scan SPRITE for spin density imaging of short relaxation time porous materials

The single-point ramped imaging with T1 enhancement (SPRITE) imaging technique has proven to be a very robust and flexible method for the study of a wide range of systems with short signal lifetimes. As a pure phase encoding technique, SPRITE is largely immune to image distortions generated by susceptibility variations, chemical shift and paramagnetic impurities. In addition, it avoids the line width restrictions on resolution common to time-based sampling, frequency encoding methods. The standard SPRITE technique is however a longitudinal steady-state imaging method; the image intensity is related to the longitudinal steady state, which not only decreases the signal-to-noise ratio, but also introduces many parameters into the image signal equation. A centric scan strategy for SPRITE removes the longitudinal steady state from the image intensity equation and increases the inherent image intensity. Two centric scan SPRITE methods, that is, Spiral-SPRITE and Conical-SPRITE, with fast acquisition and greatly reduced gradient duty cycle, are outlined. Multiple free induction decay (FID) points may be acquired during SPRITE sampling for signal averaging to increase signal-to-noise ratio or for T2* and spin density mapping without an increase in acquisition time. Experimental results show that most porous sedimentary rock and concrete samples have a single exponential T2* decay due to susceptibility difference-induced field distortion. Inhomogeneous broadening thus dominates, which suggests that spin density imaging can be easily obtained by SPRITE.     

Laser spectroscopic investigation of Ca atoms in liquid helium

Excitation and emission spectra of the 4s21S0 – 4s4p1P1 and 4s4p3PJ(J = 0, 1, 2) – 4s5s3S1 transitions of Ca atoms implanted not only in liquid 4He but also in liquid 3He have been measured. It has been found that the excitation spectra for liquid 3He show considerably smaller widths and peak shifts from transition wavelengths of a free Ca atom, in comparison with those for liquid 4He. These spectral characteristics have been well reproduced by our theoretical calculation based on a vibrating bubble model. This calculation has shown that, due to a significant difference in surface tension between the two liquids, the radius of a bubble formed around Ca in liquid 3He is larger than the one in 4He, and that this fact as well as the smaller number density of He atoms in liquid 3He cause weaker perturbation for Ca, resulting in the smaller peak shift and width for liquid 3He.     

I. Formation and rise of a bubble stream in a viscous liquid

The continuous emission of gas bubbles from a single ejection orifice immersed in a viscous fluid is considered. We first present a semi empirical model of spherical bubble growth under constant flow conditions to predict the bubble volume at the detachment stage. In a second part, we propose a physical model to describe the rise velocity of in-line interacting bubbles and we derive an expression for the net viscous force acting on the surrounding fluid. Experimental results for air/water-glycerol systems are presented for a wide range of fluid viscosity and compared with theoretical predictions. An imagery technique was used to determine the bubble size and rise velocity. The effects of fluid viscosity, gas flow rate, orifice diameter and liquid depth on the bubble stream dynamic were analyzed. We have further studied the effect of large scale recirculation flow and the influence of a neighbouring bubble stream on the bubble growth and rising velocity. Received: 23 July 1997 / Revised: 16 December 1997 / Accepted: 11 May 1998     

A front-tracking/ghost-fluid method for fluid interfaces in compressible flows

A front-tracking/ghost-fluid method is introduced for simulations of fluid interfaces in compressible flows. The new method captures fluid interfaces using explicit front-tracking and defines interface conditions with the ghost-fluid method. Several examples of multiphase flow simulations, including a shock–bubble interaction, the Richtmyer–Meshkov instability, the Rayleigh–Taylor instability, the collapse of an air bubble in water and the breakup of a water drop in air, using the Euler or the Navier–Stokes equations, are performed in order to demonstrate the accuracy and capability of the new method. The computational results are compared with experiments and earlier computational studies. The results show that the new method can simulate interface dynamics accurately, including the effect of surface tension. Results for compressible gas–water systems show that the new method can be used for simulations of fluid interface with large density differences.     

激光空泡在文丘里管中运动的动力学特性

以水为工作介质,在不同文丘里管入口压力下,利用YAG激光器产生的激光轰击水中的金属靶材产生空泡,借助高速摄像系统记录激光空泡在文丘里管中的运动过程,并采用流体动力学模拟对文丘里管中的流场特性和空泡的溃灭特性进行分析.结果表明:激光空泡在文丘里管中的运动,其形状变化可分为产生阶段、挤压阶段、溃灭初始阶段和溃灭阶段等四个阶段.空泡的溃灭取决于流场状态,当流动为层流时,空泡不发生溃灭;当流动为湍流时,空泡发生溃灭,且湍流程度越剧烈,溃灭现象越显著.搭建的激光空泡生成与运动系统能够实现空泡的定点溃灭.     

声场中水力空化泡的动力学特性

以水为工作介质,考虑了液体黏性、表面张力、可压缩性及湍流作用等情况,对文丘里管反应器中空化泡在声场作用下的动力学行为特性进行了数值研究.分析了超声波频率、声压及喉径比对空化泡运动特性以及空化泡崩溃时所形成泡温以及压力脉冲的影响.结果表明,超声将水力空化泡运动调制成稳态空化,有利于增强空化效果. 关键词： 超声波 水力空化 湍流 气泡动力学     

不同排列粗糙元对湍流拟序结构的影响

采用氢气泡流动显示技术,研究了壁面不同排列粗糙元对近壁湍流拟序结构的影响。实验中基于平均速度和水力直径的雷诺数分别为14300、31200、48000。氢泡丝上游布置的离散粗糙元直径为2．3、4和6mm,得到了光滑壁面和5种不同排列方式的离散粗糙元壁面湍流条带和湍流斑块特征。结果发现：相同雷诺数和相同排列方式下条带间距...