Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface roughness”

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

Sapphire Substrate Trenching to Advance Superconducting High-Coherence Quantum Devices

Achieving longer coherence times in superconducting qubits is essential for advancing quantum information processing and enabling scalable quantum computing. This work presents an innovative fabrication strategy focused on sapphire trenching of the substrate of superconducting high-coherence quantum devices, developed by the SQMS Nanofabrication Taskforce in collaboration with the materials characterization team. In this talk, we present a method to efficiently and consistently trench into the sapphire substrate, without compromising the surface roughness or profile of the etched sapphire. In the literature, trenching in the substrate has been shown to decrease losses, unwanted coupling and thermal stress. Sapphire’s strong covalent bonds make it harder to etch compared to silicon, yet it is more favorable for high-coherence quantum devices due to its lower dielectric losses. This innovative technique opens new possibilities for superconducting qubit fabrication: we fabricated high coherence qubit chips on sapphire substrates with trenching, to demonstrate that it will help pushing towards higher coherence times. By integrating this work with ongoing junction process optimization, design innovation, materials characterization and exploration, we outline a clear pathway toward transmon coherence times reaching millisecond scales and beyond.

Garattoni, S. [Fermilab]↗
DOE OSTI2026

Sapphire Substrate Trenching to Advance Superconducting High-Coherence Quantum Devices

Achieving longer coherence times in superconducting qubits is essential for advancing quantum information processing and enabling scalable quantum computing. This work presents an innovative fabrication strategy focused on sapphire trenching of the substrate of superconducting high-coherence quantum devices, developed by the SQMS Nanofabrication Taskforce in collaboration with the materials characterization team. In this talk, we present a method to efficiently and consistently trench into the sapphire substrate, without compromising the surface roughness or profile of the etched sapphire. In the literature, trenching in the substrate has been shown to decrease losses, unwanted coupling and thermal stress. Sapphire’s strong covalent bonds make it harder to etch compared to silicon, yet it is more favorable for high-coherence quantum devices due to its lower dielectric losses. This innovative technique opens new possibilities for superconducting qubit fabrication: we fabricated high coherence qubit chips on sapphire substrates with trenching, to demonstrate that it will help pushing towards higher coherence times. By integrating this work with ongoing junction process optimization, design innovation, materials characterization and exploration, we outline a clear pathway toward transmon coherence times reaching millisecond scales and beyond.

Garattoni, S. [Fermilab]↗
DOE OSTI2024

“Green” Fabrication of High-performance Transparent Conducting Electrodes by Blade Coating and Photonic Curing on PET for Perovskite Solar Cells

This study presents an innovative material processing approach to fabricate transparent conducting electrodes (TCEs) on polyethylene terephthalate (PET) substrates using blade coating and photonic curing. The hybrid TCEs consist of a multiscale Ag network, combining silver metal bus lines and nanowires, overcoated by an indium zinc oxide layer, and then photonically cured. Blade coating ensures film uniformity and thickness control over large areas. Photonic curing, a non-thermal processing method with significantly lower carbon emissions, enhances the conductivity and transparency of the coated layers. Our hybrid TCEs achieve an average transmittance of (81 ± 0.4)% referenced to air ((90 ± 0.4)% referenced to the PET substrate) in the visible range, an average sheet resistance of (11 ± 0.5) Ω sq−1, and an average surface roughness of (4.3 ± 0.4) nm. We benchmark these values against commercial PET/TCE substrates. Mechanical durability tests demonstrate <3% change in resistance after 2000 bending cycles at a 1 in radius. The scalable potential of the hybrid TCE fabrication method is demonstrated by high uniformity and excellent properties in 7 in × 8 in large-area samples and by performing the photonic curing process at 11 m min−1. Furthermore, halide perovskite solar cells fabricated on these hybrid TCEs achieve average and champion power conversion efficiencies of (10.5 ± 1.0) % and 12.2%, respectively, and significantly outperform devices made on commercial PET/TCEs. This work showcases our approach as a viable pathway for high-speed “green” manufacturing of high-performance TCEs on PET substrates for flexible optoelectronic devices.

