scholarly journals Progress on Understanding Rayleigh–Taylor Flow and Mixing Using Synergy Between Simulation, Modeling, and Experiment

2020 ◽  
Vol 142 (12) ◽  
Author(s):  
Oleg Schilling

Abstract Simultaneous advances in numerical methods and computing, theoretical techniques, and experimental diagnostics have all led independently to better understanding of Rayleigh–Taylor (RT) instability, turbulence, and mixing. In particular, experiments have provided significant motivation for many simulation and modeling studies, as well as validation data. Numerical simulations have also provided data that is not currently measurable or very difficult to measure accurately in RT unstable flows. Thus, simulations have also motivated new measurements in this class of buoyancy-driven flows. This overview discusses simulation and modeling studies synergistic with experiments and examples of how experiments have motivated simulations and models of RT instability, flow, and mixing. First, a brief summary of measured experimental and calculated simulation quantities, of experimental approaches, and of issues and challenges in the simulation and modeling of RT experiments is presented. Implicit large-eddy, direct numerical, and large-eddy simulations validated using RT experimental data are then discussed. This is followed by a discussion of modeling using analytical, modal, buoyancy–drag, and turbulent transport models of RT mixing experiments. The discussion will focus on three-dimensional RT mixing arising from multimode perturbations. Finally, this focused review concludes with a perspective on future simulation, modeling, and experimental directions for further research. Research in simulation and modeling of RT unstable flows, coupled with experiments, has made significant progress over the past several decades. This overview serves as an opportunity to both discuss progress and to stimulate future research on simulation and modeling of this unique class of hydrodynamically unstable turbulent flows.

2005 ◽  
Vol 122 (3) ◽  
pp. 42-55
Author(s):  
Jorge BARATA

The present paper presents a numerical study on evaporating droplets injected through a turbulent cross-stream. Several models have been used with more or less success to describe similar phenomena, but much of the reported work deals only with sprays in stagnant surroundings. The ultimate goal of this study is to develop an Eulerian/Lagragian approach to account for turbulent transport, dispersion, evaporation and coupling between both processes in practical spray injection systems, which usually include air flows in the combustion chamber like swirl, tumble and squish in I.C. engines or crossflow in gas turbines. In this work a method developed to study isothermal turbulent dispersion is extended to the case of an array of evaporating droplets through a crossflow, and the performance of two different evaporation models widely used is investigated. The convection terms were evaluated using the hybrid or the higher order QUICK scheme. The dispersed phase was treated using a Lagrangian reference frame. The differences between the two evaporation models and its applicability to the present flow are analysed in detail. During the preheating period of the Chen and Pereira [1] model the droplets are transported far away from the injector by the crossflow, while with the Sommerfeld [2] formulation for evaporation the droplet has a continuous variation of the diameter. This result has profound implications on the results because the subsequent heat transfer and turbulent dispersion is extremely affected by the size of the particles (or droplets). As a consequence, droplet diameter, temperature and mass fraction distributions were found to be strongly dependent on the evaporation model used. So, a new formulation that takes into account also the transport of the evaporating droplets needs to be developed if practical injection systems are to be simulated. Also, in order to better evaluate and to improve the vaporization models more detailed measurements of three-dimensional configurations are required.


2020 ◽  
Author(s):  
Rui Zhang ◽  
Chengwen Zhong ◽  
Sha Liu ◽  
Congshan Zhuo

Abstract In this paper, we introduce the incompressible discrete Maxwellian equilibrium distribution function and external forces into the two-stage third-order Discrete Unified Gas-Kinetic Scheme (DUGKS) for simulating low-speed incompressible turbulent flows with forcing term. The Wall-Adapting Local Eddy-viscosity (WALE) and Vreman sub-grid models for Large-Eddy Simulations (LES) of wall-bounded turbulent flows are coupled within the present framework. In order to simulate the three-dimensional turbulent flows associated with great computational cost, a parallel implementation strategy for the present framework is developed, and is validated by three canonical wall-bounded turbulent flows, viz., the fully developed turbulent channel flow at a friction Reynolds number (Re) about 180, the turbulent plane Couette flow at a friction Re number about 93 and three-dimensional lid-driven cubical cavity flow at a Re number of 12000. The turbulence statistics are computed by the present approach with both WALE and Vreman models, and their predictions match precisely with each other. Especially, the predicted flow physics of three-dimensional lid-driven cavity are consistent with the description from abundant literatures. While, they have small discrepancies in comparison to the Direct Numerical Simulation (DNS) due to the relatively low grid resolution. The present numerical results verify that the present two-stage third-order DUGKS-based LES method is capable for simulating inhomogeneous wall-bounded turbulent flows and getting reliable results with relatively coarse grids.


Author(s):  
Bakhtier Farouk

Abstract A three-dimensional time dependent computational fluid dynamic (CFD) study of laminar and turbulent thermal mixing of two flows entering a 90° T-junction pipe is presented. The two incoming flows (both liquids) in the T-junction enter the flow domain with different inlet velocities, and temperatures. Water flow is considered in both the horizontal pipe and the vertical pipe. Inlet temperature differences and temperature dependent thermophysical properties are considered. Large eddy simulations (LES) with sub-grid scale (SGS) modeling were considered for the simulation of the turbulent cases. The flow characteristics, and thermal mixing behaviors and detailed mixing structures were simulated, and they showed that thermal mixing of the two streams are closely affected by the inlet conditions of the two streams and the inlet thermophysical properties of the two streams.


