Search NASA⌕ Search

SEARCH · Search NASA

Results for “THERMAL 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 OSTI2021

MOSCATO Solver Development and Integration Plan

During FY21, we conducted ongoing development work for the MOSCATO (Molten Salt Chemistry and Transport) solver. The code development work primarily consisted of transitioning capabilities from the original version of the solver, which was written in OpenFOAM, into Nek5000. In doing so, a fast, highly parallelizable solver was created that is capable of complex chemistry and corrosion simulations for engineering-scale molten salt systems. The Nek5000 version of MOSCATO is now fully featured and capable of higher-fidelity simulations than were previously possible. Demonstration cases including a thermal convection loop have been simulated to test these new capabilities. Although capable of large-scale simulations, MOSCATO is not well-suited to parametric studies of complete reactor geometries. These types of simulations are instead better handled by reduced-order modeling codes such as ORNL’s Mole code. Reduced-order simulation tools like Mole, however, are dependent on high fidelity correlations to account for complex, coupled three-dimensional phenomena that they do not directly simulate. Tools such as MOSCATO must therefore be used to create these correlations, as suitable empirical relationships are not available for most molten salt systems. Toward that end, we used the Nek-derived version of MOSCATO to create new mass transfer correlations for three relevant cases including tubular, tube bank, and subchannel geometries. These new correlations are more accurate than any existing ones and can be readily integrated into any reduced-order modeling tools that are targeting full-scale MSR simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Corrosion of 316H Stainless Steel in Flowing FLiBe (LiF-BeF 2 )

This milestone was originally envisioned for completion in FY21 and was delayed due to the COVID-19 response and difficulties in fabricating/obtaining fluoride salts. A monometallic type 316H stainless steel thermal convection loop (TCL) was operated for 1000 h with flowing LiFBeF 2 (i.e. FLiBe) salt and a peak temperature of 650°C. In general, the attack was minimal. However, classic mass transfer was not clearly observed with mass losses in both the hot and cold legs. There was an issue with cleaning the loop with water and specimens may have oxidized, which was not an issue with the previous FLiNaK TCL.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Volumetric Rendering on Wavelet-Based Adaptive Grid

Numerical modeling of physical phenomena frequently involves processes across a wide range of spatial and temporal scales. In the last two decades, the advancements in wavelet-based numerical methodologies to solve partial differential equations, combined with the unique properties of wavelet analysis to resolve localized structures of the solution on dynamically adaptive computational meshes, make it feasible to perform large-scale numerical simulations of a variety of physical systems on a dynamically adaptive computational mesh that changes both in space and time. Volumetric visualization of the solution is an essential part of scientific computing, yet the existing volumetric visualization techniques do not take full advantage of multi-resolution wavelet analysis and are not fully tailored for visualization of a compressed solution on the wavelet-based adaptive computational mesh. Our objective is to explore the alternatives for the visualization of time-dependent data on space-time varying adaptive mesh using volume rendering while capitalizing on the available sparse data representation. Two alternative formulations are explored. The first one is based on volumetric ray casting of multi-scale datasets in wavelet space. Rather than working with the wavelets at the finest possible resolution, a partial inverse wavelet transform is performed as a preprocessing step to obtain scaling functions on a uniform grid at a user-prescribed resolution. As a result, a solution in physical space is represented by a superposition of scaling functions on a coarse regular grid and wavelets on an adaptive mesh. An efficient and accurate ray casting algorithm is based just on these coarse scaling functions. Additional details are added during the ray tracing by taking an appropriate number of wavelets into account based on support overlap with the interpolation point, wavelet coefficient magnitude, and other characteristics, such as opacity accumulation (front to back ordering) and deviation from frontal viewing direction. The second approach is based on complementing of wavelet-based adaptive mesh to the traditional Adaptive Mesh Refinement (AMR) mesh. Both algorithms are illustrated and compared to the existing volume visualization software for Rayleigh-Benard thermal convection and electron density data sets in terms of rendering time and visual quality for different data compression of both wavelet-based and AMR adaptive meshes.

Vezolainen, Alexei V.↗
DOE OSTI2024

