Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiscale analysis”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Uniaxial Tensile Properties of AS4 3D Woven Composites with Four Different Resin Systems: Experimental Results and Analysis: Property Computations

As a part of the NASA Composite Technology for Exploration project, eight different AS4 3D orthogonal woven composite panels were manufactured and were subjected to mechanical testing including uniaxial tension along the weaves' warp direction. Each set, with four different resin systems (KCR-IR6070, EP2400, RTM6, and RS-50), included weave architectures designed using 12K and 6K AS4 carbon fiber yarns. For the tension testing conducted at Room Temperature Ambient (RTA) conditions, the elastic modulus and strength of these eight panels (as-processed and thermally-cycled) were measured and compared while the potential evolution of micro-cracking before and after thermal cycling were monitored via optical microscopy and X-Ray Computed Tomography. The data set also included test results of the as-processed materials at Elevated Temperature Wet (ETW) conditions. In the second part of this study, efforts were made to compute elastic constants for AS4 6K/RTM6 and AS4 12K/RTM6 materials by implementing a finite element approach and the Multiscale Generalized Method of Cells (MSGMC) technique developed at NASA Glenn Research Center. Digimat-FE was used to model the weave architectures, assign properties, calculate yarn properties, create the finite element mesh, and compute the elastic properties by applying periodic boundary conditions to finite element models of each repeating unit cell. The required input data for MSGMC was generated using Matlab® from Digimat exported weave information. Experimental and computational results were compared, and the differences and limitations in correlating to the test data were briefly discussed.

Property Computations↗
DOE OSTI2026

Molecular insights into CO 2 -to-bicarbonate transformation in functionalized anion exchange ionomers for electrochemical separations

Bipolar membrane (BPM) electrochemical processes are a promising platform for carbon dioxide (CO 2 ) separations, but the molecular level thermodynamic and kinetic understanding of CO 2 -to-bicarbonate (HCO 3 − ) transformation remain poorly understood. This study employs a multiscale computational approach to systematically explore the adsorption and reactive transformation of CO 2 in five anion exchange ionomer systems. Classical molecular dynamics (MD) simulation results demonstrate that polymers with imidazolium groups significantly reduce CO 2 diffusion and enhance (OH − )–CO 2 interactions due to stronger electrostatic and π-interactions. Compared to the commonly used quaternary ammonium ionomers, imidazolium-functionalized ionomers show improved CO 2 proximity and interaction strength. Ab initio MD and density functional theory (DFT) calculations reveal that the benzyl-substituted imidazolium (IM-Ben) substantially reduces the energy barrier for HCO 3 − formation (∼72 meV lower) compared to the alkyl-substituted IM-nBu, while also mitigating imidazolium deprotonation under moderate hydration conditions. Transition state analysis shows IM-Ben forms more extensive hydrogen-bonding networks, which stabilize the transition state structure and contribute to a lower energy barrier for bicarbonate formation. These findings highlight the advantage of the adjacent benzyl moiety in enabling efficient CO 2 -to-bicarbonate transformation via hydrated hydroxide ion counterions, offering mechanistic insights and clear molecular design principles for optimizing anion exchange ionomers at bipolar membrane interfaces for electrochemical CO 2 separation applications.

Bipolar membranes, Reactive transformation of CO2,↗
NASA NTRS2008

New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.

Saether, E.↗
NASA NTRS2011

Analysis of Ultra High Resolution Sea Surface Temperature Level 4 Datasets

Sea surface temperature (SST) studies are often focused on improving accuracy, or understanding and quantifying uncertainties in the measurement, as SST is a leading indicator of climate change and represents the longest time series of any ocean variable observed from space. Over the past several decades SST has been studied with the use of satellite data. This allows a larger area to be studied with much more frequent measurements being taken than direct measurements collected aboard ship or buoys. The Group for High Resolution Sea Surface Temperature (GHRSST) is an international project that distributes satellite derived sea surface temperatures (SST) data from multiple platforms and sensors. The goal of the project is to distribute these SSTs for operational uses such as ocean model assimilation and decision support applications, as well as support fundamental SST research and climate studies. Examples of near real time applications include hurricane and fisheries studies and numerical weather forecasting. The JPL group has produced a new 1 km daily global Level 4 SST product, the Multiscale Ultrahigh Resolution (MUR), that blends SST data from 3 distinct NASA radiometers: the Moderate Resolution Imaging Spectroradiometer (MODIS), the Advanced Very High Resolution Radiometer (AVHRR), and the Advanced Microwave Scanning Radiometer ? Earth Observing System(AMSRE). This new product requires further validation and accuracy assessment, especially in coastal regions.We examined the accuracy of the new MUR SST product by comparing the high resolution version and a lower resolution version that has been smoothed to 19 km (but still gridded to 1 km). Both versions were compared to the same data set of in situ buoy temperature measurements with a focus on study regions of the oceans surrounding North and Central America as well as two smaller regions around the Gulf Stream and California coast. Ocean fronts exhibit high temperature gradients (Roden, 1976), and thus satellite data of SST can be used in the detection of these fronts. In this case, accuracy is less of a concern because the primary focus is on the spatial derivative of SST. We calculated the gradients for both versions of the MUR data set and did statistical comparisons focusing on the same regions.

