Reynolds number dependence of turbulence induced by the Richtmyer-Meshkov instability using direct numerical simulations

This paper investigates the Reynolds number dependence of a turbulent mixing layer evolving from the Richtmyer-Meshkov instability using a series of direct numerical simulations of a well-defined narrowband initial condition for a range of different Reynolds numbers. The growth rate exponent of the integral width and mixed mass is shown to marginally depend on the initial Reynolds number Re0, as does the minimum value of the molecular mixing fraction. The decay rates of turbulent kinetic energy and its dissipation rate are shown to decrease with increasing Re0, while the spatial distribution of these quantities is biased towards the spike side of the layer. The normalised dissipation rate and scalar dissipation rate are calculated and are observed to be approaching a high Reynolds number limit. By fitting an appropriate functional form, the asymptotic value of these two quantities is estimated as 1.54 and 0.66. Finally, an evaluation of the mixing transition criterion for unsteady flows is performed, showing that even for the highest Re0 case the turbulence in the flow is not yet fully developed. This is despite the observation of a narrow inertial range in the turbulent kinetic energy spectra, with a scaling close to -3/2

Atwood ratio dependence of Richtmyer-Meshkov flows under reshock conditions using large-eddy simulations

We study the shock-driven turbulent mixing that occurs when a perturbed planar density interface is impacted by a planar shock wave of moderate strength and subsequently reshocked. The present work is a systematic study of the influence of the relative molecular weights of the gases in the form of the initial Atwood ratio A. We investigate the cases A = ± 0.21, ±0.67 and ±0.87 that correspond to the realistic gas combinations air–CO_2, air–SF_6 and H_2–air. A canonical, three-dimensional numerical experiment, using the large-eddy simulation technique with an explicit subgrid model, reproduces the interaction within a shock tube with an endwall where the incident shock Mach number is ~1.5 and the initial interface perturbation has a fixed dominant wavelength and a fixed amplitude-to-wavelength ratio ~0.1. For positive Atwood configurations, the reshock is followed by secondary waves in the form of alternate expansion and compression waves travelling between the endwall and the mixing zone. These reverberations are shown to intensify turbulent kinetic energy and dissipation across the mixing zone. In contrast, negative Atwood number configurations produce multiple secondary reshocks following the primary reshock, and their effect on the mixing region is less pronounced. As the magnitude of A is increased, the mixing zone tends to evolve less symmetrically. The mixing zone growth rate following the primary reshock approaches a linear evolution prior to the secondary wave interactions. When considering the full range of examined Atwood numbers, measurements of this growth rate do not agree well with predictions of existing analytic reshock models such as the model by Mikaelian (Physica D, vol. 36, 1989, p. 343). Accordingly, we propose an empirical formula and also a semi-analytical, impulsive model based on a diffuse-interface approach to describe the A-dependence of the post-reshock growth rate

Towards a solution of the closure problem for convective atmospheric boundary-layer turbulence

We consider the closure problem for turbulence in the dry convective atmospheric boundary layer (CBL). Transport in the CBL is carried by small scale eddies near the surface and large plumes in the well mixed middle part up to the inversion that separates the CBL from the stably stratified air above. An analytically tractable model based on a multivariate Delta-PDF approach is developed. It is an extension of the model of Gryanik and Hartmann [1] (GH02) that additionally includes a term for background turbulence. Thus an exact solution is derived and all higher order moments (HOMs) are explained by second order moments, correlation coefficients and the skewness. The solution provides a proof of the extended universality hypothesis of GH02 which is the refinement of the Millionshchikov hypothesis (quasi- normality of FOM). This refined hypothesis states that CBL turbulence can be considered as result of a linear interpolation between the Gaussian and the very skewed turbulence regimes. Although the extended universality hypothesis was confirmed by results of field measurements, LES and DNS simulations (see e.g. [2-4]), several questions remained unexplained. These are now answered by the new model including the reasons of the universality of the functional form of the HOMs, the significant scatter of the values of the coefficients and the source of the magic of the linear interpolation. Finally, the closures 61 predicted by the model are tested against measurements and LES data. Some of the other issues of CBL turbulence, e.g. familiar kurtosis-skewness relationships and relation of area coverage parameters of plumes (so called filling factors) with HOM will be discussed also

