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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Performance of reanalysis and mesoscale models off the coast of Hawai'i

The eastern Hawai'i coast in the United States is characterized by considerable wind resource fuelled by persistent trade winds, making it an important area for energy research. The need is strong for reanalyses and higher-resolution regional simulations where observations have been historically limited, such as Hawai'i's offshore environments. However, studies using offshore observations in other parts of the world have shown that significant errors can occur in reanalyses and wind datasets, which can lead to inaccurate estimates of wind energy generation, payback periods, and extreme weather risks at project locations. The degree of such errors is influenced by a number of factors, including spatial resolution and the handling of processes within the planetary boundary layer (PBL). In this work, we provide a wind resource characterization from year-long lidar buoy measurements off the eastern coast of O'ahu, Hawai'i, an environment previously unobserved at the rotor level, and use the characterization to evaluate the performance of two simulation datasets. The O'ahu deployment location is meteorologically unique and less complex than land-based wind resource characterizations, being strongly characterized by trade winds with minimal land–atmosphere interaction influences. Despite the unique and fairly consistent meteorological conditions, we hypothesize that distinct simulation datasets will exhibit diverse ranges of errors similar to those that have been seen for other offshore locations. We find the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis version 5 (ERA5) to strongly underestimate observed wind speeds at the O'ahu location (bias = −1.54 m s −1 at a height of 140 m above sea level), while a regional Weather Research and Forecasting Model (WRF) simulation produced by the University of Hawai'i (UH-WRF) provides a significantly smaller wind speed bias (−0.25 m s −1 ), highlighting the value of running regional, higher-resolution simulations. The large bias noted for ERA5 is driven by significant underestimation of fast wind speeds (>9 m s −1 ), which the study site is largely characterized by, along with discontinuities in the ERA5 diurnal cycle. We also speculate that the relative sparsity of observations for data assimilation in this remote part of the world could influence the performance of ERA5 and that challenges with characterizing island effects could impact the performance of both datasets.

17 WIND ENERGY↗
DOE OSTI2023

Understanding grain boundary segregation under irradiation: A mesoscale modeling study

Grain boundaries (GBs) are an important component in the design of irradiation resistant materials as they play a crucial role as point defect sinks for non-equilibrium vacancies and self-interstitials. However, in the presence of sustained point defect fluxes, the coupled transport of atomic components with inherently different diffusivities often leads to radiation-induced segregation (RIS), with undesired effects on properties such as intergranular corrosion and stress-assisted cracking. In this talk, we investigate the fundamental mechanisms leading to GB segregation in BCC FeCr-based alloys using multiscale modeling and experiments. A phase-field model is employed to account for both thermodynamic and kinetic segregation mechanisms. The full Onsager transport matrix is derived using kinetic Monte Carlo, molecular dynamics, and kinetic cluster expansion (KineCluE) calculations. The thermodynamics of GB segregation is described via an atomic density-based modification to the CALPHAD free energy. The modeling predictions are then compared against segregation data from atom probe characterization for various temperatures, compositions, and neutron irradiation conditions. We will conclude by discussing the utility of the multiscale modeling approach to shed light on the fundamental mechanisms leading to solute-GB interaction.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development and preliminary validation of a mechanistic multiscale model for fuel-cladding chemical interaction in metallic nuclear fuels

Despite decades of fuel rod material and design improvements, fuel-cladding chemical interaction (FCCI) remains the single-most lifetime-limiting behavior for modern metallic fuel rods. Constraining fuel lifetime increases operating costs, limiting the economic viability of commercializing metallic nuclear fuel technology. A mechanistic multiscale model utilizing the finite element method-based MARMOT and BISON codes was developed to more confidently predict cladding-side FCCI and its impact on fuel performance. The new BISON model incorporates mesoscale models for the effects of fuel microstructure evolution on the transport of wastage-inducing lanthanides through the fuel and for the kinetics of cladding wastage layer growth. The mesoscale models, in turn, build on lanthanide transport property data obtained from the atomistic scale. Preliminary validation studies using wastage thickness and cladding profilometry data from four fuel rods irradiated in Experimental Breeder Reactor II experiment X447 and one fuel rod from Fast Flux Test Facility experiment IFR1 show that the new model predicts cladding wastage and its effects on cladding deformation as well as existing empirical FCCI correlations. The new model is expected to aid in the design of new metallic fuel concepts, including fuel additives, cladding liners, and sodium-free annular fuel geometries. In conclusion, future work will focus on broader validation and refinement of the model’s treatment of different fuel alloys and cladding materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Tropical Cyclone Super Resolution using conditional diffusion denoising probabilistic model from mesoscale simulation to LES

