Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer simulation”

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

Data-based filtered dissipation rate modelling for multi-modal turbulent combustion: evaluating a priori model generalizability

Manifold-based models offer a computationally efficient alternative to directly transporting the thermochemical state in computational simulations of turbulent reacting flows, projecting the high-dimensional thermochemical state-space onto a low-dimensional manifold. Recent efforts have yielded a manifold-based model applicable to multi-modal combustion, enabling reconstruction of the thermochemical state from solutions to two-dimensional manifold equations in mixture fraction and generalized progress variable that are parameterised by three scalar dissipation rates. In coarse-grained simulations such as Large Eddy Simulation (LES), closure of the multi-modal manifold equations and subfilter variances/covariance requires closure of three filtered scalar dissipation rates. Here, the present work adopts a data-based approach, providing closure for the three filtered scalar dissipation rates via deep neural networks (DNNs). High-fidelity datasets corresponding to an autoigniting n-dodecane jet flame and a bluff body swirl-stabilized confined lifted spray flame of two aviation fuels (Jet-A and C1) with different ignition propensities are leveraged to generate training data that spans a diverse range of thermodynamic conditions and combustion modes, including low- and high-temperature ignition regimes in addition to premixed and nonpremixed behaviour. A final DNN model is trained to enforce inherent physical constraints by learning nonlinear functional transformations of the three filtered scalar dissipation rates. The generalizability of this constrained DNN model is demonstrated a priori via conditional statistics evaluated on the lifted spray flame with C1–a configuration that had not been included in the training data. Excellent DNN agreement with conditional DNS statistics is observed, and integrated gradients are computed to identify the most sensitive input variables. The similarity of the marginal PDFs of the most informative input variables and outputs across configurations are quantified via the Wasserstein metric, demonstrating that data-based models may successfully generalize to unseen parametric conditions so long as the most informative input variables share similar distributions across training and testing datasets.

Data-based modelling↗
DOE OSTI2025

Elastic Bayesian Model Calibration

Functional data are ubiquitous in scientific modeling. For instance, quantities of interest are modeled as functions of time, space, energy, density, etc. Uncertainty quantification methods for computer models with functional response have resulted in tools for emulation, sensitivity analysis, and calibration that are widely used. However, many of these tools do not perform well when the computer model’s parameters control both the amplitude variation of the functional output and its alignment (or phase variation). This paper introduces a framework for Bayesian model calibration when the model responses are misaligned functional data. The approach generates two types of data out of the misaligned functional responses: (1) aligned functions so that the amplitude variation is isolated and (2) warping functions that isolate the phase variation. These two types of data are created for the computer simulation data (both of which may be emulated) and the experimental data. The calibration approach uses both types so that it seeks to match both the amplitude and phase of the experimental data. The framework is careful to respect constraints that arise, especially when modeling phase variation, and is framed in a way that it can be done with readily available calibration software. In conclusion, we demonstrate the techniques on two simulated data examples and on two dynamic material science problems: a strength model calibration using flyer plate experiments and an equation of state model calibration using experiments performed on the Sandia National Laboratories’ Z-machine.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Competition and cooperation: The plasticity of bacterial interactions across environments

Bacteria live in diverse communities, forming complex networks of interacting species. A central question in bacterial ecology is whether species engage in cooperative or competitive interactions. But this question often neglects the role of the environment. Here, we use genome-scale metabolic networks from two different open-access collections (AGORA and CarveMe) to assess pairwise interactions of different microbes in varying environmental conditions (provision of different environmental compounds). By computationally simulating thousands of environments for 10,000 pairs of bacteria from each collection, we found that most pairs were able to both compete and cooperate depending on the availability of environmental resources. This modeling approach allowed us to determine commonalities between environments that could facilitate the potential for cooperation or competition between a pair of species. Namely, cooperative interactions, especially obligate, were most common in less diverse environments. Further, as compounds were removed from the environment, we found interactions tended to degrade towards obligacy. However, we also found that on average at least one compound could be removed from an environment to switch the interaction from competition to facultative cooperation or vice versa. Together our approach indicates a high degree of plasticity in microbial interactions in response to the availability of environmental resources.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Understanding the Role of Hydroxyl Functionalization in Linear Poly(Ethylenimine) for Oxidation‐Resistant Direct Air Capture of CO 2

