Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model order reduction”

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 469 records · Page 26

Unstructured-grid large-eddy simulation of flow over an airfoil

Historically, large-eddy simulations (LES) have been restricted to simple geometries where spectral or finite difference methods have dominated due to their efficient use of structured grids. Structured grids, however, not only difficulty representing complex domains and adapting to complicated flow features, but also are rather inefficient for simulating flows at high Reynolds numbers. The lack of efficiency stems from the need to resolve the viscous sublayer, which requires very fine resolution in all three directions near the wall. Structured grids make use of a stretching to reduce the normal grid spacing but must carry the fine resolution in the streamwise and spanwise directions throughout the domain. The unnecessarily fine grid for much of the domain leads to disturbingly high grid estimates. Chapman (1979), and later Moin & Jimenez (1993), pointed out that, in order to advance the technology to airfoils at flight Reynolds numbers, structured grids must be abandoned in lieu of what are known as nested or unstructured grids. The finite element method can efficiently solve the Navier-Stokes equations on unstructured grids. Although the CPU cost per time step per element is somewhat higher than structured grid methods, this effect is more than offset by the reduction in the number of elements. The use of unstructured grids, coupled with the advances in dynamic subgrid-scale modeling such as those made by Germano et al. (1991) and Ghosal et al. (1994), make LES of an airfoil tractable. We have chosen the NACA 4412 airfoil at maximum lift as the first simulation since this flow has not been successfully simulated with the Reynolds-averaged Navier-Stokes equations.

Jansen, Kenneth↗

Calibration of the COBE FIRAS instrument

The Far-Infrared Absolute Spectrophotometer (FIRAS) instrument on the Cosmic Background Explorer (COBE) satellite was designed to accurately measure the spectrum of the cosmic microwave background radiation (CMBR) in the frequency range 1-95/cm with an angular resolution of 7 deg. We describe the calibration of this instrument, including the method of obtaining calibration data, reduction of data, the instrument model, fitting the model to the calibration data, and application of the resulting model solution to sky observations. The instrument model fits well for calibration data that resemble sky condition. The method of propagating detector noise through the calibration process to yield a covariance matrix of the calibrated sky data is described. The final uncertainties are variable both in frequency and position, but for a typical calibrated sky 2.6 deg square pixel and 0.7/cm spectral element the random detector noise limit is of order of a few times 10(exp -7) ergs/sq cm/s/sr cm for 2-20/cm, and the difference between the sky and the best-fit cosmic blackbody can be measured with a gain uncertainty of less than 3%.

Fixsen, D. J.↗

Numerical Modeling of Cavitating Venturi: A Flow Control Element of Propulsion System

In a propulsion system, the propellant flow and mixture ratio could be controlled either by variable area flow control valves or by passive flow control elements such as cavitating venturies. Cavitating venturies maintain constant propellant flowrate for fixed inlet conditions (pressure and temperature) and wide range of outlet pressures, thereby maintain constant, engine thrust and mixture ratio. The flowrate through the venturi reaches a constant value and becomes independent of outlet pressure when the pressure at throat becomes equal to vapor pressure. In order to develop a numerical model of propulsion system, it is necessary to model cavitating venturies in propellant feed systems. This paper presents a finite volume model of flow network of a cavitating venturi. The venturi was discretized into a number of control volumes and mass, momentum and energy conservation equations in each control volume are simultaneously solved to calculate one-dimensional pressure, density, and flowrate and temperature distribution. The numerical model predicts cavitations at the throat when outlet pressure was gradually reduced. Once cavitation starts, with further reduction of downstream pressure, no change in flowrate is found. The numerical predictions have been compared with test data and empirical equation based on Bernoulli's equation.

Majumdar, Alok↗

Magnetogasdynamic Power Extraction and Flow Conditioning for a Gas Turbine

