Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional modeling”

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 919 records · Page 51

Chemically reacting supersonic flow calculation using an assumed PDF model

This work is motivated by the need to develop accurate models for chemically reacting compressible turbulent flow fields that are present in a typical supersonic combustion ramjet (SCRAMJET) engine. In this paper the development of a new assumed probability density function (PDF) reaction model for supersonic turbulent diffusion flames and its implementation into an efficient Navier-Stokes solver are discussed. The application of this model to a supersonic hydrogen-air flame will be considered.

Farshchi, M.↗

Nonlinear feedback control for high alpha flight

Analytical aerodynamic models are derived from a high alpha 6 DOF wind tunnel model. One detail model requires some interpolation between nonlinear functions of alpha. One analytical model requires no interpolation and as such is a completely continuous model. Flight path optimization is conducted on the basic maneuvers: half-loop, 90 degree pitch-up, and level turn. The optimal control analysis uses the derived analytical model in the equations of motion and is based on both moment and force equations. The maximum principle solution for the half-loop is poststall trajectory performing the half-loop in 13.6 seconds. The agility induced by thrust vectoring capability provided a minimum effect on reducing the maneuver time. By means of thrust vectoring control the 90 degrees pitch-up maneuver can be executed in a small place over a short time interval. The agility capability of thrust vectoring is quite beneficial for pitch-up maneuvers. The level turn results are based currently on only outer layer solutions of singular perturbation. Poststall solutions provide high turn rates but generate higher losses of energy than that of classical sustained solutions.

Stalford, Harold↗

Development of a calibrated software reliability model for flight and supporting ground software for avionic systems

The object of this project was to develop and calibrate quantitative models for predicting the quality of software. Reliable flight and supporting ground software is a highly important factor in the successful operation of the space shuttle program. The models used in the present study consisted of SMERFS (Statistical Modeling and Estimation of Reliability Functions for Software). There are ten models in SMERFS. For a first run, the results obtained in modeling the cumulative number of failures versus execution time showed fairly good results for our data. Plots of cumulative software failures versus calendar weeks were made and the model results were compared with the historical data on the same graph. If the model agrees with actual historical behavior for a set of data then there is confidence in future predictions for this data. Considering the quality of the data, the models have given some significant results, even at this early stage. With better care in data collection, data analysis, recording of the fixing of failures and CPU execution times, the models should prove extremely helpful in making predictions regarding the future pattern of failures, including an estimate of the number of errors remaining in the software and the additional testing time required for the software quality to reach acceptable levels. It appears that there is no one 'best' model for all cases. It is for this reason that the aim of this project was to test several models. One of the recommendations resulting from this study is that great care must be taken in the collection of data. When using a model, the data should satisfy the model assumptions.

Lawrence, Stella↗

PDF turbulence modeling and DNS

The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in probability density function (pdf). A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models. The effect of Coriolis forces on compressible homogeneous turbulence is studied using direct numerical simulation (DNS). The numerical method used in this study is an eight order compact difference scheme. Contrary to the conclusions reached by previous DNS studies on incompressible isotropic turbulence, the present results show that the Coriolis force increases the dissipation rate of turbulent kinetic energy, and that anisotropy develops as the Coriolis force increases. The Taylor-Proudman theory does apply since the derivatives in the direction of the rotation axis vanishes rapidly. A closer analysis reveals that the dissipation rate of the incompressible component of the turbulent kinetic energy indeed decreases with a higher rotation rate, consistent with incompressible flow simulations (Bardina), while the dissipation rate of the compressible part increases; the net gain is positive. Inertial waves are observed in the simulation results.

Hsu, A. T.↗

Models of Plumes: Their Flow, Their Geometric Spreading, and Their Mixing with Interplume Flow

There are two types of plume flow models: (1) 1D models using ad hoc spreading functions, f(r); (2) MagnetoHydroDynamics (MHD) models. 1D models can be multifluid, time dependent, and incorporate very general descriptions of the energetics. They confirm empirical results that plume flow is slow relative to requirements for high speed wind. But, no published 1 D model incorporates the rapid local spreading at the base (fl(r)) which has an important effect on mass flux. The one published MHD model is isothermal, but confirms that if b=8*pi*p/absolute value(B)2<<l then the field is nearly potential below -70,000 km. Building on the MHD result, we apply a two scale approximation to calculate fl(r). We also compute the global spreading (fg(r)) out to 5.0 RSUN imposed by coronal hole geometry. Global MHD models provide a potent method of calculating fg(r). Unambiguous plume signatures have not yet been found in the solar wind. This is probably due to strong mixing of plume and interplume flows near the Sun. We describe a physical source for strong mixing due to the observed flows being unstable to shear instabilities that lead to rapid disruption.

