Search NASASearch

SEARCH · Search NASA

Results for “large eddy simulation (LES)”

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.

915 records · Page 12

Performance of the Satellite Test Assistant Robot in JPL's Space Simulation Facility

An innovative new telerobotic inspection system called STAR (the Satellite Test Assistant Robot) has been developed to assist engineers as they test new spacecraft designs in simulated space environments. STAR operates inside the ultra-cold, high-vacuum, test chambers and provides engineers seated at a remote Operator Control Station (OCS) with high resolution video and infrared (IR) images of the flight articles under test. STAR was successfully proof tested in JPL's 25-ft (7.6-m) Space Simulation Chamber where temperatures ranged from +85 C to -190 C and vacuum levels reached 5.1 x 10 -6 torr. STAR's IR Camera was used to thermally map the entire interior of the chamber for the first time. STAR also made several unexpected and important discoveries about the thermal processes occurring within the chamber. Using a calibrated test fixture arrayed with ten sample spacecraft materials, the IR camera was shown to produce highly accurate surface temperature data. This paper outlines STAR's design and reports on significant results from the thermal vacuum chamber test.

Douglas McAffee

Integrated Simulation of Weld Residual Stress Evolution and Crack Propagation Using XFEM

Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.

42 - ENGINEERING

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Wettability of Two-Dimensional Carbon Allotropes from Molecular Simulations

Force-field Monte Carlo and Molecular Dynamics simulations are used to compare wetting behaviors of model carbon sheets mimicking neat graphene, its saturated derivative, graphane, and related planar allotropes penta-graphene, γ-graphyne, and ψ-graphene in contact with aqueous droplets or an aqueous film confined between parallel carbon sheets. Atomistic and area-integrated surface/water potentials are found to be essentially equivalent in capturing moderate differences between the wetting free energies of tested substrates. Despite notable differences in mechanical and electric properties of distinct allotropes, the predicted allotrope/water contact angles span a narrow window of weakly hydrophilic values. Contact angles in the range of 80 ± 10° indicate modest hydration repulsion incapable of competing with van der Waals attraction between carbon particles. Poor dispersibility in neat water is hence a common feature of studied materials.

Biochemistry & Molecular Biology

In-Space Laser Forming Simulated Under Thermal Vacuum (TVAC) Conditions

NASA Marshall Space Flight Center is investigating laser forming under thermal vacuum conditions representative of the space environment through a partnership with the DARPA NOM4D program. The thermal gradient mechanism (TGM) is used to bend a metal sheet towards the laser optics. Laser equipment and process parameters are being optimized for thermal vacuum conditions. Such conditions are representative of the extreme cold and hot temperatures and high vacuum that exist on orbit or on the lunar surface. in situ photogrammetry and thermography assess the bend/counterbend and elucidate the TGM in thermal vacuum. Structured light scanning and metallography after bending reveal the macroscopic and microscopic structure after TGM.

thermal vacuum

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

Accurate prediction of aerothermodynamic loads in thermal non-equilibrium flows requires precise modeling of internal energy exchange. While previous direct simulation Monte Carlo frameworks have successfully implemented discrete rotational energy models for diatomic species, the treatment of polyatomic molecules has traditionally relied on continuous energy assumptions that break down at low temperatures and neglect critical high-temperature corrections. This study extends the discrete rotational energy models of Boyd and Gimelshein to fully encompass polyatomic molecules. The proposed framework implements quantized rotational energy level sampling for linear, spherical, and symmetric/asymmetric top rotors. Crucially, the model incorporates centrifugal distortion to address the limitations of the rigid-rotor assumption at hypersonic temperatures, and accounts for nuclear spin parity, which dictates the permissible rotational states and macroscopic specific heats at low temperatures. The model is verified through equilibrium sampling procedures, demonstrating agreement with theoretical quantum Boltzmann distributions and accurately reproducing thermophysical properties across a wide range of temperatures.

DSMC

Thermomechanical Modeling of Woven Materials With Particle-Based, Explicit-Fiber Simulations

