Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesoscale modeling”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An energy-based coupling approach to nonlocal interface problems.

Nonlocal models provide accurate representations of physical phenomena ranging from fracture mechanics to complex subsurface flows, settings in which traditional partial differential equation models fail to capture effects caused by long-range forces at the microscale and mesoscale. However, the application of nonlocal models to problems involving interfaces, such as multimaterial simulations and fluid-structure interaction, is hampered by the lack of a physically consistent interface theory which is needed to support numerical developments and, among other features, reduces to classical models in the limit as the extent of nonlocal interactions vanish. In this paper, we use an energy-based approach to develop a formulation of a nonlocal interface problem which provides a physically consistent extension of the classical perfect interface formulation for partial differential equations. Numerical examples in one and two dimensions validate the proposed framework and demonstrate the scope of our theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Mesoflow: An Open-Source Reacting Flow Solver for Catalysis at Mesoscale

We present the capabilities and software performance metrics of our open-source continuum solver for catalysis, Mesoflow, developed specifically for modeling transport and chemistry at the mesoscale. Our solver utilizes Cartesian block-structured adaptive mesh refinement to resolve complex catalyst surface morphologies directly obtained from X-ray tomography data. An immersed boundary based formulation enables rapid representation of complex geometries prevalent in most mesoporous catalyst interfaces. The solver is developed on top of open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Our flexible software framework enables integration of complex chemical mechanisms at heterogenous interfaces and time-split algorithms for circumventing highly disparate reaction and flow time-scales. Our current studies indicate a ten-fold performance gain by using graphics-processing-units (GPUs) compared to a single processor for representative problem sizes (2 million cell mesh). We will also present a brief introduction on how to build and use this software for application problems pertaining to catalytic upgrading and gas transport within porous catalyst particles.

adaptive meshing↗
DOE OSTI2020

Art and Science of the Cellular Mesoscale

Experimental information from microscopy, structural biology, and bioinformatics may be integrated to build structural models of entire cells with molecular detail. This integrative modeling is challenging in several ways: the intrinsic complexity of biology results in models with many closely packed and heterogeneous components; the wealth of available experimental data is scattered among multiple resources and must be gathered, reconciled, and curated; and computational infrastructure is only now gaining the capability of modeling and visualizing systems of this complexity. We present recent efforts to address these challenges, both with artistic approaches to depicting the cellular mesoscale, and development and application of methods to build quantitative models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Unpacking model inadequacy: The quantification of silver release from TRISO fuel by considering empirical and mechanistic approaches

Increasing adoption of the proposed tristructural isotropic (TRISO) particle fuel for both advanced and existing reactors makes it critical to assess and address any uncertainties and inadequacies of TRISO fission product release models. Model inadequacy stems from simplifications made to the computational model when compared to the experiments. The modeling and simulation efforts conducted using the BISON fuel performance code, along with the experimental campaigns carried out under the Advanced Gas Reactor Fuel Development and Qualification Program, afford a unique opportunity to conduct a rigorous modeling inadequacy assessment within the Bayesian uncertainty quantification (UQ) framework. Here, this study compares the standard Bayesian framework against the Kennedy-O'Hagan (KOH) framework, which explicitly represents modeling inadequacy, in regard to UQ for TRISO silver release models. For this purpose, both the traditional Arrhenius equation fitted to experimental data and the more advanced lower-length-scale (LLS)-informed model, which considers microstructure information, are independently considered. Applying the inverse UQ process on the AGR-2 and -3/4 datasets revealed modeling inadequacy to be the most dominant source of uncertainty. Experimental noise uncertainty is also significant; however, model parameter uncertainty can be considered negligible. Interestingly, both the Arrhenius equation and the LLS-informed model demonstrated similar levels of modeling inadequacy. For the forward predictive UQ, the KOH framework improved both the accuracy and quality of quantified uncertainties in comparison to the standard Bayesian framework. This is true for both the Arrhenius equation and the LLS-informed model. In comparing these modeling approaches, both demonstrated similar performance at the engineering scale, while the LLS-informed model expectedly outperformed the Arrhenius equation at the mesoscale. These conclusions highlight the importance of explicitly accounting for modeling inadequacy in the UQ process, and reinforce the need for continuous refinement of physics-based models in order to address the modeling inadequacy.

