Search NASA⌕ Search

SEARCH · Search NASA

Results for “convective instability”

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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Shock Ignition Laser-Plasma Interactions in Ignition-Scale Plasmas

Here, we use a subignition scale laser, the 30 kJ Omega, and a novel shallow-cone target to study laser-plasma interactions at the ablation-plasma density scale lengths and laser intensities anticipated for direct drive shock-ignition implosions at National Ignition Facility scale. Our results show that, under these conditions, the dominant instability is convective stimulated Raman scatter with experimental evidence of two plasmon decay (TPD) only when the density scale length is reduced. Particle-in-cell simulations indicate this is due to TPD being shifted to lower densities, removing the experimental back-scatter signature and reducing the hot-electron temperature. The experimental laser energy-coupling to hot electrons was found to be 1%–2.5%, with electron temperatures between 35 and 45 keV. Radiation-hydrodynamics simulations employing these hot-electron characteristics indicate that they should not preheat the fuel in MJ-scale shock ignition experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Effects of Atwood and Reynolds numbers on the evolution of buoyancy-driven homogeneous variable-density turbulence

The evolution of buoyancy-driven homogeneous variable-density turbulence (HVDT) at Atwood numbers up to 0.75 and large Reynolds numbers is studied by using high-resolution direct numerical simulations. To help understand the highly non-equilibrium nature of buoyancy-driven HVDT, the flow evolution is divided into four different regimes based on the behaviour of turbulent kinetic energy derivatives. The results show that each regime has a unique type of dependence on both Atwood and Reynolds numbers. It is found that the local statistics of the flow based on the flow composition are more sensitive to Atwood and Reynolds numbers compared to those based on the entire flow. It is also observed that, at higher Atwood numbers, different flow features reach their asymptotic Reynolds-number behaviour at different times. The energy spectrum defined based on the Favre fluctuations momentum has less large-scale contamination from viscous effects for variable-density flows with constant properties, compared to other forms used previously. The evolution of the energy spectrum highlights distinct dynamical features of the four flow regimes. Thus, the slope of the energy spectrum at intermediate to large scales evolves from -7/3 to -1 , as a function of the production-to-dissipation ratio. The classical Kolmogorov spectrum emerges at intermediate to high scales at the highest Reynolds numbers examined, after the turbulence starts to decay. Finally, the similarities and differences between buoyancy-driven HVDT and the more conventional stationary turbulence are discussed and new strategies and tools for analysis are proposed.

42 ENGINEERING↗
DOE OSTI2023

An investigation of non-canonical mixing in red giant stars using APOGEE 12 C/ 13 C ratios observed in open cluster stars

Standard stellar evolution theory poorly predicts the surface abundances of chemical species in low-mass, red giant branch (RGB) stars. Observations show an enhancement of p–p chain and CNO cycle products in red giant envelopes, which suggests the existence of non-canonical mixing that brings interior burning products to the surface of these stars. The 12 C/ 13 C ratio is a highly sensitive abundance metric used to probe this mixing. We investigate extra RGB mixing by examining: (1) how 12 C/ 13 C is altered along the RGB, and (2) how 12 C/ 13 C changes for stars of varying age and mass. Our sample consists of 43 red giants, spread over 15 open clusters from the Sloan Digital Sky Survey’s APOGEE DR17, that have reliable 12 C/ 13 C ratios derived from their APOGEE spectra. We vetted these 12 C/ 13 C ratios and compared them as a function of evolution and age/mass to the standard mixing model of stellar evolution, and to a model that includes prescriptions for RGB thermohaline mixing and stellar rotation. We find that the observations deviate from standard mixing models, implying the need for extra mixing. Additionally, some of the abundance patterns depart from the thermohaline model, and it is unclear whether these differences are due to incomplete observations, issues inherent to the model, our assumption of the cause of extra mixing, or any combination of these factors. Nevertheless, the surface abundances across our age/mass range clearly deviate from the standard model, agreeing with the notion of a universal mechanism for RGB extra mixing in low-mass stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dependencies of Simulated Convective Cell and System Growth Biases on Atmospheric Instability and Model Resolution

