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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

877 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2024

First Principles Study of Aluminum Doped Polycrystalline Silicon as a Potential Anode Candidate in Li‐ion Batteries

Addressing sustainable energy storage remains crucial for transitioning to renewable sources. While Li‐ion batteries have made significant contributions, enhancing their capacity through alternative materials remains a key challenge. Micro‐sized silicon is a promising anode material due to its tenfold higher theoretical capacity compared to conventional graphite. However, its substantial volumetric expansion during cycling impedes practical application due to mechanical failure and rapid capacity fading. A novel approach is proposed to mitigate this issue by incorporating trace amounts of aluminum into the micro‐sized silicon electrode using ball milling. Density functional theory (DFT) is employed to establish a theoretical framework elucidating how grain boundary sliding, a key mechanism involved in preventing mechanical failure is facilitated by the presence of trace aluminum at grain boundaries. This, in turn, reduces stress accumulation within the material, reducing the likelihood of failure. To validate the theoretical predictions, capacity retention experiments are conducted on undoped and Al‐doped micro‐sized silicon samples. In conclusion, the results demonstrate significantly reduced capacity fading in the doped sample, corroborating the theoretical framework and showcasing the potential of aluminum doping for improved Li‐ion battery performance.

25 ENERGY STORAGE↗

DOE OSTI2026

The effect of dogbone sample size in the tensile testing of TATB-based plastic‑bonded explosive materials

Mechanical properties are of interest for many plastic‑bonded explosive (PBX) materials with tensile properties being of particular interest. Direct tensile measurements using dogbone-shaped samples are considered the gold standard, but they are fairly large, making testing more costly and less desirable from a safety perspective. We investigated whether the measured tensile strength depends on the dogbone specimen size, which to our knowledge, has not been reported in the literature for PBX materials. Understanding this should inform the feasibility of employing smaller samples and how sample size should be considered when comparing PBX dogbone values in the literature. The TATB-based PBX dogbone sample size was varied by (a) scaling all dimensions proportionally and (b) varying only the length of the samples. It was observed that the measured tensile peak stress (strength) was a function of the sample size, and was more dependent on the diameter (cross-sectional area) than the length of the samples. Since peak stress is calculated as peak force normalized to the diameter of the sample, one might not expect an explicit diameter dependence for the peak stress. Therefore, these results suggest there may be an additional strengthening effect as the sample diameter is increased.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Multiscale and multidimensional modeling of particle acceleration and transport in solar flares

Multi-messenger, multi-viewpoint, and time-resolved observations of solar flares are now providing unprecedented constraints on particle acceleration sites, energy conversion, and energy transport. The interpretation of current observations, including microwave imaging spectroscopy from EOVSA, hard x-ray (HXR) imaging from Solar Orbiter/STIX, gamma-ray diagnostics from Fermi, and in situ measurements from Parker Solar Probe and Solar Orbiter, collectively demands modeling frameworks that go beyond traditional spatially unresolved, one-zone models or single-mechanism descriptions. This review surveys multiscale and multidimensional modeling approaches, including kinetic, magnetohydrodynamic (MHD), and macroscopic particle models, that are being developed to meet the need. Kinetic simulations reveal that three-dimensional (3D) effects, including field-line chaos and self-generated turbulence, are essential for sustained power-law particle acceleration. MHD simulations now capture flux-rope eruptions, plasmoid-unstable current sheets, and turbulent flare regions in realistic magnetic topologies. Macroscopic models coupling MHD with energetic-particle models produce spatially resolved electron distributions and synthetic HXR and microwave emissions for direct comparison with observations. Despite these advances, outstanding challenges remain in bridging kinetic and global scales, improving MHD simulations and macroscopic particle models, and achieving quantitative model-observation closure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Low-energy calibration of SuperCDMS HVeV cryogenic silicon calorimeters using Compton steps