An extension of the Russian AJAX concept to a turbojet is being explored. This magnetohydrodynamic (MHD) energy bypass engine cycle incorporating conventional gas turbine technology has MHD flow conditioning at the inlet to electromagnetically extract part of the inlet air kinetic energy. The electrical power generated can be used for various on-board vehicle requirements including plasma flow control around the vehicle or it may be used for augmenting the expanding flow in the high speed nozzle by MHD forces to generate more thrust. In order to achieve this interaction, the air needs to be ionized by an external means even up to fairly high flight speeds, and the leading candidates may be classified as electrical discharge devices. The present kinetic modeling calculations suggest that the use of electron beams with characteristics close to the commercially available e-beam systems (electron energy approx. 60 keV, beam current approx. 0.2 mA/sq cm) to sustain ionization in intermediate pressure, low-temperature (P = 0.1 atm, T = 300 K) supersonic air flows allows considerable reduction of the flow kinetic energy (up to 10 to 20 percent in M = 3 flows). The calculations also suggest that this can be achieved at a reasonable electron beam efficiency (eta approx. 5), even if the e-beam window losses are taken into account. At these conditions, the exit NO and O atom concentrations due to e-beam initiated chemical reactions do not exceed 30 ppm. Increasing the beam current up to approx. 2 mA/sq cm, which corresponds to a maximum electrical conductivity of sigma(sub max) approx. 0.8 mho/m at the loading parameter of K = 0.5, would result in a much greater reduction of the flow kinetic energy (up to 30 to 40 percent). The MHD channel efficiency at these conditions would be greatly reduced (to eta approx. 1) due to increased electron recombination losses in the channel. At these conditions, partial energy conversion from kinetic energy to heat would result in a significant total pressure loss (P(sub 0)/P(sub 0i) approx. 0.3). The total pressure loss can be reduced operating at the loading parameter closer to unity, at the expense of the reduced electrical power output. Raising the beam current would also result in the increase of the exit O atom concentrations (up to 600 ppm) and NO (up to 150 ppm).

Adamovich, Igor V.↗

A comparison of the results of dynamic wind-tunnel tests with theoretical predictions for an aeromechanical gust-alleviation system for light airplanes

Dynamic wind tunnel tests have been conducted on a 1/6-scale model of a general aviation airplane equipped with an all-mechanical gust alleviation system which uses auxiliary aerodynamic surfaces to drive the flaps. The longitudinal short period motions were studied under simulated gust conditions in order to verify the mathematical model used in a previous study to predict the performance of the full scale system and determine the amount of normal acceleration alleviation which could be attained. The model responses were measured for different configurations with the system active and without the system active for comparison. The tests confirmed the general relationships between the experimental variables and the model responses predicted by the mathematical model, but there were significant differences in the magnitudes of the responses. The experimental results for the model were used to estimate a reduction of 30 percent in the rms normal acceleration response of a similar full scale airplane in atmospheric turbulence.

Stewart, E. C.↗

Simultaneous multifrequency observations of the Seyfert 1 galaxy NGC 4051 - Constant optical-infrared emission observed during large-amplitude X-ray variability

Simultaneous infrared, optical, ultraviolet, and X-ray observations of the low-luminosity Seyfert 1 galaxy, NGC 4051 are reported. A new method of reduction was developed to correct the optical flux for atmospheric variations due to seeing and extinction. The X-ray flux varied by factors of up to 2 on time-scales of tens of minutes, while the optical flux remained steady to within 1 percent. The results rule out all models, including Compton scattering and synchrotron models, in which a single electron population is responsible for the formation of both the infrared-to-optical and the X-ray spectra. The optical emission region must be an order of magnitude larger than, or completely separate from, the X-ray source.

Done, C.↗

Optimal Sensor Selection for Health Monitoring Systems

Sensor data are the basis for performance and health assessment of most complex systems. Careful selection and implementation of sensors is critical to enable high fidelity system health assessment. A model-based procedure that systematically selects an optimal sensor suite for overall health assessment of a designated host system is described. This procedure, termed the Systematic Sensor Selection Strategy (S4), was developed at NASA John H. Glenn Research Center in order to enhance design phase planning and preparations for in-space propulsion health management systems (HMS). Information and capabilities required to utilize the S4 approach in support of design phase development of robust health diagnostics are outlined. A merit metric that quantifies diagnostic performance and overall risk reduction potential of individual sensor suites is introduced. The conceptual foundation for this merit metric is presented and the algorithmic organization of the S4 optimization process is described. Representative results from S4 analyses of a boost stage rocket engine previously under development as part of NASA's Next Generation Launch Technology (NGLT) program are presented.

Santi, L. Michael↗

