Search NASA⌕ Search

SEARCH · Search NASA

Results for “mesoscale simulations”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Sensitivity Analysis of Wind and Turbulence Predictions With Mesoscale‐Coupled Large Eddy Simulations Using Ensemble Machine Learning

Abstract Coupling between mesoscale models and large‐eddy simulation (LES) models is increasingly used to more realistically represent the wide range of scales of atmospheric motions affecting boundary layer winds and turbulence that need to be simulated accurately for applications such as wind energy. However, such mesoscale‐to‐microscale coupled modeling frameworks are potentially affected by a large number of uncertain closure parameters. Here, we investigate the sensitivity associated with six closure parameters related to a 1.5‐order subgrid‐scale turbulence closure for an ensemble of mesoscale‐coupled LES. The simulations are performed using the Weather Research and Forecasting model nested from horizontal resolutions of greater than a kilometer down to tens of meters. Closure parameters are varied to generate perturbed parameter ensembles for two case studies of highly sheared, convective boundary layers observed in the Columbia Basin of Oregon and Washington during the Second Wind Forecast Improvement Project. Machine learning algorithms are used to explore the sensitivity of LES predictions, considering the effects of the perturbed physical parameters alongside categorical factors such as the case study identity, measurement location, and LES resolution. For the conditions we examine, a single parameter, the eddy viscosity coefficient, is the dominant source of parametric sensitivity and its importance is comparable to the categorical factors for several of the simulation response variables we examine.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A GPU‐Based Ocean Dynamical Core for Routine Mesoscale‐Resolving Climate Simulations

Abstract We describe an ocean hydrostatic dynamical core implemented in Oceananigans optimized for Graphical Processing Unit (GPU) architectures. On 64 A100 GPUs, equivalent to 16 computational nodes in current state‐of‐the‐art supercomputers, our dynamical core can simulate a decade of near‐global ocean dynamics per wall‐clock day at an 8‐km horizontal resolution; a resolution adequate to resolve the ocean's mesoscale eddy field. Such efficiency, achieved with relatively modest hardware resources, suggests that climate simulations on GPUs can incorporate fully eddy‐resolving ocean models. This removes a major source of systematic bias in current IPCC coupled model projections, the parameterization of ocean eddies, and represents a major advance in climate modeling. We discuss the computational strategies, focusing on GPU‐specific optimization and numerical implementation details that enable such high performance.

