Search NASASearch

SEARCH · Search NASA

Results for “Neural ordinary differential equations”

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.

74 records · Page 3

Effect of Turbulence Models on Criticality Conditions in Swirling Flows

The critical state of vortex cores downstream of vortex breakdown has been studied. Base vortical flows were computed using the Reynolds-averaged, axisymmetric Navier-Stokes equations. Standard K - epsilon , RNG and second-order Reynolds stress models were employed. Results indicate that the return to supercriticality is highly dependent on the turbulence model. The K - epsilon model predicted a rapid return of the vortex to supercritical conditions, the location of which showed little sensitivity to changes in the swirl ratio. The Reynolds stress model predicted that the vortex remains subcritical to the end of the domain for each of the swirl ratios employed, and provided results in qualitative agreement with experimental work. The RNG model produced intermediate results, with a downstream movement in the critical location with increasing swirl. Calculations for which area reductions were introduced at the exit in a subcritical flow were also performed using the Reynolds stress model. The structure of the resulting recirculation zone was altered significantly. However, when area reductions were employed within supercritical flows as predicted using the two-equation models, no significant influence on the recirculation zone was noted.

Robert E Spall

PERSIANN-Unet: A Global Deep Learning Framework for Near-Real-Time Precipitation Estimation Using Infrared Data

Access to high-quality, high-resolution, near-real-time precipitation data is essential for hydrological and meteorological research and disaster mitigation. Traditional tools such as rain gauges and radar networks, though effective, have limitations, including sparse coverage in remote areas and high operational costs. Satellite data, with its global coverage and high spatial and temporal resolutions, mitigates limitations in coverage. Satellite precipitation products like Hydro Estimator (HE), Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (IMERG), and Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks (PERSIANN) utilize both geosynchronous thermal infrared (IR) and passive microwave (PMW) data in their operation. PMW sensors offer detailed atmospheric profiles but suffer from higher latency, whereas IR sensors provide lower latency but only capture cloud-top information. Despite this constraint, IR data remains attractive for low-latency precipitation estimation. Recent advances in deep learning, particularly convolutional neural networks (CNNs), have further improved satellite precipitation retrievals. This study introduces PERSIANN-Unet (PUnet or PERSIANN V3), a quasi-global algorithm covering 60°N–60°S that combines IR data, monthly climatology, and the UNet architecture to produce half-hourly precipitation estimates at 0.04° resolution. The product is evaluated against HE, IMERG, and PDIR-Now for 2022–2023. Results show that PUnet closely matches its training target, IMERG V07 Final, at the global scale, and performance is further evaluated against Stage IV as a reference over CONUS. Training PUnet on IMERG (2016–2021) leverages a high-quality, integrated PMW IR-gauge precipitation product while developing an IR-based framework not reliant on PMW availability. By operating on a single global image, PUnet avoids tile partitioning and blending steps, reducing edge discontinuities, and produces more spatially consistent precipitation fields across hemispheres.

Phu Nguyen

Turbo-Design: Open-Source Radial Equilibrium Turbomachinery Solver: Part I - Turbines

Advances in 3D Geometrical Designs and Cooling have played a significant role in improving the efficiency of turbomachinery. However, these advancements must be effectively translated back to the modeler. Machine learning can facilitate this transition. Specifically, machine learning–based loss models can bridge the gap between 3D and 1D designs, enabling modelers not only to predict velocity triangles but also to extract additional geometric features. Currently, the design tools used at NASA have not been updated to support such integration—until now. TurboDesign is an open-source, Python-based framework that replaces TD2 (LEW-11029-1) and AXOD2 (LEW-16323-1), both of which are radial equilibrium solvers for axial turbines. The goal of this update is to enable the integration of machine learning loss models into radial equilibrium equations. Additionally, TurboDesign is designed to support radial machines. This paper presents the governing equations, the assumptions underlying the code, the integration of legacy loss models, an example of machine learning model integration, and a validation comparison with CFD. All code, tutorials, and documentation are available at: https://www.github.com/nasa/turbo-design

Radial Equilibrium

Analysis of Impact Induced Damage and its Effect on Structural Integrity of Space Flight Composite Overwrapped Pressure Vessels

