Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Reduced Order 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 1,153 records · Page 64

Numerical experiments related to the summer 1980 U.S. heat wave

The NASA Goddard Laboratory for Atmospheres fourth-order GCM was used in a series of medium-range numerical forecast experiments in order to improve understanding of the severe summer 1980 heat wave over the U.S. The results show that the derived soil moisture anomalies in the summer of 1980 contributed positively to the model's simulation of the heat wave maintenance, and suggest that once a region of reduced soil moisture is established, it tends to persist and maintain warmer and drier conditions. The lower soil moisture values resulted in reduced evaporation, higher ground temperatures, increased sensible heat flux from ground to air, higher surface air temperature, lower sea-level pressure, and higher 500-mb height. The effects of North Pacific sea-surface-temperature anomalies were mostly opposite to those of the soil-moisture anomalies: enhanced northerly flow of cooler drier air, increased evaporation, lower ground and air temperature, higher sea level pressure, and lower 500 mb heights over the Great Plains.

Wolfson, N.

Integrated flight/propulsion control design for a STOVL aircraft using H-infinity control design techniques

Results are presented from an application of H(infinity) control design methodology to a centralized integrated flight/propulsion control (IFPC) system design for a supersonic short take-off and vertical landing (STOVL) fighter aircraft in transition flight. The emphasis is on formulating the H(infinity) control design problem such that the resulting controller provides robustness to modeling uncertainties and model parameter variations with flight condition. Experience gained from a preliminary H(infinity)=based IFPC design study performed earlier is used as the base to formulate the robust H(infinity) control design problem and improve the previous design. Detailed evaluation results are presented for a reduced-order controller obtained from the improved H(infinity) control design showing that the control design meets the specified nominal performance objectives as well as provides stability robustness for variations in plant system dynamics with chnages in aircraft trim speed within the transition flight envelope.

Garg, Sanjay

Hierarchical approach to reconstruct surfaces by using iteratively rectified imagery

A new approach to reconstruct the 3D surface of the object space from digital images is described. All the object points obtained by an automatic orientation procedure lead to a first approximation of the surface. Edges are computed for one image and matched to the other image by grey level correlation or least-squares matching through the scale space. To every discrete step in the scale space there exists the digital stereopair (image pyramid), the corresponding surface digital elevation model (DEM), and the warped images. The warped images in this discrete scale space representation correspond to digital orthophotos obtained from the DEMs that result from matching the image pyramid. The warped images will be used on every successive level in the image pyramid in order to reduce the foreshortening problems associated with any area-based matching method.

Schenk, Toni

Use of nonlinear identification in robust attitude and attitude rate estimation for SAMPEX

A method is described for obtaining optimal attitude estimation/identification algorithms for spacecraft lacking attitude rate measurement devices (rate gyros), and then demonstrated using actual flight data from the Solar, Anomalous, and Magnetospheric Particle Explorer (SAMPEX) spacecraft. SAMPEX does not have on-board rate sensing, and relies on sun sensors and a three-axis magnetometer for attitude determination. The absence of rate data normally reduces both the total amount of data available and the sampling density (in time) by a substantial fraction. In addition, attitude data is occasionally unavailable (for example, during sun occultation). As a result, the sensitivity of the estimates to model uncertainty and to measurement noise increases. In order to maintain accuracy in the attitude estimates, there is an increased need for accurate models of the rotational dynamics. The Minimum Model Error(MME)/Least Square Correlation(LSC) algorithm accurately identifies an improved model for SAMPEX to be used during periods of complete data loss or extreme noise. The model correction is determined by estimating only one orbit(the identification pass) just prior to the assumed data loss(the prediction pass). The MME estimator correctly predicted the states during the identification phase, but more importantly determines the necessary model correction trajectory, d(t). The LSC algorithm is then used to find this trajectory's functional form, H(x(t)). The results show significant improvement of the new corrected model's attitude estimates as compared to the original uncorrected model's estimates. The possible functional form of the correction term is limited at this point in the study to functions strictly of the estimated states. The results, however, strongly suggest that functions based on the relative position of the satellite may also be possible candidates for future consideration.

