Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equation recovery”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

UAS Well Clear Recovery Against Non-Cooperative Intruders Using Vertical Maneuvers

This paper documents a study that drove the development of a mathematical expression in the detect-and-avoid (DAA) minimum operational performance standards (MOPS) for unmanned aircraft systems (UAS). This equation describes the conditions under which vertical maneuver guidance should be provided during recovery of DAA well clear separation with a non-cooperative VFR aircraft. Although the original hypothesis was that vertical maneuvers for DAA well clear recovery should only be offered when sensor vertical rate errors are small, this paper suggests that UAS climb and descent performance should be considered-in addition to sensor errors for vertical position and vertical rate-when determining whether to offer vertical guidance. A fast-time simulation study involving 108,000 encounters between a UAS and a non-cooperative visual-flight-rules aircraft was conducted. Results are presented showing that, when vertical maneuver guidance for DAA well clear recovery was suppressed, the minimum vertical separation increased by roughly 50 feet (or horizontal separation by 500 to 800 feet). However, the percentage of encounters that had a risk of collision when performing vertical well clear recovery maneuvers was reduced as UAS vertical rate performance increased and sensor vertical rate errors decreased. A class of encounter is identified for which vertical-rate error had a large effect on the efficacy of horizontal maneuvers due to the difficulty of making the correct left/right turn decision: crossing conflict with intruder changing altitude. Overall, these results support logic that would allow vertical maneuvers when UAS vertical performance is sufficient to avoid the intruder, based on the intruder's estimated vertical position and vertical rate, as well as the vertical rate error of the UAS' sensor.

UAS in the NAS↗
DOE OSTI2022

A modified Susceptible-Infected-Recovered model for observed under-reported incidence data

Fitting Susceptible-Infected-Recovered (SIR) models to incidence data is problematic when not all infected individuals are reported. Assuming an underlying SIR model with general but known distribution for the time to recovery, this paper derives the implied differential-integral equations for observed incidence data when a fixed fraction of newly infected individuals are not observed. The parameters of the resulting system of differential equations are identifiable. Using these differential equations, we develop a stochastic model for the conditional distribution of current disease incidence given the entire past history of reported cases. We estimate the model parameters using Bayesian Markov Chain Monte-Carlo sampling of the posterior distribution. We use our model to estimate the transmission rate and fraction of asymptomatic individuals for the current Coronavirus 2019 outbreak in eight American Countries: the United States of America, Brazil, Mexico, Argentina, Chile, Colombia, Peru, and Panama, from January 2020 to May 2021. Our analysis reveals that the fraction of reported cases varies across all countries. For example, the reported incidence fraction for the United States of America varies from 0.3 to 0.6, while for Brazil it varies from 0.2 to 0.4.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2017

UAS Well Clear Recovery Against Non-Cooperative Intruders Using Vertical Maneuvers

This paper documents a study that drove the development of a mathematical expression in the minimum operational performance standards (MOPS) of detect-and-avoid (DAA) systems for unmanned aircraft systems (UAS). This equation describes the conditions under which vertical maneuver guidance could be provided during recovery of well clear separation with a non-cooperative VFR aircraft in addition to horizontal maneuver guidance. Although suppressing vertical maneuver guidance in these situations increased the minimum horizontal separation from 500 to 800 feet, the maximum severity of loss of well clear increased in about 35 of the encounters compared to when a vertical maneuver was preferred and allowed. Additionally, analysis of individual cases led to the identification of a class of encounter where vertical rate error had a large effect on horizontal maneuvers due to the difficulty of making the correct left-right turn decision: crossing conflict with intruder changing altitude. These results supported allowing vertical maneuvers when UAS vertical performance exceeds the relative vertical position and velocity accuracy of the DAA tracker given the current velocity of the UAS and the relative vertical position and velocity estimated by the DAA tracker. Looking ahead, these results indicate a need to improve guidance algorithms by utilizing maneuver stability and near mid-air collision risk when determining maneuver guidance to regain well clear separation.

detect and avoid↗
NASA NTRS2015

Cost Scaling of a Real-World Exhaust Waste Heat Recovery Thermoelectric Generator: A Deeper Dive

