SEARCH · Search NASA

Results for “baroclinic instabilities”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Hydrostatic and Non‐Hydrostatic Baroclinic Instability in the Dynamical Core of the DOE Global Climate Model

Abstract The dynamical core of the Department of Energy global climate model is used to understand the role of non‐hydrostatic dynamics in the simulation of dry and moist baroclinic waves. To reduce computational cost, the Diabatic Acceleration and REscaling approach is adopted. Scale analysis and numerical simulations suggest that the model solution is not distorted by the change of spherical metric terms due to the change of Earth radius, neither are the inter‐scale interactions strongly altered as indicated by the analysis of spectral flux. Compared with hydrostatic simulations, the onset of baroclinic instability is delayed under non‐hydrostatic dynamics as the associated weaker vertical motions tend to increase the critical wavelength and narrow the range of waves that can be baroclinically unstable. During the development of baroclinic waves, non‐hydrostatic dynamics tends to induce vertical motions in the upper troposphere, accelerating the eastward propagation of upper‐level ridges through their impact on local vorticity tendency. These processes reduce the westward tilt of the vertical ridge axes and suppress the conversion of mean flow available potential energy to eddy kinetic energy, leading to weaker baroclinic eddies than in hydrostatic simulations. The contrasts between hydrostatic and non‐hydrostatic settings hold in supplemental experiments that vary the background flow Rossby number, the amount of water vapor content and vertical resolution. We also find that mesoscale and smaller‐scale activities can be considerably under‐represented when the vertical resolution is limited to that typically used in global climate models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Numerical Mixing Suppresses Submesoscale Baroclinic Instabilities Over Sloping Bathymetry

Abstract The impacts of spurious numerical salinity mixing on the larger‐scale flow and tracer fields are characterized using idealized simulations. The idealized model is motivated by realistic simulations of the Texas‐Louisiana shelf and features oscillatory near‐inertial wind forcing. can exceed the physical mixing from the turbulence closure in frontal zones and within the mixed layer. This suggests that simulated mixing processes in frontal zones are driven largely by . Near‐inertial alongshore wind stress amplitude is varied to identify a base case that maximizes the ratio of to in simulations with no prescribed horizontal mixing. We then test the sensitivity of the base case with three tracer advection schemes (MPDATA, U3HC4, and HSIMT) and conduct ensemble runs with perturbed bathymetry. Instability growth is evaluated using the volume‐integrated eddy kinetic energy and available potential energy . While all schemes have similar total mixing, the HSIMT simulations have over double the volume‐integrated and 20% less relative to other schemes, which suppresses the release of and reduces the by roughly 25%. This results in reduced isohaline variability and steeper isopycnals, evidence that enhanced suppresses instability growth. Differences in and between the MPDATA and U3HC4 simulations are marginal. However, the U3HC4 simulations have 25% more . Experiments with variable horizontal viscosity and diffusivity coefficients show that small amounts of prescribed horizontal mixing improve the representation of the ocean state for all advection schemes by reducing the and increasing the .

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Baroclinic Control of Southern Ocean Eddy Upwelling Near Topography

In the Southern Ocean, mesoscale eddies contribute to the upwelling of deep waters along sloping isopycnals, helping to close the upper branch of the meridional overturning circulation. Eddy energy (EE) is not uniformly distributed along the Antarctic Circumpolar Current (ACC). Instead, “hotspots” of EE that are associated with enhanced eddy-induced upwelling exist downstream of topographic features. This study shows that, in idealized eddy-resolved simulations, a topographic feature in the ACC path can enhance and localize eddy-induced upwelling. However, the upwelling systematically occurs in regions where eddies grow through baroclinic instability, rather than in regions where EE is large. Across a range of parameters, along-stream eddy growth rate is a more reliable indicator of eddy upwelling than traditional parameterizations such as eddy kinetic energy, eddy potential energy, or isopycnal slope. Ocean eddy parameterizations should consider metrics specific to the growth of baroclinic instability to accurately model eddy upwelling near topography.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Enhanced Cross-Shelf Exchange by the Eddies Associated With Plume Front

