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 127 records · Page 7

Onset of cavitation and vapor bubble development over hydrophilic and hydrophobic surfaces

Cavitation, the formation of vapor bubbles as the liquid pressure is reduced below the saturated vapor pressure, often requires a substantial negative relative pressure in a pure liquid. Classical nucleation theory (CNT) provides an estimate for the rate of cavitation but there is often a disconnect between the predictions at the molecular scale compared to observations at the macroscale. We report on mesoscale simulations of cavitation based on many-body dissipative particle dynamics (mDPD), a coarse-grained molecular dynamics (MD), which bridges the two scales. A liquid layer is confined between smooth planar walls at a constant temperature, while the pressure is reduced slowly by expanding the wall-bounded domain. The wetting properties of the liquid are determined by the parameters of the interaction potentials. With hydrophilic walls, homogeneous nucleation is observed in the liquid bulk. As a bubble forms and grows, it creates a strong pressure pulse and oscillations that cause other bubbles that may have formed slightly later to collapse. For a nearly neutral wall with a contact angle close to 90 ° , heterogeneous nucleation occurs at the walls at a smaller negative pressure and generates weaker pressure oscillations. With hydrophobic walls or seed particles, heterogeneous nucleation readily occurs, where fluctuations and the merger of transient surface bubbles are significant.

Science & Technology - Other Topics↗

Non-Gaussian distribution of displacements for Lennard-Jones particles in equilibrium

Most mesoscale simulation methods assume Gaussian distributions of velocity-like quantities. These quantities are not true velocities, however, but rather time-averaged velocities or displacements of particles. Here, we show that there is a large range of coarse-graining scales where the assumption of a Gaussian distribution of these displacements fails, and a more complex distribution is required to adequately express these distribution functions of displacements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Discrete ion stochastic continuum overdamped solvent algorithm for modeling electrolytes

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Prediction of Probabilistic Shock Initiation Thresholds of Energetic Materials Through Evolution of Thermal-Mechanical Dissipation and Reactive Heating

The ignition threshold of an energetic material (EM) quantifies the macroscopic conditions for the onset of self-sustaining chemical reactions. The threshold is an important theoretical and practical measure of material attributes that relate to safety and reliability. Historically, the thresholds are measured experimentally. In this work, we present a new Lagrangian computational framework for establishing the probabilistic ignition thresholds of heterogeneous EM out of the evolutions of coupled mechanical-thermal-chemical processes using mesoscale simulations. Furthermore, the simulations explicitly account for microstructural heterogeneities, constituent properties, and interfacial processes and capture processes responsible for the development of material damage and the formation of hotspots in which chemical reactions initiate. The specific mechanisms tracked include viscoelasticity, viscoplasticity, fracture, post-fracture contact, frictional heating, heat conduction, reactive chemical heating, gaseous product generation, and convective heat transfer. To determine the ignition threshold, the minimum macroscopic loading required to achieve self-sustaining chemical reactions with a rate of reactive heat generation exceeding the rate of heat loss due to conduction and other dissipative mechanisms is determined. Probabilistic quantification of the processes and the thresholds are obtained via the use of statistically equivalent microstructure sample sets (SEMSS). The predictions are in agreement with available experimental data.

36 MATERIALS SCIENCE↗

Evaluation and Uncertainty of Radioxenon Transport with a Mesoscale Model after the February 2013 Underground Test in North Korea

Abstract The transport of radioxenon released from the February 2013 underground nuclear weapons test in North Korea was analyzed at two receptors—one at the Comprehensive Nuclear Test Ban Treaty site Rn58 in Russia (400 km downwind) and a second at Rn38 in Japan (1000 km downwind). Transport was modeled with two ensembles of mesoscale simulations, one generated with varying initial and lateral boundary conditions taken from the Global Forecasting System uncertainty ensemble, and a second created from different parameterizations and surface conditions. The wind variability was similar for the two ensembles and consistent with observations at 925 mb (1 mb = 1 hPa) but not at the surface. Biases in calculated surface winds and the radioxenon concentration in the ensembles were attributed mainly to poor simulation of the sea breeze at both locations and mountain lee affects at Rn38 in Japan. These wind regimes affected the timing of the surface radioxenon plume at Rn58 and its duration at Rn38. Surface wind variability induced by terrain and land–sea contrast (the sea breeze) also had a significant effect on the surface winds and plume dynamics, including blocking of flow approaching elevated terrain near Vladivostok and the west side of Japan. Increased plume uncertainty was seen at night because of surface wind variability. Measured surface chemical variability was larger than found in the first European Tracer Experiment in central Europe. The study found that horizontal model resolution contributes to uncertainty but not as much as vertical resolution, boundary layer parameterizations, and assimilation of surface meteorological data near the receptor.