Bonner, Justin C↗
DOE OSTI2024

When and Where Lithium Plating Occurs, Its Correlation with Microstructure Heterogeneity, and the Mechanisms That Initiate and Self-Regulate Electrochemical Heterogeneity (A02-0444)

A microstructure scale electrochemical LIB model was used to investigate lithium plating onset, material non-uniform utilization, and in-plane heterogeneities for an NMC-graphite full cell. Model predicts active material particle surface roughness and size distribution (respectively, non-uniform curvature within and between particles) initiate in-plane heterogeneity, and that particle size heterogeneity at the separator interface controls the lithium plating preferential deposition ("Where"). These in-plane heterogeneities are then exacerbated by through-plane heterogeneities induced at fast charge as electrolyte depletion occurs and concentrates intercalation reaction near the anode-separator interface. Also, magnitude and occurrence of lithium plating is controlled by effective, or macroscale, microstructure parameters ("When"). As local states of charge start to diverge between nearby active material regions, overpotential differences induced by OCP difference kick in and contribute to reduce these SOC local heterogeneities. However, for staged materials such as graphite, with OCP profile alternating between plateaus and varying regions, this balancing mechanism is, respectively, inactive and active. This leads to a dynamic, non-monotonic, in-plane heterogeneity time evolution for state of charge and Faraday current density, for which their respective in-plane heterogeneity magnitude alternates. Such behavior has been modeled both for the whole electrode at the microstructure scale and at the particle scale. In-plane heterogeneities are usually considered to be detrimental, as they result in material non-uniform utilization (i.e., under and over stressed regions) and earlier degradations. However, this work provides a more granular approach as it discriminates between a harmful in-plane heterogeneity (non-uniform curvature) that triggers SOC in-plane heterogeneity, and a beneficial in-plane heterogeneity (Faraday current density) that contributes to reduce SOC in-plane heterogeneity. This work comprehensively explains the mechanisms that initiate, exacerbate, and regulate heterogeneity at the microstructure scale, while providing some design suggestions to reduce both in-plane and through-plane heterogeneities, as summarized in the graphical abstract.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Designing Track for Electrospinning Unit and Cost-Effective Laser Scanning System

Nanofibers are produced in the Targeted Systems Department (TSD) by applying a large voltage to the nanofiber fluid and the collection apparatus known as the electrospinner. The fiber is then shot out of nozzles and collected onto the electrospinner into a nanofiber mat. The problem is that since the nozzle heads are stationary, there is non-uniform deposition of nanofiber on the collection drum leading to variation in thickness of the produced nanofiber mat. To mitigate this issue, a reciprocating mechanism is designed to move the nozzle along the drum collector along with the ability to vary its stroke length. This design was realized using Siemens NX and several parts were printed using a 3D printer, but is, as of writing, untested. Second project involved providing an alternative method to scan a large sample that had undergone certain surface deformation by beam interaction and detect out of plane deformation such as micro-scale swelling and surface roughness. Systems to carry out this scanning already exist; however, they are expensive and produce many files that need to be stitched together and are cumbersome to deal with. The solution is to create a scanning system using an already purchased scanning laser and linear stage motor. The laser and stage motor were combined to produce length results that aligned with a digital microscope but differing height results.

Ruffolo, Leopoldo↗
DOE OSTI2024

Designing a track for an electrospinning unit and cost-effective laser scanning system

Nanofibers are produced in the Targeted Systems Department (TSD) by applying a large voltage to the nanofiber fluid and the collection apparatus known as the electrospinner. The fiber is then shot out of nozzles and collected onto the electrospinner into a nanofiber mat. The problem is that since the nozzle heads are stationary, there is non-uniform deposition of nanofiber on the collection drum leading to variation in thickness of the produced nanofiber mat. To mitigate this issue, a reciprocating mechanism is designed to move the nozzle along the drum collector along with the ability to vary its stroke length. This design was realized using Siemens NX and several parts were printed using a 3D printer, but is, as of writing, untested. Second project involved providing an alternative method to scan a large sample that had undergone certain surface deformation by beam interaction and detect out of plane deformation such as micro-scale swelling and surface roughness. Systems to carry out this scanning already exist; however, they are expensive and produce many files that need to be stitched together and are cumbersome to deal with. The solution is to create a scanning system using an already purchased scanning laser and linear stage motor. The laser and stage motor were combined to produce length results that aligned with a digital microscope but differing height results.