Broadened width of high chlorophyll concentration band with wavy structures, patches, and filaments are often observed along the western coastal next to the Pearl River Estuary over the northern South China Sea shelf during the transition period from winter to summer monsoon. Whereas, there is no such wide band in other seasons. By using a high-resolution numerical model, we reveal that the complex structure and wider band of high coastal chlorophyll concentration results from the smaller scale eddies (about 20–50 km in diameter) associated with buoyant plume-induced salinity front and density fronts, which are roughly along the 30 and 50 m isobaths, respectively. Two trains of eddies are formed along the fronts by the baroclinic instability triggered by frequently alternating winds over the fronts during the period of monsoon transition. The influences of these two trains of eddies are extended in the cross-shelf direction by their interactions, and they can temporally enhance the cross-shelf flow and material exchange. In conclusion, they serve as an efficient pathway to link the inner shelf toward the continental slope.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Global Teleconnections: Atmospheric Interactions Across Scales [Slides]

This project aims at an improved understanding of non-hydrostatic dynamics in atmospheric modeling using the DOE’s next generation climate model. One aspect that we have been focusing on is baroclinic instability which is closely related to mid-latitude cyclone genesis. The other aspect that we recently developed is E3SM simulations under Radiative-Convective Equilibrium (RCE), a prototype of tropical climate.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Growth of Mesoscale Convective Systems in Observations and a Seasonal Convection-Permitting Simulation over Argentina

Abstract A 6.5-month, convection-permitting simulation is conducted over Argentina covering the Remote Sensing of Electrification, Lightning, And Mesoscale/Microscale Processes with Adaptive Ground Observations and Clouds, Aerosols, and Complex Terrain Interactions (RELAMPAGO-CACTI) field campaign and is compared with observations to evaluate mesoscale convective system (MCS) growth prediction. Observed and simulated MCSs are consistently identified, tracked, and separated into growth, mature, and decay stages using top-of-the-atmosphere infrared brightness temperature and surface rainfall. Simulated MCS number, lifetime, seasonal and diurnal cycles, and various cloud-shield characteristics including growth rate are similar to those observed. However, the simulation produces smaller rainfall areas, greater proportions of heavy rainfall, and faster system propagations. Rainfall area is significantly underestimated for long-lived MCSs but not for shorter-lived MCSs, and rain rates are always overestimated. These differences result from a combination of model and satellite retrieval biases, in which simulated MCS rain rates are shifted from light to heavy, while satellite-retrieved rainfall is too frequent relative to rain gauge estimates. However, the simulation reproduces satellite-retrieved MCS cloud-shield evolution well, supporting its usage to examine environmental controls on MCS growth. MCS initiation locations are associated with removal of convective inhibition more than maximized low-level moisture convergence or instability. Rapid growth is associated with a stronger upper-level jet (ULJ) and a deeper northwestern Argentinean low that causes a stronger northerly low-level jet (LLJ), increasing heat and moisture fluxes, low-level vertical wind shear, baroclinicity, and instability. Sustained growth corresponds to similar LLJ, baroclinicity, and instability conditions but is less sensitive to the ULJ, large-scale vertical motion, or low-level shear. Growth sustenance controls MCS maximum extent more than growth rate.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

The effect of initial conditions on mixing transition of the Richtmyer–Meshkov instability

We investigate the late-time Richtmyer–Meshkov instability (RMI) growth of sinuous perturbations on an air/sulphur hexafluoride interface (Atwood number, A ~ 0.67 ) subjected to a Mach 1.2 planar shock wave at Los Alamos National Laboratory's vertical shock tube facility. Interface perturbations are established using a novel membraneless technique where cross-flowing air and SF 6 separated by an oscillating splitter plate create a perturbed density interface. The interface formed has multi-modal features and residual small perturbations, however, a dominant mode is still noticeable. The late-time perturbation growths scale with $ka_0$ initial conditions (where $k$ is the wavenumber and $a_0$ is the initial amplitude of the dominant mode) as measured at the pre-shock interface. Past nonlinear models based on potential-flow theory, heuristic/interpolation approaches, Padé approximants and numerical simulations are evaluated against present experimental results. Accounting for an explicit $ka_0$ dependence in Sadot et al.'s model, we propose an empirical rational function that captures the asymptotic behaviour of perturbation growth for a broad range of initial conditions ( 0.30 $≤ ka_0 ≤$ 0.86 ). The onset of mixing transition and its initial condition dependence are investigated with respect to the minimum state criterion ($Re =$ 1.6 $\times$ 10 5 ) for unsteady flows by Zhou. Earlier mixing transitions for higher ka0 initial conditions are noted from local and global Reynolds number estimates which are corroborated by the existence of an inertial sub-range and formation of mixing regions indicating the physical significance of the minimum state criterion in RMI flows. The transition is accompanied by the increasing teapot-like appearance of joint probability density functions of $p–q$ (invariants of the reduced velocity gradient tensor), establishing the technique as a useful tool for turbulence detection in two-dimensional diagnostics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Floquet stability analyses of stratified oscillating boundary layers on adiabatic slopes