Mook, D. Joseph

The oscillations of superfluid neutron stars

The effects of superfluid hydrodynamics on the oscillations of neutron stars are investigated. The equations describing the small-amplitude pulsations of a neutron star's neutron-proton superfluid mixture are reduced to a system of three second-order equations for three scalar potentials. A variational principle is developed from which the frequencies of the modes of these superfluid oscillations may be estimated. These pulsations are studied by finding analytical solutions to the equations for simple uniform models of nonrotating neutron stars and numerical solutions for realistic models (including rotation). These solutions reveal that the lowest frequency modes are almost indistinguishable from those based on ordinary-fluid hydrodynamics. The analytical solutions also reveal the existence of a new set of modes having no ordinary-fluid counterpart.

Lindblom, Lee

Launch Vehicle Design Process Description and Training Formulation

A primary NASA priority is to reduce the cost and improve the effectiveness of launching payloads into space. As a consequence, significant improvements are being sought in the effectiveness, cost, and schedule of the launch vehicle design process. In order to provide a basis for understanding and improving the current design process, a model has been developed for this complex, interactive process, as reported in the references. This model requires further expansion in some specific design functions. Also, a training course for less-experienced engineers is needed to provide understanding of the process, to provide guidance for its effective implementation, and to provide a basis for major improvements in launch vehicle design process technology. The objective of this activity is to expand the description of the design process to include all pertinent design functions, and to develop a detailed outline of a training course on the design process for launch vehicles for use in educating engineers whose experience with the process has been minimal. Building on a previously-developed partial design process description, parallel sections have been written for the Avionics Design Function, the Materials Design Function, and the Manufacturing Design Function. Upon inclusion of these results, the total process description will be released as a NASA TP. The design function sections herein include descriptions of the design function responsibilities, interfaces, interactive processes, decisions (gates), and tasks. Associated figures include design function planes, gates, and tasks, along with other pertinent graphics. Also included is an expanded discussion of how the design process is divided, or compartmentalized, into manageable parts to achieve efficient and effective design. A detailed outline for an intensive two-day course on the launch vehicle design process has been developed herein, and is available for further expansion. The course is in an interactive lecture/workshop format to engage the participants in active learning. The course addresses the breadth and depth of the process, requirements, phases, participants, multidisciplinary aspects, tasks, critical elements,as well as providing guidance from previous lessons learned. The participants are led to develop their own understanding of the current process and how it can be improved. Included are course objectives and a session-by-session outline of course content. Also included is an initial identification of visual aid requirements.

Atherton, James

How Ionospheric Ions Populate the Magnetosphere during a Magnetic Storm

Ionospheric oxygen ions have been observed throughout the magnetosphere, from the plasma sheet to the ring current region. I t has been found that the O+ /H+ density ratio in the magnetosphere increases with geomagnetic activity and varies with storm phases. During the magnetic storm in late September to earIy October 2002, Cluster was orbiting in the plasma sheet and ring current regions. At prestorm time, Cluster observed high H+ density and low O+ density in the plasma sheet and lobes. During the storm main phase, 0+ density has increased by 10 times over the pre-storm level. Strong field-aligned beams of O+ were observed in the lobes. O+ fluxes were significantly reduced in the central plasma sheet during the storm recovery. However, 0+ was still evident on the boundaries of the plasma sheet and in the lobes. In order to interpret the Cluster observations and to understand how O+ ions populate the magnetosphere during a magnetic storm, we model the storm in early October 2002 using our global ion kinetic simulation (GIK). We use the LFN global simulation model to produce electric and magnetic fields in the outer magnetosphere, the Strangeway outflow scaling with Delcourt ion trajectories to include ionospheric outflows, and the Fok inner magnetospheric model for the plasmaspheric and ring current response to all particle populations. We find that the observed composition features are qualitatively reproduced by the simulations, with some quantitative differences that point to future improvements in the models.

