Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational physics”

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.

631 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Numerical simulation of frost formation and heat transfer on fin-and-tube heat exchangers in turbulent cross-flow

Frost formation in fin-and-tube heat exchangers in turbulent cross-flow presents significant challenges in industrial refrigeration applications, affecting heat transfer efficiency and operational reliability. The purpose of this work is to investigate frost deposition and growth on a staggered bank of a fin-and-tube freezer coil under turbulent forced convection conditions. The focus here is on investigating conditions that closely replicate real-world scenarios in large walk-in industrial freezers. Using a direct numerical simulation approach, we examine the flow dynamics and thermal behaviour in the presence of frost, considering turbulent regimes characterized by a Reynolds number in the range 1050 ≤ R e D , avg ≤ 4800 , with the characteristic length being the outer diameter of the tube and the velocity being the bulk fluid velocity between the plates (fins). Computational fluid dynamics simulations are employed to resolve the interactions between turbulent airflow and the frost layer. Our approach incorporates a modified immersed boundary method and a slow-time acceleration technique to address the complex dynamic interface between the continuously evolving frost layer and the flowing air stream. Our findings indicate that frost forms more on the sides of the finned surfaces (plates) and less on the tubes themselves. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗
DOE OSTI2025

Customizable wave tailoring nonlinear materials enabled by bilevel inverse design

Abstract Passive wave transformation via nonlinearity is ubiquitous in settings from acoustics to optics and electromagnetics. It is well known that different nonlinearities yield different effects on propagating signals, which raises the question of “what precise nonlinearity is the best for a given wave tailoring application?” In this work, considering a one-dimensional spring-mass chain connected by polynomial springs (a variant of the Fermi-Pasta-Ulam-Tsingou system), we introduce a bilevel inverse design method which couples the shape optimization of structures for tailored constitutive responses with reduced-order nonlinear dynamical inverse design. We apply it to two qualitatively distinct problems—minimization of peak transmitted kinetic energy from impact, and pulse shape transformation—demonstrating our method’s breadth of applicability. For the impact problem, we obtain two fundamental insights. First, small differences in nonlinearity can drastically change the dynamic response of the system, from severely under- to outperforming a comparative linear system. Second, the oft-used strategy of impact mitigation via “energy locking” bistability can be significantly outperformed by our optimal nonlinearity. We validate this case with impact experiments and find excellent agreement. This study establishes a framework for broader passive nonlinear mechanical wave tailoring material design, with applications to computing, signal processing, shock mitigation, and autonomous materials.

Science & Technology - Other Topics↗
DOE OSTI2025

Cyclic Peptides for Lanthanide Binding

Lanthanide ions are difficult to separate from one another due to their similar chemical properties. The discovery of lanthanide-binding peptides and proteins in nature has led to an increased interest in the possibility of utilizing the strong binding of peptides to lanthanide ions for their separations; as such, there has been an effort to identify or design peptides with improved lanthanide binding and selectivity toward particular lanthanide ions. Here, in this study, we designed and characterized lanthanide-binding cyclic peptides (LBCPs) with molecular dynamics simulations, electronic structure calculations, and emission spectroscopy. Luminescent decay measurements were done to determine the number of water molecules coordinated to the Eu 3+ ion in Eu-LBCP complexes and compare to the predicted number of water molecules by computation to assess the lanthanide-binding affinity of LBCPs. Measured stability constants show binding of the LBCPs to the Eu 3+ ion with stronger than micromolar affinity. We were able to identify multiple peptides that selectively bind to middle lanthanides. We describe the structural basis of the lanthanide-binding selectivity trend with strongest binding to the middle lanthanides, followed by the heavier lanthanides, and finally to the lighter ions.

ions↗
DOE OSTI2026

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE↗
DOE OSTI2026

Multitarget Rydberg gates via spatial blockade engineering

