Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative 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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Asymptotic Energies and QED Shifts for Rydberg States of Helium

This paper reviews progress that has been made in obtaining essentially exact solutions to the nonrelativistic three-body problem for helium by a combination of variational and asymptotic expansion methods. The calculation of relativistic and quantum electrodynamic corrections by perturbation theory is discussed, and in particular, methods for the accurate calculation of the Bethe logarithm part of the electron self energy are presented. As an example, the results are applied to the calculation of isotope shifts for the short-lived 'halo' nucleus He-6 relative to He-4 in order to determine the nuclear charge radius of He-6 from high precision spectroscopic measurements carried out at the Argonne National Laboratory. The results demonstrate that the high precision that is now available from atomic theory is creating new opportunities to create novel measurement tools, and helium, along with hydrogen, can be regarded as a fundamental atomic system whose spectrum is well understood for all practical purposes.

Drake, G.W.F.↗
NASA NTRS1987

Estimation of geopotential from satellite-to-satellite range rate data: Numerical results

A technique for high-resolution geopotential field estimation by recovering the harmonic coefficients from satellite-to-satellite range rate data is presented and tested against both a controlled analytical simulation of a one-day satellite mission (maximum degree and order 8) and then against a Cowell method simulation of a 32-day mission (maximum degree and order 180). Innovations include: (1) a new frequency-domain observation equation based on kinetic energy perturbations which avoids much of the complication of the usual Keplerian element perturbation approaches; (2) a new method for computing the normalized inclination functions which unlike previous methods is both efficient and numerically stable even for large harmonic degrees and orders; (3) the application of a mass storage FFT to the entire mission range rate history; (4) the exploitation of newly discovered symmetries in the block diagonal observation matrix which reduce each block to the product of (a) a real diagonal matrix factor, (b) a real trapezoidal factor with half the number of rows as before, and (c) a complex diagonal factor; (5) a block-by-block least-squares solution of the observation equation by means of a custom-designed Givens orthogonal rotation method which is both numerically stable and tailored to the trapezoidal matrix structure for fast execution.

Thobe, Glenn E.↗
NASA NTRS1992

Numerical computation of steady-state acoustic disturbances in flow

Two time domain methods for computing two dimensional steady-state acoustic disturbances propagating through internal subsonic viscous flow fields in the presence of variable area are investigated. The first method solves the Navier-Stokes equations for the combined steady and acoustic field together and subtracts the steady flow to obtain the acoustic field. The second method solves a system of perturbation equations to obtain the acoustic disturbances, making use of a separate steady flow computation as input to the system. In each case the periodic steady-state acoustic fluctuations are obtained numerically on a supercomputer using a second order unsplit explicit MacCormack predictor-corrector method. Results show that the first method is not very effective for computing acoustic disturbances of even moderate amplitude. It appears that more accurate steady flow algorithms are required for this method to succeed. On the other hand, linear and nonlinear acoustic disturbances extracted from the perturbation approach are shown to exhibit expected behavior for the problems considered. It is also found that inflow boundary conditions for an equivalent uniform duct can be successfully applied to a nonuniform duct to obtain steady-state acoustic disturbances.

Watson, W. R.↗
NASA NTRS1959

Source Distribution Method for Unsteady One-Dimensional Flows With Small Mass, Momentum, and Heat Addition and Small Area Variation

A source distribution method is presented for obtaining flow perturbations due to small unsteady area variations, mass, momentum, and heat additions in a basic uniform (or piecewise uniform) one-dimensional flow. First, the perturbations due to an elemental area variation, mass, momentum, and heat addition are found. The general solution is then represented by a spatial and temporal distribution of these elemental (source) solutions. Emphasis is placed on discussing the physical nature of the flow phenomena. The method is illustrated by several examples. These include the determination of perturbations in basic flows consisting of (1) a shock propagating through a nonuniform tube, (2) a constant-velocity piston driving a shock, (3) ideal shock-tube flows, and (4) deflagrations initiated at a closed end. The method is particularly applicable for finding the perturbations due to relatively thin wall boundary layers.

