Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lunar Surface Thermal Modeling”

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.

776 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Lunar accelerometer network gravitational observatory (LANGO)

With ground-based interferometers detecting hundreds of gravitational-wave (GW) events, GW astronomy has continued to blossom. U.S. and European scientists are developing plans to construct third-generation ground-based interferometers. ESA has proceeded through the mission formulation phase of space-based interferometer in a lower-frequency band, 10 –4 –0.1 Hz. Despite all these exciting developments, there is still a missing frequency band, 0.1–10 Hz. This mid-frequency band is rich with interesting astrophysical events. Coalescence and merger of intermediate-mass black holes (IMBHs) will occur in this frequency band. Coalescing stellar-mass BHs will pass through this frequency band days before they reach the frequency band of LIGO and Virgo. Detection of such signals would enable a mid-frequency detector to issue an advance notice to the high-frequency GW detectors, as well as to optical, x-ray and γ-ray telescopes. We propose Lunar Accelerometer Network Gravitational Observatory (LANGO) to detect GWs in this frequency band. In the first phase (LANGO 1), we propose to deploy four ambient-temperature (250 K) accelerometers in the tetrahedral or in a square configuration on the hemisphere facing the Earth. After successful operation at 250 K, LANGO would be upgraded to a cryogenic (4 K) version with over two orders of magnitude increased sensitivity (LANGO 2), with coherent rejection of seismic noise implemented. LANGO is a full-tensor detector, capable of determining the source direction and wave polarization. Each test mass (TM) is suspended as a pendulum with resonance frequency ∼ 0.01 Hz, thus is only weakly coupled to the lunar surface horizontally. LANGO is designed to detect the relative motion of globally separated, nearly free, TMs by using the Moon as a large quiet platform. LANGO 1 and 2 accelerometers aim at sensitivities ⩽ 10 –11 m s –2 Hz –1/2 and ⩽ 10 –13 m s –2 Hz –1/2 in the horizontal axes over the frequency band of 1 mHz–10 Hz, which yield GW sensitivities 1.1 x 10 -21 Hz -1/2 and 3.8 x 10 -24 Hz -1/2 at 1 Hz, respectively. The LANGO accelerometers will be 10 3 –10 5 times more sensitive than Apollo seismometers. With such sensitivity, LANGO will also make great contribution to the advancement of lunar geophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Modeling Radiolysis and Chemical Reactions during Dry Storage of Aluminum-clad Spent Nuclear Fuel

After aluminum-clad spent nuclear fuel (ASNF) is removed from the reactor, it is initially stored in spent fuel pools, which are specially designed water-filled basins that provide temporary cooling to reduce the temperature of the fuel assemblies and provide radiation shielding. ASNF continues to generate heat due to the radioactive decay of elements within the fuel, which persists for many years post-shutdown as the residual radioactive products decay into more stable elements. During the wet storage period, an oxyhydroxide layer composed of boehmite/bayerite forms on the surfaces of the aluminum cladding from exposure to water in the pools. Road-ready packaging for long-term disposition of the ASNF involves dry storage in helium backfilled DOE standard canisters (DSCs). When the ASNF is removed from water storage and dried, most of the water is removed, but some physisorbed and chemisorbed water remains in the oxyhydroxide layers. This residual water can produce hydrogen when exposed to radiation from the ASNF during dry storage. Predicting hydrogen accumulation over time in the DSCs is critical for long-term storage considerations. Previous modeling efforts have developed coupled computational fluid dynamics (CFD)-chemical models to simulate temperature, pressure, and gas phase concentrations within the DSCs. These models use the thermal field predicted by CFD as input to a radiolysis model for the gas phase and the surface oxyhydroxide layer chemistry. Given the long storage period of the DSCs and the impracticality of long-term experiments, a simulation-based approach is necessary to assess chemical evolution within the canisters. This study advances the development of a modeling framework designed to simulate the chemical evolution of spent fuel canisters. Both thermal and radiation-driven reactions are considered, with radiation kinetics quantified using G-values. Sensitivity analysis identifies key parameters influencing species composition. Reaction pathway diagrams offer insight into dominant species formation routes, enabling more effective comparisons between model predictions and experimental observations, particularly regarding the production of hydrogen. Results show that the model predicts significant hydrogen gas production with minimal oxygen generation, primarily due to hydrogen formation via boehmite pathways. These findings underscore the importance of accurately characterizing surface-bound species and radiolysis kinetics. A deeper understanding of these mechanisms is critical for evaluating the long-term safety of nuclear waste storage.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2026

