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

The Three-dimensional Collapse of a Rapidly Rotating 16 M ⊙ Star

I report on the three-dimensional (3D) hydrodynamic evolution of a rapidly rotating 16 M ⊙ star to iron core collapse. For the first time, I follow the 3D evolution of the angular momentum (AM) distribution in the iron core and convective shell burning regions for the final 10 minutes up to and including gravitational instability and core collapse. In 3D, convective regions show efficient AM transport that leads to an AM profile that differs in shape and magnitude from MESA within a few shell convective turnover timescales. For different progenitor models, such as those with tightly coupled Si/O convective shells, efficient AM transport in 3D simulations could lead to a significantly different AM distribution in the stellar interior affecting estimates of the natal neutron star or black hole spin. The results suggest that 3D AM transport in convective and rotating shell burning regions are critical components in models of massive stars and could qualitatively alter the explosion outcome and inferred compact remnant properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Gravitational-wave signal of a core-collapse supernova explosion of a 15M ⊙ star

In this paper we report on the gravitational-wave signal computed in the context of a three-dimensional simulation of a core-collapse supernova explosion of a 15 M ⊙ star. The simulation was performed with our neutrino hydrodynamics code chimera . We detail the gravitational wave strains as a function of time, for both polarizations, and discuss their physical origins. We also present the corresponding spectral signatures. Gravitational wave emission in our model has two key features: low-frequency emission (less than 200 Hz) emanates from the gain layer as a result of neutrino-driven convection and the standing accretion shock instability (SASI), and high-frequency emission (greater than 600 Hz) emanates from the proto–neutron star due to Ledoux convection within it. The high-frequency emission dominates the gravitational wave emission in our model and emanates largely from the convective layer itself, not from the convectively stable layer above it, due to convective overshoot. Moreover, the low-frequency emission emanates from the gain layer itself, not from the proto–neutron star, due to accretion onto it. We provide evidence of the SASI in our model and demonstrate that the peak of our low-frequency gravitational wave emission spectrum corresponds to it. Given its origin in the gain layer, we classify the SASI emission in our model as p-mode emission and assign a purely acoustic origin, not a vortical–acoustic origin, to it. We compare the results of our three-dimensional model analysis with those obtained from the model’s two-dimensional counterpart and find a significant reduction in the strain amplitudes in the former case, as well as significant reductions in all related quantities. Our dominant proto–neutron star gravitational wave emission is not well characterized by emission from surface g modes, complicating the relationship between peak frequencies observed and the mass and radius of the proto–neutron star expressed by analytic estimates under the assumption of surface g-mode emission. We present our frequency normalized characteristic strain along with the sensitivity curves of current- and next-generation gravitational wave detectors. This simple analysis indicates that the spectrum of gravitational wave emission between approximately 20 Hz and approximately 1 kHz, stemming from neutrino-driven convection, the SASI, accretion onto the proto–neutron star, and proto–neutron star convection, will be accessible for a Galactic event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Theory of the magnetothermal instability in coronal plasma flows

The theory of the magnetothermal instability (MTI) is revisited through the lens of the stability of uniform systems. The linear stability analysis includes flow advection and Nernst transport. The instability criteria derived distinguish between the convective and the absolute nature of the perturbation growth. It is proven that, in the region where the Nernst and plasma blowoff velocities cancel, the MTI can be absolute and wave-packet perturbations grow in situ. This instability is mediated by the internal feedback between the Biermann battery and Righi–Leduc terms. Here, the analysis is extended to derive the dispersion relation for short-wavelength perturbations developing in nonuniform profiles with the application to coronal plasmas. It is found that the condition for MTI requires the net B-field convection velocity to be small at the isothermal sonic section, and the plasma conditions in this section govern the dynamics of the instability. Analysis of hydro-equivalent implosions suggests that unstable perturbations undergo more e-foldings of growth in larger-size targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Geometric optics model of Thomson scattering enhanced by parametric coupling

