Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximation methods”

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

Aerodynamic analysis of hypersonic waverider aircraft

The purpose of this study is to validate two existing codes used by the Systems Analysis Branch at NASA ARC, and to modify the codes so they can be used to generate and analyze waverider aircraft at on-design and off-design conditions. To generate waverider configurations and perform the on-design analysis, the appropriately named Waverider code is used. The Waverider code is based on the Taylor-Maccoll equations. Validation is accomplished via a comparison with previously published results. The Waverider code is modified to incorporate a fairing to close off the base area of the waverider configuration. This creates a more realistic waverider. The Hypersonic Aircraft Vehicle Optimization Code (HAVOC) is used to perform the off-design analysis of waverider configurations generated by the Waverider code. Various approximate analysis methods are used by HAVOC to predict the aerodynamic characteristics, which are validated via a comparison with experimental results from a hypersonic test model.

Sandlin, Doral R.↗

A generalized Zel'dovich approximation to gravitational instability

The orbits of particles undergoing gravitational instability are parameterized by generalizing the Zel'dovich approximation to a series expansion of arbitrary accuracy in the nonlinear regime. The coefficients of this series are determined from an action principle, or, more generally, from moments of the equation of motion. It is shown that the series is more rapidly convergent than previous nonlinear approximations. The method therefore provides a practical means of determining particle orbits, even for highly nonlinear perturbations. As an alternative, we also outline how the nonlinear dynamics may be computed as a field theory in which the evolution of the density and the velocity is determined in fixed comoving Eulerian coordinates.

Giavalisco, M.↗

Linear, multivariable robust control with a mu perspective

The structured singular value is a linear algebra tool developed to study a particular class of matrix perturbation problems arising in robust feedback control of multivariable systems. These perturbations are called linear fractional, and are a natural way to model many types of uncertainty in linear systems, including state-space parameter uncertainty, multiplicative and additive unmodeled dynamics uncertainty, and coprime factor and gap metric uncertainty. The structured singular value theory provides a natural extension of classical SISO robustness measures and concepts to MIMO systems. The structured singular value analysis, coupled with approximate synthesis methods, make it possible to study the tradeoff between performance and uncertainty that occurs in all feedback systems. In MIMO systems, the complexity of the spatial interactions in the loop gains make it difficult to heuristically quantify the tradeoffs that must occur. This paper examines the role played by the structured singular value (and its computable bounds) in answering these questions, as well as its role in the general robust, multivariable control analysis and design problem.

Packard, Andy↗

Analysis of Ground Effect on the Lifting Airscrew

A study is presented of ground effect as applied to the lifting airscrew of the type used in modern gyroplane and helicopters. The mathematical analysis of the problem has been verified by tests made of three rotor models in the presence of a large circular "ground plane." A method of approximating the decelerating effect during vertical approach to the ground has been devised. The results of the study are presented in the form of convenient charts.

Knight, Montgomery↗

An aerodynamic analysis of the autogiro rotor with a comparison between calculated and experimental results

This report presents an extension of the autogiro theory of Glauert and Lock in which the influence of a pitch varying with the blade radius is evaluated and methods of approximating the effect of blade tip losses and the influence of reversed velocities on the retreating blades are developed. A comparison of calculated and experimental results showed that most of the rotor characteristics could be calculated with reasonable accuracy, and that the type of induced flow assumed has a secondary effect upon the net rotor forces, although the flapping motion is influenced appreciably. An approximate evaluation of the effect of parasite drag on the rotor blades established the importance of including this factor in the analysis.

Wheatley, John B↗

The longitudinal stability of elastic swept wings at supersonic speed

The longitudinal stability characteristics of elastic swept wings of high aspect ratio experiencing bending and torsional deformations are calculated for supersonic speed by the application of linearized lifting-surface theory. A parabolic wing deflection curve is assumed and the analysis is simplified by a number of structural approximations. The method is thereby limited in application to wings of high aspect ratio for which the root effects are small. Expressions for the lift, pitching-moment, and span load distribution characteristics are derived in terms of the elastic properties of the wing; namely, the design stress, the modulus of elasticity, the shearing modulus, and the maximum design load factor. The analysis applies to wings with leading edges swept behind the Mach lines. In all cases, however, the trailing edge is sonic or supersonic. Application of the method of analysis to wings with leading edges swept ahead of the Mach lines is discussed.