Accurate modeling of tropical cyclone wind fields is essential for the design, risk assessment, and operational planning of offshore energy infrastructure. While mesoscale simulations are widely used thanks to their computational efficiency, they lack the necessary resolution to capture key features such as wind shear and veer profiles as well as the distribution turbulent kinetic energy (TKE). High-fidelity large-eddy simulation (LES) models on the other hand, can resolve turbulent structures and provide a more accurate representation of the complex wind field, albeit at a higher computational cost. To address this modeling gap, we introduce a two-part generative framework to enhance the resolution and physics-capturing ability of mesoscale simulations. First, a reduced-order model based on Karhunen–Loève (KL) decomposition is used to extract dominant spatial modes from one-dimensional mean wind profiles. A multilayer perceptron (MLP) is trained to map mesoscale mode weights to their LES counterparts, enabling accurate reconstruction of vertical velocity profiles. Second, a conditional Diffusion Denoising Probabilistic Model (DDPM) is developed to super-resolve coarse and low-fidelity mesoscale velocity fields, recovering fine-scale turbulence structures and stress distributions. The framework is evaluated across different tropical cyclone intensity categories defined by the Saffir–Simpson scale and demonstrates strong performance in both interpolation and extrapolation tasks. The generated fields accurately reproduce spatial coherence, stress distributions, and spectral energy characteristics observed in LES data. By bridging the fidelity gap between mesoscale and LES outputs, this approach offers a scalable, data-driven solution for enhancing the representation of tropical cyclone wind fields, enabling more robust offshore energy infrastructure systems design in tropical-cyclone-prone areas.

17 WIND ENERGY↗
DOE OSTI2024

An atomistic survey of shear coupling in asymmetric tilt grain boundaries and interpretation using the disconnections framework

Grain Boundaries (GB) play an important role in determining the behavior of polycrystalline materials. While the mechanisms of motion and associated shear response for symmetric tilt grain boundaries (STGBs) are well studied, the same is not true for asymmetric tilt grain boundaries (ATGBs) despite their greater prevalence in polycrystals. Here, this study aims to investigate the shear response of a large collection of asymmetric tilt grain boundaries (ATGBs) using molecular dynamics (MD) simulations and interpret the data using a discrete disconnections model that works remarkably well for STGBs. MD simulations of shear-driven ATGBs show that the plastic shear (shear coupling factor) in the region swept by a GB exhibits a complex dependence on the inclination angle, and this dependence changes with the misorientation of the GB. In addition, the shear response was observed to be highly sensitive to the applied shear rate and temperature. Recognizing the spatial and temporal scale limitations of MD simulations, we extended the discrete disconnections mesoscale model of Khateret al. (2012) to calculate the nucleation barriers of disconnection modes and predict the effective shear coupling of an ATGB. We observed that the mesoscale model’s predictions of the shear coupling factor of ATGBs do not agree with those observed in MD simulations. Finally, we examine the hypotheses of our mesoscale model that contribute to disagreements between MD simulations and the mesoscale model and propose improvements to the mesoscale model for future work.

36 MATERIALS SCIENCE↗
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 OSTI2020

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets—Part I: Mesoscale simulations

Recently a mesoscale model was developed that includes thermomechanically consistent constitutive equations with dissipation due to porous compaction, inelastic distortional deformation rate, and slipping on persistent dry joint surfaces. In this Part I of a two-part paper, we describe enhancements of the mesoscale model which include: normal compliance with hysteresis, non-persistent joints with finite area, joint sets with angular deviation, and the weakening effect of fluid pressure in the joints. Examples study the influences of these effects on average strength. These computationally expensive simulations are used to inspire functional forms in a continuum model, described in Part II, which is needed for large scale simulations of real applications.

42 ENGINEERING↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Inert mesoscale simulations