Advanced reactors↗
DOE OSTI2023

Grid Spacing Sensitivities of Simulated Mid-Latitude and Tropical Mesoscale Convective Systems in the Convective Gray Zone

The main objective of this study is to observationally constrain processes in tropical and midlatitude mesoscale convective systems (MCSs), and to use these constraints for model evaluation. To accomplish this, we leverage MCS observations collected at the U.S. DOE Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site in Oklahoma and ARM's mobile GoAmazon2014/15 site in Manaus, Brazil (MAO). We simulate 13 and 11 of these observed MCSs at the SGP and MAO site, respectively, using the Weather Research and Forecasting model at 12-, 4-, 2-, and 1-km horizontal grid spacing. Observations from radiosondes, surface meteorology, and radar wind profilers are used to characterize MCS properties, such as MCS timing and location, cold pools, and convective drafts, and evaluate these simulations. SGP cases are found in better agreement with observations than MAO cases, and when simulated at 2 km, outperform simulations at 1 km regarding the timing of MCS overpass and the accuracy of surface variable trends. MAO simulations suggest a consistent improvement in model accuracy with increasing model resolution in depicting the downdraft structure, the timing of MCSs, and the surface variables changes, except for the latter two metrics at 2 km. Deficiencies are still evident at km-scales, suggesting the need for higher resolution to simulate tropical MCSs. Overall, location-dependent improvements in MCS representation are obtained with the increasing model resolution, prompting the evaluation of sub-km scale simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A mechanistic, multiscale model for predicting Pd penetration in TRISO fuels using BISON

TRistructural ISOtropic (TRISO) particles use silicon carbide (SiC) as the primary structural member and barrier against metallic fission product (FP) release. palladiums (PDs), produced by fission in the fuel kernel, can diffuse to and chemically interact with the SiC layer, degrading its structural integrity and ability to contain radioactive FPs. Existing temperature-dependent correlations for predicting Pd penetration rely on experimental data with significant scatter due to varying conditions, potentially complicating ongoing fuel qualification and licensing efforts for advanced reactors that would subject TRISO fuels to operating conditions outside of those examined in the experiments. A mechanistic model of Pd production, transport, and reaction is developed in this work to better understand and predict PD attack of SiC in TRISO particles. molecular dynamicss (MDs) simulations are utilized to calculate the Pd diffusivity in SiC grain bulk and grain boundaries. A mesoscale phase-field diffusion model, informed by the MD diffusivities, is used to develop a reduced order model (ROM) for the effect of SiC microstructure on the Pd penetration rate. The engineering scale BISON model calculates Pd production and transport, utilizing the ROM to predict penetration rates consistent with experimental data. This novel mechanistic ROM captures the effect of temperature, microstructure, and irradiation history on Pd penetration. In conclusion, these new capabilities are expected to support ongoing qualification and licensing efforts associated with near-term TRISO-fueled reactor applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Summertime Near-Surface Temperature Biases Over the Central United States in Convection-Permitting Simulations

Convection-Permitting Model (CPM) simulations of the Central United States climate for the summer of 2011 are studied to understand the causes of warm biases in 2-m air temperature (T 2m ) and related underestimates of precipitation including that from mesoscale convective systems (MCSs). Based on 10 CPM simulations and 9 coarser-resolution model simulations, we quantify contributions from evaporative fraction (EF) and radiation to the T 2m bias with both types of models overestimating T 2m largely because they underestimate EF. The performance of CPMs in capturing MCS characteristics (frequency, rainfall, propagation) varies. The pre-summer precipitation bias has large correlation with mean summertime T 2m bias but the relationship between summertime MCS mean rainfall bias and T 2m bias is non-monotonic. Analysis of lifting condensation level deficit and convective available potential energy suggests that models with T 2m warm biases and low EF have too dry and stable boundary layers, inhibiting the formation of clouds, precipitation and MCSs. Among the CPMs with differing model formulations (e.g., transpiration, infiltration, cloud macrophysics and microphysics), evidence suggests that altering the land-surface model is more effective than altering the atmospheric model in reducing T 2m biases. In conclusion, these results demonstrate that land-atmosphere interactions play a very important role in determining the summertime climate of the Central United States.

