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 325 records · Page 18

Learning two-phase microstructure evolution using neural operators and autoencoder architectures

Abstract Phase-field modeling is an effective but computationally expensive method for capturing the mesoscale morphological and microstructure evolution in materials. Hence, fast and generalizable surrogate models are needed to alleviate the cost of computationally taxing processes such as in optimization and design of materials. The intrinsic discontinuous nature of the physical phenomena incurred by the presence of sharp phase boundaries makes the training of the surrogate model cumbersome. We develop a framework that integrates a convolutional autoencoder architecture with a deep neural operator (DeepONet) to learn the dynamic evolution of a two-phase mixture and accelerate time-to-solution in predicting the microstructure evolution. We utilize the convolutional autoencoder to provide a compact representation of the microstructure data in a low-dimensional latent space. After DeepONet is trained in the latent space, it can be used to replace the high-fidelity phase-field numerical solver in interpolation tasks or to accelerate the numerical solver in extrapolation tasks.

36 MATERIALS SCIENCE↗

Resolved gravity waves in the tropical stratosphere: Impact of horizontal resolution and deep convection parametrization

Convectively generated gravity waves (CGGWs) are important for numerical weather prediction due to their effect on the quasi-biennial oscillation (QBO) in the stratosphere. Here, using global ECMWF IFS simulations at TCo7999 (or 1.25 km), TCo2559 (or 3.9 km) and TCo1279 (or 7.8 km) horizontal resolutions, sensitivity of resolved CGGWs to the horizontal resolution and to the explicit versus parametrized representation of deep convection is elucidated during the westerly shear phase of the QBO. Parametrized deep convection is found to inhibit CGGWs, resulting in a twofold reduction in CGGW forcing. When deep convection is explicitly resolved, the total CGGW forcing is almost unchanged across the horizontal resolutions. However, the contribution of long and mesoscale CGGWs (with horizontal wavelengths 100 km ≤ $λ_h$ < $1,900$ km) to the total CGGW forcing decreases and the contribution of smaller-scale CGGWs (with $λ_h$ < $100$ km) increases as the horizontal resolution increases. At the maximum CGGW forcing altitude, at TCo7999 resolution 43% of the total CGGW forcing is due to long and mesoscale waves, whereas at TCo2559 and TCo1279 resolutions their contribution is 70% and 90%, respectively. While CGGW forcing from long and mesoscale waves is similar at TCo7999 resolution with explicit deep convection and at TCo1279 resolution with parametrized deep convection, CGGW forcing from these waves is artificially enhanced at TCo1279 and TCo2559 resolutions with explicit deep convection. This is due to the explicit deep convection being too strong and having too much variance for 100 km ≤ $λ_h$ < $1,900$ km. Therefore, parametrizations of deep convection and CGGWs (to account for forcing from waves with $λ_h$ < $100$ km) are required even at TCo2559 resolution. Additionally, resolved CGGW forcing at TCo7999 resolution is examined for the easterly shear phase of the QBO; similar to the westerly shear phase, the smaller-scale waves contribute 55% to the total CGGW forcing at the maximum CGGW forcing altitude.

54 ENVIRONMENTAL SCIENCES↗

The Impacts of Rotational Mixing on the Precipitation Simulated by a Convection Permitting Model

With increased availability of computational resources, regional and global scale convection-permitting model (CPM, Δx ~ 1–10 km) simulations are becoming more common. CPMs have improved accuracy in their representation of deep convection and mesoscale convective systems (MCSs) compared to coarser resolution models. However, CPMs still exhibit convective cloud and precipitation biases relative to observations, notably a lesser frequency of light precipitation rates and greater frequency of heavy precipitation rates. In this work we hypothesize that these CPM biases are related to under-resolved mixing between convective updrafts and their surrounding environment. To test this hypothesis, we introduce a parameterization to the Weather Research and Forecasting model (WRF) that adds a small angular rotation of the grid-scale flow about the axis perpendicular to the plane of convective drafts. This rotated flow is then allowed to alter advection of moisture and hydrometeors. The effects of such mixing on precipitation characteristics are evaluated in month-long 4-km grid spacing simulations over the Amazon. The enhanced mixing transports moisture and condensate from convective cores to other areas including downdrafts. This increases the frequency of low-precipitable water and light precipitation. It also decreases the frequency of intense precipitation from isolated deep convection and MCSs, increases cloud top temperatures, reduces radar echo-top heights, and increases overall precipitation by altering the relationship of precipitation with precipitable water, in better agreement with observations. The results suggest when optimized using multiple observations, such an approach may provide a path toward more accurate representation of convection and precipitation statistics in convection-permitting simulations.

