Modelling of a single drop impact onto liquid film using particle method

The process of single liquid drop impact on thin liquid surface is numerically simulated with moving particle semi‐implicit method. The mathematical model involves gravity, viscosity and surface tension. The model is validated by the simulation of the experimental cases. It is found that the dynamic processes after impact are sensitive to the liquid pool depth and the initial drop velocity. In the cases that the initial drop velocity is low, the drop will be merged with the liquid pool and no big splash is seen. If the initial drop velocity is high enough, the dynamic process depends on the liquid depth. If the liquid film is very thin, a bowl‐shaped thin crown is formed immediately after the impact. The total crown subsequently expands outward and breaks into many tiny droplets. When the thickness of the liquid film increases, the direction of the liquid crown becomes normal to the surface and the crown propagates outward. It is also found that the radius of the crown is described by a square function of time: r_C = [c(t ? t₀)]^0.5. When the liquid film is thick enough, a crown and a deep cavity inside it are formed shortly after the impact. The bottom of the cavity is initially oblate and then the base grows downward to form a sharp corner and subsequently the corner moves downward. Copyright © 2004 John Wiley & Sons, Ltd.     

Hydrodynamics of contact of a falling drop with a liquid free surface

A fine structure of the flows developing during primary contact of freely falling drops with a deep quiescent fluid is studied using the macrophotography and high-speed video filming methods. Water drops falling in water, alcohol, and oil, as well as drops of oil, petroleum, and aqueous solutions of salt or alcohol falling in water are investigated. The work is focused on the visualization of the finespray scattering from the primary contact area. The collisions of small droplets with the surface of the submerging drop are first recorded. The direction of the spray and streamer scattering is determined by the surface tension coefficients of the coalescing liquids. The conditions under which the spray droplets collide with the drop surface are determined.     

An Analysis of the Movement of Wetting and Nonwetting Fluids in Homogeneous Porous Media

The movement of wetting and nonwetting fluid flow in columns packed with glass beads is used to understand the more complicated flows in homogeneous porous media. The motion of two immiscible liquids (oil and water) is observed with different surfactants. Through dimensional analyses, fluid velocity is well correlated with interfacial tension and less dependent on particle size. In water–oil (W/O) experiments, finger pattern flows are observed if water is the displacing fluid that flows in an oil-filled porous media, whereas oil ganglia tend to form if oil is the displacing fluid in the water-wetted porous media. The results are well described by a simple model based on an earlier theory of flow in a tube.     

Pore-Scale Monitoring of Wettability Alteration by Silica Nanoparticles During Polymer Flooding to Heavy Oil in a Five-Spot Glass Micromodel

It is well known that the oil recovery is affected by wettability of porous medium; however, the role of nanoparticles on wettability alteration of medium surfaces has remained a topic of debate in the literature. Furthermore, there is a little information of the way dispersed silica nanoparticles affect the oil recovery efficiency during polymer flooding, especially, when heavy oil is used. In this study, a series of injection experiments were performed in a five-spot glass micromodel after saturation with the heavy oil. Polyacrylamide solution and dispersed silica nanoparticles in polyacrylamide (DSNP) solution were used as injected fluids. The oil recovery as well as fluid distribution in the pores and throats was measured with analysis of continuously provided pictures during the experiments. Sessile drop method was used for measuring the contact angles of the glass surface at different states of wettability after coating by heavy oil, distilled water, dispersed silica nanoparticles in water (DSNW), polyacrylamide solution, and DSNP solution. The results showed that the silica nanoparticles caused enhanced oil recovery during polymer flooding by a factor of 10%. The distribution of DSNP solution during flooding tests in pores and throats showed strong water-wetting of the medium after flooding with this solution. The results of sessile drop experiments showed that coating with heavy oil, could make an oil-wet surface. Coating with distilled water and polymer solution could partially alter the wettability of surface to water-wet and coating with DSNW and DSNP could make a strongly water-wet surface.     

Impact dynamics of drops on thin films of viscoelastic wormlike micelle solutions