Aminopolymer-based adsorbents are a prominent class of materials being used for direct air capture of CO 2 at the industrial scale. However, improving their working lifetime, specifically by increasing their resilience to oxidative degradation, remains an ongoing challenge. Toward this end, functionalization of aminopolymers with non-amine functionalities such as hydroxyls has emerged in recent years as a promising strategy toward improving adsorbent lifetime. Although there is a growing body of work demonstrating the effectiveness of this approach and investigating the origin of this improved stability, studies to date have primarily focused on branched aminopolymer systems such as branched poly(ethylenimine). In this work, hydroxyl-functionalized linear poly(ethylenimine) is used to continue to probe the underlying protective mechanism of this strategy. A combination of thermogravimetric analysis, NMR relaxometry, differential scanning calorimetry, and computational simulations is used to better understand the relationship between the extent of chemical functionalization, physical properties, and adsorbent performance.

amine-based adsorbent↗
DOE OSTI2025

Experimental Pathways for Detecting Double Superionicity in Planetary Ices

The ice giant planets Uranus and Neptune are assumed to contain large amounts of planetary ices such as water, methane, and ammonia. The properties of mixtures of such ices at the extreme pressures and temperatures of planetary interiors are not yet well understood. Ab initio computer simulations have predicted that a number of ices exhibit a hydrogen superionic state and a doubly superionic state. Since the latter state has not yet been generated with experiments, we outline here two possible pathways for reaching and detecting such a state with dynamic compression experiments. Here, we suggest X-ray diffraction as the principal tool for detecting when the material becomes doubly superionic and the sublattice of one of the heavy nuclei melts. That would require a temperature of ~3500 K and pressures greater than ~200 GPa for H 3 NO 4 , which we use as an example material here. Such conditions can be reached with experiments that employ an initial shock that is followed by a ramp compression wave. Alternatively, one may use triple-shock compression because a single shock does not yield sufficiently high densities.

Neptune↗
DOE OSTI2024

Investigation of the notch sensitivity of tailorable long fiber discontinuous prepreg composite laminates

Tailorable discontinuous fiber composite laminates provide relative formability beyond that of continuous fiber laminates, while achieving improved mechanical performance over comparable stochastic systems. Here, in this work, the notch sensitivity of engineered prepreg platelet molded composite (PPMC) laminates is investigated using the open-hole tension (OHT) test and compared to available data for stochastic PPMCs and continuous fiber laminates made with the same material. The press-formed thermoplastic composites (AS4/PEKK) were molded with a quasi-isotropic stacking sequence. The discontinuous PPMC laminate was found to be notch insensitive with OHT strengths ranging from 145.4 MPa (CV $=$ 7%) for d/w $=$ 0.5 to 229.3 MPa (CV $=$ 9%) for d/w $=$ 0.25. The highly ordered meso-structure of the engineered PPMC laminate yields comparatively excellent mechanical properties for relatively thin laminates in contrast to stochastic systems. Both net- and gross-section failures were observed for d/w $=$ 0.25, which suggests that the engineered PPMC laminates studied here maintain a degree of inherent, internal stress concentrations that compete with those caused by geometric features such as a circular hole. Computational simulations of the OHT tests with explicitly represented platelets were found to be in good agreement with experimental measurements. The progressive failure analysis was used to conduct a numerical investigation of the stacking sequence and platelet meso-morphology.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spontaneous sodium ion storage behaviors of reduced graphene oxide anodes exceeding 100% Coulombic efficiency by modulated ion solvation