Mesoscale modeling of shock waves in Ni+Al multilayers poses significant challenges that are due, in part, to shock-induced chemical reactions. Current modeling approaches utilize reactive molecular dynamics (MD), but they are limited to resolving domains of only a few hundred nanometers. In contrast, actual multilayer superlattices can be tens of micrometers thick, and they exhibit non-ideal (i.e., wavy) interfaces. The second part of our research builds upon previous work developing physically based, thermodynamically complete equations of state for various Ni and Al intermetallic compositions. Here, we introduce a novel workflow for high-fidelity mesoscale simulations of Ni+Al multilayers using a continuum hydrocode. By increasing the simulation domain size beyond MD limitations (e.g., 2 × 6 μm 2 ) and incorporating explicit interfacial roughness, we investigate the shock response of Ni+Al multilayers at previously unexplored scales. Our experimental design encompasses nine multilayer geometries with varying roughness amplitudes and tilt angles (θ = 15°, 30°, and 45°), alongside 19 flyer impact velocities ranging from 0.3 to 3.0 km/s, resulting in a total of 171 high-fidelity simulations. The bulk shock state from inert 2D mesoscale simulations aligns with the law of mixtures, while temperature and pressure fluctuations strongly correlate with multilayer geometry types. A new metric dubbed the “hot spot probability integral” shows a greater dependence on a tilt angle than interfacial roughness.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2024

Lithium Plating Characteristics in High Areal Capacity Li-Ion Battery Electrodes

Li-ion battery degradation and safety events are often attributed to undesirable metallic lithium plating. Since their release, Li-ion battery electrodes have been made progressively thicker to provide a higher energy density. However, the propensity for plating in these thicker pairings is not well understood. Herein, we combine an experimental plating-prone condition with robust mesoscale modeling to examine electrode pairings with capacities ranging from 2.5 to 6 mAh/cm 2 and negative to positive (N/P) electrode areal capacity ratio from 0.9 to 1.8 without the need for extensive aging tests. Using both experimentation and a mesoscale model, we identify a shift from conventional high state-of-charge (SOC) type plating to high overpotential (OP) type plating as electrode thickness increases. Further, these two plating modes have distinct morphologies, identified by optical microscopy and electrochemical signatures. We demonstrate that under operating conditions where these plating modes converge, a high propensity of plating exists, revealing the importance of predicting and avoiding this overlap for a given electrode pairing. Further, we identify that thicker electrodes, beyond a capacity of 3 mAh/cm 2 or thickness >75 μm, are prone to high OP, limiting negative electrode (NE) utilization and preventing cross-sectional oversizing the NE from mitigating plating. Here, it simply contributes to added mass and volume. The experimental thermal gradient and mesoscale model either combined or independently provide techniques capable of probing performance and safety implications of mild changes to electrode design features.

25 ENERGY STORAGE↗
DOE OSTI2022

Mesoscale simulations of pressure-shear loading of granular tungsten carbide

Numerical simulations of pressure-shear loading of a granular material are performed using the shock physics code CTH. A simple mesoscale model for the granular material is used that consists of a randomly packed arrangement of solid circular or spherical grains of uniform size separated by vacuum. The grain material is described by a simple shock equation of state, elastic perfectly plastic strength model, and fracture model with baseline parameters for WC taken from previous mesoscale modeling work. Simulations using the baseline material parameters are performed at the same initial conditions of pressure-shear experiments on dry WC powders. Except for some localized flow regions appearing in simulations with an approximate treatment of sliding interfaces among grains, the samples respond elastically during shear, which is in contrast to experimental observations. By extending the simulations to higher shear wave amplitudes, macroscopic shear failure of the simulated samples is observed with the shear strength increasing with increasing stress confinement. The shear strength is also found to be strongly dependent on the grain interface treatment and on the fracture stress of the grains, though the variation in shear strength due to fracture stress decreases with increasing stress confinement. At partial compactions, the transverse velocity histories show strain-hardening behavior followed by formation of a shear interface that extends through the transverse dimensions of the sample. Near full compaction, no strain hardening is observed and, instead, the sample transitions sharply from an elastic response to formation of an internal shear interface. Agreement with experiment is shown to worsen with increasing confinement stress with simulations overpredicting the shear strengths measured in experiment. The source of the disagreement can be ultimately attributed to the Eulerian nature of the simulations, which do not treat contact and fracture realistically.

