Search NASA⌕ Search

SEARCH · Search NASA

Results for “eddying model”

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

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

At least 181 records · Page 10

High-fidelity forced convection simulations of the University of Wisconsin–Madison air-cooled reactor cavity cooling system

Among the various types of advanced nuclear technologies that exist, high-temperature gas-cooled reactors (HTGRs) stand out for their co-generation capabilities and exceptional passive safety systems. HTGRs rely on a reactor cavity cooling system (RCCS) to remove decay heat from the reactor pressure vessel (RPV) during the course of long-term transients. Several studies have focused on experimental and numerical modeling of the RCCS; however, the literature reveals a notable gap when it comes to high-fidelity simulations such as direct numerical simulations (DNS) or large eddy simulation (LES) models of the RCCS. An LES model was developed for forced convection inside the RCCS so as to obtain a high-fidelity representation of the RCCS at a scale never before attempted. Furthermore, the LES model was compared against three different Reynolds-averaged Navier–Stokes (RANS) models (Realizable k − ϵ, k − ω SST, and Reynolds stress model [RSM]). This comparison verified the robustness and accuracy of the turbulence modeling approach of the RANS models for the RCCS. Finally, based on the comparison between the LES and the RANS models, the k − ω SST RANS model of the experimental facility was developed and compared against the available experimental data.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Sensitivities of Large Eddy Simulations of Aerosol Plume Transport and Cloud Response

Abstract Cloud responses to surface‐based sources of aerosol perturbation partially depend on how turbulent transport of the aerosol to cloud base affects the spatial and temporal distribution of aerosol. Here, scenarios of plume injection below a marine stratocumulus cloud are modeled using large eddy simulations coupled to a prognostic bulk aerosol and cloud microphysics scheme. Both passive plumes, consisting of an inert tracer, and active plumes are investigated, where the latter are representative of saltwater droplet plumes such as have been proposed for marine cloud brightening. Passive plume scenarios show higher in‐plume cloud brightness (relative to out‐of‐plume) due to the predominant transport of the passive plume tracer from the near‐surface to the cloud layer within updrafts. These updrafts rise into brighter areas within the cloud deck, even in the absence of an aerosol perturbation associated with an active plume. Comparing albedo at in‐plume to out‐of‐plume locations associates the inert plume with the brightest cloud locations, without any causal effect of the plume on the cloud. Numerical sensitivities are first assessed to establish a suitable model configuration. Then sensitivity to particle injection rate is investigated. Trade‐offs are identified between the number of injected particles and the suppressive effect of droplet evaporation on plume loft and spread. Furthermore, as the near‐field in‐plume brightening effect does not depend significantly on injection rate given a suitable definition of perturbed versus unperturbed regions of the flow, plume area is a key controlling factor on the overall cloud brightening effect of an aerosol perturbation.

54 ENVIRONMENTAL SCIENCES↗

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗

Large-Eddy Simulation of a wind turbine using a Filtered Actuator Line Model

When dealing with multirotor devices such as quadcopters or wind farms, the cost of blade-resolved large-eddy simulation (LES) becomes prohibitive. Combining LES with a family of lower-fidelity models, called actuator line models (ALMs), has grown in popularity in the past decade. ALM replaces full blade resolution with an array of actuator points or lines parameterized by aerodynamic lift/drag polar plots along the blades. Body forces computed based on these actuator points are then projected onto the LES flow mesh, mimicking the effect of rotating blades on the flow. However, the optimal projection radius and the associated LES grid size is often too restrictive for multirotor simulations. Recently, a new tip-correction-based filtered ALM (F-ALM) was proposed by Martinez-Tossas and Meneveau (2019), which allows coarser-than-optimal grids by avoiding the associated overprediction of thrust. In this work, F-ALM is implemented into a high-order, in-house LES code to simulate National Renewable Energy Laboratory Phase VI wind turbine. It is then followed by a comparison between the baseline ALM and the newly implemented F-ALM in terms of instantaneous and time-averaged flow fields and blade loads, revealing the advantages of F-ALM in preventing the overprediction of power on coarse grids. Finally, this encourages accurate and affordable simulations of multirotor devices in the future.

17 WIND ENERGY↗

Data-driven enhancement of coherent structure-based models for predicting instantaneous wall turbulence

