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 1,279 records · Page 71

A Non-Intrusive Algorithm for Sensitivity Analysis of Chaotic Flow Simulations

We demonstrate a novel algorithm for computing the sensitivity of statistics in chaotic flow simulations to parameter perturbations. The algorithm is non-intrusive but requires exposing an interface. Based on the principle of shadowing in dynamical systems, this algorithm is designed to reduce the effect of the sampling error in computing sensitivity of statistics in chaotic simulations. We compare the effectiveness of this method to that of the conventional finite difference method.

Blonigan, Patrick J.↗

The Third General Scientific Assembly of the International Association of Geomagnetism and Aeronomy - Special sessions of auroral processes

Methods of timing magnetic substorms, the rapid fluctuations of aurorae, electromagnetic and electrostatic instabilities observed on the field lines of aurorae, the auroral microstructure, and the relationship of currents, electric field and particle precipitation to auroral form are discussed. Attention is given to such topics as D-perturbations as an indicator of substorm onset, the role of the magnetotail in substorms, spectral information derived from imaging data on aurorae, terrestrial kilometric radiation, and the importance of the mirror force in self-consistent models of particle fluxes, currents and potentials on auroral field lines.

Russell, C. T.↗

Gravity wave propagation in a diffusively separated atmosphere with height-dependent collision frequencies

Numerical calculations of gravity wave propagation through a two-component thermosphere with vertically-varying inter-species collision frequencies are performed using a direct integration method. Reflection of upgoing gravity wave energy into downgoing gravity waves appears to be small though nonnegligible at typical thermospheric periods and wavelengths. Coupling into diffusion waves is completely insignificant in most cases and may be determined in part by the relative abundance of the minor species. Temperature perturbation differences between species can be greater than 10% above 300 km altitude even though coupling into diffusion waves is small. Diffusive dissipation of gravity wave energy can be significant in the lower thermosphere and may be comparable to dissipation by viscosity and heat conduction below about 225 km altitude. Heating by diffusive dissipation may also be comparable to direct heating by solar radiation in the lower thermosphere. As a result of dissipative filtering of long-period and short-wavelength waves and the sensitivity of He-N2 density perturbation phase differences to diffusion, the AE-C satellite wave observations can only be fit by internal gravity waves with periods no greater than about 25 min and horizontal wavelengths no less than about 200 km.

Del Genio, A. D.↗

The RAPT control method for nonlinear systems

The Rapid Advancement Preview Tracking (RAPT) control method uses both a nonlinear system model and a linearized system model to provide a form of feedback/preview control for nonlinear systems. It uses measurements of the states to calculate the current preview control signal from a local linearization of the state equations, and then iteratively calculates a perturbation correction for that control using the nonlinear equations. The corrected control signal is then applied to the system, driving the end-effector to the next position along a specified trajectory. Computer simulations of tracking control for two-link planar manipulators indicate that this algorithm may be fast enough for real-time control.

Long, T. W.↗

Damage detection in elastic structures using vibratory residual forces and weighted sensitivity

A methodology is presented for detecting structural damage in elastic structures by nondestructive means. Measured modal test data along with a correlated analytical structural model are used to locate potentially damaged regions using residual modal force vectors and to conduct a weighted sensitivity analysis to assess the extent of mass and/or stiffness variations, where damage is characterized as a stiffness reduction. The current approach is unique among other approaches in that it accounts for (1) variations in system mass, system stiffness, and mass center (locations), (2) perturbations of both the natural frequencies and modal vectors, and (3) statistical confidence factors for the structural parameters and potential experimental instrumentation error. Moreover, this procedure can be used with either full or reduced models. A wide variety of numerical examples are presented that show that the current method provides a precise indication of both the location and the extent of structural damage.

Ricles, J. M.↗

On the application of transonic similarity rules to wings of finite span