Mirels, Harold↗
NASA NTRS1994

A Monte Carlo study of photochemical model uncertainties affecting the calculation of subsonic perturbations in the upper troposphere

Monte Carlo methods are frequently applied to the evaluation of uncertainties in models with multiple inputs that themselves have associated imprecisions. In the case of photochemical models used to evaluate changes in O3 or OH, inputs analyzed include rate coefficients measured in the laboratory, and chemical and physical constituents measured in the atmosphere. The Monte Carlo method was used with the 1-dimensional GSFC tropospheric photochemical model to examine uncertainty propagation to calculation of NO(x) and the major oxidants (O3, OH) in the upper troposphere. In all cases chemical kinetics inputs are varied and the NO(x) perturbation of a subsonic fleet is simulated. A series of model runs is used to explore sensitivities of model-computed parameters to other parameters: heterogeneous processes, the uncertainty in upper tropospheric H2O vapor measurements, aircraft emissions at different latitudes.

Thompson, Anne M.↗
NASA NTRS1990

Thermal instabilities in protogalactic clouds

The means by which a protogalaxy can fragment to form the first generation of stars and globular clusters remains an important problem in astrophysics. Gravitational instabilities grow on timescales too long to drive fragmentation before the background density grows by many orders of magnitude (see Murray and Lin 1989a, and references therein). Thermal instability provides a much more likely mechanism. After its initial collapse, a protogalactic cloud is expected to be shock heated to its virial temperature approx. 10(exp 6) K. Cooling by H and He+ below 10(exp 6) K has a negative slope, so that the cloud is subject to strong thermal instabilities. Density enhancements may then grow rapidly, fragmenting the protogalaxy as it cools to lower temperatures. The role of dynamical effects upon the growth of perturbations is considered here. The method used is similar to that used in Murray and Lin (1989a; see also the Erratum to appear September 15), which examined the growth of thermal instabilities with a one-dimensional Lagrangian hydrodynamics code, written for spherical symmetry. Perturbed regions therefore take the form of shells. The dynamical variables are integrated explicitly, while the temperature, ionization fraction, and molecular fraction are integrated implicitly, and account is taken for non-equilibrium values of these quantities.

Murray, S. D.↗
NASA NTRS1989

Feasibility study of a single, elliptical heliocentric Earth-Mars trajectory

The initial intent of this design project was to evaluate the existence and feasibility of a single elliptical heliocentric Earth/Mars trajectory. This trajectory was constrained to encounter Mars twice in its orbit, within a time interval of 15 to 180 Earth days between encounters. The single ellipse restriction was soon found to be prohibitive for reasons shown later. Therefore, the approach taken in the design of the round-trip mission to Mars was to construct single-leg trajectories which connected two planets on two prescribed dates. Three methods of trajectory design were developed. Method 1 is an eclectic approach and employs Gaussian Orbit Determination (Method 1A) and Lambert-Euler Preliminary Orbit Determination (Method 1B) in conjunction with each other. Method 2 is an additional version of Lambert's Solution to orbit determination, and both a coplanar and a noncoplanar solution were developed within Method 2. In each of these methods, the fundamental variables are two position vectors and the time between the position vectors. In all methods, the motion was considered Keplerian motion and the reference frame origin was located at the sun. Perturbative effects were not considered in Method 1. The feasibility study of round-trip Earth/Mars trajectories maintains generality by considering only heliocentric trajectory parameters and planetary approach conditions. The coordinates and velocity components of the planets, for the standard epoch J2000, were computed from an approximate set of osculating elements by the procedure outlined in an ephemeris of coordinates.

Blake, M.↗
NASA NTRS1976

A wave-envelope technique for wave propagation in non-uniform ducts

