Air-assisted atomization of liquid jets in varying levels of turbulence

2014 ◽  
Vol 764 ◽  
pp. 95-132 ◽  
Author(s):  
A. Kourmatzis ◽  
A. R. Masri

AbstractAir-assisted primary atomization is investigated in a configuration where liquid is injected in a turbulent gaseous jet flow both within as well as outside of the potential core. Cases are studied where the injection point is moved within the flow to maintain a range of constant gaseous mean velocities but changing local fluctuating velocity root-mean-square (r.m.s.) levels. Over a range of mean conditions, this allows for a systematic understanding of both the effects of gas-phase turbulence and mean shear on primary break-up independently. Extensive data is obtained and analysed from laser Doppler anemometry/phase Doppler anemometry, high-speed microscopic backlit imaging and advanced image processing. It is found that the ratio of the turbulent Weber number $\mathit{We}^{\prime }$ to the mean Weber number $\mathit{We}$ is a relevant parameter as is the turbulence intensity. The primary break-up length is found to be heavily influenced not only by the mean velocity, but also by the turbulence level and the mass fuel to air ratio. Above a particular threshold intensity level the break-up time changes in proportion to the change in the integral time scale of the flow. In addition, it is found that regardless of diameter and turbulent flow conditions at the liquid jet, the final size of ligaments converges to a value which is of the order of the measured primary instability wavelength (${\it\lambda}_{1}$). In contrast, cases of different turbulence intensity show the mean of droplet sizes diverging as the spray is advected downstream and this is because droplets are generated from ligaments, the latter of which are subjected both to Rayleigh–Taylor instabilities and turbulent fluctuations. This contribution, for the first time, examines the theoretical applicability of the Rayleigh–Taylor instability in flows where the turbulence is substantial with respect to the mean flow. It is shown that for high turbulence intensities a full theoretical reconstruction of the measured final droplet size distribution is possible from a probability density function of model Rayleigh–Taylor wavelengths (${\it\lambda}_{RT}$). In agreement with the literature (Varga et al. J. Fluid Mech., vol. 497, 2003, pp. 405–434), mean droplet sizes are found to be equal to a mean theoretical Rayleigh–Taylor wavelength normalized by a particular constant value. This, however, is only true for local turbulence intensities less than ${\sim}25\,\%$, or for ratios of the turbulent Weber number to mean Weber number ($\mathit{We}^{\prime }/\mathit{We}$) of less than ${\sim}6\,\%$. Above this, the normalization value is no longer constant, but increases with $\mathit{We}^{\prime }/\mathit{We}$. Finally, the instability wavelengths can be used as part of an approximation that estimates the total number of objects formed after break-up, where the object number is found to be dictated by a balance of both mean flow conditions and local turbulence.

Energies ◽  
2021 ◽  
Vol 14 (12) ◽  
pp. 3625
Author(s):  
Jon Hardwick ◽  
Ed B. L. Mackay ◽  
Ian G. C. Ashton ◽  
Helen C. M. Smith ◽  
Philipp R. Thies

Numerical modeling of currents and waves is used throughout the marine energy industry for resource assessment. This study compared the output of numerical flow simulations run both as a standalone model and as a two-way coupled wave–current simulation. A regional coupled flow-wave model was established covering the English Channel using the Delft D-Flow 2D model coupled with a SWAN spectral wave model. Outputs were analyzed at three tidal energy sites: Alderney Race, Big Roussel (Guernsey), and PTEC (Isle of Wight). The difference in the power in the tidal flow between coupled and standalone model runs was strongly correlated to the relative direction of the waves and currents. The net difference between the coupled and standalone runs was less than 2.5%. However, when wave and current directions were aligned, the mean flow power was increased by up to 7%, whereas, when the directions were opposed, the mean flow power was reduced by as much as 9.6%. The D-Flow Flexible Mesh model incorporates the effects of waves into the flow calculations in three areas: Stokes drift, forcing by radiation stress gradients, and enhancement of the bed shear stress. Each of these mechanisms is discussed. Forcing from radiation stress gradients is shown to be the dominant mechanism affecting the flow conditions at the sites considered, primarily caused by dissipation of wave energy due to white-capping. Wave action is an important consideration at tidal energy sites. Although the net impact on the flow power was found to be small for the present sites, the effect is site specific and may be significant at sites with large wave exposure or strong asymmetry in the flow conditions and should thus be considered for detailed resource and engineering assessments.


Author(s):  
P. A. Beau ◽  
T. Me´nard ◽  
R. Lebas ◽  
A. Berlemont ◽  
S. Tanguy ◽  
...  