Cryogenic calorimeters for low-mass dark matter searches have achieved sub-eV energy resolutions, driving advances in both low-energy calibration techniques and our understanding of detector physics. The energy deposition spectrum of gamma rays scattering off target materials exhibits step-like features, known as Compton steps, near the binding energies of atomic electrons. Here, we demonstrate a successful use of Compton steps for sub-keV calibration of cryogenic silicon calorimeters, utilizing four SuperCDMS High-Voltage eV-resolution detectors operated with 0 V bias across the crystal. This new calibration at 0 V is compared with the established high-voltage calibration using optical photons. The comparison indicates that the detector response at 0 V is about 30% weaker than expected, highlighting challenges in detector response modeling for low-mass dark matter searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Stabilizer Scars

Quantum many-body scars are eigenstates in nonintegrable isolated quantum systems that defy typical thermalization paradigms, violating the eigenstate thermalization hypothesis and quantum ergodicity. Here, we identify exact analytic scar solutions in a 2+1 dimensional lattice gauge theory in a quasi-1D limit as zero-magic resource stabilizer states. Our results also highlight the importance of magic resources for gauge theory thermalization, revealing a connection between computational complexity and quantum ergodicity.

eigenstate thermalization↗

NASA NTRS2022

Optical properties of morphologically complex black carbon aerosols: Effects of coatings

Optical properties are computed for fractal-like aggregate black carbon (BC) aerosols coated with different substances. Two models are used for these aerosols: (i) the coated aggregate model (model I), where coating is added to the voids and surroundings of BC fractal-like aggregates; and (ii) the closed-cell aggregate model (model II), where coating is added concentrically to each monomer of BC fractal-like aggregates. Our results favor choosing the coated aggregate model (model I) to simulate scattering and absorption by coated BC aerosols because this model is morphologically more realistic, and because this model yields mass absorption cross section (MAC) and backscattering linear depolarization ratio (LDR) values that are consistent with field measurements. Moreover the corresponding computed degree of linear polarization (DoLP) and LDR values are very sensitive to changes in the coating volume fraction (fvol) and the coating refractive index (m). With the same absorbing BC core, the MAC value increases steadily with increasing fvol and m values. For example, using a density of 1.8 g/cm 3 and BC refractive index 1.95 + i0.79, the calculated MAC values for uncoated BC aggregates range from 6.2 to 6.8 m 2 /g at 0.55 µm. When coating material is applied to the BC aggregates, the calculated MAC values for model I particles increase to between 9.8 and 13.2 m 2 /g (depending on m) when fvol = 87.5%. The backscattering LDR values also tend to increase with the increasing m values for the shapes and sizes considered in our study. For model I particles, the backscattering LDR values span a wide range of 4.2–27.8% at a wavelength of 0.35 µm at fvol = 87.5% when m increases from 1.33 to 1.55. Our results are relevant to analyses of polarimetric and lidar observations of smoke particles, especially when these particles undergo hygroscopic growth.

Black carbon aerosols↗

DOE OSTI2026

The impact of aerosol mixing state on immersion freezing: insights from classical nucleation theory and particle-resolved simulations

Immersion freezing, initiated by ice-nucleating particles (INPs) in supercooled aqueous droplets, plays an important role in the formation of ice crystals within clouds. The efficiency of immersion freezing depends strongly on INP composition and, crucially, on the mixing state – how chemical species are distributed across the particle population. Here, we quantify the impact of aerosol mixing state on immersion freezing using a combined theoretical and particle-resolved modeling approach. We derive analytical expressions for the frozen fraction of internally and externally mixed INP populations based on classical nucleation theory, showing that the frozen fraction is sensitive to whether ice-active species are present in all particles or only in a subset of the population. We introduce a multi-species immersion freezing scheme into the particle-resolved model PartMC, using the water activity-based immersion freezing model (ABIFM) to compute freezing probabilities for mixed-composition particles. To improve computational efficiency, we implement a Binned Tau-Leaping algorithm and demonstrate an order-of-magnitude speedup with minimal accuracy loss. Simulations reproduce the analytical trends in limiting cases and extend the analysis to more general aerosol populations, where mixing state continues to exert a substantial control on frozen fraction. Sensitivity analyses across particle size, species type, and cooling condition reveal that the mixing state effect is most pronounced when small amounts of highly efficient INPs are mixed with less efficient materials. These findings underscore the need to represent aerosol mixing state explicitly in models of heterogeneous ice nucleation to reduce uncertainty in cloud-phase partitioning.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Alumina–Titania Nanolaminate Condensers for Hot Programmable Catalysis