Computational Methods for Multi-Physics Simulation of Melting in Steelmaking

Iron and steel production accounts for approximately 8% of global carbon dioxide (CO) emissions. Pathways to decarbonize include replacing fossil fuels in iron ore reduction and electrifying other steelmaking processes. Iron pellets produced by hydrogen, called Hydrogen Direct Reduced Iron (HDRI), have property differences from those produced using conventional DRI processes. These differences may impact melting in electric arc furnaces (EAF) and other downstream processes. The physical properties of iron pellets vary significantly with temperature during heating, complicating predictions of their behavior. In this project, we seek to develop an integrated simulation, including the fluid flow and convective thermal transport around the pellet particle. We also examine conduction and phase changes within the particle as they impact the melting process. We use adaptive mesh refinement (AMR) to resolve both the changing size of the particle and the complex physics of the interaction between the pellet and the surrounding fluid. We base our simulations on the AMReX-incflo module, which allows large-scale Navier-Stokes simulation while resolving the changing particle size during melting. As we advance our numerical tools, we anticipate an improved understanding of the dynamics of HDRI melting, which will, in turn, accelerate the adoption of low-carbon technologies in the steelmaking industry.

AMReX↗
DOE OSTI2021

Nondimensional convection numbers modeling thermally stratified storage tanks: Richardson's number and hot-water tanks

Thermally stratified storage tank studies have spanned over 50 years to increase the thermal storage efficiency and accurate prediction of the outlet temperature particularly for solar applications. The studies have reviewed and modeled the jet and plume flow phenomena inside the tank due to the inlet mixing and stratification level. Kelvin–Helmholtz and Rayleigh–Taylor instabilities are the major drivers of the mixing in these tanks. Momentum jets deflecting off walls at the bottom of the tank also create significant mixing. Reviewing Richardson models shows that the categorization was based on the range of Reynolds numbers at the inlet. Unfortunately, the use of superficial velocity in calculating the Richardson number results in critical values in the literature ranging from below 0.25 to 100. The most used length scale associated with these flows is an inertial scale based on the tank height or diameter although the mixing can occur at a relatively smaller scale. The various inlet devices and a large span of flow rates experienced in thermally stratified storage tanks requisite the use of the Reynolds number in combination with a convection number for accurate one dimensional models that predict performance over the long-term. The evaluation of peak shifting of electric loads leveraging renewable sources for applications, including residential hot-water tanks, commercial water tanks, and large-scale chilled water storage tanks, require these models. Finally, this paper is focused on establishing the significance of the convection numbers in conjunction with the Reynolds number for modeling the thermal stratification in storage tanks.

14 SOLAR ENERGY↗
DOE OSTI2025

Thermomechanics coupling to Monte Carlo particle transport on unstructured mesh geometries using Cardinal

Geometry deformation due to thermal expansion influences neutron transport in many systems. Studying this phenomenon involves coupling models for neutronics, thermal hydraulics, and solid mechanics. To enable high fidelity modeling of these coupled physics, new capabilities were introduced in Cardinal, coupling OpenMC Monte Carlo particle transport models with MOOSE thermomechanical physics on unstructured moving-mesh geometries. In this work, we present a fully open-source capability leveraging on-the-fly mesh skinning to automatically regenerate OpenMC geometry, which allows multiphysics feedback from temperature, density, and geometry changes. The new capability is verified using an analytic benchmark slab problem, which couples S 2 neutron transport with thermal conduction, convective boundary conditions, Doppler-broadened cross sections, and nonlinear thermal expansion effects along the heated slab. Cardinal reproduces the analytic solutions for the neutron flux, heating, k eff , and temperature with demonstrated convergence in various error terms including mesh resolution and cross section temperature library spacing. For the nominal benchmark conditions and with a fine mesh, maximum relative errors for neutron flux, temperature, and heating are lower than 1%, while errors in integral quantities such as k eff and slab length are within 1 pcm and 48 µm, respectively. This work (i) presents a new numerical approach to thermomechanics coupling with OpenMC models, (ii) is the first (to our knowledge) to utilize a mechanical partial differential equation (PDE) solution to solve the (Griesheimer and Kooreman, 2022) analytic benchmark, and (iii) develops this verified capability within an open-source package.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Formation of a Large Cold Groundwater Mantle Helium Anomaly and High Temperature Geothermal Resources in Response to Bimodal Magmatism Near Roosevelt Hot Springs and Utah FORGE, Milford Valley, Southwest Utah