2-m air temperature↗
DOE OSTI2020

Volumetric analysis and mesh generation of real and artificial microstructural geometries

Producing a viable finite element mesh of realistic microstructural structural geometry is a critical step in analyzing the thermo-mechanical behavior of complex multi-material composites. Advancements in imaging technology such as micro computed tomography have allowed modelers to access high resolution mesoscale geometries for direct numerical simulation. However, converting from voxel based 3D images to usable finite element meshes has been challenging. A robust method including algorithms and software scripts for generating finite element meshes from 3D imaged microstructures is presented in this paper. It includes a routine for inserting cohesive elements around material interfaces to enable modeling of interface properties including delamination and damage. The algorithms and procedures presented in this method leverage currently available software packages for processing surface based geometry into volume based meshes. In addition to converting real geometry from physical imaging systems, algorithms for producing numerically generated and statistically equivalent microstructural geometry are also included. These artificial microstructures can be a valuable resource for modelers when physical specimens do not exist or are limited in quantity: Method establishes a workflow from voxel data to viable finite element mesh including interface information, and Includes a method for synthetic geometry generation based on metrics from real microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Interactions Between Clouds and Wind-Driven Surface Heat Exchanges over Land

Earth system model experiments show that increasing horizontal resolution fundamentally alters the simulated soil-moisture-precipitation feedback. Kilometer-scale simulations often produce weaker or even negative feedback compared to coarse-resolution models. A key difference of kilometer-scale models is that they resolve mesoscale secondary circulations, including boundary layer horizontal rolls and cellular structures, in addition to cold pools and downdrafts associated with convective precipitation. However, because the relevant processes occur on yet-smaller scales, these circulations are often poorly resolved. This project demonstrated that boundary layer secondary circulations significantly affect surface heat exchanges and wind gusts, and that current model parameterizations can misrepresent these processes at kilometer-scale resolution. Using DOE Atmospheric Radiation Measurement (ARM) observations and targeted experiments with the DOE Energy Exascale Earth System Model (E3SM), we identified physically unrealistic wind gust and surface flux responses to secondary circulations, diagnosed a systematic overestimation of wind shear in convective cold pools, and uncovered a multivariate relationship between land surface fluxes and the scales of updrafts that form shallow cumulus clouds. These findings provide observation-based recommendations for improving parameterizations of surface fluxes and wind gusts in high-resolution Earth system models, thereby reducing uncertainty in convective storm prediction and land-atmosphere feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

PyFLEXTRKR: a flexible feature tracking Python software for convective cloud analysis

Abstract. This paper describes the new open-source framework PyFLEXTRKR (Python FLEXible object TRacKeR), a flexible atmospheric feature tracking software package with specific capabilities to track convective clouds from a variety of observations and model simulations. This software can track any atmospheric 2D objects and handle merging and splitting explicitly. The package has a collection of multi-object identification algorithms, scalable parallelization options, and has been optimized for large datasets including global high-resolution data. We demonstrate applications of PyFLEXTRKR on tracking individual deep convective cells and mesoscale convective systems from observations and model simulations ranging from large-eddy resolving (∼100s m) to mesoscale (∼10s km) resolutions. Visualization, post-processing, and statistical analysis tools are included in the package. New Lagrangian analyses of convective clouds produced by PyFLEXTRKR applicable to a wide range of datasets and scales facilitate advanced model evaluation and development efforts as well as scientific discovery.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Efficacy of the Cell Perturbation Method in Large-Eddy Simulations of Boundary Layer Flow over Complex Terrain