54 ENVIRONMENTAL SCIENCES↗

Core energies of dislocations in bcc metals

Accurate methods and an efficient workflow for computing and documenting dislocation core energies are developed and applied to 1 2 ( 111 ) and ( 100 ) dislocations in five body-centered cubic (bcc) metals W, Ta, V, Mo, and α -Fe represented by 13 model interatomic potentials. For each dislocation type, dislocation core energies are extracted for a large number of dislocation characters thoroughly sampling the entire 2-space of crystallographic line orientations of the bcc lattice. Of particular interest, core energies of the 1 2 ( 111 ) { 110 } dislocations are found to be distinctly asymmetric with respect to the sign of the character angle, whereas core energies of ( 100 ) { 110 } junction dislocations exhibit marked cusps for line orientations vicinal to the closed-packed ( 111 ) directions. Overall, our findings furnish substantial insights for developing accurate models of dislocation core energies employed in mesoscale dislocation dynamics simulations of crystal plasticity.

36 MATERIALS SCIENCE↗

The Role of Deep Convection and Large-scale Circulation in Driving Model Spread in Low Cloud Feedback and Equilibrium Climate Sensitivity

This project aims to advance the understanding of the processes that drive the large uncertainties in climate change projections, use observations to constrain model physics and reduce the inter-model spread in equilibrium climate sensitivity (ECS). There are three major goals: 1) Characterize the representation of the physical pathways that link deep convection, large-scale circulation and low cloud feedback in CMIP6 model simulations and determine the relative contribution of each pathway to the CMIP6 model spread in low cloud feedback and ECS; 2) Use process-oriented diagnostics and multiple observations to evaluate CMIP6 model performance in capturing the observed cloud-circulation relation and deep convection characteristics including convective transition statistics and the bulk properties of mesoscale convective systems (MCSs). Error decomposition in CMIP6 models will be performed. 3) Conduct E3SM short-range hindcasts following the DOE Cloud-Associated Parameterizations Testbed (CAPT) protocol to pinpoint specific model parameters/processes that are crucial to the representation of deep convection, circulation, clouds and the pathways that connect them. We will modify convective parameters in E3SM and analyze the perturbed physics experiments (PPEs) to isolate model parameters that are critical to the uncertainty of ECS.

54 ENVIRONMENTAL SCIENCES↗

Dislocation transmission across Σ3{112} incoherent twin boundary: a combined atomistic and phase-field study

Grain boundaries (GBs) in polycrystalline materials act as impediments to dislocation motion and result in strengthening. Understanding slip transmission through GBs, specifically twin boundaries, is essential to understand the plastic deformation behavior of polycrystalline fcc materials. Here the interaction between a glide dislocation and Σ3{112} incoherent twin boundary (ITB) in copper is investigated using a combined atomistic and mesoscale approach. The material parameters and structure of the GB in the mesoscale phase field dislocation dynamics (PFDD) model are informed from Molecular Statics (MS) simulations. The structural unit of the ITB consists of an array of three partial dislocations. The interaction between a glide dislocation impinging on each of the GB partial dislocations is investigated using both PFDD and Molecular Dynamics (MD) with two boundary conditions. Transmission planes predicted by both PFDD and MD (NVT) are in agreement, and show that not all transmission events are direct. Critical transmission stresses predicted by PFDD are in the range of 276 MPa to 1380 MPa, while MD predictions are in the range from 100 MPa to 700 MPa. The PFDD and MD predictions of slip transmission are explained using dislocation theory based on isotropic linear elasticity.

36 MATERIALS SCIENCE↗

A improved equation of state for Xe gas bubbles in gamma U-Mo fuels