The impact dynamics of water drops on thin films of viscoelastic wormlike micelle solutions is experimentally studied using a high-speed digital video camera at frame rates up to 4000 frame/s. The composition and thickness of the thin film is modified to investigate the effect of fluid rheology on the evolution of crown growth, the formation of satellite droplets and the formation of the Worthington jet. The experiments are performed using a series of wormlike micelle solutions composed of a surfactant, cetyltrimethylammonium bromide (CTAB), and a salt, sodium salicylate (NaSal), in deionized water. The linear viscoelastic shear rheology of the wormlike micelle solutions is well described by a Maxwell model with a single relaxation time while the steady shear rheology is found to shear thin quite heavily. In transient homogeneous uniaxial extension, the wormlike micelle solutions demonstrate significant strain hardening. The size and velocity of the impacting drop is varied to study the relative importance of Weber, Ohnesorge, and Deborah numbers on the impact dynamics. The addition of elasticity to the thin film fluid is found to suppress the crown growth and the formation of satellite drops with the largest effects observed at small film thicknesses. A new form of the splashing threshold is postulated which accounts for the effects of viscoelasticity and collapses the satellite droplet data onto a single master curve dependent only on dimensionless film thickness and the underlying surface roughness. Additionally, a plateau is observed in the growth of the maximum height of the Worthington jet height with increasing impact velocity. It is postulated that the complex behavior of the Worthington jet growth is the result of a dissipative mechanism stemming from the scission of wormlike micelles.     

回转体高速倾斜入水的流场特性及结构响应

回转体高速入水过程涉及液体和固体的耦合作用，是一个复杂的非线性、非定常过程。为研究回转体高速入水的结构动响应及流场演变规律，本文中基于STAR-CCM+和ABAQUS平台，建立了回转体高速入水的双向流固耦合数值模型，开展了不同入水速度的回转体高速倾斜入水流固耦合数值计算。结果表明:数值计算的入水速度、位移曲线和空泡形态与实验结果良好吻合，验证了流固耦合方法的有效性；回转体倾斜高速入水的载荷先集中在触水部分边缘处，后集中于回转体底部中心处；流固耦合方法的入水冲击载荷峰值小于刚体的，弹性回转体的载荷曲线产生明显波动；撞水阶段，回转体空泡呈现不对称形态，随着入水加深，空泡不对称性变弱；入水速度60 m/s下，空泡发生表面闭合，回转体入水初速度越快，空泡表面闭合越晚；冲击载荷与入水速度有关，入水速度越大，峰值出现越早，震荡越明显，速度超过100 m/s时，回转体产生塑性形变。     

Viscous oil–water flows in a microchannel initially saturated with oil: Flow patterns and pressure drop characteristics

Immiscible viscous liquid–liquid two-phase flow patterns and pressure drop characteristics in a circular microchannel have been investigated. Water and silicone oil with a dynamic viscosity of 863 mPa s were injected into a fused silica microchannel with an inner diameter of 250 μm. As the microchannel was initially filled with the silicone oil, an oil film was found to always form and remain on the microchannel wall. Different flow patterns were observed and classified over a wide range of water and oil flow rates. A flow pattern map is presented in terms of Re, Ca, and We numbers. Two-phase pressure drop data have also been collected and analyzed to develop a simple correlation for slug, annular and annular-droplet flow patterns in terms of superficial water and oil velocities.     

开放空腔壳体倾斜入水运动特性试验研究

基于高速摄像试验方法,研究了开放空腔壳体的倾斜入水运动特性,重点分析了开放空腔结构引起的空泡流动特征和壳体运动规律. 通过试验数据分析了开放空腔内气体运动将引起独特的空泡流动和阶段性的运动规律,探讨了初始入水速度、入水姿态对入水弹道和空泡形态等运动特征的影响. 结果表明:开放空腔壳体入水空泡出现阶段波动演化现象,并先后经历两次闭合;入水空泡演化改变流体动力分布,直接影响壳体运动方式,进而改变水下弹道特征;空腔内部形成相对独立流场环境和开放端周期性流动,在重力作用下液体对空腔内下侧壁面作用力较大,加剧壳体偏转,从而改变入水运动过程的稳定性;随着入水速度的增大,空泡波动特征逐渐明显,闭合时间延迟,非对称深闭合引起的横向位移减小,但偏转角度与入水速度无关;随着初始姿态倾角减小,空泡波动程度减弱、闭合时间延迟,偏转角速度增大,闭合引起的横向位移增大.      

Flow induced on a salt waterbody due to the impingement of a freshwater drop or a water source

The particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF) techniques are used to study the flow induced on the surface of a body of saltwater when a drop impinges on its surface or when a source is present on the surface. The measurements show that the impingement of a fresh water drop causes a strong axisymmetric solutocapillary flow about the vertical line passing through the center of impact. The fluid directly below the center of impact rises upward, and near the surface it moves away from the center of impact. The flow, which develops within a fraction of second after the impact, persists for several seconds. In comparison, when a freshwater drop falls on a body of freshwater, the flow induced on the surface is much weaker and persists for a relatively shorter duration of time and the volume of water circulated is two orders of magnitude smaller. Similarly, when a fresh water source is present on a body of saltwater there is a solutocapillary flow which on the surface is away from the source and below the surface is towards the source.     

