Impact of a strip of compressible liquid on a barrier

The exact solution of the plane problem of the impact of a finite liquid strip on a rigid barrier is obtained in the linearized formulation. The velocity components, the pressure and other elements of the flow are determined by means of a velocity potential that satisfies a two-dimensional wave equation. The final expressions for them are given in terms of elementary functions that clearly reflect the wave nature of the motion. The exact solution has been thoroughly analyzed in numerous particular cases. It is shown directly that in the limit the solution of the wave problem tends to the solution of the analogous problem of the impact of an incompressible strip obtained in [1]. A logarithmic singularity of the velocity parallel to the barrier in the corner of the strip is identified. A one-dimensional model of the motion, which describes the behavior of the compressible liquid in a thin layer on impact and makes it possible to obtain a simple solution averaging the exact wave solution, is proposed. Inefficient series solutions are refined and certain numerical data on the impact characteristics for a semi-infinite compressible liquid strip, previously considered in [2–4] in connection with the study of the earthquake resistance of a dam retaining water in a semi-infinite basin, are improved. The solution obtained can be used to estimate the forces involved in the collision of solids and liquids. It would appear to be useful for developing correct and reliable numerical methods of solving the nonlinear problems of fluid impact on solids often examined in the literature [5].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 138–145, November–December, 1990.The results were obtained by the author under the scientific supervision of B. M. Malyshev (deceased).     

Investigation on the generation process of impact-driven high-speed liquid jets using a CFD technique

High-speed liquid jets have been applied to many fields of engineering, science and medicine. It is therefore of benefit to all these areas to investigate their characteristics by modern and inexpensive methods using a computational fluid dynamics (CFD) technique. Previously, high-speed liquid jets have been studied experimentally using a momentum exchange method, called the “impact driven method (IDM)”, by which the impact of a high-velocity projectile on the liquid package contained in the nozzle cavity produced the jet. The shock pulse reflections in the cavity caused by the impact then drove a multiple pulsed jet from the nozzle exit. In this study, a two-fluid simulation consisting of liquid and air can be successfully calculated by using a two-phase flow mixture model and a moving mesh for the projectile motion. The CFD results show good agreement to the results of previous experimental studies, both quantitatively and qualitatively. For the first time, the wave propagation within the liquid in the nozzle has been captured and analyzed, thereby demonstrating the dynamic characteristics of multiple pulsed high-speed liquid jets initiated by the IDM. This provides a breakthrough in the simulation of the supersonic injection of a liquid into air by using a well-known and user-friendly CFD software. It is useful fundamental knowledge for future studies of high-speed injection with applications in all its related fields.     

Longitudinal waves at a micropolar fluid/solid interface

The possibility of plane wave propagation in a micropolar fluid of infinite extent has been explored. The reflection and transmission of longitudinal elastic wave at a plane interface between a homogeneous micropolar fluid half-space and a micropolar solid half-space has also been investigated. It is found that there can exist four plane waves propagating with distinct phase speeds in an infinite micropolar fluid. All the four waves are found to be dispersive and attenuated. The reflection and transmission coefficients are found to be the functions of the angle of incidence, the elastic properties of the half-spaces and the frequency of the incident wave. The expressions of energy ratios have also been obtained in explicit form. Frequency equation for the Stoneley wave at micropolar solid/fluid interface has also been derived in the form of sixth-order determinantal expression, which is found in full agreement with the corresponding result of inviscid liquid/elastic solid interface. Numerical computations have been performed for a specific model. The dispersion curves and attenuation of the existed waves in micropolar fluid have been computed and depicted graphically. The variations of various amplitudes and energy ratios are also shown against the angle of incidence. Results of some earlier workers have been deduced from the present formulation.     

Impact on the surface of a compressible liquid

