Search NASASearch

SEARCH · Search NASA

Results for “astrophysical 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.

638 records · Page 16

Swamp Works Regolith Compaction Technologies

While the level of compaction below the lunar surface increases quickly after only a few cm of depth, in many cases during a construction mission there will be a need to excavate and transport regolith to a new location for cut-and-fill or horizontal construction of structures such as berms. In these cases to achieve high levels of bulk density, compaction must be per-formed. Additionally, in some cases surface technolo-gies such as systems that sinter/melt the surface may desire the maximum possible compaction at the sur-face to improve melting/heating performance and the final material strength properties. Kennedy Space Center’s (KSC) Swamp Works has developed two means of compaction, lunar and mar-tian compaction. Planetary Autonomous Compaction Technology (PACT) which is part of the Multifunction End Effector for Regolith Compaction Acquisition and Transfer (MEERCAT [1]) robotic arm end effector system’s capabilities and the Site Preparation Tooling for Operations on Mobility Platforms (STOMP [2]) vibratory roller compactor. PACT on MEERCAT has been demonstrated to a TRL 5 and STOMP to a TRL 4 in ambient testing. The results of PACT on MEERCAT and STOMP testing will be shared with results for various simulants including BP-1, ICN-LHT-1G (aka CSM-LHT-1G), RDW-LHT-1GH (a simulant developed for the Mason Tipping Point to match characteristics of ICN-LHT-1G), and Exolith LHS-1E. This will also include discus-sions on methods used to verify relative density before and after compaction and means to verify density effects below depth. To calculate relative density, maximum and minimum densities for simulants were taken from literature [3] [4] [5] and additional lab test-ing (publication in work).

vibration

Modeling of a Four-Stage Linear Ionization Cooling Channel for a Muon Collider in G4Beamline

A previous study of an eight-stage rectilinear ionization cooling channel in the ICOOL software demonstrated a five-order-of-magnitude reduction in a muon beam’s 6D emittance. In this study, we look to compare the ways ICOOL and Muons, Inc.’s g4Beamline software model ionization cooling by comparing their modeling of the first four stages to this optimized cooling channel constructed in ICOOL. We begin by identifying the parameters used to construct the optimized ionization cooling channel in ICOOL. We then reconstruct this beam in g4Beamline with identical parameter specifications and simulate the cooling of an identical input beam. Finally, we compare the two simulations based on their beam transmission, longitudinal emittance, and transverse emittance along the channel length. Through this process, we demonstrate that G4Beamline accurately reproduces transverse cooling results but predicts systematically different longitudinal emittance evolution while maintaining similar overall cooling performance, reproducing a 97.9% reduction in 6D emittance over four stages.