Fok, Mei-Ching

Biological Response to SPE Exposures

It has long been recognized that a single solar particle event (SPE) can produce, over a short period of time, exposures on the order of LD50 for humans under normal physiological conditions. It is further recognized that recovery from injury over the period of exposure would greatly increase the chances of survival (dose rate effects) although such effects were left unquantified. In the present report we use the bioresponse model derived from a broad range of animal and human exposure data for evaluation of troop readiness in tactical nuclear warfare to evaluate the biological risk posed by the solar event of 4 August 1972. The astronaut blood forming organ (BFO) exposure in deep space would have been 2.2 Sv (1.6 Gy) in a space suit, 1.8 Sv (1.3 Gy) in an aluminum pressure vessel, and 0.7 Sv (0.5 Gy) in an equipment room compared to an X-ray mortality threshold of 1.5 Gy (assuming high dose rate). We find BFO dose rate effectiveness factors for this SPE on the order of 3 to 4, greatly reducing the mortality risks for this event. There is an approximate 3 percent chance that an even larger event may occur for which exposures could be 2-4 times higher. Assured survival of the astronaut requires added shelter shielding and a warning system for this event. The required mass of the shelter shield can be greatly reduced by using hydrogenous materials such as polymers, water, food, and other biological materials in its construction. Limitations of the current bioresponse model arise from the exposures taking place in the microgravity environment wherein the immune system is already challenged and the effective mortality threshold may be reduced by a factor of two. Such microgravity effects could greatly affect astronaut risks.

Wilson, J. W.

Mixing and non-equilibrium chemical reaction in a compressible mixing layer

The effects of compressibility, chemical reaction exothermicity, and non-equilibrium chemical modeling in a reacting plane mixing layer were investigated by means of two dimensional direct numerical simulations. The chemical reaction was irreversible and second order of the type A + B yields Products + Heat. The general governing fluid equations of a compressible reacting flow field were solved by means of high order finite difference methods. Physical effects were then determined by examining the response of the mixing layer to variation of the relevant non-dimensionalized parameters. The simulations show that increased compressibility generally results in a suppressed mixing, and consequently a reduced chemical reaction conversion rate. Reaction heat release was found to enhance mixing at the initial stages of the layer growth, but had a stabilizing effect at later times. The increased stability manifested itself in the suppression or delay of the formation of large coherent structures within the flow. Calculations were performed for a constant rate chemical kinetics model and an Arrhenius type kinetic prototype. The choice of the model was shown to have an effect on the development of the flow. The Arrhenius model caused a greater temperature increase due to reaction than the constant kinetic model. This had the same effect as increasing the exothermicity of the reaction. Localized flame quenching was also observed when the Zeldovich number was relatively large.

Steinberger, Craig J.

QuadTune version 1: a regional tuner for global atmospheric models

When a new, better-formulated physical parameterization is introduced into a global atmospheric model, aspects of the global model solutions are sometimes degraded. Then, in order to use the new global model to address science questions, there is an incentive to restore its accuracy. Oftentimes this restoration is achieved by tuning of model parameter values. Unfortunately, the retuning process is expensive because characterizing the parameter dependence requires numerous time-consuming global simulations. To reduce the cost of tuning, this manuscript introduces a “poor man's” model tuner, “QuadTune”. QuadTune carves the globe into regions and approximates the model parameter dependence through the use of an uncorrelated quadratic emulator (i.e., response surface). The simplicity of the emulator reduces the required number of global model simulations and aids explainability of tuner behavior. Tuning removes parametric error but leaves behind model structural error. Structural error manifests itself as regional residual biases, such as stubborn biases and tuning trade-offs. To visualize these residual biases, QuadTune's software includes a set of diagnostic plots. This paper illustrates the use of the plots for characterizing residual biases with an example tuning problem.

Environmental sciences

An on-line equivalent system identification scheme for adaptive control