The main objective of our work is to develop direct numerical simulation tools for the primary break up of a jet. Results can help to determine closure relation in the ELSA model [1] which is based on a single-phase Eulerian model and on the transport equation for the mean liquid/gas interface density in turbulent flows. DNS simulations are carried out to obtain statistical information in the dense zone of the spray where nearly no experimental data are available. The numerical method should describe the interface motion precisely, handle jump conditions at the interface without artificial smoothing, and respect mass conservation. We develop a 3D code [2], where interface tracking is ensured by Level Set method, Ghost Fluid Method [3] is used to capture accurately sharp discontinuities, and coupling between Level Set and VOF methods is used for mass conservation [4]. Turbulent inflow boundary conditions are generated through correlated random velocities with a prescribed length scale. Specific care has been devoted to improve computing time with MPI parallelization. The numerical methods have been applied to investigate physical processes that are involved in the primary break up of an atomizing jet. The chosen configuration is close as possible of Diesel injection (Diameter D = 0.1 mm, Velocity = 100m/s, Liquid density = 696kg/m3, Gas density = 25kg/m3). Typical results will be presented. From the injector nozzle, the turbulence initiates some perturbations on the liquid surface, that are enhanced by the mean shear between the liquid jet and the surrounding air. The interface becomes very wrinkled and some break-up is initiated. The induced liquid parcels show a wide range of shapes. Statistics are carried out and results will be provided for liquid volume fraction, liquid/gas interface density, and turbulent correlations.


Water ◽  
2019 ◽  
Vol 12 (1) ◽  
pp. 118 ◽  
Author(s):  
Hossein Hamidifar ◽  
Alireza Keshavarzi ◽  
Paweł M. Rowiński

Trees have been used extensively by river managers for improving the river environment and ecology. The link between flow hydraulics, bed topography, habitat availability, and organic matters is influenced by vegetation. In this study, the effect of trees on the mean flow, bed topography, and bed shear stress were tested under different flow conditions. It was found that each configuration of trees produced particular flow characteristics and bed topography patterns. The SR (single row of trees) model appeared to deflect the maximum velocity downstream of the bend apex toward the inner bank, while leading the velocity to be more uniformly distributed throughout the bend. The entrainment of sediment particles occurred toward the area with higher values of turbulent kinetic energy (TKE). The results showed that both SR and DR (double rows of trees) models are effective in relieving bed erosion in sharp ingoing bends. The volume of the scoured bed was reduced up to 70.4% for tests with trees. This study shows the effectiveness of the SR model in reducing the maximum erosion depth.


2019 ◽  
Vol 22 (8) ◽  
pp. 1977-1987
Author(s):  
Xu Wang ◽  
Huaqiang Li ◽  
Zengshun Chen ◽  
Yuanhao Qian ◽  
Yanru Wang ◽  
...  

During landfall of Typhoon Haikui in Eastern China in 2012, ground level wind data were recorded using a smart monitoring system installed on JiuBao Bridge in Hangzhou, China. This article documents the mean flow and turbulence characteristics from data recorded during the storm. The results show that both turbulence intensity and gust factor decrease with the increase in the mean wind velocity. However, as the mean wind velocity increases, this trend gradually attenuates. The peak factor distribution with gust averaging time duration derived with the Typhoon Haikui data agrees well with the Durst curve. However, the longitudinal gust factor derived from the typhoon wind-speed record in this study is higher compared with the curves proposed by Durst and Krayer-Marshall. Analyses of the gust factor distribution with the turbulence intensity during the passage of the storm reveal a similarity to the empirical curves of Ishizaki and Choi. Results show that the relationship between lateral turbulence and gust factors can be well represented by a quadratic polynomial. Turbulence scale increased with mean wind velocity. The values of autocorrelation coefficients in longitudinal direction are larger than those in lateral direction. There has no obvious dependency of cross-correlation coefficients with mean wind velocity. In general, the wind characteristics in this study are shown to be very similar to those of winds under normal circumstance.


Author(s):  
C. J. Lea ◽  
A. P. Watkins

A study is made here of the application of a differential stress model (DSM) of turbulence to flows in two model reciprocating engines. For the first time this study includes compressive effects. An assessment between DSM and k-ɛ results is made comparing with laser Doppler anemometry experimental data of the mean flow and turbulence intensity levels during intake and compression strokes. A well-established two-dimensional finite-volume computer code is employed. Two discretization schemes are used, namely the HYBRID scheme and the QUICK scheme. The latter is found to be essential if differentiation is to be made between the turbulence models. During the intake stroke the DSM results are, in general, similar to the k-ɛ results in comparison to the experimental data, except for the turbulence levels, which the DSM seriously underpredicts. This is in contrast to a parallel set of calculations of steady in-flow, which showed significant gains from using the DSM, particularly at the turbulence field level. The increased number of grid lines employed in those calculations contribute to this apparent difference between steady and unsteady flows, but cycle- to-cycle variations are more likely to be the primary cause, resulting in too high levels of turbulence intensity being measured. However, during the compression stroke the DSM returns vastly superior results to the k-ɛ model at both the mean flow and turbulence intensity levels. This is because the DSM generates an anisotropic shear stress field during the early stages of compression that suppresses the main vortical structure, in line with the experimental data.


