Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer simulations”

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 541 records · Page 30

Evaluating the Importance of Conformers for Understanding the Vacuum-Ultraviolet Spectra of Oxiranes: Experiment and Theory

Vacuum-ultraviolet (VUV) absorption spectroscopy enables electronic transitions that offer the unambiguous identification of molecules. As target molecules become more complex, multifunctional species present a great challenge to both experimental and computational spectroscopy. This research reports both experimental and theoretical studies of oxiranes. Computationally, the nuclear ensemble approach has been used to accurately predict experimental spectra for a variety of molecules. However, this approach incurs great computational cost, as ensembles generally consist of thousands of geometries. The present study aims to drastically reduce the ensemble by evaluating the significance of the conformers to the predicted spectra. This approach was applied to 11 substituted oxiranes using the Conformer Rotamer Ensemble Sampling Tool (CREST) of Grimme to generate an ensemble of unique conformers determined by their Boltzmann populations. Five TD-DFT functionals (BMK, CAM-B3LYP, M06-2X, MN15, ωB97X-D) and EOM-CCSD were used to simulate the spectrum of each substituted oxirane ensemble. Computed spectra were then compared to the experiment using both qualitative and quantitative metrics. Based on these metrics, it was observed that certain conformers may not be necessary to characterize this set of oxiranes despite the temperature (323 K) of the experiment. A single conformer can then be used with TD-DFT and EOM-CCSD to replicate the experimental spectra of these medium-sized combustion species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine Learning Interatomic Potentials for Modeling Framework Flexibility and Water Uptake in NbOFFIVE-1-Ni Metal–Organic Framework

Metal–organic frameworks (MOFs), with their distinctive porous structures and tunable chemical properties, have shown immense promise in the separation and storage of gases. Currently, the accurate simulation of their adsorptive properties remains challenging, especially for systems where the molecules fit very tightly into the pores. Traditional simulation methods often approximate the frameworks as rigid and do not account for the framework flexibility seen in materials such as NbOFFIVE-1-Ni. First-principles molecular dynamics (FPMD) simulations offer the desired accuracy in modeling this flexibility but are limited by their extensive computational demands, rendering them impractical for long simulations. Conversely, classical force field-based simulations offer computational efficiency but lack the necessary accuracy. Here, to break this accuracy-efficiency trade-off, we have developed machine learning interatomic potentials trained on energies and forces from FPMD to model the framework flexibility of NbOFFIVE-1-Ni in the presence of water over nanosecond time scales. Furthermore, by integrating MLIP-driven molecular dynamics (MLIP-MD) with grand canonical Monte Carlo (GCMC) simulations, we further incorporated framework flexibility into adsorption predictions, yielding water adsorption isotherms that better align with experimental data compared to those of conventional GCMC simulations. These advances offer new opportunities for the design and optimization of MOFs in gas storage and separation applications.

adsorption↗

Fully quantum algorithm for mesoscale fluid simulations with application to partial differential equations

Fluid flow simulations marshal our most powerful computational resources. In many cases, even this is not enough. Quantum computers provide an opportunity to speed up traditional algorithms for flow simulations. We show that lattice-based mesoscale numerical methods can be executed as efficient quantum algorithms due to their statistical features. This approach revises a quantum algorithm for lattice gas automata to reduce classical computations and state preparation at every time step. For this, the algorithm approximates the qubit relative phases and subtracts them at the end of each time step. Phases are evaluated using the iterative phase estimation algorithm and subtracted using single-qubit rotation phase gates. Further, this method optimizes the quantum resource required and makes it more appropriate for near-term quantum hardware. We also demonstrate how the checkerboard deficiency that the D1Q2 scheme presents can be resolved using the D1Q3 scheme. The algorithm is validated by simulating two canonical partial differential equations: the diffusion and Burgers' equations on different quantum simulators. We find good agreement between quantum simulations and classical solutions for the presented algorithm.