Ruffolo, Leopoldo↗
DOE OSTI2025

Computational and Experimental Investigation of Discrete “W” Turbulators for sCO2 Turbine Internal Cooling

Direct-fired supercritical CO2 (sCO2 ) cycles are of particular interest for power generation due to the potential for higher thermal efficiencies with inherent CO2 capture. Direct fired sCO2 turbines must withstand high temperatures and the density of the working fluid is much higher than conventional Brayton cycles. Based on conventional Brayton cycle applications, internal serpentine cooling channels with ribbed turbulator strips are a common approach to achieve the required cooling, while minimizing coolant utilization. Most of the previous studies of coolant channel heat transfer enhancements have focused on serpentine channels using air as the cooling media at Reynold’s numbers typical of Brayton cycles. This paper investigates the Nusselt number augmentation of coolant channels with no features, ribbed turbulators, and discrete “W” turbulators with surface roughness effects included and CO2 as the coolant. The numerical results are compared to experimental results using CO2 as the cooling media at Reynold’s numbers ranging from 100,000 to 300,000. The results show the Nusselt number augmentation plateaus with increasing Reynolds number, in contrast to the decay observed in prior publications with air as the coolant media. Additionally, the CFD results indicate the discrete “W” geometry provides a 10% increase in Nusselt number augmentation relative to the angled ribbed turbulators.

computational fluid dynamics (CFD)↗
DOE OSTI2025

Phenomena Identification and Ranking Table (PIRT) for Heat Pipes

This Phenomena Identification and Ranking Table (PIRT) report provides an evaluation of key phenomena affecting the performance and operational regimes of heat pipes, particularly in the context of heat pipe microreactors (HPMRs). Heat pipes are advanced passive thermal management devices that utilize phase change and capillary action to achieve efficient heat transfer. However, due to the complexity of the phenomena coupled in the heat pipe, including phase change, turbulent transition, and compressibility effects, among others, there is high uncertainty in identifying and ranking the important phenomena affecting the operation of heat pipes and the current knowledge for their modeling and simulation and experimental measurements and instrumentation. This PIRT exercise, conducted as a collaborative effort involving the Department of Energy (DOE) Microreactor Program (MRP), the Nuclear Regulatory Commission (NRC), and university partners systematically identifies, reviews, and prioritizes critical phenomena affecting the operation of heat pipes based on their importance and knowledge levels. The report analyzes phenomena with high importance and low knowledge, such as wick de-wetting, critical heat flux, contact angles, and pressure dynamics, discussing challenges and future research directions for improving their modeling and simulation and experimental measurements. Additionally, the report addresses phenomena with low knowledge that could impact heat pipe operation during non-normal or transient operation, including frozen startup, laminar to turbulent transition, geysering, wick priming, underfilling conditions, surface roughness of the wick, NCGs trapped in the wick, and the timescales of startup and shutdown. This comprehensive evaluation serves as a valuable resource for guiding future research and development efforts, supporting the successful integration of heat pipes into critical applications such as nuclear reactors, and contributing to the advancement of heat pipe technologies in safety-critical industries.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2025

Competition between roughness and strength for scale-dependent surfaces

