Search NASA⌕ Search

SEARCH · Search NASA

Results for “circulation model”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Measuring scale-dependent shape anisotropy by coarse-graining: Application to inhomogeneous Rayleigh-Taylor turbulence

Here, we generalize the “filtering spectrum” to probe scales along different directions by spatial coarse-graining. This multidimensional filtering spectrum quantifies the spectral content of flows that are not necessarily homogeneous. From multidimensional spectral information, we propose a simple metric for shape anisotropy at various scales. The method is applied to simulations of two-dimensional (2D) and 3D Rayleigh-Taylor (RT) turbulence, which is inhomogeneous and anisotropic. We show that 3D RT has clear shape anisotropy at large scales with approximately 4:3 vertical to horizontal aspect ratio but tends toward isotropy at small scales as expected. In sharp contrast, we find that RT in 2D simulations, which are still the main modeling framework for many applications, is isotropic at large scales and its shape anisotropy increases at smaller scales where structures tend to be horizontally elongated. While this may be surprising, it is consistent with recent results that large-scale isotropy in 2D RT is due to the generation of a large-scale overturning circulation via an upscale cascade, while small-scale anisotropy is due to the stable stratification resultant from such overturning and the inefficient mixing in 2D.

42 ENGINEERING↗
DOE OSTI2020

Thermally Anisotropic Composites for Heat Redirection and Thermal Management in Building Envelopes

In building envelope systems, thermal management is important from energy conservation and thermal comfort perspectives, and is typically done via insulation, thermal mass, solar control and shading, etc. Thermal anisotropy is commonly utilized in electronics for heat dissipation. This article investigates the feasibility of utilizing thermally anisotropic composites (TACs) to reduce unwanted heat gains and losses through the building envelope by redirecting the heat transfer to or from appropriate heat sinks and sources.Numerical simulations of building envelope components using COMSOL Multiphysics and whole building simulations using EnergyPlus (E+) were performed to estimate the energy savings potential of TACs used in external walls of residential buildings. Application of TACs was compared with baseline walls with cavity insulation as well as with and without exterior continuous insulation (CI). The COMSOL simulations indicated 75% or higher reductions in heat gains when using a TAC coupled with a heat sink compared to an exterior CI case. In terms of energy savings in a hot climate, E+ models estimated cooling energy savings of 19% with a "TAC + heat sink" compared to a baseline wall with only cavity insulation. Further, a small-scale test sample with a TAC-heat sink system was built and tested in a multi-transducer heat flow meter (HFM). The heat sink was a copper tube circulating cold water. The experiments verified the re-direction of heat to a sink and significant reduction in heat transfer through the test specimen. The experiments showed that the heat redirection only occurred in the presence of both thermal anisotropy and water circulation; if either element was missing, the system behaved like a traditional envelope component.

Biswas, Kaushik↗
DOE OSTI2020

Pacific decadal oscillation remotely forced by the equatorial Pacific and the Atlantic Oceans

Abstract The Pacific Decadal Oscillation (PDO), the leading mode of Pacific decadal sea surface temperature variability, arises mainly from combinations of regional air-sea interaction within the North Pacific Ocean and remote forcing, such as from the tropical Pacific and the Atlantic. Because of such a combination of mechanisms, a question remains as to how much PDO variability originates from these regions. To better understand PDO variability, the equatorial Pacific and the Atlantic impacts on the PDO are examined using several 3-dimensional partial ocean data assimilation experiments conducted with two global climate models: the CESM1.0 and MIROC3.2m. In these partial assimilation experiments, the climate models are constrained by observed temperature and salinity anomalies, one solely in the Atlantic basin and the other solely in the equatorial Pacific basin, but are allowed to evolve freely in other regions. These experiments demonstrate that, in addition to the tropical Pacific’s role in driving PDO variability, the Atlantic can affect PDO variability by modulating the tropical Pacific climate through two proposed processes. One is the equatorial pathway, in which tropical Atlantic sea surface temperature (SST) variability causes an El Niño-like SST response in the equatorial Pacific through the reorganization of the global Walker circulation. The other is the north tropical pathway, where low-frequency SST variability associated with the Atlantic Multidecadal Oscillation induces a Matsuno-Gill type atmospheric response in the tropical Atlantic-Pacific sectors north of the equator. These results provide a quantitative assessment suggesting that 12–29% of PDO variance originates from the Atlantic Ocean and 40–44% from the tropical Pacific. The remaining 27–48% of the variance is inferred to arise from other processes such as regional ocean-atmosphere interactions in the North Pacific and possibly teleconnections from the Indian Ocean.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Winter Warming in North America Induced by Urbanization in China