The evolution of a weak, nearly plane shock wave produced by the impact on the plane boundary of a compressible liquid is considered. At the initial moment the liquid is at rest and occupies the lower half-plane. Then the points of its boundary get instantly velocities directed into the liquid domain. This leads to the formation of a shock wave the intensity of which is non-uniform due to a non-uniform distribution of the impact velocities. Initially the shock wave is plane but then it bends due to the non-linear effects and can later be focused. To analyze the liquid flow, the method of matched asymptotic expansions is used. For finite times the flow and the evolution of the shock wave are described within the framework of the acoustic approximation. For large times the flow becomes non-linear, and the form of the shock front depends essentially on the characteristics of the liquid flow behind it. If the non-uniformity of the impact velocity distribution is slight then the focusing of the shock wave is shown not to occur. The influence of viscosity of the liquid on the structure of its motion is discussed.     

Propagation of Normal Axisymmetric Waves in Compliant Elastic Cylinders Filled with and Immersed in a Fluid

The paper studies the relationship between the physical characteristics of a cylinder and the properties of normal axisymmetric waves in elastic–liquid waveguides. The cylinder is made of a compliant material in which the velocity of shear waves is less than the sonic velocity in a perfect compressible liquid. The complete system of dynamic elasticity equations and the wave equation are used to describe the wave fields in the elastic cylinder and fluid, respectively. This approach allows obtaining the dispersion characteristics of coupled normal waves in compound waveguides over wide ranges of frequencies and wavelengths. The curves of real, imaginary, and complex wave numbers versus frequency are plotted for specific pairs of waveguide materials. Computations are carried out for a thick-walled cylinder filled with a fluid and immersed in either vacuum or a fluid. It is found out that compliant and rigid materials of the cylinder affect differently the wave interaction process in elastic–liquid waveguides     

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.     

Dynamics of a rigid cylinder near a plane boundary in the radiation field of an acoustic wave

An approach is described for investigation of the interaction between a rigid body and a viscous fluid boundary under acoustic wave propagation. The influence of the liquid on the rigid body is determined as a mean force, which is a constant in the time component of the hydrodynamic force. This enables the use of a previously developed technique for calculation of pressure in a compressible viscous liquid. The technique takes into account the second-order terms with respect to the wave field parameters and is based on investigation of a system of initially nonlinear hydromechanics equations that can be simplified with respect to the wave motion parameters of the liquid. It has proven possible to retain the second-order terms for determination of stresses in the liquid without having to solve the system of nonlinear equations. The stresses can be expressed in terms of parameters found in the solution of the linearized equations of the compressible viscous liquid. In this way, the solution of linearized equations is expressed in terms of a scalar and vector potentials. The problem statement is derived for a rigid cylinder located near a rigid flat wall under the effects of a wave propagating perpendicular to the wall. The solution for this particular example is obtained.     

基于SPH方法的湿润性固壁边界条件模拟研究

固体壁面由于表面特殊结构和材料属性,时常表现出对交界面上水体的吸附作用,而这一特征对微小水体作用尤为明显。本文提出了一种湿润性固壁边界条件的计算方法,即假设壁面粒子的亲水性以及毛细吸附作用统一表现为对支持域内流体粒子的吸附力。基于光滑粒子流体动力学（SPH）方法,模拟了静态液滴在不同湿润性壁面上的变形至稳定过程。模拟了液滴撞击疏水壁面的过程,将液滴的运动过程分为碰撞、铺展、回缩和回弹四个阶段,分析各阶段壁面受力分布情况。研究表明：根据模拟液滴静态接触角的变化特点,本文湿润性固壁边界条件可以较好的反映出壁面湿润性;液滴撞击输水表面的模拟数据与试验结果趋势上吻合良好;壁面压力波伴随着液滴的铺展和回缩传播并衰减;只有在回弹后期液滴即将脱离壁面时壁面拉力起主导作用,其余各时刻壁面均以压力为主。     

基于LS-DYNA的热成型钢断裂失效预测研究

