Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixed convection”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Benchmarking CO₂ storage simulations: Results from the 11 th Society of Petroleum Engineers Comparative Solution Project

The 11 th Society of Petroleum Engineers Comparative Solution Project (shortened SPE11 herein) benchmarked simulation tools for geological carbon dioxide (CO 2 ) storage. A total of 45 groups from leading research institutions and industry across the globe signed up to participate, with 18 ultimately contributing valid results that were included in the comparative study reported here. This paper summarizes the SPE11 results. A comprehensive introduction and qualitative discussion of the submitted data are provided, together with an overview of online resources for accessing the full depth of data. A global metric for analyzing the relative distance between submissions is proposed and used to conduct a quantitative analysis of the submissions. This analysis attempts to statistically resolve the key aspects influencing the variability between submissions. The study shows that the major qualitative variation between the submitted results is related to thermal effects, dissolution-driven convective mixing, and resolution of facies discontinuities. Moreover, a strong dependence on grid resolution is observed across all three versions of the SPE11. However, our quantitative analysis suggests that the observed variations are predominantly influenced by factors not documented in the technical responses provided by the participants. We therefore identify that unreported variations due to human choices within the process of setting up, conducting, and reporting on the simulations underlying each SPE11 submission are at least as impactful as the computational choices reported.

