Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equation discovery”

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.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Learning interpretable surface elasticity properties from bulk properties via neural network equation learners

Surface elasticity is central to understanding the mechanics and stability of surfaces and interfaces. It is characterized by quantities such as surface tension, residual surface stress, and surface stiffness. However their analytical expressions are typically difficult to derive from atomistic data, and depend strongly on modeling choices. This work presents a neural network-based equation learner which combines customized activation functions and connection-based pruning to discover parsimonious, closed-form equations for surface elasticity from atomistic simulations. Applying the method to seven face-centered cubic (FCC) metals, our equation learner uncovers interpretable equations that describe both low-Miller index and high-Miller index surface properties, capturing long-tail property distributions accurately. The discovered expressions are decoupled into two components: a universal, geometry-driven orientation function, and material-specific baseline coefficients. We find that lower-order properties such as surface tension are fundamentally geometry dependent, while higher-order properties such as surface stress and elasticity show more complex geometry and material dependence. We also relate material dependent coefficients to bulk properties, forming a clear map from bulk material properties to surface elasticity. Overall, this approach demonstrates that interpretable neurosymbolic machine learning can bridge the gap between atomistic simulations and physical laws, enabling the discovery of generalizable structure–property relationships for materials science phenomena such as surface elasticity.

Equation learning↗
NASA NTRS1982

Interaction of heat production, strain rate and stress power in a plastically deforming body under tensile test

At high strain rates the heat produced by plastic deformation can give rise to a rate dependent response even if the material has rate independent constitutive equations. This effect has to be evaluated when interpreting a material test, or else it could erroneously be ascribed to viscosity. A general thermodynamic theory of tensile testing of elastic-plastic materials is given in this paper; it is valid for large strain at finite strain rates. It enables discovery of the parameters governing the thermodynamic strain rate effect, provides a method for proper interpretation of the results of the tests of dynamic plasticity, and suggests a way of planning experiments in order to detect the real contribution of viscosity.

Paglietti, A.↗
NASA NTRS1979

Dynamics of liquid films and thin jets

The theory of liquid films and thin jets as one- and two-dimensional continuums is examined. The equations of motion have led to solutions for the characteristic speeds of wave propagation for the parameters characterizing the shape. The formal analogy with a compressible fluid indicates the possibility of shock wave generation in films and jets and the formal analogy to the theory of threads and membranes leads to the discovery of some new dynamic effects. The theory is illustrated by examples.

Zak, M.↗
DOE OSTI2024

Physics-Informed Active Learning With Simultaneous Weak-Form Latent Space Dynamics Identification

The parametric greedy latent space dynamics identification (gLaSDI) framework has demonstrated promising potential for accurate and efficient modeling of high-dimensional nonlinear physical systems. However, it remains challenging to handle noisy data. Here, to enhance robustness against noise, we incorporate the weak-form estimation of nonlinear dynamics (WENDy) into gLaSDI. In the proposed weak-form gLaSDI (WgLaSDI) framework, an autoencoder and WENDy are trained simultaneously to discover intrinsic nonlinear latent-space dynamics of high-dimensional data. Compared with the standard sparse identification of nonlinear dynamics (SINDy) employed in gLaSDI, WENDy enables variance reduction and robust latent space discovery, therefore leading to more accurate and efficient reduced-order modeling. Furthermore, the greedy physics-informed active learning in WgLaSDI enables adaptive sampling of optimal training data on the fly for enhanced modeling accuracy. The effectiveness of the proposed framework is demonstrated by modeling various nonlinear dynamical problems, including viscous and inviscid Burgers' equations, time-dependent radial advection, and the Vlasov equation for plasma physics. With data that contains 5%–10% Gaussian white noise, WgLaSDI outperforms gLaSDI by orders of magnitude, achieving 1%–7% relative errors. Compared with the high-fidelity models, WgLaSDI achieves 121 to 1779x speed-up.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