A prime obstacle to the widespread use of adaptive control is the degradation of performance and possible instability resulting from the presence of unmodeled dynamics. The approach taken is to explicitly include the unstructured model uncertainty in the output error identification algorithm. The order of the compensator is successively increased by including identified modes. During this model building stage, heuristic rules are used to test for convergence prior to designing compensators. Additionally, the recursive identification algorithm as extended to multi-input, multi-output systems. Enhancements were also made to reduce the computational burden of an algorithm for obtaining minimal state space realizations from the inexact, multivariate transfer functions which result from the identification process. A number of potential adaptive control applications for this approach are illustrated using computer simulations. Results indicated that when speed of adaptation and plant stability are not critical, the proposed schemes converge to enhance system performance.

Sliwa, S. M.

About the coupling of turbulence closure models with averaged Navier-Stokes equations

The MacCormack implicit predictor-corrector model (1981) for numerical solution of the coupled Navier-Stokes equations for turbulent flows is extended to nonconservative multiequation turbulence models, as well as the inclusion of second-order Reynolds stress turbulence closure. A scalar effective pressure turbulent contribution to the pressure field is defined to approximate the effects of the Reynolds stress in strongly sheared flows. The Jacobian matrices of the transport equations are diagonalized to reduce the required computer memory and run time. Techniques are defined for including turbulence in the diagonalization. Application of the method is demonstrated with solutions generated for transonic nozzle flow and for the interaction between a supersonic flat plate boundary layer and a 12 deg compression-expansion ramp.

Vandromme, D.

On state estimation for an orbiting single tether system

The effects of instrumentation accuracy and configuration on estimation error are studied for the small expendable-tether deployment system (SEDS) using a continuous-discrete extended Kalman filter (CDEKF) state estimator. A twelfth order model that incorporates the rigid body modes of the tether as well as the satellite attitude dynamics is developed. Simulation results using the model and the estimator indicate that the originally planned instrumentation package could not estimate the state vector adequately. Recommendations are made and results presented that reduce the estimation error by adding instruments and increasing selected measurement accuracies.

Greene, Michael E.

Probabilistic Study of Fluid Structure Interaction

Probabilistic CFD design is needed because we are asked to do more with less. To cost effectively accomplish the design task, we need to formally quantify the effect of uncertainties (variables) in the design. Probabilistic design is one effective method to formally quantify the effect of uncertainties. Our objective is to establish a revolutionary new early design process, by developing non-deterministic physics-based probabilistic design tools, which will include all the life cycle processes. Breakthroughs will be sought in speed, accuracy, intelligence, and usability of the system. This paper is concerned with the usefulness of parametric optimization method coupled with a Navier-Stokes analysis code for the aero-thermodynamic design of turbomachinery combustor liner. The interconnection between the CFD code and NESSUS codes facilitated the coupling between the thermal profiles and structural design. We have developed new concepts for reducing the computational cost of unsteady, three-dimensional, compressible aerodynamic analyses for multistage turbomachinery flows. The flow was modeled by the three-dimensional Favre-Reynolds-averaged Navier-Stokes equations using the k-epsilon turbulence closure, which was integrated using an implicit third-order upwind solver. The methodology developed in this paper is expected to lead to the design optimization of turbomachinery blades.

Gorla, Rama S. R.

Probabilistic Unsteady Aerodynamic Analysis

Probabilistic CFD design is needed because we are asked to do more with less. To cost effectively accomplish the design task, we need to formally quantify the effect of uncertainties (variables) in the design. Probabilistic design is one effective method to formally quantify the effect of uncertainties. Our objective is to establish a revolutionary new early design process, by developing non-deterministic physics-based probabilistic design tools, which will include all the life cycle processes. This work was concerned with the usefulness of parametric optimization method coupled with a Navier-Stokes analysis code for the aero-thermodynamic design of turbomachinery combustor liner. The interconnection between the CFD code and NESSUS codes will facilitate the coupling between the thermal profiles and structural design. We have developed new concepts for reducing the computational cost of unsteady, three-dimensional, compressible aerodynamic analyses for multistage turbomachinery flows. The flow was modeled by the three-dimensional Favre-Reynolds-averaged Navier-Stokes equations using the k-E turbulence closure, which was integrated using an implicit third-order upwind solver. The methodology developed in this work is expected to lead to the design optimization of turbomachinery blades.

