Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Non-linear Simulation”

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.

842 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Spike-Nosed Bodies and Forward Injected Jets in Supersonic Flow

The paper contains new numerical simulation and experimental test results of blunt body drag reduction using thin spikes mounted in front of a body and one- or two-phase jets injected against a supersonic flow. Numerical simulations utilizing the NASA CFL3D code were conducted at the Hampton University Fluid Mechanics and Acoustics Laboratory (FM&AL) and experimental tests were conducted using the facilities of the IM/MSU Aeromechanics and Gas Dynamics Laboratory. Previous results were presented at the 37 t h AIAA/ASME/SAE/ASEE Joint Propulsion Conference. Those results were based on some experimental and numerical simulation tests for supersonic flow around spike-nosed or shell-nosed bodies, and numerical simulations were conducted only for a single spike-nosed or shell-nosed body at zero attack angle, alpha=0 o . In this paper, experimental test results of gas, liquid and solid particle jet injection against a supersonic flow are presented. In addition, numerical simulation results for supersonic flow around a multiple spike-nosed body with non-zero attack angles and with a gas and solid particle forward jet injection are included. Aerodynamic coefficients: drag, C D , lift, C L , and longitudinal momentum, M z , obtained by numerical simulation and experimental tests are compared and show good agreement.

M Gilinsky↗
NASA NTRSDate not supplied

Spike-Nosed Bodies and Forward Injected Jets in Supersonic Flow

The paper contains new numerical simulation and experimental test results of blunt body drag reduction using thin spikes mounted in front of a body and one- or two-phase jets injected against a supersonic flow. Numerical simulations utilizing the NASA CFL3D code were conducted at the Hampton University Fluid Mechanics and Acoustics Laboratory (FM&AL) and experimental tests were conducted using the facilities of the IM/MSU Aeromechanics and Gas Dynamics Laboratory. Previous results were presented at the 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conference. Those results were based on some experimental and numerical simulation tests for supersonic flow around spike-nosed or shell-nosed bodies, and numerical simulations were conducted only for a single spike-nosed or shell-nosed body at zero attack angle, alpha = 0 degrees. In this paper, experimental test results of gas, liquid and solid particle jet injection against a supersonic flow are presented. In addition, numerical simulation results for supersonic flow around a multiple spike-nosed body with non-zero attack angles and with a gas and solid particle forward jet injection are included. Aerodynamic coefficients: drag, C (sub D), lift, C(sub L), and longitudinal momentum, M(sub z), obtained by numerical simulation and experimental tests are compared and show good agreement.

M Gilinsky↗
DOE OSTI2026

Experimental investigation of a closed vapour box module for a divertor-like configuration in Magnum-PSI

Efficient management of extreme heat fluxes in the divertor region to extend the lifetime of the components remains a critical challenge for the realization of nuclear fusion-based power plants. Among the alternative concepts explored for the divertor region, the use of liquid metals, particularly lithium, is of interest due its ability to dissipate the incoming plasma heat flux through the vapour shielding effect (VS). In this work, we experimentally investigated a ‘closed’ configuration of a dedicated Vapour Box Module (VBM) in the linear plasma device Magnum-PSI. The goal of the experiments is to simulate the vapour box divertor environment conditions and assess its performance in terms of power mitigation and redistribution and lithium confinement. Initial testing without Li demonstrated the efficacy of a closed VBM structure in inducing detachment via neutral gas accumulation. Apertures which enabled non-condensing gas to be effectively pumped while ensuring lithium condensed on the inner surfaces were therefore added. With a lithium capillary porous structure target used, lithium is directly vaporized by the plasma, forming a dense lithium vapour cloud that interacts with the incoming plasma. This resulted in a significant reduction of the target temperature of at least 48%, together with a temperature locking effect, a phenomenon typically observed in the VS regime. Lithium vapour confinement within the VBM was strongly correlated with the wall temperature. Relatively cold walls promoted Li re-condensation and therefore improved Li confinement, although with the expected trade-off of increased hydrogenic retention on lithium-wetted surfaces. As the wall temperature increased, the confinement efficiency decreased, consistent with reduced Li re-condensation and thermally activated Li–H chemistry and remobilization at the walls. Diagnostic measurements through embedded thermocouples and calorimetry revealed that lithium vaporization and re-condensation processes also playedsignificant roles in plasma power dissipation. The results advance the case for a closed divertor chamber with direct lithium evaporation from the strike-points as a viable method to manage divertor heat fluxes in future fusion reactors.