固体壁面由于表面特殊结构和材料属性,时常表现出对交界面上水体的吸附作用,而这一特征对微小水体作用尤为明显。本文提出了一种湿润性固壁边界条件的计算方法,即假设壁面粒子的亲水性以及毛细吸附作用统一表现为对支持域内流体粒子的吸附力。基于光滑粒子流体动力学（SPH）方法,模拟了静态液滴在不同湿润性壁面上的变形至稳定过程。模拟了液滴撞击疏水壁面的过程,将液滴的运动过程分为碰撞、铺展、回缩和回弹四个阶段,分析各阶段壁面受力分布情况。研究表明：根据模拟液滴静态接触角的变化特点,本文湿润性固壁边界条件可以较好的反映出壁面湿润性;液滴撞击输水表面的模拟数据与试验结果趋势上吻合良好;壁面压力波伴随着液滴的铺展和回缩传播并衰减;只有在回弹后期液滴即将脱离壁面时壁面拉力起主导作用,其余各时刻壁面均以压力为主。     

The evolution of small rotational perturbations to a uniform rotation of a viscoelastic fluid contained between parallel plates

The transfer mechanism of momentum from the boundaries into the interior of a viscoelastic fluid is examined in the case when the boundary consists of two parallel plates.Special emphasis is given to the case of a Maxwellian fluid. It is seen that an increase in the azimuthal angular velocity at the boundaries progresses into the interior of the liquid in the form of a velocity discontinuity wave. Simultaneously, an axial circulatory flow pumps liquid from the low rotation region near the middle plane into the high rotation region near the boundaries. Detailed calculations are given for the intial stage of evolution.     

Impact of a surfacewave on an elastic hull

In order to model a ship hull’s response to the impact of surface waves, the two-dimensional problem of wave impact on an elastic beam whose ends are connected by springs with a rigid structure uniformly submerged in a fluid is considered. The fluid is assumed to be ideal and incompressible and its flow symmetric; the lateral bending of the beam is described by the Euler equation. The fluid flow and the size of the wetted region are determined simultaneously with the calculation of the the beam deflection within the framework of the Wagner approach which takes into account the reshaping of the free surface of the fluid on interacting with a body. The stresses and strains arising in the beam and at its ends during impact are found. The numerical algorithm developed makes it possible to analyze the elastic effects in fluid impacts on thin-walled structures of finite length. Moreover, as the stiffness of the connecting springs tends to zero, the solution of this problem describes the impact of an elastic beam with free ends on a weakly curved fluid surface.     

An anti‐diffusive volume of fluid method for interfacial fluid flows

A volume of fluid (VOF) method is developed combining a first‐order limited downwind scheme with higher order accurate schemes. The method is characterized by retaining a sharp fluid interface and a reduction in numerical diffusion near the interface, but avoids complicated geometrical reconstruction as occurs in most volume tracing algorithms. To demonstrate the accuracy and robustness of the method, a selection of numerical experiments are presented involving a pure advection problem, a water wave impact caused by a dam breaking and liquid sloshing in a partially filled tank. Copyright © 2011 John Wiley & Sons, Ltd.     

Waves on the interface between a nematic liquid crystal and an ideal isotropic fluid

The problem of surface wave propagation over the interface between a nematic liquid crystal and an ideal isotropic fluid is considered. For the nematic liquid crystal the Frank-Oseen model with an isotropic viscous stress tensor is used. Anisotropic surface tension is described by the Rapini model. In this formulation, for the problem of harmonic small-amplitude surface wave propagation, in the case of infinite depths of both phases, an analytical solution is obtained. The dispersion relation is derived and its properties are investigated.     

Numerical solution of the two‐phase expansion of a metastable flashing liquid jet using the dispersion‐controlled dissipative scheme

This paper presents a study of the stationary phenomenon of superheated or metastable liquid jets, flashing into a two‐dimensional axisymmetric domain, while in the two‐phase region. In general, the phenomenon starts off when a high‐pressure, high‐temperature liquid jet emerges from a small nozzle or orifice expanding into a low‐pressure chamber, below its saturation pressure taken at the injection temperature. As the process evolves, crossing the saturation curve, one observes that the fluid remains in the liquid phase reaching a superheated condition. Then, the liquid undergoes an abrupt phase change by means of an oblique evaporation wave. Across this phase change the superheated liquid becomes a two‐phase high‐speed mixture in various directions, expanding to supersonic velocities. In order to reach the downstream pressure, the supersonic fluid continues to expand, crossing a complex bow shock wave. The balance equations that govern the phenomenon are mass conservation, momentum conservation, and energy conservation, plus an equation‐of‐state for the substance. A false‐transient model is implemented using the shock capturing scheme: dispersion‐controlled dissipative (DCD), which was used to calculate the flow conditions as the steady‐state condition is reached. Numerical results with computational code DCD‐2D v1 have been analyzed. Copyright © 2009 John Wiley & Sons, Ltd.     