Keeler, Dominic [Purdue U., West Lafayette] (ORCID

Development of a silver-zinc battery system

Summary report is described of historical documentation and detailed design data for development of silver-zinc battery for use on Surveyor spacecraft. Electrical and physical characteristics of battery models are included, along with data on qualification, acceptance, solar-thermal-vacuum, mission simulation testing, and actual flight performance.

Moses, A. J.

CAE for Thermal Management of Aerospace Electronic Boards Using the BETAsoft Program

Aerospace electronic boards require special attention to thermal management due to constraints such as their need to be light, small, and maintain high power densities. Also, cooling is mainly through conductive and radiative modes with minor or negligible convective cooling. Due to these particular requirements, thermal design has become an integrated part of the electronic design process in order to avoid expensive repeat prototyping and to ensure high reliability. To achieve high speed simulations, the BETAsoft code uses semi-empirical formulations and an advanced finite difference scheme that incorporates local adaptive grids. Detailed conduction, convection and radiation heat transfer is considered. Various benchmark verifications of the software simulation compared to infrared images typically prove to be within 10% of each other. The thermal analysis of a sample avionic card in a natural convection environment is shown. Then, the individual effects of attaching metal screws to the casing, increasing radiative emissivities of the casing, increasing the conductance of the wedge lock, adding an aluminum core to the board, adding metal strips in board layers, inserting conduction pads under components, and adding heat sinks to components are demonstrated.

Kimberly Bobish

Preliminary Assessment of the Impact on the V-Band Oxygen Channels From Satellite Communication Uplinks

We calculate the percentage of time that an ATMS-like instrument [1] will be illuminated by the uplink beam of one of the proposed V-band communication system and estimate the damage resulting from such exposure. Using a combination of openly available information and educated guesses about the location and characteristics of the up/down link terminals, we constructed the ground segment of a hypothetical high-speed communication network. The space segment of the network was constructed from the orbital data of the existing Starlink constellation [2] of 6223 communication satellites (comsats) which is used as strawman to represent any other possible constellation of communication satellites. It is shown that without a very delicate balance of frequency allocations (science vs telecommunications), coupled with extremely steep and deep bandpass-defining filters, and strict adherence to the agreed limits (i.e. no out-of-band transmissions) the deployment of the telecommunication network leads to almost-complete loss of some important geophysical data. For the analysis we use the spectral characteristics of the ATMS instrument with the ephemeris for the NOAA-21 satellite [3]. The analysis is conducted for the USA and the simulation covers 8 consecutive days in July 2024. Effective and accurate microwave remote sensing of the atmosphere depends on the availability of interference-free spectrum windows at frequencies which are prescribed by physical processes [e.g. 4]. The family of resonant lines of the oxygen molecule near 60 GHz provides a unique opportunity to sample the vertical distribution of temperature and density from space, and it has been exploited for weather and climate studies from polar-orbiting satellites since 1978 (MSU on TIROS-N [5]). It remains a staple in the payloads operated by Russia, China, USA, Japan, France, India, UK, Ukraine [6] which are built around a common blueprint: a few wide-band (hundreds of MHz) channels around 50 GHz to sample the atmosphere and the surface while several more channels with high spectral resolution (few MHz) sample the individual resonant lines. Accurate retrieval of the environmental parameters depends upon the data provided by both sets of channels, and the their location in frequency space is not arbitrary and cannot be altered at will [7, 8]. The introduction of 5G technology in 2019 has driven telecommunication companies to request more bandwidth to be dedicated to their devices. This additional bandwidth is only available in spectral regions traditionally reserved for environmental and astrophysical research, such as the V-band between 50 and 60 GHz for up/downlink between satellites in low-earth orbits and terminals connected to fiberoptics network for distribution to high-speed local internet services. The power broadcast by the uplink communication leg is many orders of magnitude greater than the natural thermal signal emitted from the Earth scene. If the ground antenna were to perfectly align with the passive instrument’s antenna, the spaceborne receiver would suffer permanent, irreparable damage. While a direct boresight-to-boresight conjunction is extremely unlikely (even with a large constellation of satellites the fraction of the celestial sphere occupied by the satellites remains minuscule) the finite size of the ground station’s antenna beam in the sky suggests that the ATMS will be in the near background (as seen from the ground station) of one of the communication satellites and will be illuminated by either the main lobe or the near sidelobes of the uplink antenna more often than it is desirable. For our analysis we first calculate the position of the ATMS with respect to each of the ground stations at a resolution of 0.2 sec, then calculate the position of each of the comsats which are at least 25 deg above the station’s local horizon; finally we calculate the angle between the line-of-sight of the ATMS and the line-of-sight of the comsat. We assume that the gain pattern of the ground station is circularly symmetric; the angle-off-station-boresight then provides an attenuation of the uplink power which we use to assess the likely effect upon the passive instrument’s operations. We assume that each ground station can communicate with all the comsats in its field of view; this implies that, on average, a ground station can engage with 46 comsats simultaneously. The analysis is repeated for the case when the uplink broadcast within the ATMS passive channels (in-band scenario) and for the case when the uplink is limited to frequencies adjacent to the ATMS channels (out-of-band scenario). The antenna of the ground station is modelled as having a HPBW (Half-Power Beam Width) of 0.16 deg and EIRP (Equivalent Isotropic Radiated Power) of 70 dBW. We account for the geometric dissipation of the signal caused by the satellite orbital altitude, the attenuation induced by atmospheric gasses at 51 GHz and the mismatch between the circular polarization of the ground-based transmitting antenna and the linear polarization of the satellite-borne receiving antenna. The damages on ATMS are estimated from bench-level measurement conducted at the ATMS’ manufacturer facilities [unpublished].

passive microwave

Controlled Parametric Forcing During Directional Solidification of a Bulk Organic Alloy Under Microgravity

The response of dendritic microstructures to step-like pulling velocity conditions is investigated using microgravity directional solidification experiments conducted on DECLIC-DSI combined with phase-field simulations. Under a constant pulling velocity of 1.5 µm/s, the evolution toward steady-state growth is characterized in terms of primary spacing, dendrite drift, and tip dynamics. For the first time, side-view observations enabled direct measurement of tip radius and sidebranching frequency. When step-like oscillations of the pulling velocity are imposed, the dendritic array exhibits a strongly period-dependent response: short periods lead to rapid tip adaptation, whereas longer periods induce a phase lag between tip position and morphology, resulting in progressive tip flattening and, above a critical period, interface destabilization and dendrite splitting. Quantitative phase-field simulations, including a realistic thermal field and stochastic noise, reproduce the experimental observations and provide insight into the governing mechanisms, highlighting the role of characteristic relaxation times, sequence-dependent effects, and the irreversible reorganization of the microstructure following splitting.

Microgravity

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)