A monolithic fuel design based on a U-Mo alloy has been selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). An issue with U-Mo monolithic fuel is the large amount of swelling that takes place during operation. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic properties into mesoscale and continuum level fuel performance modeling codes. However, the thermodynamic properties of the fission gas bubbles (such as the relationship among bubble size, equilibrium Xe concentration, and bubble pressure) are not well known. This work studies Xe bubbles in ?U-Mo from a diameter of 3 nm up to 8.5 nm and from 400 K up to 700 K. The energetic relationship of Xe bubbles with regard to voids and Xe substitutional atoms is described. The transition is also determined for when a bubble becomes over-pressurized. Finally, an equation of state is fit to the pressure as a function of molar volume and temperature.

Beeler, Benjamin↗

Radiation driven diffusion in γU-Mo

A monolithic fuel design based on a U-Mo alloy has been selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). A critical phenomenon of interest with U-Mo monolithic fuel is the large amount of swelling that takes place during operation, particularly at high fission densities. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic and kinetic properties into mesoscale and continuum level fuel performance modeling codes. One such property where there exists incomplete data is the diffusion of relevant species under irradiation. Fuel performance swelling predictions rely on an accurate representation of diffusion in order to determine the rate of fission gas swelling and the local microstructural evolution. In this work, we present molecular dynamics simulations of the radiation driven diffusion of U, Mo and Xe in U-Mo nuclear fuels. Diffusion coefficients for each species are determined over a range of temperatures and compositions. In this work, updated diffusion coefficients are presented that are applicable under irradiation that incorporate both intrinsic and radiation driven diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Ocean fronts and eddies force atmospheric rivers and heavy precipitation in western North America

Atmospheric rivers (ARs) are responsible for over 90% of poleward water vapor transport in the mid-latitudes and can produce extreme precipitation when making landfall. However, weather and climate models still have difficulty simulating and predicting landfalling ARs and associated extreme precipitation, highlighting the need to better understand AR dynamics. Here, using high-resolution climate models and observations, we demonstrate that mesoscale sea-surface temperature (SST) anomalies along the Kuroshio Extension can exert a remote influence on landfalling ARs and related heavy precipitation along the west coast of North America. Inclusion of mesoscale SST forcing in the simulations results in approximately a 40% increase in landfalling ARs and up to a 30% increase in heavy precipitation in mountainous regions and this remote impact occurs on two-week time scales. The asymmetrical response of the atmosphere to warm vs. cold mesoscale SSTs over the eddy-rich Kuroshio Extension region is proposed as a forcing mechanism that results in a net increase of moisture flux above the planetary boundary layer, prompting AR genesis via enhancing moisture transport into extratropical cyclones in the presence of mesoscale SST forcing.

54 ENVIRONMENTAL SCIENCES↗

Thermal and electric field driven rf breakdown precursor formation on metal surfaces

The phenomenon of electric breakdown poses serious challenges to the design of devices that operate in high electric field environments. Experimental evidence points toward breakdown events that are accompanied by elevated temperatures and dark current spikes, which is attributed to high-asperity nanostructure formation that enhances the local electric field and triggers a runaway process. However, the exact mechanistic origin of such nanostructures under typical macroscopic operational conditions of electric field and magnetic-field-mediated heating remains poorly understood. In this work, we simulate the evolution of a copper surface under the combined action of the electric fields and elevated temperatures. Using a mesoscale curvature-driven surface evolution model, we show how a copper surface can undergo a type of dynamical instability that naturally leads to the formation of sharp asperities in realistic experimental conditions. Exploring the combined effect of fields and temperature rise, we identify the critical regimes that allow for the formation of breakdown precursors. The results show that thermoelastic stresses, while not essential, can significantly lower the critical electric field required for runaway surface instability, which is consistent with experimental observations that thermal effects can increase breakdown rates. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

mesoflow [SWR-22-56]

Mesoflow is a continuum scale simulation tool developed specifically for modeling transport and chemistry at the mesoscale. Our solver utilizes Cartesian block-structured adaptive mesh refinement to resolve complex surface morphologies (of catalysts/biomass particles among others) directly obtained from X-ray tomography data. An immersed boundary based formulation enables rapid representation of complex geometries prevalent in most mesoporous 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 (GPU) compared to a single processor for representative problem sizes (2 million cell mesh).

Sitaraman, Hariswaran↗

Phase-field modeling of radiation-induced segregation for multicomponent alloys