The objective of this research work has been to provide analytical background and support to the ongoing experimental program at NASA, White Sands Test Facility, involving testing composite overwrapped pressure vessels (COPV) for impact damage and cyclic pressurization. Preliminary theoretical basis, including the governing equations for a shallow shell subjected to internal pressure, has been established. Effects of the Griffith type cracks on the structural integrity of the cylindrical vessel were evaluated by methods of Fracture Mechanics. The results indicate that the effective mass of the pressure vessel is an important factor influencing the response to impact events. We also have found that the material properties of the target, contained in the constitutive equations of the composite attached to the Aluminum liner, dominate the impact event in the low velocity range, the material properties become less important, while the target mass distribution and the impactor mass become more significant as the velocity of the impactor increases. Therefore, at high-velocity impact it is not only the kinetic energy of the impactor but also its mass which has a significant effect on the dynamics of the event, and consequently on the induced damage. This work also suggests a methodology for an assessment of the rate of loading effects on the degradation of the material toughness associated with a high-velocity impact where the rate effects become significant. To model the rate dependence of the material response a viscoelastic-plastic constitutive equations were assumed, and on this basis predictions are made regarding the rate dependent material resistance curve. Other dynamic phenomena associated with the impact event have been treated in the framework of the Computational Mechanics using the courtesy of Prof. P. Guebelle and his graduate student at University of Illinois at Urbana-Champaign who have an access to a super-fast computer located on their campus. Finally, the guidelines for a follow-up research program are provided in the body of this report. They address three major areas: theoretical research, numerical studies, and further experimental work.

Michael P Wnuk

Form Factors, Grey Bodies and Radiation Conductances (Radks)

With today's analysis tools, large, complex thermal radiation problems are easily solved; But, as with any analytical tool, lack of an understanding of the fundamental equations and technique limitations may leave you with the wrong answer; Whether you are a new engineer or a seasoned veteran, an understanding of the techniques employed by these powerful analysis tools is crucial.

Fluids

Improvements to RANS Modeling for Aeroheating Predictions on Blunt Bodies

Accurate predictions of aeroheating are critical for designing thermal protection systems for planetary entry vehicles. For larger vehicles, turbulence in the boundary layer can substantially increase convective heating. This turbulence must be accurately modeled to ensure the thermal protection system is sufficient. The majority of hypersonic turbulence model development and validation focuses on boundary layers developing over flat-plates or sharp cones; these cases are substantially different than the boundary layer that develops over the heatshield of a blunt body traveling at hypersonic speeds. Planetary missions often use blunt body geometries, such as the 70-degree sphere-cone favored by Mars missions or the 45-degree sphere-cone planned for the upcoming DAVINCI mission. Due to smaller vehicle size and the lower velocities in the stagnation region, planetary entry vehicles have relatively low Reynolds numbers. Surface curvature and high enthalpy gradients create additional challenges. These difficulties must be addressed to obtain high accuracy needed for the ambitious planetary missions in the upcoming decade. This work focuses on both assessing and improving Reynolds-averaged Navier-Stokes (RANS) turbulence models for blunt-body geometries typical of planetary entry vehicles, with a focus on one-equation and two-equation formulations.

Mars2020

Third Annual Workshop on Space Operations Automation and Robotics (SOAR 1989)