Nordbotten, Jan M. [Univ. of Bergen (Norway); Norw↗
DOE OSTI2023

A unified non-equilibrium phase change model for injection flow modeling

The homogenous relaxation model (HRM) is one of the most widely used models to describe the liquid- gas phase transition. However, in its original formulation, it is unable to handle multispecies vapor-liquid equilibrium (VLE), which limits its applicability to single-component fluids. In this work, a unified non-equilibrium phase change model that considers the VLE of multicomponent mixtures is proposed building upon the HRM's structure. A time factor is introduced to mimic the effect of different phase change timescales due to different mechanisms, e.g., cavitation, flash-boiling, and evaporation. Here to assess the model's performance, computational fluid dynamics simulations of the internal and near-nozzle injection flow of the Engine Combustion Network's Spray G injector were performed using the nine-component PACE-20 fuel with both the unified model and the original HRM. The predicted fuel density in the near-nozzle region matched well with X-ray tomography measurements. The simulation results indicated that, whereas the HRM failed to capture the vaporization due to convective mixing between the fuel and ambient gas, the unified model performed well in predicting the mixing-driven vaporization and the corresponding evaporative cooling. Further comparisons using the nine-component fuel formula and a single-component fuel surrogate demonstrated the unified model's ability to predict preferential vaporization, which affects the predictions of local mixture composition and rate of vaporization. Finally, it is shown that the unified model is capable of representing multiple phase change mechanisms, and the relaxation time factor plays an important role in determining the degree of phase change due to the different mechanisms.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Perfect Crystals: microgravity capillary counterdiffusion crystallization of human manganese superoxide dismutase for neutron crystallography

The NASA mission Perfect Crystals used the microgravity environment on the International Space Station (ISS) to grow crystals of human manganese superoxide dismutase (MnSOD)—an oxidoreductase critical for mitochondrial vitality and human health. The mission’s overarching aim is to perform neutron protein crystallography (NPC) on MnSOD to directly visualize proton positions and derive a chemical understanding of the concerted proton electron transfers performed by the enzyme. Large crystals that are perfect enough to diffract neutrons to sufficient resolution are essential for NPC. This combination, large and perfect, is hard to achieve on Earth due to gravity-induced convective mixing. Capillary counterdiffusion methods were developed that provided a gradient of conditions for crystal growth along with a built-in time delay that prevented premature crystallization before stowage on the ISS. Here, we report a highly successful and versatile crystallization system to grow a plethora of crystals for high-resolution NPC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Experiments and gyrokinetic simulations of the nonlinear interaction between spinning magnetized plasma pressure filaments

A set of experiments using controlled, skin depth-sized plasma pressure filaments in close proximity have been carried out in a large linear magnetized plasma device. Two- and three-filament configurations have been used to determine the scale of cross field nonlinear interaction. When the filaments are separated by a distance of approximately five times the size of a single filament or less, a significant transfer of charge and energy occurs, leading to the generation of inter-filament electric fields. This has the effect of rotating the filaments and influencing the merging dynamics. Nonlinear gyrokinetic simulations using seeded filaments confirm the presence of unstable drift-Alfvén modes driven by the steep electron temperature gradient. When the filaments are within a few collisionless electron skin depths (separations twice the size of a single filament), the unstable perturbations drive the convective mixing of the density and temperature and rearrange the gradients such that they maximize in the region surrounding the filament bundle.

Sydora, R. D. (ORCID:0000000192543149)↗
DOE OSTI2022

Toward full simulations for a liquid metal blanket: part 2. Computations of MHD flows with volumetric heating for a PbLi blanket prototype at Ha ~10 4 and Gr ~10 12

On the pathway toward full simulations for a liquid metal (LM) blanket, this part 2 extends a previous study of purely magnetohydrodynamic (MHD) flows in a DCLL blanket in reference Chen et al (2020 Nucl. Fusion 60 076003) to more general conditions when the MHD flow is coupled with heat transfer. The simulated prototypic blanket module includes all components of a real LM blanket system, such as supply ducts, inlet and outlet manifolds, multiple poloidal ducts and a U-turn zone. Volumetric heating generated by fusion neutrons is added to simulate thermal effects in the flowing lead–lithium (PbLi) breeder. The MHD flow equations and the energy equation are solved with a DNS-type finite-volume code ‘MHD-UCAS’ on a very fine mesh of 470 × 10 6 cells. The applied magnetic field is 5 T (Hartmann number Ha ~ 10 4 ), the PbLi velocity in the poloidal ducts is 10 cm s –1 (Reynolds number Re ~10 5 ), whereas the maximum volumetric heating is 30 MW m –3 (Grashof number Gr ~ 10 12 ). Four cases have been simulated, including forced- and mixed-convection flows, and either an electrically conducting or insulating blanket structure. Various comparisons are made between the four computed cases and also against the purely MHD flows computed earlier in reference Chen et al (2020 Nucl. Fusion 60 076003) with regards to the (1) MHD pressure drop, (2) flow balancing, (3) temperature field, (4) flows in particular blanket components, and (5) 3D and turbulent flow effects. The strongest buoyancy effects were found in the poloidal ducts. In the electrically non-conducting blanket, the buoyancy forces lead to significant modifications of the flow structure, such as formation of reverse flows, whereas their effect on the MHD pressure drop is relatively small. In the electrically conducting blanket case, the buoyancy effects on the flow and MHD pressure drop are almost negligible.

Physics↗
DOE OSTI2020

Modelling R Coronae Borealis stars: effects of He-burning shell temperature and metallicity

The R Coronae Borealis (RCB) stars are extremely hydrogen-deficient carbon stars that produce large amounts of dust, causing sudden deep declines in brightness. They are believed to be formed primarily through white dwarf mergers. In this paper, we use mesa to investigate how post-merger objects with a range of initial He-burning shell temperatures from 2.1 to 5.4 × 10 8 K with solar and subsolar metallicities evolve into RCB stars. The most successful model of these has subsolar metallicity and an initial temperature near 3 × 10 8 K. We find a strong dependence on initial He-burning shell temperature for surface abundances of elements involved in the CNO cycle, as well as differences in effective temperature and radius of RCBs. Elements involved in nucleosynthesis present around 1 dex diminished surface abundances in the 10 per cent solar metallicity models, with the exception of carbon and lithium that are discussed in detail. Models with subsolar metallicities also exhibit longer lifetimes than their solar counterparts. In addition, we find that convective mixing of the burned material occurs only in the first few years of post-merger evolution, after which the surface abundances are constant during and after the RCB phase, providing evidence for why these stars show a strong enhancement of partial He-burning products.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A spectrophotometric analysis of cool white dwarfs in the Gaia and pan-STARRS footprint

We present a spectrophotometric analysis of 2880 cool white dwarfs within 100 pc of the Sun and cooler than T eff ∼ 10 000 K, with grizy Pan-STARRS photometry and Gaia trigonometric parallaxes available. We also supplement our data sets with near-infrared JHK photometry, when available, which is shown to be essential for interpreting the coolest white dwarfs in our sample. We perform a detailed analysis of each individual object using state-of-the-art model atmospheres appropriate for each spectral type, including DA, DC, DQ, DZ, He-rich DA, and the so-called IR-faint white dwarfs. We discuss the temperature and mass distributions of each subsample, as well as revisit the spectral evolution of cool white dwarfs. We find little evidence in our sample for the transformation of a significant fraction of DA stars into He-atmosphere white dwarfs through the process of convective mixing between T eff = 10 000 and ∼6500 K, although the situation changes drastically in the range T eff = 6500–5500 K where the fraction of He-atmosphere white dwarfs reaches ∼45 per cent. However, we also provide strong evidence that at even cooler temperatures (T eff ≲ 5200 K), most DC white dwarfs have H atmospheres. We discuss a possible mechanism to account for this sudden transformation from He- to H-atmosphere white dwarfs involving the onset of crystallization and the occurrence of magnetism. Finally, we also argue that DQ, DZ, and DC white dwarfs may form a more homogeneous population than previously believed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Temporal Forecasting of Distributed Temperature Sensing in a Thermal Hydraulic System With Machine Learning and Statistical Models

We benchmark performance of long-short term memory (LSTM) network machine learning model and autoregressive integrated moving average (ARIMA) statistical model in temporal forecasting of distributed temperature sensing (DTS). Data in this study consists of fluid temperature transient measured with two co-located Rayleigh scattering fiber optic sensors (FOS) in a forced convection mixing zone of a thermal tee. We treat each gauge of a FOS as an independent temperature sensor. We first study prediction of DTS time series using Vanilla LSTM and ARIMA models trained on prior history of the same FOS that is used for testing. The results yield maximum absolute percentage error (MaxAPE) and root mean squared percentage error (RMSPE) of 1.58% and 0.06% for ARIMA, and 3.14% and 0.44% for LSTM, respectively. Next, we investigate zero-shot forecasting (ZSF) with LSTM and ARIMA trained on history of the co-located FOS only, which is advantageous when limited training data is available. The ZSF MaxAPE and RMSPE values for ARIMA are comparable to those of the Vanilla use case, while the error values for LSTM increase. We show that in ZSF, performance of LSTM network can be improved by training on most correlated gauges between the two FOS, which are identified by calculating the Pearson correlation coefficient. The improved ZSF MaxAPE and RMSPE for LSTM are 4.4% and 0.33%, respectively. Performance of ZSF LSTM can be further enhanced through transfer learning (TL), where LSTM is re-trained on a subset of the FOS that is the target of forecasting. We show that LSTM pre-trained on correlated dataset and re-trained on 30% of testing target dataset achieves MaxAPE and RMSPE values of 2.32% and 0.28%, respectively.

ARIMA↗
DOE OSTI2021

Physics-Informed Neural Networks for Heat Transfer Problems

Abstract Physics-informed neural networks (PINNs) have gained popularity across different engineering fields due to their effectiveness in solving realistic problems with noisy data and often partially missing physics. In PINNs, automatic differentiation is leveraged to evaluate differential operators without discretization errors, and a multitask learning problem is defined in order to simultaneously fit observed data while respecting the underlying governing laws of physics. Here, we present applications of PINNs to various prototype heat transfer problems, targeting in particular realistic conditions not readily tackled with traditional computational methods. To this end, we first consider forced and mixed convection with unknown thermal boundary conditions on the heated surfaces and aim to obtain the temperature and velocity fields everywhere in the domain, including the boundaries, given some sparse temperature measurements. We also consider the prototype Stefan problem for two-phase flow, aiming to infer the moving interface, the velocity and temperature fields everywhere as well as the different conductivities of a solid and a liquid phase, given a few temperature measurements inside the domain. Finally, we present some realistic industrial applications related to power electronics to highlight the practicality of PINNs as well as the effective use of neural networks in solving general heat transfer problems of industrial complexity. Taken together, the results presented herein demonstrate that PINNs not only can solve ill-posed problems, which are beyond the reach of traditional computational methods, but they can also bridge the gap between computational and experimental heat transfer.

Engineering↗
DOE OSTI2022

Characterizing warm atmospheric boundary layer over land by combining Raman and Doppler lidar measurements

PBL plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, for PBL characterization and process study, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary. Here, new approaches are developed to provide reliable PBLH and MLH to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) user facility. Compared with widely used lidar aerosol measurements for PBLH, WVMR is a better trace for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to re-construct 2-D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm seasons’ RL and ML measurements at the SGP site for five years to study warm season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable PBLH and MLH dataset will be valuable for PBL process study, model evolution, and PBL parameterization improvement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High-Resolution Simulations of Geological CO 2 Injection: Application to the SPE11 Benchmark

Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.