First-principles DFT modeling of nitrobenzene adsorption on the Ag(111) surface at varying monolayer coverages

Here, the adsorption behavior of nitrobenzene on the Ag(111) surface as a function of coverage is investigated using density functional theory. Adsorption energies and optimized geometries are analyzed together with isolated intermolecular interaction calculations and electronic structure analysis, including Bader charge partitioning and projected density of states, to disentangle the roles of adsorbate–surface bonding and through-space adsorbate–adsorbate interactions. At low and intermediate coverages (θ = 1/9 and θ = 2/9), bidentate adsorption at top sites is favored due to strong adsorbate–surface interactions, with additional stabilization at θ = 2/9 arising from favorable intermolecular separations. At higher coverage (θ = 1/3), bidentate adsorption is destabilized by strong intermolecular repulsion, and monodentate adsorption becomes energetically preferred as rotational freedom allows more favorable intermolecular spacing. Charge density differences, Bader charge, and density of states analyses show that charge transfer from the Ag surface is localized primarily on the nitro group and increases with coverage and adsorption denticity, although this increase does not directly correlate with adsorption strength at high coverage due to competing intermolecular interactions. These results demonstrate that surface coverage can induce a transition in preferred adsorption denticity driven by intermolecular interactions, highlighting the importance of adsorbate packing in organic molecule adsorption at metal interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First in operando transient grating spectroscopy measurements on tungsten during high flux plasma operation in PISCES-RF

This work details the construction and first in operando transient grating spectroscopy measurements conducted on PISCES-RF during plasma operation. A preliminary study on a tungsten sample with varied plasma species (Ar and D 2 ), ion flux (4 × 10 21 to 6 × 10 22 ion/m 2 /s), ion energy (5 to 80 eV), and surface temperature (20 to 200 °C) is presented. Prior to plasma exposure, the thermal diffusivity and surface acoustic wave (SAW) speed, (6.9 ± 0.2) × 10 −5 m 2 /s and (2.66 ± 0.02) km/s, were measured in situ with a thermal grating period of 12.5 μm. During Ar plasma exposure, these two quantities vary according to the sample surface temperature. At low flux and varied ion energy, D 2 plasma exposure also results in thermal diffusivity and SAW speed varying with surface temperature alone. At high flux, D 2 plasma exposure results in no convincing change to the thermal diffusivity, but the SAW speed is reduced. After returning to 20 °C, post plasma exposures, the thermal diffusivity and SAW speed almost recover to the initial value prior to exposure.

Surface acoustic wave SAW↗
DOE OSTI2024

WRF-Comfort: simulating microscale variability in outdoor heat stress at the city scale with a mesoscale model

Abstract. Urban overheating and its ongoing exacerbation due to global warming and urban development lead to increased exposure to urban heat and increased thermal discomfort and heat stress. To quantify thermal stress, specific indices have been proposed that depend on air temperature, mean radiant temperature (MRT), wind speed, and relative humidity. While temperature and humidity vary on scales of hundreds of meters, MRT and wind speed are strongly affected by individual buildings and trees and vary on the meter scale. Therefore, most numerical thermal comfort studies apply microscale models to limited spatial domains (commonly representing urban neighborhoods with building blocks) with resolutions on the order of 1 m and a few hours of simulation. This prevents the analysis of the impact of city-scale adaptation and/or mitigation strategies on thermal stress and comfort. To solve this problem, we develop a methodology to estimate thermal stress indicators and their subgrid variability in mesoscale models – here applied to the multilayer urban canopy parameterization BEP-BEM within the Weather Research and Forecasting (WRF) model. The new scheme (consisting of three main steps) can readily assess intra-neighborhood-scale heat stress distributions across whole cities and for timescales of minutes to years. The first key component of the approach is the estimation of MRT in several locations within streets for different street orientations. Second, mean wind speed and its subgrid variability are downscaled as a function of the local urban morphology based on relations derived from a set of microscale LES and RANS simulations across a wide range of realistic and idealized urban morphologies. Lastly, we compute the distributions of two thermal stress indices for each grid square, combining all the subgrid values of MRT, wind speed, air temperature, and absolute humidity. From these distributions, we quantify the high and low tails of the heat stress distribution in each grid square across the city, representing the thermal diversity experienced in street canyons. In this contribution, we present the core methodology as well as simulation results for Madrid (Spain), which illustrate strong differences between heat stress indices and common heat metrics like air or surface temperature both across the city and over the diurnal cycle.