Free surface flow impact on a vertical wall: a numerical assessment

The sudden impact of a free surface flow upon a solid wall is a common occurrence in many situations in nature and technology. The design of marine structures is probably the most obvious example, but also river and dam hydraulics as well as the necessity of understanding flood and debris flow-induced damage have led to theoretical and experimental work on the mechanism of fluid slamming loads. This is therefore a very old and rich research field, which has not yet reached full maturity, so that semi-empirical methods in design practice are still the rule in many sectors. Up-to-date CFD technology with both Eulerian and Lagrangian approaches is employed to investigate highly non-stationary fluid impact on a solid wall. The development of the pressure wave produced by the impact is examined as it propagates and interacts with the fluid boundaries, as well as the subsequent build-up of high-pressure gradients of high fluid velocities. The geometry and the velocity field of the problem considered are very simple, but the results seem to provide new insight, in particular, into the connection between phenomena with different timescales.     

Initiation of explosive boiling of a droplet with a shock wave

The role of incident shock waves in the initiation of vapor explosions in volatile liquid hydrocarbons has been investigated. Experiments were carried out on single droplets (1–2 mm diameter) immersed in a host fluid and heated to temperatures at or near the limit of superheat. Shocks generated by spark discharge were directed at previously nonevaporating drops as well as at drops boiling stably at high pressure. Explosive boiling is triggered in previously nonevaporating drops only if the drop temperature is above a threshold temperature that is near the superheat limit. Interaction of a shock with a stably boiling drop immediately causes a transition to violent unstable boiling in which fine droplets are torn from the evaporating interface, generating a two-phase flow downstream. On the previously nonevaporating interface between the drop and the host liquid, multiple nucleation sites appear which grow rapidly and coalesce. Overpressures generated in the surrounding fluid during bubble collapse may reach values on the same level as the pressure jump across the shock wave that initiated the explosive boiling. A simple calculation is given, which suggests that shock focusing may influence the location at which unstable boiling is initiated.     

Propagation of a pressure wave in two-phase flow with very high void fraction

An experimental and theoretical study of a finite amplitude pressure wave propagating through a two-phase media of about 0.9999–0.99999 void fraction is performed. This two-phase media consists of many parallel liquid films in a gas. The films are perpendicular to the wave propagation direction and result in a two-phase fluid of extremely high void fraction. Experiments are done in a vertical shock tube and show that the shock wave is broken down into an initial sharply rising wave and a second gradually rising wave. The velocity of the first wave agrees well with the theoretical prediction assuming an adiabatic thermal equilibrium change, which approaches the gas sonic velocity in the two-phase flow in the low void fraction region. The second wave is caused by the complex reflection and destruction of the waves.     

Simulations of two-phase flow distribution in communicating parallel channels for a PEM fuel cell

Numerical simulations utilizing computational fluid dynamics (CFD) with a volume of fluid (VOF) method has been employed to investigate two-phase flow distribution in inter-connected parallel flow channels. The interconnections resemble gas and liquid communications in fuel cell flow fields due to the inherent or artificial structures of gas diffusion layers (GDLs). The simulation results showed that communication between parallel channels could have a great impact on the two-phase flow pattern, gas and water distribution and flow maldistribution. Wide communication channels provide a pathway for gas to short-circuit the liquid, leading to a worsened gas flow distribution. However, when the communication channels are narrow enough, they are helpful for mitigating the flow maldistribution by redistributing the liquid among the parallel flow channels through the communication channels. The simulation results were also verified by comparing the predicted and measured normalized pressure drop and the gas flow ratios at the entrance section of experimental parallel channels.