Nanolaminates composed of thin alternating layers of Al2O 3 and TiO 2 (ATO) were engineered by using atomic layer deposition as the dielectric material for a Pt-on-carbon catalytic condenser. Investigation assessed synthesis parameters including the deposition temperature, Al 2 O 3 and TiO 2 layer thicknesses, total number of layers, and a capping Al 2 O 3 layer on the maximum charge accumulation in the Pt catalyst. The highest capacitance ATO configuration demonstrated a specific capacitance of ∼1200 nF/cm 2 with working voltages of ±5 V, enabling the storage of 4 × 10 13 electrons or holes per cm 2 at room temperature. The ATO devices exhibited enhanced capacitance at elevated temperatures of up to 400 °C, suggesting the suitability of these materials for high-temperature applications. Adsorption of carbon monoxide on the Pt/C-ATO device characterized by grazing incidence infrared spectroscopy showed changes in the surface binding energy of 13.1 ± 0.8 kJ/mol for an applied external voltage bias of ±1 V.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Ramp-release experiments for strength measurements: Strain-rate dependence

This paper presents an enhanced analysis method for investigating material properties at high strain rates, extending the capability of established experimental techniques to gain more information. The ramp-release method has been applied to many experiments reported at high (≈10 5 − 10 6 s −1 ) strain-rates. More recent data gathered at the National Ignition Facility (NIF) has enabled higher (≈10 8 s −1 ) strain-rates to be studied. Here, we present an initial application of ramp-release analysis to NIF ramp-compression data, illustrating both the opportunities and the practical challenges of extending these methods to laser-driven platforms. The higher strain-rates accessed at the NIF mean that there is more strain-rate enhancement to strength, and the experimental configuration means that this enhancement is more readily seen in the data. This is enabled by the capability of avoiding peak-compression attenuation through the sample thickness with a designed hold period made possible by the pulse-shaping capability of NIF. We propose that this combination of experimental conditions and an enhanced analysis method enables the strain-rate enhancement to strength to be studied, and potentially for this to inform physics models at smaller scales than the continuum.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Hybrid Data‐Driven Discovery of High‐Performance Silver Selenide‐Based Thermoelectric Composites

Optimizing material compositions often enhances thermoelectric performances. However, the large selection of possible base elements and dopants results in a vast composition design space that is too large to systematically search using solely domain knowledge. To address this challenge, a hybrid data-driven strategy that integrates Bayesian optimization (BO) and Gaussian process regression (GPR) is proposed to optimize the composition of five elements (Ag, Se, S, Cu, and Te) in AgSe-based thermoelectric materials. Data is collected from the literature to provide prior knowledge for the initial GPR model, which is updated by actively collected experimental data during the iteration between BO and experiments. Within seven iterations, the optimized AgSe-based materials prepared using a simple high-throughput ink mixing and blade coating method deliver a high power factor of 2100 µW m −1 K −2 , which is a 75% improvement from the baseline composite (nominal composition of Ag 2 Se 1 ). In conclusion, the success of this study provides opportunities to generalize the demonstrated active machine learning technique to accelerate the development and optimization of a wide range of material systems with reduced experimental trials.

36 MATERIALS SCIENCE↗

NASA NTRS1985

Nickel-cadmium Battery Cell Reversal from Resistive Network Effects

During the individual cell short-down procedures often used for storing or reconditioning nickel-cadmium (Ni-Cd) batteries, it is possible for significant reversal of the lowest capacity cells to occur. The reversal is caused by the finite resistance of the common current-carrying leads in the resistive network that is generally used during short-down. A model is developed to evaluate the extent of such a reversal in any specific battery, and the model is verified by means of data from the short-down of a f-cell, 3.5-Ah battery. Computer simulations of short-down on a variety of battery configurations indicate the desirability of controlling capacity imbalances arising from cell configuration and battery management, limiting variability in the short-down resistors, minimizing lead resistances, and optimizing lead configurations.

Zimmerman, A. H.↗

DOE OSTI2024

