Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CONVECTION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

The Grell-Freitas Convection Parameterization: Recent Developments and Applications Within the NASA GEOS Global Model

We implemented and began to evaluate an alternative convection parameterization for the NASA Goddard Earth Observing System (GEOS) global model. The parameterization is based on the mass flux approach with several closures, for equilibrium and non-equilibrium convection, and includes scale and aerosol awareness functionalities. Recently, the scheme has been extended to a tri-modal spectral size approach to simulate the transition from shallow, mid, and deep convection regimes. In addition, the inclusion of a new closure for non-equilibrium convection resulted in a substantial gain of realism in model simulation of the diurnal cycle of convection over the land. Here, we briefly introduce the recent developments, implementation, and preliminary results of this parameterization in the NASA GEOS modeling system.

GEO↗

Polarimetric Radar Characteristics of Simulated and Observed Intense Convective Cores for a Midlatitude Continental and Tropical Maritime Environment

This study contrasts midlatitude continental and tropical maritime deep convective cores using polarimetric radar observables and retrievals from selected deep convection episodes during the MC3E and TWPICE field campaigns. The continental convective cores produce stronger radar reflectivities throughout the profiles, while maritime convective cores produce more positive differential reflectivity Z(dr) and larger specific differential phase K(dp) above the melting level. Hydrometeor identification retrievals revealed the presence of large fractions of rimed ice particles (snow aggregates) in the continental (maritime) convective cores, consistent with the Z(dr) and K(dp) observations. The regional cloud-resolving model simulations with bulk and size-resolved bin microphysics are conducted for the selected cases, and the simulation outputs are converted into polarimetric radar signals and retrievals identical to the observational composites. Both the bulk and the bin microphysics reproduce realistic land and ocean (L-O) contrasts in reflectivity, polarimetric variables of rain drops, and hydrometeor profiles, but there are still large uncertainties in describing Zdr and Kdp of ice crystals associated with the ice particle shapes/orientation assumptions. Sensitivity experiments are conducted by swapping background aerosols between the continental and maritime environments, revealing that background aerosols play a role in shaping the distinct L-O contrasts in radar reflectivity associated with raindrop sizes, in addition to the dominant role of background thermodynamics.

Toshihisa Matsui↗

QBO Responses to Increased CO2: The Roles of Convection and SSTs

Long-term changes to the Quasi-Biennial Oscillation (QBO) remain uncertain. Recent multimodel investigations (Richter et al. 2020a, Butchart et al. 2020) suggest that QBO amplitude will weaken in the future, but trends in QBO period and stability are unclear. Additionally, the parameterizations that generate the QBO in CMIP6 models vary greatly, and about half of CMIP6 models do not represent the QBO at all (Richter et al. 2020b). Tropical convection is likely to play a role in QBO trends, yet a majority of gravity wave parameterizations in QBO-resolving models are fixed and cannot respond to changes in convection. We consider future changes of the QBO using the NASA Goddard Institute for Space Studies (GISS) Model E2.2, which has a realistic QBO produced partly through use of a non-orographic gravity wave drag scheme that explicitly links lower stratospheric momentum fluxes to convection (Rind et al., 2014, 2020; Orbe et al. 2020). The experiments include historical integrations and abrupt increases in CO which we use to probe the roles of convection and SST feedbacks. Overall, both the QBO period and amplitude decrease in response to CO2, and our results link the period response to changes in convection and the amplitude response to changes in upwelling.

Quasi-Biennial Oscillation (QBO)↗

Characteristics and Origins of Frontal Convection During IMPACTS on 25 January 2020