The presence of a no-slip, impermeable, adiabatic, sloped boundary in an otherwise quiescent, stably stratified, Boussinesq flow generates baroclinic vorticity within a diffusive boundary layer. Such conditions are typical of the oscillating boundary layers on adiabatic abyssal slopes, sloped lake bathymetry and sloped coastal bathymetry in the absence of high-wavenumber internal waves, mean flows, far-field turbulence on larger scales, and resonant tidal–bathymetric interaction. We investigate the linear stability of the oscillating flow within non-dimensional parameter space typical of the $M_2$ tide and hydraulically smooth, middle-latitude abyssal slopes through Floquet linear stability analysis. The flow dynamics depends on three non-dimensional variables: the Reynolds number for Stokes’ second problem, the Prandtl number, and a frequency ratio that accounts for the resonance conditions ($C$, criticality) of the buoyant restoring force and the tidal forcing. The Floquet analysis results suggest that oscillating laminar boundary layers on adiabatic abyssal slopes are increasingly unstable as Reynolds number, criticality parameter and/or spanwise disturbance wavenumber are increased. We also show that the two-dimensional Floquet linear instability necessarily generates three-dimensional baroclinic vorticity, which suggests that the evolution of the gravitational instabilities may be nonlinear as $t→∞$.

42 ENGINEERING↗

DOE OSTI2024

Radiation and heat transport in divergent shock–bubble interactions

Shock–bubble interactions (SBIs) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal resolution, making simulations a vital tool in understanding these systems. In this study, we benchmark several radiation and thermal transport models in the xRAGE hydrodynamic code against experimental images of a collapsing mesoscale void during the passage of a 300 GPa shock. We also quantitatively examine the role of transport physics in the evolution of the SBI. This allows us to understand the dynamics of the interaction at timescales shorter than experimental imaging framerates. We find that all radiation models examined reproduce empirical shock velocities within experimental error. Radiation transport is found to reduce shock pressures by providing an additional energy pathway in the ablation region, but this effect is small (∼1% of total shock pressure). Employing a flux-limited Spitzer model for heat conduction, we find that flux limiters between 0.03 and 0.10 produce agreement with experimental velocities, suggesting that the system is well-within the Spitzer regime. Higher heat conduction is found to lower temperatures in the ablated plasma and to prevent secondary shocks at the ablation front, resulting in weaker primary shocks. Finally, we confirm that the SBI-driven instabilities observed in the HED regime are baroclinically driven, as in the low energy case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

High pressure/high temperature multiphase simulations of dodecane injection to nitrogen: Application on ECN Spray-A

The present work investigates the complex phenomena associated with pressure/high temperature dodecane injection for the Engine Combustion Network (ECN) Spray-A case, employing more elaborate thermodynamic closures, to avoid well known deficiencies concerning density and speed of sound prediction using traditional cubic models. A tabulated thermodynamic approach is proposed here, based on log 10 (p)-T tables, providing very high accuracy across a large range of pressures, spanning from 0 to 2500 bar, with only a small number of interpolation points. The tabulation approach is directly extensible to any thermodynamic model, existing or to be developed in the future. Here NIST REFPROP properties are used, combined with PC-SAFT Vapor-Liquid-Equilibrium to identify the liquid in mixtures penetration, hence avoiding the use of an arbitrary threshold for mass fraction. Identified liquid and vapour penetration are compared against experimental data from the ECN database showing a good agreement, within approximately 3–8% for axial penetration of liquid, 2% for vapor axial penetration and within experimental uncertainty for radial distribution of mass fraction. Analysis of the vortex evolution indicates that driving mechanisms behind the jet break-up are vortex tilting/stretching, then baroclinic torque, leading to Rayleigh-Taylor instabilities, closely followed by vortex dilation and finally viscous effects.