Rapid urbanization in China in the past several decades has been shown to affect surface air temperature (SAT) locally, but its global impact is unclear. Using a global climate model, we have investigated the potential impact of urbanization in China on SAT at the global scale. The expansion of urban areas in China from 1985 to 2017 leads to winter surface air warming by as much as 0.8°C in North America (NA), largely caused by warm temperature advection due to global-scale atmospheric circulation changes. The East Asian jet stream moves northward in eastern and northeastern China where the expansion of urban areas has increased SAT and sensible heat flux, forcing change of the atmospheric stationary waves and so the structure of the jet streams.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Projected End-of-Century Changes in the South American Monsoon in the CESM Large Ensemble

Projected changes in the South American monsoon system by the end of the twenty-first century are analyzed using the Community Earth System Model Large Ensemble (CESM-LENS). The wet season is shorter in LENS when compared to observations, with the mean onset occurring 19 days later and the mean retreat date 21 days earlier in the season. Despite a precipitation bias, the seasonality of rainfall over South America is reproduced in LENS, as well as the main circulation features associated with the development of the South American monsoon. Both the onset and retreat of the wet season over South America are delayed in the future compared to current climate by 3 and 7 days, respectively, with a slightly longer wet season. Central and southeastern Brazil are projected to get wetter as a result of moisture convergence from the strengthening of the South Atlantic low-level jet and a weaker South Atlantic subtropical high. The Amazon is projected to get drier by the end of the century, negatively affecting rain forest productivity. During the wet season, an increase in the frequency and intensity of extreme precipitation events is found over most of South America, and especially over northeastern and southern Brazil and La Plata. Meanwhile, during the dry season an increase in the maximum number of consecutive dry days is found over northeastern Brazil and the northern Amazon.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Biphasic decay of intact SHIV genomes following initiation of antiretroviral therapy complicates analysis of interventions targeting the reservoir

The latent reservoir for HIV-1 in resting CD4 + T cells persists despite antiretroviral therapy (ART) and precludes cure. Reservoir-targeting interventions are evaluated in ART-treated macaques infected with simian immunodeficiency virus (SIV) or simian-human immunodeficiency virus (SHIV). Efficacy is determined by reservoir measurements before and after the intervention. However, most proviruses persisting in the setting of ART are defective. In addition, intact HIV-1 and SIV genomes undergo complex, multiphasic decay observable when new infection events are blocked by ART. Intervention-induced elimination of latently infected cells must be distinguished from natural decay. Here, we address these issues for SHIV. We describe an intact proviral DNA assay that allows digital counting of SHIV genomes lacking common fatal defects. We show that intact SHIV genomes in circulating CD4 + T cells undergo biphasic decay during the first year of ART, with a rapid first phase (t 1/2 = 30.1 d) and a slower second phase (t 1/2 = 8.1 mo) that is still more rapid that the slow decay observed in people with HIV-1 on long-term ART (t 1/2 = 3.7 y). In SHIV models, most interventions are tested during 2nd phase decay. Natural 2nd phase decay must be considered in evaluating interventions as most infected cells present at this time do not become part of the stable reservoir. In addition, for interventions tested during 2nd phase decay, a caveat is that the intervention may not be equally effective in people with HIV on long-term ART whose reservoirs are dominated by latently infected cells with a slower decay rate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

First operation of the KATRIN experiment with tritium

The determination of the neutrino mass is one of the major challenges in astroparticle physics today. Direct neutrino mass experiments, based solely on the kinematics of β -decay, provide a largely model-independent probe to the neutrino mass scale. The Karlsruhe Tritium Neutrino (KATRIN) experiment is designed to directly measure the effective electron antineutrino mass with a sensitivity of 0.2 eV (90 % CL). In this work we report on the first operation of KATRIN with tritium which took place in 2018. During this commissioning phase of the tritium circulation system, excellent agreement of the theoretical prediction with the recorded spectra was found and stable conditions over a time period of 13 days could be established. These results are an essential prerequisite for the subsequent neutrino mass measurements with KATRIN in 2019.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Assessing Global Ocean Circulation Sensitivity to Changes in Antarctic Ice Discharge