The transonic aerodynamic characteristics of wings of finite span are discussed from the point of view of a unified small perturbation theory for subsonic, transonic, and supersonic flows about thin wings. This approach avoids certain ambiguities which appear if one studies transonic flows by means of equations derived under the more restrictive assumption that the local velocities are everywhere close to sonic velocity. The relation between the two methods of analysis of transonic flow is examined, the similarity rules and known solutions of transonic flow theory are reviewed, and the asymptotic behavior of the lift, drag, and pitching-moment characteristics of wings of large and small aspect ratio is discussed. It is shown that certain methods of data presentation are advantageous for the effective display of these characteristics.

Spreiter, John R↗

Effect of the acoustic boundary layer on the wave propagation in ducts

An analysis is presented for the wave propagation in two-dimensional and circular lined ducts taking into account the effects of viscosity in both the mean and the acoustic problems. The method of composite expansions is used to express each acoustic flow quantity as the sum of an inviscid part and a boundary layer part insignificant outside a thin layer next to the wall. The problem is reduced to solving a second-order ordinary differential equation for the pressure perturbation as in the inviscid acoustic case but with a modified specific wall admittance. An analytic expression is presented for the variation of the modified admittance with the wall and flow parameters, such as the acoustic boundary layer thickness, the mean velocity and temperature gradients at the wall, the frequency of oscillation, and the wavelength.

Nayfeh, A. H.↗

Remarks on the theory of aerodynamic noise

The accuracy of approximations employed in the Lighthill theory in current aerodynamic noise research is critically evaluated. Based on the method of matched asymptotic expansions, the full Navier-Stokes equations are expanded for small Mach numbers. To the first order, the near-field is generally nonisentropic. The pressure field (pseudo-sound) is generated by the incompressible Reynolds stresses in the turbulent flow and the velocity, pressure perturbation, and their derivatives on the boundaries. In the far-field, the first-order pressure (acoustic) field satisfying a linear wave equation is obtained by matching with the pseudo-sound field. A uniformly valid solution for the pressure field in a stationary or uniformly moving medium is obtained. The solution shows that the generation of the first-order aerodynamic noise does not depend on the viscous, thermal, or entropy effects in the adiabatic flow nor on the shear stress on a smooth rigid boundary.

Pan, Y. S.↗

A model study of diurnal behavior of the ionosphere and the protonosphere coupling

A new method to study dynamic behavior of the ionosphere-protonosphere coupling within the plasmasphere is developed and used to calculate ion distributions above 500 km using observed electron densities at 500 km. The method is based on the relation between the total magnetic flux tube content of H(+) above some reference height (e.g., 3000 km) and the H(+) flux at that height, which is uniquely determined by the coupled momentum and continuity equations for the quasi-steady state. The O(+) profile is perturbed from diffusive equilibrium only by ion drag with H(+). The time dependency is taken into account by applying self-consistent boundary conditions which are obtained from the solution of the time-dependent equation for the tube content of H(+). Calculations are carried out for two models. The first model deals with the ion dynamics in a tube of force which rotates with the earth. In the second model, effects of the cross-L plasma drifts are considered. Both models produce diurnal variations of H(+) and O(+) densities near 1000 km similar to observational results.

Marubashi, K.↗

Inversion of stratospheric aerosol and gaseous constituents from spacecraft solar extinction data in the 0.38-1.0-micron wavelength region

The paper discusses a possible data retrieval technique for the spaceborne Stratospheric Aerosol and Gas Experiment (SAGE). The SAGE instrument has four radiometric channels located at selected intervals in the 0.38-1.0-micron wavelength range. A data reduction procedure is described for minimization of experimental errors on the basis of a detailed simulation of the measurement sequence. An efficient and accurate inversion method is then used for the retrieval of all the constituent vertical profiles. Also, the effects of horizontally inhomogeneous distributions of the constituent vertical profiles are studied based on available data of their global distributions. A simple horizontally inhomogeneous model of stratospheric aerosol and ozone is employed to estimate the perturbation on the retrieval accuracies.