This presentation will discuss the characteristics of frontal convection and its estimation from the ER-2 multifrequency nadir reflectivity and Doppler velocity from a warm-occlusion frontal system that traversed across the Northeast U.S. on 25 January 2020. This system that occurred during the IMPACTS field deployment, brought widespread clouds and precipitation to the Northeast. The system brought primarily rain in most areas of more intense rainfall associated with embedded convection near the frontal band and some freezing rain and/or snowfall in the northern regions or higher terrain due to the pre-existing low-level cold air. The NASA ER-2 instrumented with multifrequency radar and the NASA P-3 with in situ microphysics measurements both collected data from this system. The focus of this presentation is on the estimation of vertical velocity in convective and non-convective regions using the ER-2 nadir Doppler measurements. Vertical velocities are inherently difficult to estimate in frontal systems at mesoscale and below other than from direct in situ measurement because of their small magnitudes and errors in the measurements. Estimates from nadir Doppler measurements such as from the ER-2 are challenging because of both Doppler velocity uncertainties and the need to estimate particle fallspeed since the Doppler measurement is the vertical velocity plus the hydrometeor fallspeed. In the presentation, the vertical velocity estimation, its associated errors, and a possible mitigation scheme are described. The characteristics of the convective and stratiform regions on 25 January 2020 are discussed. The vertical velocity estimated from the nadir Doppler measurements are compared with P-3 in situ measurements. Estimation of vertical velocity is one of the goals for upcoming spaceborne missions such as EarthCare and ACCP.

Gerald M. Heymsfield↗

Simulation of Radon-222 with the GEOS-Chem Global Model: Emissions, Seasonality, and Convective Transport

Radon-222 (Rn-222) is a short-lived radioactive gas naturally emitted from land surfaces and has long been used to assess convective transport in atmospheric models. In this study, we simulate Rn-222 using the GEOS-Chem chemical transport model to improve our understanding of Rn-222 emissions and surface concentration seasonality and characterize convective transport associated with two Goddard Earth Observing System (GEOS) meteorological products, the Modern-Era Retrospective analysis for Research and Applications (MERRA) and GEOS Forward Processing (GEOS-FP). We evaluate four global Rn-222 emission scenarios by comparing model results with observations at 51 surface sites. The default emission scenario in GEOS-Chem yields a moderate agreement with surface observations globally (68.9 % of data within a factor of 2) and a large underestimate of winter surface Rn-222 concentrations at Northern Hemisphere midlatitudes and high latitudes due to an oversimplified formulation of Rn-222 emission fluxes (1 atom cm−2 s−1 over land with a reduction by a factor of 3 under freezing conditions). We compose a new global Rn-222 emission scenario based on Zhang et al. (2011) and demonstrate its potential to improve simulated surface Rn-222 concentrations and seasonality. The regional components of this scenario include spatially and temporally varying emission fluxes derived from previous measurements of soil radium content and soil exhalation models, which are key factors in determining Rn-222 emission flux rates. However, large model underestimates of surface Rn-222 concentrations still exist in Asia, suggesting unusually high regional Rn-222 emissions. We therefore propose a conservative upscaling factor of 1.2 for Rn-222 emission fluxes in China, which was also constrained by observed deposition fluxes of 210Pb (a progeny of Rn-222). With this modification, the model shows better agreement with observations in Europe and North America (> 80 % of data within a factor of 2) and reasonable agreement in Asia (close to 70 %). Further constraints on Rn-222 emissions would require additional concentration and emission flux observations in the central United States, Canada, Africa, and Asia. We also compare and assess convective transport in model simulations driven by MERRA and GEOS-FP using observed Rn-222 vertical profiles in northern midlatitude summer and from three short-term airborne campaigns. While simulations with both GEOS products are able to capture the observed vertical gradient of Rn-222 concentrations in the lower troposphere (0–4 km), neither correctly represents the level of convective detrainment, resulting in biases in the middle and upper troposphere. Compared with GEOS-FP, MERRA leads to stronger convective transport of Rn-222, which is partially compensated for by its weaker large-scale vertical advection, resulting in similar global vertical distributions of Rn-222 concentrations between the two simulations. This has important implications for using chemical transport models to interpret the transport of other trace species when these GEOS products are used as driving meteorology.

Bo Zhang↗

Investigating Convective Gustiness of Winds and Fluxes Using Microwave Remote Sensing