Gorla, Rama S. R.

Measuring Massive Black Hole Binaries with LISA

The coalescence of two massive black holes produces gravitational waves (GWs) which can be detected by the space-based detector LISA. By measuring these waves, LISA can determine the various parameters which characterize the source. Measurements of the black hole masses and spins will provide information about the growth of black holes and their host galaxies over time. Measurements of a source's sky position and distance may help astronomers identify an electromagnetic counterpart to the GW event. The counterpart's redshift, combined with the GW-measured luminosity distance, can then be used to measure the Hubble constant and the dark energy parameter $w$. Because the potential science output is so high, it is useful to know in advance how well LISA can measure source parameters for a wide range of binaries. We calculate expected parameter estimation errors using the well-known Fisher matrix method. Our waveform model includes the physics of spin precession, as well as subleading harmonics. When these higher-order effects are not included, strong degeneracies between some parameters cause them to be poorly determined by a GW measurement. When precession and subleading harmonics are properly included, the degeneracies are broken, reducing parameter errors by one to several orders of magnitude.

Lang, Ryan N.

Soil Moisture Active Passive Mission L4_SM Data Product Assessment (Version 2 Validated Release)

During the post-launch SMAP calibration and validation (Cal/Val) phase there are two objectives for each science data product team: 1) calibrate, verify, and improve the performance of the science algorithm, and 2) validate the accuracy of the science data product as specified in the science requirements and according to the Cal/Val schedule. This report provides an assessment of the SMAP Level 4 Surface and Root Zone Soil Moisture Passive (L4_SM) product specifically for the product's public Version 2 validated release scheduled for 29 April 2016. The assessment of the Version 2 L4_SM data product includes comparisons of SMAP L4_SM soil moisture estimates with in situ soil moisture observations from core validation sites and sparse networks. The assessment further includes a global evaluation of the internal diagnostics from the ensemble-based data assimilation system that is used to generate the L4_SM product. This evaluation focuses on the statistics of the observation-minus-forecast (O-F) residuals and the analysis increments. Together, the core validation site comparisons and the statistics of the assimilation diagnostics are considered primary validation methodologies for the L4_SM product. Comparisons against in situ measurements from regional-scale sparse networks are considered a secondary validation methodology because such in situ measurements are subject to up-scaling errors from the point-scale to the grid cell scale of the data product. Based on the limited set of core validation sites, the wide geographic range of the sparse network sites, and the global assessment of the assimilation diagnostics, the assessment presented here meets the criteria established by the Committee on Earth Observing Satellites for Stage 2 validation and supports the validated release of the data. An analysis of the time average surface and root zone soil moisture shows that the global pattern of arid and humid regions are captured by the L4_SM estimates. Results from the core validation site comparisons indicate that "Version 2" of the L4_SM data product meets the self-imposed L4_SM accuracy requirement, which is formulated in terms of the ubRMSE: the RMSE (Root Mean Square Error) after removal of the long-term mean difference. The overall ubRMSE of the 3-hourly L4_SM surface soil moisture at the 9 km scale is 0.035 cubic meters per cubic meter requirement. The corresponding ubRMSE for L4_SM root zone soil moisture is 0.024 cubic meters per cubic meter requirement. Both of these metrics are comfortably below the 0.04 cubic meters per cubic meter requirement. The L4_SM estimates are an improvement over estimates from a model-only SMAP Nature Run version 4 (NRv4), which demonstrates the beneficial impact of the SMAP brightness temperature data. L4_SM surface soil moisture estimates are consistently more skillful than NRv4 estimates, although not by a statistically significant margin. The lack of statistical significance is not surprising given the limited data record available to date. Root zone soil moisture estimates from L4_SM and NRv4 have similar skill. Results from comparisons of the L4_SM product to in situ measurements from nearly 400 sparse network sites corroborate the core validation site results. The instantaneous soil moisture and soil temperature analysis increments are within a reasonable range and result in spatially smooth soil moisture analyses. The O-F residuals exhibit only small biases on the order of 1-3 degrees Kelvin between the (re-scaled) SMAP brightness temperature observations and the L4_SM model forecast, which indicates that the assimilation system is largely unbiased. The spatially averaged time series standard deviation of the O-F residuals is 5.9 degrees Kelvin, which reduces to 4.0 degrees Kelvin for the observation-minus-analysis (O-A) residuals, reflecting the impact of the SMAP observations on the L4_SM system. Averaged globally, the time series standard deviation of the normalized O-F residuals is close to unity, which would suggest that the magnitude of the modeled errors approximately reflects that of the actual errors. The assessment report also notes several limitations of the "Version 2" L4_SM data product and science algorithm calibration that will be addressed in future releases. Regionally, the time series standard deviation of the normalized O-F residuals deviates considerably from unity, which indicates that the L4_SM assimilation algorithm either over- or under-estimates the actual errors that are present in the system. Planned improvements include revised land model parameters, revised error parameters for the land model and the assimilated SMAP observations, and revised surface meteorological forcing data for the operational period and underlying climatological data. Moreover, a refined analysis of the impact of SMAP observations will be facilitated by the construction of additional variants of the model-only reference data. Nevertheless, the “Version 2” validated release of the L4_SM product is sufficiently mature and of adequate quality for distribution to and use by the larger science and application communities.