The principal goal of this research is to assess how projected increases in meltwater discharge and iceberg calving from the Antarctic Ice Sheet will impact climate in the next 250 years. In short, the sensitivity of the global ocean circulation and climate, and in particular the Atlantic Meridional Overturning Circulation (AMOC), to large increases in iceberg calving and meltwater discharges to the Southern Ocean are poorly understood. The requirement to investigate this topic is heightened by growing evidence that the West Antarctic Ice Sheet (WAIS) is vulnerable to rapid retreat and collapse through marine-ice-sheet-instability (MISI) on multidecadal-to-centennial timescales. Future model simulations of the AIS show the potential for rapid ice sheet retreat in the next 50 – 300 years, with such a collapse associated with the discharge of enormous volumes of ice and meltwater to the Southern Ocean. This project uses a combination of an ice sheet/shelf model and a fully coupled (ocean-atmosphere-seaice) model (CESM 1.2) to explore the sensitive of global climate to future changes in runoff from the Antarctic continent under future greenhouse gas emissions scenarios IPCC RCP4.5 and 8.5. Overall, accounting for Antarctic meltwater reduces the projected global mean anthropogenic greenhouse warming by up to 2 deg C. during peak WAIS collapse and by ~1 deg. for most of the simulation. In addition, the increased meltwater contribution raises sub-surface ocean temperatures at the Antarctic ice sheet margin by more than 1ºC, such that the ice shelf basal melt rates increase beyond those currently projected. In contrast, the model predicts a strong cooling of both the surface air and ocean temperatures over the Southern Ocean that results in a massive expansion and thickening of the perennial sea-ice. The slow-down in anthropogenic warming also results in the onset of ice-free conditions in the Arctic being delayed in the model by several decades, and the projected slowdown in the AMOC to reduce.

58 GEOSCIENCES↗
DOE OSTI2026

Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior (Final Technical Report)

Enhanced Geothermal Systems (EGS) produce clean energy by circulating fluid through hot rock deep underground and bringing that heat to the surface to generate electricity. For this process to work reliably, fluids must be able to move efficiently through networks of natural or engineered fractures in the rock. Enhancing and maintaining subsurface permeability over time is essential for sustainable energy production. However, fluid injection changes the underground temperature, pressure, rock stress, and chemistry, which can alter permeability and sometimes trigger earthquakes. Predicting these interconnected processes remains a key challenge. To address this, we combined high-temperature laboratory experiments with high-fidelity simulations to better understand how fractures in geothermal reservoirs evolve over time. Our experiments measured how fractures respond to stress, slip, slip rate, and chemical reactions under geothermal conditions. These data were integrated into coupled thermal-hydrological-mechanical-chemical and earthquake (THMC+E) models tailored to the Utah FORGE site. The validated modeling framework improves predictions of reservoir performance and seismic response and helps guide operational decisions. This work reduces technical risk and strengthens the scientific foundation needed to make geothermal energy a reliable and scalable clean energy resource.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Determining Spatial Scales of Soil Moisture—Cloud Coupling Pathways Using Semi-Idealized Simulations

The spatial scales of land-atmosphere interactions over heterogeneous, realistic soil moisture patterns are critical to improve our understanding of land-atmosphere interactions. However, these scales are poorly quantified. We use high-resolution numerical experiments to quantify the spatial scales for a selected day during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) campaign in Oklahoma. A semi-idealized approach (removing background winds and surface topography while retaining the real-world data for other boundary conditions) enables examination of the local influence from the surface. Simulations that apply uniform surface, variable land cover with uniform soil moisture, and variable land cover with variable soil moisture are then compared. Land cover variations, especially water and urban surfaces, create two variance peaks in the surface sensible heat flux in the scales >10 km (meso-β) and <6 km (meso-γ). Both peaks are enhanced by the observation-based soil moisture distribution that has most of its variance in the meso-β scale. Atmospheric response in the meso-β scale is hydrostatic with weak vertical motion. The land cover variations induce secondary circulations in the scales of 2.5–6 km. The additional soil moisture variability intensifies extreme circulations over dry soil patches and sharp sensible heat flux gradients, slightly extending the characteristic scales to 2.5–9 km. While their areal coverage is small, soil moisture variations on such “hotspots” lead to substantially larger convective cloud structure and domain-mean rainfall, calling for a consideration of their non-linear effect in large-scale models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Vegetation forcing modulates global land monsoon and water resources in a CO 2 -enriched climate