Geology↗
DOE OSTI2025

The glass-ceiling convective regime and the origin and diversity of coronae on Venus

Venus and Earth are rocky planets of roughly the same size and bulk density, yet their surface volcanic and tectonic features appear substantially different. On Venus, the coexistence of large volcanic highlands—interpreted as the surface expression of long-lived mantle plumes—alongside coronae, smaller features thought to be caused by transient thermal diapirs, remains enigmatic. Using two-dimensional numerical models of mantle convection with sharp and broad mineral phase transitions for pyrolite, we show that both scales of upwellings can be generated in a stagnant lid planet with an interior temperature 250 to 400 K warmer than Earth’s. The smaller plumes originate from a ~600 km deep internal layer that exists as a consequence of the different sequence of mineral phase transitions that occur in warmer mantles less processed and differentiated by partial melting and volcanism. Future models that include melting will provide further tests of our hypothesis.

Science & Technology - Other Topics↗
DOE OSTI2026

Demonstration of gold nanorod systems for enhanced total efficiency: Experimental and numerical analysis

This study investigates the photothermal performance of gold nanorods engineered to exhibit longitudinal plasmon resonances at 695 nm, 780 nm, and 970 nm. The work combines synthesis, structural characterization, extinction measurements, numerical modeling, and controlled temperature experiments to quantify how nanorod geometry, resonance tuning, concentration, and chamber shape jointly influence heat generation. Transmission electron microscopy confirms that increasing nanorod aspect ratio systematically shifts the longitudinal plasmon peak toward the near-infrared region. Extinction measurements show strong agreement with theoretical predictions, validating the numerical model across two independent datasets. Three chamber geometries were tested under laser excitation at 640 nm, 808 nm, and 980 nm: an ascending stepped base, a flat base, and a descending stepped base. Without nanorods, the ascending geometry produced the highest efficiency due to enhanced natural convection. After introducing gold nanorods, all geometries exhibited substantial thermal enhancement, with total efficiencies exceeding 20%. The strongest improvement was obtained for nanorods resonant at 780 nm with a mass concentration of 4.6 mg/mL implemented on the descending stepped-base geometry. This performance resulted from the combined effect of spectral overlapping with the 808 nm laser, the highest nanorod concentration, and localized heat accumulation that intensified buoyancy-driven flow. The findings demonstrate that total efficiency is governed by a synergistic interplay between optical resonance, nanoparticle concentration, and macroscopic chamber design, revealing the system-level coupling between nanoscale plasmonic absorption and macroscale heat-transfer phenomena. The results provide a validated framework for tuning nanoscale plasmonic absorbers and optimizing thermal systems for applications requiring efficient light-to-heat conversion.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Hierarchical Multi-agent Large Language Model Reasoning for Autonomous Heterogeneous Catalyst Discovery

Artificial intelligence is reshaping scientific exploration, but most methods automate procedural tasks without engaging in scientific reasoning, limiting autonomy in discovery. We demonstrate that hierarchical agentic large language model reasoning can efficiently drive simulation and scientific exploration. Across two chemical applications, CO adsorption on Cu surface transition metal adatoms and on M–N–C catalysts, reasoning-guided exploration reduces required atomistic simulations by up to 90% relative to heuristic or random selection. Comparisons across single-agent, multi-agent, and stochastic baselines show that hierarchical strategies yield more coherent and information-efficient search trajectories. Reasoning traces reveal chemically grounded decisions that cannot be explained by semantic bias or stochastic sampling. We realize these agentic reasoning strategies in Materials Agents for Simulation and Theory in Electronic-structure Reasoning (MASTER), a multimodal system that translates natural language into density functional theory workflows. Altogether, multi-agent collaboration accelerates heterogeneous catalyst discovery and marks a step toward more autonomous, reasoning-guided scientific exploration.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