Chu, W. P.↗

High-frequency measurements of multilayer ceramic capacitors

A resonant coaxial transmission line, short circuited at one end and open circuited at the other, whose fundamental resonant frequency and Q factor are known, is perturbed with a test capacitor connected either in series at the shorted end of the line, or in shunt at the open end. Measuring the Q factor of the system with the delta f technique yields the effective series resistance, capacitance, and the Q factor of the test specimen. This method of measurement has the advantage that there are no adjustable elements to alter circuit conditions in an unprescribed way, the only variable is the frequency which can be measured with an uncertainty of less than 1 ppm, the loss of the line as a function of frequency is quite predictable, and the Q factor of the line can be made sufficiently high to support accurate measurements of low loss capacitors.

Lafferty, R. E.↗

Analytic computation of energy derivatives - Relationships among partial derivatives of a variationally determined function

Formulas are presented relating Taylor series expansion coefficients of three functions of several variables, the energy of the trial wave function (W), the energy computed using the optimized variational wave function (E), and the response function (lambda), under certain conditions. Partial derivatives of lambda are obtained through solution of a recursive system of linear equations, and solution through order n yields derivatives of E through order 2n + 1, extending Puley's application of Wigner's 2n + 1 rule to partial derivatives in couple perturbation theory. An examination of numerical accuracy shows that the usual two-term second derivative formula is less stable than an alternative four-term formula, and that previous claims that energy derivatives are stationary properties of the wave function are fallacious. The results have application to quantum theoretical methods for the computation of derivative properties such as infrared frequencies and intensities.

King, H. F.↗

Accurate collision-induced line-coupling parameters for the fundamental band of CO in He - Close coupling and coupled states scattering calculations

The first accurate theoretical values for off-diagonal (i.e., line-coupling) pressure-broadening cross sections are presented. Calculations were done for CO perturbed by He at thermal collision energies using an accurate ab initio potential energy surface. Converged close coupling, i.e., numerically exact values, were obtained for coupling to the R(0) and R(2) lines. These were used to test the coupled states (CS) and infinite order sudden (IOS) approximate scattering methods. CS was found to be of quantitative accuracy (a few percent) and has been used to obtain coupling values for lines to R(10). IOS values are less accurate, but, owing to their simplicity, may nonetheless prove useful as has been recently demonstrated.

Green, Sheldon↗

Dust particles from comets and asteroids collected at the Earth's orbit: Parent-daughter relationships

The relative contributions of comets and asteroids to the reservoir of dust in the interplanetary medium is not well known. There are direct observations of dust released from comets and there is evidence to associate the IRAS dust bands with possible collisions of Asteroids in the main belt. It is believed that one may combine lab analysis of the physics and chemistry of captured particles with orbital data in order to identify comet and asteroid parent bodies. It is possible to use the collected orbits of the dust to connect with its source in two ways. One is to consider the long time orbit evolution of the dust under Poynting-Robertson drag. The other is to look at the prompt orbit change of dust from comets onto trajectories that intersect the earth's orbit. In order to characterize the orbits of dust particles evolved over a long period of time, a study of its orbital evolution was undertaken. Various parameters associated with these dust orbits as they cross the Earth's orbit were considered in order to see if one may discriminate between particles evolved from comets and asteroids. The method was to calculate by a numerical procedure the orbits of dust particles after they left their parent bodies. It appears that as the particles pass the Earth's orbit, asteroidal grains and cometary grains can be differentiated on the basis of their measured orbital eccentricities even after much planetary perturbation. Broad parent daughter associations can be made on this basis from measurement of their trajectories intercepted in earth orbit.

Jackson, A. A.↗

Use of source distributions for evaluating theoretical aerodynamics of thin finite wings at supersonic speeds