Tropical atmospheric convection is known to be associated with enhanced winds over the oceans that can drive increased rates of mass and energy transfer through turbulent latent and sensible heat fluxes. However, convection and its associated rainfall usually obscures the surface wind signals observed using passive microwave observations. Air-sea flux estimates based on the use of microwave imagers typically mask in areas of rainfall and thus are likely underrepresenting the true turbulent flux transfer. However, both monostatic and bistatic radar scatterometer systems are allowing for more detailed analysis of the convective near-surface environment. While rain still impact C-band scatterometer observations, the L-band observables from the Cyclone Global Navigation Satellite System (CYGNSS) are providing improved observing capabilities in all-weather conditions. In this study, we will present ongoing research to leverage CYGNSS observations in addition to other scatterometer, in situ, and satellite-based precipitation information to provide an improved characterization of surface wind enhancement within tropical oceanic convection. The results will specifically highlight current capabilities and limitations to assess convective wind enhancement and wind gustiness and their relationship to the surface turbulent fluxes of heat and water.

J Brent Roberts↗

Moist Convection in the Giant Planet Atmospheres

The outer planets of our Solar System display a myriad of interesting cloud features, of different colors and sizes. The differences between the types of observed clouds suggest a complex interplay between the dynamics and chemistry at play on these atmospheres. Particularly, the stark difference between the banded structures of Jupiter and Saturn, vs the sporadic clouds on the ice giants highlights the varieties in dynamic, chemical and thermal processes that shape these atmospheres. Since the early explorations of these planets by spacecrafts, such as Voyager and Voyager 2, there are many outstanding questions about the long term stability of the observed features. One hypothesis is that the internal heat generated during the formation of these planets is transported to the upper atmosphere through latent heat release from convecting clouds (i.e., moist convection). In this review, we present evidences of moist convective activity on the gas giant atmospheres of our Solar System from remote sensing data, both from ground- and space-based observations. We detail the processes that drive moist convective activity, both in terms of the dynamics as well as the microphysical processes that shape the resulting clouds. Finally, we also discuss the effects of moist convection on shaping the large scale dynamics (such as jet structures on these planets).

Planetary Science↗

Diagnosing added value of convection-permitting regional models using precipitation event identification and tracking

Dynamical downscaling with high-resolution regional climate models may offer the possibility of realistically reproducing precipitation and weather events in climate simulations. As resolutions fall to order kilometers, the use of explicit rather than parametrized convection may offer even greater fidelity. However, these increased resolutions both allow and require increasingly complex diagnostics for evaluating model fidelity. Here we focus on precipitation evaluation and analyze five 2-month-long dynamically downscaled model runs over the continental United States that employ different convective and microphysics parameterizations, including one high-resolution convection-permitting simulation. All model runs use the Weather Research and Forecasting Model driven by National Center for Environmental Prediction reanalysis data. We show that employing a novel rainstorm identification and tracking algorithm that allocates essentially all rainfall to individual precipitation events allows new insights into model biases. Results include that, at least in these runs, model wet bias is driven by excessive areal extent of individual precipitating events, and that the effect is time-dependent, producing excessive diurnal cycle amplitude. This amplified cycle is driven not by new production of events but by excessive daytime enlargement of long-lived precipitation events. Here, we further show that in the domain average, precipitation biases appear best represented as additive offsets. Of all model configurations evaluated, convection-permitting simulations most consistently reduced biases in precipitation event characteristics.

54 ENVIRONMENTAL SCIENCES↗

Efficient In‐Cloud Removal of Aerosols by Deep Convection

Abstract Convective systems dominate the vertical transport of aerosols and trace gases. The most recent in situ aerosol measurements presented here show that the concentrations of primary aerosols including sea salt and black carbon drop by factors of 10 to 10,000 from the surface to the upper troposphere. In this study we show that the default convective transport scheme in the National Science Foundation/Department of Energy Community Earth System Model results in a high bias of 10–1,000 times the measured aerosol mass for black carbon and sea salt in the middle and upper troposphere. A modified transport scheme, which considers aerosol activation from entrained air above the cloud base and aerosol‐cloud interaction associated with convection, dramatically improves model agreement with in situ measurements suggesting that deep convection can efficiently remove primary aerosols. We suggest that models that fail to consider secondary activation may overestimate black carbon's radiative forcing by a factor of 2.