SMAP L4_SM

Calibrating the Helium Pressurization System for the Space Shuttle Liquid-Hydrogen Tank

Analysis of the results from the STS-114 tanking tests and subsequent launch called into question existing thermal and mass models of helium pressurization of the liquid hydrogen tank. This hydrogen tank, which makes up the bottom two-thirds of the External Tank, is pressurized prior to launch to avoid cavitation in the Shuttle Main Engine pumps. At about 2 minutes prior to launch, the main vent valve is closed, and pressurized helium flows into the tank ullage space to achieve set point pressure. As the helium gas cools, its pressure drops, calling for additional helium. Subsequent helium flows are provided in short, timed pulses. The number of pulses is taken as a rough leak indicator. An analysis of thermal models by Marshall Space Flight Center showed considerable uncertainty in the pressure-versus-time behavior of the helium ullage space and the ability to predict the number of pulses normally expected. Kennedy Space Center proposed to calibrate the dime-sized orifice, which together with valves, controls the helium flow quantity (Figure 1). Pressure and temperature sensors were installed to provide upstream and downstream measurements necessary to compute flow rate based on the orifice discharge coefficient. An assessment of flow testing with helium indicated an extremely costly use of this critical resource. In order to reduce costs, we proposed removing the orifices from each Mobile Launcher Platform (MLP) and asking Colorado Engineering Experiment Station Inc. (CEESI) to calibrate the flow. CEESI has a high-pressure air flow system with traceable flow meters capable of handling the large flow rates. However, literature research indicated that square-edged orifices of small diameters often exhibit significant hysteresis and nonrepeatability in the vicinity of choked or sonic flow. Fortunately, the MLP orifices behaved relatively well in testing (Figure 2). Using curve fitting of the air-flow data, in conjunction with ASME orifice modeling equations, a method of relating the helium mass flow to measured air flow data was obtained. This analysis showed that the highest uncertainty in flow occurred in the vicinity of the choking pressure ratio, as would be expected. In addition, analysis of typical flow pulses showed that most of the helium flow occurred either well below or well above this uncertain area. The final result is the ability to provide postlaunch estimates of helium mass flows that are within 1.5 percent of the actual value.

Source record