Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “New York State”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Three-dimensional high-content imaging of unstained soft tissue with subcellular resolution using a laboratory-based X-ray microscope

With increasing interest in studying biological systems across spatial scales—from centimeters down to nanometers—histology continues to be the gold standard for tissue imaging at cellular resolution, providing an essential bridge between macroscopic and nanoscopic analysis. However, its inherently destructive and two-dimensional nature limits its ability to capture the full three-dimensional complexity of tissue architecture. Here, we show that phase-contrast X-ray microscopy can enable three-dimensional virtual histology with subcellular resolution. This technique provides direct quantification of electron density without restrictive assumptions, allowing for direct characterization of cellular nuclei in a standard laboratory setting. By combining high spatial resolution and soft tissue contrast, with automated segmentation of cell nuclei, we demonstrated virtual Hematoxylin and Eosin (H&E) staining using machine learning-based style transfer, yielding volumetric datasets compatible with existing histopathological analysis tools. Furthermore, by integrating electron density and the sensitivity to nanometric features of the dark field contrast channel, we achieve stain-free, high-content imaging capable of distinguishing nuclei and extracellular matrix.

3D virtual histology↗
DOE OSTI2025

Probing Effects of Electron Addition to Polycyclic Conjugated Hydrocarbons Containing Four‐Membered Rings

The rectangular cyclobutadiene (CBD, C 4 H 4 ) is a unique moiety for building nonbenzenoid polycyclic conjugated hydrocarbons with interesting electron‐accepting properties. Herein, the investigation on chemical reduction of several CBD‐containing polycyclic hydrocarbons with increasing conjugation length is reported: biphenylene (C 12 H 8 ), dimethyl[2]naphthalene (C 22 H 16 ), and tetramethyl‐dibenzo‐[3]phenylene (C 30 H 22 ). The two‐step sequential reduction is first demonstrated by in situ spectroscopic investigation and then confirmed by the isolation of single crystals of the reduced products. The X‐ray crystallographic analysis reveals the formation of several mono‐ and doubly reduced products in solvent‐separated and complexed forms. The crystal structures for both neutral parents and corresponding reduced products unravel the changes in bond alternation in each ring of the fused systems. Density functional theory (DFT) and nucleus‐independent chemical shift (NICS) scan calculations reveal that the two‐electron addition reduces the aromatic character in the benzenoid rings but has minor influence on the antiaromatic CBD rings.