Abstract This study evaluates convective cell properties and their relationships with convective and stratiform rainfall within a season‐long convection‐permitting weather research and forecasting simulation over central Argentina using radar, satellite, and radiosonde measurements from the RELAMPAGO‐CACTI field campaign. The simulation slightly underestimates radar‐estimated rainfall over the ∼3.5‐month evaluation period but underestimates stratiform rainfall by 46% and overestimates convective rainfall by 43%. As convective available potential energy (CAPE) increases, the convective rainfall overestimation decreases, but the stratiform rainfall underestimation increases such that the contribution of convective to total rainfall remains constantly high biased by ∼26%. Overestimated convective rainfall arises from the simulation generating 2.6 times more precipitating convective cells (14,299) than observed by radar (5,662) despite similar observed and simulated cell growth processes, with relatively wide cells contributing mostly to excessive convective rainfall. Relatively shallow cells, typically reaching heights of 4–7 km, contribute most to the cell number bias. This cell number bias increases as CAPE decreases, potentially because cells and their updrafts become narrower and more under‐resolved as CAPE decreases. The gross overproduction of precipitating shallow cells leads to overly efficient precipitation and inadequate detrainment of ice aloft, thereby diminishing the formation of robust stratiform rainfall regions. Decreasing model horizontal grid spacing from 3 to 1 or 0.333 km for low (<300 J kg −1 ) and high CAPE (>1,000 J kg −1 ) cases results in minimal change to cell number, depth, and convective‐to‐stratiform partitioning biases. This suggests that improving prediction of these convective properties depends on factors beyond solely increasing model resolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Convective-storm environments in subtropical South America from high-frequency soundings during RELAMPAGO-CACTI

During the RELAMPAGO-CACTI field experiments in 2018-19, an unprecedented number of balloon-borne soundings were collected in Argentina. Radiosondes were launched from both fixed and mobile platforms, yielding 2712 soundings during the period 15 October 2018-30 April 2019. Approximately 20% of these soundings were collected by highly mobile platforms, strategically positioned for each intensive observing period, and launching approximately once per hour. The combination of fixed and mobile soundings capture both the overall conditions characterizing the RELAMPAGO-CACTI campaign, as well as the detailed evolution of environments supporting the initiation and upscale growth of deep convective storms, including some that produced hazardous hail and heavy rainfall. Episodes of frequent convection were characterized by sufficient quantities of moisture and instability for deep convection, along with deep-layer vertical wind shear supportive of organized or rotating storms. Eleven soundings showed most-unstable convective available potential energy (MUCAPE) exceeding 6000 J kg-1, comparable to the extreme instability observed in other parts of the world with intense deep convection. Parameters used to diagnose severe-storm potential showed that conditions were often favorable for supercells and severe hail, but not for tornadoes, primarily owing to insufficient low-level wind shear. High-frequency soundings also revealed the structure and evolution of the boundary layer leading up to convection initiation, convectively generated cold pools, the South American Low-Level Jet (SALLJ), and elevated nocturnal convection. This sounding dataset will enable improved understanding and prediction of convective storms and their surroundings in subtropical South America, as well as comparisons with other heavily studied regions such as the central United States that have not previously been possible.

Schumacher, Russ S.↗
DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

On the Dependence of Simulated Convection on Domain Size in CRMs

Abstract We present a heuristic model to explain the suppression of deep convection in convection‐resolving models (CRMs) with a small number of grid columns, such as those used in super‐parameterized or multi‐scale modeling framework (MMF) general circulation models (GCM) of the atmosphere. Domains with few grid columns require greater instability to sustain convection because they force a large convective fraction, driving strong compensating subsidence warming. Updraft dilution, which is stronger for reduced horizontal grid spacing, enhances this effect. Thus, suppression of deep convection in CRMs with few grid columns can be reduced by increasing grid spacing. Radiative‐convective equilibrium simulations using standalone CRM simulations with the System for Atmospheric Modeling (SAM) and using GCM‐coupled CRM simulations with the Energy Exascale Earth System Model (E3SM)‐MMF confirm the heuristic model results.

CRM↗
DOE OSTI2022

Turbulent Rotating Rayleigh–Bénard Convection

Rotation with thermally induced buoyancy governs many astrophysical and geophysical processes in the atmosphere, ocean, sun, and Earth's liquid-metal outer core. Rotating Rayleigh–Bénard convection (RRBC) is an experimental system that has features of rotation and buoyancy, where a container of height H and temperature difference Δ between its bottom and top is rotated about its vertical axis with angular velocity Ω. The strength of buoyancy is reflected in the Rayleigh number (~ H 3 Δ) and that of the Coriolis force in the Ekman and Rossby numbers (~Ω -1 ). Rotation suppresses the convective onset, introduces instabilities, changes the velocity boundary layers, modifies the shape of thermal structures from plumes to vortical columns, affects the large-scale circulation, and can decrease or enhance global heat transport depending on buoyant and Coriolis forcing. RRBC is an extremely rich system, with features directly comparable to geophysical and astrophysical phenomena. Here we review RRBC studies, suggest a unifying heat transport scaling approach for the transition between rotation-dominated and buoyancy-dominated regimes in RRBC, and discuss non-Oberbeck–Boussinesq and centrifugal effects.