The method of variation of parameters is used to analyze wave propagation in variable area, plane ducts with no mean flow. The method has an explicit representation of the fast axial variation of the acoustic modes, and numerical integration is required only for the slower axial variations of the mode amplitudes and phases. Results are presented which demonstrate the numerical advantages of the method. Comparison of the results with those of a small perturbation theory are given. The relationship between this method and the method of weighted residuals is discussed.

Kaiser, J. E.↗
NASA NTRS2011

Predicting Spacecraft Trajectories by the WeavEncke Method

A combination of methods is proposed of predicting spacecraft trajectories that possibly include multiple maneuvers and/or perturbing accelerations, with greater speed, accuracy, and repeatability than were heretofore achievable. The combination is denoted the WeavEncke method because it is based on unpublished studies by Jonathan Weaver of the orbit-prediction formulation of the noted astronomer Johann Franz Encke. Weaver evaluated a number of alternatives that arise within that formulation, arriving at an orbit-predicting algorithm optimized for complex trajectory operations. In the WeavEncke method, Encke's method of prediction of perturbed orbits is enhanced by application of modern numerical methods. Among these methods are efficient Kepler s-equation time-of-flight solutions and self-starting numerical integration with time as the independent variable. Self-starting numerical integration satisfies the requirements for accuracy, reproducibility, and efficiency (and, hence, speed). Self-starting numerical integration also supports fully analytic regulation of integration step sizes, thereby further increasing speed while maintaining accuracy.

Weaver, Jonathan K.↗
NASA NTRS1987

Shape design sensitivity analysis and optimization of three dimensional elastic solids using geometric modeling and automatic regridding

An automatic regridding method and a three dimensional shape design parameterization technique were constructed and integrated into a unified theory of shape design sensitivity analysis. An algorithm was developed for general shape design sensitivity analysis of three dimensional eleastic solids. Numerical implementation of this shape design sensitivity analysis method was carried out using the finite element code ANSYS. The unified theory of shape design sensitivity analysis uses the material derivative of continuum mechanics with a design velocity field that represents shape change effects over the structural design. Automatic regridding methods were developed by generating a domain velocity field with boundary displacement method. Shape design parameterization for three dimensional surface design problems was illustrated using a Bezier surface with boundary perturbations that depend linearly on the perturbation of design parameters. A linearization method of optimization, LINRM, was used to obtain optimum shapes. Three examples from different engineering disciplines were investigated to demonstrate the accuracy and versatility of this shape design sensitivity analysis method.

Yao, Tse-Min↗
NASA NTRS1976

Rotor Dynamic State and Parameter Identification from Hovering Transients

State and parameter identifications based on a form of the maximum likelihood method are applied to the problem of extracting linear perturbation models, including rotor dynamic inflow effects, from transient blade flapping measurements. The estimation method is first studied in computer simulations and then applied to cyclic pitch stirring transients generated with a four-bladed rotor model operating in hovering trim conditions. The analytical perturbation models extracted from the transient test results are compared with transient and frequency response tests not used in the state and parameter identification. The identified analytical perturbation model is also compared with a simple theory. The method that is applicable both to small scale and full scale dynamic rotor testing is being extended to perturbations from forward flight trim conditions.

Hohenemser, K. H.↗
NASA NTRS1997

Computational Investigation of Experimental Interaction Impedance Obtained by Perturbation for Helical Traveling-Wave Tube Structures

Conventional methods used to measure the cold-test interaction impedance of helical slow-wave structures involve perturbing a helical circuit with a cylindrical dielectric rod placed on the central axis of the circuit. It has been shown that the difference in resonant frequency or axial phase shift between the perturbed and unperturbed circuits can be related to the interaction impedance. However, because of the complex configuration of the helical circuit, deriving this relationship involves-several approximations. With the advent of accurate three- dimensional helical circuit models, these standard approximations can be fully investigated. This paper addresses the most prominent approximations made in the analysis for measured interaction impedance by Lagerstrom and investigates their accuracy using the three-dimensional simulation code MAFIA. It is shown that a more accurate value of interaction impedance can be obtained by using three-dimensional computational methods rather than performing costly and time consuming experimental cold-test measurements.