A challenge to simulating turbulent flow in multiscale atmospheric applications is the efficient generation of resolved turbulence motions over an area of interest. One approach is to apply small perturbations to flow variables near the inflow planes of turbulence-resolving simulation domains nested within larger mesoscale domains. While this approach has been examined in numerous idealized and simple terrain cases, its efficacy in complex terrain environments has not yet been fully explored. Here, we examine the benefits of the stochastic cell perturbation method (CPM) over real complex terrain using data from the 2017 Perdigão field campaign, conducted in an approximately 2-km wide valley situated between two nearly parallel ridges. Following a typical configuration for multiscale simulation using nested domains within the Weather Research and Forecasting (WRF) model to downscale from the mesoscale to a large-eddy simulation (LES), we apply the CPM on a domain with horizontal grid spacing of 150 m. At this resolution, spurious coherent structures are often observed under unstable atmospheric conditions with moderate mean wind speeds. Results from such an intermediate resolution grid are often nested down for finer, more detailed LES, where these spurious structures adversely affect the development of turbulence on the subsequent finer grid nest. We therefore examine the impacts of the CPM on the representation of turbulence within the nested LES domain under moderate mean flow conditions in three different stability regimes: weakly convective, strongly convective, and weakly stable. In addition, two different resolutions of the underlying terrain are used to explore the role of the complex topography itself in generating turbulent structures. We demonstrate that the CPM improves the representation of turbulence within the LES domain, relative to the use of high-resolution complex terrain alone. During the convective conditions, the CPM improves the rate at which smaller-scales of turbulence form, while also accelerating the attenuation of the spurious numerically generated roll structures near the inflow boundary. During stable conditions, the coarse mesh spacing of the intermediate LES domain used herein was insufficient to maintain resolved turbulence using CPM as the flow develops downstream, highlighting the need for yet higher resolution under even weakly stable conditions, and the importance of accurate representation of flow on intermediate LES grids.

17 WIND ENERGY↗
DOE OSTI2021

Decadal Variability of Southeast US Rainfall in an Eddying Global Coupled Model

Abstract Ocean variability is a dominant source of remote rainfall predictability, but in many cases the physical mechanisms driving this predictability are not fully understood. This study examines how ocean mesoscales (i.e., the Gulf Stream SST front) affect decadal Southeast US (SEUS) rainfall, arguing that the local imprint of large‐scale teleconnections is sensitive to resolved mesoscale features. Based on global coupled model experiments with eddying and eddy‐parameterizing ocean, we find that a resolved Gulf Stream improves localized rainfall and remote circulation response in the SEUS. The eddying model generally improves the air‐sea interactions in the Gulf Stream and the North Atlantic Subtropical High that modulate SEUS rainfall over decadal timescales. The eddy‐parameterizing simulation fails to capture the sharp SST gradient associated with the Gulf Stream and overestimates the role of tropical Pacific SST anomalies in the SEUS rainfall.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Monsoonal MCS Initiation, Rainfall, and Diurnal Gravity Waves over the Bay of Bengal: Observation and a Linear Model

Abstract Previous observational studies have indicated that mesoscale convective systems (MCSs) contribute the majority of precipitation over the Bay of Bengal (BoB) during the summer monsoon season, yet their initiation and propagation remain incompletely understood. To fill this knowledge gap, we conducted a comprehensive study using a combination of 20-yr satellite observations, MCS tracking, reanalysis data, and a theoretical linear model. Satellite observations reveal clear diurnal propagation signals of MCS initiation frequency and rainfall from the west coast of the BoB toward the central BoB, with the MCS rainfall propagating slightly slower than the MCS initiation frequency. Global reanalysis data indicate a strong association between the offshore-propagating MCS initiation frequency/rainfall and diurnal low-level wind perturbations, implying the potential role of gravity waves. To verify the hypothesis, we developed a 2D linear model that can be driven by realistic meteorological fields from reanalysis. The linear model realistically reproduces the characteristics of offshore-propagating diurnal wind perturbations. The wind perturbations, as well as the offshore propagation signals of MCS initiation frequency and rainfall, are associated with diurnal gravity waves emitted from the coastal regions, which in turn are caused by the diurnal land–sea thermal contrast. The ambient wind speed and vertical wind shear play crucial roles in modulating the timing, propagation, and amplitude of diurnal gravity waves. Using the linear model and satellite observations, we further show that the stronger monsoonal flows lead to faster offshore propagation of diurnal gravity waves, which subsequently control the offshore propagation signals of MCS initiation and rainfall. Significance Statement Rainfall over the Bay of Bengal (BoB) is primarily contributed by large and organized rainfall systems in the summer monsoon season. During this season, these systems are commonly observed over the east coast of India around midnight, the western BoB in the morning, and the central BoB in the afternoon. This eastward rainfall propagation is confirmed by observations, reanalysis data, and a theoretical model to have a strong association with atmospheric diurnal gravity waves. These waves are caused by land–sea thermal contrast and can trigger rainfall systems over the offshore regions. We also found that the diurnal gravity waves, as well as their triggered rainfall systems, can be greatly modulated by the large-scale monsoonal flows. All the above findings improve our understanding of diurnal rainfall cycle over the BoB and may contribute to the future improvement of rainfall forecast over the region.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Enabling Particulate Materials Processing Science for High-Consequence, Small-Lot Precision Manufacturing