Abstract A large mantle helium anomaly and separate domains of high heat flow are the predominant manifestations of bimodal magmatic activity in the Milford valley. The mantle helium anomaly (1.9–2.6 R/Ra) covers 270 km 2 and is subdivided into two separated domains: a cold shallow groundwater regime and high temperature hydrothermal activity. The zone of anomalous heat flow covers >100 km 2 and is also subdivided into two adjacent domains, comprising hydrothermal activity at Roosevelt Hot Springs (RHS) (3–7 W/m 2 ) and conductive heat flow (100–180 mW/m 2 ). While the transfer of heat and mantle helium at RHS are coupled, heat and helium transfer are decoupled in the adjacent cold groundwater regime to the west. Both the mantle helium and geothermal anomalies are attributed to recent mafic‐felsic magmatic intrusions of >400 km 3 , however, the absence of volcanic eruptions <500,000 years indicates magmas stall before rising to shallow crustal level <10 km depth. Deep level magmatism produces a felsic composition melt, which is inferred to be responsible for the widespread and near uniform range of diluted mantle helium values. A thick and impermeable mass of crystalline granitic basement rock at the mid‐crustal level divides the ascent of mantle helium into separate flow paths. It may also impede the rise of buoyant magma trapping thermal energy that facilitates partial melting, slow cooling, and development of a thick thermal aureole. Partitioning of convective and conductive thermal regimes and independent flow paths supplying deeply derived helium characterize the development of a large long‐lived magma‐related geothermal system.

Geochemistry & Geophysics↗
DOE OSTI2021

Simulation of flame acceleration and deflagration-to-detonation transition with heat transfer in HE product and fractured HE

A simple heat transfer model has been developed for high explosives. For the gaseous high explosive (HE) product, the Smargorinsky eddy viscosity model has been adopted for estimating the turbulent thermal conductivity. Convective heat transfer in the fractured solid HE is modeled as the effective thermal conduction in the porous media, which is equivalent to the enthalpy flux carried by the velocity determined by Darcy’s law with Forchheimer term.Other features of the model include a multi-reaction chain, equation-of-state for each species involved, and a two-phase hotspot model. Test calculations are performed with various initial disturbances. Onset and progress of flame acceleration (FA) and chemical reactions are strongly affected by the heat transfer and initial disturbance. Staged development of FA also affirms the necessity of modeling based on a multi-reaction chain. When conditions warrant,e.g., stronger initial disturbance, FA evolves to deflagration-to-detonation transition (DDT),demonstrating automatic capture of the DDT.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
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 OSTI2022

Analytical benchmark solution for 1-D neutron transport coupled with thermal conduction and material expansion

In this work we present fully analytical solutions for a class of finite, homogeneous, 1-D slab benchmark problems with nonlinear temperature feedback effects. The proposed class of benchmarks include multiplicative 1-D neutron transport (limited to quasistatic S{sub 2} with μ = ±1) coupled with thermal conduction, convection, Doppler broadening, and expansion effects along the length of the slab. This class of benchmark models, along with the corresponding analytical solutions, are valuable for validating multiphysics analysis frameworks that support coupled neutronics/thermal/structural calculations. Analytical solutions for the benchmarks are obtained by introducing an ansatz that the equilibrium flux and temperature distributions in the slab have the same shape. Specific values for the total microscopic cross section, σ{sub t,0}, and conductive heat transfer coefficient, h, that satisfy the assumed ansatz are then determined. A discussion of the procedure for generating benchmark models and analytical solutions is provided, along with numerical results for an example set of model parameters and thoughts on practical applications of the benchmarks for multiphysics code validation. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Collaborative Proposal: Improving understanding of the internal structure and dynamics of deep convection using ARM observations and large eddy simulations