A risk based performance evaluation of plate-and-frame heat exchangers

Plate-and-frame heat exchangers (PHEs) operating in process industries are fouled to a greater or lesser extent depending on surface temperature, surface condition, material of construction, fluid velocity, flow geometry and fluid composition. This fouling phenomenon is time-dependent and will result in a decrease in the overall heat transfer coefficient and increase in the pressure drop of the PHE. Once the overall heat transfer coefficient decreases to a minimum acceptable level, cleaning of the equipment becomes necessary to restore the performance. In this paper, we present a simple probabilistic approach to characterize various fouling models that are commonly encountered in many industries. These random fouling growth models are then used to investigate the impact on risk based thermal effectiveness, overall heat transfer coefficient and the hot- and cold-fluid outlet temperatures of a PHE. All the results are presented in a generalized form in order to demonstrate the generality of the risk-based procedure discussed in this paper.     

Characterization of the critical transition from annular to wavy-stratified flow for oil–water mixtures in horizontal pipes

The transition from annular to wavy-stratified oil–water adiabatic flow within horizontal pipes is experimentally analyzed, and a semiempirical model is proposed. The transition is referred to as critical because it occurs suddenly, giving rise to a sharp and strong increase in the pressure drop due to the contact of the high-viscosity oil with the pipe wall. This could lead to a dangerous accident in pipelines. Experimental runs were performed on eight test sections of both Plexiglas^® and Pyrex^® pipes with internal diameters ranging from 21.5 to 50 mm, using tap water and oil with viscosity about 880 times higher than that of water. On the basis of pressure drop measurement and flow pattern visualization, the transition boundary between annular and wavy-stratified flow was analytically determined and compared with flow pattern maps.     

Crown behavior and bubble entrainment during a drop impact on a liquid film

Physical and mathematical models are established to simulate a single liquid drop impinging onto a liquid film using the coupled level set and volume of fluid method. The crown liquid sheet after impact is obtained, which coincides well with the experimental results in literatures. Influence of Weber number, Reynolds number and the dimensionless film thickness on the crown diameter and height is discussed quantitatively. Results indicate that the crown diameter is independent of the two non-dimensional numbers, while it can be increased by reducing the dimensionless film thickness. The crown height increases with the increasing of Weber number, but Reynolds number has small effect on it. Mechanism of the jet formation process is revealed by analyzing pressure distribution and velocity field in the liquid. It is found that both pressure difference in the neck region and velocity discontinuity can greatly affect the jet formation. Besides, the bubble entrainment phenomenon during a liquid drop impact on a liquid film is successfully captured with this numerical method. It is found that the increase in both impact Weber number and the drop diameter contributes to the emerging of bubble rings.     

Collision between immiscible drops with large surface tension difference: diesel oil and water

The collision outcomes of immiscible drops with large surface tension difference, namely, a water drop and a diesel oil drop, were observed experimentally. In a near head-on collision between immiscible drops with large surface tension difference, an “overlaying” action for the drop of the smaller surface tension, i.e., the diesel oil drop, to go around the surface of the drop of the larger surface tension, i.e., the water drop, occurs during the collision. This overlaying action reduces the reflex energy for head-on collisions, making reflex separation more difficult to occur. At the same time, due to the immiscibility, the liquid bridge during stretching separation becomes narrower, which makes stretching separation easier to happen. No coalescence could be observed for a collision of Weber number greater than 60. In addition, compound drops are produced frequently.     

Symmetry Properties of the Elastic Energy of a Woven Fabric with Bending and Twisting Resistance

A woven fabric can be described as a surface made of two families of fibers: in this work we study how the geometry of the weave pattern affects the symmetry properties of the elastic energy of the surface. Four basic symmetry classes of weave patterns are possible, depending on the angle between the fibers and their material properties. The properties of the pattern determine the material symmetry group of the network, under which the elastic energy is invariant. We derive representations for the energy of a woven fabric that are invariant under the symmetry group of the network, and discuss the relation of these invariants with the curvature and twist of the fibers.      

Scattered fracture of porous materials with brittle skeleton