Institute for Computational Mechanics in Propulsion (ICOMP)

The Institute for Computational Mechanics in Propulsion (ICOMP) is a combined activity of Case Western Reserve University, Ohio Aerospace Institute (OAI) and NASA Lewis. The purpose of ICOMP is to develop techniques to improve problem solving capabilities in all aspects of computational mechanics related to propulsion. The activities at ICOMP during 1991 are described

The 1999 Center for Simulation of Dynamic Response in Materials Annual Technical Report

Introduction: This annual report describes research accomplishments for FY 99 of the Center for Simulation of Dynamic Response of Materials. The Center is constructing a virtual shock physics facility in which the full three dimensional response of a variety of target materials can be computed for a wide range of compressive, ten- sional, and shear loadings, including those produced by detonation of energetic materials. The goals are to facilitate computation of a variety of experiments in which strong shock and detonation waves are made to impinge on targets consisting of various combinations of materials, compute the subsequent dy- namic response of the target materials, and validate these computations against experimental data

Pressure Gain Combustion: Fuel Spray and Shockwave Interaction

Pressure gain combustion can attain higher thermodynamic cycle efficiency in gas turbine power systems, resulting in the reduction of specific fuel consumption/fuel burn and Carbon dioxide emissions.There are many ways to achieve pressure gain and the present research investigates pressure gain through shock bubble (gas and liquid bubble) interaction (SBI) using computational fluid dynamics (CFD) simulations. The numerical simulations have been performed in 2D and 3D representations of the shock tube to depict the interaction of a planar shock wave with distinct gas and liquid inhomogeneities. The three scenarios considered cover the interaction of a planar shock wave in air with: spherical helium bubble (Mach number, Ma = 1.25); cylindrical helium bubble (Ma = 1.22) and cylindrical water bubble (Ma = 1.47). To perform these simulations, the Unsteady Reynolds-Averaged Navier-Stokes (URANS) mathematical model and the coupled level set and VOF method within the commercial CFD code, ANSYS FLUENT, have been applied. A finite volume method (FVM) is also employed to solve the governing equations. For the spherical and cylindrical gas bubble cases, various quantitative analyses are presented and compared to the experimental work of Haas and Sturtevant (1987). These include: refracted wave, transmitted wave, upstream interface, downstream interface, jet, vortex filament, non-dimensional bubble, and vortex velocities. The predicted non-dimensional bubble and vortex velocities have also been compared with experimental data, a simple model of shock- induced Rayleigh Taylor (RT) instability and other theoretical models. Comparisons are also shown between the predicted bubble length/width and the experimentally measured results to elucidate changes in the shape and size of the 2D and 3D bubbles. Additional quantitative analyses are also presented for the spherical bubble involving the size estimation of the vortex pair as well as their spacing. For the shock cylindrical water bubble interaction case, the quantitative predictions include: displacement/drift, acceleration, distortion in the lateral direction, distortion in flow direction, area variation from bubble distortion, as well as drag coefficient and are compared to the experimental measurements of Igra et al. (2002). It has been demonstrated that 3D simulations compare very well with the experimental data, suggesting that 3D simulations are necessary to capture SBI process accurately. Finally, comprehensive flow visualization has been used to elucidate the shock-bubble interaction (SBI) process from bubble compression to the formation of the vortex filaments (cylindrical helium bubble), vortex rings (spherical helium bubble), vortices (cylindrical water bubble) as well as the production and distribution of vorticity. It is demonstrated for the first time that turbulence is generated at the early phase of the SBI process, with the maximum turbulence intensity reaching about 20% around the vortex filaments/vortex rings regions for the cylindrical/spherical helium bubble cases respectively and about 22% for the cylindrical water bubble case at the later phase of the interaction process

