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 361 records · Page 20

High-Pressure Equation of State of 1,3,5-triamino-2,4,6-trinitrobenzene: Insights into the Monoclinic Phase Transition, Hydrogen Bonding, and Anharmonicity

The high-pressure equation of state (EOS) of energetic materials (EMs) is important for continuum and mesoscale models of detonation performance and initiation safety. Additionally, obtaining a high-fidelity EOS of the insensitive EM 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) has proven to be difficult because of challenges in experimental characterization at high pressures (HPs). In this work, powder X-ray diffraction patterns were fitted using the recently discovered monoclinic I2/a phase above 4 GPa, which shows that TATB is less compressible than when indexed with the triclinic $P\bar{1}$ phase. First-principles calculations were performed with Perdew–Burke–Ernzerhof (PBE) and PBE0 functionals including thermal effects using the $P\bar{1}$ phase. PBE0 improves the description of hydrogen bonding and thus predicts accurate planar a and b lattice parameters under ambient conditions. However, discrepancies in the predicted lattice parameters above 4–10 GPa compared with experimental measurements indexed with $P\bar{1}$ are further evidence of a structural modification at high pressure. Layer sliding defects are formed during molecular dynamics simulations, which induces an anharmonic effect on the thermal expansion of the c lattice parameter. In short, the results provide several insights into determining high-fidelity EOS parameters for TATB and other molecular crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synchrotron Imaging of Pore Formation in Li Metal Solid-State Batteries Aided by Machine Learning

High-rate capable, reversible lithium metal anodes are necessary for next generation energy storage systems. In situ tomography of Li/LLZO/Li cells is carried out to track morphological transformations in Li metal electrodes. Machine learning enables tracking the lithium metal morphology during galvanostatic cycling. Nonuniform lithium electrode kinetics are observed at both electrodes during cycling. Hot spots in lithium metal are correlated with microstructural anisotropy in LLZO. Mesoscale modeling reveals that regions with lower effective properties (transport and mechanical) are nuclei for failure. Advanced visualization combined with electrochemistry represents an important pathway toward resolving non-equilibrium effects that limit rate capabilities of solid-state batteries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Understanding the Influence of Li 7 La 3 Zr 2 O 12 Nanofibers on Critical Current Density and Coulombic Efficiency in Composite Polymer Electrolytes

Composite polymer electrolytes (CPEs)are attractive materials for solid-state lithium metal batteries, owing to their high ionic conductivity from ceramic ionic conductors and flexibility from polymer components. As with all lithium metal batteries, however, CPEs face the challenge of dendrite formation and propagation. Not only does this lower the critical current density (CCD) before cell shorting, but the uncontrolled growth of lithium deposits may limit Coulombic efficiency (CE) by creating dead lithium. Here, we present a fundamental study on how the ceramic components of CPEs influence these characteristics. CPE membranes based on poly-(ethylene oxide) and lithium bis-(trifluoro methane sulfonyl)-imide(PEO-LiTFSI) with Li 7 La 3 Zr 2 O 12 (LLZO) nanofibers were fabricated with industrially relevant roll-to-roll manufacturing techniques. Galvanostatic cycling with lithium symmetric cells shows that the CCD can be tripled by including 50 wt % LLZO, but half-cell cycling reveals that this comes at the cost of CE. Varying the LLZO loading shows that even a small amount of LLZO drastically lowers the CE, from 88% at 0 wt % LLZO to 77% at just 2 wt % LLZO. Mesoscale modeling reveals that the increase in CCD cannot be explained by an increase in the macroscopic or microscopic stiffness of the electrolyte; only the microstructure of the LLZO nanofibers in the PEO-LiTFSI matrix slows dendrite growth by presenting physical barriers that the dendrites must push or grow around. This tortuous lithium growth mechanism around the LLZO is corroborated with mass spectrometry imaging. In conclusion, this work highlights important elements to consider in the design of CPEs for high-efficiency lithium metal batteries.

36 MATERIALS SCIENCE↗

Chemical and Structural Insights into Solid Electrolyte Interphase Evolution for Sodium Metal Electrodes