Predictions of the spatial representation of instantaneous wall-bounded flows, via coherent structure-based models, are highly sensitive to the geometry of the representative structures employed by them. In this study, we propose a methodology to extract the three-dimensional (3-D) geometry of the statistically significant eddies from multi-point wall-turbulence datasets, for direct implementation into these models to improve their predictions. The methodology is employed here for reconstructing a 3-D statistical picture of the inertial wall coherent turbulence for all canonical wall-bounded flows, across a decade of friction Reynolds number (Re T ). These structures are responsible for the Re T -dependence of the skin-friction drag and also facilitate the inner-outer interactions, making them key targets of structure-based models. The empirical analysis brings out the geometric self-similarity of the large-scale wall-coherent motions and also suggests the hairpin packet as the representative flow structure for all wall-bounded flows, thereby aligning with the framework on which the attached eddy model (AEM) is based. The same framework is extended here to also model the very-large-scaled motions, with a consideration of their differences in internal versus external flows. Implementation of the empirically-obtained geometric scalings for these large structures into the AEM is shown to enhance the instantaneous flow predictions for all three velocity components. Finally, an active flow control system driven by the same geometric scalings is conceptualized, towards favourably altering the influence of the wall coherent motions on the skin-friction drag.

42 ENGINEERING↗

Characterizing Performance of Freshwater Wetland Methane Models Across Time Scales at FLUXNET-CH 4 Sites Using Wavelet Analyses

Process-based land surface models are important tools for estimating global wetland methane (CH 4 ) emissions and projecting their behavior across space and time. So far there are no performance assessments of model responses to drivers at multiple time scales. In this study, we apply wavelet analysis to identify the dominant time scales contributing to model uncertainty in the frequency domain. We evaluate seven wetland models at 23 eddy covariance tower sites. Our study first characterizes site-level patterns of freshwater wetland CH 4 fluxes (FCH 4 ) at different time scales. A Monte Carlo approach was developed to incorporate flux observation error to avoid misidentification of the time scales that dominate model error. Our results suggest that (a) significant model-observation disagreements are mainly at multi-day time scales (<15 days); (b) most of the models can capture the CH 4 variability at monthly and seasonal time scales (>32 days) for the boreal and Arctic tundra wetland sites but have significant bias in variability at seasonal time scales for temperate and tropical/subtropical sites; (c) model errors exhibit increasing power spectrum as time scale increases, indicating that biases at time scales <5 days could contribute to persistent systematic biases on longer time scales; and (d) differences in error pattern are related to model structure (e.g., proxy of CH 4 production). Our evaluation suggests the need to accurately replicate FCH 4 variability, especially at short time scales, in future wetland CH 4 model developments.

58 GEOSCIENCES↗

Non-locality of mean scalar transport in two-dimensional Rayleigh–Taylor instability using the macroscopic forcing method

The importance of non-locality of mean scalar transport in two-dimensional Rayleigh–Taylor Instability (RTI) is investigated. The macroscopic forcing method is utilized to measure spatio-temporal moments of the eddy diffusivity kernel representing passive scalar transport in the ensemble averaged fields. Presented in this work are several studies assessing the importance of the higher-order moments of the eddy diffusivity, which contain information about non-locality, in models for RTI. First, it is demonstrated through a comparison of leading-order models that a purely local eddy diffusivity is insufficient to capture the mean field evolution of the mass fraction in RTI. Therefore, higher-order moments of the eddy diffusivity operator are not negligible. Models are then constructed by utilizing the measured higher-order moments. It is demonstrated that an explicit operator based on the Kramers–Moyal expansion of the eddy diffusivity kernel is insufficient. An implicit operator construction that matches the measured moments is shown to offer improvements relative to the local model in a converging fashion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Large-eddy simulation and Reynolds-averaged Navier-Stokes modeling of three Rayleigh-Taylor mixing configurations with gravity reversal