sea surface temperature SST)↗
DOE OSTI2026

Prelithiated SiO x /Graphite-NMC811 Cells: Capacity Loss, Impedance Rise and Hidden Degradation Pathways Revealed Using Three-Electrode Diagnostics

Electrodes containing SiO x /graphite (Gr) materials are attractive as anodes for high-energy lithium-ion batteries. However, their mechanical deformation, electrochemical response, and impedance evolution during long-term cycling are strongly coupled, complicating accurate diagnosis of performance fade mechanisms. In this work, the behavior of electrochemically prelithiated SiO x /Gr anodes paired with NMC811 cathodes is systematically investigated using techniques that include in-situ dilatometry, three-electrode electrochemistry, and multiscale post-cycling microscopy. The SiO x /Gr electrode exhibits a maximum expansion of 49% upon lithiation to 10 mV vs Li + /Li, with 84% of the expansion and 91% of the capacity being reversible. In full cells, relatively stable cycling with only 12% capacity fade over 500 cycles is observed. Three-electrode experiments reveal cell-level impedance growth, which is dominated by the NMC811 cathode: the SiO x /Gr anode exhibits minimal net impedance rise and an initial impedance decrease at low potentials. Despite this apparent electrochemical stability, cross-sectional SEM, PFIB tomography, and cryo-STEM reveal irreversible anode thickening caused by the accumulation of an inorganic-rich solid electrolyte interphase (SEI) permeating the anode bulk. Electrode potential-shift analysis further demonstrates that Li + ions released from lithium reservoirs in the prelithiated anode mask true lithium inventory loss during aging. These results demonstrate that low-expansion SiO x /Gr anodes can simultaneously exhibit favorable cycling and impedance metrics while undergoing substantial, hidden degradation, underscoring the importance of electrode-resolved diagnostics for evaluating prelithiated silicon-based anodes.

25 ENERGY STORAGE↗
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

Structural complexity of snapshots of two-dimensional Fermi-Hubbard systems

The development of quantum gas microscopy for two-dimensional optical lattices has provided an unparalleled tool to study the Fermi-Hubbard model (FHM) with ultracold atoms. Spin-resolved projective measurements, or snapshots, have played a significant role in quantifying correlation functions, theory verification, and thus the uncovering of underlying physical phenomena such as antiferromagnetism at commensurate filling on bipartite lattices and other charge and spin correlations, as well as dynamical properties at various densities. Here we employ a recent concept, the multiscale structural complexity, and show that when computed for the snapshots (of single spin species, local moments, or total density) it can provide a theory-free property, immediately accessible to experiments. Specifically, after benchmarking results for Ising and $XY$ models, we study the structural complexity of snapshots of the repulsive FHM in the two-dimensional square lattice as a function of doping and temperature. We generate projective measurements using determinant quantum Monte Carlo and compare their complexities against those from the experiment. We demonstrate that these complexities are linked to relevant physical observables such as the entropy and double occupancy. Their behaviors capture the development of correlations and relevant length scales in the system. Furthermore, we provide an open-source code in python which can be implemented into data analysis routines in experimental settings for the square lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Tracking animal movements via collaborative acoustic telemetry networks: Multiscale habitat use, phenology, and management insights