Suess, Steven T.↗

Vertical Structure and Optical Properties of Titans Aerosols from Radiance Measurements Made Inside and Outside the Atmosphere

Prompted by the detection of stratospheric cloud layers by Cassini's Composite Infrared Spectrometer (CIRS; see Anderson, C.M., Samuelson, R.E. [2011]. Icarus 212, 762-778), we have re-examined the observations made by the Descent Imager/Spectral Radiometer (DISR) in the atmosphere of Titan together with two constraints from measurements made outside the atmosphere. No evidence of thin layers (<1 km) in the DISR image data sets is seen beyond the three previously reported layers at 21 km, 11 km, and 7 km by Karkoschka and Tomasko (Karkoschka, E., Tomasko, M.G. [2009]. Icarus 199, 442-448). On the other hand, there is evidence of a thicker layer centered at about 55 km. A rise in radiance gradients in the Downward-Looking Visible Spectrometer (DLVS) data below 55 km indicates an increase in the volume extinction coefficient near this altitude. To fit the geometric albedo measured from outside the atmosphere the decrease in the single scattering albedo of Titan's aerosols at high altitudes, noted in earlier studies of DISR data, must continue to much higher altitudes. The altitude of Titan's limb as a function of wavelength requires that the scale height of the aerosols decrease with altitude from the 65 km value seen in the DISR observations below 140 km to the 45 km value at higher altitudes. We compared the variation of radiance with nadir angle observed in the DISR images to improve our aerosol model. Our new aerosol model fits the altitude and wavelength variations of the observations at small and intermediate nadir angles but not for large nadir angles, indicating an effect that is not reproduced by our radiative transfer model. The volume extinction profiles are modeled by continuous functions except near the enhancement level near 55 km altitude. The wavelength dependence of the extinction optical depth is similar to earlier results at wavelengths from 500 to 700 nm, but is smaller at shorter wavelengths and larger toward longer wavelengths. A Hapke-like model is used for the ground reflectivity, and the variation of the Hapke single scattering albedo with wavelength is given. Fits to the visible spectrometers looking upward and downward are achieved except in the methane bands longward of 720 nm. This is possibly due to uncertainties in extrapolation of laboratory measurements from 1 km-am paths to much longer paths at lower pressures. It could also be due to changes in the single scattering phase functions at low altitudes, which strongly affect the path length through methane that the photons travel. We demonstrate the effects on the model fits by varying each model parameter individually in order to illustrate the sensitivity of our determination of each model parameter.

Titan↗

Lunar Phase Function at 1064 Nm from Lunar Orbiter Laser Altimeter Passive and Active Radiometry

We present initial calibration and results of passive radiometry collected by the Lunar Orbiter Laser Altimeter onboard the Lunar Reconnaissance Orbiter over the course of 12 months. After correcting for time- and temperature-dependent dark noise and detector responsivity variations, the LOLA passive radiometry measurements are brought onto the absolute radiance scale of the SELENE Spectral Profiler. The resulting photometric precision is estimated to be 5%. We leverage the unique ability of LOLA to measure normal albedo to explore the 1064 nm phase function's dependence on various geologic parameters. On a global scale, we find that iron abundance and optical maturity (quantified by FeO and OMAT) are the dominant controlling parameters. Titanium abundance (TiO2), surface roughness on decimeter to decameter scales, and soil thermo- physical properties have a smaller effect, but the latter two are correlated with OMAT, indicating that exposure age is the driving force behind their effects in a globally-averaged sense. The phase function also exhibits a dependence on surface slope at approximately 300 m baselines, possibly the result of mass wasting exposing immature material and/or less space weathering due to reduced sky visibility. Modeling the photometric function in the Hapke framework, we find that, relative to the highlands, the maria exhibit decreased backscattering, a smaller opposition effect (OE) width, and a smaller OE amplitude. Immature highlands regolith has a higher backscattering fraction and a larger OE width compared to mature highlands regolith. Within the maria, the backscattering fraction and OE width show little dependence on TiO2 and OMAT. Variations in the phase function shape at large phase angles are observed in and around the Copernican-aged Jackson crater, including its dark halo, a putative impact melt deposit. Finally, the phase function of the Reiner Gamma Formation behaves more optically immature than is typical for its composition and OMAT, suggesting the visible-to-near-infrared spectrum and phase function respond differently to the unusual regolith evolution and properties at this location. The phase difference map revealed additional geologically-influenced variations in the phase function's shape. In particular, variations were observed associated with the dark halo around Jackson crater, the impact ejecta of Copernicus and Giordano Bruno, and the Reiner Gamma Formation. For the latter, we 915 found that the phase function behaves more optically immature than the global phase function for its composition and OMAT, suggesting a difference in how the visible-to-near-IR spectrum and phase function respond to the unusual regolith evolution and properties at this location.