Silvestri, Simone [Massachusetts Institute of Tech↗
DOE OSTI2020

Development of a Time-Height Profile Assimilation Technique for Large-Eddy Simulation

Mesoscale-to-microscale coupling (MMC) aims to address the limited scope of traditional large-eddy simulations by driving the microscale flow with information concerning large-scale weather patterns provided by mesoscale models. We present a new offline MMC technique for horizontally homogeneous microscale flow conditions, in which internal forcing terms are computed based on mesoscale time–height profiles of mean-flow quantities. The advantage of such an approach is that it can be used to drive a microscale simulation with either mesoscale or observational data, and that it does not rely on specific terms in the mesoscale budget equations, which are typically not part of the default output of a mesoscale solver. The performance of the proposed profile assimilation technique is assessed based on the simulation of a typical diurnal cycle over the Scaled Wind Farm Technology site in west Texas. Results indicate that simple data assimilation techniques lead to unphysically high levels of shear and turbulence caused by the algorithm’s inability to cope with inaccuracies in the mesoscale time–height profiles. Modifying the algorithm to account for vertical coherence in the mesoscale source terms gives the microscale solver a greater ability to correct the provided mesoscale time–height profiles, leading to improved predictions of shear and turbulence statistics. The resulting turbulence statistics are in good agreement with meteorological tower observations and simulation results obtained with state-of-the-art coupling techniques using mesoscale budget components.

data assimilation↗
DOE OSTI2020

An Adaptive-Mesh-Refinement Based Computational Tool for Simulating Catalysis at Mesoscale

In this work, we present a computational tool for mesoscale applications using open-source exascale- computing compatible adaptive-mesh-refinement (AMR) library, AMReX [2]. AMReX is software library that enables development of application solvers with block-structured Cartesian AMR. Our tool has capabilities to include realistic geometry representation, chemical species transport, reactions and thermodynamics that are critical for capturing mesoscale physics. A significant achievement is the ability of our solver to automatically import electron microscopy data in the form of a stereolithography (STL) or pixelated file format (mrc, tiff) without undergoing the tedious task of unstructured mesh generation. This feature allows for rapid simulation of catalyst particles with complex morphologies using an immersed-boundary formulation. The use of AMR allows for higher resolutions at catalyst surface interfaces, which in turn provides an accurate description of surface reactions and transport. Our solver uses a hybrid distributed and shared memory parallelism (OpenMP/GPU-based) with which strong scaling up to 10,000 processors for realistic catalyst particle simulations have been demonstrated.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Turbulent kinetic energy over large offshore wind farms observed and simulated by the mesoscale model WRF (3.8.1)

Abstract. Wind farms affect local weather and microclimates;hence, parameterizations of their effects have been developed for numerical weather prediction models.While most wind farm parameterizations (WFPs) include drag effects of wind farms, models differ on whether or not an additionalturbulent kinetic energy (TKE) source should be included in theseparameterizations to simulate the impact of wind farms on the boundary layer.Therefore, we use aircraft measurements above large offshore wind farms in stable conditions to evaluateWFP choices. Of the three case studies we examine, we find the simulated ambient backgroundflow to agree with observations of temperature stratification and winds. This agreement allows us toexplore the sensitivity of simulated wind farm effects with respect to modeling choicessuch as whether or not to include a TKE source, horizontal resolution, vertical resolutionand advection of TKE. For a stably stratified marine atmospheric boundary layer (MABL), a TKE source and a horizontal resolutionon the order of 5 km or finer are necessary to represent the impact of offshore wind farms on the MABL.Additionally, TKE advection results in excessively reduced TKE over the wind farms, which in turn causesan underestimation of the wind speed deficit above the wind farm. Furthermore, using fine verticalresolution increases the agreement of the simulated wind speed with satellite observations of surface wind speed.

17 WIND ENERGY↗
DOE OSTI2023

Non-Empirical and Self-Interaction Corrections for DFTB: Towards Accurate Quantum Simulations for Large Mesoscale Systems (Final Report)

This project was comprised of two complementary (but parallel) thrusts: (1) implementing massively-parallelized computing hardware (with new computational hardware that may replace GPUs) in the density functional tight binding (DFTB) approach and (2) developing new capabilities in DFTB to calculate the electronic structure and dynamics of large chemical systems. While classical molecular dynamics can handle hundreds of thousands of atoms, it cannot provide a first-principles based description of chemical systems at the quantum level. At the other extreme, conventional Kohn-Sham DFT methods can probe the true quantum mechanical nature of chemical systems; however, these methods cannot tackle the large sizes and length scales relevant to dynamics simulations of realistic systems. The DFTB formalism utilized in this project provides a viable approach for probing these large systems at a quantum mechanical level of detail. However, to utilize the DFTB approach for accurate calculations of electronic properties, it is crucial to incorporate quantum-based corrections in DFTB since exchange-correlation effects can still remain very strong in these large systems. At the same time, enhancing the computational efficiency of DFTB is also essential since optimal computational performance is required for addressing the large size scales associated with realistic chemical systems. As such, the new non-empirical corrections and computing hardware enhancements implemented in this project will enable accurate and computationally efficient approaches to directly probe electronic properties in these large, complex systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dislocation loop bias and void swelling in irradiated α-iron from mesoscale and atomistic simulations

Abstract Dislocation loops are ubiquitous in irradiated materials, and dislocation loop bias plays a critical role in void swelling. However, due to complicated interactions between dislocation loops and point defects, it is challenging to evaluate the bias factors of dislocation loops. Here, we determine the bias of sessile < 100 > loops in α-iron using a recently developed atomistic approach based on the lifetime of point defects. We establish a mechanistic understanding of the loop interaction based on the diffusion tendency of point defects near the loop core region. Mobile self-interstitial atoms tend to be absorbed from the edge of the loop, and a trapping region perpendicular to the habit plane of the loop exists. The dislocation loop bias is found to be substantially lower than those of straight dislocations in α-iron and should be included in swelling rate estimates. With the obtained sink strength and bias values, agreement is achieved with experimental results for both absolute values and temperature dependence.

Yu, Ziang↗
DOE OSTI2020

Benchmark problems for the Mesoscale Multiphysics Phase Field Simulator (MEMPHIS)

This report details the current benchmark results to verify, validate and demonstrate the capabilities of the in-house multi-physics phase-field modeling framework Mesoscale Multiphysics Phase Field Simulator (MEMPHIS) developed at the Center for Integrated Nanotechnologies (CINT). MEMPHIS is a general phase-field capability to model various nanoscience and materials science phenomena related to microstructure evolution. MEMPHIS has been benchmarked against a suite of reported classical phase-field benchmark problems to verify and validate the correctness, accuracy and precision of the models and numerical methods currently implemented into the code.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Modeling Offshore Wind Farm Performance in Coastal Low-Level Jets Using Coupled Mesoscale-Microscale Large Eddy Simulations

Accurately predicting wind farm reliability under complex offshore atmospheric conditions remains a key challenge, particularly during noncanonical meteorological events such as coastal low-level jets (LLJs). LLJs, characterized by strong nonmonotonic vertical shear and directional veer, depart significantly from the simplified inflow assumptions embedded in conventional design standards, low-fidelity engineering models, and microscale large eddy simulations of the atmospheric boundary layer. In this work, we use the virtual wind farm framework—an exascale, graphics processing unit–accelerated large eddy simulation platform coupled with high-fidelity aeroservoelastic turbine models and advanced mesoscale-microscale coupling via the ExaWind software stack—to investigate turbine responses under realistic LLJ forcing. Simulations are performed over the U.S. North Atlantic offshore domain with the use of meteorological inputs from New York State Energy Research and Development Authority buoy data, focusing on a representative LLJ case impacting the International Energy Agency 15 MW reference turbine. Our results show that LLJs can cause up to 50% power deficits in downstream turbine rows and significantly amplify low-speed shaft and tower loads through nonlinear coupling between complex inflow characteristics and turbine structural dynamics. Two primary mechanisms drive these load amplifications: (1) unique LLJ inflow features—including veer and vertical/lateral shear—and (2) the downstream evolution of the flow under stable thermal stratification, which suppresses turbulence mixing and alters wake recovery. These mechanisms produce streamwise variations in turbine loading not captured by standard hub height–based metrics or existing design load case (DLC) definitions. This study highlights the critical role of rotor-scale flow gradients in driving fatigue and system-level aeroelastic responses, challenging current DLC and control strategies. We advocate the integration of full-flow field, environment-aware wind inputs into load modeling and control algorithms. By leveraging exascale computing to resolve mesoscale-microscale coupling, this work lays the groundwork for next-generation offshore wind turbine design and operation in meteorologically complex marine environments.

17 WIND ENERGY↗
DOE OSTI2025

Toward Understanding the Differences between Mesoscale and Large-Eddy Simulations of Tropical Cyclones

In this work, we investigate the ability of mesoscale and large-eddy simulation (LES) model configurations to predict the mean wind speed profile within the boundary layer of tropical cyclones (TCs). To this end, we perform idealized simulations of five hypothetical intense storms ranging from categories 1 to 5 on the Saffir–Simpson scale and extract time-averaged quantities near the eyewall region. We compare the model-generated data against mean wind speed profiles compiled from dropsondes launched from reconnaissance aircraft operating in the North Atlantic basin. Our analysis shows that mesoscale- and LES-generated mean wind fields display important differences in the boundary layer, including the magnitude of shear as well as the height where their low-level wind speed maxima are located. In addition, a comparison between the two model configurations with the dropsonde data shows that both modeling approaches are unable to capture the typical structure of mean winds in the lower part of the TC boundary layer (10–500 m), calling into question the use of simulations of near-axisymmetric storms for investigating the wind structure of past events. To better understand these differences, we conduct a momentum-budget analysis and show that modeled turbulent fluxes are underestimated in the mesoscale boundary layer parameterization compared to the LES model. Based on the analysis of the horizontal turbulent fluxes and their potential impact on mean flow quantities, a TC-specific boundary layer parameterization may be needed.

17 WIND ENERGY↗
DOE OSTI2022

Mesoscale evaluation of AMPS using AWARE radar observations of a wind and precipitation event over the Ross Island region of Antarctica

Surface, upper-air, and radar observations are used to assess the performance of the Antarctic Mesoscale Prediction System (AMPS) in simulating the mesoscale aspects of a wind and precipitation event over the Ross Island region of Antarctica that spanned January 16–20, 2016. The observations, collected during the Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE), provide a unique dataset for evaluating AMPS, especially the radar observations that facilitate a three-dimensional depiction of winds and precipitation. Comparisons of AMPS forecast data with surface meteorology, balloon-sounding, and profiling radar observations at and above sites near McMurdo Station reveal a mixture of similarities and differences. A generally southerly flow is evident at low levels in both the AMPS simulations and observed Doppler radial velocities. AMPS winds are comparable to those observed at the surface and aloft in terms of magnitude, direction, and timing but the strongest simulated southerly flow is displaced eastward relative to the observations. AMPS-simulated reflectivity over the broader Ross Island region is more limited in areal extent and smaller in magnitude than observed by a scanning Doppler radar. Three episodes of surface precipitation are observed near McMurdo Station over the five-day event with peak rates of ~3 mm h -1 and a total accumulation of ~22 mm. However, AMPS produces no surface precipitation at that location over the five-day event due to a low-level dry bias in the forecasts. Herein, the results show the first observationally based three-dimensional understanding of meteorology in the Ross Island region.