The electrostatics of a dusty plasma

Voyager 1 and 2 observations of the Saturnian ring system have led to the discovery of several interesting phenomena associated with its fine dust component. The dust grains are immersed in a plasma. The present paper is concerned with the electrostatics of a dusty plasma, giving particular attention to the capacitance and charge of an individual grain in the presence of neighboring grains and the surrounding plasma. The Poisson equation and gauge considerations are discussed along with a solution of the Poisson equation for several models, taking into account impermeable grains, permeable grains, and a spherical capacitor model. A comparison of potential shapes for three models is conducted, and the effect of neighbors on charging currents is explored.

Whipple, E. C.↗
NASA NTRS1998

Radar Studies in the Solar System

We aid in study of the solar system by means of ground-based radar. We have concentrated on: (1) developing the ephemerides needed to acquire radar data at Arecibo Observatory and (2) analyzing the resultant data to: test fundamental laws of gravitation; determine the size , shape, topography, and spin vectors of the targets; and study the surface properties of these objects, through their scattering law and polarization characteristics. We are engaged in radar observations of asteroids and comets, both as systematically planned targets and as "targets of opportunity." In the course of the program, we have prepared ephemerides for about 80 asteroids and three comets, and the radar observations have been made or attempted at the Arecibo Observatory, in most cases successfully, and in some cases on more than one apparition. The results of these observations have included echo spectra for the targets and, in some cases, delay - Doppler images and measurements of the total round-trip delay to the targets. Perhaps the most dramatic of these results are the images obtained for asteroids (4179) Toutatis and 1989PB (Castalia), which were revealed to be double-lobed objects by the radar images. Besides these direct results, the radar observations have furnished information on the sizes and shapes of the targets through analysis of the Doppler width of the echoes as a function of time, and on the surface properties (such as composition, bulk density, and roughness) through analysis of the reflectivity and of the polarization state of the echoes. We have also refined the orbits of the observed asteroids as a result of the Doppler (and in some cases delay) measurements from the radar observations. Although the orbits of main-belt asteroids accessible to ground-based radar are quite well known from the available optical data, some near-Earth objects have been seen by radar very soon after their optical discovery (for example, 199OMF, just eight days after discovery). In such cases. the radar results ensure that the object in question can be anticipated and identified at the next apparition. We have also participated in radar studies of the terrestrial planets. The results of these studies have included both planetary topography profiles from the analysis of round-trip delays to points along the target Doppler equator and determinations of the target spin state. The latter is of special interest in the case of Venus, which is very close to, but not on, a multi-body spin-orbit resonance such that Venus rotates 12 times for every 8 Earth orbits and 13 Venus orbits. As a result, Venus presents nearly the same face toward Earth at each inferior conjunction. Our latest results confirm that the spin state of Venus is slightly off the resonance. The delay measurements from planetary 2 ranging have also been used in combination with other types of range data in testing general relativity with increasing accuracy. We have also been engaged in radar studies of planetary satellites. Using our ephemerides, Arecibo made radar observations of the Galilean satellites of Jupiter and of Mars' satellite Phobos during the favorable opposition seasons (1988-1992 for Jupiter and 1990 for Mars). An attempt was also made to observe Deimos, but without detecting an echo. In 1997, an attempt was made to observe Saturn's satellite Titan, using the newly upgraded Arecibo radar system for transmitting and the Goldstone radar for receiving, but no echo was detected. The study of satellites by radar is in many ways similar to that of asteroids. The results from these observations have included characterization of the surface properties from the reflectivity and polarization ratio, as well as (in the case of the large satellites of Jupiter) the variation of reflectivity with incidence angle.

Shaprio, Irwin I.↗
DOE OSTI2025

Probing a New Regime of Neutrino Self-Interactions with Astrophysical Neutrinos and the Relativistic Cosmic Neutrino Background