Structural alloys under irradiation are known to undergo radiation-induced solute redistribution (RIS) at grain boundaries, leading to detrimental effects such as intergranular corrosion and stress-assisted cracking. To better understand the phenomenon, improved models of RIS applicable to concentrated, multicomponent alloys, and mesoscale microstructures are needed. In this talk, we present a novel grand-potential-based phase-field model to account for the complete set of multicomponent kinetic and thermodynamic couplings between atoms and point defects in the Onsager transport equations. We demonstrate multiscale modeling capability by deriving the Onsager coefficient matrix from atomistic-based Kinetic Monte Carlo simulations. Model predictions and validations of RIS and the effect of defect production, grain boundary sink strength and density are demonstrated for a model FCC FeCrNi system. We also demonstrate the novel capability to describe RIS in the presence of equilibrium segregation described using a density-based CALPHAD thermodynamic framework. This modeling approach overcomes certain limitations in conventional RIS models and provides a step closer towards multiscale modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Simulating Catalysis with Realistic Pellet Geometries Using Mesoflow: A Case Study of Catalytic Propane Dehydrogenation

We present a case study of catalytic propane dehydrogenation with our open-source multiphysics solver, Mesoflow. The solver was developed to simulate reactive flow coupled to heterogeneous catalytic reactions and deactivation in the context of complex, mesoscale geometry. The method leverages cartesian block-structured adaptive mesh refinement to capture realistic catalyst microstructural features acquired directly from X-ray computed tomography data. A kinetic model for propane dehydrogenation and catalyst deactivation was developed based on temporal analysis of products (TAP) reactor experiments. The TAP reactor experiments allow for precise characterization of intrinsic kinetic reaction steps which are implemented into Mesoflow simulations to model the spatial and temporal evolution of reactants, products, and catalyst active sites. The short-term and long-term deactivation behavior is studied by using XCT data collected from fresh and aged catalyst pellets, which exhibit different microstructural features. This study employs time-splitting algorithms to connect disparate reaction and flow timescales, enabling the simulations to achieve realistic deactivation timescales on the order of minutes while the flow time-scales for small particles (100 microns) are several milliseconds. We also introduce a flexible automated python script that writes the necessary files to construct a Mesoflow simulation from user-created chemical mechanisms. We will also introduce a few new features that are added to Mesoflow such as higher order schemes, implicit chemistry integrators and the ability to run on AMD and NVIDIA graphics-processing-units.

AMReX↗

Multiscale investigation of the microstructural mechanisms driving ratchet growth in PBX 9502

The high explosive PBX 9502 undergoes irreversible expansion during thermal cycling (“ratchet growth”). Recent innovations in thermomechanical modeling via homogenization strategies are beginning to incorporate mesoscale information such as grain size, total porosity, and spatial distribution of voids and cracks. To generate a complete experimental data set to challenge and inform these models, PBX 9502 pellets were thermally cycled, cross-sectioned using ion polishing, and imaged in high resolution with scanning electron microscopy. Ratchet growth was found to drive expansion through microcracking. Microcracks were affected by agglomeration of crystals within the PBX. Virgin material showed greater ratchet growth than recycled material.

36 MATERIALS SCIENCE↗

Development of a Model for Irradiation-Induced Grain Growth in UO2 thin films

In this work, we develop a model of irradiation-induced grain growth in UO2 using the MARMOT mesoscale nuclear materials simulation tool. We couple the existing thermally activated grain growth model with a heat conduction model that includes random heat sources representing thermal spikes. We compare the results with the irradiation data on UO2 thin films.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Use of Remote Sensing and In-Situ Observations to Develop and Evaluate Improved Representations of Convection and Clouds for the ACME Model