54 ENVIRONMENTAL SCIENCES↗

Ascending Columns, WTG, and Convective Aggregation

Analytic solutions are derived for a convecting atmosphere with mean ascent using a zero-buoyancy bulk-plume approximation for moist convection. It has been suggested that such solutions should serve as a model for the relationship between humidity, instability, and precipitation in the tropics, but it is shown here that this interpretation is incompatible with the observed weak temperature gradient (WTG). Instead, the solutions can be used to understand the atmospheric state averaged over all tropical convecting regions. Further, using the analytic solutions in this way, they predict the changes in humidity, instability, and precipitation as a function of the size of the moist patch in a convectively aggregated state.

54 ENVIRONMENTAL SCIENCES↗

A Vertically Resolved MSE Framework Highlights the Role of the Boundary Layer in Convective Self-Aggregation

Abstract Convective self-aggregation refers to a phenomenon in which random convection can self-organize into large-scale clusters over an ocean surface with uniform temperature in cloud-resolving models. Previous literature studies convective aggregation primarily by analyzing vertically integrated (VI) moist static energy (MSE) variance. That is the global MSE variance, including both the local MSE variance at a given altitude and the covariance of MSE anomalies between different altitudes. Here we present a vertically resolved (VR) MSE framework that focuses on the local MSE variance to study convective self-aggregation. Using a cloud-resolving simulation, we show that the development of self-aggregation is associated with an increase of local MSE variance, and that the diabatic and adiabatic generation of the MSE variance is mainly dominated by the boundary layer (BL; the lowest 2 km). The results agree with recent numerical simulation results and the available potential energy analyses showing that the BL plays a key role in the development of self-aggregation. Additionally, we find that the lower free troposphere (2–4 km) also generates significant MSE variance in the first 15 days. We further present a detailed comparison between the global and local MSE variance frameworks in their mathematical formulation and diagnostic results, highlighting their differences.

54 ENVIRONMENTAL SCIENCES↗

An Investigation of LES Wall Modeling for Rayleigh–Bénard Convection via Interpretable and Physics-Aware Feedforward Neural Networks with DNS

Abstract The traditional approach of using the Monin–Obukhov similarity theory (MOST) to model near-surface processes in large-eddy simulations (LESs) can lead to significant errors in natural convection. In this study, we propose an alternative approach based on feedforward neural networks (FNNs) trained on output from direct numerical simulation (DNS). To evaluate the performance, we conduct both a priori and a posteriori tests. In the a priori (offline) tests, we compare the statistics of the surface shear stress and heat flux, computed from filtered DNS input variables, to the stress and flux obtained from the filtered DNS. Additionally, we investigate the importance of various input features using the Shapley additive explanations value and the conditional average of the filter grid cells. In the a posteriori (online) tests, we implement the trained models in the System for Atmospheric Modeling (SAM) LES and compare the LES-generated surface shear stress and heat flux with those in the DNS. Our findings reveal that vertical velocity, a traditionally overlooked flow quantity, is one of the most important input features for determining the wall fluxes. Increasing the number of input features improves the a priori test results but does not always improve the model performance in the a posteriori tests because of the differences in input variables between the LES and DNS. Last, we show that physics-aware FNN models trained with logarithmic and scaled parameters can well extrapolate to more intense convection scenarios than in the training dataset, whereas those trained with primitive flow quantities cannot. Significance Statement The traditional near-surface turbulence model, based on a shear-dominated boundary layer flow, does not represent near-surface turbulence in natural convection. Using a feedforward neural network (FNN), we can construct a more accurate model that better represents the near-surface turbulence in various flows and reveals previously overlooked controlling factors and process interactions. Our study shows that the FNN-generated models outperform the traditional model and highlight the importance of the near-surface vertical velocity. Furthermore, the physics-aware FNN models exhibit the potential to extrapolate to convective flows of various intensities beyond the range of the training dataset, suggesting their broader applicability for more accurate modeling of near-surface turbulence.