AMPS↗
DOE OSTI2021

In-Situ Ion Irradiation to Add Irradiation Assisted Grain Growth to the Marmot Tool

In the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the MARMOT mesoscale fuel performance tool is used to inform the development of mechanistic materials models for the BISON fuel performance tool. The grain size of the fuel has a large impact on its performance, directly impacting heat conduction, fission gas release, creep, and fracture. Thus, atomistic and mesoscale MARMOT simulations have been used to investigate grain boundary migration and grain growth in UO 2 . However, the current grain growth models in MARMOT does not consider the effect of irradiation on causing grain growth. This research project experimentally studied and implemented irradiation effects on grain growth of UO 2 into MARMOT.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Evaluation of Mesoscale Convective Systems in Climate Simulations: Methodological Development and Results from MPAS-CAM over the U.S.

In this study, a process-oriented approach is developed to evaluate warm-season mesoscale convective system (MCS) precipitation and their favorable large-scale meteorological patterns (FLSMPs) over the U.S. This approach features a novel observation-driven MCS-tracking algorithm using infrared brightness temperature and precipitation feature at 12, 25 and 50 km resolution and metrics to evaluate the model large-scale environment favorable for MCS initiation. The tracking algorithm successfully reproduces the observed MCS statistics from a reference 4-km radar MCS database. To demonstrate the utility of the new methodologies in evaluating MCS in climate simulations with mesoscale resolution, the process-oriented approach is applied to two climate simulations produced by the Variable-Resolution Model for Prediction Across Scales coupled to the Community Atmosphere Model physics, with refined horizontal grid spacing at 50 km and 25 km over North America. With the tracking algorithm applied to simulations and observations at equivalent resolutions, the simulated number of MCS and associated precipitation amount, frequency and intensity are found to be consistently underestimated in the Central U.S., particularly from May to August. The simulated MCS precipitation shows little diurnal variation and lasts too long, while MCS precipitation area is too large and intensity is too weak. The model is able to simulate four types of observed FLSMP associated with frontal systems and low-level jets (LLJ) in spring, but the frequencies are underestimated because of low-level dry bias and weaker LLJ. Precipitation simulated under different FLSMPs peak during daytime, in contrast to the observed nocturnal peak. Implications of these findings for future model development and diagnostics are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets - Part II: Continuum modeling