The solid electrolyte interphase (SEI) critically governs the reversibility of sodium metal batteries, through dynamically mediating ion transport and interfacial reactions. However, its kinetic evolution under operating conditions, and how it influences interfacial stability, remains poorly understood. Here, in this study, we reveal that the SEI undergoes coupled chemical and mechanical changes during sodium plating and stripping, leading to spatial and temporal heterogeneity that drives interfacial degradation. Synchrotron operando grazing-incidence wide-angle X-ray scattering and soft X-ray absorption spectroscopy capture the sequential formation and dissolution of inorganic SEI phases (NaF, NaH, NaOH, Na 2 PO 3 F), accompanied by depth-dependent alterations in organic SEI components. Mesoscale modeling connects this evolving SEI heterogeneity to localized current density fluctuations and stress accumulation at the Na interface, identifying pathways to electrically isolated sodium formation. These findings show that SEI instability fundamentally limits reversibility in sodium metal batteries, and that controlling SEI chemistry–mechanics coupling is essential to achieving its durability.

36 MATERIALS SCIENCE↗

Substitution-Mediated Calcination of Nickel-Based Cathodes: Decoupling Lithiation and Crystallization

Nickel-based layered cathodes such as LiNiO 2 offer high energy density for lithium-ion batteries, yet improvements in cycling performance and safety are required for practical use–often achieved through manganese and cobalt substitution as in LiNi 0.80 Mn 0.10 Co 0.10 O 2 (NMC811). However, how such substitution impacts calcination, the key process that governs lithiation, structural ordering, crystallization, and ultimately the resulting material properties, remains unclear. Here, in this study, we investigate substitution-mediated calcination dynamics in NMC811 compared to LiNiO 2 using multiscale-correlated in situ spectroscopy and atomistic-to-mesoscale modeling. While both systems progress through the same sequence of intermediates toward the thermodynamically favored layered phase, NMC811 exhibits an earlier onset of layering, concurrent with hydroxide decomposition followed by sluggish crystallization. Modeling reveals that Mn and Co lower the energy barrier for lithium incorporation and ordering but increase the penalty for interlayer gliding, thereby slowing crystal growth at elevated temperatures. This substitution-mediated decoupling of lithiation and crystallization explains the fine-grained microstructure observed in NMC811 versus coarsened particles in LiNiO 2 and establishes a mechanistic framework for predictive microstructure engineering of Ni-based cathodes.

Calcination↗

Simulating Mixed‐Phase Open Cellular Clouds Observed During COMBLE: Evaluation of Parameterized Turbulence Closure

Marine cold-air outbreaks, or CAOs, are airmass transformations whereby relatively cold boundary layer (BL) air is transported over relatively warm water. To more deeply understand BL and mixed-phase cloud properties during CAO conditions, the Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) took place from late 2019 into early 2020. During COMBLE, the U.S. Department of Energy's first Atmospheric Radiation Measurement Mobile Facility (AMF1) was deployed to Andenes, Norway, far downstream (~1,000 km) from the Arctic pack ice. This study examines the two most intense CAOs sampled at the AMF1 site. The observed BL structures are open cellular with high (~3–5 km) and cold (–30 to –50 °C) cloud tops, and they often have pockets of high liquid water paths (LWPs; up to ~1,000 g m –2 ) associated with strong updrafts and enhanced turbulence. We use a high-resolution mesoscale model to explore how well four turbulence closure methods represent open cellular clouds. After applying a radar simulator to model outputs for direct evaluation, cloud top properties agree well with AMF1 observations (within ~10%), but radar reflectivity and LWP agreement is more variable. Results suggest that the turbulent Prandtl number may play an important role for the simulated BL and cloud properties. All simulations produce enhanced precipitation rates that are well-correlated with a cloud transition. Finally, the eddy-diffusivity/mass-flux approach produces the deepest cloud layer and therefore the largest and most coherent cellular structures. Furthermore, we recommend the use of a non-local turbulence closure approach to better capture turbulent processes in intense CAOs.

54 ENVIRONMENTAL SCIENCES↗

Linking electronic structure calculations to generalized stacking fault energies in multicomponent alloys