42 ENGINEERING↗

DOE OSTI2020

Dependence of Enstrophy Transport and Mixed Mass on Dimensionality and Initial Conditions in the Richtmyer–Meshkov Instability Induced Flows

Here we present a comparative study of the enstrophy budget and mixed mass between two- and three-dimensional flows induced by Richtmyer–Meshkov instability (RMI). Specifically, the individual contributions to the enstrophy budget due to the production from baroclinicity and from vortex stretching (which vanishes in two-dimensional (2D) flow) are delineated. This is enabled by a set of two- and three-dimensional computations at Atwood 0.5 having both narrow- and broad-band perturbations. A further three-dimensional (3D) computation is conducted at Atwood 0.9 using an identical narrowband perturbation to the Atwood 0.5 case to examine the sensitivity to density ratio. The mixed mass is also considered with the goal to obtain insight on how faithfully a simplified calculation performed in two dimensions can capture the mixed mass for an inertial confinement fusion (ICF) or other practical application. It is shown that the late time power law decay of variable density enstrophy is substantially different in two and three dimensions for the narrowband initial perturbation. The baroclinic production term is negligible in three dimensions (aside from the initial shock interaction), as vortex stretching is larger by two orders of magnitude. The lack of vortex stretching considerably reduces the decay rate in both narrowband and broadband perturbations in two dimensions. In terms of mixed mass, the lack of vortex stretching reduces the mixed mass in two dimensions compared to three in all cases. In the broadband cases, the spectral bandwidth in the 2D case is wider; hence, there is a longer time period of sustained linear growth which reduces the normalized mixed mass further.

42 ENGINEERING↗

DOE OSTI2021

Vortex-sheet modeling of hydrodynamic instabilities produced by an oblique shock interacting with a perturbed interface in the HED regime

In this work, we consider hydrodynamic instabilities produced by the interaction of an oblique shock with a perturbed material interface under high-energy-density (HED) conditions. During this interaction, a baroclinic torque is generated along the interface due to the misalignment between the density and pressure gradients, thus leading to perturbation growth. Our objective is to understand the competition between the impulsive acceleration due to the normal component of the shock velocity, which drives the Richtmyer–Meshkov instability, and the shear flow across the interface due to the tangential component of the shock velocity, which drives the Kelvin–Helmholtz instability, as well as its relation to perturbation growth. Since the vorticity resulting from the shock-interface interaction is confined to the interface, we describe the perturbation growth using a two-dimensional vortex-sheet model. We demonstrate the ability of the vortex-sheet model to reproduce roll-up dynamics for non-zero Atwood numbers by comparing to past laser-driven HED experiments. We determine the dependence of the interface dynamics on the tilt angle and propose a time scaling for the perturbation growth at early time. Eventually, this scaling will serve as a platform for the design of future experiments. This study is the first attempt to incorporate into a vortex-sheet model the time-dependent interface decompression and the deceleration (as well as the corresponding Rayleigh–Taylor instability) arising from laser turn-off.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

The influence of surface heat fluxes on the growth of idealized monsoon depressions

This study explores the effect of surface sensible and latent heat fluxes on monsoon depressions using a series of idealized convection-permitting simulations. Each experiment is initialized with a small-amplitude wave that is allowed to grow within an environment representative of the South Asian monsoon. Comparing experiments with and without interactive surface heat fluxes, it is found that these fluxes enhance the growth of the simulated vortices. Without interactive surface fluxes, the strengthening period is short and the vortices fail to reach intensities characteristic of stronger monsoon depressions. Using a large set of experiments in which the vertical and meridional shear are systematically varied, it is found that surface heat fluxes enhance intensity the most when the upper-level shear is weak, the lower-level shear and associated moist static energy (MSE) gradient are sufficiently steep, and the lower-level meridional shear is strong. Here, these experiments reveal two different regimes of convection-coupled monsoon depression growth: one in which convection is driven by MSE advection and one in which it is driven by surface heat fluxes and quasigeostrophic forcing for ascent. Both regimes require sufficiently strong meridional shear to achieve initial growth by barotropic instability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Shock-driven dispersal of a corrugated finite-thickness particle layer

