Buoyancy Driven Effects in Formation of Mixing Zones

Author(s):  
Arindam Banerjee

The release of a mass of hydrogen fuel (gas) into the ambient atmosphere results in the transient formation of flammable mixture zones that represent potential fire, explosion and toxic hazards. The formation of mixing zones of air and hydrogen for this simple geometry follows the classical Rayleigh Taylor (R-T) instability, which is induced when a heavy fluid is placed over a light fluid in a gravitational field. Buoyancy driven mixing in such flow configurations is studied by using the Boussinesq approximation and considering the flow to be laminar. However, this approximation is valid only at low Atwood numbers (non-dimensional density differences). As Atwood number increases (>0.1, i.e. large density differences) the Boussinesq approximation is no longer valid and a distinct bubble and spike geometry of Rayleigh-Taylor buoyant plumes is formed. Aside from asymmetry in the flow the Atwood number also affects key parameters such as the growth constants and molecular mix. The effect of initial conditions on the growth rate of turbulent Rayleigh-Taylor (RT) mixing has been studied using carefully formulated numerical simulations. A monotone integrated large-eddy simulation (MILES) using a finite-volume technique was employed to solve the three-dimensional incompressible Euler equations with numerical dissipation for air and hydrogen mixing at Atwood number 0.875. The study also provides preliminary guidelines for reducing the fire and explosion hazards in enclosures where such situations are present.

2011 ◽  
Vol 121-126 ◽  
pp. 3657-3661
Author(s):  
Dun Zhang ◽  
Yuan Zheng ◽  
Ying Zhao ◽  
Jian Jun Huang

Numerical simulation of three-dimensional transient turbulent flow in the whole flow passage of a Francis turbine were based upon the large eddy simulation(LES) technique on Smargorinsky model and sliding mesh technology. The steady flow data simulated with the standard k-εmodel was used as the initial conditions for the unsteady simulation. The results show that LES can do well transient turbulent flow simulation in a Francis turbine with complex geometry. The computational method provides some reference for exploring the mechanism of eddy formation in a complex turbulent of hydraulic machinery.


Atmosphere ◽  
2021 ◽  
Vol 12 (9) ◽  
pp. 1107
Author(s):  
Akshay A. Gowardhan ◽  
Dana L. McGuffin ◽  
Donald D. Lucas ◽  
Stephanie J. Neuscamman ◽  
Otto Alvarez ◽  
...  

Fast and accurate predictions of the flow and transport of materials in urban and complex terrain areas are challenging because of the heterogeneity of buildings and land features of different shapes and sizes connected by canyons and channels, which results in complex patterns of turbulence that can enhance material concentrations in certain regions. To address this challenge, we have developed an efficient three-dimensional computational fluid dynamics (CFD) code called Aeolus that is based on first principles for predicting transport and dispersion of materials in complex terrain and urban areas. The model can be run in a very efficient Reynolds average Navier–Stokes (RANS) mode or a detailed large eddy simulation (LES) mode. The RANS version of Aeolus was previously validated against field data for tracer gas and radiological dispersal releases. As a part of this work, we have validated the Aeolus model in LES mode against two different sets of data: (1) turbulence quantities measured in complex terrain at Askervein Hill; and (2) wind and tracer data from the Joint Urban 2003 field campaign for urban topography. As a third set-up, we have applied Aeolus to simulate cloud rise dynamics for buoyant plumes from high-temperature explosions. For all three cases, Aeolus LES predictions compare well to observations and other models. These results indicate that Aeolus LES can be used to accurately simulate turbulent flow and transport for a wide range of applications and scales.


Author(s):  
Malcolm J. Andrews ◽  
Farzaneh F. Jebrail ◽  
Arindam Banerjee

High Atwood number (non-dimensional density difference) effects in buoyancy driven flows are discussed. Buoyancy driven (natural convection) flows may be treated as Boussinesq for small Atwood number, but as Atwood number increases (>0.1, i.e. large temperature differences) the Boussinesq approximation is no longer valid and a distinct "bubble" and "spike" geometry of Rayleigh-Taylor buoyant plumes is formed. Aside from asymmetry in the flow the Atwood number also affects key turbulent mix parameters such as the molecular mix, and heat transfer coefficients. This paper presents recent experimental work being performed in the buoyancy driven mix laboratory at Texas A&M University, using air/helium as mixing components (upto At ~ 0.5). Corresponding numerical simulations of the experiments performed at Los Alamos is also presented, and future directions for the research discussed.


2010 ◽  
Vol 654 ◽  
pp. 99-139 ◽  
Author(s):  
B. THORNBER ◽  
D. DRIKAKIS ◽  
D. L. YOUNGS ◽  
R. J. R. WILLIAMS

