Search NASA⌕ Search

SEARCH · Search NASA

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 OSTI2024

Interdependence of Support Wettability ‐ Electrodeposition Rate‐ Sodium Metal Anode and SEI Microstructure

Abstract This study examines how current collector support chemistry (sodiophilic intermetallic Na 2 Te vs. sodiophobic baseline Cu) and electrodeposition rate affect microstructure of sodium metal and its solid electrolyte interphase (SEI). Capacity and current (6 mAh cm −2 , 0.5–3 mA cm −2 ) representative of commercially relevant mass loading in anode‐free sodium metal battery (AF‐SMBs) are analyzed. Synchrotron X‐ray nanotomography and grazing‐incidence wide‐angle X‐ray scattering (GIWAXS) are combined with cryogenic ion beam (cryo‐FIB) microscopy. Highlighted are major differences in film morphology, internal porosity, and crystallographic preferred orientation e.g. (110) vs. (100) and (211) with support and deposition rate. Within the SEI, sodium fluoride (NaF) is more prevalent with Te−Cu versus sodium hydride (NaH) and sodium hydroxide (NaOH) with baseline Cu. Due to competitive grain growth the preferred orientation of sodium crystallites depends on film thickness. Mesoscale modeling delineates the role of SEI (ionic conductivity, morphology) on electrodeposit growth and onset of electrochemical instability.

Lo, Chang‐An↗
DOE OSTI2024

Alumina – Stabilized SEI and CEI in Potassium Metal Batteries

Aluminum oxide (Al 2 O 3 ) nanopowder is spin-coated onto both sides of commercial polypropene separator to create artificial solid-electrolyte interphase (SEI) and artificial cathode electrolyte interface (CEI) in potassium metal batteries (KMBs). Here, this significantly enhances the stability, including of KMBs with Prussian Blue (PB) cathodes. For example, symmetric cells are stable after 1,000 cycles at 0.5 mA/cm 2 –0.5 mAh/cm 2 and 3.0 mA/cm 2 –0.5 mAh/cm 2 . Alumina modified separators promote electrolyte wetting and increase ionic conductivity (0.59 vs. 0.2 mS/cm) and transference number (0.81 vs. 0.23). Cryo-stage focused ion beam (cryo-FIB) analysis of cycled modified anode demonstrates dense and planar electrodeposits, versus unmodified baseline consisting of metal filaments (dendrites) interspersed with pores and SEI. Alumina-modified CEI also suppresses elemental Fe crossover and reduces cathode cracking. Mesoscale modeling of metal – SEI interactions captures crucial role of intrinsic heterogeneities, illustrating how artificial SEI affects reaction current distribution, conductivity and morphological stability.

25 ENERGY STORAGE↗
DOE OSTI2024

Interdependence of Support Wettability - Electrodeposition Rate- Sodium Metal Anode and SEI Microstructure

Here, this study examines how current collector support chemistry (sodiophilic intermetallic Na 2 Te vs. sodiophobic baseline Cu) and electrodeposition rate affect microstructure of sodium metal and its solid electrolyte interphase (SEI). Capacity and current (6 mAh cm -2 , 0.5–3 mA cm -2 ) representative of commercially relevant mass loading in anode-free sodium metal battery (AF-SMBs) are analyzed. Synchrotron X-ray nanotomography and grazing-incidence wide-angle X-ray scattering (GIWAXS) are combined with cryogenic ion beam (cryo-FIB) microscopy. Highlighted are major differences in film morphology, internal porosity, and crystallographic preferred orientation e.g. (110) vs. (100) and (211) with support and deposition rate. Within the SEI, sodium fluoride (NaF) is more prevalent with Te-Cu versus sodium hydride (NaH) and sodium hydroxide (NaOH) with baseline Cu. Due to competitive grain growth the preferred orientation of sodium crystallites depends on film thickness. Mesoscale modeling delineates the role of SEI (ionic conductivity, morphology) on electrodeposit growth and onset of electrochemical instability.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Impact of Salt Chemistry and Interlayer Architecture on Solid Electrolyte Interphase Composition and Performance

Reversible lithium metal anodes require interphases that are chemically stable, mechanically robust, and compatible with scalable processing. Here, we use slot-die-coated, salt-rich PEO interlayers to design SEI nanolayers on Cu current collectors and isolate how salt anion chemistry controls interphase composition and lithium plating/stripping behavior in architectures relevant to Li plating and anode-free cells. Mesoscale modeling shows that dendrite initiation at SEI defects is governed by a tradeoff among interphase thickness, ionic conductivity, and stiffness through defect-driven current focusing, motivating controlled formation from thin precursor layers. Using UHV-transferred XPS and ToF-SIMS depth profiling after controlled Li deposition, we find that fluorinated sulfonylimide salts (LiTFSI, LiFSI) form mixed inorganic/organic interphases containing LiF and PEO-derived alkoxides, with LiFSI producing a more LiF-dominant Li-facing surface than LiTFSI under identical conditions. In contrast, LiNO 3 -containing interlayers yield a comparatively thick nitrate-derived interphase that reduces upon Li contact to an N-rich inorganic layer (Li–N/N 3− , LiN x O y ). In Cu| |Li cells, these designed interlayers reduce interfacial resistance and improve Coulombic efficiency (CE) and critical current density (CCD) relative to uncoated controls, with LiFSI- and LiNO 3 -based interlayers providing the highest CE/CCD. Overall, the results demonstrate a manufacturing-relevant approach to engineer SEI nanolayers by salt chemistry, providing pathways to improve reversible lithium metal anodes while reducing reliance on LiTFSI.