In this work, high-fidelity large-eddy simulation (LES) is performed of Rayleigh-Taylor (RT) mixing in three different configurations involving gravity reversal. In each configuration, LES results are compared with one-dimensional Reynolds-averaged Navier-Stokes (RANS) results, and a deficiency in a commonly used transport equation for the mass-flux velocity, $a_j$, is identified. In the first configuration, a classical two-component RT mixing layer is allowed to develop before it is subjected to rapid acceleration reversal. In the second configuration, a three-component RT mixing layer with an intermediate density layer is allowed to develop before being subjected to rapid acceleration reversal. Finally, in the third configuration, a light layer is interposed between two heavy layers; in this configuration, only one interface is RT-unstable at a time as it undergoes rapid acceleration reversal. In all cases, a commonly used buoyancy production closure in the $a_j$ transport equation is shown to lead to significant over-prediction of mixing layer growth after gravity reversal. An alternative formulation for this closure is then presented which is shown to more accurately capture the stabilization effect of gravity reversal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of turbulence-radiation coupling in pressurized oxy-combustion

Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emission in power generation. Due to the strong impact of pressure on thermal radiation, heat transfer in POC differs significantly from the situation of conventional atmospheric pressure combustion. Thus thermal radiation in POC needs to be investigated to aid new combustor development. The present computational work is a step in this direction, initiating a systematic analysis of thermal radiation and heat transfer in a pilot-scale POC combustor, which has been developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. While most previous computational works on flame radiation focused on atmospheric pressure condition, this work considered a pressurized flame. Specifically, a 15-bar POC combustor of power 50 kW is modeled employing the Ansys FLUENT commercial platform, using both Reynolds-averaged Navier-Stokes (RANS) modeling and large-eddy simulation (LES). A recently published and validated global radiation model is used to predict the radiative property of the flue gas. The discrete ordinates (DO) radiation model is chosen to solve the radiative transport equation. Incident radiation on the walls of the combustor is identified and investigated. It is revealed that for this pilot-scale, pressurized combustor, thermal radiation exhibits similar patterns in the RANS and LES models. However, the magnitudes of radiation are different in both models, presumably because interactions between radiation and turbulence are embedded in the LES model but not in the RANS model. This difference thereby underlines the major goal of this research: investigating the impact of such turbulence-to-radiation coupling on thermal radiation in a pilot-scale POC setting. Most previous studies on turbulence-to-radiation interactions focused on bench-scale, atmospheric pressure flames. This work is the first effort to extend this study to a pilot-scale, pressurized flame.

Li, Lei↗

Bayesian Cloud Property Retrievals from ARM Active and Passive Measurements

The optimum use of the continuous measurements of thermodynamics, radiation, aerosols, clouds and precipitation from the DOE Atmospheric Radiation Measurement (ARM) program is key to achieve the DOE Atmospheric System Research (ASR)’s objectives. One of the key mission requirements is to retrieve cloud and precipitation properties, as well as vertical motion parameters, along the vertical cross- section defined by the profiling active sensors. Such retrievals are challenging to perform continuously in the entire spectrum of cloud and precipitation conditions due to the large natural microphysical and dynamical variability, the often-limited information content in the measurements, and the lack of proper characterization of measurement quality and uncertainty. Today, the acquisition of new remote and in-situ sensors by the ARM program creates opportunities to address the microphysical retrieval problem by exploiting new, more robust retrieval techniques and integrating various scattered advancements in both sensor techniques and retrieval algorithms. During this project, we constructed a robust Bayesian Markov chain Monte Carlo (MCMC) cloud property retrieval algorithm that includes a state of the art radar forward model. Our MCMC-based retrieval produces both the best estimate of height-resolved cloud and precipitation properties in the radar profile, as well as an estimate of the in-cloud vertical motion and turbulence. In addition, the MCMC algorithm automatically produces robust and flexible estimates of retrieval uncertainty. We tested the algorithm on several synthetic cloud profiles obtained from large eddy simulation (LES) models with bin-resolved microphysics.

54 ENVIRONMENTAL SCIENCES↗

What Controls Crystal Diversity and Microphysical Variability in Cirrus Clouds?

Abstract Variability of ice microphysical properties like crystal size and density in cirrus clouds is important for climate through its impact on radiative forcing, but challenging to represent in models. For the first time, recent laboratory experiments of particle growth (tied to crystal morphology via deposition density) are combined with a state‐of‐the‐art Lagrangian particle‐based microphysics model in large‐eddy simulations to examine sources of microphysical variability in cirrus. Simulated particle size distributions compare well against balloon‐borne observations. Overall, microphysical variability is dominated by variability in the particles' thermodynamic histories. However, diversity in crystal morphology notably increases spatial variability of mean particle size and density, especially at mid‐levels in the cloud. Little correlation between instantaneous crystal properties and supersaturation occurs even though the modeled particle morphology is directly tied to supersaturation based on laboratory measurements. Thus, the individual thermodynamic paths of each particle, not the instantaneous conditions, control the evolution of particle properties.

54 ENVIRONMENTAL SCIENCES↗

A priori analysis on deep learning of subgrid-scale parameterizations for Kraichnan turbulence