This Laboratory Directed Research and Development project developed and applied closely coupled experimental and computational tools to investigate powder compaction across multiple length scales. The primary motivation for this work is to provide connections between powder feedstock characteristics, processing conditions, and powder pellet properties in the context of powder-based energetic components manufacturing. We have focused our efforts on multicrystalline cellulose, a molecular crystalline surrogate material that is mechanically similar to several energetic materials of interest, but provides several advantages for fundamental investigations. We report extensive experimental characterization ranging in length scale from nanometers to macroscopic, bulk behavior. Experiments included nanoindentation of well-controlled, micron-scale pillar geometries milled into the surface of individual particles, single-particle crushing experiments, in-situ optical and computed tomography imaging of the compaction of multiple particles in different geometries, and bulk powder compaction. In order to capture the large plastic deformation and fracture of particles in computational models, we have advanced two distinct meshfree Lagrangian simulation techniques: 1.) bonded particle methods, which extend existing discrete element method capabilities in the Sandia-developed , open-source LAMMPS code to capture particle deformation and fracture and 2.) extensions of peridynamics for application to mesoscale powder compaction, including a novel material model that includes plasticity and creep. We have demonstrated both methods for simulations of single-particle crushing as well as mesoscale multi-particle compaction, with favorable comparisons to experimental data. We have used small-scale, mechanical characterization data to inform material models, and in-situ imaging of mesoscale particle structures to provide initial conditions for simulations. Both mesostructure porosity characteristics and overall stress-strain behavior were found to be in good agreement between simulations and experiments. We have thus demonstrated a novel multi-scale, closely coupled experimental and computational approach to the study of powder compaction. This enables a wide range of possible investigations into feedstock-process-structure relationships in powder-based materials, with immediate applications in energetic component manufacturing, as well as other particle-based components and processes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modeling shock-induced void collapse in single-crystal Ta systems at the mesoscales

Understanding the role of microstructural heterogeneities on the shock wave propagation and defect evolution behavior is essential to predicting the dynamic response of metals. Heterogeneities, such as voids, provide challenges to understanding the wave propagation behavior as the shock-void interaction can collapse the void and result in large plastic strains and significant heating (hot spot formation) in the metal. Accurate modeling of this phenomenon requires predicting the void collapse mechanisms and the related heat generation and dissipation mechanisms (hotspot formation) that determine the microstructure evolution. While molecular dynamics (MD) simulations can model the void collapse behavior, the time/length scale capabilities pose a challenge to connect with continuum models or the experimental scales. Here, this study presents the capability of the newly developed quasi-coarse-grained dynamics (QCGD) method that extends the MD simulations to larger system sizes and longer times to model this phenomenon. This study uses QCGD simulations to investigate the mechanisms of shock wave interactions with pre-existing voids in single-crystal Ta microstructures for variations in shock pressures, void size, and loading orientations. For a given orientation, the plasticity contributions and rates of void collapse are observed to vary with shock pressures and void size. The larger void sizes and higher pressures result in increased temperatures (hot spots) and faster void collapse rates and unravel the variations in the plasticity contributions. In addition, QCGD simulations investigate the post-collapse microstructure evolution as a release wave travels through the hot spot region. The simulations reveal that the reduction in temperatures due to heat dissipation initiates the dynamic recrystallization behavior in the hotspot regions.

36 MATERIALS SCIENCE↗