A research area emerging in the multiphase flow community is the study of shock-driven multiphase instability (SDMI), a gas–particle analog of the traditional fluid-fluid Richtmyer–Meshkov instability (RMI). In this work, we study the interaction of planar air shocks with corrugated glass particle curtains through the use of numerical simulations with an Eulerian–Lagrangian approach. One objective of this study is to compare the simulated particle curtains to a comparable set of shock tube experiments performed to analyze traditional RMI of a gas curtain. The simulations are set to match the experimental shock Mach numbers and perturbation wavelengths (3.6 and 7.2 mm) while also matching the Atwood number of the experiments to the multiphase Atwood number of the simulations. Varying particle diameters are tested in the simulations to explore the impact of particle diameter on the evolution of the particle curtain. This simulation setup allows for a one-to-one comparison between RMI and SDMI under comparable conditions while also allowing for a separate study into the validity of the use of the multiphase Atwood number to compare the single-phase and multiphase instabilities. In particular, we show that the comparison depends on the diameter of the particles (thus, dependent on the Stokes number of the flow). A second objective of this study is to analyze the effect of the initial particle volume fraction on the evolution of the curtain and the behavior of the instability. This is done through analyzing the effect of the multiphase terms of the vorticity evolution equation on the vorticity deposition in SDMI. Also discussed is the effect of the particle diameter on the multiphase generation terms as well as in the baroclinic vorticity generation term in SDMI as the shock passes over the curtain.

42 ENGINEERING↗

DOE OSTI2025

An Observational Evaluation of RKW Theory over the U.S. Southern Great Plains

The theory of Rotunno et al. (“RKW” theory) addresses the behavior of squall-line cold pools in vertically sheared flows. It predicts that, within a given thermodynamic environment, a balance between baroclinic vorticity generation by the cold pool and low-level environmental vertical wind shear induces an upright updraft along the gust front that maximizes the initiation of new convective cells. Although this theory has been evaluated numerically, its applicability to observed systems remains unclear and is limited by a lack of critical measurements, including high-frequency thermodynamic and wind profiles across the gust front. Herein, observations from the Atmospheric Radiation Measurement Southern Great Plains (ARM-SGP) observatory near Lamont, Oklahoma, are used to evaluate RKW theory for 10 well-observed squall lines over a 11-yr period. For this evaluation, RKW parameters including cold-pool intensity (c), low-level ambient, line-normal vertical shear (ΔV n ), subcloud and cloud-layer updraft tilts, and multiple measures of system intensity are estimated. Furthermore, the c estimates rely on thermodynamic retrievals from the Atmosphere Emitted Radiance Interferometer (AERI), which are uncertain but verify reasonably well against independent observations. As predicted by the theory, for c/ΔV n ≥ 1, c/ΔV n correlates positively with updraft tilt and negatively with system intensity, but these results are not always statistically significant and are also sensitive to the method by which ΔV n is evaluated. Specifically, ΔV n evaluations that extend above the cold-pool top yield greater consistency with RKW predictions. Also, some measures of intensity correlate more strongly with standard moist instability metrics than with RKW parameters.

Cold pools↗

DOE OSTI2024

Analysis of single-mode Richtmyer–Meshkov instability using high-order incompressible vorticity-streamfunction and shock-capturing simulations

In this study, two- and three-dimensional simulation results obtained using a new high-order incompressible, variable-density vorticity–streamfunction (VS) method and data from previous ninth-order weighted essentially nonoscillatory (WENO) shock-capturing simulations are used to investigate the nonlinear dynamics of single-mode Richtmyer–Meshkov instability using a model of a Mach 1.3 air(acetone)/SF6 shock tube experiment. A comparison of the density fields from both simulations with the experimental images demonstrates very good agreement in the large-scale structure with both methods but differences in the small-scale structure. The WENO method captures the small-scale disordered structure observed in the experiment, while the VS method partially captures such structure and yields a strong rotating core. The perturbation amplitude growth from the simulations generally agrees well with the experiment. The simulation bubble and spike amplitudes agree well at early times. At later times, the WENO bubble amplitude is smaller than the VS amplitude and vice versa for the spike amplitude. The predictions of nonlinear single-mode instability growth models are shown to agree with the simulation amplitudes at early-to-intermediate times but underpredict the amplitudes at later times in the nonlinear regime. Visualizations of the mass fraction and enstrophy isosurfaces, velocity and vorticity fields, and baroclinic vorticity production and vortex stretching terms from the three-dimensional simulations indicate that, with the exception of the small-scale structure within the rollups, the VS and WENO results are in good agreement.