36 MATERIALS SCIENCE↗
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↗
DOE OSTI2022

The Effect of Surface Heating Heterogeneity on Boundary Layer Height and Its Dependence on Background Wind Speed

The planetary boundary layer height (PBLH) is a fundamental variable of the planetary boundary layer (PBL). The effect of surface heterogeneity is challenging in the PBL parameterization for numerical weather and climate models. Here we use large eddy simulation data to investigate the impact of surface heating heterogeneity on the PBLH spatial variations and mean values over a domain of typical mesoscale model grid size. It is found that this impact depends on background wind condition. Variance decomposition and Fourier spectra both reveal that for the fully developed sheared convective PBL, surface heterogeneity contributes little to the spatial variation of PBLH, except for the case with heterogeneity scale of 14.4 km. This suggests that surface heterogeneity with length scales less than ~10 km, which is typical for mesoscale modeling, does not induce significant subgrid spatial variation of PBLH. The domain-averaged PBLH is found to generally increase as surface heterogeneity scale decreases, until reaching the scale on the order of PBLH. The findings can be used to inform the parameterization development of PBL processes over heterogeneous surface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Bicrystallography-informed Frenkel–Kontorova model for interlayer dislocations in strained 2D heterostructures

In recent years, van der Waals (vdW) heterostructures and homostructures, which consist of stacks of two-dimensional (2D) materials, have risen to prominence due to their association with exotic quantum phenomena originating from correlated electronic states harbored by them. Atomistic scale relaxation effects play an extremely important role in the electronic scale quantum physics of these systems, providing means of manipulation of these materials and allowing them to be tailored for emergent technologies. We investigate such structural relaxation effects in this work using atomistic and mesoscale models, within the context of twisted bilayer graphene — a well-known heterostructure system that features moiré patterns arising from the lattices of the two graphene layers. For small twist angles, atomic relaxation effects in this system are associated with the natural emergence of interface dislocations or strain solitons, which result from the cyclic nature of the generalized stacking fault energy (GSFE), that measures the interface energy based on the relative movement of the two layers. Here, in this work, we first demonstrate using atomistic simulations that atomic reconstruction in bilayer graphene under a large twist also results from interface dislocations, although the Burgers vectors of such dislocations are considerably smaller than those observed in small-twist systems. To reveal the translational invariance of the heterointerface responsible for the formation of such dislocations, we derive the translational symmetry of the GSFE of a 2D heterostructure using the notions of coincident site lattices (CSLs) and displacement shift complete lattices (DSCLs). The workhorse for this exercise is a recently developed Smith normal form bicrystallography framework. Next, we construct a bicrystallography-informed and frame-invariant Frenkel–Kontorova model, which can predict the formation of strain solitons in arbitrary 2D heterostructures, and apply it to study a heterostrained, large-twist bilayer graphene system. Our mesoscale model is found to produce results consistent with atomistic simulations. We anticipate that the model will be invaluable in predicting structural relaxation and for providing insights into various heterostructure systems, especially in cases where the fundamental unit cell is large and therefore, atomistic simulations are computationally expensive.

2D heterostructures↗
DOE OSTI2026

Nonlocal Effect of Percolated Particle Networks on Viscoelasticity of Polymer–Filler Nanocomposites: A Mesoscale Simulation Study

With nanoparticles (NPs) as fillers, polymer nanocomposites (PNCs) usually exhibit enhanced mechanical properties. However, a direct connection between the microscopic structural relaxation and macroscopic mechanical properties of PNCs remains to be established. To investigate the micro-to-macro connection, we develop a mesoscale model, in which the NPbridging polymer chains are represented by a dynamic bonded interaction between NPs, and the bulk polymer matrix is implicitly modeled by overdamped Langevin dynamics. Extensive equilibrium simulations are performed to quantify the microscopic dynamics of model PNCs. Systematic analyses of modified Rouse dynamics, dynamic structure factor, and relaxation modulus uncover that the microscopic relaxation dynamics of PNCs are significantly decelerated across different length scales because of nonlocal effects of percolated particle networks a phenomenon that has not been adequately captured in prior simulation studies. We find that NPvolume- fraction and NP-bonding-energy barrier are the two critical variables that affect bulk viscoelasticity the most. The proposed mesoscale model is versatile and provides a powerful framework for studying structure−property relations of different PNCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