Numerical simulation of liquefied fuel spills: I. Instantaneous release into a confined area

The shallow-water equations in radial symmetry are solved numerically to simulate the collapse of a cylindrical liquid column into an area surrounded by a concentric dike. The following three subcases of this problem are considered: a liquid column collapsing onto a layer of the same liquid, a liquid column collapsing onto a solid surface, and a column of lighter liquid collapsing onto a heavier liquid (i.e. liquefied natural gas (LNG) spilled onto water). The results for the three categories are compared and the differences and similarities between them are analysed.     

Numerical modeling of wave interactions with coastal structures by a constrained interpolation profile/immersed boundary method

A high‐order difference method based multiphase model is proposed to simulate nonlinear interactions between water wave and submerged coastal structures. The model is based on the Navier–Stokes equations using a constrained interpolation profile (CIP) method for the flow solver, and employs an immersed boundary method (IBM) for the treatment of wave–structure interactions. A more accurate interface capturing scheme, the volume of fluid/weighed line interface calculation (VOF/WLIC) scheme, is adopted as the interface capturing method. A series of computations are performed to verify the application of the model for simulations of fluid interaction with various structures. These problems include flow over a fixed cylinder, water entry of a circular cylinder and solitary waves passing various submerged coastal structures. Computations are compared with the available analytical, experimental and other numerical results and good agreement is obtained. The results of this study demonstrate the accuracy and applications of the proposed model to simulate the nonlinear flow phenomena and capture the complex free surface flow. Copyright © 2015 John Wiley & Sons, Ltd.     

Dynamic stabilization of Rayleigh capillary instability

The effect of high-frequency vibrations on the stability of a cylindrical liquid interface is studied. In the absence of external disturbances the interface will be unstable if the length of the liquid cylinder exceeds the length of the normal boundary section. It is shown that vibrations circularly polarized in the plane of the normal section can suppress the development of instability however great the length of the liquid cylinder. The effect of the density ratio of the liquids and the dimensions of the rigid outer shell on the stability of the system is investigated. It is shown that vibrations can stabilize the cylindrical interface only if the radius of the shell is not too great as compared with the radius of the liquid cylinder. The critical value of the radius ratio is approximately equal to 1.58 and does not depend on the density ratio of the liquids.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 3–7, November–December, 1991.     

Numerical simulation of the flow with contact lines

The paper proposes a physical model for the motion of the contact line and the gas-liquid interface. The local motion of the contact line at the solid wall is assumed and the interface between gas and liquid is traced by a level function. The numerically. The motion of the water column in a vertical pipe is computed and the results are in good agreement with experimental data.     

固体火箭燃气射流驱动液柱过程的CFD分析

固体火箭燃气射流驱动液柱过程会产生一个复杂的非稳态多相流场,为了研究液柱对固体火箭发动机工作过程中射流流场的降温效果,并揭示燃气冲击液柱的流动演化和气水之间的相互作用,利用FLUENT软件中耦合了液态水汽化方程的VOF多相流计算模型对燃气与液柱之间的耦合流动及相变过程进行了数值模拟,并与无液柱情况下射流流场的计算结果进行了对比分析。计算结果表明,当有液柱平衡体时射流流场中的压力、温度、速度波动幅度均减小,减弱了射流流场中的湍流脉动强度;液柱与燃气之间的汽化以及液柱的阻碍作用减小了射流流场的轴向发展位移,尾管后的完全发展射流流场核心区域内的压力峰值降低了0.9 MPa,温度峰值降低了503 K,速度峰值降低了291 m/s,验证了实验中液柱对燃气射流流场的降温效果。     

Experimental investigation of the influence of column scale,gas density and liquid properties on gas holdup in bubble columns

Measurements of gas holdups in bubble columns of 0.16, 0.30 and 0.33 m diameter were carried out. These columns were operated in co-current flow of gas and liquid phases and in semibatch mode. The column of 0.33 m diameter was operated at elevated pressures of up to 3.6 MPa. Nitrogen was employed as the gas phase and deionized water, aqueous solutions of ethanol and acetone and pure acetone and cumene as the liquid phase. The effects of differing liquid properties, gas density (due to elevated pressure), temperature, column diameter and superficial liquid velocity on gas holdup were studied. The gas holdup measurements were utilized by differential pressure measurements at different positions along the height of the bubble columns which allowed for the identification of axial gas holdup profiles. A decrease of gas holdup with increasing column diameter and an increase of gas holdup with increasing pressure was observed. The effect of a slightly decreasing gas holdup with increasing liquid velocity was found to exist at smaller column diameters. The use of organic solvents as the liquid phase resulted in a significant increase in gas holdup compared to deionized water. It is found that published gas holdup models are mostly unable to predict the results obtained in this study.     