We formulate a theory of Thomson scattering (TS) from inhomogeneous plasmas using a ray-tracing approach. The scattered radiation is described within the eikonal approximation, including convective gain arising from scattering instabilities driven by the TS probe laser. In the limit of a weak probe and vanishing gain, the theory recovers the standard TS cross section and scattered power due to volumetric particle noise. The ray-tracing formalism is implemented in a TS fitting framework, enabling accurate analysis of inhomogeneous plasmas whose background profiles can be provided by large-scale hydrodynamic simulations or experimental diagnostics. Here, several examples of synthetic TS spectra are presented, and further extensions of the theory are discussed.

Carleton, D. E. [University of Alberta, Edmonton, ↗
DOE OSTI2022

Laser melting modes in metal powder bed fusion additive manufacturing

In laser powder bed fusion additive manufacturing of metals, extreme thermal conditions create many highly dynamic physical phenomena such as vaporization and recoil, Marangoni convection, and protrusion and keyhole instability. Collectively however, the full set of phenomena is too complicated for practical applications and, in reality, the melting modes are used as a guideline for printing. With increasing local material temperature beyond the boiling point, the mode can change from conduction to keyhole. These mode designations ignore laser-matter interaction details but in many cases are adequate to determine the approximate microstructures and hence the properties of the build. To date, no consistent, common, and coherent definitions have been agreed upon because of historic limitations in melt pool and vapor depression morphology measurements. Here, we distinguish process-based definitions of different melting modes from those based on postmortem evidence. The latter are mainly derived from the transverse cross-sections of the fusion zone, whereas the former come directly from time-resolved x-ray imaging of melt pool and vapor depression morphologies. These process-based definitions are more strict and physically sound, and they offer new guidelines for laser additive manufacturing practices and create new research directions. Further, we highlight the significance of the keyhole, which substantially enhances the laser energy absorption by the melt pool. Recent studies strongly suggest that stable-keyhole laser melting enables efficient, sustainable, and robust additive manufacturing. The realization of this scenario demands the development of multiphysics models, signal translations from morphology to other feasible signals, and in-process metrology across platforms and scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Convection-Permitting Ensembles of an Isolated Mountain Thunderstorm during RELAMPAGO/CACTI

Abstract The north–south-oriented Sierras de Córdoba (SDC) ridge in central Argentina is noted for initiating thunderstorms that may grow into intense mesoscale convective systems (MCSs). It also initiates more isolated, shorter-lived cells under weaker synoptic forcing. These cells are less impactful than MCSs but may be difficult to predict in convective-scale numerical weather prediction (NWP) due to their strong sensitivities to subgrid and partially resolved processes. To study the mechanisms and predictability of such cells, convection-permitting ensemble simulations were conducted of an isolated, diurnally forced SDC thunderstorm during Cloud, Aerosol, and Complex Terrain Interactions (CACTI)/Remote Sensing of Electrification, Lightning, and Mesoscale/Microscale Processes with Adaptive Ground Observations (RELAMPAGO). The rich observational data facilitated detailed ensemble verification, where dry biases in the surface energy balance and soil moisture were identified. These biases promoted rapid removal of convective inhibition and an early onset of precipitating cells over the SDC that were shallower and weaker than the observed cell. Correction, and then overcorrection, of the soil moisture bias in two successive ensembles was required to rectify the surface energy balance and improve the representation of the SDC cell. Nevertheless, substantial ensemble variability in convective precipitation was found, with some members producing more widespread convection than observed and others producing no deep convection at all. This variability was largely explained by a combination of thermodynamic and dynamic mechanisms, dominated by a positive sensitivity of convective precipitation to preconvective moist instability over the ridge. Secondary sensitivities were found to low-level upward mass flux and midlevel cross-barrier winds, the latter of which caused gravity waves with elevated downdrafts that tended to suppress incipient clouds.

Lopez, Andres [Department of Atmospheric and Ocean↗
DOE OSTI2020

Electrodeposition of Zinc in Aqueous Electrolytes Containing High Molecular Weight Polymers

Strong polarization of the ion distribution in liquid electrolytes subjected to potential differences exceeding the thermal voltage, V T = kT/e, produces a hydrodynamic instability termed electroconvection at ion-selective interfaces. Electroconvection is desirable in some situations (e.g., electrodialysis) because it promotes mixing in a stagnant electrolyte layer, enhancing the ion flux at a fixed potential difference. It is undesirable in others (e.g., electrodeposition of metals) where early experiments show that convective fluid rolls associated with hydrodynamic instability in bounded electrolytes produce preferential metal deposition at localized regions on an electrode. Such localized deposition drives the family of morphological instabilities loosely termed dendritic electrodeposition. We experimentally investigate the effect of ultrahigh molecular weight polymer additives on the onset conditions and physical characteristics of both instabilities. Direct observations of electrodeposit morphology and tracer particle motions used in tandem with indirect electrokinetic measurements reveal that even at moderate concentrations, the polymer additives have a large beneficial effect in extending the range of electric potentials where stable electrodeposition is observed. Additionally, we report that at polymer concentrations above the entanglement threshold, high molecular weight polymers impart elasticity to liquid electrolytes which dampen electroconvective flow at a cation-selective interface but have at most a minimal effect on the bulk ionic conductivity of the liquid.

25 ENERGY STORAGE↗
DOE OSTI2020

The Role of Multiscale Interaction in Tropical Cyclogenesis and Its Predictability in Near-Global Aquaplanet Cloud-Resolving Simulations

Tropical cyclogenesis (TCG) is a multiscale process that involves interactions between large-scale circulation and small-scale convection. A near-global aquaplanet cloud-resolving model (NGAqua) with 4-km horizontal grid spacing that produces tropical cyclones (TCs) is used to investigate TCG and its predictability. This study analyzes an ensemble of three 20-day NGAqua simulations, with initial white-noise perturbations of low-level humidity. TCs develop spontaneously from the northern edge of the intertropical convergence zone (ITCZ), where large-scale flows and tropical convection provide necessary conditions for barotropic instability. Zonal bands of positive low-level absolute vorticity organize into cyclonic vortices, some of which develop into TCs. A new algorithm is developed to track the cyclonic vortices. A vortex-following framework analysis of the low-level vorticity budget shows that vertical stretching of absolute vorticity due to convective heating contributes positively to the vorticity spinup of the TCs. A case study and composite analyses suggest that sufficient humidity is key for convective development. TCG in these three NGAqua simulations undergoes the same series of interactions. The locations of cyclonic vortices are broadly predetermined by planetary-scale circulation and humidity patterns associated with ITCZ breakdown, which are predictable up to 10 days. Whether and when the cyclonic vortices become TCs depend on the somewhat more random feedback between convection and vorticity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The Roles of Westward-Propagating Waves and the QBO in Limiting MJO Propagation

Abstract A recent study categorized the Madden–Julian oscillation (MJO) during boreal winter season into four types called stand, jump, slow, and fast MJO. This study focuses on the stand and jump MJO. Based on whether their convection penetrates the Maritime Continent (MC), stand and jump MJOs are seen as non-penetrating (NP) MJOs, while the other two are seen as eastward-penetrating (EP) MJOs. Results reveal the relative roles of the westward-propagating wave (WPW), as well as the QBO and ENSO, in limiting MJO propagation. Lack of the premoistening over the southern sea surface of the MC stops NP MJO from penetrating the MC. The active convection of the WPWs hinders the descending branch of the NP MJO circulation and therefore leads to the insufficient meridional advective moistening over the southern sea surface of the MC. The independent convection over the Pacific for jump MJOs is influenced by a combined effect of the QBO and ENSO. The tropopause instability induced by the MJO is found to significantly decouple from its convection over the Pacific in westerly QBO (QBOW) winters more than in easterly QBO (QBOE) winters. For jump MJOs, the independent convection over the central Pacific comes from local WPWs whose amplification and further development into deep convection are correlated to jump the MJOs’ decoupled tropopause instability. For stand MJOs, however, the seasonal-mean La Niña–like cool SST anomalies weaken the WPW activity over the central Pacific and confine WPWs within the western Pacific. Therefore, the decoupled tropopause instability of stand MJOs is out phase of WPWs and fails to induce an independent convection over the central Pacific.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Core collapse supernova gravitational wave emission for progenitors of 9.6, 15, and 25 M ⊙

Here, we present gravitational wave emission predictions based on three core collapse supernova simulations corresponding to three different progenitor masses. The masses span a large range, between 9.6 and 25M ⊙ , are all initially nonrotating, and are of two metallicities: zero and solar. We compute both the temporal evolution of the gravitational wave strains for both the plus and the cross polarizations, as well as their spectral decomposition and characteristic strains. The temporal evolution of our zero metallicity 9.6M ⊙ progenitor model is distinct from the temporal evolution of our solar metallicity 15M ⊙ progenitor model and our zero metallicity 25M ⊙ progenitor model. In the former case, the high-frequency gravitational wave emission is largely confined to a brief time period ~75 ms after bounce, whereas in the latter two cases high-frequency emission does not commence until ~125 ms after bounce or later. The excitation mechanisms of the high-frequency emission in all three cases correspond to proto-neutron star convection and accretion onto the proto-neutron star from the convective gain layer above it, with the former playing the dominant role for most of the evolution. The low-frequency emission in all three models exhibits very similar behavior. At frequencies below ~250 Hz, gravitational waves are emitted by neutrino-driven convection and the standing accretion shock instability (SASI). This emission extends throughout the simulations when a gain region is present. In all three models, explosion is observed at ~125, ~500, and ~250 ms after bounce in the 9.6, 15, and 25M ⊙ progenitor models, respectively. At these times, the low-frequency gravitational wave emission is joined by very low-frequency emission, below ~10 Hz. These very low-frequency episodes are the result of explosion and begin at the above designated explosion times in each of our models. Our characteristic strains tell us that, in principle, all three gravitational wave signals would be detectable by current-generation detectors for a supernova at a distance of 10 kpc. However, our 9.6M ⊙ progenitor model is a significantly weaker source of gravitational waves, with strain amplitudes approximately 5–10 times less than in our other two models. The characteristic strain for this model tells us that such a supernova would be detectable only within a much more narrow frequency range around the maximum sensitivity of today’s detectors. Finally, in our 9.6M ⊙ progenitor model, we see very high-frequency gravitational radiation, extending up to ~2000 Hz. This feature results from the interaction of shock- and deleptonization-induced convection with perturbations introduced in the progenitor by nuclear burning during core collapse. While unique to the 9.6M ⊙ progenitor model analyzed here, this very high-frequency emission may, in fact, be a generic feature of the predictions for the gravitational wave emission from all core collapse supernova models when simulations are performed with three-dimensional progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Flow Channels and the Generation of Alfvénic Turbulence Along Storm‐Time Inner Magnetospheric Field‐Lines

Abstract A feature of Earth's storm‐time magnetosphere outside the plasmapause is the occurrence of broad‐spectrum Alfvénic fluctuations. In this letter observations from the Van Allen Probes are compared with 3‐D fluid‐kinetic simulations of an evolving convective flow channel to investigate the mechanisms generating the observed spectrum. It is shown how narrow channels of fast convection are unstable to the Kelvin‐Helmholtz instability which on closed field‐lines initiates a cascade to small scales. Sustained driving of the flow combined with reflection from the topside ionosphere leads to the generation of an intensified spectrum of electromagnetic structures having similar spectral and morphological characteristics to those observed. This process couples enhanced magnetospheric convection to kinetic scale electromagnetic fluctuations that drive particle transport, scattering and energization through the outer radiation belt and ring current during geomagnetic storms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Symmetry-breaking-induced rare fluctuations in a time-delay dynamic system

Inspired by the experimental and numerical findings, we study the dynamic instabilities of two coupled nonlinear delay differential equations that are used to describe the coherent oscillations between the top and bottom boundary layers in turbulent Rayleigh–Bénard convection. By introducing two sensitivity parameters for the instabilities of the top and bottom boundary layers, we find three different types of solutions, namely in-phase single-period oscillations, multi-period oscillations and chaos. The chaos solution contains rare but large amplitude fluctuations. The statistical properties of these fluctuations are consistent with those observed in the experiment for the massive eruption of thermal plumes, which causes random reversals of the large-scale circulation in turbulent Rayleigh–Bénard convection. Finally, our study thus provides new insights into the origin of rare massive eruptions and sudden changes of large-scale flow pattern that are often observed in convection systems of geophysical and astrophysical scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

The Physics behind Precipitation Onset Bias in CMIP6 Models: The Pseudo-Entrainment Diagnostic and Trade-Offs between Lapse Rate and Humidity

Conditional instability and the buoyancy of plumes drive moist convection but have a variety of representations in model convective schemes. Vertical thermodynamic structure information from Atmospheric Radiation Measurement (ARM) sites and reanalysis (ERA5), satellite-derived precipitation (TRMM3b42), and diagnostics relevant for plume buoyancy are used to assess climate models. Previous work has shown that CMIP6 models represent moist convective processes more accurately than their CMIP5 counterparts. However, certain biases in convective onset remain pervasive among generations of CMIP modeling efforts. Here we diagnose these biases in a cohort of nine CMIP6 models with sub-daily output, assessing conditional instability in profiles of equivalent potential temperature θ e and saturation equivalent potential temperature θ es in comparison to a plume model with different mixing assumptions. Most models capture qualitative aspects of the θ es vertical structure, including a substantial decrease with height in the lower free troposphere associated with the entrainment of subsaturated air. We define a “pseudo-entrainment” diagnostic that combines subsaturation and a θ es measure of conditional instability similar to what entrainment would produce under the small-buoyancy approximation. This captures the trade-off between larger θ es lapse rates (entrainment of dry air) and small subsaturation (permits positive buoyancy despite high entrainment). This pseudo-entrainment diagnostic is also a reasonable indicator of the critical value of integrated buoyancy for precipitation onset. Models with poor θ e /θ es structure (those using variants of the Tiedtke Scheme) or low entrainment runs of CAM5, and models with low subsaturation, such as NASA-GISS, lie outside the observational range in this diagnostic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Land–Atmosphere Coupling at the U.S. Southern Great Plains: A Comparison on Local Convective Regimes between ARM Observations, Reanalysis, and Climate Model Simulations

Using the 9-yr warm-season observations at the Atmospheric Radiation Measurement Southern Great Plains site, in this work we assess the land–atmosphere (LA) coupling in the North American Regional Reanalysis (NARR) and two climate models: hindcasts with the Community Atmosphere Model version 5.1 by Cloud-Associated Parameterizations Testbed (CAM5-CAPT) and nudged runs with the Energy Exascale Earth System Model Atmosphere Model version 1 Regionally Refined Model (EAMv1-RRM). We focus on three local convective regimes and diagnose model behaviors using the local coupling metrics. NARR agrees well with observations except a slightly warmer and drier surface with higher downwelling shortwave radiation and lower evaporative fraction. On clear-sky days, it shows warmer and drier early-morning conditions in both models with significant underestimates in surface evaporation by EAMv1-RRM. On the majority of the ARM-observed shallow cumulus days, there is no or little low-level clouds in either model. When captured in models, the simulated shallow cumulus shows much less cloud fraction and lower cloud bases than observed. On the days with late-afternoon deep convection, models tend to present a stable early-morning lower atmosphere more frequently than the observations, suggesting that the deep convection is triggered more often by elevated instabilities. Generally, CAM5-CAPT can reproduce the local LA coupling processes to some extent due to the constrained early-morning conditions and large-scale winds. EAMv1-RRM exhibits large precipitation deficits and warm and dry biases toward mid-to-late summers, which may be an amplification through a positive LA feedback among initial atmosphere and land states, convection triggering and large-scale circulations.

58 GEOSCIENCES↗
DOE OSTI2019

Validation of RELAP-7 for Forced Convective and Natural Circulation Reactor Flows

New reactor designs such as small modular reactors and newly proposed passive safety systems rely on natural circulation two-phase flow cooling of the reactor core under accident conditions. Without a pump to drive the flow, the circulation rate is a strong function of void fraction. Under these conditions where the flow rate, pressure drop, and heat transfer are all closely coupled through the mixture density, the accurate prediction of void fraction is pivotal. Furthermore, while these passive systems boast a more reliable cooling option, the sensitivity of the system to two-phase flow instabilities is more likely than in forced convective cooling and can result in dangerous density wave or flow excursion conditions. Improvement of one-dimensional system codes and further development of three-dimensional code capability are all important efforts in addressing these issues. The proposed work is a critical step towards providing confidence through validation of the capability of thermal-hydraulic tools to meet the emerging challenges in reactor simulation. Historically, system analysis codes develop a stiffness over time due to fine-tuning of adjustable code parameters that result in a reluctance or inability to adopt improved models. Therefore, a rigorous effort should be invested in preparation for RELAP-7 to identify required modeling improvements of constitutive relations and validation of the code applicability to the wide range of possible reactor flow conditions. Part of this validation involves the uncertainty quantification which is a critical step in the Code Scaling, Applicability, and Uncertainty (CSAU) methodology and Best-Estimate Plus Uncertainty (BEPU) licensing method. BEPU allows recovery of safety margins and potentially optimizing plant performance, resulting in economic benefit for the reactor operators. Beginning with the collection of past experimental data covering prototypic LWR conditions, this document reviews three critical datasets available for the validation of system codes. The experimental methods, procedures, facilities, instrumentation, and supplementary information are consolidated to aid the validation efforts of this work and future projects. These three datasets (824 conditions) are chosen due to their wide acceptance and demonstrated application for code validation. Two experimental datasets from the Nuclear Power Engineering Corporation (NUPEC) span prototypic conditions of Boiling Water Reactors (BWR) and Pressurized Water Reactors (PWR). These datasets are labeled as the BFBT (BWR Full-size Fine-mesh Bundle Test) and the PSBT (PWR Subchannel and Bundle Test) datasets respectively. The third dataset is taken in an internally heated annulus that is scaled based on a BWR subassembly and includes two-phase flows important to post-blowdown analysis. Together, the datasets consider the full range of pressure and mass flux expected in LWR forced convective flows from steady-state to long term cooling conditions. Recognizing the deficiency of available data under low-pressure, low-flow, a new natural circulation dataset consisting of 107 conditions is obtained. With system pressures ranging between 145 and 950 kPa, natural circulation is driven by a heat flux up to 275 kW/m2 in a 3 m long heated region of the annulus test section which is followed by a 2 m long unheated region before entering a condenser. Stable natural circulation flow rates in this facility were measured to be between 182 and 590 kg/m2-s. Pressure and temperature are measured at ten locations in the facility. In addition, measurements of void fraction, interfacial area concentration, and gas velocity are obtained at five axial locations in the annulus test section. The RELAP5-3D code demonstrates good accuracy for void fraction predictions in BFBT and PSBT benchmarks, but shows relatively large discrepancy in the annulus experiments at lower pressures. In two-phase natural circulation, large oscillations are observed in many system parameters during the simulation of 19 experiment cases which were conducted at low pressures, and the model shows good predictive performance with the remaining 88 experiment cases. However, some discrepancies in void fraction and temperature predictions are found which may require improved constitutive modeling under these low-pressure, low-flow condition. The quantified uncertainty in void fraction predictions in three forced convection benchmarks due to uncertainties in boundary conditions and geometries are relatively small, except for a few cases in the annulus experiment. The output uncertainty of the natural circulation model is quantified by considering the measurement error in boundary conditions, and it is shown that the output uncertainty ranges are relatively small for the range of conditions considered. Through this work, a clear pathway is provided for validation of RELAP7 once the code is available.

42 ENGINEERING↗
DOE OSTI2021

Heat transfer and pseudo phase transition for low-Reynolds, mixed-convection channel flow in the supercritical thermodynamic regime

Fluids at supercritical thermodynamic conditions are inherently complex due to the large variations in thermodynamic and transport properties. Multiple recent numerical and experimental investigations indicate ongoing interest for these fluids, especially supercritical CO and supercritical water, for a variety of applications. The potential heat transfer benefits in this regime are directly influenced by the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. Limited studies have been conducted for low-Reynolds/intermediate-Rayleigh numbers mixed-convection channel flow of supercritical water, which is the focus of this study. To investigate the thermally driven hydrodynamic instabilities in this regime, we use a high-order fully-implicit numerical method. As expected, in this channel configuration and flow regime, a competition arises between the forced and the natural convection. Buoyancy, driven by the thermal gradient at the wall, pushes the fluid upwards in the channel, while inertial forces drive the fluid from inlet to outlet. We observe significant differences in the heat transfer coefficient between the top and bottom walls of the channel flow, due to natural convection, and unstable plumes. We further quantify the buoyancy forces along the channel flow through evaluation of the Richardson number. Simulations show that with a higher negative Richardson number, there are more unstable structures along the channel, resulting in enhanced mixing. Furthermore, we demonstrate the existence of gas-like, liquid-like, and the possibility of pseudo-two-phase regimes in the supercritical fluid. The pseudo-two-phase regime occurs when the range of fluid temperatures is crossing the pseudo-critical line, also known as the Widom line. The density of the fluid varies significantly near the Widom line, and further influences the natural convection and heat transfer in the flow.

42 ENGINEERING↗
DOE OSTI2025

A Decadal Hybrid GCM Simulation Using Deep‐Learning‐Based Cloud and Convection Parameterization Generalized to a Warm Climate

A critical challenge for machine‐learning (ML) parameterization in global climate models (GCMs) is to achieve stable, accurate simulations under climates not seen during training. Previous studies have demonstrated promising offline performance and year‐long online stability in aquaplanet simulations but have encountered difficulties in real geography and under climate warming. Here we report that a GCM with real geography configuration using neural‐network‐based cloud and convection parameterization, trained exclusively with present‐day climate data, successfully performs a stable, decade‐long simulation of a warm climate with +4 K sea surface temperature (SST). The neural network (NN) is based on Han et al. (2023, https://doi.org/10.1029/2022ms003508 ) with additional inputs. The simulation captures the global precipitation distribution, surface temperatures, vertical atmospheric structures, and extreme precipitation very well, closely matching simulations from both the superparameterized CAM (SPCAM) and the conventional CAM5 in the warm climate without accuracy degradation compared to those in the baseline climate. Moreover, it produces a climate response to +4 K SST in atmospheric thermodynamic states and circulations similar to those from SPCAM and CAM5. Prognostic ablation tests on NN input variables show that the NN without convective memory as input suffers from numerical instability, and the NN without considering radiative variables and land fraction as input, or with reduced training samples produce less accurate results. To our knowledge, this is the first time an ML parameterization successfully achieves online extrapolation to a warm climate without using additional warm‐climate data for training. It demonstrates the potential of ML‐driven parameterizations for credible long‐term climate projections.

Atmosphere model↗
DOE OSTI2023

Experimental investigation of power transient flow boiling

To advance the mechanistic understanding of power transients in flow boiling, experimental investigations are performed to capture the thermal response of four different kinds of cladding materials, including FeCrAls, to Fuchs reactivity initiated accident transients and linear ramp power transients. The transition boiling regime is phenomenologically reflected in power-transient flow boiling, which is different from steady-state boiling. In light of this, the critical heat flux that is usually featured with temperature overshooting does not indicate thermal safety margin of cladding materials as conservatively as the maximum heat flux, which is greater than critical heat flux in the power transient boiling curve. This is physically attributed to the thermal energy deposition on the cladding wall. In addition, the thermal energy released from the cladding wall results in different boiling heat transfer coefficient for the decreasing and increasing power stages respectively during the Fuchs power transient. The maximum heat flux difference gap, which appears due to changes in cladding materials and transient time scales, is appreciable under an intermediate heat convection regime. However, this difference gap is gradually weakened by the progressive increasing of mass flux and/or inlet subcooling due to the enhanced dominance of heat convection over heat conduction. Finally, the small difference gap of maximum heat flux is found to be insignificant under the weak heat convection because of power transient induced annular flow instability.

42 ENGINEERING↗