Part I of this two-part paper describes enhancements of a mesoscale model to include: many non-persistent fluid saturated joints with finite area, normal compliance with hysteresis and the weakening effects of angular deviation in the joint sets and of fluid pressure in the joints. The objective of this Part II is to use the results of the mesoscale model as simulated data to inspire functional forms in a continuum model for rock masses with multiple compliant fluid saturated joint sets which include these enhancements. Also, an algorithm is proposed for reducing the number of joint sets being modeled. Examples are presented here which show that the continuum model can predict similar anisotropic response to that observed in mesoscale simulations. In contrast with the mesoscale model, the continuum model can be used for large scale simulations of real applications, smoothly transitioning from coupled inelastic thermomechanical response near an explosive source to orthotropic elastic response far from the source.

58 GEOSCIENCES↗
DOE OSTI2021

Mesoscale Convective Systems in a Superparameterized E3SM Simulation at High Resolution

Accurately representing mesoscale convective systems (MCSs) is crucial to simulating the energy and water cycles in global climate models. Using a novel MCS identification and tracking algorithm applied to observations and model simulations, we evaluate how well the Energy Exascale Earth System Model simulates MCSs over the central US. Simulations performed by E3SM at 25 km grid spacing with (SP-E3SM) and without (E3SM) superparameterization are compared and evaluated against observations using multiple metrics highlighting important MCS characteristics. Compared to E3SM, SP-E3SM better simulates MCS number and MCS precipitation amount, diurnal cycle, propagation, and the probability distribution of precipitation rate in both spring and summer. The improvement from SP is partly contributed by improvement in simulating the large-scale environments, featuring enhanced atmospheric low-level moisture and larger moisture transport to the central US relative to E3SM. However, SP-E3SM still underestimates MCS precipitation amount, particularly in summer. This underestimation is closely related to the negative bias in MCS precipitation intensity, although the drier environments simulated during summer also contributes to underestimation of MCS frequency. Without SP, the larger bias in MCS precipitation amount is closely related to the negative bias in MCS frequency, while the substantial dry bias in the large-scale environments also contributes to underestimation of MCS intensity. Our results suggest that SP improves MCS simulation by improving modeling of the large-scale environments and convection initiation, which are both major limiting factors in E3SM even at 25 km grid spacing where deep convection is represented by a cumulus parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The Role of Local Inhomogeneities on Dendrite Growth in LLZO-Based Solid Electrolytes

The majority of the ceramic solid electrolytes (LLZO, LATP) demonstrate polycrystalline grain/grain-boundary (G/GB) microstructure. Higher lithium (Li) concentration and lower mechanical stiffness result in current focusing at the GBs. Growth of Li dendrites through local inhomogeneities and subsequent short circuit of the cell is a major concern. Recent studies have revealed that bulk Li metal is a viscoplastic material that has low (~0.3 MPa) and high (~1.0 MPa) yield strength during deformation at smaller and larger rates of strain, respectively. It has been argued that during deposition at smaller current densities, due to its lower yield strength, Li metal should demonstrate plastic flow against stiff ceramic electrolytes, and Li dendrites will be prevented from penetrating through solid electrolytes. In this manuscript, a multiscale modeling framework has been developed for predicting properties of GBs and the bulk of ceramic electrolytes using atomistic calculations for input to mesoscale models. Using the parameters obtained from the atomistic simulations, the mesoscale model reveals that, given enough time, even at low charge rates, lithium dendrites can grow through the GBs of LLZO. The present multiscale model results also provide information regarding the dendrite growth velocity through LLZO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