A series of publications on the source-distribution methods for evaluating the aerodynamics of thin wings at supersonic speeds is summarized, extended, and unified. Included in the first part are the deviations of: (a) the linearized partial-differential equation for unsteady flow at a substantially constant Mach number. b) The source-distribution solution for the perturbation-velocity potential that satisfies the boundary conditions of tangential flow at the surface and in the plane of the wing; and (c) the integral equation for determining the strength and the location of sources to describe the interaction effects (as represented by upwash) of the bottom and top wing surfaces through the region between the finite wing boundary and the foremost Mach wave. The second part deals with steady-state thin-wing problems. The third part of the report approximates the integral equation for unsteady upwash and includes a solution of approximate equation. Expressions are then derived to evaluate the load distributions for time-dependent finite-wing motions.

Evvard, John C↗

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 2: Homogeneous medium at Mach 3 with gamma = 5/3

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing spheres of varying sizes (from 10 down to 0.01 accretion radii) move at Mach 3 relative to a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3). To accommodate the long-range gravitational forces, the extent of the computational volume is 32(exp 3) accretion radii. We examine the influence of numerical procedure on physical behavior. The hydrodynamics is modeled by the 'piecewise parabolic method.' No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (5-10) grids around the sphere, each finer grid being a factor of 2 smaller in zone dimension that the next coarser grid. The largest dynamic range (ratio of size of the largest grid to size of the finest zone) is 16,384. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid, while at the same time evolving the gas on the coarser grids. Initially (at time t = 0-10), a shock front is set up, a Mach cone develops, and the accretion column is observable. Eventually the flow becomes unstable, destroying axisymmetry. This happens approximately when the mass accretion rate reaches the values (+/- 10%) predicted by the Bondi-Hoyle accretion formula (factor of 2 included). However, our three-dimensional models do not show the highly dynamic flip-flop flow so prominent in two-dimensional calculations performed by other authors. The flow, and thus the accretion rate of all quantities, shows quasi-periodic (P approximately equals 5) cycles between quiescent and active states. The interpolation formula proposed in an accompanying paper is found to follow the collected numerical data to within approximately 30%. The specific angular momentum accreted is of the same order of magnitude as the values previously found for two-dimensional flows.

Ruffert, Maximilian↗

Convection, nucleosynthesis, and core collapse

We use a piecewise parabolic method hydrodynamics code (PROMETHEUS) to study convective burning in two dimensions in an oxygen shell prior to core collapse. Significant mixing beyond convective boundaries determined by mixing-length theory brings fuel (C-12) into the convective regon, causing hot spots of nuclear burning. Plumes dominate the velocity structure. Finite perturbations arise in a region in which O-16 will be explosively burned to Ni-56 when the star explodes; the resulting instabilities and mixing are likely to distribute Ni-56 throughout the supernova envelope. Inhomogeneities in Y(sub e) may be large enough to affect core collapse and will affect explosive nucleosynthesis. The nature of convective burning is dramatically different from that assumed in one-dimensional simulations; quantitative estimates of nucleosynthetic yields, core masses, and the approach to core collapse will be affected.

Bazan, Grant↗

Expendable Launch Vehicle Studies

Analytical support studies of expendable launch vehicles concentrates on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane linearized dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed. The design of a robust LQR controller based on the reduced order system is accomplished using the parameter perturbation technique. The ELV modeling and analysis team has spent the past three years working on the theoretical development and application of sensitivity analysis to solve eigenvalue problems associated with structural dynamics. Specific application areas include stochastic vibrations, viscously damped vibrations and nonlinear dynamics. In stochastic linear vibrations, sensitivity analysis methods have been developed which determined the expectation and variance of the response for given probability density functions of various stochastic parameters. Examples worked include beam vibrations with up to six stochastic parameters; and the eigensensitivity results differ little from those obtained, with much effect, from Monte Carlo techniques. For viscously damped vibrations, eigensensitivity analysis has given excellent results for homogeneous beams, modeled by either quadratic or quartic eigenvalue equations. Specific applications include Kevin and Maxwell-type viscoelastic beams.

Bainum, P. M.↗