Neutrino self-interactions beyond the standard model have profound implications in astrophysics and cosmology. In this Letter, we study an uncharted scenario in which one of the three neutrino species has a mass smaller than the temperature of the cosmic neutrino background. This results in a relativistic component that significantly broadens the absorption feature on the astrophysical neutrino spectra, in contrast to the sharply peaked absorption expected in the extensively studied scenarios assuming a fully nonrelativistic cosmic neutrino background. By solving the Boltzmann equations for neutrino absorption and regeneration, we demonstrate that this mechanism provides novel sensitivity to sub-keV mediator masses, well below the traditional ∼1–100 MeV range. Future observations of the diffuse supernova neutrino background with Hyper-Kamiokande could probe coupling strengths down to 𝑔 ∼10 −8 , surpassing existing constraints by orders of magnitude. These findings open new directions for discoveries and offer crucial insights into the interplay between neutrinos and the dark sector.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2022

Assurance Equations: A Cost and Criticality Model for Optimizing Quality Assurance Surveillance

The cost of quality vs cost of failure correction has been a long-running topic of discussion within the Aerospace community. It leads directly to concepts of “risk tolerance”, and risk-based decision-making. It would be valuable if there was a way to compute the optimal investment in customer-executed quality assurance activities using defect significance with respect to performance objectives, the activities’ defect detection effectiveness, and the cost-penalty for late discovery of impactful defects. This optimization is particularly of interest to projects whose budget constraints significantly limit their risk management options.The cost to fix defects (i.e., failure correction) escalates as the project matures. There have been studies attempting to determine the relative cost of fixing defects discovered during various phases of a project life cycle with important implications, all of which suggest growth factors are large. The commonly referred to 1:10:100 rule represents a cost multiplier for repair/rework across the Design to Fab to Test hardware development phases. Cost premiums for QA activities also accumulate when they are treated as mandatory (due to schedule drag) or are performed later than their assigned phase.This paper describes the modeling of development phase -dependencies in the conduct of typical customer-executed quality assurance activities. Our initial modeling encompasses:• Distinct phases of the production lifecycle• Multiple kinds of Defects, each with some a-priori likelihood of being present• Each defect’s impact on performance Objectives for a type of hardware• The cost and efficacy of assurance techniques at detecting such Defects• The costs of fixing those Defects detected in a given phase of the production lifecycleThe model captures assurance activities’ abilities to Detect defects. Upon detection it is assumed that the Defect is immediately fixed. Defects that “escape” detection by some activity may thereafter be detected by a later activity, but by then the cost of fixing the Defect may have escalated. Defects are related to the performance Objectives they would detract from, were those Defects to remain present in the operating system.We have constructed and are exploring, a model that relates the importance of hardware system elements to mission objectives, the impact of types of Defects on those hardware types, the cost of customer-executed assurance activities (i.e., supplier controls) and their effectiveness towards reducing an impactful quality escape, and the cost of Defect correction across production phase. We describe the approach taken to select the key model aspects, why they are relevant to our NASA mission, and our efforts to populate it with relevant and contemporary data. We use a notional example to illustrate model design and function.

Plante, Jeannette↗
DOE OSTI2026

Stoichiometrically-informed symbolic regression for extracting chemical reaction mechanisms from data

A data-driven computational method is introduced to extract chemical reaction mechanisms from time series chemical concentration data. It is realized through the use of dynamic symbolic regression in which a sparse analytical form for a dynamical system is discoverable from the underlying data. We specifically develop the stoichiometrically-informed symbolic regression (SISR) method to address a standing challenge in complex chemical reaction networks: given a time-series dataset of concentrations of several components, what is the mechanism and the associated rate constants? SISR finds the optimal mechanism, kinetic equations and rate constants by combining differential optimization with a genetic optimization approach that searches a symbolic space of possible reaction mechanisms. Use of SISR in several paradigmatic examples spanning linear and nonlinear reaction schemes results in excellent agreement between true and predicted mechanisms, including when the method is applied to noisy data. The advantages of a stoichiometrically-informed approach such as SISR to address reaction discovery is illustrated through comparison with the use of generic state-of-the-art data-driven approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Optimization techniques in self-similar compressible flow