54 ENVIRONMENTAL SCIENCES↗

BISON microstructure-based pulverization criterion in high burnup structure

To improve the economics of commercial nuclear power production, utilities are seeking to increase the allowable burnup limit of UO$_2$ fuel. One of the main factors that contributes to the current burnup limit of 62 GWd/MTU in commercial light water reactors (LWRs) is the risk of fine fragmentation or pulverization during a loss of coolant accident (LOCA). Pulverization primarily occurs at high burnups, especially when the high burnup structure (HBS) has formed. To allow the industry to pursue increased burnup and develop mitigation strategies, it is essential to have improved capability to predict the onset of pulverization. However, the mechanism of pulverization is not well understood, and the existing predictive capabilities implemented in the BISON fuel performance code are empirical in nature. In this report, mesoscale simulations are used to improve understanding of the formation mechanism of the HBS and how it responds during a LOCA transient, and inform development of a BISON pulverization criterion. A phase-field model was used to simulate the evolution of bubble pressure as a result of HBS formation. The simulations showed that gas atoms diffuse from grain interiors to the new grain boundaries created during HBS formation, and diffuse rapidly along these grain boundaries to reach existing bubbles. This causes an increase in bubble pressure in existing bubbles, leading to bubble growth during steady-state operation. To simulate the response of HBS bubbles to a LOCA transient, a newly developed phase-field model was used; in agreement with preliminary results from FY20, bubble size did not change significantly during the duration of the transient. A phase-field fracture model was used to study fragmentation patterns in the HBS, including using input from the phase-field model as initial conditions. Phase-field fracture simulations were used to determine a pulverization criterion for BISON. A function for the critical pressure for pulverization to occur was fit to data from the phase-field fracture simulations, and this function was implemented as a material property in BISON. For comparison, an analytical criterion for pulverization was developed and implemented within the same material property.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Energy Research and Forecasting (ERF) v1

ERF is a next-generation regional atmospheric modeling code for simulation mesoscale atmospheric flows. Its capabilities include the ability to model advection, diffusion, turbulence, terrain and moisture physics.

Almgren, Ann↗

Comparison of Gravity Wave Temperature Variances from Ray-Based Spectral Parameterization of Convective Gravity Wave Drag with AIRS Observations

The realism of ray-based spectral parameterization of convective gravity wave drag, which considers the updated moving speed of the convective source and multiple wave propagation directions, is tested against the Atmospheric Infrared Sounder (AIRS) onboard the Aqua satellite. Offline parameterization calculations are performed using the global reanalysis data for January and July 2005, and gravity wave temperature variances (GWTVs) are calculated at z = 2.5 hPa (unfiltered GWTV). AIRS-filtered GWTV, which is directly compared with AIRS, is calculated by applying the AIRS visibility function to the unfiltered GWTV. A comparison between the parameterization calculations and AIRS observations shows that the spatial distribution of the AIRS-filtered GWTV agrees well with that of the AIRS GWTV. However, the magnitude of the AIRS-filtered GWTV is smaller than that of the AIRS GWTV. When an additional cloud top gravity wave momentum flux spectrum with longer horizontal wavelength components that were obtained from the mesoscale simulations is included in the parameterization, both the magnitude and spatial distribution of the AIRS-filtered GWTVs from the parameterization are in good agreement with those of the AIRS GWTVs. The AIRS GWTV can be reproduced reasonably well by the parameterization not only with multiple wave propagation directions but also with two wave propagation directions of 45 degrees (northeast-southwest) and 135 degrees (northwest-southeast), which are optimally chosen for computational efficiency.

convective gravity wave drag↗

Hot Structure Concepts and Materials for Space Exploration Project

Hot structures are multifunctional material systems that integrate the functions of thermal protection systems (TPS) and primary structure to increase mass efficiency in aerospace vehicles. The goal of this project was to improve our capabilities for developing lightweight, reusable, and reliable hot structures concepts to enable advanced space exploration missions. Specifically, current hot structure carbon/carbon (C/C) composite material systems have low interlaminar strength, which can lead to material failures at low load levels, and large material property variability in as-fabricated parts. Two parallel paths were followed to address these shortcomings. First, analysis models were developed to identify the conditions under which excessive interlaminar stresses occur and to evaluate the potential of stitching the C/C laminae to improve the interlaminar thermal-structural response. This effort included the development of high-fidelity, fracture-based thermal-structural finite element analysis models to evaluate critical features needed to produce efficient hot structure concepts. The second path focused on developing atomistic and mesoscale simulation tools to model the fabrication process in an effort to reduce the time and cost of production while improving the quality of the resulting parts. Aspects of the process currently being investigated include the infusion of the matrix precursor into the fiber preform, the cure of the precursor molecules into their polymeric form, and the carbonization and graphitization of the polymer to form the carbon matrix material. Results of the 6-month study are presented herein.