反射激波作用下三维凹气柱界面演化的数值研究

激波与气柱相互作用是Richtmyer-Meshkov不稳定性研究的经典案例.单次激波与二维气柱相互作用已得到广泛关注,但是反射激波再次冲击气柱(尤其是三维气柱)的研究较少,相关演化规律和机理尚不清楚.反射激波再次冲击演化中的气柱界面会产生新的斜压涡量,影响涡量的输运和分布,从而影响界面的演化.本文采用自主开发的HOW...     

Effect of pulsations on the stability of a gas column

The effect of radial pulsations on the stability of a compressible cylindrical gas column surrounded by an ambient liquid is discussed. In the absence of pulsations, the stationary interface is susceptible to the Rayleigh capillary instability, which promotes the growth of longitudinal waves whose wave length is larger than 2 times the column radius, irrespective of the Reynolds number. A Floquet stability analysis for potential flow shows that the pulsations further destabilize the interface by extending the range of unstable wave numbers to a sequence of islands. A similar stability analysis for Stokes flow shows that the pulsations also have a destabilizing influence, though the presence of an insoluble surfactant has a competing stabilizing influence that may cause an overall reduction in the range of unstable wave numbers.     

A numerical scheme for strong blast wave driven by explosion

After the detonation of a solid high explosive, the material has extremely high pressure keeping the solid density and expands rapidly driving strong shock wave. In order to simulate this blast wave, a stable and accurate numerical scheme is required due to large density and pressure changes in time and space. The compressible fluid equations are solved by a fractional step procedure which consists of the advection phase and non‐advection phase. The former employs the Rational function CIP scheme in order to preserve monotone signals, and the latter is solved by interpolated differential operator scheme for achieving the accurate calculation. The procedure is categorized into the fractionally stepped semi‐Lagrangian. The accuracy of our scheme is confirmed by checking the one‐dimensional plane shock tube problem with 10³ times initial density and pressure jump in comparison with the analytic solution. The Sedov–Taylor blast wave problem is also examined in the two‐dimensional cylindrical coordinate in order to check the spherical symmetry and the convergence rates. Two‐ and three‐dimensional simulations for the blast waves from the explosion in the underground magazine are carried out. It is found that the numerical results show quantitatively good agreement with the experimental data. Copyright © 2006 John Wiley & Sons, Ltd.     

The EHD-driven fluid flow and deformation of a liquid jet by a transverse electric field

The aim of this study is to investigate the effect of a uniform transverse electric field on the steady-state behavior of a liquid cylinder surrounded by another liquid of infinite extent. The governing electrohydrodynamic equations are solved for Newtonian and immiscible fluids in the framework of leaky-dielectric theory and in the limit of small electric field and fluid inertia. A detailed analysis of the electrical and hydrodynamic stresses acting on the interface separating the two fluids is presented, and an expression is found for the interface deformation for small distortions from a circular shape. The electrical stresses acting on the interface of two leaky-dielectric liquids are compared with those acting on an interface separating a perfect dielectric or infinitely conducting core fluid cylinder from a surrounding perfect dielectric fluid. A comparison is made between the results of this study and those of a similar study for fluids with permeable interfaces and the classical results for liquid drops.     

Self-similar solution to the problem of the expansion of a cylindrical column of conducting gas in a longitudinal magnetic field

The nonstationary radial motion of a long cylindrical column of conducting gas in a time-varying longitudinal magnetic field is considered. Exact solutions are found by the method of separating the variables for the system of equations of magnetohydrodynamics on the assumption that the statistical pressure of the plasma at the boundary of the column is proportional to the external magnetic pressure. Some numerical computations are performed and the energetic characteristics of the interaction process are calculated. The ratio of the useful work done by the gas over an infinite time interval to the initial energy of the column is given as a function of the magnetic Reynolds number. We note that a similar method was applied in [1], where not only was the average temperature taken over the cross section, but the inertia of the medium was also neglected. When the inertia is taken into account, we have the additional requirement that the statistical pressure be proportional to the magnetic pressure at the boundary of the column.A physically similar model may be interpreted, for example, as the expansion of a compressible conducting gas column in a nonconducting incompressible fluid situated in a permeable cylinder of some radius R infinite along the axis of symmetry. The requirement that the statistical pressure be proportional to the magnetic pressure reduces to the condition that the external pressure on the boundary of the permeable cylinder of radius R should vary according to a specific law, which may easily be determined.     