Fiber-based materials are extensively used to protect spacecraft during entry. Insulative fibers, often in a fiber network or woven, provide rigidity, strength, and control of material anisotropy and density. Woven thermal protection materials, such as ADEPT (Adaptable, Deployable Entry and Placement Technology), 3D-MAT (3-Dimensional Multifunctional Ablative Thermal Protection), and 3MDCP (3D Woven Mid-Density Carbon Phenolic), enable missions with stronger and denser materials for entry profiles with high shear and heat flux. Vulnerabilities to woven thermal protection materials include manufacturing-induced material property variation, and impact from micrometeoroids. Simulating woven materials under these conditions require models that can resolve hierarchal structures, thermomechanical behavior, and failure. To address this, we simulate weave thermal conduction and mechanical deformation. We simulate the full weave with a coarse-grained yarn model is presented. The model combines a validated, high-resolution single 3MDCP yarn model and phenolic resin model. Instead of modeling every fiber, each yarn ply with order 10, instead of order 1000, fibers. The discrete element bonded particle model (DEM-BPM) of fibers captures the thermal and mechanical behavior within and between fibers. We study the proportion of heat transfer and stress via the contact network, fiber bonds, and overall weave geometry.

bonded particle

Tensorized Interior Radiative Heat Transfer for a Scalable and Calibrated Building Energy Simulator

Building energy simulation is a critical tool for developing and testing advanced control strategies, such as Reinforcement Learning (RL), to provide demand flexibility and affordable energy costs. The recently introduced Smart Buildings Control Suite (sbsim) provides a lightweight, scalable, and data-calibrated simulation environment based on a 2D finite-difference model. However, the initial model primarily focused on conductive and convective heat transfer, neglecting the significant impact of long-wave radiative heat exchange between interior surfaces. This paper presents a significant extension to the sbsim framework by incorporating a physically-grounded model for interior radiative heat transfer. Our primary contribution is the development and integration of a fully tensorized radiative heat transfer module, which preserves the computational efficiency and scalability of the original simulator. This was achieved by developing a pipeline for view factor calculation, including an algorithm to identify directly seeing surfaces within complex floor plans, and formulating the net radiation equations for efficient execution on modern hardware accelerators. We validate the numerical accuracy of our tensorized implementation by comparing its results against a traditional iterative approach, demonstrating identical outcomes. This enhancement increases the physical fidelity of sbsim, enabling more accurate training of RL agents for building energy optimization.

Ham, Sang woo

Relativistic Magnetohydrodynamic Simulations of Giant Magnetar Bursts

Gradual crustal deformation can generate strongly twisted magnetic fields around magnetars, potentially triggering giant flares with total energies exceeding 10 44 erg. In this letter, we present the first relativistic magnetohydrodynamic simulation of a surface shear-driven magnetar eruption, capturing reconnection-driven plasma heating, the ejection of relativistically hot plasma, and the formation of a hot fireball confined within the inner magnetosphere. We find that magnetic reconnection in the equatorial current sheet launches a hot trailing outflow capable of powering the initial spike observed in giant flares, while simultaneously leaving behind a thermally stratified fireball with sufficient thermal energy to produce the pulsating, decaying tail. Together, these features provide a self-consistent physical framework for understanding the observed energetics of magnetar giant flares. The eruption also expels a magnetically dominated giant plasmoid carrying up to ∼9% of the magnetosphere’s total magnetic energy. Furthermore, our simulation demonstrates how the plasmoid drives the formation of a blast wave—an important ingredient in models linking magnetar eruptions to fast radio bursts.

79 ASTRONOMY AND ASTROPHYSICS

Simulating the Lunar Thermal Environment for the Surface-Deployed LEMS Artemis III TVAC Test

The Lunar Environment Monitoring Station (LEMS) is an autonomous, survive-the lunar-night seismic suite to be deployed on the Lunar surface by the Artemis III crew and designed to operate continuously for two years. It will see the extremes of the Lunar south pole thermal environment where surface temperatures range from -200°C to +20°C and where nighttime duration is at least 354 hours. Given power and mass constraints, the thermal system is limited to 2.5 Watts of heat during the lunar night. Presented is the system thermal vacuum test plan and configuration to validate the thermal control system performance in a flight-like environment. The test will measure the nighttime heat leaks for the bus and 18 seismometers to ensure lunar night survival and operability.

Thermal