Abstract The generalized stacking fault energy is a key ingredient to mesoscale models of dislocations. Here we develop an approach to quantify the dependence of generalized stacking fault energies on the degree of chemical disorder in multicomponent alloys. We introduce the notion of a “configurationally-resolved planar fault” (CRPF) energy and extend the cluster expansion method from alloy theory to express the CRPF as a function of chemical occupation variables of sites surrounding the fault. We apply the approach to explore the composition and temperature dependence of the unstable stacking fault energy (USF) in binary Mo–Nb alloys. First-principles calculations are used to parameterize a formation energy and CRPF cluster expansion. Monte Carlo simulations show that the distribution of USF energies is significantly affected by chemical composition and temperature. The formalism is broadly applicable to arbitrary crystal structures and alloy chemistries and will enable the development of rigorous models for deformation mechanisms in high-entropy alloys.

Natarajan, Anirudh Raju (ORCID:0000000202529480)↗

Conformal hetero-electrolyte interface between soft oxyhalides and garnet enables low-pressure lithium-reservoir-free solid-state batteries

The operation of solid-state batteries with a lithium metal anode and a high voltage cathode requires solid electrolytes (SEs) that are chemically stable with lithium, have a wide electrochemical window, and accommodate volume changes in the electrodes. Unfortunately, no SE has exhibited satisfactory mechanical and electrochemical properties that fit these requirements to date. Dual solid electrolyte systems that use a different SE for the anolyte and catholyte present a viable solution. Here, we focus on oxyhalide SEs that demonstrate superior ionic conductivity and cathode compatibility, where their lithium metal reactivity and poor reduction stability can be resolved using a lithium garnet (Li 6.5 La 3 Zr 1.5 Ta 0.5 O 12 , LLZTO) separator. Nonetheless, this imposes a new hetero-electrolyte (H-E) interface at the anolyte|catholyte contact that defines the ion transport across the boundary. We report its promising properties, which are deconvoluted from the electrical measurements of bilayer symmetric cells, for three representative oxyhalide catholytes, LiNbOCl 4 , LiTaOCl 4 , and Li 3 Al 3 O 2 Cl 8 . Pressure-dependent measurements reveal that the relative softness of the oxyhalides (hardness ≤0.4 GPa) enables H-E resistances lower than 150 Ω cm 2 at 2–3 MPa. Mesoscale modelling reveals that the transfer-active contact area of oxyhalides with LLZTO is about 2–3-fold higher than that of argyrodite, Li 6 PS 5 Cl. The low H-E resistance of the LiNbOCl 4 |Li 6.5 La 3 Zr 1.5 Ta 0.5 O 12 dual electrolyte enables the cycling of a Li|LiNi 0.82 Mn 0.07 Co 0.11 O 2 full cell with a high discharge capacity (200 mA h g −1 ) at 60 °C and ∼7 MPa. Importantly, we demonstrate a Li-reservoir-free full cell with high Coulombic efficiency (>99.5%) and capacity at 1 MPa using this approach coupled with a garnet-silver interlayer.

