scholarly journals Large eddy simulation of high atwood number rayleigh-taylor mixing

2019 ◽  
Vol 128 ◽  
pp. 08001
Author(s):  
Ilyas Yilmaz

Large eddy simulation of Rayleigh-Taylor instability at high Atwood numbers is performed using recently developed, kinetic energy-conserving, non-dissipative, fully-implicit, finite volume algorithm. The algorithm does not rely on the Boussinesq assumption. It also allows density and viscosity to vary. No interface capturing mechanism is requried. Because of its advanced features, unlike the pure incompressible ones, it does not suffer from the loss of physical accuracy at high Atwood numbers. Many diagnostics including local mole fractions, bubble and spike growth rates, mixing efficiencies, Taylor micro-scales, Reynolds stresses and their anisotropies are computed to analyze the high Atwood number effects. The density ratio dependence for the ratio of spike to bubble heights is also studied. Results show that higher Atwood numbers are characterized by increasing ratio of spike to bubble growth rates, higher speeds of bubble and especially spike fronts, faster development in instability, similarity in late time mixing values, and mixing asymmetry.

Author(s):  
F. F. Grinstein ◽  
A. A. Gowardhan ◽  
J. R. Ristorcelli

Under-resolved computer simulations are typically unavoidable in practical turbulent flow applications exhibiting extreme geometrical complexity and a broad range of length and time scales. An important unsettled issue is whether filtered-out and subgrid spatial scales can significantly alter the evolution of resolved larger scales of motion and practical flow integral measures. Predictability issues in implicit large eddy simulation of under-resolved mixing of material scalars driven by under-resolved velocity fields and initial conditions are discussed in the context of shock-driven turbulent mixing. The particular focus is on effects of resolved spectral content and interfacial morphology of initial conditions on transitional and late-time turbulent mixing in the fundamental planar shock-tube configuration.


Author(s):  
Yunfei Wang ◽  
Huanlong Chen ◽  
Huaping Liu ◽  
Yanping Song ◽  
Fu Chen

An in-house large eddy simulation (LES) code based on three-dimensional compressible N-S equations is used to research the impact of incoming wakes on unsteady evolution characteristic in a low-pressure turbine (LPT) cascade. The Mach number is 0.4 and Reynolds number is 0.6 × 105 (based on the axial chord and outlet velocity). The reduced frequency of incoming wakes is Fred = 0 (without wakes), 0.37 and 0.74. A detailed analysis of Reynolds stresses and turbulent kinetic energy inside the boundary layer has been carried out. Particular consideration is devoted to the transport process of incoming wakes and the intermittent property of the unsteady boundary layer. With the increase of reduced frequency, the inhibiting effect of wakes on boundary layer separation gradually enhances. The separation at the rear part of the suction side is weakened and the separation point moves downstream. However, incoming wakes lead to an increase in dissipation and aerodynamic losses in the main flow area. Excessive reduced frequency ( Fred = 0.74) causes the main flow area to become one of the main source areas of loss. An optimal reduced frequency exists to minimize the aerodynamic loss of the linear cascade.


2005 ◽  
Vol 44 (5) ◽  
pp. 571-590 ◽  
Author(s):  
Ronald Calhoun ◽  
Frank Gouveia ◽  
Joseph Shinn ◽  
Stevens Chan ◽  
Dave Stevens ◽  
...  

Abstract A field program to study atmospheric releases around a complex building was performed in the summers of 1999 and 2000. The focus of this paper is to compare field data with a large-eddy simulation (LES) code to assess the ability of the LES approach to yield additional insight into atmospheric release scenarios. In particular, transient aspects of the velocity and concentration signals are studied. The simulation utilized the finite-element method with a high-fidelity representation of the complex building. Trees were represented with a canopy term in the momentum equation. Inflow and outflow conditions were used. The upwind velocity was constructed from a logarithmic law fitted to velocities obtained on two levels from a tower equipped with a 2D sonic anemometer. A number of different kinds of comparisons of the transient velocity and concentration signals are presented—direct signal versus time, spectral, Reynolds stresses, turbulent kinetic energy signals, and autocorrelations. It is concluded that the LES approach does provide additional insight, but the authors argue that the proper use of LES should include consideration of cost and may require an increased connection to field sensors; that is, higher-resolution boundary and initial conditions need to be provided to realize the full potential of LES.


Author(s):  
Jonathan K. Lai ◽  
Elia Merzari ◽  
Marilyn Delgado ◽  
Samuel J. Lee ◽  
Saya Lee ◽  
...  

The helical coil steam generator (HCSG) is a compact heat exchanger that can have high heat transfer even when the pressure drop is low. This makes it advantageous in small modular reactors and high-temperature reactor designs. In order to investigate the fluid phenomena around these helical banked tubes, a test section was built at Texas A&M University to represent flow across two half-rods within HCSG. This study focuses on the validation of large eddy simulation (LES) for this particular geometry. Pressure tap and particle image velocimetry (PIV) measurements have been recorded at an inlet Reynolds number of 8643, and both mean and fluctuating data is compared with the numerical results. The highly scalable spectral-element code Nek5000 has been used to produce the LES calculations. First, simulations of varying polynomial order expansions are made to determine the spatial resolution required to capture the turbulent scales. Then, simulations with different inlet conditions are compared with experimental data. The pressure drop shows good agreement with pressure tap measurements while velocity shows similar characteristics with PIV. Furthermore, the components of the Reynolds stresses and modes from proper orthogonal decomposition have been developed to validate the physics captured.