Geosciences↗
DOE OSTI2024

Modulation of regional carbon uptake by AMOC and alkalinity changes in the subpolar North Atlantic under a warming climate

The slowdown of the Atlantic Meridional Overturning Circulation (AMOC) and associated consequences on ocean carbon uptake could have large implications for the Earth's climate system and its global carbon cycle. This study analyzes ten Earth System Models from the Coupled Model Intercomparison Project Phase 6 (CMIP6) and reveals that the regional carbon uptake in the subpolar North Atlantic under a high CO 2 emission scenario moderately correlates with the decline in AMOC at 40°N. AMOC transports warm and salty subtropical waters to the subpolar regions. Models with stronger AMOC slowdown generally exhibit weaker surface warming and larger decline of surface salinity and alkalinity. We consider two plausible mechanisms linking the AMOC slowdown to the decline of regional CO 2 uptake: the reduction in surface alkalinity and diminished subduction. The decline of surface salinity and alkalinity reduces the ocean's capacity to buffer acids leading to a reduced CO 2 uptake. This important contribution is unique to the North Atlantic. Diminished convective mixing and subduction of surface water can further decrease the downward transport of anthropogenic carbon, as also shown in previous research. The centennial trends of p CO 2 are decomposed into four components driven by temperature, salinity, alkalinity and dissolved inorganic carbon, revealing that alkalinity and dissolved inorganic carbon are both significant contributors. The alkalinity-driven p CO 2 essentially follows surface salinity, establishing the linkage between AMOC slowdown and alkalinity decline. Our results indicate that alkalinity changes are important for the interplay between AMOC and the regional carbon sequestration ability across the late 20 th and the entirety of the 21 st century in the subpolar North Atlantic.