In the present study, we investigate different data-driven parameterizations for large eddy simulation of two-dimensional turbulence in the a priori settings. These models utilize resolved flow field variables on the coarser grid to estimate the subgrid-scale stresses. We use data-driven closure models based on localized learning that employs a multilayer feedforward artificial neural network with point-to-point mapping and neighboring stencil data mapping, and convolutional neural network fed by data snapshots of the whole domain. The performance of these data-driven closure models is measured through a probability density function and is compared with the dynamic Smagorinksy model (DSM). The quantitative performance is evaluated using the cross-correlation coefficient between the true and predicted stresses. We analyze different frameworks in terms of the amount of training data, selection of input and output features, their characteristics in modeling with accuracy, and training and deployment computational time. We also demonstrate computational gain that can be achieved using the intelligent eddy viscosity model that learns eddy viscosity computed by the DSM instead of subgrid-scale stresses. Here, we detail the hyperparameters optimization of these models using the grid search algorithm.

42 ENGINEERING↗

Turbulence modeling to aid tidal energy resource characterization in the Western Passage, Maine, USA

Numerical models combined with field measurements are regularly used to characterize tidal energy resources at potential energetic sites. However, most existing works only focus on the tidal hydrodynamic characteristics, and turbulence parameters are often not reported because of the lack of high-quality turbulence measurements and the limitations of numerical models in resolving turbulent eddies. In this study, we used FVCOM - a hydrostatic primitive equation (HPE) model - to characterize the tidal energy resource in the Western Passage, Maine, USA, by taking care of the essential macro-scale turbulence properties. We observed an excellent model performance using the Mellor-Yamada Level 2.5 Turbulence Model; estimating the spatial and vertical distribution of the turbulent kinetic energy and intensity added a new perspective to the site ranking for tidal energy converter (TEC) deployments. In addition, we also examined the role of channel geometry and bathymetry, such as headlands and underwater sills, in enhancing turbulent eddies around potential TEC siting locations. Ultimately, the detailed analysis of the turbulent flow characteristics has changed the site-ranking results and demonstrated that the regional-scale HPE models could be used for the relative understanding of more or less turbulent sites for a refined resource assessment.

16 TIDAL AND WAVE POWER↗

Methods to Evaluate Subcolumn Profiles Based on Two-Point Diagnostics

In atmospheric models, stochastic generation of subgrid-scale profiles or “subcolumns” has been used for a variety of purposes. Such subcolumns can be generated from subgrid probability density functions (PDFs) at different vertical levels, when such PDFs are available. To do so, the generator needs to decide how strongly points should be correlated in the vertical, that is, how much the values should be overlapped. This is sometimes called “PDF overlap.” To assess vertical correlation in a simplified, observable setting, here the vertical correlation of vertical velocity in subcloud layers is examined. Doppler lidar is used to evaluate the vertical profiles of vertical velocity produced by a large-eddy simulation (LES) model and the Subgrid Importance Latin Hypercube Sampler (SILHS) subcolumn generator. In order to diagnose unrealistic features in subcolumn profiles, various statistical diagnostics are examined here, including the bivariate PDF of vertical velocity at two separated points (i.e., altitudes), the two-point velocity correlation, the integral correlation length, the PDF of two-point velocity differences, and the skewness and kurtosis of two-point velocity differences. The profiles produced by LES match lidar well, except that they are too smooth at small scales. The profiles produced by SILHS exhibit sharp jumps from updraft to downdraft that are not observed in the lidar data. To reduce the generation of these unrealistically sharp jumps, the SILHS sampling method is revised. The diagnostics confirm that the revised sampling method reduces the overprediction of sharp jumps.

54 ENVIRONMENTAL SCIENCES↗

Mesh and model requirements for capturing deep-stall aerodynamics in low-Mach-number flows

Here, this work presents a comprehensive computational fluid dynamics investigation of the effects of grid resolution and turbulence-model choice for capturing the unsteady three-dimensional aerodynamic performance of NACA 0012 and 0021 airfoils, with specific focus on the deep-stall regime. At high angles of attack (α), wind turbine blades routinely experience vortex-induced vibrations, which can cause significant structural damages. Accurate predictions of post-stall aerodynamics can identify the frequencies at which such vibrations maybe triggered. In this context, the NACA 0012 airfoil simulations are conducted at a chord-based Reynolds number, Re c =2×10 6 , with the k-ω Shear-Stress Transport Reynolds-Averaged Navier-Stokes (RANS) and Improved Delayed Detached Eddy Simulation (IDDES) hybrid RANS-Large Eddy Simulation turbulence models. The effect of mesh resolution both in the wall-normal and spanwise directions is investigated. Only the IDDES model with a minimum spanwise resolution of 24 cells per chord length correctly predicts the aerodynamic forces. Spectral analysis shows the peak primary shedding frequency at α=30°, which signifies the end of the stall region. In the post-stall regime, both lift and drag frequencies drop asymptotically with increasing α. The Strouhal number, based on normalised chord length, remains nearly constant in this region. Based on this study, NACA 0021 airfoil runs are performed with IDDES for Re c =2.7×10 5 and 2.0×10 6 on the finest wall-normal mesh and three spanwise grids. Simulations conducted on the finer spanwise grids demonstrate grid independence and show good agreement with experiments. The effect of varying Rec on the airfoil frequency statistics is investigated. Additionally, comparison studies are presented to investigate the impact of airfoil thickness on the frequency content at Re c =2.0×10 6 . The results from the study provide guidance on the choice of mesh resolution with the IDDES model to accurately capture aerodynamic quantities for complex industrial applications.