anion chemistry↗
DOE OSTI2026

Machine learning insights into microstructural origins of transport and mechanical properties in porous microstructures

Multifunctional porous materials are increasingly needed across various fields, but their complex microstructures create significant challenges due to the intricate microstructure-property relationships. This complexity, combined with limitations of traditional analysis methods, hinders efforts to understand and optimize microstructure–property relationships. Here, to address this, we integrate physics-based mesoscale modeling with interpretable machine learning (ML) to uncover how microstructural features govern effective diffusivity and elastic modulus. At constant porosity, we show diffusivity varies by over 150 × and modulus by ∼50 ×, highlighting the power of microstructure engineering. Statistical analysis reveals bimodal behavior in diffusivity and unimodal in modulus. ML identifies connectivity as the dominant factor, while modulus is also sensitive to domain size and feature interactions. Controlled simulations further highlight domain shape as a critical feature for modulus. This framework enables efficient exploration of microstructure-property correlations, offering new insights to guide the design of advanced porous materials.

Bicontinuous microstructure↗
DOE OSTI2026

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↗
DOE OSTI2026

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↗
DOE OSTI2024

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↗
DOE OSTI2026

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↗
DOE OSTI2026

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↗
DOE OSTI2024

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↗
DOE OSTI2025

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↗
DOE OSTI2024

Realistic Precipitation Diurnal Cycle in Global Convection-Permitting Models by Resolving Mesoscale Convective Systems

Accurately representing the precipitation diurnal cycle has long been a challenge for global climate models (GCMs). Here we evaluate the precipitation diurnal cycle in the DYAMOND global convection-permitting models (CPMs) and CMIP6 HighResMIP models. Comparison of the high- (25–50 km) and low-resolution (100–250 km) models with parameterized convection in HighResMIP shows that simply increasing model resolution does not noticeably improve the precipitation diurnal cycle. In contrast, CPMs can better capture the observed amplitude and timing of precipitation diurnal cycle. However, the simulated spatial variation of timing in CPMs is smaller than observed, leading to an exaggeration of the spatially averaged diurnal amplitude. The better-simulated precipitation diurnal cycle in the CPMs is tied to mesoscale convective systems (MCSs), which contribute about half of the total precipitation. The observed life cycle of MCSs, including initiation and mature stages, is well captured in the CPMs, leading to a more realistic precipitation diurnal cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sensitivities of the West Greenland Current to Greenland Ice Sheet Meltwater in a Mesoscale Ocean/Sea Ice Model

Abstract Meltwater from the Greenland Ice Sheet can alter the continental shelf/slope circulation and cross-shelf freshwater fluxes and limit deep convection in adjacent basins through surface freshening. We explore the impacts on the West Greenland Current and eastern Labrador Sea with different vertical distributions of the meltwater forcing. In this study, we present the results from global coupled ocean/sea ice simulations, forced with atmospheric reanalysis, that are mesoscale eddy-active (∼2–3-km horizontal spacing) and eddy-permitting (∼6–7-km horizontal spacing) in the study region. We compare the West Greenland Current in mesoscale eddy-active and eddy-permitting without meltwater to highlight the role of small-scale features. The mesoscale eddy-active configuration is then used to assess the change in the eastern Labrador Sea when meltwater is added to the surface or vertically distributed to account for mixing within fjords. In both simulations with meltwater, the West Greenland and West Greenland Coastal Currents are faster than in the simulation with no meltwater; their mean surface speeds are the highest in the vertical distribution case. In the latter case, there is enhanced baroclinic conversion at the shelf break compared to the simulation with no meltwater. When meltwater is vertically distributed, there is an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy. In addition, in the eastern Labrador Sea, the salinity is lower and the meltwater volume is greater when meltwater is vertically distributed. Therefore, the West Greenland Current is sensitive to how meltwater is added to the ocean with implications for the freshening of the Labrador Sea. Significance Statement Our goal is to understand how the flux of freshwater across the West Greenland continental slope into the Labrador Sea is modified by meltwater from the Greenland Ice Sheet. We compare the simulations of the ocean that capture key dynamics along the West Greenland continental slope that have no meltwater, meltwater added to the ocean surface, and meltwater distributed vertically to represent the mixing within fjords. When meltwater is added, the currents along the continental slope are faster, with the greatest increase when meltwater is vertically distributed. In that case, there is enhanced freshening of the Labrador Sea because modified density gradients generate more eddies. Proper representation of the vertical structure of meltwater is important for projecting the impact of freshwater on the subpolar North Atlantic.