1996 ◽  
Vol 2 (3) ◽  
pp. 209-218 ◽  
Author(s):  
Richard B. Rivir ◽  
Mingking K. Chyu ◽  
Paul K. Maciejewski

Hot-wire measurements of the mean flow, turbulence characteristics, and integral scale in a square channel roughened with transverse ribs mounted on two opposing sidewalls are presented for three rib configurations: single rib, in-line multiple ribs, and staggered multiple ribs. Test conditions for multiple ribs use p/H = 10, H/D 0.17, andRe⁡D23,000. Measured results highlight the spatial distribution and evolution of turbulence intensity and integral scale from the flow entrance of the first period to the developed regime near the exit of the third period. The highly turbulent, shear layer initiated near the trailing upper-edge of a rib elevates the turbulence level in the mainstream of the channel. The magnitude of turbulence intensity in the channel core rises from 0.7% in the approaching flow to about 20–25% near the exit of the third period. The integral scale dominating the mainstream flow increases from approximately one-half the rib-height, 0.5H, in the approaching flow to 1.5-2.5H behind the first rib and further downstream.


2020 ◽  
Vol 5 (10) ◽  
pp. 1199-1203
Author(s):  
Md. Mosharrof Hossain ◽  
Muhammed Hasnain Kabir Nayeem ◽  
Dr. Md Abu Taher Ali

In this investigation experiment was carried out in 80 mm diameter swirling pipe jet, where swirl was generated by attaching wedge-shaped helixes in the pipe. All measurements were taken at Re 5.3e4. In the plain pipe jet the potential core was found to exist up to x/D=5 but in the swirling jet there was no existence of potential core. The mean velocity profiles were found to be influenced by the presence of wedge-shaped helixes in the pipe. The velocity profiles indicated the presence of sinusoidal flow field in the radial direction existed only in the near field of the jet. This flow field died out after x/D=3 and the existence of jet flow diminished after x/D=5.


2019 ◽  
Vol 863 ◽  
pp. 454-493 ◽  
Author(s):  
Qian-Cheng Wang ◽  
Zhen-Guo Wang ◽  
Ming-Bo Sun ◽  
Rui Yang ◽  
Yu-Xin Zhao ◽  
...  

Direct numerical simulation is conducted to uncover the response of a supersonic turbulent boundary layer to streamwise concave curvature and the related physical mechanisms at a Mach number of 2.95. Streamwise variations of mean flow properties, turbulence statistics and turbulent structures are analysed. A method to define the boundary layer thickness based on the principal strain rate is proposed, which is applicable for boundary layers subjected to wall-normal pressure and velocity gradients. While the wall friction grows with the wall turning, the friction velocity decreases. A logarithmic region with constant slope exists in the concave boundary layer. However, with smaller slope, it is located lower than that of the flat boundary layer. Streamwise varying trends of the velocity and the principal strain rate within different wall-normal regions are different. The turbulence level is promoted by the concave curvature. Due to the increased turbulence generation in the outer layer, secondary bumps are noted in the profiles of streamwise and spanwise turbulence intensity. Peak positions in profiles of wall-normal turbulence intensity and Reynolds shear stress are pushed outward because of the same reason. Attributed to the Görtler instability, the streamwise extended vortices within the hairpin packets are intensified and more vortices are generated. Through accumulations of these vortices with a similar sense of rotation, large-scale streamwise roll cells are formed. Originated from the very large-scale motions and by promoting the ejection, sweep and spanwise events, the formation of large-scale streamwise roll cells is the physical cause of the alterations of the mean properties and turbulence statistics. The roll cells further give rise to the vortex generation. The large number of hairpin vortices formed in the near-wall region lead to the improved wall-normal correlation of turbulence in the concave boundary layer.


2004 ◽  
Vol 126 (2) ◽  
pp. 257-265 ◽  
Author(s):  
Qiang Zhang ◽  
Sang Woo Lee ◽  
Phillip M. Ligrani

The effects of surface roughness on the aerodynamic performance of turbine airfoils are investigated with different inlet turbulence intensity levels of 0.9%, 5.5% and 16.2%. Three symmetric airfoils, each with the same shape and exterior dimensions, are employed with different rough surfaces. The nonuniform, irregular, 3-D roughness is characterized using the equivalent sand grain roughness size. Mach numbers along the airfoil range from 0.4 to 0.7. Chord Reynolds numbers based on inlet and exit flow conditions are 0.54×106 and 1.02×106, respectively. The contributions of varying surface roughness and turbulence intensity level to aerodynamic losses, Mach number profiles, normalized kinetic energy profiles, and Integrated Aerodynamics Losses (IAL) are quantified. Results show that effects of changing the surface roughness condition on IAL values are substantial, whereas the effects of different inlet turbulence intensity levels are generally relatively small.


Sign in / Sign up

Export Citation Format

Share Document