Zhu, Yikun [Department of Chemistry University at ↗
DOE OSTI2026

Mortality among workers at the Rocky Flats Plant, 1951–2017

The Rocky Flats (RFs) Plant operated from 1951–1989 as part of the U.S. Department of Energy (DOE) nuclear complex. Its primary mission was weapons component fabrication, whereby workers were potentially exposed to radioactive and non-radioactive hazards. RF worker mortality was compared to the general population, and dose-response relationships between mortality and radiation organ doses were examined. RF workers first employed between 1951 and 1979 for ⩾30 d were identified (n = 9397). Vital status was determined using national and state death records up to 2017. Organ doses from external photons and neutrons irritation and internalised plutonium (Pu), americium (Am), and uranium (U) were modelled as cumulative lagged total doses per year. Beryllium exposure was evaluated as an effect modifier using data from the DOE Nationwide Beryllium Medical Program. Statistical analyses included standardised mortality ratios (SMRs), Cox proportional hazard models, and excess relative risk (ERR) models. Approximately 53.2% of workers were deceased by the end of the study. Nearly 90% were monitored for radiation exposure, with a mean weighted absorbed dose of 59.0 mGy for the lungs. Nearly 45% of workers had intakes of alpha-particle emitting radionuclides, and 46.7% were monitored for neutrons. Leading causes of death included ischemic heart disease (n = 999) and lung cancer (n = 361). The highest SMRs were observed for berylliosis (SMR: 176.9; 95% CI: 76.2, 348.7; n < 10) and asbestosis (SMR: 4.65; 95% CI: 2.23, 8.55; n = 10). Dose-response analyses showed no statistical increase in risk from low-dose radiation including lung cancer (ERR per 100 mGy: −0.02; 95% CI: −0.11, 0.08; n = 361) and Parkinson’s disease (ERR per 100 mGy: 0.13; 95% CI: −0.26, 0.31; n = 57). Approximately 45% of workers were monitored for beryllium, with a weak non-significant indication of effect modification for lung cancer risk. The RF cohort showed no evidence of a statistically significant increase in mortality from occupational radiation exposure. However, this study was limited by low statistical power, which inhibits the ability to detect effects. Future pooling of Million Person Study (MPS) cohorts will provide further insights, particularly regarding Pu as a carcinogen.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2025

Tailoring Mechanical Properties of a Multi-Principle-Element Alloy by a Multi-Length-Scale Approach

A multiscale strengthening approach was employed to tailor the mechanical properties of a NiCoCr-based multi-principal element alloy (MPEA). Initially, the addition of Mo was utilized to strengthen the alloys by achieving grain size refinement, severe lattice distortion, and a reduction in stacking fault energy (SFE). Subsequently, the Suzuki segregation phenomenon was leveraged by carefully controlling stacking fault (SF) density via pre-straining and a subsequent 500°C heat treatment. Consequently, the Mo-added MPE alloys that were pre-strained and subsequently heat-treated exhibited a true yield stress of 2.3 GPa, with a corresponding fracture elongation of 12%. The Mo promotes complex SF interactions and the SFs formed during prior straining serve as sites for Suzuki segregation during subsequent heat treatment, varying the local SFE with Cr segregation, creating a roughened landscape and hindering dislocation dynamics. This multiscale strengthening approach presents a promising route for achieving desired mechanical properties in MPEAs.

mechanical properties↗
DOE OSTI2026

Co 2 P-Pt Heterostructure Interfaces for Electrocatalytic Hydrogen Evolution

Pt-based electrocatalysts are effective for the hydrogen evolution reaction (HER); however, their limited ability to facilitate water dissociation and suboptimal hydrogen binding energy (H BE ) in alkaline electrolytes result in slow reaction kinetics, which hinders their cost-efficiency and practical applications. This study reports the synthesis of Co 2 P-Pt heterostructure nanorods using a seed-mediated growth method, producing a high density of Co 2 P-Pt interfacial sites. Density functional theory (DFT) calculations indicate that electronic interactions at these interfaces optimize H BE on Pt, while the interfacial sites promote water dissociation. The Co 2 P-Pt nanorods demonstrate an overpotential of 14 mV at 10 mA cm −2 for the HER, highlighting the potential of precisely engineered metal-metal phosphide interfaces for enhancing electrocatalytic efficiency.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Twin Polaritons: Classical versus Quantum Features in Polaritonic Spectra

Understanding whether a polaritonic phenomenon is fundamentally quantum or classical is essential for building accurate theoretical models and guiding experimental design. Here, in this work, we address this question in the context of polaritonic spectra and report an intriguing new feature: the twin polariton, an additional splitting beyond the primary resonant polariton splitting originating from vacuum field fluctuations. We show that the twin polariton persists in the many-molecule limit under permutationally symmetrical initial-state constraint and that it follows the same linear dependence on coupling strength as the primary polariton splitting. This establishes a novel mechanism by which a quantum feature (the twin polariton) can be tuned through a classical one (the primary polariton), offering new opportunities to probe and control the fundamental nature of polaritonic systems.

classical optics↗
DOE OSTI2025

Tuning Solvation Dynamics of Electrolytes at Their Eutectic Point Through Halide Identity

Deep eutectic solvents (DESs) are regarded as highly promising solvent systems for redox flow batteries. DESs, composed of choline halides (ChX, X = F − , Cl − , Br − , I − ) and ethylene glycol (EG), exhibit distinct physicochemical properties at their eutectic points, including halide-dependent phase behavior, viscosity, polarity, conductivity, and solvation dynamics. In this study, we investigate the effects of the halide identity on the solvation properties of ChX:EG mixtures at varying mol % of ChX salt content. The solvatochromic polarity based on ET(30) measurements indicates higher polarity for larger halides (I − > Br − ) than for smaller halides (Cl − > F − ), which exhibit larger compensating solvation shells. The ionic conductivity follows the trend of the solvent fluidity (the inverse of the viscosity), namely ChCl > ChBr > ChI > ChF, influenced by the ion mobility and solvodynamic radii. Measurements of the liquidus temperatures (T L ) reveal that the system with ChCl exhibits the deepest eutectic point (at ~20 mol % ChCl), while ChBr and ChI have shallower minima at ~10 mol % ChBr and ~3 mol % ChI, respectively. ChF does not display a eutectic transition but instead appears to readily supercool at salt concentrations above 30 mol % ChF. Consistent with the phase transition measurements, femtosecond transient absorption spectroscopy shows that in the ChCl system, the solvation dynamics become faster with an increasing salt concentration up to ~16.67 mol %, after which the dynamics slow down with further increases in the salt content. The ChF-based system exhibits similar behavior, though with slower dynamics. In contrast, the solvation dynamics of the systems containing ChBr and ChI monotonously slow down with an increasing salt concentration, in agreement with the phase transition measurements, which show that the eutectic points occur at low salt concentrations. These measurements suggest that the solvent composition and, in particular, the identity of the halide anion play a significant role in the solvation behavior of these ethylene-glycol-based DESs, offering a foundation for tuning the DES properties for specific applications.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Connected Traffic Signal Coordination Optimization Framework through Network-Wide Adaptive Linear Quadratic Regulator–Based Control Strategy

Traffic congestion in metropolitan areas causes several significant challenges, such as longer travel times, decreased productivity, increased fuel consumption and vehicle emissions, and even severe injuries during crashes. Traffic signal control is a management approach to reduce traffic congestion and allocate the appropriate right of way for safety and mobility efficiency, both in temporal and spatial domains. Here, this study proposes a network-wide adaptive signal control coordination optimization framework based on the linear quadratic regulator algorithm. The traffic flow conditions driven by signal control inputs are formulated based on their network-wide state-space representation. After modeling traffic control regulation constraints, an adaptive linear quadratic regulator algorithm is designed to maximize the network-wide total throughput under the current conditions. Optimal signal control split time durations for multiple intersections in the network are derived by solving the algebraic Riccati equation. Furthermore, the recursive least square parameter estimation method is employed to quantify dynamic traffic condition changes. To verify the effectiveness of this proposed signal control framework, both simulation and real-world experimental tests are conducted for multiple intersections in downtown Chattanooga, Tennessee, United States. In preparation for real-world experimental tests, pipelines for real-time data processing implementation and historical traffic flow data analysis are conducted. The test results demonstrate that the proposed control framework achieves a decrease in travel time by up to 19.4%, total time spent (TTS) by up to 11.9%, and relative queue balance (RQB) by up to 15.6%. The research findings indicate that the proposed signal control framework can be generalized to handle large scale signal control optimization network-wide.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Corrigendum to “The Emissions Model Intercomparison Project (Emissions-MIP): quantifying model sensitivity to emission characteristics” published in Atmos. Chem. Phys., 23, 14779–14799, 2023

Anthropogenic emissions of aerosols and precursor compounds are known to significantly affect the energy balance of the Earth–atmosphere system, alter the formation of clouds and precipitation, and have a substantial impact on human health and the environment. Global models are an essential tool for examining the impacts of these emissions. In this study, we examine the sensitivity of model results to the assumed height of SO 2 injection, seasonality of SO 2 and black carbon (BC) particulate emissions, and the assumed fraction of SO 2 emissions that is injected into the atmosphere as particulate phase sulfate (SO 4 ) in 11 climate and chemistry models, including both chemical transport models and the atmospheric component of Earth system models. We find large variation in atmospheric lifetime across models for SO 2 , SO 4 , and BC, with a particularly large relative variation for SO 2 , which indicates that fundamental aspects of atmospheric sulfur chemistry remain uncertain. Of the perturbations examined in this study, the assumed height of SO 2 injection had the largest overall impacts, particularly on global mean net radiative flux (maximum difference of -0.35?W?m-2), SO2 lifetime over Northern Hemisphere land (maximum difference of 0.8?d), surface SO 2 concentration (up to 59?% decrease), and surface sulfate concentration (up to 23?% increase). Emitting SO 2 at height consistently increased SO 2 and SO 4 column burdens and shortwave cooling, with varying magnitudes, but had inconsistent effects across models on the sign of the change in implied cloud forcing. The assumed SO 4 emission fraction also had a significant impact on net radiative flux and surface sulfate concentration. Because these properties are not standardized across models this is a source of inter-model diversity typically neglected in model intercomparisons. These results imply a need to ensure that anthropogenic emission injection height and SO 4 emission fraction are accurately and consistently represented in global models.

Ahsan, Hamza [Pacific Northwest National Laborator↗
DOE OSTI2025

Ultra-Fast Non-Volatile Resistive Switching Devices with Over 512 Distinct and Stable Levels for Memory and Neuromorphic Computing

Low-current multilevel programmability with inherent non-volatility and high stability of resistance states is required for both multi-bit memory storage and deep learning accelerators but is difficult to achieve. Here, in a resistive switching system, this work realizes >512 (>9 bits) distinct non-volatile conductance levels with stable retention for each state with current levels down to the nanoampere range, highly promising for potential integration with small processing nodes with ultra-low power consumption requirements. This is achieved by demonstrating a new thin film design concept that encompasses three key features: an ultra-thin epitaxial oxygen ionic switching layer that provides a tunable energy barrier at the bottom electrode, an overcoat amorphous layer that acts as an ion migration barrier for stable state retention, and a partial conductive filament as a localized electronic transport channel to the epitaxial switching layer. A large dynamic resistance range of up to seven orders of magnitude is achieved with reset-free transitions among intermediate states, and programmability is demonstrated with ultra-fast (20 ns) pulses. Artificial neural network (ANN) simulations, based on the experimental performance and its non-idealities, demonstrate close-to-ideal inference accuracies for various Modified National Institute of Standards and Technology (MNIST) data sets.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nucleation and Antiphase Twin Control in Bi 2 Se 3 via Step‐Terminated Al 2 O 3 Substrates

The epitaxial synthesis of high-quality 2D layered materials is an essential driver of both fundamental physics studies and technological applications. Bi 2 Se 3 , a prototypical 2D layered topological insulator, is sensitive to defects imparted during the growth, either thermodynamically or due to the film-substrate interaction. Here, in this study, it is shown that step-terminated Al 2 O 3 substrates with a high miscut angle (3°) can effectively suppress a particular hard-to-mitigate defect, the antiphase twin. Systematic investigations across a range of growth temperatures and substrate miscut angles confirm that atomic step edges act as preferential nucleation sites, stabilizing a single twin domain. First-principles calculations suggest that there is a significant energy barrier for twin boundary formation at step edges, supporting the experimental observations. Detailed structural characterization indicates that this twin-selectivity is lost through the mechanism of the 2D layers overgrowing the step edges, leading to higher twin density as the thickness increases. These findings highlight the complex energy landscape unique to 2D materials that is driven by the interplay between substrate properties, nucleation dynamics, and defect formation, and overcoming and controlling these are critical to improve material quality for quantum and electronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Constructing Highly Porous Low Iridium Anode Catalysts Via Dealloying for Proton Exchange Membrane Water Electrolyzers

Iridium (Ir) is the most active and durable anode catalyst for the oxygen evolution reaction (OER) for proton exchange membrane water electrolyzers (PEMWEs). However, their large-scale applications are hindered by high costs and scarcity of Ir. Lowering Ir loadings below 1.0 mgcm -2 causes significantly reduced PEMWE performance and durability. Therefore, developing efficient low Ir-based catalysts is critical to widely commercializing PEMWEs. Herein, an approach is presented for designing porous Ir metal aerogel (MA) catalysts via chemically dealloying IrCu alloys. In this study, the unique hierarchical pore structures and multiple channels of the Ir MA catalyst significantly increase electrochemical surface area (ECSA) and enhance OER activity compared to conventional Ir black catalysts, providing an effective solution to design low-Ir catalysts with improved Ir utilization and enhanced stability. An optimized membrane electrode assembly (MEA) with an Ir loading of 0.5 mg Ir cm -2 generated 2.0 A cm -2 at 1.79 V, higher than the Ir black at a loading of 2.0 mg Ir cm -2 (1.63 A cm -2 ). The low-Ir MEA demonstrated an acceptable decay rate of ≈40 µV h -1 during durability tests at 0.5 (>1200 h) and 2.0 A cm -2 (400 h), outperforming the commercial Ir-based MEA (175 µV h -1 at 2.0 mg Ir cm -2 ).

36 MATERIALS SCIENCE↗
DOE OSTI2026

Residual symmetries and scalar multiplet vacuum alignment in non-Abelian flavour models

We demonstrate that, upon minimizing a renormalizable, single-scalar potential invariant under a non-Abelian symmetry, special orientations in the associated vacuum alignment of the scalar multiplet correspond to the preservation of a discrete residual flavour symmetry in the broken phase of the theory. Conversely, we show that these special scalar alignments are perturbed when additional Lagrangian operators (e.g. renormalizable, multi-flavon operators and/or effective, higher-dimensional operators) are present that break said residual symmetry, leading to a vacuum reorientation and phenomenological consequences. We therefore construct a one-to-one correspondence principle between broken residual symmetries and vacuum alignment corrections, providing a mechanism to identify (and correct) a subtle but persistent form of phenomenologically relevant fine-tuning embedded in — but often ignored by — most successful non-Abelian flavour models. We first establish this correspondence in a set of toy models based on the S4 permutation symmetry, and then apply the lessons learned to the more realistic A4 Altarelli-Feruglio and ∆(27) Universal Texture Zero models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Convergence analysis for a nonlocal gradient descent method via directional Gaussian smoothing

We analyze the convergence of a nonlocal gradient descent method for minimizing a class of high-dimensional non-convex functions, where a directional Gaussian smoothing (DGS) is proposed to define the nonlocal gradient (also referred to as the DGS gradient). The method was first proposed in [Zhang et al., Enabling long-range exploration in minimization of multimodal functions, UAI 2021], in which multiple numerical experiments showed that replacing the traditional local gradient with the DGS gradient can help the optimizers escape local minima more easily and significantly improve their performance. However, a rigorous theory for the efficiency of the method on nonconvex landscape is lacking. In this work, we investigate the scenario where the objective function is composed of a convex function, perturbed by deterministic oscillating noise. We provide a convergence theory under which the iterates exponentially converge to a tightened neighborhood of the solution, whose size is characterized by the noise wavelength. Here, we also establish a correlation between the optimal values of the Gaussian smoothing radius and the noise wavelength, thus justifying the advantage of using moderate or large smoothing radii with the method. Furthermore, if the noise level decays to zero when approaching the global minimum, we prove that DGS-based optimization converges to the exact global minimum with linear rates, similarly to standard gradient-based methods in optimizing convex functions. Several numerical experiments are provided to confirm our theory and illustrate the superiority of the approach over those based on the local gradient.

Tran, Hoang [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2026

Fluorinated ionic liquids as gas chromatographic stationary phases for the separation of volatile per- and polyfluoroalkyl substances

Background Here, the production of fluorinated organic compounds in the manufacturing, semiconductor, and pharmaceutical industries has increased exponentially over the past decade. This rapid growth has created an urgent need for efficient chromatographic platforms capable of selectively separating these compounds from complex mixtures, not only to support industrial quality control and waste management practices, but also to enable reliable environmental monitoring of volatile fluorinated contaminants. Conventional GC stationary phases lack the fluorophilic interactions needed for highly fluorinated analytes. Consequently, there is a clear demand for specialized stationary phases designed to improve chromatographic retention and selectivity for these compounds. Results Three stationary phases composed of fluorinated ionic liquids (ILs) with varied extent of fluorination were prepared to study fluorophilic interactions with fluorinated/non-fluorinated probe molecules by gas chromatography (GC). IL stationary phases featuring linear and branched perfluoroalkyl moieties, as well as a branched alkyl moiety, were systematically investigated. Chromatographic performance was examined using fluorinated compounds and their hydrocarbon analogs, including CF 3 -substituted aromatics, aliphatic alcohols, fluorotelomer alcohols (FTOHs), and perfluoroalkenes. Measurements on 5 m and 20 m columns revealed that the IL possessing branched alkyl provided stronger dispersive and hydrogen bonding interactions toward non-fluorinated aromatic and long-chain alcohols, whereas the fluorinated ILs enhanced retention of highly fluorinated FTOHs and perfluorodecene. Comprehensive two-dimensional GC (GC × GC), using a nonpolar primary column coupled with secondary columns featuring cross-bonded poly(trifluoropropylmethyl siloxane) (Rtx-200 ms), the branched fluorinated IL, or the branched non-fluorinated IL, highlighted complementary selectivity with the branched fluorinated IL providing the strongest interactions with fluorinated analytes. Significance These results demonstrate that fluorinated IL stationary phases are promising alternatives to conventional polysiloxane stationary phases for improving the separation of per- and polyfluoroalkyl substances and related fluorinated compounds. By correlating IL structure with fluorophilic interactions, this work establishes design principles for GC stationary phases that enable enhanced selectivity for highly fluorinated analytes while maintaining complementary interactions with non-fluorinated compounds.

Comprehensive two-dimensional GC↗
DOE OSTI2025

GrainPaint: A multi-scale diffusion-based generative model for microstructure reconstruction of large-scale objects

Simulation-based approaches to microstructure generation can suffer from a variety of limitations, such as high memory usage, long computational times, and difficulties in generating complex geometries. Generative machine learning models present a way around these issues, but they have previously been limited by the fixed size of their generation area. Here, we present a new microstructure generation methodology leveraging advances in inpainting using denoising diffusion models to overcome this generation area limitation. We show that microstructures generated with the presented methodology are statistically similar to grain structures generated with a kinetic Monte Carlo simulator, SPPARKS.

36 MATERIALS SCIENCE↗
DOE OSTI2025

In-situ monitoring of atomically dispersed Pt sites supported on OMS-2 during CO 2 activation

Atomically dispersed catalysts have drawn great interest lately, as they showcase a high density of active sites, selectivity, and high turnover frequencies in oxidation chemistry due to labile oxygen activation. In contrast, the applications of these catalysts have lagged in reduction reactions due to the ambiguity caused by the sintering and restructuring of active sites. To bridge this gap, the evolution of Pt 4+ isomorphically substituted into an octahedral molecular sieve structure (OMS-2) under reductive conditions was correlatively characterized using multiple in-situ analytical techniques such as ambient pressure X-ray photoelectron spectroscopy, environmental transmission electron microscopy, and solid-state nuclear magnetic resonance. The surface dynamics of the Pt single atoms were revealed during the Reverse Water Gas Shift (RWGS) reaction, where the active sites were identified as two-coordinated platinum single atoms. Under reaction, we show nonbinding atoms adjacent to the single atoms restructured the motif of the single atoms to Pt 2+ via ion mobility of potassium, increasing the activation energy by 25.6 kJ/mol. Here, this work also highlights the potential for increased stability of the single atom sites via isomorphic substitution of the metal oxide support, since the Pt-OMS-2 catalyst retained activity for about 33 h before deactivation, after which nanoparticles were observed in TEM images. This work offers a new perspective in single atom synthesis using the metal oxide as the host for the single atom site, instead of adatoms on the surface

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