Hydrodynamic Phenomena During High-Speed Collision Between Liquid Droplet and Rigid Plane

1973 ◽  
Vol 95 (2) ◽  
pp. 276-292 ◽  
Author(s):  
Yen C. Huang ◽  
F. G. Hammitt ◽  
W-J Yang

The dynamics of high-speed impact between a compressible water droplet and a rigid solid surface is investigated analytically. The purpose of the study is to examine the mechanism leading to the erosion of a material due to liquid impingement. A Compressible-Cell-and-Marker (ComCAM) numerical method is developed to solve the differential equations governing the unsteady, two-dimensional liquid-solid impact phenomena. The method is designed to solve this unsteady portion up until the flow reasonably approaches the steady-state solution. The validity of the method is confirmed by comparing its numerical results with the idealized exact solution for the classical one-dimensional liquid impact problem. The accuracy of the numerical reresults is found to be very good in that only slight numerical oscillations occur. Viscosity and surface tension are neglected as seems resaonable with the relatively large drops and high velocities considered. Pressure and velocity distributions are solved as a function of time. The deformation of a drop is also recorded for three different shapes: cylindrical, spherical, and a combination of the two, which may more closely model the actual droplet shapes to be encountered in such impacts. Typical liquid impact Mach numbers of 0.2 and 0.5 (sonic velocity referred to water) were studied. Thus impact velocities of about 980 and 2450 fps are considered. Compression predominates during the early stages of the impact, while rarefaction governs later, during which time the radial lateral flow velocity exceeds the initial impact velocity. The reflection of compression waves and the lateral flow leads to the possibility of cavitation within the drop, due to the consequent generation of negative pressures, exists. The maximum pressure calculated in this two-dimensional liquid impact problem is found to be less than the one-dimensional maximum pressure for all three different droplets in various degrees. It is found that droplet shape impact angle and liquid impact Mach number are the only important parameters of the problem for a flat fully-rigid target surface. As more time elapses, i.e., up to 2–3 μsec for a 2.0 mm-dia drop, the maximum pressure shifts from the center of the contact area radially outward, while the pressure at the center attenuates rapidly toward conventional stagnation pressure.

In the initial stage of liquid-drop impact, the contact region expands faster than the wave speed in the liquid. This causes compressible behaviour in the liquid, and high transient pressures. High-velocity jetting results when the wave motion in the liquid overtakes the expanding contact edge and moves up the free surface of the drop. The detailed pressure fields in this early time history of impact have been calculated by Lesser ( Proc . R . Soc . Lond . 377, 289 (1981)) for both two and three-dimensional liquid masses and for targets of finite admittance. An important result is that the edge pressures exceed the central ‘water-hammer’ pressure 3ρ 0 CU i and at the time of shock-detachment approach ca . 3ρ 0 CU i . At this stage the edge pressures, for both spherical drops and two-dimensional liquid wedges, depend only on the impact velocity and the instantaneous angle between the liquid and solid surfaces. This suggests that the essential features of the early stage of liquid impact can be usefully studied by producing impacts with two-dimensional liquid wedges, and predicted data for pressures, shock angles and velocities are presented. Experiments are described for producing impacts with preformed shapes by using water-gelatine mixtures and observing the impact events with high-speed photography. The results confirm the main features of the model and give information on edge pressures, jetting, cavitation in the liquid and the effect of the admittance of the solid. The relevance of the results to the damage and erosion of materials subjected to liquid impact is discussed. In particular, it is possible to explain the apparently low damage-threshold of some materials, the form of damage and its development with repeated impact. The study highlights the importance of the detailed surface geometry in the region of contact.


Author(s):  
Xintian Liu ◽  
Yang Qu ◽  
Xiaobing Yang ◽  
Yongfeng Shen