Two-Phase Flow and the Capillarity Phenomenon in Porous Solids – A Continuum Thermomechanical Approach

Two-phase flow and capillarity phenomenon in porous solids, well known in physics and engineering, are treated from a rigorous continuum thermomechanical point of view for the first time. A ternary model, consisting of a porous solid phase, a liquid phase, and a gas phase, is investigated within the framework of thermodynamics. The main result of the evaluation of the entropy principle turns out to be that the interaction forces between the solid, gas, and liquid phases are dependent on the free Helmholtz energy functions of the corresponding phases and on the gradient of the liquid density. The classical result for the driving volume force for raising a water column in a narrow tube against the force of gravity is contained in the general investigation.     

An efficient shock-capturing algorithm for compressible flows in a duct of variable cross-section

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a gas in a single spatial co-ordinate. This includes flow in a duct of variable cross-section as well as flow with slab, cylindrical or spherical symmetry and can prove useful when testing codes for the two-dimensional equations governing compressible flow of a gas. The resulting scheme requires an average of the flow variables across the interface between cells and for computational efficiency this average is chosen to be the arithmetic mean, which is in contrast to the usual ‘square root’ averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and a comparison is made in the cylindrical case with results from a two-dimensional problem with no sources.     

Hydrodynamics and energy consumption studies in a three-phase liquid circulating three-phase fluid bed contactor

The hydrodynamics and energy consumption have been studied in a cold flow, bubbling and turbulent, pressurized gas–liquid–solid three-phase fluidized bed (0.15 m ID × 1 m height) with concurrent gas–liquid up flow is proposed with the intention of increasing the gas hold up. The hydrodynamic behaviour is described and characterised by some specific gas and liquid velocities. Particles are easily fluidized and can be uniformly distributed over the whole height of the column. The effect of parameters like liquid flow rate, gas flow rate, particle loading, particle size, and solid density on gas hold up and effect of gas flow rate, solid density and particle size on solid hold up, energy consumption and minimum fluidization velocity has been studied. At the elevated pressures a superior method for better prediction of minimum fluidization velocity and terminal settling velocities has been adopted. The results have been interpreted with Bernoulli’s theorem and Richardson–Zaki equation. Based on the assumption of the gas and liquid as a pretend fluid, a simplification has been made to predict the particle terminal settling velocities. The Richardson–Zaki parameter n′ was compared with Renzo’s results. A correlation has been proposed with the experimental results for the three-phase fluidization.     

Flow pattern transition in liquid-liquid flows with a transverse cylinder

The effect of a cylindrical bluff body on the interface characteristics of stratified two-phase, oil-water, pipe flows is experimentally investigated with high speed Particle Image Velocimetry (PIV). The motivation was to study the feasibility of flow pattern map actuation by using a transverse cylinder immersed in water in the stratified pattern, and particularly the transition from separated to dispersed flows. The cylinder has a diameter of 5 mm and is located at 6.75 mm from the bottom of the pipe in a 37 mm ID acrylic test section. Velocity profiles were obtained in the middle plane of the pipe. For reference, single phase flows were also investigated for Reynolds numbers from 1550 to 3488. It was found that the flow behind the cylinder was similar to the two dimensional cases, while the presence of the lower pipe wall diverted the vorticity layers towards the top. In two-phase flows, the Froude number (from 1.4 to 1.8) and the depth of the cylinder submergence below the interface affected the generation of waves. For high Froude numbers and low depths of submergence the counter rotating von Karman vortices generated by the cylinder interacted with the interface. In this case, the vorticity clusters from the top of the cylinder were seen to attach at the wave crests. At high depths of submergence, a jet like flow appeared between the top of the cylinder and the interface. High speed imaging revealed that the presence of the cylinder reduced to lower mixture velocities the transition from separated to dual continuous flows where drops of one phase appear into the other.     

Water Evaporation Versus Condensation in a Hygroscopic Soil

The liquid/vapour phase change of water in soil is involved in many environmental geotechnical processes. In the case of hygroscopic soils, the liquid water is strongly adsorbed on the solid phase and this particular thermodynamic state can highly influence the phase change kinetics. Based on the linear Thermodynamic of Irreversible Processes ideas, the non-equilibrium phase change rate is written as a linear function of the water chemical potential difference between the liquid and vapour state. In this relation, the system is characterized by a phenomenological coefficient that depends on the state variables. Using an original experimental set-up able to analyze the response of a porous medium subjected to non-equilibrium conditions, the phase change coefficient is determined in various configurations. This paper focuses on the influence of the gas phase pressure and underlines that a low gas pressure decreases the phase change kinetics. Then, evaporation and condensation processes are compared showing an asymmetric behaviour. These experimental results are interpreted from a microscopic point of view by relying on recent works dealing with molecular dynamics numerical simulation of the liquid/gas interface.     