17 WIND ENERGY↗

Emergence of a nocturnal low-level jet from a broad baroclinic zone

An analytical model is presented for the generation of a Blackadar-like nocturnal low-level jet in a broad baroclinic zone. The flow is forced from below (flat ground) by a surface buoyancy gradient and from above (free atmosphere) by a constant pressure gradient force. Diurnally-varying mixing coefficients are specified to increase abruptly at sunrise and decrease abruptly at sunset. With attention restricted to a surface buoyancy that varies linearly with a horizontal coordinate, the Boussinesq-approximated equations of motion, thermal energy, and mass conservation reduce to a system of one-dimensional equations that can be solved analytically. Sensitivity tests with southerly jets suggest that (i) stronger jets are associated with larger decreases of the eddy viscosity at sunset (as in Blackadar theory), (ii) the nighttime surface buoyancy gradient has little impact on jet strength, and (iii) for pure baroclinic forcing (no free-atmosphere geostrophic wind), the nighttime eddy diffusivity has little impact on jet strength, but the daytime eddy diffusivity is very important and has a larger impact than the daytime eddy viscosity. The model was applied to a jet that developed in fair weather conditions over the Great Plains from southern Texas to northern South Dakota on 1 May 2020. The ECMWF Reanalysis v5 (ERA5) for the afternoon prior to jet formation showed that a broad north-south-oriented baroclinic zone covered much of the region. The peak model-predicted winds were in good agreement with ERA5 winds and lidar data from the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) central facility in north-central Oklahoma.

54 ENVIRONMENTAL SCIENCES↗

RCEMIP-II: mock-Walker simulations as phase II of the radiative–convective equilibrium model intercomparison project

Abstract. The radiative–convective equilibrium (RCE) model intercomparison project (RCEMIP) leveraged the simplicity of RCE to focus attention on moist convective processes and their interactions with radiation and circulation across a wide range of model types including cloud-resolving models (CRMs), general circulation models (GCMs), single-column models, global cloud-resolving models, and large-eddy simulations. While several robust results emerged across the spectrum of models that participated in the first phase of RCEMIP (RCEMIP-I), two points that stand out are (1) the strikingly large diversity in simulated climate states and (2) the strong imprint of convective self-aggregation on the climate state. However, the lack of consensus in the structure of self-aggregation and its response to warming is a barrier to understanding. Gaining a deeper understanding of convective aggregation and tropical climate will require reducing the degrees of freedom with which convection can vary. Therefore, we propose phase II of RCEMIP (RCEMIP-II) that utilizes a prescribed sinusoidal sea surface temperature (SST) pattern to provide a constraint on the structure of convection and move one critical step up the model hierarchy. This so-called “mock-Walker” configuration generates features that resemble observed tropical circulations. The specification of the mock-Walker protocol for RCEMIP-II is described, along with example results from one CRM and one GCM. RCEMIP-II will consist of five required simulations: three simulations with the same three mean SSTs as in RCEMIP-I but with an SST gradient and two additional simulations at one of the mean SSTs with different values of the SST gradients. We also test the sensitivity to the imposed SST gradient and the domain size. Under weak SST gradients, unforced self-aggregation emerges across the entire domain, similar to what was found in RCEMIP. As the SST gradient increases, the convective region narrows and is more confined to the warmest SSTs. At warmer mean SSTs and stronger SST gradients, low-frequency variability in the convective aggregation emerges, suggesting that simulations of at least 200 d may be needed to achieve robust equilibrium statistics in this configuration. Simulations with different domain sizes generally have similar mean statistics and convective structures, depending on the value of the SST gradient. The prescribed SST boundary condition is the only difference in the set-up between RCEMIP-II and RCEMIP-I, which enables comparison between the two; however, we also welcome participation in RCEMIP-II from models that did not participate in RCEMIP-I.

Geology↗