Background:: In the process of high-speed driving, the wheel hub is constantly subjected to the impact load from the ground. Therefore, it is important to estimate the fatigue life of the hub in the design and production process. Objective:: This paper introduces a method to study the fatigue life of car hub based on the road load collected from test site. Methods:: Based on interval analysis, the distribution characteristics of load spectrum are analyzed. The fatigue life estimation of one - dimensional and two - dimensional load spectra is compared by compiling load spectra. Results:: According to the S-N curve cluster and the one-dimensional program load spectrum, the estimated range fatigue life of the hub is 397,100 km to 529,700 km. For unsymmetrical cyclic loading, each level means and amplitude of load were obtained through the Goodman fatigue empirical formula, and then according to S-N curve clusters in the upper and lower curves and two-dimensional program load spectrum, estimates the fatigue life of wheel hub of the interval is 329900 km to 435200 km, than one-dimensional load spectrum fatigue life was reduced by 16.9% - 17.8%. Conclusion:: This paper lays a foundation for the prediction of fatigue life and the bench test of fatigue durability of auto parts subjected to complex and variable random loads. At the same time, the research method can also be used to estimate the fatigue life of other bearing parts or high-speed moving parts and assemblies.


2000 ◽  
Author(s):  
Songbin Wei ◽  
Imin Kao

Abstract In wiresaw manufacturing process where thin wire moving at high speed is pushed onto ingot to produce slices of wafer, the wire is constrained by two wafer walls as it slices into the ingot. In this paper, we investigate the vibration of such wire under the constraints of wafer walls. To address this problem, the model for wire vibration with impact to wafer walls is developed. The equation of motion is discretized using the Galerkin’s method. The principle of impulse and momentum is utilized to solve the impact problem. The results of analysis and simulation indicate that the response under a pointwise sinusoidal excitation is neither periodical nor symmetric with respect to the horizontal axis, due to the excitation from the impact. The wire vibration behavior is affected dramatically by the wafer wall constraints.


Author(s):  
V. Vlasenko ◽  
A. Shiryaeva

New quasi-two-dimensional (2.5D) approach to description of three-dimensional (3D) flows in ducts is proposed. It generalizes quasi-one-dimensional (quasi-1D, 1.5D) theories. Calculations are performed in the (x; y) plane, but variable width of duct in the z direction is taken into account. Derivation of 2.5D approximation equations is given. Tests for verification of 2.5D calculations are proposed. Parametrical 2.5D calculations of flow with hydrogen combustion in an elliptical combustor of a high-speed aircraft, investigated within HEXAFLY-INT international project, are described. Optimal scheme of fuel injection is found and explained. For one regime, 2.5D and 3D calculations are compared. The new approach is recommended for use during preliminary design of combustion chambers.


Fluids ◽  
2020 ◽  
Vol 5 (4) ◽  
pp. 168 ◽  
Author(s):  
Agostino Lauria ◽  
Giancarlo Alfonsi ◽  
Ali Tafarojnoruz

Ski jump spillways are frequently implemented to dissipate energy from high-speed flows. The general feature of this structure is to transform the spillway flow into a free jet up to a location where the impact of the jet creates a plunge pool, representing an area for potential erosion phenomena. In the present investigation, several tests with different ski jump bucket angles are executed numerically by means of the OpenFOAM® digital library, taking advantage of the Reynolds-averaged Navier–Stokes equations (RANS) approach. The results are compared to those obtained experimentally by other authors as related to the jet length and shape, obtaining physical insights into the jet characteristics. Particular attention is given to the maximum pressure head at the tailwater. Simple equations are proposed to predict the maximum dynamic pressure head acting on the tailwater, as dependent upon the Froude number, and the maximum pressure head on the bucket. Results of this study provide useful suggestions for the design of ski jump spillways in dam construction.


2018 ◽  
Vol 5 (12) ◽  
pp. 181101 ◽  
Author(s):  
Y. Tatekura ◽  
M. Watanabe ◽  
K. Kobayashi ◽  
T. Sanada

The prime objective of this study is to answer the question: How large is the pressure developed at the instant of a spherical liquid droplet impact on a solid surface? Engel first proposed that the maximum pressure rise generated by a spherical liquid droplet impact on a solid surface is different from the one-dimensional water-hammer pressure by a spherical shape factor (Engel 1955 J. Res. Natl Bur. Stand. 55 (5), 281–298). Many researchers have since proposed various factors to accurately predict the maximum pressure rise. We numerically found that the maximum pressure rise can be predicted by the combination of water-hammer theory and the shock relation; then, we analytically extended Engel’s elastic impact model, by realizing that the progression speed of the contact between the gas–liquid interface and the solid surface is much faster than the compression wavefront propagation speed at the instant of the impact. We successfully correct Engel’s theory so that it can accurately provide the maximum pressure rise at the instant of impact between a spherical liquid droplet and solid surface, that is, no shape factor appears in the theory.