Optimization of a 10–Step Reduced CH₄/O₂ Combustion Mechanism for RDRE Applications

A 10-step, 9-species reduced mechanism for methane/oxygen combustion is evaluated and optimized for use in rotating detonation rocket engine (RDRE) simulations at an operating pressure of 10 bar. The mechanism is assessed against GRI-Mech 3.0 for ignition delay times, adiabatic equilibrium temperature and species, Chapman-Jouguet (CJ) detonation parameters, and laminar flame speed. A sequential Nelder–Mead optimization procedure is applied to the rate constants of reactions R1, R4, and R5, targeting ignition delay agreement over the equivalence ratio range φ = 1.0–1.4 and temperature range T = 1400–2000 K. The optimized mechanism achieves ignition delay agreement within 28% of GRI-Mech 3.0 for the temperature range relevant to RDRE operation (τ ≥ 5 μs). Adiabatic equilibrium temperatures and CJ detonation parameters are reproduced within 0.15% of GRI-Mech 3.0 values, confirming thermodynamic consistency of the mechanism. Laminar flame speed predictions deviate significantly from GRI-Mech 3.0, which is attributed to the absence of HO₂ chemistry in the reduced mechanism; this limitation does not affect RDRE simulation fidelity since laminar flame propagation is of secondary importance in RDREs.

Computational Fluid Dynamics

Solar-Wind Bombardment of a Surface in Space

The solar wind is described and its effects on surfaces are reviewed. Experimental results are given for materials subjected to simulated solar-wind bombardment. Among the effects discussed are sputtering-type erosion, chemical reaction with paint vehicles and reduction of oxide pigments, and production of lunar-surface optical characteristics in a basalt powder.

Solar Wind

A Radiometric Consistent Spectral Fingerprinting Algorithm for Continuity Products of Hyperspectral Sounders

A radiometric consistent climate fingerprinting methodology has been developed to derive long-term temperature, water vapor, cloud, trace gases, and surface skin temperature anomaly time series from the hyper-spectral sounder measurements of multiple platforms. The spectral fingerprinting methodology requires the use of radiative kernels that are radiometrically consistent with observations. Radiative kernels are built using space-time averaged Jacobians that are physically retrieved from observations under all sky conditions. The physical retrieval algorithm uses the Principal Component based Radiative Transfer Model (PCRTM) for the forward simulation. The incorporation of multiple scattering simulation in PCRTM allows the direct radiative relationship between single field-of-view (FOV) radiance observations and corresponding thermal dynamic variables including cloud properties to be established. Therefore, radiance ?closure? can be achieved under all-sky conditions by the fingerprinting scheme. This methodology has been used to derive climate anomalies from the space-time averaged spectra of AIRS/AMSU and CrIS/ATMS. The use of a consistent fingerprinting scheme provides an effective mean of generating continuity product by merging observations from different platforms and therefore facilitating the long-term climate trend study.

Wan Wu

General Electric 20-ampere hour nickel-cadmium battery

The interaction, effect, and controllability of the performance parameters of the General Electric 20-ampere-hour, 24-cell nickel cadmium battery are investigated. The battery was cycled under simulated orbit conditions. The acquired data was analyzed and evaluated in terms of battery parameters and performance characteristics. Conclusions and tests results are presented along with recommendations for further study.

Kirsch, W. W.