Cost is equally important to power density or efficiency for the adoption of waste heat recovery thermoelectric generators (TEG) in many transportation and industrial energy recovery applications. In many cases the system design that minimizes cost (e.g., the $/W value) can be very different than the design that maximizes the system's efficiency or power density, and it is important to understand the relationship between those designs to optimize TEG performance-cost compromises. Expanding on recent cost analysis work and using more detailed system modeling, an enhanced cost scaling analysis of a waste heat recovery thermoelectric generator with more detailed, coupled treatment of the heat exchangers has been performed. In this analysis, the effect of the heat lost to the environment and updated relationships between the hot-side and cold-side conductances that maximize power output are considered. This coupled thermal and thermoelectric treatment of the exhaust waste heat recovery thermoelectric generator yields modified cost scaling and design optimization equations, which are now strongly dependent on the heat leakage fraction, exhaust mass flow rate, and heat exchanger effectiveness. This work shows that heat exchanger costs most often dominate the overall TE system costs, that it is extremely difficult to escape this regime, and in order to achieve TE system costs of $1/W it is necessary to achieve heat exchanger costs of $1/(W/K). Minimum TE system costs per watt generally coincide with maximum power points, but Preferred TE Design Regimes are identified where there is little cost penalty for moving into regions of higher efficiency and slightly lower power outputs. These regimes are closely tied to previously-identified low cost design regimes. This work shows that the optimum fill factor Fopt minimizing system costs decreases as heat losses increase, and increases as exhaust mass flow rate and heat exchanger effectiveness increase. These findings have profound implications on the design and operation of various thermoelectric (TE) waste heat 3 recovery systems. This work highlights the importance of heat exchanger costs on the overall TEG system costs, quantifies the possible TEG performance-cost domain space based on heat exchanger effects, and provides a focus for future system research and development efforts.

cost analysis↗
DOE OSTI2020

MBMS1.0: An Open-Source Code for Modeling and Simulation of Membrane-Based Dehumidification and Energy Recovery

Membrane-based dehumidification is currently being considered as a promising solution for the building application due to its low cost and very limited energy consumption. Developing a simple and efficient open-source code simulation tool is important for boosting the optimization and evaluation of such device in HVAC community. This paper reports a first-order physics based model which accounts for the fundamental heat and mass transfer of humid-air vapor at feed side to flow stream at permeate side. The current model comprises two membrane mass transfer submodels (i.e. microstructure model and performance map model); and it adopts a segment-by-segment methodology for discretizing heat and mass transfer governing equations. The model is capable of simulating both dehumidifiers and energy recovery ventilators with parallel-flow cross-flow, and counter-flow configurations. The model was validated with the measurements at appropriate device. The practices in dehumidification and energy recovery exchangers are also discussed. The model and open-source codes are expected to become a solid fundament for developing a more comprehensive and accurate membrane-based dehumidification in the future.

Gao, Zhiming↗
DOE OSTI2026

Discovery of Probabilistic Dirichlet-to-Neumann Maps on Graphs

Dirichlet-to-Neumann maps enable the coupling of multiphysics simulations across computational subdomains by ensuring continuity of state variables and fluxes at artificial interfaces. We present a novel method for learning Dirichlet-to-Neumann maps on graphs using Gaussian processes, specifically for problems where the data obey a conservation law arising from an underlying partial differential equation. Our approach combines discrete exterior calculus and nonlinear optimal recovery to infer relationships between vertex and edge values. This framework yields data-driven predictions with uncertainty quantification across the entire graph, even when observations are limited to a subset of vertices and edges. By minimizing the reproducing kernel Hilbert space norm while penalizing kernel complexity through maximum likelihood estimation, our method ensures that the resulting surrogate strictly enforces conservation laws without overfitting. We demonstrate our method on two representative applications: subsurface flow in fracture networks and arterial blood flow. Finally, the results demonstrate that the method maintains high accuracy and well-calibrated uncertainty estimates even under severe data scarcity, highlighting its potential for scientific applications where limited data and reliable uncertainty quantification are critical.

Dirichlet-to-Neumann map↗
NASA NTRS1980