42 ENGINEERING↗
DOE OSTI2025

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗
DOE OSTI2021

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Synoptic and Mesoscale Modulation of Dynamic and Thermodynamic Impacts on Central Arctic Sea Ice During MOSAiC

The primary objectives of this proposed project are to determine 1) to what extent dynamic and thermodynamic air‐ice interactions in the Central Arctic are driven by atmospheric synoptic and mesoscale phenomena, and 2) the key processes involved in these interactions. These phenomena include large‐scale features such as fronts, high‐ and low‐pressure systems, low‐level jets, and cloud bands, but may also include smaller‐scale phenomena such as mesoscale instabilities (e.g., symmetric instability or slantwise convection), transient meso‐lows and gravity waves. Dynamic air‐ice interactions include direct atmospheric stress and atmospheric stress divergence/deformation; thermodynamic air‐ice interactions involve the net surface energy budget, especially in association with non‐stratocumulus clouds generally occurring with stronger dynamical forcing.

58 GEOSCIENCES↗
DOE OSTI2026

Onset of Rayleigh-Bénard convection at fuel-ablator interfaces in inertial confinement fusion simulations

Here, we demonstrate the onset of Rayleigh-Bénard convection (RBC) near the fuel-ablator interface of inertial confinement fusion (ICF) implosion simulations for the scenario in which Rayleigh-Taylor or Richtmyer-Meshkov instabilities are suppressed. This convective heating within the cold fuel can degrade the adiabat and therefore performance of ICF implosions even without mixing of fuel and ablator. Our work suggests additional considerations to realize design adiabats in otherwise ideal ICF implosions: first, via control of the temperature differences between ablator and fuel and, second, control of the fuel acceleration which we show has a complex interplay with the RBC dynamics.

Physics - Plasma physics↗
DOE OSTI2020

On the character of turbulent-like flows in self-consistent models of core-collapse supernovae

Neutrino-driven convection plays a crucial role in the development of core-collapse supernova (CCSN) explosions. However, the complex mechanism that triggers the shock revival and the subsequent explosion has remained inscrutable for many decades. Multidimensional simulations suggest that the growth of fluid instabilities, especially development of convection, will determine the morphology of the explosion. We have performed 3D simulations using spherical-polar coordinates covering a reduced angular extent (90° computational domain), and with angular resolutions of 2°, 1°, 1/2°, and 1/4°, to study the development of turbulent-like flows in core-collapse supernova explosions on a time scale of order 100 ms. We have employed the multi-physics Chimera code that includes detailed nuclear physics and spectral neutrino transport. Coarse resolution models do not develop an inertial range, presumably due to the bottleneck effect, such that the energy is prevented from cascading down to small scales and tends to accumulate at large scales. High-resolution models instead, start to recover the k −5/3 scaling of Kolmogorov's theory. Stochasticity and few simulation samples limit our ability to predict the development of explosions. Over the simulated time period, our models show no clear trend in improving (or diminishing) conditions for explosion as the angular resolution is increased. However, we find that disordered flow provides an effective pressure, as characterized by the Reynolds stress, behind the shock. This contribution, ~40%–50% of the thermal pressure, aids shock revival and thus the development of the explosion. Finally, we show that the kinetic energy power spectrum of reduced angular extent and full 4π models are consistent, thus indicating that a 90° computational domain is an adequate configuration to study the character of turbulence in CCSNe.

Casanova Bustamante, Jordi↗
DOE OSTI2020

Laser-plasma instabilities in long scale-length plasmas relevant to shock-ignition

We present 2D particle-in-cell simulations of laser-plasma instabilities (LPI) performed using conditions relevant to experiments on the OMEGA and NIF laser facilities. The laser intensity used is 2×10 15 W cm -2 , which allows comparison with previous experiments and simulations investigating shock ignition, though this is a relatively low intensity for the scheme. We find that the large convective gain of the instabilities leads to dynamics that are largely controlled by pump depletion and to scattering or absorption by LPIs predominantly occurring away from the quarter-critical density. The two-plasmon decay (TPD) instability is dominant for OMEGA-type conditions and remains important at the NIF-scale. In both cases, most absorption from TPD occurs near its Landau cutoff density, and we examine its nonlinear dynamics in this region, focusing on the ion-acoustic waves responsible for saturation. For NIF-type conditions, stimulated Raman scattering (SRS) and stimulated Brillouin scattering divert significant fractions of laser energy, and we find that SRS backscatter can occur with high gain at low densities due to kinetic effects. An in-depth analysis of hot-electron production is performed to identify the main sources and their characteristic electron temperature. For both simulations presented, the overall hot-electron temperature lies between 30 and 35 keV. This is relatively low and will likely increase shock pressure; however, the distributions also contain a significant number of high-energy (E k >100 keV) electrons that would likely cause unacceptable preheat. We, therefore, suggest some strategies that may be utilized to minimize this high-energy component. I. INTROD

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Wall modes and the transition to bulk convection in rotating Rayleigh-Bénard convection