Valorization of consolidated bioprocessing residues for bioplastics

This study demonstrates an organic solvent-free processing strategy to valorize consolidated bioprocessing (CBP) residues, from switchgrass and poplar biomass, into functional poly(butylene succinate) (PBS)-based biocomposites using high-shear homogenization (HSH). HSH transformed the switchgrass and poplar CBP residues (CBP-R) into fine, uniformly distributed particles and microfibers. The composites of PBS with homogenized switchgrass residues (H-CBP-R-SG) or homogenized poplar residues (H-CBP-R-P) at a 70/30 weight ratio exhibited improved processability and mechanical integrity, with the Young's modulus for the PBS/H-CBP-R-SG and PBS/H-CBP-R-P nearly doubling to 0.66 ± 0.07 GPa and 0.65 ± 0.04 GPa, respectively, compared to neat PBS (0.36 ± 0.02 GPa). Dynamic Mechanical Analysis (DMA) reveals a significant suppression of the tan δ peak magnitude, indicating that HSH-mediated physical activation facilitates stress transfer in composites typical of covalent chemical grafting systems. While the transition to a stiffness-dominated profile reduces ductility, the resulting composites exhibit the dimensional stability and resistance to thermal warping required for high-fidelity FDM 3D printing and injection molding. Beyond material performance, comprehensive techno-economic analysis (TEA) and life cycle assessment (LCA) confirmed that diverting CBP residues into composite products can improve the economic viability of the biorefinery without substantially increasing biorefinery global warming potential (GWP). At a 30 wt% blend ratio, incorporating residuals into PBS yielded a minimum selling price for the composite of $\$4.07$ per kg compared to the conventional bioplastic price of $\$5.00$ per kg. This approach aligns with circular bioeconomy principles by converting waste streams into value-added products. Furthermore, this innovative strategy addresses key challenges in bioplastic development, including cost, compatibility, and performance, while simultaneously advancing waste minimization strategies for sustainable manufacturing systems.

09 BIOMASS FUELS

Socio-Techno-Economic Feasibility of Deep Geothermal with EGS for Residential District Heating in the Northeastern United States

In this paper, we present socio-techno-economic simulation results for enhanced geothermal system (EGS) district heating systems in the northeastern United States. We applied the geospatial Distributed Geothermal Market Demand (dGeo) tool coupled with GEOPHIRES to evaluate the feasibility of geothermal deep direct-use, utilizing EGS reservoirs for providing residential heating with district heating systems. The results were combined with JEDI results to assess social impacts. Simulation results assuming doublet EGS reservoirs with target production temperature of 80 degrees C and newly constructed district heating systems indicate a wide range in system levelized cost of heat (LCOH) values ranging from ~$15/MMBtu to over $1,000/MMBtu, with most the attractive regions urban areas with high total thermal demand and high thermal demand density. We estimate about 60 GWth in potential installed capacity with LCOH values under $50/MMBtu and 100 GWth in potential capacity with LCOH values under $100/MMBtu, respectively.

15 GEOTHERMAL ENERGY

A Probabilistic Approach to Load Modeling for Central HVAC Systems in Large Commercial Buildings for Retrofit Decisions Under Uncertainty

Retrofitting central HVAC systems in large commercial buildings with advanced technologies like heat recovery chillers (HRCs) offers a significant opportunity to enhance energy efficiency. However, analyzing these retrofits is challenging with traditional whole-building simulation tools, which require intensive calibration and struggle to model innovative system configurations and controls. To overcome these limitations, this study proposes a load profilebased retrofit analysis framework that provides better decisions under uncertainty. The main focus of this paper is the development of a probabilistic load profile model that can be used in the framework by using exploratory data analysis (EDA) of measured building data to properly quantify its inherent variability. A non-parametric Gaussian Process (GP) model was employed to capture the time- and weather-dependent characteristics of the heating load while explicitly modeling its uncertainty. The model's effectiveness is demonstrated through strong predictive performance on unseen data and physically interpretable insights into load behavior. This data-driven, probabilistic load profile serves as a robust and flexible input for subsequent system simulations, enabling a more confident and statistically sound analysis of retrofit potential.

Ham, S W

AT2025ulz and S250818k: Leveraging DESI Spectroscopy in the Hunt for a Kilonova Associated with a Subsolar-mass Gravitational-wave Candidate