Abstract Estuaries support diverse fish and invertebrate communities, including resident species that rely on estuarine habitats year‐round and transient migratory species. The unique movement patterns of these animals connect habitats within and far beyond the estuary and are integrally linked to fisheries management objectives. With a focus on Chesapeake Bay, this study leveraged data from collaborative acoustic telemetry networks in the northwest Atlantic to assess habitat use and phenology of movements for seven species of fish (cownose rays, dusky sharks, smooth dogfish, alewife, striped bass, common carp, and blue catfish) and one invertebrate (horseshoe crabs). A total of 288 acoustically tagged individuals were detected >3.2 million times (6,743 to 2,095,717 detections per species) on receivers across ~20.5 degrees of latitude spanning the North American Atlantic seaboard from Florida, USA, to New Brunswick, Canada. Common metrics of movement and phenology grouped these species as resident (common carp, blue catfish, horseshoe crabs), primarily resident in estuaries (juvenile striped bass), and coastal migrant (cownose rays, dusky sharks, smooth dogfish, alewife); maximum distance traveled varied by three orders of magnitude among these species. Further analysis of phenology for coastal migrants elucidated the timing and duration of these species' use of Chesapeake Bay. Collectively, movements linked habitats within Chesapeake Bay and connected the estuary to coastal ecosystems both to the north (e.g., alewife) and south (e.g., cownose rays), creating networks of fisheries management jurisdictions that varied in complexity and identified opportunities for enhancement to current management or co‐management of some species. Our results elucidate the importance of estuaries to species with diverse movement behaviors, identify scales and pathways of habitat connectivity via animal movements, and highlight the utility of collaborative acoustic telemetry networks for quantifying movements relevant to both ecological research and fisheries management.

Livernois, Mariah C.↗
DOE OSTI2024

Dynamic data-driven multiscale modeling for predicting the degradation of a 316L stainless steel nuclear cladding material

Here, we have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ionic Interdiffusion at Cathode|Solid-Electrolyte Interface: A Machine Learning–Assisted Multiscale Investigation and Mitigation Strategies

Future lithium batteries are expected to use solid electrolytes to achieve higher energy density and fast charge capabilities. However, most solid electrolytes are thermodynamically unstable against layered oxide cathodes. In this study, the stability of LiCoO2 (LCO) cathode with Li10GeP2S12 (LGPS) solid electrolyte is investigated using ab initio molecular dynamics (AIMD) and machine learning molecular dynamics (MLMD). The propensity of ionic interdiffusion, formation of a passivating interphase layer, and corresponding decay in cell performance is addressed using a continuum model. Large-scale MLMD simulations confirm that the LCO|LGPS interface permits interdiffusion of cobalt (Co) and other ionic species, leading to the formation and growth of a resistive interphase and to dramatic capacity fade even in the first cycle. We examine the literature evidence that incorporating a thin layer of LiNb0.5Ta0.5O3 (LNTO) between LCO and LGPS prevents the interdiffusion of ions. Atomistic simulations suggest that substituting lithium (Li) in LNTO with Co is thermodynamically unfavorable, thereby inhibiting ionic interdiffusion. The stable Nb5+/Ta5+ states form a rigid metal-oxide framework, which consequently also prevents the substitution of niobium (Nb) or tantalum (Ta). However, continuum-level analysis suggests that the higher mechanical stiffness of LNTO can lead to interfacial delamination between the LCO and LNTO. This phenomenon reduces the effectiveness of the protective layer. This paper, therefore, highlights the need to develop novel interlayers that balance low ionic interdiffusion with low mechanical stiffness.

Ncube, Musawenkosi K.↗
NASA NTRS2020

An Algorithm for Characterization of Fiber Aggregation in Composite Microstructures

Composite structures are susceptible to localized flaws, or variability, that drive global failure. To capture this variance, a multiscale model needs to be introduced that not only accurately represents the statistical nature of the composite microstructure but is also efficient. At the microscale, fiber aggregation creates local stress concentrations, where failure is likely to occur sooner than expected. A method of rapid microstructure generation, in conjunction with cluster characterization, is examined to develop accurate reproductions of 2D composite cross-sections. Parameterization of shape and size of fiber clusters is used to characterize representative volume elements. The discrete element method is used to generate pseudo-microstructures. The clustering parameters from the pseudo- and actual microstructure arrangements (obtained from micrographs) can be compared to determine the validity of the representative volume element generated with the discrete element method. Results show a promising approach to evaluating randomness of fiber distributions and how to accurately recreate microstructures for strength analysis.

Carbon fiber↗
DOE OSTI2025

Implications of AI and Data Center Service Demand on Energy, Water, and Economic Systems

This talk focuses on the challenges posed by the rapid growth of data centers (DCs), which rely heavily on electricity, water, land, and mineral resources, and have significant interactions with economic growth. The U.S. leads in DC deployment, but understanding how resource demands and constraints impact DC expansion requires strategic multisector analysis. To address this, the talk introduces a new conceptual framework and modeling capability within GCAM-USA, allowing integrated analysis of DC service demand, electricity, and direct and indirect water needs under varied scenarios. Insights include state-level water and electricity trends for DCs, uncertainty bounds, and the broader implications for resources and economics.

Data Centers↗
NASA NTRS2021