54 ENVIRONMENTAL SCIENCES↗

Towards Unifying the Planetary Boundary Layer and Shallow Convection in CAM5 with the Eddy-Diffusivity/Mass-Flux Approach

The modular structure of the boundary layer and convection parameterizations in atmospheric models have long been affecting the numerical representation of subgrid-scale motions and their mutual interactions. A promising alternative, the eddy-diffusivity/mass-flux approach (EDMF), has the potential for unifying the existing formulations into a consistent scheme and improving some of the long-standing issues. This study documents a step towards developing such a unified approach by implementing a stochastic multi-plume EDMF scheme into the Community Atmosphere Model (Version 5.0). Its performance in single-column mode is evaluated against the control parameterization and large-eddy simulation (LES) for two benchmark cases: marine and continental shallow convection. Overall, the results for the two parameterizations agree well with each other and with LES in terms of mean profiles of moist conserved variables and their vertical fluxes, as well as the updraft properties. However, systematic differences between the two schemes, especially for transient continental convection, are also documented. Using EDMF helps improve some of the parameterized features of shallow convection. In particular, for the highest tested vertical resolution, the EDMF cloud base and top heights and the vertical fluxes of energy and water are remarkably close to LES.

54 ENVIRONMENTAL SCIENCES↗

Critical heat flux in natural convection cooled TRIGA reactors with hexagonal bundle

A three-rod bundle Critical Heat Flux (CHF) study at low flow, low pressure, and natural convection condition has been conducted, simulating TRIGA reactors with the hexagonally configured core. The test section is a custom-made trefoil shape tube with three identical fuel pin heater rods located symmetrically inside. The full scale fuel rod is electrically heated with a chopped-cosine axial power profile. CHF experiments were carried out with the following conditions: inlet water subcooling from 30 K to 95 K; pressure from 110 kPa to 230 kPa; mass flux up to 150 kg/m{sup 2}s. About 50 CHF data points were collected and compared with a few existing CHF correlations whose application ranges are close to the testing conditions. Some tests were performed with the forced convection to identify the potential difference between the CHF under the natural convection and forced convection. The relevance of the CHF to test parameters is investigated. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The influence of water table depth and the free atmospheric state on convective rainfall predisposition

Here, we report a mechanistic model for the soil-plant system is coupled to a conventional slab representation of the atmospheric boundary layer (ABL) to explore the role of groundwater table (WT) variations and free atmospheric (FA) states on convective rainfall predisposition (CRP) at a Loblolly pine plantation site situated in the lower coastal plain of North Carolina. Predisposition is quantified using the crossing between modeled lifting condensation level (LCL) and convectively grown ABL depth. The LCL-ABL depth crossing is necessary for air saturation but not sufficient for cloud formation and subsequent convective rainfall occurrence. However, such crossing forms the main template for which all subsequent dynamical processes regulating the formation (or suppression) of convective rainfall operate on. If the feedback between surface fluxes and FA conditions is neglected, a reduction in latent heat flux associated with reduced WT levels is shown to enhance the ABL-LCL crossing probability. When the soil-plant system is fully coupled with ABL dynamics thereby allowing feedback with ABL temperature and humidity, FA states remain the leading control on CRP. However, vegetation water stress plays a role in controlling ABL-LCL crossing when the humidity supply by the FA is within an intermediate range of values. When FA humidity supply is low, cloud formation is suppressed independent of surface latent heat flux. Similarly, when FA moisture supply is high, cloud formation can occur independent of surface latent heat flux. In an intermediate regime of FA moisture supply, the surface latent heat flux controlled by soil water availability can supplement (or suppress) the necessary water vapor leading to reduced LCL and subsequent ABL-LCL crossing. Lastly, it is shown that this intermediate state corresponds to FA values around the mode in observed humidity lapse rates γ w (between -2.5 × 10 -6 and -1.5 × 10 -6 kg kg -1 m -1 ), suggesting that vegetation water uptake may be controlling CRP at the study site.