Full-coverage film cooling. I - Three-dimensional measurements of turbulence structure. II - Prediction of the recovery-region hydrodynamics

Hydrodynamic measurements of turbulence structure were performed with a triaxial hot wire in the full coverage and the recovery regions following an array of injection holes under isothermal conditions at ambient temperature and pressure for blowing ratios of 0.9 and 0.4. High levels of turbulence kinetic energy (TKE) were determined for low blowing, and low TKE levels were found for the high blowing levels; in the recovery region, the flow can be represented by a model with an outer boundary layer and a 2-dimensional inner boundary layer. Recovery region hydrodynamics can be modelled by considering that a new boundary layer started to grow immediately after the end of blowing; the Prandtl mixing length distributions calculated from the values of mean velocity and turbulent shear stresses were consistent with the presence of a dual boundary layer structure in the recovery region. The program used here contains a one-equation model of turbulence, using turbulence kinetic energy with an algebraic mixing length; this 2-dimensional, finite difference program can predict the mean velocity and turbulence kinetic energy profiles based on initial values, boundary conditions, and a closure condition.

Yavuzkurt, S.↗
NASA NTRS1983

A theory of rotating stall of multistage axial compressors

A theoretical analysis was made of rotating stall in axial compressors of many stages, finding conditions for a permanent, straight-through traveling disturbance, with the steady compressor characteristic assumed known, and with simple lag processes ascribed to the flows in the inlet, blade passages, and exit regions. For weak disturbances, predicted stall propagation speeds agree well with experimental results. For a locally-parabolic compressor characteristic, an exact nonlinear solution is found and discussed. For deep stall, the stall-zone boundary is most abrupt at the trailing edge, as expected. When a complete characteristic having unstalling and reverse-flow features is adopted, limit cycles governed by a Lienard's equation are found. Analysis of these cycles yields predictions of recovery from rotating stall; a relaxation oscillation is found at some limiting flow coefficient, above which no solution exists. Recovery is apparently independent of lag processes in the blade passages, but instead depends on the lags originating in the inlet and exit flows, and also on the shape of the given characteristic diagram. Small external lags and tall diagrams favor early recovery. Implications for future research are discussed.

Moore, F. K.↗
NASA NTRS1990

A viscoplastic model with application to LiF-22 percent CaF2 hypereutectic salt

A viscoplastic model for class M (metal-like behavior) materials is presented. One novel feature is its use of internal variables to change the stress exponent of creep (where n is approximately = 5) to that of natural creep (where n = 3), in accordance with experimental observations. Another feature is the introduction of a coupling in the evolution equations of the kinematic and isotropic internal variables, making thermal recovery of the kinematic variable implicit. These features enable the viscoplastic model to reduce to that of steady-state creep in closed form. In addition, the hardening parameters associated with the two internal state variables (one scalar-valued, the other tensor-valued) are considered to be functions of state, instead of being taken as constant-valued. This feature enables each internal variable to represent a much wider spectrum of internal states for the material. The model is applied to a LiF-22 percent CaF2 hypereutectic salt, which is being considered as a thermal energy storage material for space-based solar dynamic power systems.

Freed, A. D.↗
NASA NTRS1990

A fully Sinc-Galerkin method for Euler-Bernoulli beam models

A fully Sinc-Galerkin method in both space and time is presented for fourth-order time-dependent partial differential equations with fixed and cantilever boundary conditions. The Sinc discretizations for the second-order temporal problem and the fourth-order spatial problems are presented. Alternate formulations for variable parameter fourth-order problems are given which prove to be especially useful when applying the forward techniques to parameter recovery problems. The discrete system which corresponds to the time-dependent partial differential equations of interest are then formulated. Computational issues are discussed and a robust and efficient algorithm for solving the resulting matrix system is outlined. Numerical results which highlight the method are given for problems with both analytic and singular solutions as well as fixed and cantilever boundary conditions.

Smith, R. C.↗
DOE OSTI2022

Dynamic properties of FeCrMnNi, a high entropy alloy