Upper ocean mixing, surface heat fluxes, and heat content variability in the upper 150 m during Hurricane Laura (2020)

Observations (from thermistor sensors and a hurricane glider) and model-based temperature were used to answer questions including the following: 1) What physical and environmental conditions explain the rapid intensification of Hurricane Laura? 2) How did the pre-existing warm mixed layer at Stone mooring (StM) influence the degree of cooling in the mixed layer? 3) How did vertical mixing, surface heat fluxes, and advective processes collectively shape the mixed layer heat evolution at StM? Observations measured by thermistor sensors showed that Hurricane Laura induced a cooling of -1.2oC in the mixed layer when Hurricane Laura’s wind speed increased to 47 m/s on 26 August. This contrasts with the -1.04oC mixed layer temperature change estimated from model-based temperature. At StM, the presence of a 31oC warm mixed layer and elevated heat content (60-80kJ/cm2) effectively preconditioned the upper ocean ahead of Hurricane Laura’s passage. After Hurricane Laura passed, turbulent mixing associated with the storm transported the heat anomalies downward into subsurface layers. In addition, this study shows that a surface heat flux of 5.04 kJ/cm2 supported the intensity of Hurricane Laura as it traversed over the StM. Using a mixed layer heat balance model, this study shows that entrainment flux and surface flux mostly contributed to the observed change in the mixed layer temperature.

Tsai, Senam↗
DOE OSTI2027

Low thermal budget dopant activation in shallow junctions of implanted Si via Tunable plasma enhanced annealing

Low thermal budget (LTB) annealing has become increasingly critical for shallow junctions as semiconductor devices are scaled down. Plasma Enhanced Annealing (PEA) has emerged as a promising LTB annealing process, however its mechanisms remain unclear. In this study, arsenic and boron implanted silicon wafers were subjected to helium or argon ion bombardment generated by annealing plasmas under controlled ion energies and doses in a home-built annealing reactor. The structural and electrical changes were characterized by four-point probe measurements, secondary ion mass spectrometry, Raman spectroscopy and spectroscopic ellipsometry. A pronounced reduction in sheet resistance, accompanied by a decrease in the damage layer thickness demonstrates the surface selective annealing capability of PEA process. Comparative studies using He and Ar plasmas and related ion bombardment reveal a dependence of activation efficacy on both ion species and the dopant implantation profile. Furthermore, a multi-step annealing mechanism is proposed, consisting of an initial (low dose) structural recovery stage followed by the generation of “extended” defects and concluded at large doses by an enhanced dopant activation. The proposed annealing kinetics are thought to be controlled by ion flux, dose and energy. Finally, these results highlight PEA as an effective, surface-selective, LTB approach for shallow-dopant activation and near surface annealing, and provide mechanistic insights into PEA processes.

dopant activation↗
DOE OSTI2022

In Situ Tracking of Nonthermal Plasma Etching of ZIF-8 Films

Surface characterization is critical for understanding the processes used for preparing catalysts, sorbents, and membranes. Nonthermal plasma (NTP) is a process that achieves high reactivity at low temperatures and is used to tailor the surface properties of materials. In this work, we combine the capabilities of infrared reflection absorption spectroscopy (IRRAS) with NTP for the in situ interrogation of zeolitic imidazolate framework-8 (ZIF-8) thin films to probe modifications in the material induced by oxygen and nitrogen plasmas. The IRRAS measurements in oxygen plasma reveal etching of organic ligands with sequential removal of the methyl group and imidazole ring and with the formation of carbonyl moieties (C═O). In contrast, nitrogen plasma induces mild etching and grafting of nitrile groups (−C≡N). Scanning electron microscopy imaging shows that oxygen plasma, at prolonged times, significantly degrades the ZIF-8 film at the grain boundaries. Treatment of ZIF-8 membranes using mild plasma conditions yields a fivefold enhancement for H 2 /N 2 and CO 2 /CH 4 ideal selectivities and an eightfold enhancement for CO 2 /N 2 ideal selectivity. Additionally, the new tools described here can be used for spectroscopic in situ tracking of plasma-induced chemistry on thin films in general.