Recent observational and large eddy simulation (LES) modeling studies have nearly unanimously supported the view of deep cumulus convection being composed of a series of quasi-spherical bubbles of buoyant air, known as moist thermals. Despite the prevalence of moist thermals in deep convection, a comprehensive theory for the dynamics of these structures is lacking. Most current conceptual models for cumulus convection are based on canonical scaling theories for dry thermals or plumes; however, there is considerable evidence that the behavior of moist thermals differs markedly from these theories. Furthermore, the theoretical basis for most cumulus parameterizations originates from the plume conceptual model, and therefore these parameterizations are inconsistent with the real structure of moist convection. Motivated by the aforementioned knowledge gaps, this “end-to-end” research effort use theory, observations, numerical simulations, and direct improvements to the Zhang-McFarlane (ZM) convection scheme in the global climate Community Atmosphere Model (CAM) to address the following research questions: What key environmental parameters determine whether or not shallow convection will transition into deep convection, in the context of thermal-like updrafts? What factors regulate the size of thermals within cumulus updrafts? How does vertical wind shear influence thermal behavior, and as a consequence, vertical velocity and mass flux profiles and the shallow-to-deep convective transition? What are the critical processes that determine updraft vertical velocities and their connection to the vertical mass flux profile for thermal-like updrafts? Idealized LES modeling will be used in conjunction with theoretical models for the core properties of thermal-like updrafts to better understand key processes that regulate thermal ascent rates and entrainment properties. Thermal-tracking procedures will be used to characterize the behavior of thermals within the LES, and recently developed direct measures of entrainment and detrainment will be used to quantify entrainment/detrainment rates. Building from these results, we will analyze the structure of moist thermals from hemispheric range-height indicator scans taken during the Atmospheric Radiation Measurement Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign, and from “real case” LES of CACTI events. This combined modeling and observational analysis will provide essential validation for the existing body of research on moist thermal dynamics, which is based primarily on modeling studies. With the insight gained from the aforementioned activities, we will modify the Zhang-McFarlane convection scheme to improve its representation of updraft vertical velocity and entrainment rate profiles. These process-level changes will be tested in the Community Atmosphere Model to assess the impact on global climate simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Aerosol-deep convection interaction based on joint cell-thermal tracking in Large Eddy Simulations during the TRACER campaign

In cumulus clouds, aerosol concentrations control cloud droplet concentrations, modifying cloud radiative properties, precipitation processes, and cloud electrification. However, mechanisms of aerosol-deep convection interactions are not well understood due to complex cloud dynamics and microphysics. We investigate the interaction of aerosols with isolated deep convection using Large Eddy Simulations of two cases during the TRacking Aerosol Convection interactions ExpeRiment (TRACER) near Houston, Texas, using a joint cell-thermal tracking algorithm. Cumulus thermals are droplet generators, since supersaturation and droplet nucleation coincide with thermal centers, where the strongest updrafts occur. Primary ice crystal formation does not take place inside thermals, but at layers where previous thermals detrained moisture. As subsequent thermals containing supercooled droplets penetrate these layers, hail and graupel form at or near these thermals. Higher aerosol concentrations result in higher droplet concentrations that suppress drizzle, delay warm rain processes, and transport more moisture aloft. This increases snow and ice amount, as well as graupel and hail, leading to more lightning. Polluted thermals initiate at slightly higher altitudes, and are slightly larger and faster, suggesting a weak invigoration. We also find more thermals per cell, but fewer isolated cells, since convection is more aggregated and intense, especially near the end of the 24 h simulation. Non-linear mesoscale feedback likely triggered by temperature and moisture responses to aerosol-thermal interactions causes the aggregation. Time-lagged aerosol-reinitialization experiments show that the mesoscale response is the predominant forcing for the invigoration. These changes happen within one day, on a smaller scale than previously suggested.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Collaborative Proposal: Improving Understanding of the Internal Structure and Dynamics of Deep Convection Using ARM Observations and Large Eddy Simulations (Final progress report)

