Results for “Neural ordinary differential equations” — Page 2

NASA NTRS1992

Comparison between the PISO algorithm and preconditioning methods for compressible flow

Two widely used family of algorithms, pressure-based and density-based methods, have been developed for computational fluid dynamics (CFD) problems over the years. Pressure-based methods (such as SIMPLE and PISO) use a Poisson-like equation for updating pressure instead of the continuity equation, while density-based methods use the continuity equation to update density (an equation of state is used to provide density in pressure based schemes and pressure in density based schemes). Pressure-based methods were developed originally for incompressible flows at low Reynolds numbers and were then extended to high Reynolds numbers and compressible applications. On the other hand, density based methods were originally developed for transonic flows and have been extended down to low Mach numbers through the use of preconditioning techniques. We compare these two very different approaches to solving the Navier-Stokes equations in order to gain an understanding of their similarities and differences. Specifically, we consider the PISO scheme as a representative pressure-based method and contrast it with a recently developed preconditioning scheme. We also compare the relative performance of the PISO algorithm with a Euler implicit algorithm that is employed to solve the preconditioned equations by means of a vector stability analysis.

Charles L Merkle↗

NASA NTRS1977

The Evolution of the Moon and the Terrestrial Planets

The thermal evolutions of the Moon, Mars, Venus, and Mercury are calculated theoretically starting from cosmochemical condensation models. An assortment of geological, geochemical, and geophysical data are used to constrain both the present-day temperatures and the thermal histories of the planets' interiors. Such data imply that the planets were heated during or shortly after formation and that all the terrestrial planets started their differentiations early in their history. Initial temperatures and core formation play the most important roles in the early differentiation. The size of the planet is the primary factor in determining its present-day thermal state. The Moon, smallest in size, is characterized as a differentiated body with a crust, a thick solid mantle, and an interior region which may be partially molten. It is presently cooling rapidly and is relatively inactive tectonically. Mercury, which probably has a large core, may have a 500-km-thick solid lithosphere and a partially molten core, if it is assumed that some heat sources exist in the core. If this is not the case, the planet's interior temperatures are everywhere below the melting curve for iron. The thermal evolution is dominated by the core separation and the high conductivity of iron that makes up the bulk of Mercury. Mars, intermediate in size, is assumed to have differentiated an Fe-FeS core. While the formation of an early crust is evident, large-scale melting and differentiation of the mantle silicates has occurred at least up until 1 b.y. ago. Present-day temperature profiles indicate moderate tectonic activity at the present time. Venus is characterized as a planet not unlike the Earth in many respects. Core formation has occurred probably during the first billion years after the formation. The present-day temperatures indicate a partially molten upper mantle overlain by a 100-km-thick lithosphere and a molten Fe-Ni core. We can expect that today Venus may have tectonic processes similar to the Earth's.

M Nafi Toksöz↗

NASA NTRS1995

A Navier-Stokes Boundary Element Solver