AMOC↗
DOE OSTI2024

Modeling the Progenitor Stars of Observed Type IIP Supernovae

Abstract Type IIP supernovae (SNe IIP) are thought to originate from the explosion of massive stars >10M ⊙ . Their luminosity is primarily powered by the explosion energy and the radioactive decay energy of 56 Co, with the photosphere location regulated by hydrogen recombination. However, the physical connections between SNe IIP and their progenitor stars remain unclear. This paper presents a comprehensive study of SNe IIP and their progenitor stars by using the one-dimensional stellar evolution code,MESA. Our model grids consider the effects of stellar metallicity, mass, and rotation in the evolution of massive stars, as well as the explosion energy and 56 Ni production in modeling supernovae. To elucidate the observed SNe IIP and their origins, we compare their light curves (LCs) with our models. Furthermore, we investigate the impact of stellar parameters on LCs by considering stellar mass, metallicity, rotation, explosion energy, and 56 Ni production. We find that more massive stars exhibit longer plateaus due to increased photon diffusion time caused by massive ejecta. Higher metallicity leads to increased opacity and mass loss of progenitor stars. Rapid rotation affects internal stellar structures, enhancing convective mixing and mass loss, potentially affecting the plateau’s brightness and duration. Higher explosion energy results in brighter but shorter plateaus due to faster-moving ejecta. 56 Ni mass affects late-time luminosity and plateau duration, with larger masses leading to slower declines.

Astronomy & Astrophysics↗
DOE OSTI2025

ARM SGP PBLH and MLH datasets from Raman lidar and Doppler lidar

The planetary boundary layer (PBL) plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including the daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary for studying PBL characterization and processes. Here, new approaches are developed to provide reliable PBLH and MLH estimates to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) User Facility. Compared to widely used lidar aerosol measurements for PBLH, WVMR is a better tracer for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to reconstruct 2D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm-season RL and ML measurements at the SGP site for five years to study warm-season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in a warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable data set of PBLH and MLH values will be valuable for studying PBL processes, model evolution, and PBL parameterization improvements. The MLH dataset includes the MLH in values of km above ground level. The PBLH dataset includes the PBLH in values of km above ground level, along with a flag ("situation_PBLH") to determine the state of the PBL (1 = Cloudy Condition, 2 = Stable Layer, 3 = Multi-layer WVMR structure, 4 = Well-Mixed PBL, 5 = A de-coupled layer, 6 = Other).

mixing layer height↗
DOE OSTI2020

Coarse grained simulation of convectively driven turbulent mixing, transition, and turbulence decay

Accurate predictions with quantifiable uncertainties are essential to many practical turbulent flow applications exhibiting extreme geometrical complexity and broad ranges of length and time scales. Under-resolved computer simulations are typically unavoidable in such applications, and implicit large-eddy simulation (ILES) often becomes the effective strategy. In this work, we focus on assessing ILES initialized with well-characterized 25 6 3 homogeneous isotropic turbulence datasets generated with direct numerical simulation (DNS). ILES solutions based on the LANL xRAGE code are studied as function of resolution for 6 4 3 , 12 8 3 , 25 6 3 , and 51 2 3 grids. ILES performance of new directionally-unsplit high-order numerical hydrodynamics algorithms in xRAGE is examined in this context. Compared to the initial 25 6 3 DNS, we find longer inertial subranges and higher turbulence Reynolds number (Re) for 25 6 3 and 51 2 3 xRAGE — attributed to having linked DNS (Navier–Stokes based) solutions to nominally inviscid (higher Re) Euler based ILES solutions. For fixed 25 6 3 resolution, we find that significantly larger simulated turbulence Re can be achieved with the higher-order unsplit (vs. split) discretizations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Impacts of Stochastic Mixing in Idealized Convection-Permitting Simulations of Squall Lines