Jig Twist Optimization of Mach 0.8 Transonic Truss-Braced Wing Aircraft

This paper presents a jig twist optimization study of Mach 0.8 Transonic Truss-Braced Wing (TTBW) aircraft using FUN3D. The solver has been validated against wind tunnel test data. An aeroelastic model deformation correction method is developed for the wind tunnel model to improve the validation. The flight condition corresponding to Mach 0.8 and design lift coefficient 0.695 is selected for the TTBW aircraft jig twist optimization to reduce the drag coefficient. The static aeroelastic simulation module in FUN3D is used for the aeroelastic simulation. Three different orders of the Chebyshev polynomials are used as shape function to represent the variations of the jig twist distribution along the wing span. The preliminary optimization results show that the Mach 0.8 TTBW aircraft with the optimized jig twist achieves a drag reduction of about 5.7 counts at the design condition when 5 𝑡ℎ order Chebyshev polynomials are used as shape function.

TTBW↗

Modeling of Carbon Nanotube Schottky Barrier Reduction for Holes in Air

A model is proposed for the recently observed lower Schottky barrier for holes in air than in vacuum at a metallic electrode - semiconducting carbon nanotube (CNT) junction. In air, the CNT is oxidized. The oxygen molecules on the CNT surface will attract negative charges and will be negatively charged. The negative charges on the solid surface generally enhance the surface dipole and provide stronger electron confinement within the solid. Thus the CNT work function will increase in air. Then the Schottky barrier for holes will have to increase in air according to the standard band-alignment theory, but this is against the experiment. In order to overcome this difficulty, we propose a new Schottky barrier model, assuming that there is a transition region between the electrode and the CNT, where an appreciable potential can drop. We assume that electrodes are placed at the side of the CNT and the contact is formed via the van der Waals interaction. Since the electrode surface is not atomically flat in practical situations, the contact will not be intimate and there will be a transition region before the usual metallic bulk structure starts. In this view, the Schottky barrier is determined not only by the work function difference of the electrodes and the CNT, but also by the potential drop in the transition region. The oxidation increases the potential drop and leads to a lower Schottky barrier for holes. This relation will be demonstrated in a graphical manner by expressing the potential drop as a function of CNT band bending. The potential drop is comparable to the CNT band gap when, e.g., the transition region is about a few Angstroms thick, the oxygen molecule coverage is about 10 %, and an oxygen molecule is charged by about 0.01-0.1 of the unit charge (1-10 % of molecules are charged). The proposed mechanism prevails for both p - and n-CNTs. The model consistently explains the CNT oxidation experiments reported so far.

Yamada, Toshishige↗

Plate measurement techniques and reduction methods used by the West German satellite observers, and resulting consequences for the observation

The West German tracking stations are equipped with ballistic cameras. Plate measurement and plate reduction must therefore follow photogrammetric methods. Approximately 100 star positions and 200 satellite positions are measured on each plate. The mathematical model for spatial rotation of the bundle of rays is extended by including terms for distortion and internal orientation of the camera as well as by providing terms for refraction which are computed for the measured coordinates of the star positions on the plate. From the measuring accuracy of the plate coordinates it follows that the timing accuracy for the exposures has to be about one millisecond, in order to obtain a homogeneous system.

Deker, H.↗

Image Registration: A Necessary Evil

Registration of test and reference images is a key component of nearly all PSP data reduction techniques. This is done to ensure that a test image pixel viewing a particular point on the model is ratioed by the reference image pixel which views the same point. Typically registration is needed to account for model motion due to differing airloads when the wind-off and wind-on images are taken. Registration is also necessary when two cameras are used for simultaneous acquisition of data from a dual-frequency paint. This presentation will discuss the advantages and disadvantages of several different image registration techniques. In order to do so, it is necessary to propose both an accuracy requirement for image registration and a means for measuring the accuracy of a particular technique. High contrast regions in the unregistered images are most sensitive to registration errors, and it is proposed that these regions be used to establish the error limits for registration. Once this is done, the actual registration error can be determined by locating corresponding points on the test and reference images, and determining how well a particular registration technique matches them. An example of this procedure is shown for three transforms used to register images of a semispan model. Thirty control points were located on the model. A subset of the points were used to determine the coefficients of each registration transform, and the error with which each transform aligned the remaining points was determined. The results indicate the general superiority of a third-order polynomial over other candidate transforms, as well as showing how registration accuracy varies with number of control points. Finally, it is proposed that image registration may eventually be done away with completely. As more accurate image resection techniques and more detailed model surface grids become available, it will be possible to map raw image data onto the model surface accurately. Intensity ratio data can then be obtained by a "model surface ratio," rather than an image ratio. The problems and advantages of this technique will be discussed.