Rocks famously have scale-dependent strength, yet the actual dependence is notoriously hard to measure or incorporate into any theoretical framework. Natural rough surfaces present an opportunity to solve the problem. Surfaces sliding in shear evolve as protrusions collide. These asperities can deform or break, thus creating a new surface shape. In particular, natural surfaces have roughness at all scales as well as scale-dependent strength. Based on a scaling analysis, we have previously suggested that the scale-dependent aspect ratio of steady-state surfaces should be proportional to the scale-dependent shear strain at yield. If true, scale-dependent strength could easily be inferred from natural surfaces. Thus, moving beyond the scaling argument to a rigorous treatment of scale-dependent strength for multiscale rough surfaces in shear is important. However, analytic frameworks for analyzing multiscale problems are challenging, as conventional continuum mechanics typically involves a single value for a material property across scales. Here, in this work, we build on the formalism of Persson (2001) that presents a method to compute contact area for rough surfaces with a prescribed topographic spectrum using a stochastic differential equation. The Persson formalism allows for plastic yield under normal loading of otherwise elastic materials and leaves open the possibility of scale-dependent yield stress. In this study, we pursue this route to develop a theory and numerical results for the yielding of a rough, elastoplastic surface with scale-dependent yield stress. Here, we examine surfaces for which the power spectrum of the topography 𝐶 and yield stress 𝑌 follow power laws as a function of scale 𝜆, such that 𝐶∼𝜆 −𝑚 and 𝑌∼𝜆 −𝑛 , respectively. In this formal treatment of the problem, we focus on surfaces in contact and the resulting yield and do not impose shear. Numerical solutions show that the deviation from the elastic scaling solution is bounded as expected by the prior 1D heuristic scaling argument that anticipates the Hurst exponent as 1−𝑛. We also show that the plasticity is expected to erode the contacts if 𝑚 is lower than 𝑛−3, which corresponds to a Hurst exponent lower than 1−𝑛/2. This result is rigorously sound for 2D, i.e., realistic surfaces, and quantitatively different than the prior scaling argument. The theory now permits a correspondingly quantitative approach to interpreting natural surfaces.

elasticity↗
DOE OSTI2024

Experimental and numerical investigation of fracture conductivity between non-smooth rock surfaces with and without proppant

The enhancement of fracture conductivity is vital for the efficient recovery of subsurface resources, such as geothermal energy and petroleum hydrocarbons. Proppants, granular materials injected into hydraulic fractures to maintain their conductivity, have been studied primarily in the context of smooth fractures (i.e., fractures between smooth rock surfaces). However, non-smooth fractures (i.e., fractures between rough rock surfaces) are common in geoenergy reservoirs and thus require further investigations. In this study, we conducted laboratory measurements of fracture conductivity on shale slabs with non-smooth surfaces and carried out numerical simulation using the lattice Boltzmann (LB) method, which aimed to investigate the conductivity of non-smooth fractures with and without proppants placement. When ceramic proppant with an areal concentration of 2 lb/ft 2 was placed in the fracture, the conductivity was enhanced by roughly 3-8 times compared to fractures without proppant. In fractures with proppant, gas-measured conductivity was higher than that measured with water due to proppant embedment caused by water. The experiments demonstrate the advantages of using proppant in fractures, even if the rock surface roughness can provide certain fracture conductivity via the self-propping mechanism. For fractures without proppants, high rock surface roughness is not necessarily favorable for enhancing fracture conductivity because the self-propping mechanism requires shear slip along the fracture surface. If there is no shear slip, high rock surface roughness can cause a detrimental effect on the fracture conductivity due to the interlocking effect. Utilizing advanced experimental equipment and LB modeling, this research explores the interplays between proppant placement, fracture geometry, and stress conditions to develop a comprehensive understanding of the productivity in non-smooth fractures. Further, the outcomes of this investigation indicate the importance of creating fractures with surface roughness during hydraulic fracturing and will contribute to the development of more efficient stimulation techniques for subsurface energy extraction.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Surface Variability Mapping and Roughness Analysis of the Moon Using a Coarse–Graining Decomposition