97 MATHEMATICS AND COMPUTING↗

A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method

Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.

Schwarz alternating method↗

Which way does the dendrite grow? Competition among epitaxy, preferred growth direction, and thermal gradients in powder bed fusion additive manufacturing

The as-processed microstructure of metal alloy parts manufactured through laser powder bed fusion (LPBF) is heavily derived from the cellular dendritic solidification. The growth direction of dendrites within the melt pool is determined through competition among epitaxial growth, preferred growth directions, and maximum thermal gradients. However, the dominant factor and the specific role of each in developing melt pool microstructures remain unknown. Here, in this study, we performed single laser track scans on an SS316L single crystal substrate and combined experimental characterization of microstructure and crystal orientations with Computational Fluid Dynamics simulations of thermal gradients to evaluate the role of each factor in determining dendritic growth direction and evolution. Our results reveal that epitaxial growth dominates microstructure development by preferentially growing along a single 〈100〉 variant of the single crystal substrate adjacent to the melt pool boundary. Under LPBF’s highly curved and rapidly evolving thermal field, this preferential dendrite variant selection and its continued growth from the melt pool boundary to the centerline are governed by the local temperature gradient magnitude at the solid-liquid interface, rather than by the instantaneous maximum temperature gradient direction alone. Using these findings, we successfully predict changes in the dendrite growth direction with changing laser scan direction on a single crystal substrate, and show that the geometric melt pool centerline can deviate from the microstructural centerline because asymmetric local temperature gradient magnitudes transiently limit growth, resulting in different dendrite travel distances on each side of the melt pool.

36 MATERIALS SCIENCE↗

Ultra-low thermal resistance and pressure drop copper and copper-tungsten diamond-shaped pin fin cold plates for liquid cooling of electronics

Modern and future data centers face increasing cooling challenges due to increasing chip thermal design power and die size, along with the need to reduce energy consumption used for cooling. High performance cooling solutions that maintain a low chip junction temperature are needed to ensure electronics reliability. This work develops an ultra-low thermal resistance and low pressure drop 75 mm × 75 mm cold plate, intended for next-generation electronics cooling. The cold plate features an array of diamond-shaped pin fins and integrated copper tungsten heat spreader, selected for its low coefficient of thermal expansion which reduces thermomechanical deformation and allows for closer integration of the cold plate with silicon dies. Starting with 300 candidate designs, three-dimensional computational fluid dynamics simulations predict the thermal-hydraulic performance of cold plate subsections. The highest performing geometries are evaluated with high fidelity simulations. Four cold plates are manufactured for experiments: three with diamond-shaped pin fins and one with straights fins for comparison purposes. The cold plates are fabricated from copper-tungsten (CuW), copper (Cu), or aluminum-silicon-magnesium alloy (AlSi10Mg). The diamond-shaped pin fins achieve a roughly 15 % lower thermal resistance compared to the conventional straight fin microchannel. The highest performing design achieves a chip-to-coolant (including thermal interface material) thermal resistance of 9.0 K/kW in CuW and 6.9 K/kW in Cu under a 1 kW heat load with an inlet-to-outlet pressure drop of 9.0 kPa and water as the working fluid. This work demonstrates ultra-low thermal resistance and pressure drop cold plates for large die, high heat load applications, and shows that CuW is an attractive cold plate material for improved reliability in next generation data center cooling.

Coefficient of thermal expansion↗

Metal additively manufactured wavy fin cold-plate architecture for improved thermal-hydraulic performance