Frick, C W↗

On the determination of the take-off characteristics of a seaplane

The present paper presents an attempt to coordinate the available theoretical and experimental data on planing surfaces so as to develop an approximate analytical method for the determination of the water resistance of a seaplane without any preliminary towing tests in the tank.

Perelmuter, A↗

Interlaminar stress analysis of dropped-ply laminated plates and shells by a mixed method

A mixed method of approximation based on Reissner's variational principle is developed for the linear analysis of interlaminar stresses in laminated composites, with special interest in laminates that contain terminated internal plies (dropped-ply laminates). Two models are derived, one for problems of generalized plane deformation and the other for the axisymmetric response of shells of revolution. A layerwise approach is taken in which the stress field is assumed with an explicit dependence on the thickness coordinate in each layer. The dependence of the stress field on the thickness coordinate is determined such that the three-dimensional equilibrium equations are satisfied by the approximation. The solution domain is reduced to one dimension by integration through the thickness. Continuity of tractions and displacements between layers is imposed. The governing two-point boundary value problem is composed of a system of both differential and algebraic equations (DAE's) and their associated boundary conditions. Careful evaluation of the system of DAE's was required to arrive at a form that allowed application of a one-step finite difference approximation. A two-stage Gauss implicit Runge-Kutta finite difference scheme was used for the solution because of its relatively high degree of accuracy. Patch tests of the two models revealed problems with solution accuracy for the axisymmetric model of a cylindrical shell loaded by internal pressure. Parametric studies of dropped-ply laminate characteristics and their influence on the interlaminar stresses were performed using the generalized plane deformation model. Eccentricity of the middle surface of the laminate through the ply drop-off was found to have a minimal effect on the interlaminar stresses under longitudinal compression, transverse tension, and in-plane shear. A second study found the stiffness change across the ply termination to have a much greater influence on the interlaminar stresses.

Harrison, Peter N.↗

Kinematics of the ion tail of comet P/Swift-Tuttle

We have obtained long-slit high resolution spectra of the H2O(+) 6199 A complex in the near tail of comet P/Swift-Tuttle. The observations were made using the Hamilton echelle spectrometer fed by the Lick Observatory 0.6 m coude auxiliary telescope. For most of our observations, the spectral slit was aligned along the Sun-tail axis and the cometary nucleus was placed at one end of the slit, giving us spectra having the spatial and spectral resolution needed to measure the radial velocity and velocity dispersion continuously down the cometary tail out to a distance of 4X10(exp 5) km. The radial velocities confirm the earlier more restricted observations by Rauer & Jockers (1993) and by Wyckoff & Lindholm (1994) showing that the tail motions are indeed bulk flows in the antisolar direction. Out to 3X10(exp 5) km in the tail typical bulk flows are at a speed of approximately 30 km/s. The velocity dispersion, (sigma(sub r)), of the H2O(+) lines follows a pattern that is quite systematic; sigma(sub r) is smallest near the cometary nucleus, and steadily increases down the tail. The highest velocity dispersions are found ahead of the nucleus and off the tail axis. These velocity dispersions are equivalent to ion temperatures ranging from 10(exp 5) to 10(exp 6) K. We note a clear anticorrelation between the H2O(+) line intensities (related to the ion density) and the bulk flow and dispersion velocities; direct mass loading of the solar wind by the observed water ions may be responsible. We discuss several approximate equipartition methods used to infer local magnetic fields induced by the interaction of the cometary ions with the solar wind particle/field stream. Typical fields derived are near 50 nT. The measured tailward accelerations are consistent with this order of magnitude B field.

Spinrad, Hyron↗

Experimental Investigation of Impact in Landing on Water