The global monsoon is characterised by transitions between pronounced dry and wet seasons, affecting food security for two-thirds of the world’s population. Rising atmospheric CO 2 influences the terrestrial hydrological cycle through climate-radiative and vegetation-physiological forcings. How these two forcings affect the seasonal intensity and characteristics of monsoonal precipitation and runoff is poorly understood. Here we use four Earth System Models to show that in a CO 2 -enriched climate, radiative forcing changes drive annual precipitation increases for most monsoon regions. Further, vegetation feedbacks substantially affect annual precipitation in North and South America and Australia monsoon regions. In the dry season, runoff increases over most monsoon regions, due to stomatal closure-driven evapotranspiration reductions and associated atmospheric circulation change. Our results imply that flood risks may amplify in the wet season. However, the lengthening of the monsoon rainfall season and reduced evapotranspiration will shorten the water resources scarcity period for most monsoon regions.

58 GEOSCIENCES↗
DOE OSTI2022

Concomitant immunity to M. tuberculosis infection

Some persistent infections provide a level of immunity that protects against reinfection with the same pathogen, a process referred to as concomitant immunity. To explore the phenomenon of concomitant immunity during Mycobacterium tuberculosis infection, we utilized HostSim, a previously published virtual host model of the immune response following Mtb infection. By simulating reinfection scenarios and comparing with data from non-human primate studies, we propose a hypothesis that the durability of a concomitant immune response against Mtb is intrinsically tied to levels of tissue resident memory T cells (Trms) during primary infection, with a secondary but important role for circulating Mtb-specific T cells. Further, we compare HostSim reinfection experiments to observational TB studies from the pre-antibiotic era to predict that the upper bound of the lifespan of resident memory T cells in human lung tissue is likely 2–3 years. To the authors’ knowledge, this is the first estimate of resident memory T-cell lifespan in humans. Our findings are a first step towards demonstrating the important role of Trms in preventing disease and suggest that the induction of lung Trms is likely critical for vaccine success.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

What Controls the Entrainment Rate of Dry Buoyant Thermals with Varying Initial Aspect Ratio?

Here, this study uses theory and numerical simulations to analyze the nondimensional spreading rate a (change in radius with height) of buoyant thermals as they rise and entrain surrounding environmental fluid. A focus is on how a varies with initial thermal aspect ratio A r , defined as height divided by width of the initial buoyancy perturbation. An analytic equation for thermal ascent rate w t that depends on a is derived from the thermal-volume-averaged momentum bud get equation. The thermal top height when w t is maximum, defining a critical height z c , is inversely proportional to a. The height z c also corresponds to the thermal top height when buoyant fluid along the thermal’s vertical axis is fully replaced by entrained nonbuoyant environmental fluid rising from below the thermal. The time scale for this process is controlled by the vertical velocity of parcels rising upward through the thermal’s core. This parcel vertical velocity is approximated from Hill’s analytic spherical vortex, yielding an analytic inverse relation between a and A r . Physically, this a–A r relation is connected to changes in circulation as A r is modified. Numerical simulations of thermals with A r varied from 0.5 to 2 give a values close to the analytic theoretical relation, with a factor of ~3 decrease in a as A r is increased from 0.5 to 2. The theory also explains why a of initially spherical thermals from past laboratory and modeling studies is about 0.15. Overall, this study provides a theoretical underpinning for understanding the entrainment behavior of thermals, relevant to buoyantly driven atmospheric flows.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Machine learning-assisted elucidation of CD81–CD44 interactions in promoting cancer stemness and extracellular vesicle integrity