Rapid growth in artificial intelligence and data center workloads demands high-performance liquid cooling to manage increasing chip power. This study presents two metal-additive-manufactured cold plates with sinusoidal fins, constant-amplitude wavy fins and linearly variable-amplitude wavy fins and compares them against metal-additive-manufactured straight fins using experiments conducted at 1 kW heat dissipation as well as high-fidelity 3D conjugate computational fluid dynamic simulations. The cold plates were printed in AlSi10Mg material and underwent design using a Python-automated workflow prior to manufacture and testing. The experiments show that wavy fins reduce the normalized thermal resistance by 35 to 45 % at water flow rates from 1 to 4 LPM. At a fixed 20 kPa pressure drop, the variable-waviness design lowered peak surface temperature by 9 °C and thermal resistance by 51 %, while edge-channel maldistribution in the constant wavy fin design limited gains. A thermal resistance breakdown revealed that 55–63 % of the total thermal resistance in wavy designs comes from base heat conduction, 27–33 % from fin heat conduction, and 9–13 % from fin heat convection, indicating the need to address conduction bottlenecks. Parametric sweeps identify a 3 mm fin pitch as optimal, and that horizontal inlet/outlet manifolds further reduce pressure drop by 30–60 % and thermal resistance by 9–16 % relative to vertical inlet-outlet manifolds. The results yield comprehensive guidelines for fin geometry, manifold alignment, material selection and additive-manufacturing constraints to realize high-performance liquid-cooled cold plates for power-dense electronics.

3d printing↗

A performant energy-conserving particle reweighting method for Particle-in-Cell simulations

A new particle-based reweighting method is developed and demonstrated in the Aleph Particle-in-Cell with Direct Simulation Monte Carlo (PIC-DSMC) program. Novel splitting and merging algorithms ensure that modified particles maintain physically consistent positions and velocities. This method allows a single reweighting simulation to efficiently model plasma evolution over orders of magnitude variation in density, while accurately preserving energy distribution functions (EDFs). Demonstrations on electrostatic sheath and collisional rate dynamics show that reweighting simulations achieve accuracy comparable to fixed weight simulations with substantial computational time savings. This highly performant reweighting method is recommended for modeling plasma applications that require accurate resolution of EDFs or exhibit significant density variations in time or space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Molten pool dynamics and humping suppression in high-speed laser welding via tailored beam configurations

High-speed laser welding is essential for increasing the production rate of fuel cell fabrication. However, when the welding speed exceeds a critical limit, humping occurs and reduces the weld quality. In this study, two tailored beam configurations, including an adjustable ring mode and a dual-beam configuration, were employed to suppress humping. Computational fluid dynamics simulations were performed to elucidate the underlying suppression mechanisms. Here, the results show that, in the adjustable ring mode, humping mitigation arises from a reduced backward cross-sectional melt flow rate and a more stable molten pool. In the dual-beam configuration, humping suppression is attributed to the deceleration of melt flow, the conduction-mode behavior of the trailing beam, and the widening of the molten pool induced by the trailing laser. Furthermore, because the dual-beam configuration directly modifies the trailing molten pool dynamics, it achieves more effective humping suppression, extending the welding speed limit to 1.50 m/s, compared with 1.00 m/s for the adjustable ring mode.

08 HYDROGEN↗

Temperature sensitivity of the equilibrium neutronics and accident analysis of the HTR-10

Pebble-Bed High-Temperature Gas-cooled Reactors (PB-HTGR) are moderated by the graphite in the fuel pebbles and the graphite reflector surrounding the pebble-bed. Because graphite is by far the most abundant material in PB-HTGRs and the primary moderator, accurate modeling of the graphite material, including density, impurities, and temperatures, is crucial for accurate computational modeling and simulation of these reactors. While main characteristics of the graphite components are often known, the local temperature is less well known and often averaged over all components. Here, this work studies the impact of considering accurate temperature profiles in the graphite material on the generation of a small PB-HTGR model at the state of equilibrium operation and on short-term accident progression. The fuel compositions for the PB-HTGR were determined using a jump-in equilibrium modeling method, the Axial Radial Zone Equilibrium Modeling (AR-ZEM) method. In contrast to previous work, the AR-ZEM method was used considering thermal-hydraulic feedback from the MELCOR code to determine temperatures of the fuel pebbles and the surrounding graphite reflector. The consideration of an axial and radial temperature profile in the core and reflector, as opposed to uniform material temperatures, had an impact of almost 1,300 pcm on the equilibrium core eigenvalue and caused significant differences in the discharged plutonium fuel inventory with up to 4.9% and 11.0% for Pu-239 and Pu-242, respectively. To assess the impact on short-term accident progression, two Anticipated Transient Without SCRAM (ATWS) events, a Pressurized Loss of Forced Coolant (PLOFC) and a Control Rod Withdrawal (CRW) with loss of flow, were simulated with MELCOR. The use of temperature profiles in the equilibrium core models did not reveal a significant impact on the temperature, power, or reactivity responses during the transients. In conclusion, a need for consideration of accurate temperature profiles, in particular for the graphite reflector, was found for the generation of equilibrium PB-HTGRs core models using jump-in methods, but detailed temperature profiles may not be necessary when performing conservative transient analysis.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Autonomous sensor suite for evaluating fish-turbine interactions and environmental impacts in marine renewable energy and hydropower