moon↗

Lunar Phase Function at 1064 Nm from Lunar Orbiter Laser Altimeter Passive and Active Radiometry

We present initial calibration and results of passive radiometry collected by the Lunar Orbiter Laser Al- timeter onboard the Lunar Reconnaissance Orbiter over the course of 12 months. After correcting for time- and temperature-dependent dark noise and detector responsivity variations, the LOLA passive radiometry measurements are brought onto the absolute radiance scale of the SELENE Spectral Profiler. The resulting photometric precision is estimated to be approximately 5%. We leverage the unique ability of LOLA to measure normal albedo to explore the 1064 nm phase function's dependence on various geologic parameters. On a global scale, we find that iron abundance and optical maturity (quantified by FeO and OMAT) are the dominant controlling parameters. Titanium abundance (TiO2 ), surface roughness on decimeter to decameter scales, and soil thermophysical properties have a smaller effect, but the latter two are correlated with OMAT, indicating that exposure age is the driving force behind their effects in a globally-averaged sense. The phase function also exhibits a dependence on surface slope at approximately 300 m baselines, possibly the result of mass wasting exposing immature material and/or less space weathering due to reduced sky visibility. Modeling the photometric function in the Hapke framework, we find that, relative to the highlands, the maria exhibit decreased backscattering, a smaller opposition effect (OE) width, and a smaller OE amplitude. Immature highlands regolith has a higher backscattering fraction and a larger OE width compared to mature highlands regolith. Within the maria, the backscattering fraction and OE width show little dependence on TiO2 and OMAT. Variations in the phase function shape at large phase angles are observed in and around the Copernican-aged Jackson crater, including its dark halo, a putative impact melt deposit. Finally, the phase function of the Reiner Gamma Formation behaves more optically immature than is typical for its composition and OMAT, suggesting the visible-to-near-infrared spectrum and phase function respond differently to the unusual regolith evolution and properties at this location.

LOLA↗

Development of microgravity, full body functional reach envelope using 3-D computer graphic models and virtual reality technology

In microgravity conditions mobility is greatly enhanced and body stability is difficult to achieve. Because of these difficulties, optimum placement and accessibility of objects and controls can be critical to required tasks on board shuttle flights or on the proposed space station. Anthropometric measurement of the maximum reach of occupants of a microgravity environment provide knowledge about maximum functional placement for tasking situations. Calculations for a full body, functional reach envelope for microgravity environments are imperative. To this end, three dimensional computer modeled human figures, providing a method of anthropometric measurement, were used to locate the data points that define the full body, functional reach envelope. Virtual reality technology was utilized to enable an occupant of the microgravity environment to experience movement within the reach envelope while immersed in a simulated microgravity environment.

Lindsey, Patricia F.↗

Theoretical calculations of ion acceleration in the vicinity of comet Giacobini-Zinner