Tumor-initiating cells with reprogramming plasticity or stem-progenitor cell properties (stemness) are thought to be essential for cancer development and metastatic regeneration in many cancers; however, elucidation of the underlying molecular network and pathways remains demanding. Combining machine learning and experimental investigation, here we report CD81, a tetraspanin transmembrane protein known to be enriched in extracellular vesicles (EVs), as a newly identified driver of breast cancer stemness and metastasis. Using protein structure modeling and interface prediction-guided mutagenesis, we demonstrate that membrane CD81 interacts with CD44 through their extracellular regions in promoting tumor cell cluster formation and lung metastasis of triple negative breast cancer (TNBC) in human and mouse models. In-depth global and phosphoproteomic analyses of tumor cells deficient with CD81 or CD44 unveils endocytosis-related pathway alterations, leading to further identification of a quality-keeping role of CD44 and CD81 in EV secretion as well as in EV-associated stemness-promoting function. CD81 is coexpressed along with CD44 in human circulating tumor cells (CTCs) and enriched in clustered CTCs that promote cancer stemness and metastasis, supporting the clinical significance of CD81 in association with patient outcomes. Our study highlights machine learning as a powerful tool in facilitating the molecular understanding of new molecular targets in regulating stemness and metastasis of TNBC.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Two Large-Scale Meteorological Patterns are Associated with Short-Duration Dry Spells in the Northeastern United States

Large-scale meteorological pattern (LSMP)–based analysis is used in a novel way to understand meteorological conditions before and during short-duration dry spells over the northeastern United States. These LSMPs are useful to assess models and select better-performing models for future projections. Dry-spell events are identified from histograms of consecutive dry days below a daily precipitation threshold. Events lasting 12 days or longer, which correspond to ~10% of dry-spell events, are examined. The 500-hPa stream-function anomaly fields for the first 12 days of each event are time averaged, and k -means clustering is applied to isolate the dry-spell-related LSMPs. The first cluster has a strong low pressure anomaly over the Atlantic Ocean, southeast of the region, and is more common in winter and spring. The second cluster has strong high pressure over east-central North America and is most common during autumn. Over the region, both clusters have negative specific humidity anomalies, negative integrated vapor transport from the north, and subsidence associated with a midlatitude jet stream dipole structure that reinforces upper-level convergence. Subsidence is supported by cold-air advection in the first cluster and the location on the east side of the lower-level high pressure in the second cluster. Extratropical cyclone storm tracks are generally shifted southward of the region during the dry spells. Individual events lie on a continuum between two distinct clusters. These clusters have similar local, but different remote, properties. Although dry spells occur with greater frequency during drought months, most dry spells occur during nondrought months. Significance Statement: This study examines the large-scale weather patterns and meteorological conditions associated with dry-spell events lasting at least 2 weeks while affecting the northeastern United States. A statistical approach groups events together on the basis of similar atmospheric features. We find two distinct sets of patterns that we call large-scale meteorological patterns. These patterns reduce moisture, foster localized sinking, and shift the storm track southward along the Atlantic seaboard, all of which reduce precipitation. Besides greater understanding, knowing the meteorological patterns during short-term dryness in the region provides an important tool to assess how well atmospheric models reproduce these specific patterns. More dry spells occur in non-drought months than in drought months, which means that dry spells can occur without preexisting drought conditions.

58 GEOSCIENCES↗
DOE OSTI2022

Local-thermal-gradient and large-scale-circulation impacts on turbine-height wind speed forecasting over the Columbia River Basin

Abstract. We investigate the sensitivity of turbine-height wind speed forecast to initial condition (IC) uncertainties over the Columbia River Gorge (CRG) and Columbia River Basin (CRB) for two typical weather phenomena, i.e., local-thermal-gradient-induced marine air intrusion and a cold frontal passage. Four types of turbine-height wind forecast anomalies and their associated IC uncertainties related to local thermal gradients and large-scale circulations are identified using the self-organizing map (SOM) technique. The four SOM types are categorized into two patterns, each accounting for half of the ensemble members. The first pattern corresponds to IC uncertainties that alter the wind forecast through a modulating weather system, which produces the strongest wind anomalies in the CRG and CRB. In the second pattern, the moderate uncertainties in local thermal gradient and large-scale circulation jointly contribute to wind forecast anomaly. We analyze the cross section of wind and temperature anomalies through the gorge to explore the evolution of vertical features of each SOM type. The turbine-height wind anomalies induced by large-scale IC uncertainties are more concentrated near the front. In contrast, turbine-height wind anomalies induced by the local IC thermal uncertainties are found above the surface thermal anomalies. Moreover, the wind forecast accuracy in the CRG and CRB is limited by IC uncertainties in a few specific regions, e.g., the 2 m temperature within the basin and large-scale circulation over the northeast Pacific around 140∘ W.