The extent of agreement of the theoretical impact computations with the actual phenomenon has not as yet been fully clarified. There is on the one hand a certain imperfection in the theory (simplifying assumptions made) and on the other an insufficiency in the experimental data available. The object of our present paper is to show how far test results agree with the available approximate computation methods, to investigate in greater detail the physical nature of impact on water, and to perfect the experimental method of studying the phenomenon.

Kreps, R. L.↗

Charge Resolution of the Silicon Matrix of the ATIC Experiment

ATIC (Advanced Thin Ionization Calorimeter) is a balloon borne experiment designed to measure the cosmic ray composition for elements from hydrogen to iron and their energy spectra from approx.50 GeV to near 100 TeV. It consists of a Si-matrix detector to determine the charge of a CRT particle, a scintillator hodoscope for tracking, carbon interaction targets and a fully active BGO calorimeter. ATIC had its first flight from McMurdo, Antarctica from 28/12/2000 to 13/01/2001. The ATIC flight collected approximately 25 million events. The silicon matrix of the ATIC spectrometer is designed to resolve individual elements from proton to iron. To provide this resolution careful calibration of each pixel of the silicon matrix is required. Firstly, for each electronic channel of the matrix the pedestal value was subtracted taking into account its drift during the flight. The muon calibration made before the flight was used then to convert electric signals (in ADC channel number) to energy deposits in each pixel. However, the preflight muon calibration was not accurate enough for the purpose, because of lack of statistics in each pixel. To improve charge resolution the correction was done for the position of Helium peak in each pixel during the flight . The other way to set electric signals in electronics channels of the Si-matrix to one scale was correction for electric channel gains accurately measured in laboratory. In these measurements it was found that small different nonlinearities for different channels are present in the region of charge Z > 20. The correction for these non-linearities was not done yet. In linear approximation the method provides practically the same resolution as muon calibration plus He-peak correction. For searching a pixel with the signal of primary particle an indication from the cascade in the calorimeter was used. For this purpose a trajectory was reconstructed using weight centers of energy deposits in BGO layers. The point of intersection of this trajectory with Si-matrix and its RMS was determined. The pixel with maximal signal in 3sigma region was taken as sought. The signal in this pixel was corrected by trajectory zenith angle. The preliminary results on charge resolution of the Si-matrix in the range from protons to iron are presented.

Zatsepin, V. I.↗

Evaluation of Methods to Estimate the Surface Downwelling Longwave Flux during Arctic Winter

Surface longwave radiation fluxes dominate the energy budget of nighttime polar regions, yet little is known about the relative accuracy of existing satellite-based techniques to estimate this parameter. We compare eight methods to estimate the downwelling longwave radiation flux and to validate their performance with measurements from two field programs in thc Arctic: the Coordinated Eastern Arctic Experiment (CEAREX ) conducted in the Barents Sea during the autumn and winter of 1988, and the Lead Experiment performed in the Beaufort Sea in the spring of 1992. Five of the eight methods were developed for satellite-derived quantities, and three are simple parameterizations based on surface observations. All of the algorithms require information about cloud fraction, which is provided from the NASA-NOAA Television and Infrared Observation Satellite (TIROS) Operational Vertical Sounder (TOVS) polar pathfinder dataset (Path-P): some techniques ingest temperature and moisture profiles (also from Path-P): one-half of the methods assume that clouds are opaque and have a constant geometric thickness of 50 hPa, and three include no thickness information whatsoever. With a somewhat limited validation dataset, the following primary conclusions result: (1) all methods exhibit approximately the same correlations with measurements and rms differences, but the biases range from -34 W sq m (16% of the mean) to nearly 0; (2) the error analysis described here indicates that the assumption of a 50-hPa cloud thickness is too thin by a factor of 2 on average in polar nighttime conditions; (3) cloud-overlap techniques. which effectively increase mean cloud thickness, significantly improve the results; (4) simple Arctic-specific parameterizations performed poorly, probably because they were developed with surface-observed cloud fractions; and (5) the single algorithm that includes an estimate of cloud thickness exhibits the smallest differences from observations.