We investigate the one-dimensional (1D) inviscid compressible flow equations for an ideal gas through the lens of optimization techniques. It is the case that, to our knowledge, optimization analysis applied to the so-called “linear velocity” solutions of the Euler compressible flow equations has not been previously conducted. Through both gradient-based and variational techniques, new variants of well-studied flow scenarios, i.e., self-similar, 1D, linear velocity solution class to idealized inviscid compressible flow equations, are determined, as encoded in both the kinematic and thermodynamic properties of this self-similar solution class. With the kinematics of the said solutions being driven by a self-similar “scale radius” and the thermodynamics being driven separately through the appearance of an arbitrary function, a myriad of new solution classes is possible. Acting as a guide to more realistic physical circumstances as well as discovery, it is the hope that the presented cases serve as the framework for future investigations into the intersection of self-similarity and optimization techniques. Fields of study that may find this work to be of interest include aerodynamic design, flow control, inertial confinement fusion, physics-informed neural networks, and other related areas of interest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1978

Cosmic-ray streaming perpendicular to the mean magnetic field. II - The gyrophase distribution function

The distribution function of cosmic rays streaming perpendicular to the mean magnetic field in a turbulent medium is reexamined. Urch's (1977) discovery that in quasi-linear theory, the flux is due to particles at 90 deg pitch angle is discussed and shown to be consistent with previous formulations of the theory. It is pointed out that this flux of particles at 90 deg cannot be arbitrarily set equal to zero, and hence the alternative theory which proceeds from this premise is dismissed. A further, basic inconsistency in Urch's transport equation is demonstrated, and the connection between quasi-linear theory and compound diffusion is discussed.

Forman, M. A.↗
NASA NTRS1980

A global biogeocenotical biosphere simulation

This model of the D. Forrester type, constructed in differential equations, predicts the food and mineral supply for the factors biosphere population, depending on two socio-economic factors, until about the year 2500. If contemporary rates of natural resources utilization are maintained and there is no management of the environment, food resources will begin to limit human population growth after 2200, and mineral resources will after 2300. A decrease in the biosphere pollution, increase in effective agricultural production, and discovery of new energy sources may forestall or completely avert the onset of a crisis situation. Conservation measures, according to the model, are to a considerable extent realizable only if carried out simultaneously in both areas.

Moiseyev, N. N.↗
NASA NTRS1988

Chemseal 3808-A2 penetration into small leak path

A possible fix to a leak in the oxidizer system of the Space Shuttle Discovery's attitude control system was proposed by MSFC. This fix involved the passing of a shuttlecock past the leaking Dynaflow fitting and sealing the vent tube containing the fitting with Chemseal 3808-A2. The question of whether the Chemseal 3808-A2 can flow into the leak path and provide a better seal was addressed analytically and by experiment to verify the analytical formula used. The results show that the equations are applicable and that the Chemseal will flow into the expected leak path and seal.

Carruth, M. R., Jr.↗
DOE OSTI2025

Robust and tunable oxide nanoscrolls for solar-driven H 2 generation and storage