Sustainable Production of Biomass‐Derived Graphite and Graphene Conductive Inks from Biochar

Abstract Graphite is a commonly used raw material across many industries and the demand for high‐quality graphite has been increasing in recent years, especially as a primary component for lithium‐ion batteries. However, graphite production is currently limited by production shortages, uneven geographical distribution, and significant environmental impacts incurred from conventional processing. Here, an efficient method of synthesizing biomass‐derived graphite from biochar is presented as a sustainable alternative to natural and synthetic graphite. The resulting bio‐graphite equals or exceeds quantitative quality metrics of spheroidized natural graphite, achieving a RamanI D /I G ratio of 0.051 and crystallite size parallel to the graphene layers (L a ) of 2.08 µm. This bio‐graphite is directly applied as a raw input to liquid‐phase exfoliation of graphene for the scalable production of conductive inks. The spin‐coated films from the bio‐graphene ink exhibit the highest conductivity among all biomass‐derived graphene or carbon materials, reaching 3.58 ± 0.16 × 10 4 S m −1 . Life cycle assessment demonstrates that this bio‐graphite requires less fossil fuel and produces reduced greenhouse gas emissions compared to incumbent methods for natural, synthesized, and other bio‐derived graphitic materials. This work thus offers a sustainable, locally adaptable solution for producing state‐of‐the‐art graphite that is suitable for bio‐graphene and other high‐value products.

Chemistry↗

DOE OSTI2024

Trace benzene capture by decoration of structural defects in metal–organic framework materials

Abstract Capture of trace benzene is an important and challenging task. Metal–organic framework materials are promising sorbents for a variety of gases, but their limited capacity towards benzene at low concentration remains unresolved. Here we report the adsorption of trace benzene by decorating a structural defect in MIL-125-defect with single-atom metal centres to afford MIL-125-X (X = Mn, Fe, Co, Ni, Cu, Zn; MIL-125, Ti 8 O 8 (OH) 4 (BDC) 6 where H 2 BDC is 1,4-benzenedicarboxylic acid). At 298 K, MIL-125-Zn exhibits a benzene uptake of 7.63 mmol g −1 at 1.2 mbar and 5.33 mmol g −1 at 0.12 mbar, and breakthrough experiments confirm the removal of trace benzene (from 5 to <0.5 ppm) from air (up to 111,000 min g −1 of metal–organic framework), even after exposure to moisture. The binding of benzene to the defect and open Zn(II) sites at low pressure has been visualized by diffraction, scattering and spectroscopy. This work highlights the importance of fine-tuning pore chemistry for designing adsorbents for the removal of air pollutants.

Chemistry↗

DOE OSTI2026

Ultrafast One-Dimensional Peierls-Distortion Dynamics in 1⁢T′−Re⁢S 2 Revealed by 4D Electron Microscopy

Rhenium disulfide (ReS 2 ), a prototypical 2D semiconductor with an anisotropic 1⁢T′ structure due to the pronounced Peierls distortion, has demonstrated great potential for polarization-dependent optoelectronics and photonics. Here, we report an ultrafast phase transition occurring within 1 ps, accompanied by a one-dimensional Peierls-distortion relaxation in 1⁢T′−ReS 2 revealed with combined 4D electron microscopy and time-dependent density functional theory calculations. Upon femtosecond laser pulse excitation, the Re-Re dimerization morphology rapidly transforms from a diamond cluster to zigzag chains and persists for several nanoseconds. Here, this ultrafast Peierls-distortion relaxation is further verified by transient changes in optical anisotropy via polarization-dependent transient absorption spectroscopy. Time-dependent density functional theory calculations attribute this novel phase transition to strong correlation between Peierls distortion and impulsive photoexcited carrier doping, predicting a transient band-gap collapse. This Letter opens up an exciting avenue for ultrafast control of Peierls-distortion-induced anisotropy and the metal-insulator transition in 1⁢T′−ReS 2 .

1T’-ReS2↗

NASA NTRSDate not supplied

A GPU-based Approach for Turbomachinery Application

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Gas Turbine↗

DOE OSTI2024

The effect of gamma ray irradiation on few layered MoSe2: A material for nuclear and space applications