This study investigates impacts of altering subgrid-scale mixing in “convection-permitting” kilometer-scale horizontal-grid-spacing (Δ h ) simulations by applying either constant or stochastic multiplicative factors to the horizontal mixing coefficients within the Weather Research and Forecasting Model. In quasi-idealized 1-km Δ h simulations of two observationally based squall-line cases, constant enhanced mixing produces larger updraft cores that are more dilute at upper levels, weakens the cold pool, rear-inflow jet, and front-to-rear flow of the squall line, and degrades the model’s effective resolution. Reducing mixing by a constant multiplicative factor has the opposite effect on all metrics. Completely turning off parameterized horizontal mixing produces bulk updraft statistics and squall-line mesoscale structure closest to an LES “benchmark” among all 1-km simulations, although the updraft cores are too undilute. The stochastic mixing scheme, which applies a multiplicative factor to the mixing coefficients that varies stochastically in time and space, is employed at 0.5-, 1-, and 2-km Δ h . It generally reduces midlevel vertical velocities and enhances upper-level vertical velocities compared to simulations using the standard mixing scheme, with more substantial impacts at 1- and 2-km Δ h compared to 0.5-km Δ h . Further, the stochastic scheme also increases updraft dilution to better agree with the LES for one case, but has less impact on the other case. Stochastic mixing acts to weaken the cold pool but without a significant impact on squall-line propagation. It also does not affect the model’s overall effective resolution unlike applying constant multiplicative factors to the mixing coefficients.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Evaluation of a Stochastic Mixing Scheme in the Deep Convective Gray Zone Using a Tropical Oceanic Deep Convection Case Study

Abstract A stochastic horizontal subgrid‐scale mixing scheme is evaluated in ensemble simulations of a tropical oceanic deep convection case using a horizontal grid spacing (Δ h ) of 3 km. The stochastic scheme, which perturbs the horizontal mixing coefficient according to a prescribed spatiotemporal autocorrelation scale, is found to generally increase mesoscale organization and convective intensity relative to a non‐stochastic control simulation. Perturbations applied at relatively short autocorrelation scales induce differences relative to the control that are more systematic than those from perturbations applied at relatively long scales that yield more variable outcomes. A simulation with mixing enhanced by a constant factor of 4 significantly increases mesoscale organization and convective intensity, while turning off horizontal subgrid‐scale mixing decreases both. Total rainfall is modulated by a combination of mesoscale organization, areal coverage of convection, and convective intensity. The stochastic simulations tend to behave more similarly to the constant enhanced mixing simulation owing to greater impacts from enhanced mixing as compared to reduced mixing. The impacts of stochastic mixing are robust, ascertained by comparing the stochastic mixing ensembles with a non‐stochastic mixing ensemble that has grid‐scale noise added to the initial thermodynamic field. Compared to radar observations and a higher resolution Δ h = 1 km simulation, stochastic mixing seemingly degrades the simulation performance. These results imply that stochastic mixing produces non‐negligible impacts on convective system properties and evolution but does not lead to an improved representation of convective cloud characteristics in the case studied here.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Hydrodynamical Simulations of Proton Ingestion Flashes in Type I X-Ray Bursts

We perform the first multidimensional fluid simulations of thermonuclear helium ignition underneath a hydrogen-rich shell. This situation is relevant to Type I X-ray bursts on neutron stars that accrete from a hydrogen-rich companion. Using the low-Mach number fluid code MAESTROeX, we investigate the growth of the convection zone due to nuclear burning, and the evolution of the chemical abundances in the atmosphere of the star. We also examine the convective boundary mixing processes that cause the evolution to differ significantly from previous one-dimensional simulations that rely on mixing-length theory. We find that the convection zone grows outward as penetrating fluid elements cool the overlying radiative layer, rather than directly from the increasing entropy of the convection zone itself. Simultaneously, these flows efficiently mix composition, carrying carbon out of and protons into the convection zone even before contact with the hydrogen shell. We discuss the implications of these effects for future modeling of these events and observations.

astrophysical fluid dynamics↗