IRRAS↗
DOE OSTI2026

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)↗
DOE OSTI2026

Simple device modification simultaneously enhances indoor passive evaporative cooling power and duration

Passive evaporative cooling from a damp cloth or vessel can be used for localized cooling of perishables such as food and medicine, drinking water, and even people. Such applications have two key objectives, which are normally in tension with each other: minimizing the water consumption rate and maximizing the cooling (steady-state temperature swing between the cool cloth and the warm surroundings). Here, we present a simple device that delivers higher performance in both metrics simultaneously by adding a perforated aluminum foil sheet suspended over a stagnant air layer to thermally shield the evaporating surface while also facilitating mass transfer. Experimental results demonstrate a simultaneous 10%–25% increase in temperature swing and 2–3× reduction in water consumption rate as compared to conventional evaporative cooling. This improvement is achieved through careful modeling to optimize the coupled heat and mass transfer through the airgap (thermally insulating and vapor permeable) and perforated foil (infrared reflective and vapor permeable).

Chen, Sarah Mizuno↗
DOE OSTI2026

Observed Land Surface Influence on Atmospheric Heat and Moisture Profiles During Interstorms

Land-atmospheric (L-A) feedbacks have historically been studied using models whose structure and parameterizations influence outcomes and insights. The representation of L-A feedbacks based on observations alone remains an ongoing challenge for understanding boundary layer development and precipitation. To address this gap, we use ground-based passive remote sensing and in-situ observations to present an analysis of the atmosphere during 103 interstorm soil moisture drydown events spanning nine warm seasons (2016–2024) in the U.S. Southern Great Plains region. By separating events based on local L-A coupling signals and characterizing the profiles of atmospheric heat and moisture to surface energy flux behavior, we investigate the physical mechanisms linking land surface processes to boundary layer development. We find that during interstorm drydowns, the atmospheric column follows a consistent pattern: moisture increases within the boundary layer, peaks near its top, and declines rapidly above, while warming occurs through the depth. Drydowns that shift toward evaporation produce stronger and deeper thermodynamic responses than cases dominated by sensible heating, which are weaker and shallower. Additionally, moisture is accumulated faster within the boundary layer during shorter drydowns, with longer drydowns representing slower, moisture-limited growth. Drydowns with wetter initial soil moisture will sustain stronger moistening within and above the boundary layer, accelerating buoyancy growth and convective potential toward the next storm. These results provide observational evidence linking surface flux evolution to boundary layer thermodynamics and offer a process-level benchmark for evaluating coupled L-A representations in models and demonstrating the influence of soil moisture on short-term weather forecasting skill.

Zhang, M. S. [Massachusetts Inst. of Technology (M↗
DOE OSTI2026

Mechanical, Electrochemical & Thermal Modeling of Electric Vehicle Batteries for Crash Simulation (CRADA CRD-19-00811 Final Report)

Under the proposed effort in partnership with Hyundai Motor Company (HMC), the National Laboratory of the Rockies (NLR) will cooperate with HMC to develop mathematical models for battery cells and modules for simulating abuse response in batteries subject to the type of mechanical crushing that can occur in a full motor vehicle crash.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Mechanical, Electrochemical & Thermal Models - Training (CRADA CRD-19-00798 Final Report)

Under the proposed effort in partnership with Hyundai Motor Company (HMC), the National Renewable Energy Laboratory (NLR) will host Dr. Jaeyoung Lim from HMC for a period of one year to jointly develop mathematical models for battery cells and modules subject to mechanical crush. NLR will assist with the development of mathematical models that Dr. Lim will incorporate into his research effort on new concepts of mobility with electric vehicles.

33 ADVANCED PROPULSION SYSTEMS↗