In recent years, emerging two-dimensional (2D) materials, such as molybdenum diselenide (MoSe2), have been at the center of attention for many researchers. This is due to their unique and fascinating physicochemical properties that make them attractive in space and defense applications that include shielding harsh irradiation environments. In this study, we examined the effects of gamma (γ) rays at various doses on the structural, chemical, and optical properties of MoSe2 layers. After the samples were exposed to intense gamma radiation (from a 60Co source) with various exposure times to vary the total accumulated dosage (up to 100 kGy), Raman and photoluminescence spectroscopies were used to study and probe radiation-induced changes to the samples. When compared to pristine materials, very few changes in optical properties were typically observed, indicating good robustness with little sensitivity, even at relatively high doses of gamma radiation. The imaging using scanning electron microscopy revealed a number of nano-hillocks that were connected to substrate alterations. X-ray photoelectron spectroscopies revealed that Mo’s binding energies remained the same, but Se’s binding energies blueshifted. We associated this shift with the decrease in Se vacancies that occurred after irradiation as a result of Mo atoms creating adatoms next to Se atoms. When compared to pristine materials, very few changes in optical, chemical, and structural properties were typically observed. These findings highlight the inherent resilience of MoSe2 in hostile radioactive conditions, which spurs additional research into their optical, electrical, and structural characteristics as well as exploration for potential space, energy, and defense applications.

Materials Science↗

DOE OSTI2024

Salt‐Assisted Vapor–Liquid–Solid Growth of 1D van der Waals Materials

The method of salt-assisted vapor–liquid–solid (VLS) growth is introduced to synthesize 1D nanostructures of trichalcogenide van der Waals (vdW) materials, exemplified by niobium trisulfide (NbS 3 ). The method uses a unique catalyst consisting of an alloy of Au and an alkali metal halide (NaCl) to enable rapid and directional growth. High yields of two types of NbS 3 1D nanostructures, nanowires and nanoribbons, each with sub-ten nanometer diameter, tens of micrometers length, and distinct 1D morphology and growth orientation are demonstrated. Strategies to control the location, size, and morphology of growth, and extend the growth method to synthesize other transition metal trichalcogenides, NbSe 3 and TiS 3 , as nanowires are demonstrated. Finally, the role of the Au–NaCl alloy catalyst in guiding VLS synthesis is described and the growth mechanism based on the relationships measured between structure (growth orientation, morphology, and dimensions) and growth conditions (catalyst volume and growth time) is discussed. These results introduce opportunities to expand the library of emerging 1D vdW materials to make use of their unique properties through controlled growth at nanoscale dimensions.

1D van der Waals material↗

DOE OSTI2023

Effect of the Nature of Both Cation and Anion Substitution on the Structural Symmetry of Li‐Rich 3 d ‐Metal Chalcogenide Electrodes

Abstract Li‐rich layered chalcogenides have recently led to better understanding of the anionic redox process and its associated high capacity while providing ways to overcome its practical limitations of voltage fade and irreversibility. This study reports on the feasibility of triggering anionic activity in Li 2 TiS 3 , through anionic substitution (Se for S) or cationic substitution (Fe for Ti). Herein, the chalcogenide chemical space is further explored to prepare mono‐substituted Li 1.7 Ti 0.85 Mn 0.45 Ch 3 (Ch = S/Se) and doubly substituted cationic and anionic phases (Li 1.7 Ti 0.85 Fe 0.45 S 3‐z Se z ) which crystallize either in the O3‐ or O1‐type structures depending upon substituents. All series show a bell‐shape capacity variation as function of the transition metal (TM) substitution degree with values up to 240 mAh g −1 . For specific compositions, a structural O3 to O1 phase transition is observed upon Li removal, which is not reversible upon Li re‐insertion due to kinetic limitations and negatively affects long‐term cycling performance. Density functional theory (DFT) calculations confirm the O3/O1 relative stability along the different series and point subtle electronic differences in the TM‐doping, rationalizing the structural and electrochemical behaviors of these phases upon cycling. These findings provide further insights into the link between structural and electronic stability, which is of key importance for designing chalcogenide‐based anionic redox compounds.

Chemistry↗