Multi-target gates offer the potential to reduce gate depth in syndrome extraction for quantum error correction. Although neutral-atom quantum computers have demonstrated native multi-qubit gates, existing approaches that avoid additional control or multiple atomic species have been limited to single-target gates. We propose single-control-multi-target CZ^n gates on a single-species neutral-atom platform that require no extra control and have gate durations comparable to standard CZ gates. Our approach leverages tailored interatomic distances to create an asymmetric blockade between the control and target atoms. Using a GPU-accelerated pulse synthesis protocol, we design smooth control pulses for CZZ and CZZZ gates, achieving fidelities of up to 99.55% and $99.24\%$, respectively, even in the presence of simulated atom placement errors and Rydberg-state decay. Our approach is most effective for N=2 (CZZ) and N=3 targets (CZZZ); for larger N, increasing spatial crowding of the targets introduces significant challenges for maintaining the required blockade asymmetry. This work presents a practical path to implementing low-overhead multi-target gates in single-species neutral-atom systems, significantly reducing the resource overhead for syndrome extraction. To motivate the impact of these gates, we apply a greedy scheduling algorithm and we demonstrate that our proposed gates can reduce the number of atom reconfiguration costs by up to 50% for color code syndrome extraction of code distances greater than 5.

Stein, Samuel A.↗
DOE OSTI2026

A transport-based framework for solidification cracking in Ni-based superalloys processed by laser powder bed fusion

Solidification cracking remains a persistent barrier to laser powder bed fusion (LPBF) processing of solid-solution-strengthened (SSS) Ni-based superalloys and is commonly attributed to carbide formation, liquation-type failure, or elemental segregation. In this work, the cracking behavior of three SSS Ni-based superalloys—Inconel 625 (625), Inconel 617 (617), and Haynes 230 (230)—is examined under comparable LPBF conditions to evaluate the origins of their cracking susceptibility. Carbides were present in both 625 and 230, yet cracking behavior differed significantly. MC-type carbides in 625 remain discrete and preserve liquid connectivity, whereas Cr-rich M23C6 carbides in 230 form at cellular triple junctions that bottleneck the interdendritic liquid network, influencing cracking through liquid connectivity rather than as an independent cause. Alloy 617 exhibited anomalous segregation without resolvable secondary phases yet cracked mildly, indicating that extensive carbide formation is not a prerequisite for cracking. To rationalize these observations, a transport-limited framework is proposed in which cracking occurs when terminal-liquid feeding cannot accommodate solidification-induced strain; among the mechanisms evaluated, this framework is most consistent with the physical constraints imposed by LPBF solidification. Transport-based crack-susceptibility indexes incorporating permeability, viscosity, and solidification kinetics distinguish the alloy hierarchy; permeability ratios Krat increased from approximately 7 (625) to 11.5 (617) and 14.5 (230), with corresponding increases in the μ-weighted mushy zone risk index (CSIMZRM). These results are consistent with transport-limited liquid feeding as a rate-limiting contributor to solidification cracking in LPBF-processed SSS Ni-based superalloys, providing a physically grounded basis for alloy and process design.

Hyer, Holden [ORNL] (ORCID:0000000343915561)↗
DOE OSTI2026

Superconducting phase diagram of multilayer square-planar nickelates

The discovery of superconductivity in square-planar nickelates has offered a rich materials platform to explore the origins of high-temperature superconductivity. However, experimental investigations have largely been limited to the infinite-layer R NiO 2 ( R , rare earth) nickelates. For this work we constructed a phase diagram of multilayer square-planar Nd n+1 Ni n O 2n+2 compounds and found signatures of superconductivity for dimensionality n = 4 to 8. Upon decreasing n , the superconducting anisotropy evolves owing to 4ƒ electron effects, and electronic structure characteristics approach cuprate-like behavior. Magnetic fluctuations persist from within the superconducting regime and into the overdoped, nonsuperconducting regime. The superconducting regime overlaps with that of chemically doped infinite-layer nickelates, demonstrating underlying commonalities as well as differences across varying structural realizations of square-planar nickelates. Our work establishes this layered template for creating new nickel-based superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Study of wrap mode impact on Pseudomonas aeruginosa motion in the chemotactic field of fungi

An agent-based discrete computational model biologically calibrated to Pseudomonas aeruginosa migration is used to explore the impacts of bacterial reversals and wrap mode on the efficiency of motion in different environments, both with and without chemotaxis. It is first shown that wrap mode increases the exploration of continuous multimodal chemotactic profiles such as those produced by biologically relevant fungal networks. For cells undergoing a run-reverse pattern, it is shown that the bacteria are likely to remain at the first local chemoattractant maximal production site on a hypha they find. However, with wrap mode, the bacteria can more easily escape these local sites to further explore their neighboring environment along the fungi, suggesting that wrap mode may be beneficial for migration along the fungi in liquid. In a different set of simulations of bacterial motion close to an isolated chemotactic source, wrap mode is shown to increase the ability of a bacterium to reorient toward the source while reducing the overall motion required for similar chemotactic efficiency as a run-reverse strategy, suggesting a potential metabolic benefit. In contrast, model simulations show that wrap mode can increase the rate of dispersal of P. aeruginosa in a nonchemotactic environment.