Chiacchio, Marc↗

A Physical Model to Estimate Snowfall over Land using AMSU-B Observations

In this study, we present an improved physical model to retrieve snowfall rate over land using brightness temperature observations from the National Oceanic and Atmospheric Administration's (NOAA) Advanced Microwave Sounder Unit-B (AMSU-B) at 89 GHz, 150 GHz, 183.3 +/- 1 GHz, 183.3 +/- 3 GHz, and 183.3 +/- 7 GHz. The retrieval model is applied to the New England blizzard of March 5, 2001 which deposited about 75 cm of snow over much of Vermont, New Hampshire, and northern New York. In this improved physical model, prior retrieval assumptions about snowflake shape, particle size distributions, environmental conditions, and optimization methodology have been updated. Here, single scattering parameters for snow particles are calculated with the Discrete-Dipole Approximation (DDA) method instead of assuming spherical shapes. Five different snow particle models (hexagonal columns, hexagonal plates, and three different kinds of aggregates) are considered. Snow particle size distributions are assumed to vary with air temperature and to follow aircraft measurements described by previous studies. Brightness temperatures at AMSU-B frequencies for the New England blizzard are calculated using these DDA calculated single scattering parameters and particle size distributions. The vertical profiles of pressure, temperature, relative humidity and hydrometeors are provided by MM5 model simulations. These profiles are treated as the a priori data base in the Bayesian retrieval algorithm. In algorithm applications to the blizzard data, calculated brightness temperatures associated with selected database profiles agree with AMSU-B observations to within about +/- 5 K at all five frequencies. Retrieved snowfall rates compare favorably with the near-concurrent National Weather Service (NWS) radar reflectivity measurements. The relationships between the NWS radar measured reflectivities Z(sub e) and retrieved snowfall rate R for a given snow particle model are derived by a histogram matching technique. All of these Z(sub e)-R relationships fall in the range of previously established Z(sub e)-R relationships for snowfall. This suggests that the current physical model developed in this study can reliably estimate the snowfall rate over land using the AMSU-B measured brightness temperatures.

Kim, Min-Jeong↗

Estimating Shock Spectra: Extensions beyond GEVS

Shock response spectra (SRS) are the standard description of some vibration environments on spacecraft for equipment qualification. For shock events produced by pyrotechnic devices, SRS can have significant frequency content as high as 10 kHz. It is difficult to construct and analyze finite element models that can resolve dynamic behavior at such high frequencies. GEVS provides simple, empirically based methods for approximating the SRS for a wide variety of shock events. It begins with a base SRS according to the type of pyrotechnic device, and then provides attenuation relations to adjust this SRS according to distance from the shock source, the type of structural frame and the properties of any structural joints between the source and equipment. In our paper we extend GEVS to include more detailed information about the spacecraft structure. To retain the general framework of GEVS, we begin with a base SRS and adjust this SRS using attenuation relations. We use modal and traveling wave concepts to derive the attenuation relations for simple canonical structures. Then we show how these concepts can be used to analyze more complex structures using finite element mode shapes to explicitly calculate the attenuation factors. Since the low- to mid-frequency finite element modal information is extrapolated to obtain the low- to high-frequency attenuation relations, the resulting attenuated SRS is formulated as an upper bound rather than as mean predicted values. We illustrate the extended GEVS approach by analyzing the impact response of composite tubes and the shock response of the STEREO spacecraft.

Igusa, Takeru↗

Sensing a Changing Chemical Mixture Using an Electronic Nose