This is the final technical report for the DOE Atmospheric System Research funded project entitled "Collaborative Proposal: Improving Understanding of the Internal Structure and Dynamics of Deep Convection Using ARM Observations and Large Eddy Simulations". This project addressed several science questions related to the structure and behavior of deep convective clouds in the atmosphere. It applied knowledge gained from this work to improve the deep cumulus convection parameterization in the Community Earth System Model (CESM). Specific outcomes related to the project include an improved understanding of: (1) how convective cloud environments regulate the size of thermals reaching the upper troposphere, and hence the shallow-to-deep convective transition; (2) the role of dynamic pressure forcing in thermal and convective cloud evolution; and (3) the role of vertical wind shear on the behavior of moist thermals and the transition from shallow to deep convection. Although the primary motivation of this project was to improve basic understanding of thermal dynamics in atmospheric convection, an important component of this work was using this improved understanding to modify the representation of entrainment and vertical velocity profiles in the Zhang-McFarlane cumulus parameterization used in CESM. This work also had a strong observational component, with analyses of Atmospheric Radiation Measurement (ARM) observations to characterize thermal characteristics and behavior. Observations collected during the ARM-supported CACTI field project were used with high-resolution simulations to explore environmental controls on deep convection initiation. Data from the ARM-supported MC3E field project were also used to evaluate the representation of updraft vertical velocities in the Zhang-McFarlane cumulus parameterization. This project directly supported the training of two post-doctoral researchers and led to 16 peer-reviewed published papers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Dynamic Behavior of Oval-Twisted Helical Tube Heat Exchanger: Numerical Study with RELAP5-3D

Convective heat transfer characteristics and theoretical thermal stress behaviors are numerically calculated using RELAP5-3D for the helical-coiled once-through steam generator (H-OTSG) and the novel heat exchanger design known as the oval-twisted helically coiled heat exchanger (OTHCHX) under (1) fluctuating wall temperature conditions, (2) square-wave pulsating flow conditions, and (3) the combined effects of fluctuating wall temperature and square-wave pulsating flow conditions. Heat transfer coefficient models for the H-OTSG and OTHCHX were developed based on existing data and implemented into RELAP5-3D, successfully capturing the N⁢uavg behavior within 8% to 10% of the reported data. Under fluctuating wall temperature conditions, the OTHCHX displayed higher N⁢u avg behavior than the H-OTSG. As 𝑓 increased, the $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ decreased. The $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ was higher for the OTHCHX than for the H-OTSG under fluctuating wall temperature conditions. Under pulsating flow conditions, the H-OTSG and OTHCHX displayed much higher 𝑁⁢𝑢 𝑎𝑣𝑔 than under constant flow conditions. The H-OTSG displayed a higher $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ over the OTHCHX. Under combined fluctuating wall temperature and pulsating flow conditions, the augmented heat transfer behavior from the pulsating flow was counteracted by the wall temperature fluctuations, producing slightly higher 𝑁⁢𝑢 𝑎𝑣𝑔 over constant wall temperature, constant flow conditions, but much lower than only constant pulsating flow under constant wall temperature conditions. The effects of simultaneous wall temperature fluctuations and square-wave pulsating flow caused higher $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ than that of only wall temperature fluctuations or pulsating flow. As the Reynolds number (Re) increased, $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ increased. However, when 𝑓=𝑓$_{\dot{m}}$, the $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ showed decreasing values as Re increased. In conclusion, the results indicate that thermal-fluid resonance can help mitigate thermal stresses.

Thermal stress↗
DOE OSTI2020

Numerical Simulations of Orographic Convection across Multiple Gray Zones