Rechargeable batteries are essential energy storage devices that power portable devices and electrical vehicles throughout the world. In general, it is thought that the electrochemical performance of rechargeable batteries is mostly determined by the electrodes within them and that the electrolyte plays a relatively passive role. However, ion transport and storage can be greatly influenced by the electrolyte solution structure, specifically, ion solvation within the bulk and ion desolvation across the electrode/electrolyte interfaces. Herein, we studied the role of the electrolyte as an active component of electrochemical energy storage devices. We found that with an appropriate electrolyte formulation, ion storage in disordered carbonaceous anode materials can occur spontaneously without externally supplied electrical energy. Reduced graphene oxide (RGO) in an ether-based electrolyte demonstrates ‘spontaneous' ion storage behaviors of adsorbing and inserting the solvated ions utilizing facilitated permeability and wettability of RGO, which results in Coulombic efficiency of ~145% due to additional charging capacity of ~180 mAh g -1 during electrochemical processes. The unexpected spontaneous ion storage behavior was extensively investigated using a combination of electrochemical analyses and diagnostics, advanced characterizations, and computational simulation. In conclusion, we believe the spontaneous ion storage behavior offers a new way to further improve the energy efficiency of practical rechargeable batteries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The Henkin-McGill critical temperature test for explosives revisited: An overview of historical and modern method development

Many overarching standards, regulations, or other requirements necessitate that the critical temperature (T crit ) of an energetic material be known or estimate prior to operations such as heating. However, they rarely, if ever, provide a specific method for doing so. While other methods exist to calculate T crit , such as the Frank-Kamenetskii equation, computer simulation, etc., these cannot be readily utilized for energetics that are physical mixtures (i.e. most pyrotechnics) or materials where detailed material properties required for such calculations are lacking. Here, in this study, a COTS (commercial off-the-shelf) SBAT (Simulated Bulk Auto-ignition Test) apparatus is modified to perform Henkin cookoff tests. This creates a simple, efficient, and cost-effective solution to estimate required T crit values for energetic materials. After a historical overview of the evolution of the Henkin test, several common energetic materials (PETN, RDX, etc.) for which T crit has been readily calculated, and historical data is available, were analyzed for comparison and verification purposes. This was followed by a variety of pyrotechnic mixtures and other materials to where calculation methods cannot be readily used. The modified apparatus, as well as an updated sample shell sealing method, produced results that aligned well with historical data. It also produced reasonable T crit estimates for those materials where the Frank-Kamenetskii equation cannot be applied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Tunable Nanostructures from Inverse Surfactants

Hierarchical materials in the natural world are often made through the self-assembly of amphiphilic molecules. Achieving similar structural complexity in synthetic materials requires understanding how various molecular parameters affect assembly behavior. In recent years, inverse surfactants─molecules with hydrophobic head groups and hydrophilic macromolecular tails─have been shown to self-assemble into supramolecular assemblies in aqueous solutions that show promise for a number of applications, including drug delivery. Here, we build an understanding of the morphological phase diagram of inverse surfactants using insights from scattering experiments, computer simulations, and statistical mechanics. The scattering and simulation results reveal that changing the headgroup size is an important molecular knob in controlling morphological transitions. The molecular size ratio of the hydrophobic group to the hydrophilic group emerges as a crucial dimensionless quantity in our theory and plays a determining role in setting the micelle structure and the transition from mesoscale to macroscale aggregates. Our minimal theory is able to qualitatively explain the key features of the morphological phase diagram, including the prevalence of fiber-like structures in comparison to spherical and planar micelles. Together, these findings provide a more complete picture of the molecular dependencies of assemblies of inverse surfactants, which we hope may aid in the de novo design of supramolecular structures.