On 2025 August 18, the LIGO–Virgo–KAGRA collaboration reported a subthreshold gravitational-wave candidate detection consistent with a subsolar-mass neutron star merger, denoted S250818k. An optical transient, AT2025ulz, was discovered within the localization region. AT2025ulz initially appeared to meet the expected behavior of kilonova emission, the telltale signature of a binary neutron star merger. The transient subsequently rebrightened after ∼5 days and was classified as a Type IIb supernova. In this work, we analyze the observations of its host galaxy obtained by the Dark Energy Spectroscopic Instrument (DESI). From the DESI spectrum, we obtain a secure redshift of z = 0.084840 ± 0.000006. If S250818k has an astrophysical origin, this places the transient within 2σ of the gravitational-wave distance and results in an integral overlap between the gravitational-wave alert and the transient location of log 10 $\mathcal{I} ≈ 3.9 - 4.2$. Our analysis of the host galaxy’s spectral energy distribution reveals a star-forming, dusty galaxy with stellar mass ∼10 10 M ⊙ , broadly consistent with the population of both short gamma-ray bursts and core-collapse supernova host galaxies. We also present our follow-up of DESI-selected candidate host galaxies using the Fraunhofer Telescope at the Wendelstein Observatory, and show the promise of DESI for associating or rejecting candidate electromagnetic counterparts to gravitational-wave alerts. These results emphasize the value of DESI’s extensive spectroscopic dataset in rapidly characterizing host galaxies, enabling spectroscopic host subtraction, and guiding targeted follow-up.

Hall, Xander J. [Carnegie Mellon University, Pitts

Vidyut3d: A Gpu Accelerated Fluid Solver for Non-Equilibrium Plasmas on Adaptive Grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

Sitaraman, Hariswaran

Progress on Shape Memory Alloy Actuator Development for Active Clearance Control

Results of a numerical analysis evaluating the feasibility of high-temperature shape memory alloys (HTSMA) for active clearance control actuation in the high-pressure turbine section of a modern turbofan engine has been conducted. The prototype actuator concept considered here consists of parallel HTSMA wires attached to the shroud that is located on the exterior of the turbine case. A transient model of an HTSMA actuator was used to evaluate active clearance control at various operating points in a test bed aircraft engine simulation. For the engine under consideration, each actuator must be designed to counteract loads from 380 to 2000 lbf and displace at least 0.033 in. Design results show that an actuator comprised of 10 wires 2 in. in length is adequate for control at critical engine operating points and still exhibit acceptable failsafe operability and cycle life. A proportional-integral- derivative (PID) controller with integrator windup protection was implemented to control clearance amidst engine transients during a normal mission. Simulation results show that the control system exhibits minimal variability in clearance control performance across the operating envelope. The final actuator design is sufficiently small to fit within the limited space outside the high-pressure turbine case and is shown to consume only small amounts of bleed air to adequately regulate temperature.

Jonathan DeCastro

Superdiffusion resilience in Heisenberg chains with two-dimensional interactions on a quantum processor

Superdiffusive spin transport in the one-dimensional (1D) Heisenberg model is a key theoretical discovery in nonequilibrium quantum many-body physics. Although extensively studied in 1D systems, the breakdown and sustenance of superdiffusion in two-dimensional (2D) lattices with integrability-breaking terms, as found in real materials, remains an open question. To address this, we develop a toy model that extends the 1D Heisenberg model with a representative set of 2D interaction types and tunable strengths. Our model exhibits varying degrees of superdiffusion breakdown depending on the interaction type, spanning ballistic to diffusive regimes. We establish and justify a hierarchy of 2D interactions based on their resilience against superdiffusion breakdown: Heisenberg >𝑋⁢𝑋 > Ising. This precise control over the superdiffusive behavior also enables rigorous benchmarking of quantum hardware, and our simulations on IBM's Heron devices confirm the hardware's ability to accurately capture these many-body nonequilibrium phenomena. Overall, our results are relevant not only to simulating superdiffusion in real materials, such as the 1D Heisenberg compound KCuF3, which contains modest nonintegrable 2D terms, but also to extending superdiffusive behavior to larger 2D qubit lattices and other 2D materials.

Alagarsamy Manikandan, Keerthi Kumaran [ORNL]