2011 ◽  
Vol 56 (4) ◽  
pp. 1-12 ◽  
Author(s):  
K. Richter ◽  
A. Le Pape ◽  
T. Knopp ◽  
M. Costes ◽  
V. Gleize ◽  
...  

A joint comprehensive validation activity on the structured numerical method elsA and the hybrid numerical method TAU was conducted with respect to dynamic stall applications. To improve two-dimensional prediction, the influence of several factors on the dynamic stall prediction was investigated. The validation was performed for three deep dynamic stall test cases of the rotor blade airfoil OA209 against experimental data from two-dimensional pitching airfoil experiments, covering low-speed and high-speed conditions. The requirements for spatial discretization and for temporal resolution in elsA and TAU are shown. The impact of turbulence modeling is discussed for a variety of turbulence models ranging from one-equation Spalart–Allmaras-type models to state-of-the-art, seven-equation Reynolds stress models. The influence of the prediction of laminar/turbulent boundary layer transition on the numerical dynamic stall simulation is described. Results of both numerical methods are compared to allow conclusions to be drawn with respect to an improved prediction of dynamic stall.


Author(s):  
Duohuan Wu ◽  
Jing Wang ◽  
Peiran Yang ◽  
Ton Lubrecht

In this study, the effect of oil starvation on isothermal elastohydrodynamic lubrication of an impact motion is explored with the aid of numerical techniques. During the impact process, on comparison with the fully lubricated results, the pressure and film thickness are much lower and the entrapped film shape does not happen. The rebound is delayed by the oil starvation assumption. During the rebound process, a periphery entrapment is seen in the starved film thickness distribution. Under the starved condition, the maximum pressure gradient is higher. The central film thickness and minimum film thickness exhibit different variations compared with the results by fully flooded assumption.


The type of stress pulse produced when a liquid mass strikes a solid at high velocity is first examined. Compressible behaviour, giving rise to a sharp peak of pressure, is found to occur in the initial stages of the impact. The duration of this peak depends on the dimensions and impact velocity of the liquid mass, and also on the compressible wave velocity for the liquid. A comparison is made with pulses produced by solid/solid impact and by the detonation of small quantities of explosive. Both the high-speed liquid impact and the explosive loading give intense pulses of duration only a few microseconds. A solid/solid impact has, by comparison, a much longer impact time of the order of hundreds of microseconds. The fracture of glasses and hard polymers using these three types of loading is described. The development of fracture is followed by high-speed photography. Differences in the modes of fracture are attributed to variations in the shape and duration of the applied stress pulses. Short circumferential fractures produced around the loaded area in liquid impact and explosive loading are shown to be initiated by the Rayleigh surface wave at points where flaws existed. More complex fracture patterns on the front surfaces of plates are due to the reinforcement of the surface wave with components of stress waves reflected from the back surface. A combination of impact loading and etching makes it possible to investigate the distribution and depths of flaws, their role in the fracture process, and the effect which etching has upon them. The observation on the deformation produced in solids by liquid impact has practical significance in the problem of supersonic aircraft flying through rain and in the erosion of turbine blades moving at high velocity through wet steam.


2021 ◽  
Author(s):  
Rasesh Dalal

A simplified two-dimensional numerical model of a window with a between-panes pleated blind has been developed using commercial computational fluid dynamics software. Knowledge of the effect of blinds on the free convection is important for understanding and predicting the impact of shading devices on the overall thermal performance of a window. Numerical results have been obtained for three fill gases (air, argon and krypton) and several blind geometries over the Rayleigh number range 10³ [less than or equal to]Ra[less than or equal to]10[to the exponent of 5]. The results show that pleated blinds can have a strong effect on window thermal performance. In addition, it has been shown that the data from a convection-only model can be combined with a simplified one-dimensional model to estimate the overall U-value of the enclosure. Using this procedure, the convection data can be applied to a window/blind assembly with arbitrary radiation parameters.


Sign in / Sign up

Export Citation Format

Share Document