Bell, James↗

Reaction engineering for materials processing in space: Reduction of ilmenite by hydrogen and carbon monoxide

Oxygen is a consumable material which needs to be produced continuously in most space missions. Its use for propulsion as well as life support makes oxygen one of the largest volume chemicals to be produced in space. Production of oxygen from lunar materials is of particular interest and is very attractive possibility. The kinetics and mechanism of reduction of ilmenite by carbon monoxide and hydrogen at 800 to 1100 C were investigated. The temporal profiles of conversion for carbon monoxide have a sigmoidal shape and indicate the presence of three different stages (induction, acceleration, and deceleration) during the reduction reaction. The apparent activation energy decreases from 18 kcal/mole at 10 percent conversion to 10 kcal/mole at 50 percent conversion. The reaction is first order with respect to carbon monoxide under the experimental conditions studied. Both SEM and EDX analysis show that the diffusion of Fe product away from the reaction front and through the TiO2 phase, followed by the nucleation and growth of a separate Fe phase are important steps affecting the process kinetics. The results from hydrogen reduction show that the mechanism of ilmenite reduction by hydrogen is similar to that by carbon monoxide. However, the titanium dioxide can be further reduced by hydrogen at 800 to 1000 C. The detailed comparison and theoretical modeling of both reduction processes is presented.

Zhao, Y.↗

The Operation and Evolution of the Swift X-ray Telescope

The Swift X-ray Telescope (XRT) is a CCD based X-ray telescope designed for localization, spectroscopy and long term light curve monitoring of Gamma-Ray Bursts and their X-ray afterglows. Since the launch of Swift in November 2004, the XRT has undergone significant evolution in the way it is operated. Shortly after launch there was a failure of the thermo-electric cooler on the XRT CCD, which led to the XRT team being required to devise a method of keeping the XRT CCD temperature below 50C utilizing only passive cooling by minimizing the exposure of the XRT radiator to the Earth. We present in this paper an update on how the modeling of this passive cooling method has improved in first -1000 days since the method was devised, and the success rate of this method in day-to-day planning. We also discuss the changes to the operational modes and onboard software of the XRT. These changes include improved rapid data product generation in order to improve speed of rapid Gamma-Ray Burst response and localization to the community; changes to the way XRT observation modes are chosen in order to better fine tune data aquisition to a particular science goal; reduction of "mode switching" caused by the contamination of the CCD by Earth light or high temperature effects.

Kennea, Jamie↗

A low-order global gravity field of Venus and dynamical implications

An estimate of sixth-degree and sixth-order harmonic coefficients of the global gravity field of Venus is obtained by processing the long periodic variations of the mean orbital elements of the Pioneer Venus orbiter. Approximately 220 days of data are included in this reduction, which provides almost complete longitudinal coverage. Oblateness is estimated to be -5.97 + or - 3.2 x 10 to the -6. It is noted that the amplitudes of other coefficients are similar to the predicted coefficients using Kaula's rule under equal stress assumption. Atmospheric density values as a function of altitude are obtained to help model drag perturbation. A radial acceleration map at 100 km above the Venus surface is generated, and correlation between gravity anomalies and major topographic features is observed. Spectral analysis of the harmonic model suggests that the interior density anomalies are like those on earth. The orientation angles of the principal axes of the moments of inertia are computed, and deviation of the maximum moment of inertia axis from the spin axis is observed to be small (less than 5 deg). It is found that the minimum moment of inertia axis passes through the Aphrodite and Beta regions.

Ananda, M. P.↗

A generic multi-flex-body dynamics, controls simulation tool for space station