42 ENGINEERING↗

DOE OSTI2022

Three-dimensional simulations of reshocked inclined Richtmyer-Meshkov instability: Effects of initial perturbations

The effect of initial perturbations on the evolution of the inclined Richtmyer-Meshkov turbulent mixing layer before and after reshock initiated by a shock wave with Mach number 1.55 is investigated through three-dimensional (3D) simulations using the flash code. The 3D simulations aim to reproduce both predominantly single-mode and multimode interfaces between light and heavy gases (N 2 -CO 2 , Atwood number, A≈0.22; amplitude to wavelength ratio of 0.088) which were created in an inclined shock tube facility to analyze the effects of initial conditions on mixing development in the entire flow field. The two-dimensional center slices of 3D simulations are compared with the experimental results to validate the computational code. Mixing width, mixed mass, mixed-mass thickness, and circulation in addition to concentration fields are shown to be in good agreement with the experimental data. The three-dimensional density and vorticity fields are first presented to qualitatively describe the flow behavior before and after reshock. Several measured density/velocity-related quantities indicate that the growth of the mixing material is strongly dependent on initial conditions. Before reshock and at early times after reshock, flow is clearly maintaining the memory of initial perturbations. However, at late time after reshock, although the large wavelength feature still dominates the flow motion, and the morphology of the two different interfaces indicates several differences, by breakdown of large-scale coherent structures to much finer scales, the memory of small scales of the multimode initial perturbation is not as clear as pre-reshock. Regarding three-dimensionality of the flow, before reshock in the multimode case, the baroclinic vorticity production, circulation, turbulent kinetic energy, and turbulent mass flux suggest that the small-scale roll-up features along the large inclined wavelength quickly evolves in all three dimensions. The coherent vortex tubes break down to smaller wormlike vortex structures, and turbulent fluctuations in the out-of-plane dimension are comparable to the spanwise direction. After reshock, this three-dimensionality of mixing growth was observed in the flow for both initial conditions. The results of this work represent a significant extension of previous computational studies performed on this specific topic. A different code with a different numerical method is validated through comparison with the experimental data. The initial perturbations are directly measured from the experimental results. Moreover, the entire three-dimensional experimental shock tube domain is simulated, and more quantities are investigated to understand the mixing mechanism and instability evolution in all three dimensions.

42 ENGINEERING↗

DOE OSTI2020

Effect of the numerical discretization scheme in Shock-Driven turbulent mixing simulations

In this work, we evaluate the effects of distinct numerical strategies in simulations of shock-driven turbulent mixing. An air-SF 6 -air gas curtain subjected to Mach 1.2 shock-waves is computed with Implicit Large-Eddy Simulation based on three numerical schemes: directional-split with Harten-Lax-van Leer (HLL) solver; directional-unsplit using HLL-Contact (HLLC) solver; and directional-unsplit with HLLC solver and a Low Mach number Correction (LMC). The results illustrate the importance of the numerical strategy to the accuracy of the predictions. Whereas both split and unsplit schemes result in a similar spatial development of the initial shock-driven instability, only the unsplit schemes can predict the turbulent mixing transition after reshock observed by the laboratory experiments. Such feature increases the mixing rate of the two fluids, this being particularly pronounced when the LMC is active due to i) the reduced flow Mach number, and ii) the larger effective Reynolds number. Since the selected mixing problem is driven by the deposition of vorticity at the fluids’ interface, the resultant flow physics is analyzed by investigating the contribution of distinct inviscid mechanisms to the production of vorticity: baroclinicity, stretching, and dilatation. As expected, it is observed that the production of vorticity is initially dominated by the baroclinicity mechanism. Yet, the relevance of the remainder mechanisms is enhanced after the reshock and may even surpass the baroclinicity term.

42 ENGINEERING↗