Simulating the Lunar Thermal Environment for the Surface-Deployed LEMS Artemis III TVAC Test

The Lunar Environment Monitoring Station (LEMS) is an autonomous, survive-the lunar-night seismic suite to be deployed on the Lunar surface by the Artemis III crew and designed to operate continuously for two years. It will see the extremes of the Lunar south pole thermal environment where surface temperatures range from -200°C to +20°C and where nighttime duration is at least 354 hours. Given power and mass constraints, the thermal system is limited to 2.5 Watts of heat during the lunar night. Presented is the system thermal vacuum test plan and configuration to validate the thermal control system performance in a flight-like environment. The test will measure the nighttime heat leaks for the bus and 18 seismometers to ensure lunar night survival and operability.

Thermal

Simulating the Lunar Thermal Environment for the Surface-Deployed LEMS Artemis III TVAC Test

The Lunar Environment Monitoring Station (LEMS) is an autonomous, survive-the lunar-night seismic suite to be deployed on the Lunar surface by the Artemis III crew and designed to operate continuously for two years. It will see the extremes of the Lunar south pole thermal environment where surface temperatures range from -200°C to +20°C and where nighttime duration is at least 354 hours. Given power and mass constraints, the thermal system is limited to 2.5 Watts of heat during the lunar night. Presented is the system thermal vacuum test plan and configuration to validate the thermal control system performance in a flight-like environment. The test will measure the nighttime heat leaks for the bus and 18 seismometers to ensure lunar night survival and operability.

TVAC

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Discrete Element Method

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Woven

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science

The Use of Immersion Rigs for High Temperature Hydrogen Exposure Testing Within the Nuclear Thermal Rocket Element Environmental Simulator (NTREES): NTREES-5 Hydrogen Heat Exchanger (N5HHX-01)

The Nuclear Thermal Rocket Element Environmental Simulator (NTREES) facility at NASA’s Marshall Space Flight Center (MSFC) subjects materials and components typically found on nuclear thermal propulsion (NTP) systems to non-nuclear simulations of prototypic NTP reactor thermohydraulic environments. With 1.2 MW of inductive power available, the currently employed methods to simulate nuclear heating in NTREES consist of 1) inductively heating a test specimen, which typically cannot be performed on ceramic materials, or 2) heating gas that envelops the test specimen. The latter of these has provided a basis for developing so-called “immersion rigs” to be used in NTREES. Furthermore, immersion rigs are especially useful in the case of smaller, technically challenging test specimens. One such device, the “NTREES-5 Hydrogen Heat Exchanger” (N5HHX-01), designed and deployed to test small material coupons to ~0.5 g/s of 2800 K, 500 psi Hydrogen, is discussed herein. To achieve this desired delivery temperature level, the most suitable material for the N5HHX’s construction was pure tungsten, which imposed budgetary and manufacturing challenges. Also discussed in this paper is the design of the N5HHX rig’s sample container, which became the first iteration of improving NTREES small material coupon test operations by permitting rapid changeout of different material samples.

Reactor Testing

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Ka-Band Monopulse Antenna-Pointing Systems Analysis and Simulation

NASA's Deep Space Network (DSN) has been using both 70-m and 34-m reflector antennas to communicate with spacecraft at S-band (2.3 GHz) and X-band (8.45 GHz). To improve the quality of telecommunication and to meet future mission requirements, JPL has been developing 34-m Ka-band (32-GHz) beamwave guide antennas. Presently, antenna pointing operates in either the open-loop mode with blind pointing using navigation predicts or the closed-loop mode with conical scan (conscan). Pointing accuracy under normal conscan operating conditions is in the neighborhood of 5 mdeg. This is acceptable at S- and X-bands, but not enough at Ka-band. Due to the narrow beamwidth at Ka-band, it is important to improve pointing accuracy significantly (approx. 2 mdeg). Monopulse antenna tracking is one scheme being developed to meet the stringent pointing-accuracy requirement at Ka-band. Other advantages of monopulse tracking include low sensitivity to signal amplitude fluctuations as well as single-pulse processing for acquisition and tracking. This article presents system modeling, signal processing, simulation, and implementation of Ka-band monopulse tracking feed for antennas in NASA/DSN ground stations.

V Y Lo