Marine renewable energy (MRE) harnesses ocean-based resources such as waves, tides, currents, and thermal or salinity gradients for sustainable power generation. It has the potential to complement existing renewable resources, support remote communities, and contribute to decarbonization efforts. However, understanding the hydrodynamic forces created by MRE devices and their impacts on marine life is critical for responsible deployment. Here, to address these concerns, advanced sensor devices, including the Marine Sensor Fish (MSF), Sensor Fish Mini (SF Mini), and Flexible Sensor Fish (FSF), were developed to measure interactions between aquatic organisms and MRE systems. This paper details the design, manufacturing, calibration, and field deployment of these sensor suites, highlighting their ability to capture key physical stressors such as shear forces, pressure changes, and collision impacts. The MSF successfully evaluated turbine interactions at a tidal turbine in the Salish Sea, capturing data on turbulence, collision impact, and pressure gradients. The SF Mini validated hydrodynamic conditions in scaled hydraulic models, supporting computational fluid dynamics simulations. The FSF, with its flexible silicone body, measured species-specific impacts in turbulent environments. This research demonstrates the potential of Sensor Fish technology to advance sustainable marine energy systems by reducing biological impacts and informing environmentally sustainable designs.

Ecological impacts↗

Steric Mapping, Ligand Dynamics, and Cycloisomerization Catalysis with Redox Robust Mn I/0/-I Dicarbenes

Manganese is perhaps the most electronically versatile element, yet the redox properties, reactivity, and catalytic applications of low-valent Mn 0 /Mn −I complexes remain underexplored due to the propensity for Mn 0 to dimerize, quenching highenergy metalloradicals. We report a series of redox-active monometallic Mn I , Mn 0 and Mn −I complexes containing a BH 2 - bridged dicarbene, characterized using a suite of experimental and cutting-edge computational (DFT) methods. Slow electron transfer kinetics at Mn I/0 are observed, with computations and electrochemical simulations in excellent agreement with experimental values. Despite the lack of steric bulk at the BH 2 -bridged Mn 0 , the t Bu groups at the dicarbene provide adequate steric protection to prevent dimerization, with percent buried volume (% V bur ) serving as a valuable steric ranking tool. We also show that a %V bur > 83% prevents dimerization for a diverse array of Mn 0 complexes from the literature. Ligand sterics of BPh 2 -and BH 2 -bridged complexes dictate reaction outcomes when Mn I and Mn −I are exposed to nucleophiles and electrophiles, respectively, while Mn0 facilitates the radical cycloisomerization catalysis of 6-iodo-1- hexene at room temperature. Furthermore, this work underscores the importance of ligand sterics in rationalizing reactivity patterns at Mn and provides valuable insights for designing chelating ligands that can selectively leverage Mn I/0/‑I states in redox-mediated catalytic reactions.

Ligands↗

Modeling High Current Pulsed Discharge in AA Battery Cathodes: The Effect of Localized Charging during Rest