Palmer, Max [University of California, Santa Barba↗

LiDAR Measurements of Wind Shear Exponents and Turbulence Intensity Offshore the Northeast United States

This paper presents wind speed shear exponents and turbulence intensity measurements collected from LiDARs measuring wind speeds from 40m to 200m above sea-level and provides comparisons to industry design guidelines. The high-altitude wind speed data are unique and represent some of the first measurements made offshore in this part of the country, which is actively being developed for offshore wind. The data is used to support the New England Aqua Ventus I Floating Offshore Wind Farm to be located 17km offshore the Northeast United States. Multiple LiDAR measurements were made using a DeepCLiDAR floating buoy and LiDARs located on a nearby island. Here, the measured wind speed shear exponents are compared against industry standard mesoscale model outputs and offshore design codes including the American Bureau of Shipping, American Petroleum Institute, and DNV-GL guides. Significant variation in the vertical wind speed profile occurs throughout the year which is not addressed in design standards. Additionally, turbulence intensity measurements made from the LiDAR, although not widely accepted in the scientific community, are presented and compared against industry guidelines.

42 ENGINEERING↗

Phase Field Dislocation Dynamics (PFDD) version 2.x

This disclosure is for version 2.x of a mesoscale model called Phase Field Dislocation Dynamics (PFDD). PFDD is used for investigating deformation in nanoscale (grain sizes of ~300 nm and less) materials, such as metals and alloys. This approach models the motion and interaction of individual defects, namely dislocations, in the material using scalar-valued phase field variables, also called order parameters. The system is evolved through energy minimization thus the model calculates the total energy density in terms of the phase field variables. The energy minimization is completed using the Ginzburg-Landau equation, and is implemented with explicit time integration. The total system energy can be comprised of several terms, including the strain energy (which describes dislocation-dislocation interactions), the energy due to an applied stress (dislocation interactions with the applied stress), and a core/lattice (perfect dislocations) or generalized stacking fault (partial dislocations) energy (described the dislocation core structure). The latter term in particular may vary based on the crystal structure being modeled and is typically informed using lower length scale (e.g., atomistic) approaches, although no such (atomistic) calculations are completed within the PFDD algorithm. This basic formulation was previously reviewed by Los Alamos National Laboratory and released under license number C17113. This previously reviewed version we will henceforth refer to as PFDD v1.0. PFDD v1.0 consisted of 2 codes (one parallel and one serial) plus input files, all written in the C language. This new disclosure is addressing the next versions of the PFDD, versions 2.x. There have been several enhancements of PFDD v1.0, which are described here and included in the attached code, which we will refer to as PFDD v2.0. There are also several new features described here that are either planned or already in process and are expected to be subsequent releases, i.e., v2.1, v2.2, ...v2.x.

Hunter, Abigail↗

A Multiscale Simulation Study of the Effects of Nantucket Island on Atmospheric Boundary Layer Flow

Multiscale atmospheric simulations were conducted to investigate how Nantucket Island modifies the marine atmospheric boundary layer (ABL) flow, how the mechanisms responsible are represented in mesoscale versus large-eddy simulation (LES) domains, and potential impacts on wind energy production downstream, in moderate and high wind speed cases. Here, to isolate the effects of the island, each case study was simulated (i) with and without island topography included and (ii) with observed versus spatially uniform sea surface temperature (SST) values, the latter removing the influence of spatial variability in SST, while also enabling the examination of how the island effects vary with near-surface thermodynamic stability over the surrounding ocean. Island effects were found to be strongest during daytime heating, which generates convective rolls that redistribute momentum vertically over the island and downstream. Colder SSTs and weaker winds strengthen these rolls and extend their influence. Flow statistics within a representative wind turbine rotor-swept area downstream show that while vertical momentum redistribution enhances near-surface wind speeds, those aloft decrease, resulting in minimal net change to rotor-equivalent wind speed and slightly reduced vertical shear. Comparisons of mesoscale and LES domains reveal that both capture the roughness-driven near-surface slowdown, but mesoscale models underestimate vertical momentum transport compared to LES, which explicitly resolves convective rolls. Increasing LES resolution beyond the grid spacing required to begin capturing the largest convective roll structures ∼ O(100) m yields minimal additional momentum transport but increases resolved turbulence kinetic energy, relevant for stress loading.

Atmosphere-land interaction↗

Numerical Simulations of Orographic Convection across Multiple Gray Zones

Idealized simulations are used to determine the sensitivity of moist orographic convection to horizontal grid spacing Δ h . In simulated mechanically (MECH) and thermally (THERM) forced convection over an isolated ridge, Δ h is varied systematically over both the deep-convection (Δ h ~ 10 –1 km) and turbulence (Δ h ~ 1 km–100 m) gray zones. To aid physical interpretation, a new parcel-based bulk entrainment/detrainment diagnosis for horizontally heterogeneous flows is developed. Within the deep-convection gray zone, the Δ h sensitivity is dominated by differences in parameterized versus explicit convection; the former initiates convection too far upstream of the ridge (MECH) and too early in the diurnal heating cycle (THERM). These errors stem in part from a large underprediction of parameterized entrainment and detrainment. Within the turbulence gray zone, sensitivities to Δh arise from the representation of both subcloud- and cloud-layer turbulence. As Δ h is decreased, MECH exhibits stronger cloud-layer entrainment to enhance the convective mass flux M co , while THERM exhibits stronger detrainment to suppress Mco and delay convection initiation. The latter is reinforced by increased subcloud turbulence at smaller Δ h , which leads to drying and diffusion of the central updraft responsible for initiating moist convection. Numerical convergence to a robust solution occurs only in THERM, which develops a fully turbulent flow with a resolved inertial subrange (for Δ h ≤ 250 m). In MECH, by contrast, turbulent transition occurs within the orographic cloud, the details of which depend on both physical location and Δ h .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Resolving intragranular stress fields in plastically deformed titanium using point-focused high-energy diffraction microscopy

Abstract The response of a polycrystalline material to a mechanical load depends not only on the response of each individual grain, but also on the interaction with its neighbors. These interactions lead to local, intragranular stress concentrations that often dictate the initiation of plastic deformation and consequently the macroscopic stress–strain behavior. However, very few experimental studies have quantified intragranular stresses across bulk, three-dimensional volumes. In this work, a synchrotron X-ray diffraction technique called point-focused high-energy diffraction microscopy (pf-HEDM) is used to characterize intragranular deformation across a bulk, plastically deformed, polycrystalline titanium specimen. The results reveal the heterogenous stress distributions within individual grains and across grain boundaries, a stress concentration between a low and high Schmid factor grain pair, and a stress gradient near an extension twinning boundary. This work demonstrates the potential for the future use of pf-HEDM for understanding the local deformation associated with networks of grains and informing mesoscale models. Graphical abstract

36 MATERIALS SCIENCE↗

Transmitted wave measurements in cold sprayed materials under dynamic compression

Spray-formed materials have complex microstructures which pose challenges for microscale and mesoscale modeling. To constrain these models, experimental measurements of wave profiles when subjecting the material to dynamic compression are necessary. The use of a gas gun to launch a shock into a material is a traditional method to understand wave propagation and provide information of time-dependent stress variations due to complex microstructures. This data contains information on wave reverberations within a material and provides a boundary condition for simulation. Here we present measurements of the wavespeed and wave profile at the rear surface of tantalum, niobium, and a tantalum/niobium blend subjected to plate impact. Measurements of the Hugoniot elastic limit are compared to previous work and wavespeeds are compared to longitudinal sound velocity measurements to examine wave damping due to the porous microstructure.

36 MATERIALS SCIENCE↗

Final report

The bulk of the research results supported by this grant had to do with the diurnal cycle of convection and/or convective organization over the Amazon with some analysis also done over the Southern Great Plains, and large-scale drivers of changes in these convective properties like the Amazonian low-level jet. A comprehensive suite of DOE observational and modeling tools was utilized. Overall, the research results collectively emphasize the significance of accurately observing and modeling mesoscale convective systems and the diurnal cycle of convection, particularly in the Amazon, to enhance understanding and forecasting of regional precipitation patterns and climate dynamics.

58 GEOSCIENCES↗

Development of Accelerated Kinetic Monte Carlo Code for Simulation of Helium Bubble Evolution

A mesoscale model to predict helium bubble evolution is needed for tritium applications. Such a model requires that the conventional kinetic Monte Carlo (kMC) simulations be significantly accelerated. The objective of this report is to (a) highlight the concepts and mathematical expressions of the accelerated method for defect implementation that have not been published, (b) show an example input file to run the kMC code, and (c) provide suggestions on future improvement following my retirement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Provider - University of Washington - Raw Data

GOES satellite, WSR-88D radar, sounding balloon, surface weather observation and mesoscale model data from the area within, up, and downstream of the study area for use in event logging procedures and case study investigations.

17 WIND ENERGY↗

Provider - University of Washington - Reviewed Data

GOES satellite, WSR-88D radar, sounding balloon, surface weather observation and mesoscale model data from the area within, up, and downstream of the study area for use in event logging procedures and case study investigations.

17 WIND ENERGY↗