thermal↗

Mesoscale Convective Systems Simulated by a High–Resolution Global Nonhydrostatic Model Over the United States and China

Mesoscale convective systems (MCSs) contribute a large fraction of warm-season precipitation and generate hazardous weather with substantial socio-economic impacts. Uncertainties in convection parameterizations in climate models limit our understanding of MCS characteristics and reliability of future projection. We examine the MCS simulation from the global 14-km Nonhydrostatic ICosahedral Atmospheric Model (NICAM) without cumulus parameterization against satellite observation from Global Precipitation Measurement (GPM) during 2001-2008. We focus on MCSs over the central U.S. and eastern China where MCSs prevail from March to August. A process-oriented tracking method incorporating both cloud and precipitation criteria is used to identify and track MCSs. About 140/100 MCSs initiate in the central U.S./eastern China per warm season and most of them initiate east of high mountains and in coastal regions. The frequency distribution of MCS lifetime is well captured in NICAM. But the simulated MCSs have stronger precipitation, smaller precipitation area, and larger cold cloud system than observed in both regions, which may be caused by weak entrainment as it is not well resolved at 14 km resolution. Here, the simulated MCS number is also underestimated in summer. By examining the climatological and MCS large-scale environments, the significant underestimation of MCS number in summer over the central U.S. may be attributed to the large climatological dry bias in the atmosphere. For China, mean moisture in summer is well simulated but deficiency in capturing the dynamic condition related to the coastal topography for triggering convection may have contributed to underestimation of MCS even in a sufficiently moist environment.

54 ENVIRONMENTAL SCIENCES↗

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↗

The interactive role of subsynoptic scale jet sreak and planetary boundary layer adjustments in organizing an apparently isolated convective complex

A mesoscale atmospheric simulation system is described that is being developed in order to improve the simulation of subsynoptic and mesoscale adjustments associated with cyclogenesis, severe storm development, and significant atmospheric transport processes. Present emphasis in model development is in the parameterization of physical processes, time-dependent boundary conditions, sophisticated initialization and analysis procedures, nested grid solutions, and applications software development. Basic characteristics of the system as of March 1982 are listed. In a case study, the Grand Island tornado outbreak of 3 June 1980 is considered in substantial detail. Results of simulations with a mesoscale atmospheric simulation system indicate that over the high plains subtle interactions between existing jet streaks and deep well mixed boundary layers can lead to well organized patterns of mesoscale divergence and pressure falls. The amplitude and positioning of these mesoscale features is a function of the subtle nonlinear interaction between the pre-existing jet-streak and deep well mixed boundary layers. Model results for the case study indicate that the model has the potential for forecasting the precursor mesoscale convective environment.

Kaplan, M. L.↗

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↗

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↗

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↗

Microphysical modeling of cirrus. 1: Comparison with 1986 FIRE IFO measurements

We have used a one-dimensional model of cirrus formation to study the development of cirrus clouds during the 1986 First ISCCP (International Satellite Cloud Climatology Project) Regional Experiment (FIRE) intensive field observations (IFO). The cirrus model includes microphysical, dynamical, and radiative processes. Sulfate aerosols, solution drops, ice crystals, and water vapor are all treated as interactive elements in the model. Ice crystal size distributions are fully resolved based on calculations of homogeneous freezing nucleation, growth by water vapor deposition, evaporation, coagulation, and vertical transport. We have focused on the cirrus observed on November 1, 1986. Vertical wind speed for the one-dimensional simulation is taken from a mesoscale model simulation for the appropriate time period. The mesoscale model simulation suggested that strong upward motions over Wyoming and subsequent horizontal transport of upper level moisture were responsible for the cirrus observed over Wisconsin on this date. We assumed that our one-dimensional model could be used to represent a vertical column moving from Wyoming to Wisconsin over a period of several hours. Ice crystal nucleation occurs in our model in the 8 to 10-km region as a result of the strong updrafts (and cooling) early in the simulation. Growth, coagulation, and sedimentation of these ice crystals result in a broad cloud region (5-10 km thick) with an optical depth of 1-2 after a few hours, in agreement with the FIRE measurements. Comparison with aircraft microphysical measurements made over Wisconsin indicates that the simulation generated reasonable ice water content, but the predicted ice number densities are too low, especially for radii less than about 50 microns. Sensitivity tests suggest that better agreement between simulated and observed microphysical properties is achieved if the nucleation rate is higher or stronger vertical mixing (perhaps associated with multidimensional motions) is present.

Jensen, Eric J.↗

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↗