Direct Multipoint Observations Capturing the Reformation of a Supercritical Fast Magnetosonic Shock

Using multipoint Magnetospheric Multiscale (MMS) observations in an unusual string-of-pearls configuration, we examine in detail observations of the reformation of a fast magnetosonic shock observed on the upstream edge of a foreshock transient structure upstream of Earth’s bow shock. The four MMS spacecraft were separated by several hundred km, comparable to suprathermal ion gyro-radius scales or several ion inertial lengths. At least half of the shock reformation cycle was observed, with a new shock ramp rising up out of the “foot” region of the original shock ramp. Using the multipoint observations, we convert the observed time-series data into distance along the shock normal in the shock’s rest frame. That conversion allows for a unique study of the relative spatial scales of the shock’s various features, including the shock’s growth rate, and how they evolve during the reformation cycle. Analysis indicates that: the growth rate increases during reformation, electron-scale physics play an important role in the shock reformation, and energy conversion processes also undergo the same cyclical periodicity as reformation. Strong, thin electron-kinetic-scale current sheets and large-amplitude electrostatic and electromagnetic waves are reported. Evidence is also presented of nonlinear wave decay from electromagnetic whistler-mode “lion roars” to electrostatic solitary waves in the downstream plasma regime. Results highlight the critical cross-scale coupling between electron-kinetic-and ion-kinetic-scale processes and details of the nature of nonstationarity, shock-front reformation at collisionless, fast magnetosonic shocks.

D L Turner↗
DOE OSTI2024

Investigating the Crosslinking, Degradation, and Adhesion Behavior of Photovoltaic Encapsulants Under Thermal Accelerated Aging

Degradation of photovoltaic (PV) module encapsulant characteristics that lead to mechanical embrittlement and delamination remains a cause of failure in solar installations. A multiscale reliability model connecting the encapsulant mechanical and fracture properties to the degraded molecular structure and interfacial bonding to adjacent solar cell and glass substrates was previously published. The model, developed primarily for poly(ethylene-co-vinyl acetate) acetate (EVA) encapsulants, remains to be experimentally validated. Determining the degradation and crosslinking kinetics of alternative encapsulants, such as polyolefin elastomer (POE) and EVA/POE/EVA composites (EPE), can generalize the model. In this work, we subject fully cured EVA, POE, and EPE encapsulants to accelerated thermal aging to determine how high temperatures impact reaction kinetics. An increase in gel content (crosslinking) and decrease in crystallinity of the encapsulants under hot-aerobic (90 degrees C, 22% RH) and hot-anaerobic (90 degrees C, sealed in N 2 air) aging were observed, even in the absence of UV and crosslinking initiators. Fourier transform infrared spectroscopy (FTIR)-attenuated total reflectance analysis showed insignificant encapsulant degradation, demonstrating the critical role of UV and moisture in accelerating degradation. Adhesion testing performed on coupon-level specimens (cell/encapsulant/glass laminates) showed decreases in adhesion energy, Gc , from 5000 h of hot-dry (90 degrees C, ~1% RH) and hot-humid (90 degrees C, 60% RH) aging. POE coupons demonstrated the best stability, followed by EPE then EVA. For EVA and POE, hot-humid aged coupons experienced a larger decrease in Gc due to enhanced hydrolytic degradation. Hot-dry aging condition demonstrated that thermal degradation of the interface could be significant even if the encapsulant experiences negligible degradation in the absence of UV and elevated humidity.

accelerated aging↗
NASA NTRS2016

Tool for Generation of MAC/GMC Representative Unit Cell for CMC/PMC Analysis