Kory, Carol L.↗
NASA NTRS1998

Computational Investigation of Experimental Interaction Impedance Obtained by Perturbation for Helical Traveling-Wave Tube Structures

Conventional methods used to measure the cold-test interaction impedance of helical slow-wave structures involve perturbing a helical circuit with a cylindrical dielectric rod placed on the central axis of the circuit. It has been shown that the difference in resonant frequency or axial phase shift between the perturbed and unperturbed circuits can be related to the interaction impedance. However, because of the complex configuration of the helical circuit, deriving this relationship involves several approximations. With the advent of accurate three-dimensional (3-D) helical circuit models, these standard approximations can be fully investigated. This paper addresses the most prominent approximations made in the analysis for measured interaction impedance by Lagerstrom and investigates their accuracy using the 3-D simulation code MAFIA. It is shown that a more accurate value of interaction impedance can be obtained by using 3-D computational methods rather than performing costly and time consuming experimental cold-test measurements.

Kory, Carol L.↗
NASA NTRS1976

On the Tesseral-Harmonics Resonance Problem in Artificial-Satellite Theory, Part 2

Equations were derived for the perturbations on an artificial satellite when the motion of the satellite is commensurable with that of the earth. This was done by first selecting the tesseral harmonics that contribute the most to the perturbations and then by applying Hori's method by use of Lie series. Here, are introduced some modifications to the perturbations, which now result in better agreement with numerical integration.

Romanowicz, B. A.↗
NASA NTRS1979

A higher order panel method for linearized supersonic flow

The basic integral equations of linearized supersonic theory for an advanced supersonic panel method are derived. Methods using only linear varying source strength over each panel or only quadratic doublet strength over each panel gave good agreement with analytic solutions over cones and zero thickness cambered wings. For three dimensional bodies and wings of general shape, combined source and doublet panels with interior boundary conditions to eliminate the internal perturbations lead to a stable method providing good agreement experiment. A panel system with all edges contiguous resulted from dividing the basic four point non-planar panel into eight triangular subpanels, and the doublet strength was made continuous at all edges by a quadratic distribution over each subpanel. Superinclined panels were developed and tested on s simple nacelle and on an airplane model having engine inlets, with excellent results.

Ehlers, F. E.↗
NASA NTRS1985

Perturbation-iteration theory for analyzing microwave striplines

A perturbation-iteration technique is presented for determining the propagation constant and characteristic impedance of an unshielded microstrip transmission line. The method converges to the correct solution with a few iterations at each frequency and is equivalent to a full wave analysis. The perturbation-iteration method gives a direct solution for the propagation constant without having to find the roots of a transcendental dispersion equation. The theory is presented in detail along with numerical results for the effective dielectric constant and characteristic impedance for a wide range of substrate dielectric constants, stripline dimensions, and frequencies.

Kretch, B. E.↗
NASA NTRS1978

Perturbations of non-resonant satellite orbits due to a rotating earth

The dominant perturbations of the motion of a satellite near the earth are due to atmospheric drag and the non-symmetrical gravitational field. Atmospheric drag perturbation continually pulls the satellite in and out of the different long period resonant frequencies. The result is that the resonances never become apparent and may be neglected. The tesseral harmonics have no true secular perturbation but the periodicities in the mean motion induce a secular perturbation in the mean anomaly. This secular perturbation may be determined by simply using the average mean motion instead of the osculating mean motion. The Von Ziepel method is used to determine tesseral perturbations. The solution is found first in the singular DS phi elements and then rewritten in the PS phi elements to remove singularities. The notation used in the development is described in the appendix.

Mueller, A.↗