The overachieving goal of the whole CMDV-MCS project is to improve understanding of warm season continental convection and to develop treatments of convection and microphysics capable of representing mesoscale convective systems (MCSs) features in large-scale models. Our tasks for this project contributing to the overachieving goal include: (1) Improve the ice nucleation formulation for MG2 and P3 cloud microphysics schemes; (2) Improve the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM; and (3) Test the performance of improved ice microphysics in E3SM with observation data. In this project, we have (1) Improved the ice nucleation parameterization for MG2 and P3 in E3SM by implementing two advanced empirical parameterizations with connection to aerosols. The two deterministic heterogeneous ice nucleation parameterizations (i.e., DeMott et al., 2015; Niemand et al., 2012) were merged with the MG2 and P3 cloud microphysics schemes in E3SM. Long-term simulations were conducted to examine the impacts of these new parameterizations on simulated cloud properties; (2) Improved the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM. We evaluated the double Gaussian PDF of vertical velocity simulated by the Cloud Layers Unified By Binormals (CLUBB) and the sub-column vertical velocity sampled from the Subgrid Importance Latin Hypercube Sampler (SILHS) in E3SM. We introduced the vertical velocity variance induced by topographic gravity waves for ice nucleation and droplet activation; and (3) Tested the performance of improved ice microphysics in E3SM with observation data. We tested the new treatments of ice nucleation in the single column model (SCM) mode for the stratiform mixed-phase clouds observed during 9-10 October 2004 in the DOE ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) and for the convective clouds observed on 20 May 2011 in the Midlatitude Continental Convective Clouds Experiment (MC3E). Modeled ice nucleating particles (INPs) concentrations were compared against observations collected around the globe.

54 ENVIRONMENTAL SCIENCES↗

Bridging length scales in organic mixed ionic–electronic conductors through internal strain and mesoscale dynamics

Understanding the structural and dynamic properties of disordered systems at the mesoscale is crucial. This is particularly important in organic mixed ionic-electronic conductors (OMIECs), which undergo significant and complex structural changes when operated in an electrolyte. Here, in this study, we investigate the mesoscale strain, reversibility and dynamics of a model OMIEC material under external electrochemical potential using operando X-ray photon correlation spectroscopy. Our results reveal that strain and structural hysteresis depend on the sample's cycling history, establishing a comprehensive kinetic sequence bridging the macroscopic and microscopic behaviours of OMIECs. Furthermore, we uncover the equilibrium and non-equilibrium dynamics of charge carriers and material-doping states, highlighting the unexpected coupling between charge carrier dynamics and mesoscale order. These findings advance our understanding of the structure-dynamics-function relationships in OMIECs, opening pathways for designing and engineering materials with improved performance and functionality in non-equilibrium states during device operation. Understanding mesoscale structure and dynamics in organic mixed ionic-electronic conductors is crucial. Mesoscale strain kinetics and structural hysteresis have been studied, and they uncover the coupling between charge carrier dynamics and mesoscale order in organic mixed ionic-electronic conductors.

36 MATERIALS SCIENCE↗

Contrasting characteristics of open- and closed-cellular stratocumulus cloud in the eastern North Atlantic

Abstract. Extensive regions of marine boundary layer cloud impact the radiative balance through their significant shortwave albedo while having little impact on outgoing longwave radiation. Despite this importance, these cloud systems remain poorly represented in large-scale models due to difficulty in representing the processes that drive their life cycle and coverage. In particular, the mesoscale organization and cellular structure of marine boundary clouds have important implications for the subsequent cloud feedbacks. In this study, we use long-term (2013–2018) observations from the Atmospheric Radiation Measurement (ARM) Facility's Eastern North Atlantic (ENA) site on Graciosa Island, Azores, Portugal, to identify cloud cases with open- or closed-cellular organization. More than 500 h of each organization type are identified. The ARM observations are combined with reanalysis and satellite products to quantify the cloud, precipitation, aerosol, thermodynamic, and large-scale synoptic characteristics associated with these cloud types. Our analysis shows that both cloud organization populations occur during similar sea surface temperature conditions, but the open-cell cases are distinguished by stronger cold-air advection and large-scale subsidence compared to the closed-cell cases, consistent with their formation during cold-air outbreaks. We also find that the open-cell cases were associated with deeper boundary layers, stronger low-level winds, and higher rain rates compared to their closed-cell counterparts. Finally, raindrops with diameters larger than 1 mm were routinely recorded at the surface during both populations, with a higher number of large drops during the open-cellular cases. The similarities and differences noted herein provide important insights into the environmental and cloud characteristics during varying marine boundary layer cloud mesoscale organization and will be useful for the evaluation of model simulations for ENA marine clouds.

54 ENVIRONMENTAL SCIENCES↗