The goal of this paper was to assess the dynamic properties of a new class of materials, High entropy alloys (HEAs). Specifically, gas gun experiments coupled with recovery were performed on a specific HEA, FeCrMnNi, to measure its equation of state (EOS) and its spall strength, which is a measure of the stress required to nucleate voids under shock loading conditions. While there has been a plethora of work performed to assess the mechanical properties of these HEAs under uniaxial stress conditions as a function of strain rate, investigation of its properties in dynamic extremes remains rare. The current work fills this gap in knowledge for this novel class of materials. Our results show that the measured EOS for this material at one velocity was in reasonable agreement with an estimated Hugoniot. Furthermore, while the spall strength for this material was found to be ~1.9 GPa, with some variation based on sample location, the main failure mode was brittle. This brittle failure leading to formation of cracks in the material is different than the ductile failure observed in pure iron and its alloys.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Viscoplasticity: A thermodynamic formulation

A thermodynamic foundation using the concept of internal state variables is given for a general theory of viscoplasticity, as it applies to initially isotropic materials. Three fundamental internal state variables are admitted. They are: a tensor valued back stress for kinematic effects, and the scalar valued drag and yield strengths for isotropic effects. All three are considered to phenomenologically evolve according to competitive processes between strain hardening, strain induced dynamic recovery, and time induced static recovery. Within this phenomenological framework, a thermodynamically admissible set of evolution equations is put forth. This theory allows each of the three fundamental internal variables to be composed as a sum of independently evolving constituents.

Freed, A. D.↗
NASA NTRS1991

A viscoplastic theory with thermodynamic considerations

A thermodynamic foundation using the concept of internal state variables is given for a general theory of viscoplasticity for initially isotropic materials. Three fundamental, internal, state variables are admitted; they are: a tensorial back stress for kinematic effects, and scalar drag and yield strengths for isotropic effects. All three are considered to evolve phenomenologically according to competitive processes between strain hardening, deformation-induced dynamic recovery, and thermally induced static recovery. Within this phenomenological framework, a thermodynamically admissible set of evolution equations is proposed. The theory allows each of the three internal variables to be composed as a sum of independently evolving constituents. The evolution of internal state can also include terms that vary linearly with the external variable rates, whose presence affects the energy dissipation properties of a material.

Freed, A. D.↗
NASA NTRS1991

A viscoplastic theory with thermodynamic considerations

A thermodynamic foundation using the concept of internal state variables is given for a general theory of viscoplasticity for initially isotropic materials. Three, fundamental, internal, state variables are admitted; they are: a tensorial back stress for kinematic effects, and scalar drag and yield strengths for isotropic effects. All three are considered to evolve phenomenologically according to competitive processes between strain hardening, deformation induced dynamic recovery, and thermally induced static recovery. Within this phenomenological framework, a thermodynamically admissible set of evolution equations is proposed. The theory allows each of the three internal variables to be composed as a sum of independently evolving constituents. The evolution of internal state can also include terms that vary linearly with the external variable rates, whose presence affects the energy dissipation properties of a material.

Freed, A. D.↗
NASA NTRS1993

Computational analysis of methods for reduction of induced drag

The purpose of this effort was to perform a computational flow analysis of a design concept centered around induced drag reduction and tip-vortex energy recovery. The flow model solves the unsteady three-dimensional Euler equations, discretized as a finite-volume method, utilizing a high-resolution approximate Riemann solver for cell interface flux definitions. The numerical scheme is an approximately-factored block LU implicit Newton iterative-refinement method. Multiblock domain decomposition is used to partition the field into an ordered arrangement of blocks. Three configurations are analyzed: a baseline fuselage-wing, a fuselage-wing-nacelle, and a fuselage-wing-nacelle-propfan. Aerodynamic force coefficients, propfan performance coefficients, and flowfield maps are used to qualitatively access design efficacy. Where appropriate, comparisons are made with available experimental data.

Janus, J. M.↗
NASA NTRS1994

A three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawn and dusk sides of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H+ fluxes at tens of keV, which are always over-estimated. A newly-invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm-time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, j(sub o)(1+Ay(exp n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗
NASA NTRS1995

Three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawnside and duskside of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H(+) fluxes at tens of keV, which are always overestimated. A newly invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, J(sub o)(1 + Ay(sup n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