54 ENVIRONMENTAL SCIENCES↗

The use of idealised experiments in testing a new convective parametrization: Performance of CoMorph-A

CoMorph is a new mass-flux convection parametrization under development at the Met Office designed for use within the Unified Model and its successor model, LFRic. Use of a three-dimensional idealised model enables controlled tests of the performance of the scheme across different regimes. This includes the interaction between the physical parametrizations and the resolved dynamics, allowing study of the emergent organisation of convection on the resolved scale. A selection of well-known cases is revisited here, with the purpose of documenting the extent to which CoMorph captures a range of important, but challenging, behaviour such as the diurnal cycle and sensitivity to tropospheric moisture. Simulations using CoMorph-A, a new physics package, that has been demonstrated to perform well at numerical weather prediction (NWP) and climate scales, are compared against the current global atmosphere configuration and high-resolution results. In addition to an entirely new convection scheme, the package of changes includes significant changes to the cloud, microphysics, and boundary-layer parametrizations. Recognising that CoMorph-A is the first version of a scheme that will continue to be substantially developed and to obtain good performance, compromises in tuning have had to be made. These idealised tests therefore show what works well in this configuration, and what areas will require further work. As such, it is quite a demanding testbed and could be viewed as some of the equipment required for a “convective playground”.

54 ENVIRONMENTAL SCIENCES↗

Improved natural convection heat transfer correlations for reactor cavity cooling systems of high-temperature gas-cooled reactors: From computational fluid dynamics to Pronghorn

The Reactor Cavity Cooling System (RCCS) is a common reactor safety system in High Temperature Gas Cooled Reactors (HTGR) that removes heat from the Reactor Pressure Vessel (RPV) by radiation ($\sim 80\%$) and natural convection ($\sim 20\%$). For simulation of accident scenarios of HTGRs, intermediate fidelity and system codes models must be employed for limiting the models' execution time. While accurate quantification of the radiative heat transfer is available in these models, quantification of natural convection must rely on correlations of questionable accuracy for the Nusselt number. Commonly used correlations are based in experiments performed at low Rayleigh numbers and/or using isothermal walls in simplified geometries. Here, this work improves on the accuracy of natural convection heat transfer correlations in support for HTGR designs. These correlations include both local and average Nusselt numbers as a function of the global Rayleigh number, the local Rayleigh number, and the temperature profile at the hot wall of the RCCS. In the absence of dedicated experiments and the difficulty of performing high-fidelity simulations at realistic Rayleigh numbers, the data to fit the correlations are generated with Computational Fluid Dynamics (CFD) using Reynolds Averaged Navier-Stokes (RANS) models. First, a careful selection of the RANS turbulence model is performed by comparing the results obtained with different RANS turbulence models against high fidelity simulations of natural convection at $Ra \ 1 \times 10^{11}$ in a rectangular cavity. Next, the selected model is used to perform simulations of an HTGR cavity at different high Rayleigh numbers $\in [6.1 \times 10^{11},2.9 \times 10^{13}]$ to encompass several HTGR designs, assuming an isothermal RPV wall. The results obtained are used to fit a correlation for the average and space-varying Nusselt number as a function of the global and local Rayleigh numbers via a sparsity-promoting least-squares method. The selected RANS model is then used to perform simulations of a PBMR 400 HTGR cavity with the temperature profiles at the RPV wall obtained during a PLOFC transient. We use the results obtained to fit a temperature-dependent correction to the space-varying Nusselt number with the sparsity-promoting least-squares method. The results obtained in this work, enable system-level codes, such as Pronghorn, to perform higher-fidelity simulations of the heat exchange process in the RCCS while still maintaining a low computational cost.

42 ENGINEERING↗

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↗