Christakopoulos, Panagiotis [ORNL] (ORCID:00090004↗
DOE OSTI2026

Electrolyte-Dependent, “Microscopically Irreversible” H-Atom Transfer Kinetics of Ce-Based Metal–Organic Framework, Ce-MOF-808

Redox reactions at the interface of metal oxides and protic electrolytes almost always involve protons and electrons in equal amounts. Given the stoichiometry, these proton-coupled electron transfer (PCET) reactions are thermochemically equivalent to net H-atom transfer (HAT) reactions. The correlation between the chemical nature of solid catalysts and HAT kinetics has been employed for decades as the design principle for energy-relevant reactions (e.g., reactions of 2H + /H 2 ). More recently, chemists have experimentally determined that a change in liquid electrolytes that alters the microenvironment at the redox-active sites has an equally profound impact on electrocatalysis involving PCET/HAT. Yet, precise correlations between the chemical nature of electrolytes and the PCET kinetics are, to date, rare in the literature. Herein, we report our findings using the Ce-based metal−organic framework, Ce-MOF-808, as a model system. Each Ce 6 (μ 3 −O) 4 (μ 3 − OH) 4 (OH) 6 (H 2 O) 6 node of this MOF undergoes a 1H + /1e − redox reaction. Using chronoamperometry and the Cottrell analysis, we have determined that the PCET hopping kinetics within the pores of Ce-MOF-808 can change by orders of magnitude by altering the buffer species and the proton activity of the electrolyte. Furthermore, in all buffers, reductive reactions were ∼3−10 times faster in kinetics than the reverse oxidative reaction with the same electrochemical driving force, suggesting that the system, at first glance, violates the principle of microscopic reversibility. Isothermal titration calorimetry (ITC) and computational simulations corroborated that the buffer-node binding thermodynamics are quite distinct, depending on the chemical nature of the buffer and the oxidation state of the node. Together, these results suggest that the substrate and the product during the oxidative vs reductive reaction of Ce-MOF-808 are chemically different species, which explains the apparent ‘microscopic irreversibility.’ Thus, the rational modulation of electrolytes can dramatically enhance PCET kinetics, even though the solid electrodes remain identical. Implications of these findings are contrasted with the electrochemical/electrocatalytic behavior of other redox-active MOFs, heterogeneous catalysts, and enzymatic systems at the solid−liquid interface.

Ce-based MOF↗
DOE OSTI2026

Effect of Lithium Doping on MgO Hydroxylation and Carbonation

Recovery of magnesium from brines can potentially be used to source MgO (periclase) as a CO 2 sorbent or for Mg-based cements. However, it is not clear how common impurities in brines, such as lithium, affect the resulting MgO reactivity. Here, to test the effect of lithium incorporation on MgO reactivity for hydration and carbonation, we combined computational simulations with experiments. Experimentally altered (Mg,Li)O with a low dopant concentration (0.012 ± 0.002% w/w Li) was characterized using synchrotron-based X-ray scattering and high-resolution electron microscopy to measure reaction layer formation on (Mg,Li)O. Single-crystal X-ray diffraction analysis of (Mg,Li)O demonstrates that the incorporation of lithium leads to the formation of oxygen vacancies. The presence of vacancies is likely causing faster hydroxylation rates as predicted by ab initio molecular dynamics simulations. However, the faster hydroxylation rates likely lead to faster passivation of the surface because we observe thinner reaction layers on (Mg,Li)O samples both over short time periods (30 days) and over long time periods (28 years). After 28 years, the reaction layer on the (Mg,Li)O sample was less than one-third of the thickness of that of the pure MgO sample. In addition, over 30 days, reaction layers on (Mg,Li)O samples primarily formed at steps rather than on terraces, in contrast to our previous observations on MgO. Based on our results, naturally occurring impurities in MgO modify its reactivity even at very low concentrations and need to be considered for accurate reaction rate prediction for application of MgO as a CO 2 sorbent or in cements.

MgO↗
DOE OSTI2024

Electrochemically Determined and Structurally Justified Thermochemistry of H atom Transfer on Ti-Oxo Nodes of the Colloidal Metal–Organic Framework Ti-MIL-125

Titanium dioxide (TiO 2 ) has long been employed as (photo)electrodes for reactions relevant to energy storage and renewable energy synthesis. Proton-coupled electron transfer (PCET) reactions with equimolar amounts of protons and electrons at the TiO 2 surface or within the bulk structure lie at the center of these reactions. Because a proton and an electron are thermochemically equivalent to an H-atom, these reactions are essentially H-atom transfer reactions. Thermodynamics of H-atom transfer has a complex dependence on the synthetic protocol and chemical history of the electrode, the reaction medium, and many others; together, these complications preclude the understanding of the H-atom transfer thermochemistry with atomic-level structural knowledge. Herein, we report our success in employing open-circuit potential (E OCP ) measurements to quantitatively determine the H-atom transfer thermochemistry at structurally well-defined Ti-oxo clusters within a colloidally stabilized metal–organic framework (MOFs), Ti-MIL-125. The free energy to transfer H-atom, Ti 3+ O–H bond dissociation free energy (BDFE), was measured to be 68(2) kcal mol -1 . To the best of our understanding, this is the first report on using E OCP measurements to quantify thermochemistry on any MOFs. The proton topology, the structural change upon the redox reaction, and BDFE values were further quantitatively corroborated using computational simulations. Furthermore, comparisons of the E OCP -derived BDFEs of Ti-MIL-125 to similar parameters in the literature suggest that E OCP should be the preferred method for quantitatively accurate BDFE calculations. Here, the reported success in employing E OCP for nanosized Ti- MIL-125 should lay the ground for thermochemical measurements of other colloidal systems, which are otherwise challenging. Implications of these measurements on Ti-MIL-125 as an H-atom acceptor in chemical reactions and comparisons with other MOFs/metal oxides are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Interactions of Polar and Nonpolar Groups of Alcohols in Zeolite Pores

Understanding the quantitative interactions among zeolite pore walls, Bro̷nsted acid sites, and molecules with both polar and nonpolar regions is essential for scoping out the potential of zeolites as sorbents and catalysts. Purely siliceous zeolites (MFI and Beta in the present study) are hydrophobic, whereas those containing aluminum are considered hydrophilic, preferentially adsorbing organic molecules even in aqueous environments. To characterize these interactions, we use primary alcohols of increasing molecular weight, quantifying their specific interactions in the confined pore space of the alkyl (CH x ) and OH groups. Three types of interactions were identified: (i) alkyl CH x groups interacting with the zeolite pore walls (approximately 10 kJ mol −1 per carbon), (ii) alcohol OH groups interacting with the pore walls (30−35 kJ mol −1 ), and (iii) alcohol OH groups interacting with Bro̷nsted acid sites (37 kJ mol −1 ). All three interactions were well mirrored by computational simulations. The contribution of the alkyl CH x groups was inferred from the incremental increase in sorption enthalpy with increasing molecular weight; the interaction strength of the OH groups was determined by extrapolating the global adsorption enthalpy of the alcohols to a hypothetical OH group without an alkyl group. This value was identical to the adsorption enthalpy of water. The experiments demonstrated that only water has an adsorption enthalpy on zeolite pore walls lower than its condensation enthalpy (30−35 kJ mol −1 vs 45 kJ mol −1 ), limiting the concentration of water that can be adsorbed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Size-Dependent Optical Band Gaps in Metal–Organic Framework Nanoparticles

Decades of research into size-dependent semiconductor optical gaps have focused on quantum confinement as the dominant mechanism. Emerging reports indicate that lattice strain─intentional or incidental─can impart optical shifts similar or greater in magnitude. Here, we report evidence of optical absorption and photoluminescence spectra of M(1,2,3-triazolate)2 (M = Mg, Cr, Mn, Fe, Co, Cu, Zn, or Cd) nanoparticles that blueshift from bulk values with decreasing particle sizes in a manner that defies explanation by conventional quantum confinement. Here, the phenomenon persists for particle sizes as large as 200 nm, whereas quantum confinement generally ceases beyond 20–30 nm diameters and follows a weaker dependence on the particle radius. Computational simulations and crystallographic analysis suggest that this behavior arises from size-dependent changes to metal–linker bonding that manifest in strain values comparable to literature reports of strain-induced optical shifts in other classes of materials. This behavior appears beyond this family of materials in other notable examples of metal–organic frameworks (MOFs), including the well-studied Cu3(trimesate)2 (CuBTC), where smaller sizes correlate with blueshifted optical gaps. Taken together, these results represent one of the few examples of size-dependent strain in crystalline materials and reinforce the emerging view that MOFs become softer materials when isolated as nanoparticles.

Electrical conductivity↗
DOE OSTI2024

Electrochemical formation of bis(fluorosulfonyl)imide-derived solid-electrolyte interphase at Li-metal potential

Lithium bis(fluorosulfonyl)imide-based liquid electrolytes are promising for realizing high Coulombic efficiency and long cycle life for next-generation Li-metal batteries. However, the role of anions in the formation of solid-electrolyte interphase remains unclear. Herein, we combine electrochemical analyses and X-ray photoelectron spectroscopy measurements both with and without sample washing, together with computational simulations, to propose the reaction pathways of electrolyte decomposition and correlate the interphase component solubility with the efficacy of passivation. Additionally, we discovered that not all the products derived from interphase-forming reactions were incorporated into the resulting passivation layer, with a notable portion present in the liquid electrolyte. Further, we found that the high-performance electrolytes can afford a sufficiently passivating interphase with minimized electrolyte decomposition, by incorporating more anion-decomposition products. Overall, this work presents a systematic approach of coupling electrochemical and surface analyses to paint a comprehensive picture of solid-electrolyte interphase formation, while identifying key attributes of high-performance electrolytes for guiding future designs.

Yu, Weilai [Stanford University, CA (United States↗
DOE OSTI2026

Lithium metal-mediated electrochemical reduction of per- and poly-fluoroalkyl substances

Per- and poly-fluoroalkyl substances (PFAS) have substantial environmental and health hazards. Unfortunately, current degradation routes require high temperatures or corrosive conditions and/or lead to incomplete defluorination and the generation of shorter alkyl chains. Inspired by the lithium-metal battery literature, here we develop an electrochemical degradation process that leverages reactive metals and highly reducing environments. Here, we show that electrodeposited lithium metal can enable 95% degradation and 94% defluorination of perfluorooctanoic acid to LiF without forming any shorter C 2 –C 6 PFAS as end products. Using computational simulations, we find that electron transfer from lithium to perfluorooctanoic acid leads to rapid C–F bond cleavage, fluoride formation and carbon chain fragmentation. We expand the scope to other PFAS compounds and demonstrate substantial degrees of degradation on over 22 different PFAS, plus complete mineralization to inorganic fluorides. Finally, we use the mineralized F − as a fluorine source for the synthesis of fluorinated non-PFAS compounds to complete a circular fluorine loop from waste to valuable product.

Sarkar, Bidushi [Univ. of Chicago, IL (United Stat↗
DOE OSTI2025

Space-time crystals from particle-like topological solitons

Time crystals are unexpected states of matter that spontaneously break time-translation symmetry either in a discrete or continuous manner. However, spatially mesoscale space-time crystals that break both space and time symmetries have not been reported. Here we report a continuous space-time crystal in a nematic liquid crystal driven by ambient-power, constant-intensity unstructured light. Our numerically constructed four-dimensional configurations exhibit good agreement with these experimental findings. Although meeting the established criteria to identify time-crystalline order, both experiments and computer simulations reveal a space-time crystallization phase formed by particle-like topological solitons. The robustness against temporal perturbations and spatiotemporal dislocations shows the stability and rigidity of the studied space-time crystals, which relates to their locally topological nature and many-body interactions between emergent spontaneously twisted, particle-like solitonic building blocks. Their potential technological utility includes optical devices, photonic space-time crystal generators, telecommunications and anti-counterfeiting designs, among others.

Liquid crystals↗