During high current operation, substantial heterogeneity develops within battery cathodes, particularly when their thickness is large. Heterogeneity relaxation during subsequent rest is important for understanding battery performance under pulsed conditions. Localized charge balancing phenomena within batteries at zero net current are not well understood and merit investigation. In this work, the heterogeneity within cathodes of commercial alkaline Zn–MnO 2 batteries is measured during discharge and monitored during rest using energy dispersive X-ray diffraction (EDXRD). Significant gradients in protonation form during discharge and partially relax under rest. It is demonstrated that the proton gradient relaxation is through local redox activity at zero net current, where local (de)protonation works to redistribute charge across the cathode thickness. To support this redox-based relaxation, a fundamental kinetic study on prismatic MnO 2 cathodes is conducted to determine an appropriate model to describe both discharge and charge kinetics of MnO 2 . These kinetics are incorporated into a computational model to simulate the proton gradient formation and partial relaxation under identical discharge conditions as the operando EDXRD experiments. Model and experimental data are found to be in excellent agreement, correctly predicting localized charge balancing at rest.

Batteries↗

Reactivity of Carbonate Solvent Electrolytes on Lithium Silicon Anodes

Silicon (Si) is promising for lithium-ion battery (LIB) anodes due to their high theoretical capacity and low electrochemical potential. However, significant challenges remain, including severe volumetric expansion during cycling and the electrochemical instability of electrolytes, which leads to the formation of a nonuniform solid electrolyte interphase (SEI). To investigate SEI formation mechanisms, computational molecular dynamics simulations offer valuable insights. In this work, we examine the trajectories and charge transfer behavior of lithium hexafluorophosphate (LiPF 6 ) salt with various solvent compositions using density functional theory (DFT) and ab initio molecular dynamics (AIMD). Among the tested electrolyte systems, LiPF 6 with vinylene carbonate (VC) added to ethyl methyl carbonate (EMC) exhibits the lowest reactivity with the Si anode. In contrast, the effects of fluoroethylene carbonate (FEC) and VC depend on whether the primary solvent is EMC alone or a mixture of ethylene carbonate (EC) and EMC. Moreover, we show that electrolyte reactivity varies with the degree of lithiation of the Si anode (LiSi vs Li 15 Si 4 ) and under different charge states. To decouple electrolyte reactivity from surface effects, we analyze the dissociation and formation energies of individual species from solvated configurations. Overall, these first-principles-based findings provide a strategic foundation for electrolyte design to improve cycling stability and extend calendar life in LIBs using Si anodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cold-Sprayed NMC622 Composite as a Cathode for Lithium-Ion Batteries

The growing demand for high-energy, low-cost lithium-ion batteries (LIBs) to power electric vehicles (EVs) necessitates advances in both materials and manufacturing processes. Conventional cathode fabrication methods, such as slurry casting and drying, are energy-intensive and pose challenges for scalability and environmental compliance. In this study, we propose cold-spray (CS) deposition as a solvent-free approach for fabricating LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622) composite cathodes. Powder blends of NMC622, poly(vinylidene fluoride) (PVDF), and carbon black (CB) are directly deposited onto stainless steel and Inconel substrates under varied gas temperatures, pressures (and thus velocities), and standoff distances. The effects of temperature on the deposit morphology, coating density, and volume are systematically investigated. Computational fluid dynamics simulations reveal that increasing the gas temperature enhances the particle velocity, narrows the spray angle, and reduces the mass concentration radially at the nozzle outlet. CS deposition results in a dense cathode microstructure, accompanied by fracture in polycrystalline NMC622 particles. X-ray diffraction analysis further verifies that there are no phase changes during the deposition process. The electrochemical performance of the cold-sprayed cathodes reveals an initial capacity of approximately 96 mAh g –1 for single-crystal NMC622 and 167 mAh g –1 for polycrystalline NMC622. While these values are modest compared to state-of-the-art slurry-cast cathodes, which typically exhibit 180–200 mAh g –1 under optimized conditions. The results demonstrate a competitive performance given the solvent-free nature of the CS process and compare favorably with tape-cast samples made from identical feedstock. In conclusion, he CS process enables the formation of dense, binder-integrated cathode coatings without the need for solvent processing, offering a promising pathway for scalable, energy-efficient dry electrode manufacturing of next-generation LIBs.