Papers presented at the Third Annual Workshop on Space Operations Automation and Robotics (SOAR '89), hosted by the NASA Lyndon 8. Johnson Space Center at Houston, Texas, on July 25-27, 1989, are documented herein. During the three days, approximately 100 technical papers were presented by experts from NASA, the USAF, universities, and technical companies. Also held were panel discussions on Air Force/NASA AI Overview and Expert System Verification and Validation. Tutorial sessions included Neural Networks; Theory and Application of Back Propagation; Verification and Validation of Expert Systems/ Evaluation of Expert System Tools; and Technical Environment for Modular Architectures for Robotics in Space; and are not documented herein. Technical topics addressed included intelligent systems, robotics, human factors, and environment.

Knowledge representation

Measure sea level air pressure from space to improve knowledge and forecasting of the atmospheric state

Modern numerical weather prediction (NWP) and analysis models require globally-observed meteorological data including sea-level pressure (SLP) for accurate operations. Up until now, SLP has only been measured by in situ instruments from ships, buoys, and ocean platforms. These measurements are sparse with large gaps, leaving models starved of this critical information to constrain the atmospheric state. Recent advancements in differential absorption radar (DAR) provide a path to close this critical observation gap through spaceborne observations in the coming decade, improving the analysis models relied upon for atmospheric research and the weather forecasts depended upon daily for public safety and commerce.

Matthew L Walker McLinden

Emergence of Complex Modes in Lightly Damped Structural Systems

In this short presentation format, we will provide a set of useful equations, relationships, and physical interpretations direct from the damped two-mode interaction problem. We will demonstrate that the results from the two-mode interaction problem are in direct agreement with observed test data. Practical methods for estimating the parameters driving the relations direct from modal test data are provided. We finalize with two examples demonstrating weak and strong interactions.

spaceflight hardware

Time-History Statistics of Soot Formation in A Model Gas Turbine Combustor

Soot formation is a complex dynamic and intermittent process determined by properties of the fuel, combustor design, and combustor operation. Although the major steps in soot formation (i.e., formation of precursors, inception, growth and evolution) are similar for a variety of carbonaceous fuels, applications, and operating conditions, it remains unclear when the temporal transition between these steps occurs. An engineering prediction tool coupled with computational fluid physics (CFD), therefore needs to accurately model all these complex steps. To develop such a model, we propose the time-history concept for understanding the time dependency of soot formation as a function of local properties (i.e., temperature, velocity, local fuel air ratio, etc.). We continue our previous work with modeling the DLR aero-combustor [1] with our updated in-house CFD code, Open National Combustion Code (OpenNCC), that now includes a Multiple Time-Scale Flamelet Progress Variable approach and a the semi-empirical two-equation soot model. We injected massless tracer particles upstream of the injector region of the combustor to collect time-history statistics of the solution variables. The correlations between the collected statistics with respect to the experimental soot volume fraction data showed that time-history effect of certain flow variables, including turbulent kinetic energy (TKE), and multiple species is indeed important for soot formation. We then conducted a time-history based correlation analysis to determine the key species and the concentration ranges critical for soot formation (C6H5-based nucleation, acetylene-based surface growth, and oxidation with OH and O2). Based on the time-history correlation coefficient (THCC) analysis, we propose possible modifications to improve the current two-equation model.

LES

Computation of Optimal Interplanetary Low-Thrust Trajectories With Bounded Thrust Magnitude By Means of the Generalized Newton-Raphson Method

The generalized Newton-Raphson method, an iterative procedure for solving nonlinear operator equations, has been extended in application to variational problems with bounded control variables. A minimum fuel interplanetary low thrust orbital transfer problem is worked out in detail to demonstrate the practical aspects of the algorithm as well as its computational effectiveness. The control variables are the thrust magnitude, limited from zero to some prescribed maximum value, and the thrust steering angle.

Computation

Report of the Microbial Development Working Group

In formulating ideas on the relationship of gravity to the development, growth, and reproduction of microorganisms, a rather liberal definition of microorganisms is used which includes bacteria, yeasts, protists, filamentous fungi, and single cells in culture. A principal advantage of microorganisms as experimental subjects is the rigor with which they can be defined and controlled. As single cells, each cell may be regarded as identical to the others in the population. This property applies to the morphology, physiology, and genetic parameters of the cells. The growth and development of the population is subject to precise manipulation as the nutritional requirements are known and minimal media formulations have been developed. Growth and differentiation can be manipulated in a variety of ways, such as alteration of the culture temperature and food supply, or by use of mutants. Finally, the short generation times of microorganisms provide the opportunity to conduct multigenerational studies within practical time limits and, in a similar vein, cellular responses to various stimuli or stresses are conveniently monitored because of the rapid response times of single cells.

G Nelson

Evaluation of Cycle Life and Characterization of YTP 45 Ah Li-Ion Battery for EMU

Li-ion batteries, with longer cycle life and higher energy density features, are now more and more attractive and applied in multiple fields. The YTP 45 Ah Li-ion battery has been evaluated here and may be employed in EMU in the future. Evaluations were on: (1) Cycle life tests - 500 cycles total (completed 40 cycles in simulated shuttle use mode and 460 cycles in an accelerated use mode, and recorded differential voltage of individual cell in battery); (2) Characterization test - discharge capacity measurement in environment temperature of -10, 25, 50 C before and after 500 cycles; and (3) Thermal testing - charge and discharge at 50 C and -10 C before and after 500 cycles. The battery showed less than a 9% drop of initial discharge capacity and energy within 500 cycles with 475 cycles 59% DOD plus 25 cycles 100% DOD. The EOD voltage ranged from 16.0 to 18.0 V, which fits the requirement for operating the EMU.

Deng, Yi

Determination of Local Experimental Heat-Transfer Coefficients on Combustion Side of an Ammonia-Oxygen Rocket

Local experimental heat-transfer coefficients were measured in the chamber and throat of a 2400-pound-thrust ammonia-oxygen rocket engine with a nominal chamber pressure of 600 pounds per square inch absolute. Three injector configurations were used. The rocket engine was run over a range of oxidant-fuel ratio and chamber pressure. The injector that achieved the best performance also produced the highest rates of heat flux at design conditions. The heat-transfer data from the best-performing injector agreed well with the simplified equation developed by Bartz at the throat region. A large spread of data was observed for the chamber. This spread was attributed generally to the variations of combustion processes. The spread was least evident, however, with the best-performing injector.

Curt H Liebert

Detailed examination of LDEF's frame and the A0178 thermal blankets by the meteoroid and debris special investigations group

A responsibility of the group is to define the hypervelocity particle environment encountered by LDEF during its stay in low Earth orbit. LDEF's 6061-T6 aluminum frame and the 'Teflon silver-inconel paint' thermal blankets represent large surface areas that were widely distributed around the spacecraft. The results are reported of high resolution scans of approx. 0.36 and approx. 0.31 sq m for the intercostals and thermal blankets, respectively. The findings are in qualitative agreement with existing model predictions that suggest high differential bombardment histories for surfaces pointing into specific directions relative to the velocity vector of a non-spinning platform in LEO. The production rate for craters greater than or = 50 microns in diameter in aluminum and penetration holes greater than or = 100 microns in diameter in thin foil materials differ by more than a factor of 10 between forward and rearward facing surfaces. These are substantial differences that must be considered during the design of future long duration space platforms in LEO.

Thomas H. See

Thermal Cycling and Isothermal Deformation Response of Polycrystalline NiTi: Simulations vs. Experiment

A recent microstructure-based FEM model that couples crystal-based plasticity, the B2<-> MB190 phase transformation and anisotropic elasticity at the grain scale is calibrated to recent data for polycrystalline NiTi (49.9 at.% Ni). Inputs include anisotropic elastic properties, texture and differential scanning calorimetry data, as well as a subset of recent isothermal deformation and load-biased thermal cycling data. The model is assessed against additional experimental data. Several experimental trends are captured - in particular, the transformation strain during thermal cycling monotonically increases and reaches a peak with increasing bias stress. This is achieved, in part, by modifying the martensite hardening matrix proposed by Patoor et al. [Patoor E, Eberhardt A, Berveiller M. J Phys IV 1996;6:277]. Some experimental trends are underestimated - in particular, the ratcheting of macrostrain during thermal cycling. This may reflect a model limitation that transformation-plasticity coupling is captured on a coarse (grain) scale but not on a fine (martensitic plate) scale.

Phase Transformations

Boundary Layer Analysis of Shock Tube Flows

Shock tubes offer a controlled environment to reproduce kinetic and radiative phenomena characteristic of atmospheric entry flows under ground-test conditions. The boundary layer developing behind the incident shock wave determines the available test time, influences particle residence times important for similarity scaling, and can affect radiative energy transport. In this study, we couple a quasi-1D space marcher with the compressible boundary-layer equations to numerically compute the post-shock flow in a shock frame of reference for various test gas mixtures representative of different planetary atmospheres. The solvers are individually verified against CFD simulations and analytical correlations available in the literature. Coupled solutions are computed for a finite-rate chemistry in the boundary layer and a non-catalytic isothermal wall. Results yield refined estimates of the maximum separation distance, as well as insights into concentration profiles of relevant species within the boundary layer.

Andrea Fagnani

Mechanics of Preloaded Bolt Tensile Loading With Focus on Load Introduction Factor

The bolt tensile and joint separation loads are directly influenced by the locations at which the external loads enter the clamped members of a preloaded bolted joint (PBJ) and the associated load-paths through the joint. This physical load introduction mechanism affecting the bolt tensile loading is typically represented in the bolt tensile load equation, in part, by a load introduction factor (LIF), which was shown by H.M. Lee of Marshall Spaceflight Center to be a natural product of the bolt tensile load equation using a linear spring stiffness model. This LIF, being a function of load-path stiffness, has subsequently been denoted as the stiffness-based LIF (SBLIF), providing a framework to calculate the LIF using whatever load-path stiffness approximations are appropriate. Expanding upon the work of Lee, it is shown that the SBLIF and the joint stiffness factor are functions of the stiffnesses of the same load-paths and regions within a PBJ, and thus they should not be treated as independent variables. Mathematical expressions for the SBLIF are presented. Comparisons are shown between the analytically calculated SBLIF, the analytically calculated geometric LIF (GLIF), which is a simple clamped-member thickness ratio, the experimentally derived LIF, and the LIF determined by finite element analysis (FEA). Using experiment and FEA as a benchmark, the SBLIF, using traditional load-path stiffness approximations, enables a more accurate prediction of bolt tensile loading than the GLIF, although it can be unconservative near joint separation. The GLIF generally attributes more of the externally applied tensile load to the bolt than does the SBLIF, potentially resulting in heavier and/or more costly bolted joints. Mathematical relationships between the SBLIF and the GLIF are developed. Supplemental material is provided in the appendixes where the historical practice of using the joint compressive stiffness in place of the joint tensile stiffness is evaluated. The appendixes include step-by-step examples demonstrating the calculation of the SBLIF using traditional stiffness approximations and conclude with the development of the joint diagram in terms of the SBLIF, culminating into formulas for the key features of a joint diagram, which is useful for programming.

Load Path