The lunar surface contains a wide variety of topographic shapes and features, each with different distributions and scales, and any analysis technique to objectively measure roughness must respect these qualities. Coarse-graining is a naturally scale-dependent filtering technique that preserves scale-dependent symmetries and produces coarse elevation maps that gradually erase the smaller features from the original topography. In this study of the lunar surface, we present two surface variability metrics obtained from coarse-graining lunar topography: fine elevation and coarse curvature. Both metrics are isotropic, deterministic, slope-independent, and coordinate-agnostic. Fine (detrended) elevation is acquired by subtracting the coarse elevation from the original topography and contains features that are smaller than the coarse-graining length-scale. Coarse curvature is the Laplacian of coarsened topography, and naturally quantifies the curvature at any scale and indicates whether a location is elevated or depressed relative to its neighborhood at that scale. We find that highlands and maria have distinct roughness characteristics at all length-scales. Our topographic spectra reveal four scale-breaks that mark characteristic shifts in surface roughness: 100, 300, 1,000, and 4,000 km. Comparing fine elevation distributions between maria and highlands, we show that maria fine elevation is biased toward smaller-magnitude elevations and that the maria–highland discrepancies are more pronounced at larger length-scales. Here, we also provide local examples of selected regions to demonstrate that these metrics can successfully distinguish geological features of different length-scales.

58 GEOSCIENCES↗
DOE OSTI2025

Uniting Surface Properties With Hydrodynamic Roughness in Shallow Overland Flow Models

Describing flow resistance from the properties of an underlying surface is a challenge in surface hydrology. Runoff models must specify a resistance formulation or “roughness scheme”—describing the functional relationship between flow resistance and flow depth/velocity—and its parameters. Uncertainty in runoff predictions derives from both the selected roughness scheme (e.g., Darcy Weisbach, Manning's, or laminar flow equations), and its parameterization with a roughness coefficient (e.g., Manning's n). Both choices are informed by model calibration to data, usually discharge, and, if available, velocity. In this study, a Saint Venant Equation-based runoff model is calibrated to discharge and velocity data from 112 rainfall simulator experiments. The results are used to identify the optimal roughness scheme among four widely-used options for each experiment, and to explore whether surface properties can be used to select the optimal roughness scheme and its coefficient. Among the tested roughness schemes, a transitional flow equation provided the best fit to the plurality of experiments. The most suitable roughness scheme for a given experiment was not related to measured surface properties. Regression models predicted the calibrated roughness coefficients with adjusted values between 0.48 and 0.54, depending on the roughness scheme used. Litter cover was the best predictor of the roughness coefficient, followed by soil cover and average canopy gap size. The results suggest that selection of an optimal roughness scheme based on surface properties alone remains difficult, but that once a scheme is selected, roughness coefficients can be estimated from surface properties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The Surface-Topography Challenge: A Multi-Laboratory Benchmark Study to Advance the Characterization of Topography

Surface performance is critically influenced by topography in virtually all real-world applications. The current standard practice is to describe topography using one of a few industry-standard parameters. The most commonly reported number is Ra, the average absolute deviation of the height from the mean line (at some, not necessarily known or specified, lateral length scale). However, other parameters, particularly those that are scale-dependent, influence surface and interfacial properties; for example the local surface slope is critical for visual appearance, friction, and wear. The present Surface-Topography Challenge was launched to raise awareness for the need of a multi-scale description, but also to assess the reliability of different metrology techniques. In the resulting international collaborative effort, 153 scientists and engineers from 64 research groups and companies across 20 countries characterized statistically equivalent samples from two different surfaces: a “rough” and a “smooth” surface. The results of the 2088 measurements constitute the most comprehensive surface description ever compiled. We find wide disagreement across measurements and techniques when the lateral scale of the measurement is ignored. Consensus is established through scale-dependent parameters while removing data that violates an established resolution criterion and deviates from the majority measurements at each length scale. Our findings suggest best practices for characterizing and specifying topography. The public release of the accumulated data and presented analyses enables global reuse for further scientific investigation and benchmarking.

42 ENGINEERING↗
DOE OSTI2025

Coalescence-Induced Spontaneous Shedding of Microdroplets on Superhydrophobic Surfaces Featuring Enclosed Micropillars with Hierarchical Roughness