58 GEOSCIENCES↗
DOE OSTI2025

Multiscale Modeling of the Mechanical Response of Silicon Carbide Composite Within the Accelerated Fuel Qualification Framework

The accelerated fuel qualification (AFQ) framework has been used for the initial development of multiscale modeling of silicon carbide (SiC) fiber reinforced composite (SiC-SiC). The AFQ framework provides a methodology to leverage physics-informed multiscale modeling along with a reduced set of empirical test data to reduce the time and cost of licensing and qualification of new nuclear fuel systems while maintaining the overall nuclear power plant safety case. SiC-SiC is being proposed for in-core applications, most notably fuel cladding, for current and next-generation nuclear reactors because of its high temperature stability, irradiation tolerance, and ability to withstand many accident conditions. As these composites exhibit multiscale architectures and complex microstructure-based fracture mechanics, it is an appealing use case for the AFQ methodology. While the end goal of this work is a single multiscale model that can be used for predictive in-core performance, current focus is on the individual various length scale models. Four individual models have been initially developed from microscale to engineering system level to capture key physics-based effects across different length scales. These models include a microscale homogenized tow model, a mesoscale fast Fourier transform–based weave model that integrates the homogenized tow model, a mesoscale finite element–based weave model, and a system-level BISON fuel performance model. Results of these models have undergone an initial comparison with separate-effects test data showing a good match to experimental results. By using the AFQ framework during model development, several near-term benefits have been secured including a reduction in development time for the SiC-SiC cladding, more targeted irradiation testing, and a better understanding of uncertainty.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Resource Assessment for Distributed Wind Energy: An Evaluation of Best-Practice Methods in the Continental US

Current wind resources within the United States (US) indicate a potential to profitably install nearly 1,400 gigawatts of distributed wind (DW) capacity. This amount is equivalent to over half of the United States’ current energy demand from electricity, making it enough to power millions of homes and businesses and replace countless fossil fuel-based generating plants. Despite the potential growth of DW in the US, deployments are presently hindered by a lack of confidence in resource estimation methods. One potential challenge is that smaller-scale turbines, with hub heights of 40 meters or less, are disproportionately impacted by obstacles such as buildings and vegetation. These obstacles may produce complex wake effects, best modeled with high-fidelity complex fluid dynamics (CFD) models that are too computationally expensive to use for routine siting and resource assessment. Thus, installers today make use of heuristics and simple equations to approximate the impact of obstacles while also leveraging long-term resource data from commercial or publicly available atmospheric models. This study evaluates these historical and commonly used methods alongside new lower-order obstacle models produced from CFD simulations and measurement-based bias correction. The preliminary results from this study show the importance of taking care in the choice and application of mesoscale atmospheric models and the significant value of bias correction using measurements from nearby meteorological towers. Detailed obstacle modeling provides only modest additional gains in performance and, in some cases, can add error, especially at sites where turbines have already been located to avoid obvious impact from upwind obstacles. These findings reinforce the importance of collecting in situ measurements and suggest that obstacle models may be better applied in practice to automated or computer-aided siting, rather than in economic wind resource assessments.

17 WIND ENERGY↗
DOE OSTI2025

Elucidating the Morphological Instability Mechanisms During Electrodeposition of Active Metals

The primary aim of this investigation is to elucidate the mechanisms by which a group of metals pertinent to battery technology exhibit unstable growth morphology during electrodeposition, which are crucial for the operational efficacy of batteries. The metals under consideration include alkali metals alongside zinc, all of which are characterized by high atomic mobility and reactivity with the electrolytes employed in battery systems. These materials represent promising electrode candidates for next-generation rechargeable batteries that are anticipated to deliver significantly enhanced energy density at a reduced cost. Nevertheless, the widespread occurrence of distinctly non-uniform morphologies, exemplified by the emergence of filamentous and moss-like structures, constitutes a significant impediment to their effective utilization, as it leads to accelerated capacity degradation and potential internal short circuits. This study addresses a critical question that has thus far garnered limited attention: Does the residual stress induced by the plating process play an instrumental role in inducing the morphological instability? Utilizing zinc as the primary model system, this project will be executed through three interconnected tasks: (1) the quantification of the relationship between plating stress and electroplating morphologies, (2) the three-dimensional imaging of the microstructure of the mossy structure at the nanoscale, and (3) the development of a mesoscale electroplating model to interpret empirical findings and generate additional predictive insights.

25 ENERGY STORAGE↗