17 WIND ENERGY↗
DOE OSTI2024

Identification of equine mares as reservoir hosts for pathogenic species of Leptospira

Equine leptospirosis can result in abortion, stillbirth, neonatal death, placentitis, and uveitis. Horses can also act as subclinical reservoir hosts of infection, which are characterized as asymptomatic carriers that persistently excrete leptospires and transmit disease. In this study, PCR and culture were used to assess urinary shedding of pathogenic Leptospira from 37 asymptomatic mares. Three asymptomatic mares, designated as H2, H8, and H9, were PCR-positive for lipL32, a gene specific for pathogenic species of Leptospira. One asymptomatic mare, H9, was culture-positive, and the recovered isolate was classified as L. kirschneri serogroup Australis serovar Rushan. DNA capture and enrichment of Leptospira genomic DNA from PCR-positive, culture-negative samples determined that asymptomatic mare H8 was also shedding L. kirschneri serogroup Australis, whereas asymptomatic mare H2 was shedding L. interrogans serogroup Icterohaemorrhagiae. Sera from all asymptomatic mares were tested by the microscopic agglutination test (MAT) and 35 of 37 (94.6%) were seropositive with titers ranging from 1:100 to 1:3200. In contrast to asymptomatic mares, mare H44 presented with acute spontaneous abortion and a serum MAT titer of 1:102,400 to L. interrogans serogroup Pomona serovar Pomona. Comparison of L. kirschneri serogroup Australis strain H9 with that of L. interrogans serogroup Pomona strain H44 in the hamster model of leptospirosis corroborated differences in virulence of strains. Since lipopolysaccharide (LPS) is a protective antigen in bacterin vaccines, the LPS of strain H9 (associated with subclinical carriage) was compared with strain H44 (associated with spontaneous abortion). This revealed different LPS profiles and immunoreactivity with reference antisera. It is essential to know what species and serovars of Leptospira are circulating in equine populations to design efficacious vaccines and diagnostic tests. Our results demonstrate that horses in the US can act as reservoir hosts of leptospirosis and shed diverse pathogenic Leptospira species via urine. This report also details the detection of L. kirschneri serogroup Australis serovar Rushan, a species and serotype of Leptospira, not previously reported in the US.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Exploring the impact of surface topography on Rayleigh-Bénard dry convection in the Pi cloud chamber using OpenFOAM: In cylindrical and rectangular geometries

The Pi convection-cloud chamber can generate steady-state turbulence in both rectangular and cylindrical shapes via Rayleigh-Bénard convection (RBC) by maintaining warm bottom and cold top surfaces. Although most experiments in the Pi chamber were conducted in cylindrical shapes, all previous Pi chamber simulations were conducted in a rectangular shape due to the limitations of those models to discretize a cylindrical domain when using the finite difference method therein. Here, we use OpenFOAM, an open-source finite-volume-based Computational Fluid Dynamics (CFD) software package, to conduct Large-Eddy Simulation (LES) of dry RBC in the Pi chamber at high Rayleigh numbers (10 8 to 10 9 ). Results show that large-scale circulation (LSC) direction varies in the chamber with a constant side wall temperature. Imposing a slight temperature imbalance at the side wall ranging from 0.1 to 0.7 degrees can lock the LSC, aligning better with Pi chamber observations, particularly at higher Rayleigh numbers. In addition, we examine the impact of surface topography on LSC and heat transfer in RBC systems within cylindrical and rectangular shapes under varying conditions. Results show that roughing top/bottom surfaces by adding bars of a few tens millimeters height can strengthen thermal plumes and enhance temperature fluctuations in the chamber. Furthermore, we observe that different bar height configurations lead to notable changes in LSC orientation and thermal stratification, highlighting the complex interactions between surface features and convection patterns. This finding highlights how surface topography and chamber geometry affect Rayleigh-Bénard convection, improving understanding of turbulent heat transfer and atmospheric boundary-layer processes. Direct Numerical Simulations (DNS) are also conducted to validate LES results. In conclusion, while LES effectively captures qualitative behaviors seen in DNS, it tends to underestimate velocity variances near walls, illustrating a trade-off between computational efficiency and accuracy.

54 ENVIRONMENTAL SCIENCES↗