Author(s):  
K. Ya. Yakubovskiy ◽  
◽  
A. B. Lebedev ◽  
P. D. Toktaliev ◽  
◽  
...  

The effect of initial nonuniformity and fluctuations of fuel concentration on the combustion stability and NOx and CO emission in the model combustion chamber was analyzed with the use of previously developed simple and computationally inexpensive Large Eddy Simulation (LES) methodology for simulation of three-dimensional unsteady turbulent flows with premixed combustion of methane-air mixture in low-emission combustion chamber which geometry is represented by channel with the backward facing step. Typical sizes of the combustion chamber, flow parameters, turbulence level, and method of flame front stabilization are close to those of full-sized industrial combustors.


2015 ◽  
Vol 15 (13) ◽  
pp. 7369-7389 ◽  
Author(s):  
J. Picot ◽  
R. Paoli ◽  
O. Thouron ◽  
D. Cariolle

Abstract. In this work, the evolution of contrails in the vortex and dissipation regimes is studied by means of fully three-dimensional large-eddy simulation (LES) coupled to a Lagrangian particle tracking method to treat the ice phase. In this paper, fine-scale atmospheric turbulence is generated and sustained by means of a stochastic forcing that mimics the properties of stably stratified turbulent flows as those occurring in the upper troposphere and lower stratosphere. The initial flow field is composed of the turbulent background flow and a wake flow obtained from separate LES of the jet regime. Atmospheric turbulence is the main driver of the wake instability and the structure of the resulting wake is sensitive to the intensity of the perturbations, primarily in the vertical direction. A stronger turbulence accelerates the onset of the instability, which results in shorter contrail descent and more effective mixing in the interior of the plume. However, the self-induced turbulence that is produced in the wake after the vortex breakup dominates over background turbulence until the end of the vortex regime and controls the mixing with ambient air. This results in mean microphysical characteristics such as ice mass and optical depth that are slightly affected by the intensity of atmospheric turbulence. However, the background humidity and temperature have a first-order effect on the survival of ice crystals and particle size distribution, which is in line with recent studies.


Author(s):  
N. Gourdain ◽  
F. Sicot ◽  
F. Duchaine ◽  
L. Gicquel

A better understanding of turbulent unsteady flows is a necessary step towards a breakthrough in the design of modern compressors. Owing to high Reynolds numbers and very complex geometry, the flow that develops in such industrial machines is extremely hard to predict. At this time, the most popular method to simulate these flows is still based on a Reynolds-averaged Navier–Stokes approach. However, there is some evidence that this formalism is not accurate for these components, especially when a description of time-dependent turbulent flows is desired. With the increase in computing power, large eddy simulation (LES) emerges as a promising technique to improve both knowledge of complex physics and reliability of flow solver predictions. The objective of the paper is thus to give an overview of the current status of LES for industrial compressor flows as well as to propose future research axes regarding the use of LES for compressor design. While the use of wall-resolved LES for industrial multistage compressors at realistic Reynolds number should not be ready before 2035, some possibilities exist to reduce the cost of LES, such as wall modelling and the adaptation of the phase-lag condition. This paper also points out the necessity to combine LES to techniques able to tackle complex geometries. Indeed LES alone, i.e. without prior knowledge of such flows for grid construction or the prohibitive yet ideal use of fully homogeneous meshes to predict compressor flows, is quite limited today.


Author(s):  
C. Hah

A computational procedure based on the compressible Reynolds-averaged Navier-Stokes equation has been developed for the viscous flow through an isolated compressor rotor. The numerical scheme is based on fully conservative control volume formulation and solves various conservation equations in fully elliptic form on the rotating coordinates fixed on the rotor. An algebraic Reynolds stress model is used to describe the turbulent transport terms. The numerical procedure has been applied to predict three-dimensional turbulent flows through two different isolated compressor rotors. The detailed quantitative comparisons with two sets of well-documented data show that the developed computational procedure predicts the viscous flow development over the blading and in the wake with the accuracy satisfactory for most engineering purposes; the computer code can be used for the guidance of advanced rotor design.


2001 ◽  
Vol 447 ◽  
pp. 377-408 ◽  
Author(s):  
Y.-N. YOUNG ◽  
H. TUFO ◽  
A. DUBEY ◽  
R. ROSNER

We investigate the miscible Rayleigh–Taylor (RT) instability in both two and three dimensions using direct numerical simulations, where the working fluid is assumed incompressible under the Boussinesq approximation. We first consider the case of randomly perturbed interfaces. With a variety of diagnostics, we develop a physical picture for the detailed temporal development of the mixed layer: we identify three distinct evolutionary phases in this development, which can be related to detailed variations in the growth of the mixing zone. Our analysis provides an explanation for the observed differences between two- and three-dimensional RT instability; the analysis also leads us to concentrate on the RT models which (i) work equally well for both laminar and turbulent flows, and (ii) do not depend on turbulent scaling within the mixing layer between fluids. These candidate RT models are based on point sources within bubbles (or plumes) and their interaction with each other (or the background flow). With this motivation, we examine the evolution of single plumes, and relate our numerical results (for single plumes) to a simple analytical model for plume evolution.


Sign in / Sign up

Export Citation Format

Share Document