A method of using an electronic nose to detect an airborne mixture of known chemical compounds and measure the temporally varying concentrations of the individual compounds is undergoing development. In a typical intended application, the method would be used to monitor the air in an inhabited space (e.g., the interior of a building) for the release of solvents, toxic fumes, and other compounds that are regarded as contaminants. At the present state of development, the method affords a capability for identifying and quantitating one or two compounds that are members of a set of some number (typically of the order of a dozen) known compounds. In principle, the method could be extended to enable monitoring of more than two compounds. An electronic nose consists of an array of sensors, typically made from polymer carbon composites, the electrical resistances of which change upon exposure to a variety of chemicals. By design, each sensor is unique in its responses to these chemicals: some or all of the sensitivities of a given sensor to the various vapors differ from the corresponding sensitivities of other sensors. In general, the responses of the sensors are nonlinear functions of the concentrations of the chemicals. Hence, mathematically, the monitoring problem is to solve the set of time-dependent nonlinear equations for the sensor responses to obtain the time dependent concentrations of individual compounds. In the present developmental method, successive approximations of the solution are generated by a learning algorithm based on independent-component analysis (ICA) an established information theoretic approach for transforming a vector of observed interdependent signals into a set of signals that are as nearly statistically independent as possible.

Duong, Tuan↗

A Simple Method for High-Lift Propeller Conceptual Design

In this paper, we present a simple method for designing propellers that are placed upstream of the leading edge of a wing in order to augment lift. Because the primary purpose of these "high-lift propellers" is to increase lift rather than produce thrust, these props are best viewed as a form of high-lift device; consequently, they should be designed differently than traditional propellers. We present a theory that describes how these props can be designed to provide a relatively uniform axial velocity increase, which is hypothesized to be advantageous for lift augmentation based on a literature survey. Computational modeling indicates that such propellers can generate the same average induced axial velocity while consuming less power and producing less thrust than conventional propeller designs. For an example problem based on specifications for NASA's Scalable Convergent Electric Propulsion Technology and Operations Research (SCEPTOR) flight demonstrator, a propeller designed with the new method requires approximately 15% less power and produces approximately 11% less thrust than one designed for minimum induced loss. Higher-order modeling and/or wind tunnel testing are needed to verify the predicted performance.

Patterson, Michael↗

Evaluation of Higher-order Quadrature Schemes in Improving Computational Efficiency for Orientation-averaged Single-Scattering Properties of Nonspherical Ice Particles

We evaluate several high-order quadrature schemes for accuracy and efficacy in obtaining orientation-averaged single-scattering properties (SSPs). We use the highly efficient MIDAS to perform electromagnetic scattering calculations to evaluate the gain in efficiency from these schemes. MIDAS is shown to be superior to DDSCAT, a popular discrete dipole approximation (DDA) method. This study is motivated by the fact that quality physical precipitation retrievals rely on using accurate orientation-averaged SSPs derived from realistic hydrometeors as input to radiative transfer Models (RTMs). The DDA has been a popular choice for single-scattering calculations, due to its versatility with respect to target geometry. However, being iterative-solver-based (ISB), the most used DDA codes, e.g. DDSCAT and ADDA, must solve the scattering problem for each orientation of the target separately. As the size parameter and geometric anisotropy of the hydrometeor increase, the number of orientations needed to obtain accurate orientation-averages can increase drastically and so does the computation cost incurred by the ISB-DDA methods. MIDAS is a Direct-Solver-Based (DSB) code, its decomposition of the original large matrix with a high rank into multiple more manageable smaller matrices of lower ranks makes it much more computationally efficient while maintaining excellent accuracy. In addition, direct solvers consider all requested orientations at once, giving MIDAS further advantage over popular ISB-DDA methods. MIDAS, when combined with high-order quadrature for orientation averaging, can be greater than three orders of magnitude more efficient in obtaining RTM-ready SSPs of complex-shaped hydrometeors than existing ISB-DDA methods, with the native quadrature schemes they offer.

Ines Fenni↗

Approximate line shapes for hydrogen

Two independent methods are presented for calculating radiative transport within hydrogen lines. In Method 1, a simple equation is proposed for calculating the line shape. In Method 2, the line shape is assumed to be a dispersion profile and an equation is presented for calculating the half half-width. The results obtained for the line shapes and curves of growth by the two approximate methods are compared with similar results using the detailed line shapes by Vidal et al.

Sutton, K.↗