A comparative study of direct‐forcing immersed boundary‐lattice Boltzmann methods for stationary complex boundaries

In this study, we assess several interface schemes for stationary complex boundary flows under the direct‐forcing immersed boundary‐lattice Boltzmann methods (IB‐LBM) based on a split‐forcing lattice Boltzmann equation (LBE). Our strategy is to couple various interface schemes, which were adopted in the previous direct‐forcing immersed boundary methods (IBM), with the split‐forcing LBE, which enables us to directly use the direct‐forcing concept in the lattice Boltzmann calculation algorithm with a second‐order accuracy without involving the Navier–Stokes equation. In this study, we investigate not only common diffuse interface schemes but also a sharp interface scheme. For the diffuse interface scheme, we consider explicit and implicit interface schemes. In the calculation of velocity interpolation and force distribution, we use the 2‐ and 4‐point discrete delta functions, which give the second‐order approximation. For the sharp interface scheme, we deal with the exterior sharp interface scheme, where we impose the force density on exterior (solid) nodes nearest to the boundary. All tested schemes show a second‐order overall accuracy when the simulation results of the Taylor–Green decaying vortex are compared with the analytical solutions. It is also confirmed that for stationary complex boundary flows, the sharper the interface scheme, the more accurate the results are. In the simulation of flows past a circular cylinder, the results from each interface scheme are comparable to those from other corresponding numerical schemes. Copyright © 2010 John Wiley & Sons, Ltd.     

纳米阵列中气体驱替液体的流动特征

纳米尺度下气体驱动液体流动特征在纳流控芯片及页岩气开发中具有广泛的应用前景. 利用管径规格为292.8 nm,206.2 nm,89.2 nm,67.0 nm,26.1 nm的氧化铝膜为纳米阵列,进行气驱水实验和单相气体流动实验,分析纳米尺度下气驱水流动特征. 实验表明,纳米阵列中气驱水时气体流量随驱动压力变化经历三个阶段:第一阶段流量缓慢增大,且比单相气体流量降低约一个数量级;第二阶段纳米阵列中的水被大量驱替出,流量迅速增大;第三阶段纳米阵列中的水全部被驱替出,流动特征与单相气体流动保持一致. 分析表明,气驱水第一阶段存在气液界面毛细管力的“钉扎”作用及固液界面相互作用力的影响,是产生非线性流动的主要原因;而一旦“钉扎”作用破坏,气体进入管道推动界面运动,气柱与液柱之间的毛细曲面曲率变化,毛细管力减小,气体流量急剧增大,其中毛细管力随驱替压力增大急剧变化,是造成第二阶段气体流量突变的主要原因.      

Equations of the mechanics of porous media saturated with liquid and gas

The approach proposed by Podil'chuk [1] is used to derive a system of equations of motion for saturated porous media, allowance being made for the mutual influence of the solid, liquid, and gas phases. The permeabilities of the anisotropic porous medium are assumed to depend on the direction. It is shown that when there are no gas phases and the liquid is incompressible the system of equations reduces to the general equations of the theory of elasticity of an anisotropic body with fictitious stress components. For a porous medium saturated with liquid, the relationships between the permeabilities and the anisotropy constants are obtained. The motion of liquid in an elastic porous medium in the form of an orthotropic cylindrical region with a cavity in the form of a circular cylinder is considered as an example.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 82–87, July–August, 1981.     

Unstable waves at the interface of a diamagnetic liquid column contained in a magnetic fluid

A magnetically held pipeless fluid transporting system is a potentially useful noncontact technique in science and technology. Keeping a stable interface between the diamagnetic liquid column and the surrounding magnetic fluid is essential. The present study gives an experimental study on the unstable interfacial wave leading to discontinuity of a liquid column maintained in a magnetic fluid by a non-uniform magnetic field. The experiments were carried out by setting two magnetic north poles in opposite position to produce the magnetic field, using water to form the diamagnetic liquid column, and a diluted kerosene-base magnetic fluid as the surrounding fluid. Unstable interfacial waves were measured. It was found that the magnetic pressure acting on the interface can stabilise the interface. The experimental results agree well with the theoretical stability condition derived from a two-fluids model.