Using global interpolation functions (GIF's) boundary element solutions are obtained for two-dimensional laminar flows. Two schemes are proposed for handling the convective terms. The first treats convection as a forcing function, and converts the flow equations to pseudo-Poisson equations. In the second scheme, some convective effect is incorporated into the fundamental solution used in constructing the pertinent integral equations. The lid-driven cavity flow is selected as the benchmark problem.

O Lafe↗

NASA NTRS1993

An Initial Study of Using the 34-M Antenna for Lunar Mission Support

As a part of the long-range planning of future Mars and lunar missions, a feasibility study was made of a 34-m antenna system with differentially pointed multiple beams. The performance loss mechanisms of the differentially pointed multiple-beam systems were identified and quantified. Techniques that can significantly improve the multiple-beam system performance are identified. The goal is to determine the feasibility of using the 34-m antenna to support widely separated elements associated with lunar missions.

A Cha↗

NASA NTRS1985

Battery-Charge-State Model

Charge-state model for lead/acid batteries proposed as part of effort to make equivalent of fuel gage for battery-powered vehicles. Models based on equations that approximate observable characteristics of battery electrochemistry. Uses linear equations, easier to simulate on computer, and gives smooth transitions between charge, discharge, and recuperation.

Vivian, H. C.↗

NASA NTRS2014

Analytical Model for Steady Flow through a Finite Channel with One Porous Wall with Arbitrary Variable Suction or Injection

This paper presents an exact solution of two-dimensional laminar flow through a finite length channel with one porous wall. It improves upon previous solutions by (1) satisfying the no-slip boundary condition at the channel dead end, (2) adding a turbulent term to the porous wall boundary condition, (3) allowing for arbitrary variable suction or injection across the porous wall, and (4) model validation against new cryogenic liquid hydrogen and oxygen experimental data. Of particular interest in the current work is the modeling of cryogenic propellant flow through a porous liquid acquisition device (LAD) screen and channel inside a propellant tank. First, a detailed review of the literature is presented for previously attempted solutions to channel flow with one porous wall. Next, the governing equations, boundary conditions, and model assumptions are used to derive the analytical flow solution and present general model results for pressure and velocity fields within the channel. Then, the model solution is compared with horizontal LAD channel flow data in liquid oxygen as well as vertical LAD channel flow data in an inverted outflow configuration in liquid hydrogen. Model results are used to update the static cryogenic bubble point pressure model with a dynamic bubble point term which factors in enhanced convection and cooling at the screen during propellant outflow. Convective heat transfer at the LAD screen during outflow is also quantified by comparing model and data. The new analytical flow solution with the dynamic bubble point model is shown to compare well with available cryogenic experimental data

Navier Stokes Equations↗

NASA NTRS2005

Properties and Potential of Two (Ni,Pt)Ti Alloys for Use as High-Temperature Actuator Materials

The microstructure, transformation temperatures, basic tensile properties, shape memory behavior, and work output for two (Ni,Ti)Pt high-temperature shape memory alloys have been characterized. One was a Ni 30 Pt 20 Ti 50 alloy (referred to as 20Pt) with transformation temperatures above 230 °C and the other was a Ni 20 Pt 30 Ti 50 alloy (30Pt) with transformation temperatures above 530 °C. Both materials displayed shape memory behavior and were capable of 100% (no-load) strain recovery for strain levels up to their fracture limit (3-4%) when deformed at room temperature. For the 20Pt alloy, the tensile strength, modulus, and ductility dramatically increased when the material was tested just above the austenite finish (A f ) temperature. For the 30Pt alloy, a similar change in yield behavior at temperatures above the A f was not observed. In this case the strength of the austenite phase was at best comparable and generally much weaker than the martensite phase. A ductility minimum was also observed just below the A s temperature in this alloy. As a result of these differences in tensile behavior, the two alloys performed completely different when thermally cycled under constant load. The 20Pt alloy behaved similar to conventional binary NiTi alloys with work output due to the martensite-to-austenite transformation initially increasing with applied stress. The maximum work output measured in the 20Pt alloy was nearly 9 J/cm 3 and was limited by the tensile ductility of the material. In contrast, the martensite-to-austenite transformation in the 30Pt alloy was not capable of performing work against any bias load. The reason for this behavior was traced back to its basic mechanical properties, where the yield strength of the austenite phase was similar to or lower than that of the martensite phase, depending on temperature. Hence, the recovery or transformation strain for the 30Pt alloy under load was essentially zero, resulting in zero work output.

High-Temperature Shape Memory Alloy↗

NASA NTRS2017

Modeling for Battery Prognostics

For any battery-powered vehicles (be it unmanned aerial vehicles, small passenger aircraft, or assets in exoplanetary operations) to operate at maximum efficiency and reliability, it is critical to monitor battery health as well performance and to predict end of discharge (EOD) and end of useful life (EOL). To fulfil these needs, it is important to capture the battery's inherent characteristics as well as operational knowledge in the form of models that can be used by monitoring, diagnostic, and prognostic algorithms. Several battery modeling methodologies have been developed in last few years as the understanding of underlying electrochemical mechanics has been advancing. The models can generally be classified as empirical models, electrochemical engineering models, multi-physics models, and molecular/atomist. Empirical models are based on fitting certain functions to past experimental data, without making use of any physicochemical principles. Electrical circuit equivalent models are an example of such empirical models. Electrochemical engineering models are typically continuum models that include electrochemical kinetics and transport phenomena. Each model has its advantages and disadvantages. The former type of model has the advantage of being computationally efficient, but has limited accuracy and robustness, due to the approximations used in developed model, and as a result of such approximations, cannot represent aging well. The latter type of model has the advantage of being very accurate, but is often computationally inefficient, having to solve complex sets of partial differential equations, and thus not suited well for online prognostic applications. In addition both multi-physics and atomist models are computationally expensive hence are even less suited to online application An electrochemistry-based model of Li-ion batteries has been developed, that captures crucial electrochemical processes, captures effects of aging, is computationally efficient, and is of suitable accuracy for reliable EOD prediction in a variety of operational profiles. The model can be considered an electrochemical engineering model, but unlike most such models found in the literature, certain approximations are done that allow to retain computational efficiency for online implementation of the model. Although the focus here is on Li-ion batteries, the model is quite general and can be applied to different chemistries through a change of model parameter values. Progress on model development, providing model validation results and EOD prediction results is being presented.

Prognostics↗

NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

NASA NTRS1981

Battery failure model derived from flaw theory

A previously derived failure model for battery lifetime is discussed in terms of growth rate of the flaw, distribution of flaw sizes, and number of flaws. Equations are presented for determining the failure model for a nickel cadmium battery.

Schulman, I.↗

NASA NTRSDate not supplied

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

RTV↗

NASA NTRSDate not supplied

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

silicone intumescence↗

NASA NTRSDate not supplied

Power Optimal Control Allocation for Distributed Electric Propulsion in a Series/Parallel Hybrid Powertrain

This paper describes a thrust allocation scheme that minimizes power consumption in an electrified powertrain with Distributed Electric Propulsion. It takes advantage of an observation about component efficiency maps, that efficiency is often highest at high torque/high speed conditions. This optimal approach is demonstrated to satisfy total thrust and net yaw axis torque requirements, making it suitable for utilizing differential thrust for maneuvering.

electrified aircraft propulsion↗

NASA NTRSDate not supplied

Power Optimal Control Allocation for Distributed Electric Propulsion in a Series/Parallel Hybrid Powertrain

This paper describes a thrust allocation scheme that minimizes power consumption in an electrified powertrain with Distributed Electric Propulsion. It takes advantage of an observation about component efficiency maps, that efficiency is often highest at high torque/high speed conditions. This optimal approach is demonstrated to satisfy total thrust and net yaw axis torque requirements, making it suitable for utilizing differential thrust for maneuvering.

electrified aircraft propulsion↗

NASA NTRS2026

Power Optimal Control Allocation for Distributed Electric Propulsion in a Series/Parallel Hybrid Powertrain

This paper describes a thrust allocation scheme that minimizes power consumption in an electrified powertrain with Distributed Electric Propulsion. It takes advantage of an observation about component efficiency maps, that efficiency is often highest at high torque/high speed conditions. This optimal approach is demonstrated to satisfy total thrust and net yaw axis torque requirements, making it suitable for utilizing differential thrust for maneuvering.

electrified aircraft propulsion↗

NASA NTRS2026

Benchmarking Bayesian Optimization Frameworks and Acquisition Strategies for Materials Discovery and Autonomous Laboratories

Bayesian optimization (BO) can accelerate materials discovery by guiding expensive experiments toward the most promising processing conditions. We systematically compare five BO surrogate and framework combinations (Gaussian processes in Ax, Gaussian processes and Monte-Carlo neural networks in BayBE, random forests in Lolopy, and tree-structured Parzen (TPE) estimators in Hyperopt) on three benchmarks that mimic common materials design tasks (a discrete solid-electrolyte composition space, a hybrid discrete/continuous laminate-composite design problem solved with micromechanics modeling, and the continuous Ishigami analytic function which is a standard optimization benchmark). Each BO surrogate is paired with posterior mean, probability of improvement, and expected improvement acquisition functions and run for 100 trials from randomized initial samples with uniform random search providing a control. Across five random seeds per setting, BayBE’s Gaussian-process surrogate with expected improvement consistently reached ≥95 % of the known optimum in the fewest evaluations, while Lolopy’s random forest matched or exceeded GP performance on purely categorical or mixed spaces at a higher computational cost. Posterior mean alone often stagnated at local optima, underscoring the need for exploration, whereas probability and expected improvement balanced exploration and exploitation leading to better optimization in fewer trials. Execution times ranged from milliseconds for TPE to minutes for neural-network and random-forest surrogates. These results establish baseline expectations for BO in automated materials laboratories and highlight expected improvement with Gaussian processes as a reliable first choice, with random forests offering a strong alternative when categorical variables dominate. The benchmark suite and code are released to facilitate future surrogate, acquisition, and constraint-handling research in data-driven materials optimization.

Bayesian optimization↗

NASA NTRS2005

Lunar and Planetary Science XXXVI, Part II

Some topics covered: Implications of internal fragmentation on the structure of comets; Atmospheric excitation of mars polar motion; Dunite viscosity dependence on oxygen fugacity; Cross profile and volume analysis of bahram valles on mars; Calculations of the fluxes of 10-250 kV lunar leakage gamma rays; Alluvian fans on mars; Investigating the sources of the apollo 14 high-Al mare basalts; Relationship of coronae, regional plains and rift zones on venus; and Chemical differentiation and internal structure of europa and callisto.

Lunar and Planetary Science↗

NASA NTRS1977

Deep Electromagnetic Sounding of the Moon With Lunokhod 2 Data

Results of electromagnetic sounding distinguished an outer high resistance shell about 200 km thick in the moon's structure. A preliminary petrological interpretation of the moon's layers indicated their origin as a consequence of differentiation of the initial peridotite material. Upon melting, 20% to 40% of the material melts and is removed to form a high resistance basaltic shell underlain by a layer of spinal peridotites enriched in divalent iron oxides and having a reduced resistance.

L L Van'yan↗