Hansen, Austin [University of California, Riversid↗
DOE OSTI1969

Application of Nuclear Technology to Art Identification Problems: First Annual Report

Throughout the centuries, as man has become increasingly affluent, his interest in antiquities and objects of art, and the price he has been willing to pay for them, has increased. Inevitably, one result has been to attract the unscrupulous to the production and sale of counterfeit paintings, sculptures, and other works of art. The skill of forgers varies but often is highly sophisticated. However, dealers, museum curators, and other experts have also become more skillful in the detection of forgeries. Scientific tools of increasing sensitivity and sophistication have gradually been applied to the examination of the materials of art and archaeology. Such tools frequently support and render more reliable the judgment of the experts. The quality of forgeries may improve further as their makers become acquainted with the new methods of detection and in turn learn to circumvent or confound the methods. Therefore, the development of still more advanced methods of examination has great utility in the art world. This is particularly true if the ultimate effect is to make circumvention of the methods of examination so costly that the economic incentives for producing forgeries may be substantially reduced, if not eliminated.

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

Precise Modeling of a Complex Solenoidal Magnetic Field Using a Combination of Analytic Functions and a PINN

We demonstrate an iterative approach to modeling a sparsely measured magnetic field in a large-bore solenoid. This approach uses a hybrid of traditional and machine learning techniques. The traditional technique is a linear least-squares fit using a series solution to Laplace's equation, while the machine learning technique involves the training of a physics-informed neural network (PINN) on the least-squares fit residuals. We use a newly defined activation function "DELTAsnake," a modification to the snake activation function proposed by Ziyin et al. that allows for stronger curvature and non-monotonicity. The combined model approximately obeys Maxwell's equations to a level sufficient for producing high quality physics simulations and analysis. Our approach is applied to a highly realistic calculation of the expected magnetic field in the Mu2e experiment's Detector Solenoid which includes a simple model for the expected statistical measurement uncertainties. Using ten toy measurement simulations, we demonstrate the capabilities of our model in comparison to the least-squares method alone; the least-squares method alone results in a reduced chi-squared statistic of ${2.15 \pm 0.01}$, while our approach improves the reduced chi-square to ${1.034 \pm 0.005}$. Furthermore, for an average toy simulation, we show that the range of the RMS of the three field component residuals reduces from ${0.07-0.37}$ Gauss to ${0.05-0.07}$ Gauss. We find that this novel method is robust against a realistic systematic uncertainty deriving from Hall probe calibration bias and can be used to significantly reduce the number of measurements required to achieve an accurate model.

Kampa, Cole [Caltech] (ORCID:0000000192972920)↗
DOE OSTI2026

Imaging of a van der Waals spin-orbit torque system using spin ensembles in hBN

Recently, optically active spin defects embedded in two-dimensional (2D) van der Waals (vdW) crystals have emerged as a transformative quantum sensing platform to explore cutting-edge materials science. Taking advantage of excellent solid-state integrability, this new class of spin defects can be readily arranged in nanoscale proximity to target materials, showing great promise for realizing in-situ quantum sensing of microscopic spin and charge behaviors in vdW heterostructures. Here we report hexagonal boron nitride-based quantum imaging of field-free deterministic magnetic switching and electric current distributions in an all-vdW spin-orbit torque (SOT) system. By visualizing variations of nanoscale magnetic stray field profile of room-temperature 2D magnet Fe 3 GaTe 2 under different SOT conditions, we show how the magnetic switching evolves from deterministic to stochastic behavior due to the interplay between spin orientations, anisotropy and Joule heating. Micromagnetic simulations rationalize our results well, revealing the role of field-like SOT in inhibiting thermal fluctuation driven stochastic switching and chaotic multi-domain competition. This understanding, which is otherwise difficult to access by conventional transport measurements, offers valuable insights into material design, testing, and performance evaluation of next-generation vdW spintronic devices.

Imaging techniques↗
DOE OSTI2026

Large Language Model for Validation, Optical Calibration, and Learning (VOCAL) Distributed Temperature Sensing Interface

Distributed temperature sensing (DTS) using fiber optic sensors (FOS) offers a promising method for temperature measurements in advanced reactors, such as sodium fast reactors and molten salt cooled reactors. To support the calibration and validation of DTS measurements, Argonne National Laboratory developed the Validation, Optical Calibration, and Learning (VOCAL) software package. This report describes the integration of a local large language model (LLM) with a retrieval-augmented generation (RAG) system into the VOCAL interface to serve as an interactive user assistant. The LLM framework enhances the VOCAL platform’s accessibility to users by explaining interface components, clarifying inputs and outputs, and answering user queries dynamically in real-time. The accuracy of the LLM assistant performance was evaluated with 20 queries regarding the interface and its parameters using experimental data from the Thermal Hydraulic Experimental Test Article (THETA) facility. Results demonstrate that the LLM achieved a 95% accuracy rate, with a BERTScore of 0.8816 and SBERT value of 0.7417. Furthermore, validation of the RAG system within the LLM framework showed optimal accuracy with k-values between 1 and 2 using the k-refinement convergence test. The prompt perturbation analysis demonstrated good initial consistency for the RAG system, exhibiting the highest accuracy under punctuation variations and the greatest sensitivity under query reordering. Notably, the model’s errors were limited to data retrieval failures rather than factual hallucinations, reinforcing its baseline reliability. The integration of LLM provides a highly accurate, userfriendly enhancement to the VOCAL platform without disrupting its core computational capabilities for FOS calibration and validation.

Hong, Evan↗
DOE OSTI2026

Biochar–polymer composites for 3D printing: a review

Biochar, a bio-based co-product of biofuel production via thermochemical conversion, holds potential as a filler for polymer composites to reduce costs, improve thermomechanical properties, and aid in environmental remediation. 3D-printed biochar composites have received growing interest over the past few years but have experienced difficulties such as poor layer adhesion and nozzle clogging. Currently, no literature review examines 3D-printed biochar composites and related biochar properties in-depth. This work summarizes and discusses recent studies on 3D-printed polymer and biochar composites and examines their mechanical, thermal, and additional properties that result from each study. Technical challenges in printability, such as nozzle clogging from particle size and biochar aggregation, are also discussed. Furthermore, this work discusses the variability of biochar properties resulting from the pyrolysis conditions and feedstock choice in relation to potential 3D printing outcomes. In particular, several studies reported that high lignin feedstocks could be candidates for 3D printing. The post-processing approaches of the biochar via physical and chemical methods are also introduced. Ball milling appears to hold the most promise for physical treatments due to its tunability of particle size, surface area, and functional groups, while chemical treatments with acids or alkalis are used to tailor biochar porosity and wettability. Overall, it was determined that future research needs to be done relating biochar production and post-processing methods to resulting 3D printing parameters as the number of studies is limited.

Day, Rachel [Auburn University, Auburn, AL]↗
DOE OSTI2026

Unraveling interphase-driven failure pathways in LiMn0.6Fe0.4PO4/graphite pouch cells

LiMnxFe1−xPO4 (LMFP) is a promising high-voltage, thermally stable, and earth-abundant cathode material, yet its practical application is limited by interphase instability and Mn dissolution. In this work, we systematically evaluate LiMn0.6Fe0.4PO4/graphite pouch cells using three electrolyte formulations including control carbonate electrolyte, control + 2 wt% vinylene carbonate (VC), and control + 2 wt% VC + 1 wt% 1,3,2-dioxathiolane 2,2-dioxide (DTD), to establish how electrolyte composition governs interphase chemistry and long-term degradation. Electrochemical testing shows that both additives are preferentially reduced prior to ethylene carbonate (EC) during cell formation, generating robust cathode-electrolyte interphase (CEI) and solid-electrolyte interphase (SEI) layers that suppress gas evolution and raise the first-cycle coulombic efficiency to 89.3%. Additionally, the dual-additive electrolyte delivers the most stable performance, retaining over 85% capacity after 600 cycles while minimizing impedance growth under long-term cycling at C/3 and 40 °C. Soft X-ray absorption spectroscopy confirms that VC + DTD effectively suppresses electrolyte oxidation at the cathode surface, and micro-X-ray fluorescence shows substantially reduced Mn dissolution and deposition on the graphite anode. Density functional theory simulations further provided insights into the structural and energetic influences of alkoxide species on the cathode surface, proposing a Mn2+ extraction mechanism. The combined experimental and computational findings establish a mechanistic link between electrolyte composition and interphase evolution, highlighting the effectiveness of electrolyte engineering for extending the operational lifetime of LMFP-based lithium-ion batteries.

Chak, Chanmonirath Michael↗
DOE OSTI2026

Preliminary Assessment for the Electric Load Shifting Potential of Integrating Thermal Energy Storage with Heat Pumps in Residential Buildings in Texas of United States

The widespread adoption of electric-driven heat pumps for heating and cooling is expected to significantly increase electric demand. Cooling electric demand will result in an increase in peak hours electric loads, placing additional strain on the grid during peak hours. This challenge, combined with the current rapidly increasing demand from data centers, exacerbates the electric demand duck curve problem. Integrating thermal energy storage (TES) with heat pump can shift electric use for heating and cooling from peak to off-peak hours of the electric grid, which can help flatten the daily electric demand profile of a building. This paper presents a novel design that integrates heat pump with TES (HP-TES), which uses phase change materials (PCM). Heat pump charges TES by melting or freezing PCM during off-peak hours when there is no thermal demand from the building. TES is then discharged (i.e., by freezing or melting PCM) during peak hours to provide a more favorable heat source or heat sink for the heat pump to meet the thermal demand of the building with lower electricity use than conventional air-source heat pumps. Computer simulations were developed to predict the performance of HP-TES applied to a typical single-family house in the US. The building-level simulation results were scaled up to preliminarily assess the aggregated impacts of deploying HP-TES across all single-family houses in Texas of the United States, including reduction of peak demand of the electric grid.

Anees, Fady [ORNL]↗
DOE OSTI2026

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗
DOE OSTI2025

Investigating the Combustion Performance of Dual Fuel Combustion with Diesel and Port Injected Hydrogen in a Large Bore Locomotive Engine

The heavy-duty transportation sector has primarily relied on conventional diesel combustion engines given their reliability and high thermal efficiency relative to spark ignition engines, but increased focus on reducing greenhouse gas emissions has led to investigation into alternative fuels. Gaseous hydrogen fuel has garnered a great deal of recent interest in the engine community given it has zero carbon, but hydrogen is not available at the scale and cost that petroleum fuels are currently available, and this is a barrier to adoption for industries that are looking to decarbonize their operations. Because of the fuel flexibility provided, dual fuel technology offers a pathway for some industries to adopt hydrogen as a fuel source while maintaining sufficient flexibility in times and locations where the new fuel is not yet available. This computational study investigates dual fuel combustion in a large bore locomotive engine architecture using direct injected diesel and port injected gaseous hydrogen fuel. With an optimal port fuel injection configuration from previous work, simulations of varying substitution ratio, compression ratio, manifold air temperature, diesel injection timing, and diesel injection pressure were performed to understand their effect on combustion performance. Results indicated that both increased substitution ratio and higher intake air temperature accelerates hydrogen flame propagation and can result in high peak cylinder pressures. Additionally, diesel injection timing and injection pressure were demonstrated as effective methods for controlling dual fuel combustion heat release rates.

ODonnell, Patrick Christopher↗
DOE OSTI2026

Purcell enhancement of directional edge photocurrent in a van der Waals self-cavity

Cavities provide a means to manipulate the optical and electronic responses of quantum materials by selectively enhancing light-matter interaction at specific frequencies and momenta. While cavities typically involve external structures, exfoliated flakes of van der Waals (vdW) materials can form intrinsic self-cavities due to their small finite dimensions, confining electromagnetic fields into plasmonic cavity modes, characterized by standing-wave current distributions. While cavity-enhanced phenomena are well-studied at optical frequencies, the impact of self-cavities on nonlinear electronic responses—such as directional photocurrent—remains largely unexplored, particularly in the terahertz regime, critical for emerging ultrafast optoelectronic technologies. Here, we report a self-cavity-induced Purcell enhancement of directional photocurrents in the vdW semimetal WTe 2 . Using ultrafast optoelectronic circuitry, we measured coherent near-field THz emission resulting from nonlinear photocurrents excited at the sample edges. We observed enhanced emission at finite frequencies, tunable via excitation fluence and sample geometry, which we attribute to plasmonic interference effects controlled by the cavity boundaries. We developed an analytical theory that captures the cavity resonance conditions and spectral response across multiple devices. Our findings establish WTe 2 as a bias-free, geometry-tunable THz emitter and demonstrate the potential of self-cavity engineering for controlling nonlinear, nonequilibrium dynamics in quantum materials.

condensed-matter physics↗