Ionization of cometary neutral molecules produces ions which are picked up by the solar wind. The cometary ion pickup process for comet Giacobini-Zinner is studied in two ways: (1) with a test particle method in which trajectories are numerically calculated for several thousand ions whose initial locations were chosen randomly with probability proportional to the neutral density and (2) with a quasi-linear diffusion model. The cometary ion distribution function was calculated with the test particle model at several locations upstream of the bow shock and for several types of magnetic fluctuations (or waves). These waves were allowed to propagate in both directions along the magnetic field at the Alfven speed. Both pitch angle scattering and energy diffusion are evident in the derived ion distributions. The monochromatic waves result in less ion acceleration than turbulent fluctuations with about the same amplitude. The calculated ion distribution functions are in reasonable agreement with the distributions measured in the vicinity of comet Giacobini-Zinner in 1985 by particle detectors on the ICE spacecraft when the ratio of power in sunward propagating Alfven waves to the power in antisunward propagating waves is assumed to lie between about 20% and 50%. However, the quasi-linear diffusion model results agree best with the measured distribution functions when the power ratio is only about 3 percent.

Mckenzie, M. L.↗

Model error estimation for large flexible spacecraft

The basic model obtained for spacecraft control system design applications is used to develop the control and estimation algorithms which constitute integral elements of the control system. The performance of this system is limited on account of modeling errors. The presence of such errors is inevitable in connection with truncated dynamics, parameter uncertainties, neglected nonlinearities, and external disturbances. In many cases, an approximate knowledge of the model errors can only be established by means of an estimation procedure. The present investigation has the objective to outline an approach to estimate errors inherent in a reduced-order model and to illustrate its application to modeling problems arising in spacecraft control system design. The considered procedure makes use of the principle of least-squares. In the discussed application this principle is employed to minimize a quadratic functional of the model errors.

Rodriguez, G.↗

Communication Limits Due to Photon-Detector Jitter

A theoretical and experimental study was conducted of the limit imposed by photon-detector jitter on the capacity of a pulse-position-modulated optical communication system in which the receiver operates in a photon-counting (weak-signal) regime. Photon-detector jitter is a random delay between impingement of a photon and generation of an electrical pulse by the detector. In the study, jitter statistics were computed from jitter measurements made on several photon detectors. The probability density of jitter was mathematically modeled by use of a weighted sum of Gaussian functions. Parameters of the model were adjusted to fit histograms representing the measured-jitter statistics. Likelihoods of assigning detector-output pulses to correct pulse time slots in the presence of jitter were derived and used to compute channel capacities and corresponding losses due to jitter. It was found that the loss, expressed as the ratio between the signal power needed to achieve a specified capacity in the presence of jitter and that needed to obtain the same capacity in the absence of jitter, is well approximated as a quadratic function of the standard deviation of the jitter in units of pulse-time-slot duration.

Moision, Bruce E.↗

Distributed digital signal processors for multi-body flexible structures

Multi-body flexible structures, such as those currently under investigation in spacecraft design, are large scale (high-order) dimensional systems. Controlling and filtering such structures is a computationally complex problem. This is particularly important when many sensors and actuators are located along the structure and need to be processed in real time. This report summarizes research activity focused on solving the signal processing (that is, information processing) issues of multi-body structures. A distributed architecture is developed in which single loop processors are employed for local filtering and control. By implementing such a philosophy with an embedded controller configuration, a supervising controller may be used to process global data and make global decisions as the local devices are processing local information. A hardware testbed, a position controller system for a servo motor, is employed to illustrate the capabilities of the embedded controller structure. Several filtering and control structures which can be modeled as rational functions can be implemented on the system developed in this research effort. Thus the results of the study provide a support tool for many Control/Structure Interaction (CSI) NASA testbeds such as the Evolutionary model and the nine-bay truss structure.

Lee, Gordon K. F.↗

Dynamic models for simulation of the 70-M antenna axis servos

Dynamic models for the various functional elements of the 70 m antenna axis servos are described. The model representing the digital position controller, the linear and nonlinear properties of the physical hardware, and the dynamics of the flexible antenna structure are encoded in six major function blocks. The general modular structure of the function blocks facilitates their adaptation to a variety of dynamic simulation studies. Model parameter values were calculated from component specifications and design data. A simulation using the models to predict limit cycle behavior produced results in excellent agreement with field test data from the DSS 14 70-m antenna.

Hill, R. E.↗

Local Ultrasonic Resonance Spectroscopy of Lithium Metal Batteries for Aerospace Applications