2009 ◽  
Vol 643 ◽  
pp. 279-308 ◽  
Author(s):  
D. CHUNG ◽  
D. I. PULLIN

We report direct numerical simulation (DNS) and large-eddy simulation (LES) of statistically stationary buoyancy-driven turbulent mixing of an active scalar. We use an adaptation of the fringe-region technique, which continually supplies the flow with unmixed fluids at two opposite faces of a triply periodic domain in the presence of gravity, effectively maintaining an unstably stratified, but statistically stationary flow. We also develop a new method to solve the governing equations, based on the Helmholtz–Hodge decomposition, that guarantees discrete mass conservation regardless of iteration errors. Whilst some statistics were found to be sensitive to the computational box size, we show, from inner-scaled planar spectra, that the small scales exhibit similarity independent of Reynolds number, density ratio and aspect ratio. We also perform LES of the present flow using the stretched-vortex subgrid-scale (SGS) model. The utility of an SGS scalar flux closure for passive scalars is demonstrated in the present active-scalar, stably stratified flow setting. The multi-scale character of the stretched-vortex SGS model is shown to enable extension of some second-order statistics to subgrid scales. Comparisons with DNS velocity spectra and velocity-density cospectra show that both the resolved-scale and SGS-extended components of the LES spectra accurately capture important features of the DNS spectra, including small-scale anisotropy and the shape of the viscous roll-off.


2019 ◽  
Vol 866 ◽  
pp. 503-525 ◽  
Author(s):  
Racheet Matai ◽  
Paul Durbin

Turbulent flow over a series of increasingly high, two-dimensional bumps is studied by well-resolved large-eddy simulation. The mean flow and Reynolds stresses for the lowest bump are in good agreement with experimental data. The flow encounters a favourable pressure gradient over the windward side of the bump, but does not relaminarize, as is evident from near-wall fluctuations. A patch of high turbulent kinetic energy forms in the lee of the bump and extends into the wake. It originates near the surface, before flow separation, and has a significant influence on flow development. The highest bumps create a small separation bubble, whereas flow over the lowest bump does not separate. The log law is absent over the entire bump, evidencing strong disequilibrium. This dataset was created for data-driven modelling. An optimization method is used to extract fields of variables that are used in turbulence closure models. From this, it is shown how these models fail to correctly predict the behaviour of these variables near to the surface. The discrepancies extend further away from the wall in the adverse pressure gradient and recovery regions than in the favourable pressure gradient region.


1998 ◽  
Vol 120 (2) ◽  
pp. 243-252 ◽  
Author(s):  
S. A. Jordan ◽  
S. A. Ragab

The formation and the downstream transport of the Strouhal vortices in the near wake of a circular cylinder are investigated using the large-eddy simulation (LES) method. The governing equations are formulated in curvilinear coordinates to accommodate a nonorthogonal grid with formal development of a dynamic model to account for the subgrid turbulent scales. Results were produced with and without use of the model. The focus of the investigation is at a subcritical Reynolds number of 5600. Using the dynamic model, the LES results compared best to the published experimental data in terms of both the global and local wake characteristics such as the drag and base pressure coefficients, shedding and detection frequencies, peak vorticity, and the downstream mean velocity-defect and Reynolds stresses. The results further showed streamwise filaments that connect subsequent Strouhal vortices. Qualitatively, the time-averaged Reynolds stresses of the formation region revealed similar symmetric characteristics over the range 525 ≤ Re ≤ 140,000.


Author(s):  
Gregory Rodebaugh ◽  
Zachary Stratton ◽  
Gregory Laskowski ◽  
Michael Benson

Film cooling holes with a compound angle are commonly used on high pressure turbine components in lieu of axial holes to improve effectiveness or as a result of manufacturing constraints. Whereas large eddy simulation (LES) of axial holes is becoming more common place, assessment of LES predictive ability for compound angle hole has been limited. For this study, the selected compound angle round (CAR) hole configuration has a 30 degree injection angle, a 45 degree compound angle, and a density ratio of 1.5. The geometry, flow conditions, and experimental adiabatic effectiveness validation data are from McClintic et al. [28]. The low free stream Mach number of the experiment puts the flow in the incompressible regime. Two LES solvers are evaluated, Fluent and FDL3Di, on structured meshes with a range of blowing ratios simulated for plenum, inline coolant crossflow, and counter coolant crossflow feed holes. When a steady inlet profile is used for the main flow, LES agreement with the data is poor. The inclusion of a resolved turbulent boundary layer significantly improves the predictive quality for both solvers; consequently, resolved inflow turbulence is a required aspect for CAR hole LES. The remaining differences between the simulations and IR data are partly attributed to the steady coolant inlet profiles used for the counter and inline cross feeds.


Sign in / Sign up

Export Citation Format

Share Document