This study investigated water vapor condensation on superhydrophobic surfaces (SHSs) featuring micropillars enclosed by wall lattices and having three-tier hierarchical roughness. A total of five samples were created with three (NW-J, W200-J and W400-J samples) having large micropillar depth (~6 μm) and two (NW-S and W200-S samples) having small micropillar depth (~1 μm). Two distinct condensate removal modes were observed during condensation: coalescence-induced jumping on samples with large micropillar depth and coalescence-induced shedding on samples with small micropillar depth. The results showed that the diameter of the shedding droplet on the W200-S sample having small micropillar depth could be as small as 107 μm, as compared to the theoretical critical diameter of 267 μm for gravitational shedding on the same sample. The enhanced functionality of the three-tier nanotextures on the W200-S sample could effectively suppress localized pinning of the three-phase contact line and Wenzel neck formation during the growth of condensate droplets. Consequently, during multidroplet coalescence, the released surface energy easily overcomes the solid–liquid adhesion, leading to spontaneous shedding of merged droplets. The inclusion of the wall lattice aids condensate growth by the droplet self-alignment along the walls and promoting coalescence. As a result, the W200-S sample exhibited the highest condensate collection as well. In conclusion, the proposed surface design has great potential for scaling up and implementation in heating, ventilation, and air-conditioning equipment due to the simplicity of the surface morphology and the facile spray-coating method used to achieve hierarchical roughness.

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

Atomic scale etching of diamond: insights from molecular dynamics simulations

Diamond is a promising material for multiple applications in quantum information processing and sensing as well as applications in microelectronics. However, diamond devices can be limited by surface defects that compromise charge stability and spin coherence, among others. Improved strategies in plasma etching of diamond could play an important role in minimizing or eliminating these defects. In this work, we explore plasma-assisted atomic scale etching of diamond using argon ions (Ar + ), hydrogen ions (H + ) and hydrogen atoms (H). We employ classical molecular dynamics (MD) simulations and test several interatomic potentials based on the Reactive Empirical Bond Order (REBO) form with comparisons to a variety of published experimental results. We performed MD simulations of low-energy hydrogen ($\leqslant$50 eV) and argon ( $\leqslant$200 eV) ion bombardment of diamond surfaces. Ar + bombardment can be used to locally smooth initially rough diamond surfaces via the formation of an amorphous C layer, the thickness of which increases with argon ion energy. Subsequent exposure with hydrogen ions (or fast neutrals) will selectively etch this amorphous C layer, leaving the underlying diamond layer mostly intact if the H energy is maintained below about 10 eV. The simulations suggest that combining Ar + smoothing with selective, near threshold energy H removal of amorphous C can be an effective strategy for diamond surface engineering, leading to more reliable and sensitive diamond color center devices.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Epoxy Coating of Biofiber: An Effective Modifier of Biofiber Physical and Flow Properties for Improved Tensile Behavior of Biofiber-Reinforced Biocomposite

Biocomposites combine renewable, plant-based fibers with degradable polymers and are an attractive option for sustainable, lightweight, and cost-effective materials with a low carbon footprint, especially for large-scale additive manufacturing. One of the major challenges in the widespread adoption of biocomposites is that their mechanical performance is significantly inferior to that of synthetic composites. Surface treatment is a common and effective technique to improve the mechanical properties of the biofibers used in biocomposites. This study aims to investigate the physical and flow properties of surface-treated biofibers, as well as the tensile properties of their PLA-based biocomposite, to gain insights into how surface treatment changes the fiber‘s characteristics and biocomposite‘s mechanical properties. Surface treatment was created using a two-component epoxy system by reacting poly(bisphenol A-co-epichlorohydrin) glycidyl end-capped (PBG) and dicyandiamide (DICY). The treatment was tested on two different biofibers (loblolly pine and corn stover fibers) with three different PBG/DICY molar ratios (0.25, 0.5, and 2). Results showed that surface-treated fibers improved the tensile strength and Young‘s modulus of the biocomposites. Loblolly pine biocomposites from fibers treated with a PBG/DICY ratio of 0.25 exhibited the best tensile properties. The surface treatment resulted in a more loosely dispersed fiber bulk structure, as evidenced by less fiber agglomeration into smaller particle sizes, higher fiber sphericity, and lower loose bulk density. This can enhance stress distribution and the overall mechanical performance of the biocomposites. Additionally, surface-treated fibers exhibited better dynamic flow properties.

36 MATERIALS SCIENCE↗
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↗