垂直上升气液环状流流动参数的数值计算

提出了一个新的气核-液膜耦合模型来求解垂直上升气液环状流在充分发展段的流动参数.本模型考虑了液膜、气核以及它们之间的相互影响和作用.模型中基本的气核区域和液膜区域的质量和动量方程由Fluent6.3.26进行求解,而液滴方程以及相界面上的夹带和沉积作用通过用户自定义接口函数UDF(User Defined Functi...     

Interpretation of atomization rates of the liquid film in gas-liquid annular flow

For vertical gas-liquid annular flow the fraction of the liquid in the gas is controlled by the rate of atomization of the liquid film flowing along the wall and the rate of deposition of droplets entrained in the gas. Measurements of the rate of atomization are interpreted by a Kelvin-Helmholtz mechanism. Small wavelets on the liquid film are visualized to be entrained when wave-induced variations in the gas pressure cannot be counterbalanced by surface tension effects.     

Hillslope soil erosion and runoff model for natural rainfall events

By using the momentum theorem and waterbalance principle, basic equations of slope runoff were derived, soil erosion by raindrop splash and runoff were discussed and a model was established for decribing hillslope soil erosion processes. The numerical solution of the model was obtained by adopting the Preissmann format and considering the common solution-determining conditions, from which not only the runoff and soil erosion but also their processes can be described. The model was validated by ten groups of observation data of Soil Conservation Ecological Science and Technology Demonstration Park of Jiangxi Province. Comparisons show that the maximum relative error between simulation and experimental data is about 10.98% for total runoff and 15 % for total erosion, 5.2% for runoffprocess and 6.1% for erosion process, indicating that the model is conceptually realistic and reliable and offers a feasible approach for further studies on the soil erosion process.     

The onset of gas entrainment from a flowing stratified gas–liquid regime in dual discharging branches: Part II: Critical conditions at low to moderate branch Froude numbers

This is the second of a two part investigation. Experiments were performed in a 50.8 mm diameter horizontal pipe with three 6.35 mm diameter branches located at the test section mid-span. The inlet length was 1.8 m, and three branch orientations were tested at 0° (side), 45° (inclined), and 90° (bottom) from horizontal. Water and air, operating at 206 kPa, were used and both fluids flowed co-currently within the inlet in the smooth-stratified regime. The inlet superficial velocity of the liquid phase ranged between 0.04 and 0.15 m/s while in the gas phase values of 0.3, 0.4, and 1 m/s were tested. Three different dual discharge combinations were tested and included side-inclined, side-bottom, and inclined-bottom. The tested branch flow Froude numbers were limited between low to moderate values which ranged between 1 and 23. Extensive experimental data are reported for the critical conditions at the onset of gas entrainment during dual discharge. A novel map was developed for the inclined-bottom branch configuration showing the relationship between the inlet superficial liquid velocity and branch Froude numbers. This map was used to quantify the three observed modes of gas entrainment during dual discharge. These modes were classified as onset of gas entrainment in the inclined branch only, in the bottom branch only, or both the inclined and bottom branches simultaneously. The critical height at the onset of gas entrainment results are compared to published models and data sets and poor agreement was found with studies conducted in stratified gas–liquid reservoirs.     

PTV experiments of subcooled boiling flow through a vertical rectangular channel

Time resolved Particle Tracking Velocimetry (PTV) experiments were carried out to investigate turbulent, subcooled boiling flow of refrigerant HFE-301 through a vertical rectangular channel with one heated wall. Measurements were performed with liquid Reynolds numbers (based on the hydraulic diameter) of Re = 3309, 9929 and 16,549 over a wall heat flux range of 0.0–64.0 kW/m². Turbulence statistics are inferred from PTV full-field velocity measurements. Quantities such as: instantaneous 2D velocity fields, time-averaged axial and normal velocities, axial and normal turbulence intensities, and Reynolds stresses are obtained. The present results agree well with previous studies and provides new information due to the full-field nature of the technique. This work is an attempt to provide turbulent subcooled boiling flow data for validation and improvement of two-phase flow computational models.     

Film thickness measurements in liquid–liquid slug flow regimes

At present there is significant interest in the development of small scale medical diagnostic equipment. These devices offer faster processing times and require smaller sample volumes than equivalent macro scale systems. Although significant attention has been focused upon their outputs, little attention has been devoted to the detailed fluid mechanics that govern the flow mechanisms within these devices. Conventionally, the samples in these small scale devices are segmented into distinct discrete droplets or slugs which are suspended in an organic carrier phase. Separating these slugs from the channel wall is a very thin film of the organic carrier phase.The magnitude of this film is the focus of the present study and the effects of sample slug length and carrier phase fluidic properties on the film are examined over a range of Capillary numbers. A non-intrusive optical technique was used to capture images of the flow from which the magnitude of the film was determined.The experimental results show that the film is not constant along the length of the slug; however above a threshold value for slug length, a region of constant film thickness exists. When compared with existing correlations in the literature, the experimental data showed reasonable agreement with the Bretherton model when the Capillary number was calculated based on the mean two phase flow velocity. However, significant differences were observed when the Capillary number was redefined to account for the mean velocity at the liquid interface, i.e., the mean slug velocity.Analysis of the experimental data revealed that it fell into two distinct flow regimes; a visco-capillary regime and a visco-inertial regime. A modified Taylor expression is presented to estimate the magnitude of the film for flows in the visco-capillary regime while a new model is put forward, based on Capillary and Weber numbers, for flows in the visco-inertial regime. Overall, this study provides some novel insights into parameters, such as aqueous slug length and carrier phase fluidic properties, that affect the thickness of the film in liquid–liquid slug flow regimes.     

Simulation of multiphase flows in porous media

The ability to numerically simulate single phase and multiphase flow of fluids in porous media is extremely important in developing an understanding of the complex phenomena governing the flow. The flow is complicated by the presence of heterogeneities in the reservoir and by phenomena such as diffusion, dispersion, and viscous fingering. These effects must be modeled by terms in coupled systems of nonlinear partial differential equations which form the basis of the simulator. The simulator must be able to handle both single and multiphase flows and the transition regimes between the two. A discussion of some of the aspects of modeling dispersion and viscous fingering is presented along with directions for future work.The partial differential equation models are convection-dominated and contain important local effects. An operator-splitting technique is used to address these different effects accurately. Convection is treated by time stepping along the characteristics of the associated pure convection problem, and diffusion is modeled via a Galerkin method for single phase flow and a Petrov-Galerkin technique for multiphase regimes. ELLAM (Eulerian-Lagrangian Localized Adjoint Methods) are discussed to effectively treat the advection-dominated processes. Accurate approximations of the fluid velocities needed in the Eulerian-Lagrangian time-stepping procedure are obtained by mixed finite element methods. Adaptive local grid refinement techniques are then indicated to resolve important local phenomena around wells and large heterogeneities or to resolve the moving internal boundary layers which often govern the mass transfer between phases.     

Structure of thin liquid films in gas-liquid horizontal flow

The structure of the liquid film in horizontal annular flow is studied visually using the refractive index matching technique. The liquid film is found to contain significant amount of air bubbles, which are continuously entrained, broken up and released by the rolling motion within the film. A new conceptual picture of the gas-liquid interface is presented.     

Shock-induced Near-wall Two-phase Flow Structure Over a Micron-sized Particles Bed

The present paper studies numerical modelling of near-wall two-phase flows induced by a normal shock wave moving at a constant speed, over a micron-sized particles bed. In this two-fluid model, the possibility of particle trajectory intersection is considered and a full Lagrangian formulation of the dispersed phase is introduced. The finiteness of the Reynolds and Mach numbers of the flow around a particle as well as the fineness of the particle sizes are taken into account in describing the interactions between the carrier- and dispersed-phases. For the small mass-loading ratio case, the numerical simulation of flow structure of the two phases is implemented and the profiles of the particle number density are obtained under the constant-flux condition on the wall. The effects of the shock Mach number and the particle size and material density on particle entrainment motion are discussed in detail. The obtained results indicate that interphase non-equilibrium in the velocity and temperature is a common feature for this type of flows and a local particle accumulation zone may form near the envelope of the particle trajectory family.     

不同参数下环形通道内湍流旋流流动的模拟

应用一种合理考虑湍流一旋流相互作用及湍流脉动各向异性的新的代数ＲｅｙｎｏｌｄＳ应力模型,对环形通道内的湍流旋流流动进行了数值模拟．研究了旋流数、进口轴向速度和内外半径比等参数对环形通道内湍流旋流流动的影响,以及由此产生的流场变化对强化环形通道内传热的作用．     

Grid-averaged Lagrangian LES model for multiphase turbulent flow

A grid-averaged Lagrangian (GAL) model for dispersed particle motion in multiphase turbulent flow is presented to provide a large eddy simulation (LES) model for multiphase turbulent flow in which a quite large number of particles are involved. The GAL model is based on an averaging operation for a Lagrangian-type equation of motion of a particle over a computational grid volume and a procedure of reallocation of a dispersed particle cloud with its centroid movement to each grid. The model is therefore a mixed Eulerian–Lagrangian model which can effectively reduce computational time compared with existing Lagrangian-type models, without losing the advantage of Lagrangian-type models that they can properly describe the dynamical evolution of particles. Since the GAL model adopts the grid-volume averaging operation it can easily provide an effective SGS model for LES modeling of multiphase turbulent flow. The validity of the multiphase LES model developed, which is named the GAL-LES model, is confirmed through its application to a particle plume, in which the present model is found to simulate large-eddy motion usually observed in a jet and plume, and to give good agreements with experimental data.     

内管自转且公转时环管幂律流的受力分析

数值模拟了环管中内管偏心自转且公转时由轴向压力所驱动的幂律流体充分发展层流,分析了内管上的流体作用力。结果表明,内管偏心自转时流体作用力具有推动内管作和自转同向公转的效果。当只有外力矩驱动内管自转时,由于流体的作用,随内管线密度的不同,内管能达到的受力平衡态也不同:线密度较小时内管仅能在同心自转时达到受力平衡;线密度较大时内管能在作具有不变角速度和偏心率公转时达到受力平衡,且内管线密度越大,对应的受力平衡的公转的偏心率也越大。     

环形通道内湍流旋流流动的数值模拟

本文应用一种考虑湍流－旋流相互作用及湍流脉动各向异性的新的代数Ｒｅｙｎｏｌｄｓ应力模型,对环形通道内的湍流旋流流动进行了数值模拟,研究了改主为旋流流数,进口轴向速度及半径比等参数对环形通道内湍流流动的影响,以及对强化环形通道内传热的作用。     

Towards large eddy simulation of isothermal two-phase flows: Governing equations and a priori tests

This article reports on the potential of application of LES in the calculation of turbulent two-phase flows, in the case where each phase is resolved and interfaces remain much larger than the mesh size. In comparison with single-phase flow, successful application of LES to resolve two-phase flow problems should account for the complex interaction between turbulence and interfaces. Non-linear transfers of turbulent energy across the interface have to be accurately modeled. The derivation of the complete filtered two-phase flow governing equations has been formulated to deal with turbulence at the interface in a comprehensive and practical way. Explicit filtering of 2D direct numerical simulations has been employed to evaluate the order of magnitude of the new subgrid contributions. A parametric study on the academic test case of two counter-rotative vortices and a more complex test case of phase inversion in a closed box have been utilized to perform an order of magnitude analysis of different transport mechanisms. Important features of turbulent energy transfer across the interface have been discussed. Analyses of the numerical results have been conducted to derive conclusions on the relative importance of the different subgrid scale contributions, and modeling issues and solutions are provided.     

Numerical simulation of the modulated flow pattern for vertical upflows by the phase separation concept

The multiphase heat transfer could be enhanced by creating thin liquid film on the wall. The phase separation concept is called due to the separated flow paths of liquid and gas over the tube cross section to yield thin liquid film. Our proposed heat transfer tube consists of an annular region close to the wall and a core region, interfaced by a suspending mesh cylinder in the tube. The heat transfer tube is a multiscale system with micron scale of mesh pores, miniature scale of annular region and macroscale of tube diameter and length. Great effort has been made to link from micron scale to macroscale. The Volume of Fluid (VOF) method simulates air/water two-phase flow for vertical upflow. The three-dimensional system was successfully converted to a two-dimensional one by using three equivalent criteria for mesh pores. The non-uniform base grid generation and dynamic grid adaption method capture the bubble interface. The numerical results successfully reproduce our experimental results. The numerical findings identify the following mechanisms for the enhanced heat transfer: (a) counter-current flow exists with upward flow in the annular region and downward flow in the core region; (b) void fractions are exact zero in the core region and higher in the annular region; (c) the liquid film thicknesses are decreased to 1/6–1/3 of those in the bare tube section; (d) the gas–liquid mixture travels much faster in the annular region than in the bare tube; (e) three-levels of liquid circulation exists: meter-scale bulk liquid circulation, moderate-scale liquid circulation around a single-elongated-ring-slug-bubble, and microliquid circulation following the ring-slug-bubble tails. These liquid circulations promote the fluid mixing over the whole tube length and within the radial direction. The modulated parameters of void fractions, velocities and liquid film thicknesses in the annular region and three-levels of liquid circulation are greatly beneficial for the multiphase heat transfer enhancement.     

Study of gas-sheared liquid film in horizontal rectangular duct using high-speed LIF technique: Three-dimensional wavy structure and its relation to liquid entrainment

The flow of a liquid film sheared by high velocity gas stream in a horizontal rectangular duct was investigated using a high-speed laser-induced fluorescence technique. Measurements of local film thickness were resolved in both longitudinal and transverse coordinates with high spatial and temporal resolution. It was found that the generation of fast and slow ripples by the disturbance waves was qualitatively the same as it was observed earlier in completely different conditions. The transverse size and curvature of the disturbance waves and ripples were measured. A relationship between the three-dimensional structure of ripples on top of disturbance waves and the two mechanisms of liquid entrainment, known as ‘bag break-up’ and ‘ligament break-up’, is proposed.     

Simulations of dispersed turbulent multiphase flow

Direct numerical simulations of homogeneous isotropic turbulence are used to investigate the effects of turbulence on the transport of particles in gas flows or bubbles in liquid flows. The inertia associated with the bubbles or the particles leads to locally strong concentrations of these in regions of instantaneously strong vorticity for bubbles or strain-rate for particles. This alters the average settling rates and other processes. If the mass-loading of the dispersed phase is significant a random “turbulent” flow is generated by the particle settling. A simple demonstration of this is given, showing the statistically axisymmetric character of this flow and how it can modify an ambient turbulent flow.