This document describes a recently developed analysis tool that enhances the resident capabilities of the Micromechanics Analysis Code with the Generalized Method of Cells (MAC/GMC) 4.0. This tool is especially useful in analyzing ceramic matrix composites (CMCs), where higher fidelity with improved accuracy of local response is needed. The tool, however, can be used for analyzing polymer matrix composites (PMCs) as well. MAC/GMC 4.0 is a composite material and laminate analysis software developed at NASA Glenn Research Center. The software package has been built around the concept of the generalized method of cells (GMC). The computer code is developed with a user friendly framework, along with a library of local inelastic, damage, and failure models. Further, application of simulated thermomechanical loading, generation of output results, and selection of architectures to represent the composite material have been automated to increase the user friendliness, as well as to make it more robust in terms of input preparation and code execution. Finally, classical lamination theory has been implemented within the software, wherein GMC is used to model the composite material response of each ply. Thus, the full range of GMC composite material capabilities is available for analysis of arbitrary laminate configurations as well. The primary focus of the current effort is to provide a graphical user interface (GUI) capability that generates a number of different user-defined repeating unit cells (RUCs). In addition, the code has provisions for generation of a MAC/GMC-compatible input text file that can be merged with any MAC/GMC input file tailored to analyze composite materials. Although the primary intention was to address the three different constituents and phases that are usually present in CMCs-namely, fibers, matrix, and interphase-it can be easily modified to address two-phase polymer matrix composite (PMC) materials where an interphase is absent. Currently, the tool capability includes generation of RUCs for square packing, hexagonal packing, and random fiber packing as well as RUCs based on actual composite micrographs. All these options have the fibers modeled as having a circular cross-sectional area. In addition, a simplified version of RUC is provided where the fibers are treated as having a square cross section and are distributed randomly. This RUC facilitates a speedy analysis using the higher fidelity version of GMC known as HFGMC. The first four mentioned options above support uniform subcell discretization. The last one has variable subcell sizes due to the primary intention of keeping the RUC size to a minimum to gain the speed ups using the higher fidelity version of MAC. The code is implemented within the MATLAB (The Mathworks, Inc., Natick, MA) developmental framework; however, a standalone application that does not need a priori MATLAB installation is also created with the aid of the MATLAB compiler.

Materials Engineering↗
DOE OSTI2026

Detecting and Characterizing Fracture Zones Using a Convolutional Neural Network

This project directly supports the Geothermal Technologies Office (GTO) objectives outlined in the Multi-Year Program Plan (MYPP) by advancing two key research areas: “Exploration and Characterization” and “Data, Modeling, and Analysis.” This project has successfully demonstrated a pre-drilling ability to image and characterize the distribution and connectivity of subsurface faults and fractures, key parameters for identifying permeable pathways that enable geothermal fluids to circulate and produce energy. Specifically, we developed and implemented innovative machine learning methodologies to enhance geothermal exploration. Large-scale faults were detected using a Convolutional Neural Network (CNN), while small-scale fractures were characterized using a novel Double-Beam Neural Network (DBNN). These tools have proven both technically effective and cost-efficient by reducing reliance on expensive exploratory drilling. Through collaboration with our geothermal industry partner, this research has significantly advanced techniques for identifying hidden geothermal systems and extending the productive lifespan of existing geothermal fields. We applied our methods to two geothermal fields—Soda Lake (Nevada) and Lightning Dock (New Mexico)—to identify shallow steam-charged fracture zones and characterize deep faults at depths of 1.5-2 km. The steam zone identified at the Soda Lake geothermal field showed excellent agreement with prior drilling data, validating the effectiveness of our approaches. In addition, the analysis revealed three new prospective drilling targets for further development and verification. The outcomes of this project improve our scientific understanding of geothermal reservoir behavior, enhance exploration efficiency, extend the economic life of existing geothermal plants. Ultimately, these advancements contribute to GTO’s goal of achieving more sustainable, affordable, and data-driven geothermal energy development across the United States.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗
DOE OSTI2024

Machine-learning-assisted deciphering of microstructural effects on ionic transport in composite materials: A case study of Li 7 La 3 Zr 2 O 12 -LiCoO 2

The effective diffusivity of ionic species in multiphase materials is critical for the design and function of composite materials for electrochemical energy storage. In practice, effective diffusivity depends sensitively not only on the intrinsic diffusivities of constituting materials but also on their topological arrangement; nevertheless, these coupled contributions are oversimplified in most analytical models. Here, we combine atomistically informed mesoscale modeling and machine learning (ML) analysis to unravel how such features affect effective diffusivity in two-phase composites. Using the Li 7 La 3 Zr 2 O 12 -LiCoO 2 composite solid-state battery cathode as a model system, we compute effective diffusivity for 600 distinct dense polycrystalline microstructures with different topological configurations of grains, grain boundaries, and heterointerfaces. We verify that in addition to atomic-scale variabilities, microstructural feature diversity can significantly impact effective transport properties. Across the ensemble of test microstructures, this often results in bimodal distributions of effective diffusivity that encompass two qualitatively distinct operating mechanisms, which we identify via flux analysis. An ML approach reveals that the most critical determining factors for effective diffusivity are the connectivity of bulk phases and their heterointerfaces. The role of ionic mobility at the heterointerfaces is also discussed. These insights highlight the combined importance of microstructure and interface engineering in tuning the transport properties of ionic species in composite materials. In conclusion, our framework can also be extended for understanding generic microstructure-property relationships in other complex multiphase materials.

25 ENERGY STORAGE↗