Batteries↗

Unveiling mechanisms and onset threshold of humping in high-speed laser welding

The fabrication of fuel cells relies on a rapid laser welding process. However, challenges arise with the occurrence of humping when the welding speed surpasses a critical threshold, which poses difficulties in achieving a smooth surface finish and a consistent weld strength. This study aims to elucidate the humping mechanisms by analyzing the morphology of molten pool and the characteristics of melt flow at varying welding speeds via in situ synchrotron high-speed X-ray imaging and computational fluid dynamics simulations. Our findings indicate that the short keyhole rear wall, the high backward melt velocity, and the prolonged tail of molten pool are the primary factors contributing to the onset of humping. Furthermore, a dimensionless humping index (π h ) was introduced, which successfully captured the onset threshold of humping across different literatures. This index not only provides a quantitative description of the humping formation tendency but also serves as a valuable tool for optimizing the laser welding process.

42 ENGINEERING↗

Many-body entanglement in solid-state emitters

The preparation and control of quantum states lie at the heart of quantum information science. Recent advances in solid-state quantum emitters (QEs) and nanophotonics have transformed the landscape of quantum photonic technologies, enabling scalable generation of quantum states of light and matter. A new frontier in solid-state quantum photonics is the engineering of many-body interactions between QEs and photons to achieve robust coherence and controllable many-body entanglement. These entangled states, including photonic graph and cluster states, superradiant emission and emergent quantum phases, are promising for quantum computation, sensing and simulation. However, intrinsic inhomogeneities and decoherence in solid-state platforms pose considerable challenges in realizing such complex entangled states. This Review provides an overview of fundamental many-body interactions and dynamics at the light–matter interfaces of solid-state QEs and discusses recent advances in mitigating decoherence and harnessing robust many-body coherence.

36 MATERIALS SCIENCE↗

Explainable machine learning for incipient anomaly detection in compact molten salt heat exchanger with overlapping feature distributions

High-temperature molten salt-cooled reactors (MSCRs) are a promising next-generation nuclear technology option, offering efficient power conversion and inherent safety features. However, the reliability of these systems depends on the robust operation of heat exchangers (HXs), which are susceptible to failure due to temperature gradients and channel plugging caused by fluid freezing. Conventional monitoring methods, relying on inlet and outlet measurements, lack the spatial resolution needed to detect early-stage faults. We propose a novel design of a compact salt-to-salt matrix-type HX design consisting of interleaved arrays of parallel tubes, with integrated synthetic fiber optic distributed temperature sensing (DTS) to enable localized detection of incipient faults. To evaluate performance of this design, we generate high-fidelity synthetic data using heat transfer computational modeling to simulate channel plugging, and introduce sensor noise for realistic modeling of measurements. The dataset comprises of 97% normal operation and 3% anomaly cases, with each anomaly class representing 1% of the data. These early anomalies result in overlapping temperature profiles between normal and faulty channels, producing a non-separable dataset that challenges traditional classification techniques. We benchmark eight supervised machine learning (ML) models and demonstrate that XGBoost achieves the highest performance. To improve transparency, we develop an explainability framework combining Shapley values and partially ordered sets (POSETs) to quantify and structurally analyze feature importance. This approach identifies both dominant predictors and ambiguous feature relationships, enhancing trust and interpretability. Our results highlight the potential of combining DTS and explainable ML with intelligent feature selection to improve predictive maintenance and ensure operational resilience in advanced nuclear systems.

Prantikos, Konstantinos [Argonne National Laborato↗