We investigate states of rapidly rotating Rayleigh-Bénard convection in a cylindrical cell over a range of Rayleigh numbers 3 × 10 5 ≤ Ra ≤ 5 × 10 9 and Ekman numbers 10 − 6 ≤ Ek ≤ 10 − 4 for Prandtl number Pr = 0.8 and aspect ratios 1 / 5 ≤ Γ ≤ 5 using direct numerical simulations. We characterize, for perfectly insulating sidewall boundary conditions, the first transition to convection via wall mode instability and the nonlinear growth and instability of the resulting wall mode states, including a secondary transition to time dependence. We show how the radial structure of the vertical velocity u z and the temperature T is captured well by the linear eigenfunctions of the wall mode instability where the radial width of u z is δ u z ∼ Ek 1 / 3 r / H whereas δ T ∼ e − k r ( k is the wave number of a laterally infinite wall mode state). The disparity in spatial scales for Ek = 10 − 6 means that the heat transport is dominated by the radial structure of u z since T varies slowly over the radial scale δ u z . We further describe how the transition to a state of bulk convection is influenced by the presence of the wall mode states. We use temporal and spatial scales as measures of the local state of convection and the Nusselt number Nu as representative of global transport. Our results elucidate the evolution of the wall state of rotating convection and confirm that wall modes are strongly linked with the boundary zonal flow being the robust remnant of nonlinear wall mode states. We also show how the heat transport ( Nu ) contributions of wall modes and bulk modes are related and discuss approaches to disentangling their relative contributions. Published by the American Physical Society 2024

58 GEOSCIENCES↗
DOE OSTI2024

Tornadoes in Southeast South America: Mesoscale to Planetary-Scale Environments

Abstract A multiscale analysis of the environment supporting tornadoes in southeast South America (SESA) was conducted based on a self-constructed database of 74 reports. Composites of environmental and convective parameters from ERA5 were generated relative to tornado events. The distribution of the reported tornadoes maximizes over the Argentine plains, while events are rare close to the Andes and south of Sierras de Córdoba. Events are relatively common in all seasons except in winter. Proximity environment evolution shows enhanced instability, deep-layer vertical wind shear, storm-relative helicity, reduced convective inhibition, and a lowered lifting condensation level before or during the development of tornadic storms in SESA. No consistent signal in low-level wind shear is seen during tornado occurrence. However, a curved hodograph with counterclockwise rotation is present. The Significant Tornado Parameter (STP) is also maximized prior to tornadogenesis, most strongly associated with enhanced CAPE. Differences in the convective environment between tornadoes in SESA and the U.S. Great Plains are discussed. On the synoptic scale, tornado events are associated with a strong anomalous trough crossing the southern Andes that triggers lee cyclogenesis, subsequently enhancing the South American low-level jet (SALLJ) that increases moisture advection to support deep convection. This synoptic trough also enhances vertical shear that, along with enhanced instability, sustains organized convection capable of producing tornadic storms. At planetary scales, the tornadic environment is modulated by Rossby wave trains that appear to be forced by convection near northern Australia. Madden–Julian oscillation phase 3 preferentially occurs 1–2 weeks ahead of tornado occurrence. Significance Statement The main goal of this study is to describe what atmospheric conditions (from local to global scales) are present prior to and during tornadic storms impacting southeast South America (SESA). Increasing potential for deep convection, wind shear, and potential for rotating updrafts, as well as reducing convective inhibition and cloud-base height, are predominant a few hours before and during the events in connection to low-level northerly winds enhancing moisture transport to the region. Remote convective activity near northern Australia appears to influence large-scale atmospheric circulation that subsequently triggers convective storms supporting tornadogenesis 1–2 weeks later in SESA. Our findings highlight the importance of accounting for atmospheric processes occurring at different scales to understand and predict tornado occurrences.

Meteorology & Atmospheric Sciences↗