As next-generation aircraft and vehicles continue to develop, so do their associated energy demands. Lithium metal batteries are a leading candidate to fulfill this energy requirement, but these batteries are prone to internal dendrite defects that can lead to catastrophic thermal runaway events. Current battery management systems are capable of mitigating such risks, but are unable to detect such defects until thermal runaway has already begun. Various nondestructive evaluation (NDE) techniques, particularly ultrasonic NDE, can directly monitor internal battery parameters giving them the potential to detect critical defects prior to catastrophic failure. However, most of the current battery NDE research has focused on improved battery state-of-charge (SOC) and stateof- health (SOH) monitoring with little emphasis on critical defect detection. Thus, a measurement technique sensitive to subtle battery defects is needed. In addition, the complex mechanics of ultrasound in porous, thin, multilayered batteries prompt the use of physics-based simulation to guide inspections. In this work, an ultrasonic NDE technique has been developed utilizing frequency domain analysis of local battery resonances to detect the presence of battery defects. This technique is a practical extension of local ultrasonic resonance spectroscopy (LURS) – which previously required non-contact laser ultrasonics – to measurements with piezoelectric contact and immersion scan transducers. To extend the technique to work with piezoelectric transducers, ultrasonic battery measurements were compared to a sans-battery calibration measurement. Then, a linear systems deconvolution was used to eliminate the transfer functions of extraneous factors such as the transducer and electronics, leaving only the frequency-dependent battery reflection coefficient. The LURS technique was first validated on stainless steel and aluminum plates, producing reflection coefficients in line with analytical and numerical finite element modeling (FEM) results. Functioning Li-metal pouch cells were then seeded with lithium chip defects prior to LURS measurements. The presence of these defects is shown to cause a measurable shift in the battery’s through-thickness local resonances. 2D, frequency-domain poroelastic models of ultrasonic propagation in a single-cell lithium metal pouch battery were created and corroborated these findings. Thus, this work has both extended and proven the feasibility of the LURS technique in the detection of local battery defects.

Ultrasound↗

Review of Hydrologic Models for Evaluating Use of Remote Sensing Capabilities

Hydrologic models most commonly used by federal agencies for hydrologic forecasting are reviewed. Six catchment models and one snow accumulation and ablation model are reviewed. Information on the structure, parameters, states, and required inputs is presented in schematic diagrams and in tables. The primary and secondary roles of parameters and state variables with respect to their function in the models are identified. The information will be used to evaluate the usefulness of remote sensing capabilities in the operational use of hydrologic models.

Peck, E. L.↗

Modeling preparation costs for space missions by using major cost drivers

A model is suggested for making long range planning cost estimates for Deep Space Network (DSN) support of future space missions. The model is a function of major mission cost drivers, such as maintenance and operations, downlink frequency upgrade, uplink frequency upgrade, telemetry upgrade, antenna gain/noise temperature, radiometric accuracy upgrade, radio science upgrade, and very long baseline interferometry. The model is derived from actual cost data from three space missions: Voyager (Uranus), Voyager (Neptune), and Magellan. The model allows one to estimate the total cost and the cost over time of a similar future space mission. The model was back tested against the three projects and gave cost estimates that range from 17 pct. below to 19 pct. above actual mission preparation costs. The model was also compared with two other independent projects: Mariner Jupiter/Saturn (MJS later became Voyager) and Viking. The model gave total preparation cost estimates that range from 15 pct. above to 4 pct. below actual total preparation costs for MJS and Viking, respectively.

Sherif, J. S.↗

Full Gradient Solution to Adaptive Hybrid Control

This paper focuses on the adaptation mechanisms in adaptive hybrid controllers. Most adaptive hybrid controllers update two filters individually according to the filtered-reference least mean squares (FxLMS) algorithm. Because this algorithm was derived for feedforward control, it does not take into account the presence of a feedback loop in the gradient calculation. This paper provides a derivation of the proper weight vector gradient for hybrid (or feedback) controllers that takes into account the presence of feedback. In this formulation, a single weight vector is updated rather than two individually. An internal model structure is assumed for the feedback part of the controller. The full gradient is equivalent to that used in the standard FxLMS algorithm with the addition of a recursive term that is a function of the modeling error. Some simulations are provided to highlight the advantages of using the full gradient in the weight vector update rather than the approximation.

Bean, Jacob↗