Hydrogen gas is a promising alternative to fossil fuels due to its high energy output and environmentally safe byproducts. Various morphologies of photocatalytic materials have been explored for high-efficiency H 2 production, for instance, quasi-1D nanoscroll structures that provide a larger surface-to-volume ratio. Recently, we predicted layer-by-layer formation of stable oxide nanoscrolls directly from dichalcogenide precursors, eliminating the need for costly formation of two-dimensional oxides for a roll-up synthesis of nanoscrolls. Here, in this study, we evaluate the suitability of those oxide nanoscroll materials—MoO 3 , WO 3 , PdO 2 , HfO 2 , and GeO 2 —for solar-driven photocatalytic H 2 production and storage. Using excited state theory coupled with Bethe–Salpeter equation simulations, we discern their electronic and optical properties as a function of interlayer scroll spacing and find them to be highly conducive for solar-driven photocatalysis. Additionally, using ab initio molecular dynamics simulations, we show that they are also suitable for H 2 storage as the nanoscrolls exhibit an effective trapping of hydrogen, even in the presence of defects and vacancies in the oxides. This work thus demonstrates the discovery of robust and tunable oxide nanoscrolls as materials for advancing solar-driven hydrogen technologies.

Chemical compounds↗
NASA NTRS1981

Discovery and observation of BY Draconis variables

The discovery of BY Draconis variables was efficiently accomplished by a spectroscopic survey of dK-M stars for weak H-alpha emissions, using 1-2 A resolution. The four BY Dra variables discovered are all spectroscopic binaries with P values lower than about 10 d, in light of which, it is noted that the onset of high surface activity and appearance of H-alpha emission occur sharply at v(equator) of approximately 5 km/sec. At v(equator) of about 3 km/sec, dK-M stars have low levels of surface activity. It is found that while there is a range of Ca II emission strength, and only the strongest emitters of this line are BY Dra and/or flare stars, the H-alpha feature changes abruptly to an emission feature signaling the onset of flaring and/or the BY Dra syndrome. An increase of the rotation rate above v(equator) 5 km/sec does not appear to increase the level of surface activity.

Bopp, B. W.↗
NASA NTRS1985

Energy partitioning in the Io plasma torus

The Io plasma torus, which emits substantial amounts of EUV radiation from multiply ionized oxygen and sulfur, represents one of the significant discoveries of the Voyager 1 mission. The present study is concerned with the effects of various physical processes on the density and energy partitioning of the Io plasma torus. For the ions the dominant effects, which have different relative importance in different parameter regimes, are electron impact ionization, charge exchange, electron-ion collisions, and confinement time loss. Attention is given to a physical model, the solution of the quasi-linear equations, the obtained results, and an interpretation of observations in relation to the model results.

Smith, R. A.↗
NASA NTRS2005

Alluvial Fans on Mars

Moore and Howard [1] reported the discovery of large alluvial fans in craters on Mars. Their initial survey from 0-30 S found that these fans clustered in three distinct regions and occurred at around the +1 km MOLA defined Mars datum. However, due to incomplete image coverage, Moore and Howard [1]could not conduct a comprehensive survey. They also recognized, though did not quantitatively address, gravity scaling issues. Here, we briefly discuss the identification of alluvial fans on Mars, then consider the general equations governing the deposition of alluvial fans and hypothesize a method for learning about grain size in alluvial fans on Mars.

Kraal, E. R.↗
DOE OSTI2024

Machine-assisted discovery of integrable symplectic mappings

Integrable systems possess a hidden symmetry associated with the existence of conserved quantities known as integrals of motion. These systems play an important role in understanding general dynamics in accelerators and have potential for future designs. This work will cover two automated methods for finding integrable symplectic maps of the plane. The first algorithm is based on the observation that the evolution of an integrable system in phase space is confined to a lower-dimensional submanifold of a specific type. The second algorithm relies on an analysis of dynamical variables. Both methods rediscover some of the famous McMillan-Suris integrable mappings and ultra-discrete Painlev\'e equations. Over 100 new integrable families are presented and analyzed, some of which are isolated in the space of parameters, while others are families with one parameter (or the ratio of parameters) being either continuous or discrete. In addition, the newly discovered maps are related to a general 2D symplectic map through the use of discrete perturbation theory. A method is proposed for constructing smooth near-integrable dynamical systems based on mappings with polygon invariants.

43 PARTICLE ACCELERATORS↗