Simulations of Turbulent Flows with Strong Shocks and Density Variations: Final Report

The target of this SciDAC Science Application was to develop a new capability based on high-order and high-resolution schemes to simulate shock-turbulence interactions and multi-material mixing in planar and spherical geometries, and to study Rayleigh-Taylor and Richtmyer-Meshkov turbulent mixing. These fundamental problems have direct application in high-speed engineering flows, such as inertial confinement fusion (ICF) capsule implosions and scramjet combustion, and also in the natural occurrence of supernovae explosions. Another component of this project was the development of subgrid-scale (SGS) models for large-eddy simulations of flows involving shock-turbulence interaction and multi-material mixing, that were to be validated with the DNS databases generated during the program. The numerical codes developed are designed for massively-parallel computer architectures, ensuring good scaling performance. Their algorithms were validated by means of a sequence of benchmark problems. The original multi-stage plan for this five-year project included the following milestones: 1) refinement of numerical algorithms for application to the shock-turbulence interaction problem and multi-material mixing (years 1-2); 2) direct numerical simulations (DNS) of canonical shock-turbulence interaction (years 2-3), targeted at improving our understanding of the physics behind the combined two phenomena and also at guiding the development of SGS models; 3) large-eddy simulations (LES) of shock-turbulence interaction (years 3-5), improving SGS models based on the DNS obtained in the previous phase; 4) DNS of planar/spherical RM multi-material mixing (years 3-5), also with the two-fold objective of gaining insight into the relevant physics of this instability and aiding in devising new modeling strategies for multi-material mixing; 5) LES of planar/spherical RM mixing (years 4-5), integrating the improved SGS and multi-material models developed in stages 3 and 5. This final report is outlined as follows. Section 2 shows an assessment of numerical algorithms that are best suited for the numerical simulation of compressible flows involving turbulence and shock phenomena. Sections 3 and 4 deal with the canonical shock-turbulence interaction problem, from the DNS and LES perspectives, respectively. Section 5 considers the shock-turbulence inter-action in spherical geometry, in particular, the interaction of a converging shock with isotropic turbulence as well as the problem of the blast wave. Section 6 describes the study of shock-accelerated mixing through planar and spherical Richtmyer-Meshkov mixing as well as the shock-curtain interaction problem In section 7 we acknowledge the different interactions between Stanford and other institutions participating in this SciDAC project, as well as several external collaborations made possible through it. Section 8 presents a list of publications and presentations that have been generated during the course of this SciDAC project. Finally, section 9 concludes this report with the list of personnel at Stanford University funded by this SciDAC project

Large Eddy Simulations in Astrophysics

In this review, the methodology of large eddy simulations (LES) is introduced and applications in astrophysics are discussed. As theoretical framework, the scale decomposition of the dynamical equations for neutral fluids by means of spatial filtering is explained. For cosmological applications, the filtered equations in comoving coordinates are also presented. To obtain a closed set of equations that can be evolved in LES, several subgrid scale models for the interactions between numerically resolved and unresolved scales are discussed, in particular the subgrid scale turbulence energy equation model. It is then shown how model coefficients can be calculated, either by dynamical procedures or, a priori, from high-resolution data. For astrophysical applications, adaptive mesh refinement is often indispensable. It is shown that the subgrid scale turbulence energy model allows for a particularly elegant and physically well motivated way of preserving momentum and energy conservation in AMR simulations. Moreover, the notion of shear-improved models for inhomogeneous and non-stationary turbulence is introduced. Finally, applications of LES to turbulent combustion in thermonuclear supernovae, star formation and feedback in galaxies, and cosmological structure formation are reviewed.Comment: 64 pages, 23 figures, submitted to Living Reviews in Computational Astrophysic