An order (n) multiflex body Space Station simulation tool is introduced. The flex multibody modeling is generic enough to model all phases of Space Station from build up through to Assembly Complete configuration and beyond. Multibody subsystems such as the Mobile Servicing System (MSS) undergoing a prescribed translation and rotation are also allowed. The software includes aerodynamic, gravity gradient, and magnetic field models. User defined controllers can be discrete or continuous. Extensive preprocessing of 'body by body' NASTRAN flex data is built in. A significant aspect, too, is the integrated controls design capability which includes model reduction and analytic linearization.

London, Ken W.↗

Eigenvalue error analysis of viscously damped structures using a Ritz reduction method

The efficient solution of the eigenvalue problem that results from inserting passive dampers with variable stiffness and damping coefficients into a structure is addressed. Eigenanalysis of reduced models obtained by retaining a number of normal modes augmented with Ritz vectors corresponding to the static solutions resulting from the load patterns introduced by the dampers has been empirically shown to yield excellent approximations to the full eigenvalue problem. An analysis of this technique in the case of a single damper is presented. A priori and a posteriori error estimates are generated and tested on numerical examples. Comparison theorems with modally truncated models and a Markov parameter matching reduced-order model are derived. These theorems corroborate the heuristic that residual flexibility methods improve low-frequency approximation of the system. The analysis leads to other techniques for eigenvalue approximation. Approximate closed-form solutions are derived that include a refinement to eigenvalue derivative methods for approximation. An efficient Newton scheme is also developed. A numerical example is presented demonstrating the effectiveness of each of these methods.

Chu, Cheng-Chih↗

Understanding Slat Noise Sources

Model-scale aeroacoustic tests of large civil transports point to the leading-edge slat as a dominant high-lift noise source in the low- to mid-frequencies during aircraft approach and landing. Using generic multi-element high-lift models, complementary experimental and numerical tests were carefully planned and executed at NASA in order to isolate slat noise sources and the underlying noise generation mechanisms. In this paper, a brief overview of the supporting computational effort undertaken at NASA Langley Research Center, is provided. Both tonal and broadband aspects of slat noise are discussed. Recent gains in predicting a slat s far-field acoustic noise, current shortcomings of numerical simulations, and other remaining open issues, are presented. Finally, an example of the ever-expanding role of computational simulations in noise reduction studies also is given.

Khorrami, Medhi R.↗

Observations of the Hot Horizontal Branch Stars in the Metal-Rich Bulge Globular Cluster NGC 6388

The metal-rich bulge globular cluster NGC 6388 shows a distinct blue horizontal-branch tail in its colour-magnitude diagram (Rich et al. 1997) and is thus a strong case of the well-known 2nd Parameter Problem. In addition, its horizontal branch (HB) shows an upward tilt toward bluer colours, which cannot be explained by canonical evolutionary models. Several non-canonical scenarios have been proposed to explain these puzzling observations. In order to test the predictions of these scenarios, we have obtained medium resolution spectra to determine the atmospheric parameters of a sample of the blue HB stars in NGC 6388.Using the medium resolution spectra, we determine effective temperatures, surface gravities and helium abundances by fitting the observed Balmer and helium lines with appropriate theoretical stellar spectra. As we know the distance to the cluster, we can verify our results by determining masses for the stars. During the data reduction we took special care to correctly subtract the background, which is dominated by the overlapping spectra of cool stars. The cool blue tail stars in our sample with T(sub eff) approximately 10000 K have lower than canonical surface gravities, suggesting that these stars are, on average, approximately equal to 0.4 mag brighter than canonical HB stars in agreement with the observed upward slope of the HB in NGC 6388. Moreover, the mean mass of these stars agrees well with theoretical predictions. In contrast, the hot blue tail stars in our sample with T(sub eff) greater than or equal to 12000 K show significantly lower surface gravities than predicted by any scenario, which can reproduce the photometric observations. Their masses are also too low by about a factor of 2 compared to theoretical predictions. The physical parameters of the blue HB stars at about 10,000 K support the helium pollution scenario. The low gravities and masses of the hot blue tail stars, however, are probably caused by problems with the data reduction, most likely due to remaining background light in the spectra, which would affect the fainter hot blue tail stars much more strongly than the brighter cool blue tail stars. Our study of the hot blue tail stars in NGC 6388 illustrates the obstacles which are encountered when attempting to determine the atmospheric parameters of hot HB stars in very crowded fields using ground-based observations. We discuss these obstacles and offer possible solutions for future projects.

Moehler, S.↗