This paper investigates the influence of different three-dimensional multi-mode initial conditions on the rate of growth of a mixing layer initiated via a Richtmyer–Meshkov instability through a series of well-controlled numerical experiments. Results are presented for large-eddy simulation of narrowband and broadband perturbations at grid resolutions up to 3 × 109 points using two completely different numerical methods, and comparisons are made with theory and experiment. It is shown that the mixing-layer growth is strongly dependent on initial conditions, the narrowband case giving a power-law exponent θ ≈ 0.26 at low Atwood and θ ≈ 0.3 at high Atwood numbers. The broadband case uses a perturbation power spectrum of the form P(k) ∝ k−2 with a proposed theoretical growth rate of θ = 2/3. The numerical results confirm this; however, they highlight the necessity of a very fine grid to capture an appropriately broad range of initial scales. In addition, an analysis of the kinetic energy decay rates, fluctuating kinetic energy spectra, plane-averaged volume fraction profiles and mixing parameters is presented for each case.


1995 ◽  
Vol 13 (3) ◽  
pp. 423-440 ◽  
Author(s):  
J. Hecht ◽  
D. Ofer ◽  
U. Alon ◽  
D. Shvarts ◽  
S.A. Orszag ◽  
...  

The nonlinear stage in the growth of the Rayleigh-Taylor instability in three dimensions (3D) is studied using a 3D multimaterial hydrodynamic code. The growth of a single classical 3D square and rectangular modes is compared to the growth in planar and cylindrical geometries and found to be close to the corresponding cylindrical mode, which is in agreement with a new Layzer-type model for 3D bubble growth. The Atwood number effect on the final shape of the instability is demonstrated. Calculations in spherical geometry of the late deceleration stage of a typical ICF pellet have been performed. The different late time shapes obtained are shown to be a result of the initial conditions and the high Atwood number. Finally, preliminary results of calculations of two-mode coupling and random perturbations growth in 3D are presented.


Author(s):  
Adra Benhacine ◽  
Zoubir Nemouchi ◽  
Lyes Khezzar ◽  
Nabil Kharoua

A numerical study of a turbulent plane jet impinging on a convex surface and on a flat surface is presented, using the large eddy simulation approach and the Smagorinski-Lilly sub-grid-scale model. The effects of the wall curvature on the unsteady filtered, and the steady mean, parameters characterizing the dynamics of the wall jet are addressed in particular. In the free jet upstream of the impingement region, significant and fairly ordered velocity fluctuations, that are not turbulent in nature, are observed inside the potential core. Kelvin-Helmholtz instabilities in the shear layer between the jet and the surrounding air are detected in the form of wavy sheets of vorticity. Rolled up vortices are detached from these sheets in a more or less periodic manner, evolving into distorted three dimensional structures. Along the wall jet the Coanda effect causes a marked suction along the convex surface compared with the flat one. As a result, relatively important tangential velocities and a stretching of sporadic streamwise vortices are observed, leading to friction coefficient values on the curved wall higher than those on the flat wall.


Atmosphere ◽  
2021 ◽  
Vol 12 (7) ◽  
pp. 906
Author(s):  
Ivan Bašták Ďurán ◽  
Martin Köhler ◽  
Astrid Eichhorn-Müller ◽  
Vera Maurer ◽  
Juerg Schmidli ◽  
...  

The single-column mode (SCM) of the ICON (ICOsahedral Nonhydrostatic) modeling framework is presented. The primary purpose of the ICON SCM is to use it as a tool for research, model evaluation and development. Thanks to the simplified geometry of the ICON SCM, various aspects of the ICON model, in particular the model physics, can be studied in a well-controlled environment. Additionally, the ICON SCM has a reduced computational cost and a low data storage demand. The ICON SCM can be utilized for idealized cases—several well-established cases are already included—or for semi-realistic cases based on analyses or model forecasts. As the case setup is defined by a single NetCDF file, new cases can be prepared easily by the modification of this file. We demonstrate the usage of the ICON SCM for different idealized cases such as shallow convection, stratocumulus clouds, and radiative transfer. Additionally, the ICON SCM is tested for a semi-realistic case together with an equivalent three-dimensional setup and the large eddy simulation mode of ICON. Such consistent comparisons across the hierarchy of ICON configurations are very helpful for model development. The ICON SCM will be implemented into the operational ICON model and will serve as an additional tool for advancing the development of the ICON model.


Author(s):  
P. J. Wicks

AbstractIn this paper, a model for lateral dispersion in open-channel flow is studied involving a diffusion equation which has a nonlinear term describing the effect of buoyancy. The model is used to investigate the interaction of two buoyant pollutant plumes. An approximate analytic technique involving Hermite polynomials is applied to the resulting PDEs to reduce them to a system of ODEs for the centroids and widths of the two plumes. The ODEs are then solved numerically. A rich variety of behaviour occurs depending on the relative positions, widths and strengths of the initial discharges. It is found that for two plumes of equal strength and width discharged side-by-side, the plumes move apart and the rate of spreading is inhibited by their interaction, whereas when one plume is initially much wider than the other, both plumes tend to drift to the side of the narrower plume. Finally, the PDEs are solved numerically for two sets of initial conditions and a comparison is made with the ODE solutions. Agreement is found to be good.


Sign in / Sign up

Export Citation Format

Share Document