Romano, Fabio [Dutch Institute for Fundamental Ene↗
DOE OSTI2025

Geospatial modeling of near subsurface temperatures of the contiguous United States for assessment of materials degradation

Abstract Understanding subsurface temperature variations is crucial for assessing material degradation in underground structures. This study maps subsurface temperatures across the contiguous United States for depths from 50 to 3500 m, comparing linear interpolation, gradient boosting (LightGBM), neural networks, and a novel hybrid approach combining linear interpolation with LightGBM. Results reveal heterogeneous temperature patterns both horizontally and vertically. The hybrid model performed best achieving a root mean square error of 2.61 °C at shallow depths (50–350 m). Model performance generally decreased with depth, highlighting challenges in deep temperature prediction. State-level analyses emphasized the importance of considering local geological factors. This study provides valuable insights for designing efficient underground facilities and infrastructure, underscoring the need for depth-specific and region-specific modeling approaches in subsurface temperature assessment.

Science & Technology - Other Topics↗
NASA NTRSDate not supplied

The Use of Immersion Rigs for High Temperature Hydrogen Exposure Testing within the Nuclear Thermal Rocket Element Environmental Simulator (NTREES): Thermal Soak Rig (TSR)

The Nuclear Thermal Rocket Element Environmental Simulator (NTREES) facility was purpose constructed to perform non-nuclear evaluations of nuclear thermal propulsion (NTP) system fuel materials and structures within prototypic thermochemical environments. This system has been utilized steadily in its ability to subject test specimens to thermochemical and thermohydraulic environments simulating that of an operating nuclear rocket engine. Fission heat is simulated by induction power and experiments are conducted within a ~1000 psi pressure vessel. Hydrogen is conventionally passed through the heated fuel surrogate test specimen while pressure, temperature, and gas species data are collected at various points along the experiment. In order to test fuel and material coupon samples, a class of test apparatus named “immersion rigs” are being developed and employed to more rapidly test these smaller and more technically challenging test specimen. One example of a promising potential fuel structure, Tristructural-isotropic (TRISO) particles, presents unique challenges for testing of this type. TRISO fuel micro-particles are spheroids typically on the order of 500 – 1000 μm in diameter, and exposing a batch sample to hot hydrogen requires purpose-built special test equipment. Thusly, an immersion rig was developed and successfully demonstrated to expose ~1 g of micro-particles to hydrogen gas at temperatures and pressures relevant to NTP systems for the purpose of fuel evaluation. The rig, comprised primarily of graphite and pure tungsten, houses in its core a batch of micro-particles between pucks of porous silicon carbide (SiC). This approach permits gas flow while simultaneously retaining the particles in place. Herein is a discussion of the design, analysis, fabrication, and testing of the NTREES Thermal Soak Rig (TSR).

Space Nuclear Propulsion↗
NASA NTRSDate not supplied

Spacecraft Fire Safety Predictions using Verified Saffire Model

A model developed using Fire Dynamics Simulator (FDS) that aimed to determine the effect of a fire in a spacecraft was validated by data collected during the Saffire campaign. The model used inlet and outlet temperatures and CO 2 concentrations of the Saffire payload where fire spread was taking place to determine the amount of heat and combustion products that made it into Northrop Grumman’s Cygnus vehicle. The model was then validated using six remote sensors in various places in the vehicle, as well as a far field device (FFD) in the open zenith section that was representative of average vehicle values. The current work focuses on using the model to predict fire safety scenarios. One simulation aimed to determine the fate of HCl, which sticks to surfaces. The model prediction showed that the HCl was removed from the atmosphere rapidly. This compared well against the FFD data in the Saffire VI campaign event where a bottle of 5% HCl was released into the vehicle. An additional simulation where the Environmental Control and Life Support System (ECLSS) was shut off once the FFD reached 5 ppm of HCl showed that HCl stayed in the atmosphere considerably longer. Continuing the simulation with the ECLSS activated and after temperatures returned to their initial conditions, resulted in a rapid removal of HCL similar to that observed in the original HCl release scenario model. Finally, a simulation that used the heat release rate from a lithium-ion battery test to determine the effect it would have on a spacecraft was performed. This simulation used a heat addition rate that is considerably higher than what was determined from the burning of solid fuels in the Saffire campaign and hence produced a non-trivial temperature increase in more locations within the vehicle.

Fire Safety↗
NASA NTRS2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λc, that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λc-based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

Amanda Stricklan↗
NASA NTRS1966

Non-Metallic Materials

Nonmetallic materials development - cryogenic insulation, adhesives research, and membrane diffusion theory

Materials Science↗
NASA NTRS1966

Non-magnetic battery case Patent

Nonmagnetic hermetically sealed battery case made of epoxy resin and woven glass tape for use with electrochemical cells in spacecraft

Hennigan, T. J.↗
DOE OSTI2026

Constraining primordial non-Gaussianity from DESI DR1 quasars and Planck PR4 CMB lensing

We present the first measurement of local-type primordial non-Gaussianity from the cross-correlation between 1.2 million spectroscopically confirmed quasars from the first data release (DR1) of the Dark Energy Spectroscopic Instrument (DESI) and the Planck PR4 CMB lensing reconstructions. The analysis is performed in three tomographic redshift bins covering 0.8 < z < 3.5, covering a sky fraction of ∼20%. We adopt a catalog-based pseudo-C ℓ estimator and apply linear imaging weights validated on noiseless mocks. Compared to previous analyses using photometric quasar samples, our results benefit from the high purity of the DESI spectroscopic sample, the reduced noise of PR4 lensing, and the absence of excess large-scale power in the spectroscopic quasar auto-correlation. Fitting simultaneously for the non-Gaussianity parameter f NL and the linear bias amplitude in each redshift bin, we obtain f NL = 2 +28 -34 for a response parameter p = 1.6, and f NL = 6 +20 -24 for p = 1.0. These results improve the constraints on f NL by ∼35% compared to the previous analysis based on the Legacy Imaging Survey DR9. Additionally, we derive an optimal weighting scheme to maximize the constraining power. In this case, and assuming p = 1.6, we obtain f NL = 19 +25 -31 . Our results demonstrate the statistical power of DESI quasars for probing inflationary physics, and highlight the promise of future DESI data releases.

cosmological parameters from CMBR↗
NASA NTRSDate not supplied

Material Properties and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

PATO↗
NASA NTRSDate not supplied

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG↗
NASA NTRS1998

High Temperature - Thin Film Strain Gages Based on Alloys of Indium Tin Oxide

A stable, high temperature strain gage based on reactively sputtered indium tin oxide (ITO) was demonstrated at temperatures up to 1050 C. These strain sensors exhibited relatively large, negative gage factors at room temperature and their piezoresistive response was both linear and reproducible when strained up to 700 micro-in/in. When cycled between compression and tension, these sensors also showed very little hysteresis, indicating excellent mechanical stability. Thin film strain gages based on selected ITO alloys withstood more than 50,000 strain cycles of +/- 500 micro-in/in during 180 hours of testing in air at 1000 C, with minimal drift at temperature. Drift rates as low as 0.0009%/hr at 1000 C were observed for ITO films that were annealed in nitrogen at 700 C prior to strain testing. These results compare favorably with state of the art 10 micro-m thick PdCr films deposited by NASA, where drift rates of 0.047%/hr at 1050 C were observed. Nitrogen annealing not only produced the lowest drift rates to date, but also produce the largest dynamic gage factors (G = 23.5). These wide bandgap, semiconductor strain sensors also exhibited moderately low temperature coefficients of resistance (TCR) at temperatures up to 1100 C, when tested in a nitrogen ambient. A TCR of +230 ppm/C over the temperature range 200 C < T < 500 C and a TCR of -469 ppm/C over the temperature range 600 C < T < 1100 C was observed for the films tested in nitrogen. However, the resistivity behavior changed considerably when the same films were tested in oxygen ambients. A TCR of -1560 ppm/C was obtained over the temperature range of 200 C < T < 1100 C. When similar films were protected with an overcoat or when ITO films were prepared with higher oxygen contents in the plasma, two distinct TCR's were observed. At T < 800 C, a linear TCR of -210 ppm/C was observed and at T > 800 C, a linear TCR of -2170 DDm/C was observed. The combination of a moderately low TCR and a relatively large gage factor make these semiconducting oxide films promising candidates for the active strain elements in high temperature thin film strain gages, particularly in applications where static strain measurement is desired.

Otto J Gregory↗
DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗
NASA NTRS2000

Chalcogenide-Based Non-Volatile Memory Technology

Chalcogenide is a proven phase change material used in re-writeable CDs and DVDs. This material changes phases, reversibly and quickly, between an amorphous state that is dull in appearance and electrically high in resistance, and a polycrystalline state that is highly reflective and low in resistance. The application of this commercially proven technology to create dense, high-speed, non-volatile semiconductor memories is discussed.

J Maimon↗
NASA NTRSDate not supplied

Atmospheric, Non-Combusting Droplet Sizing and Spray Imaging Results From a Non-Proprietary Aero Gas Turbine Engine Pressure Swirl Atomizer using Shadowgraphy and Planar Laser Scatter

Members from industry and federal agencies of the United States have developed a collaborating group focused on improving the understanding of the atomization processes common to gas turbine engines. The group is aimed on providing high quality data using x-ray diagnostics in the internal geometry and the near-field dense spray as well as optical diagnostics in the visible spectrum, downstream of the nozzle. This paper reports results from one of the two non-proprietary injectors designs from the group: the primary-secondary injector. In particular, we characterize the spray from only the primary, pressure-swirl atomizing circuit using water through the nozzle. We observed Sauter mean diameter decreased with increased water pressure drop across the nozzle. SMD was also collected at three axial distances and showed decreased size with distance from the nozzle. Planar laser sheet droplet scattering provided qualitative assessment of the spray. Higher water pressure drops provided better atomization than lower pressure drops. The addition of swirling air flow through the injector alleviates and improves the atomization quality, particularly at the lower water pressure drops.

optical diagnostics↗