A model of damage accumulation in a porous medium with a brittle skeleton saturated with a compressible fluid is formulated in the isothermal approximation. The model takes account of the skeleton elastic energy transformation into the surface energy of microcracks. In the case of arbitrary deformations of an anisotropic material, constitutive equations are obtained in a general form that is necessary and sufficient for the objectivity and thermodynamic consistency principles to be satisfied. We also formulate the kinetics equation ensuring that the scattered fracture dissipation is nonnegative for any loading history. For small deviations from the initial state, we propose an elastic potential which permits describing the principal characteristics of the behavior of a saturated porous medium with a brittle skeleton. We study the acoustic properties of the material under study and find their relationship with the strength criterion depending on the accumulated damage and the material current deformation. We consider the problem of scattered fracture of a saturated porous material in a neighborhood of a spherical cavity. We show that the cavity failure occurs if the Hadamard condition is violated.     

Modeling of a Dynamic Magneto-Fluid Seal in the Presence of a Pressure Drop

The ferrohydrodynamic problem, with a disconnected free surface, of the stability of an annular magneto-fluid seal under the action of centrifugal forces and a pressure drop is solved numerically. The effect of the shaft speed on the critical pressure drop is investigated in relation to the volume of the magneto-fluid plug, the magnetic field strength, the magnetic properties of the fluid, the gap width, and the shaft radius. The flow pattern and the thermal power released in the sealing layer are studied.     

Investigation of the characteristics of particulate flows through fibrous filters using the lattice Boltzmann method

A fibrous filter is one of the most common systems used to separate suspended particles from air.Two important factors(i.e.,the pressure drop and capture efficiency) are usually used to evaluate the performance of this type of filter.This study considers three two-dimensional arrangements of fibers(parallel,staggered,and random) to geometrically model fibrous media.The lattice Boltzmann method is employed to numerically simulate fluid flow through the filter.The Lagrangian form of the equation of motion of a particle is numerically solved to track the path of each particle in the flow field,where a one-way interaction between the fluid and particles is considered.The effects of pertinent parameters such as the fiber arrangement,solid volume fraction,particle-to-fiber diameter ratio,particle-to-fluid density ratio,Reynolds number,Stokes number,and size of the fibrous medium on the pressure drop and capture efficiency are studied.The obtained results are compared with existing empirical and theoretical findings and discussed.     

Investigation on heavy crude-water two phase flow and related flow characteristics

In this paper, heavy crude oil–water flows are studied in a horizontal stainless steel test section with 25.4 mm ID and overall length of 50 m. Crude oil (viscosity = 628.1 mPa s, interfacial tension with water = 10.33 mN/m at 60 °C) and water, collected from an oilfield, were used as test fluids. Visual observations, local sampling and pressure drop measurements were used to identify the flow patterns and their transitions. It was found that in all conditions studied there was a water-in-oil emulsion present. At low mixture velocities and water fractions this occupied the whole pipe cross section. As the velocity or the volume fraction increased water appeared to segregate. At high water fractions and mixture velocities annular flow appeared with the water-in-oil emulsion in the core surrounded by a water layer. The results were compared with those from a model oil with the same viscosity. At low water fractions there was a similarity between the patterns observed with the two oil systems characterized by water segregation from an oil continuous dispersion with increasing water fraction or mixture velocity. However, at high water fractions an oil-in-water dispersion formed with the model oil that was not seen with the crude oil. Pressure drop was generally higher for the crude oil system compared to the model one, while in both cases it decreased when water started to segregate and form layers in contact with the pipe wall. The differences between the two oil systems are attributed to the natural surfactants present in the heavy crude oil (such as asphaltenes and resins), which tend to accumulate on the water/oil interface, retard film drainage and maintain the stability of water drops in oil.     

Forms of partial breakdown of an oil body in a compound vortex

Over a wide range of parameters, the pattern of a two-layer (sunflower oil-water) fluid flow in a compound vortex is visualized. The vortex is formed by a disk mounted on the bottom of a cylindrical container and rotating at a constant angular velocity. The shapes of the interface of two media (water-oil, oil-air) are determined for different flow regimes, and the conditions and types of their breakdown are studied. Under supercritical conditions, the individual oil droplets, detached from the oil-water interface, form a direct water-oil emulsion near the lower edge of the body. With increase in the angular velocity of the disk, an invert emulsion is formed, i.e. cells containing water with oil shells. The dimensions of typical structural components of the flow are determined.     

Equilibrium shapes of a drop pendant in an electrostatic field

The stationary shapes of a conducting fluid drop in the gap between the plates of a plane capacitor are studied. The drop is held on the upper plate by the surface tension forces. The self-consistent problem of the determination of the drop shape and the charge distribution over its surface is solved. Estimates are obtained for the maximum volume of the stationary drop at the given fluid parameters and electric field strength.