Idealized simulations are used to determine the sensitivity of moist orographic convection to horizontal grid spacing Δ h . In simulated mechanically (MECH) and thermally (THERM) forced convection over an isolated ridge, Δ h is varied systematically over both the deep-convection (Δ h ~ 10 –1 km) and turbulence (Δ h ~ 1 km–100 m) gray zones. To aid physical interpretation, a new parcel-based bulk entrainment/detrainment diagnosis for horizontally heterogeneous flows is developed. Within the deep-convection gray zone, the Δ h sensitivity is dominated by differences in parameterized versus explicit convection; the former initiates convection too far upstream of the ridge (MECH) and too early in the diurnal heating cycle (THERM). These errors stem in part from a large underprediction of parameterized entrainment and detrainment. Within the turbulence gray zone, sensitivities to Δh arise from the representation of both subcloud- and cloud-layer turbulence. As Δ h is decreased, MECH exhibits stronger cloud-layer entrainment to enhance the convective mass flux M co , while THERM exhibits stronger detrainment to suppress Mco and delay convection initiation. The latter is reinforced by increased subcloud turbulence at smaller Δ h , which leads to drying and diffusion of the central updraft responsible for initiating moist convection. Numerical convergence to a robust solution occurs only in THERM, which develops a fully turbulent flow with a resolved inertial subrange (for Δ h ≤ 250 m). In MECH, by contrast, turbulent transition occurs within the orographic cloud, the details of which depend on both physical location and Δ h .

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

Collective effect of thermal plumes on temperature fluctuations in a closed Rayleigh–Benard convection cell

Here, we report a systematic study of the collective effect of thermal plumes on the probability density function (PDF) P(δT) of temperature fluctuations δT(t) in turbulent Rayleigh-Benard convection. By decomposing δT(t) into four basic fluctuation modes associated with single and multiple warm and cold plumes and a turbulent background, we derive an analytic form of P(δT) based on the convolutions of the five independent modes. To test the derived form of P(δT) in the multiple-plume regions, where the thermal plumes are heavily populated, we conduct time series measurements of temperature fluctuations in two convection cells; one is a vertical thin disk and the other is an upright cylinder of aspect ratio unity. For a given normalized position in most regions of the convection cell, all of the measured PDFs P(δT) for different Rayleigh numbers fall onto a single master curve, once δT is normalized by its rms value σ T . It is found that the measured P(δT / σT ) at different locations along the symmetric horizontal and vertical axes of the convection cells can all be well described by the derived form of P(δT / σT ). The fitted values of the parameters associated with the number of plumes in multiple plume clusters and their relative strengths and degrees of intermittency are closely linked to the spatial distribution of thermal plumes and local dynamics of the large-scale circulation in a closed convection cell. Our work thus provides a unified theoretical approach for understanding scalar PDFs in a turbulent field, which is very useful not only for the present study but also for the study of many turbulent mixing problems of practical interest.

42 ENGINEERING↗
DOE OSTI2024

Mechanical and Thermal Forcing for Upslope Flows and Cumulus Convection over the Sierras de Córdoba

Abstract The upslope flow processes affecting the vertical extent of orographic cumulus convection are examined using observations from the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. Specifically, clear air returns from the U.S. Department of Energy (DOE) second-generation C-band scanning Atmospheric Radiation Measurement (ARM) precipitation radar (CSAPR2) are used to characterize the structure and variability of the ridge-normal (i.e., up/downslope) flow components, which transport mass to the crest of Argentina’s Sierras de Córdoba and contribute to convective initiation. Data are compiled for the entire CACTI period (October–April), including days with clear skies, shallow cumuli, cumulus congestus, and deep convection. To examine shared variability among >70 000 radar scans, we use (i) a principal component analysis (PCA) to isolate modes of variability in the upslope flow and (ii) composite analysis based on convective outcomes, determined from GOES-16 satellite observations. These data are contextualized with observed surface sensible heat fluxes, thermodynamic profiles, and synoptic-scale analysis. Results indicate distinct thermally and mechanically forced upslope flow modes, modulated by diurnal heating and synoptic-scale variations, respectively. In some instances, there is a superposition of thermal and mechanical forcing, yielding either deeper or shallower upslope flow. The composite analyses based on satellite data show that successively deeper convective outcomes are associated with successively deeper upslope flow layers that more readily transport mass to the ridge crest in conjunction with lower lifting condensation levels, facilitating convective initiation. These results help to isolate the forcing mechanisms for orographic convection and thus provide a foundation for parameterizing orographic convective processes in coarse resolution models.

Meteorology & Atmospheric Sciences↗