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,171 records · Page 65

Substructure Versus Property-Level Dispersed Modes Calculation

This paper calculates the effect of perturbed finite element mass and stiffness values on the eigenvectors and eigenvalues of the finite element model. The structure is perturbed in two ways: at the "subelement" level and at the material property level. In the subelement eigenvalue uncertainty analysis the mass and stiffness of each subelement is perturbed by a factor before being assembled into the global matrices. In the property-level eigenvalue uncertainty analysis all material density and stiffness parameters of the structure are perturbed modified prior to the eigenvalue analysis. The eigenvalue and eigenvector dispersions of each analysis (subelement and property-level) are also calculated using an analytical sensitivity approximation. Two structural models are used to compare these methods: a cantilevered beam model, and a model of the Space Launch System. For each structural model it is shown how well the analytical sensitivity modes approximate the exact modes when the uncertainties are applied at the subelement level and at the property level.

Stewart, Eric C.↗

Sensitivity analysis for aeroacoustic and aeroelastic design of turbomachinery blades

A new method for computing the effect that small changes in the airfoil shape and cascade geometry have on the aeroacoustic and aeroelastic behavior of turbomachinery cascades is presented. The nonlinear unsteady flow is assumed to be composed of a nonlinear steady flow plus a small perturbation unsteady flow that is harmonic in time. First, the full potential equation is used to describe the behavior of the nonlinear mean (steady) flow through a two-dimensional cascade. The small disturbance unsteady flow through the cascade is described by the linearized Euler equations. Using rapid distortion theory, the unsteady velocity is split into a rotational part that contains the vorticity and an irrotational part described by a scalar potential. The unsteady vorticity transport is described analytically in terms of the drift and stream functions computed from the steady flow. Hence, the solution of the linearized Euler equations may be reduced to a single inhomogeneous equation for the unsteady potential. The steady flow and small disturbance unsteady flow equations are discretized using bilinear quadrilateral isoparametric finite elements. The nonlinear mean flow solution and streamline computational grid are computed simultaneously using Newton iteration. At each step of the Newton iteration, LU decomposition is used to solve the resulting set of linear equations. The unsteady flow problem is linear, and is also solved using LU decomposition. Next, a sensitivity analysis is performed to determine the effect small changes in cascade and airfoil geometry have on the mean and unsteady flow fields. The sensitivity analysis makes use of the nominal steady and unsteady flow LU decompositions so that no additional matrices need to be factored. Hence, the present method is computationally very efficient. To demonstrate how the sensitivity analysis may be used to redesign cascades, a compressor is redesigned for improved aeroelastic stability and two different fan exit guide vanes are redesigned for reduced downstream radiated noise. In addition, a framework detailing how the two-dimensional version of the method may be used to redesign three-dimensional geometries is presented.

Lorence, Christopher B.↗

Gaussian variational equations for osculating elements of an arbitrary separable reference orbit.

Lagrange-type equations are often used in planetary theory and sometimes in satellite theory. The equations express the variation of osculating Keplerian elements in terms of derivatives of a disturbing function or a perturbing potential. When the perturbing force is derivable from a potential, it is possible to convert the Lagrange-type equations to another form. This form, usually, attributed to Gauss, contains the perturbing forces instead of the derivatives of the potential. The case of an arbitrary separable reference orbit is discussed together with a lemma, a Keplerian check, and questions of the applicability of the equations to the spheroidal method.

Vinti, J. P.↗

Effect of the properties of a porous coating on boundary layer stability

Drawing off gas from the boundary layer is a well-known method for increasing the stability of boundary layers. The increase in stability is primarily connected with a change in the velocity profile form in the case of suction. On the basis of the assumption that the velocity perturbations on a porous slate do not equal zero, the influence of the properties of a permeable surface upon the boundary layer stability were studied.

Gaponov, S. A.↗

A mean spherical model for soft potentials: The hard core revealed as a perturbation

The mean spherical approximation for fluids is extended to treat the case of dense systems interacting via soft-potentials. The extension takes the form of a generalized statement concerning the behavior of the direct correlation function c(r) and radial distribution g(r). From a detailed analysis that views the hard core portion of a potential as a perturbation on the whole, a specific model is proposed which possesses analytic solutions for both Coulomb and Yukawa potentials, in addition to certain other remarkable properties. A variational principle for the model leads to a relatively simple method for obtaining numerical solutions.

Rosenfeld, Y.↗

Mean-spherical model for soft potentials - The hard core revealed as a perturbation

The mean-spherical approximation for fluids is extended to treat the case of dense systems interacting via soft potentials. The extension takes the form of a generalized statement concerning the behavior of the direct-correlation function c(r) and the radial-distribution function g(r). From a detailed analysis that views the hard-core portion of a potential as a perturbation on the whole, a specific model is proposed which possesses analytic solutions for both Coulomb and Yukawa potentials, in addition to certain other remarkable properties. A variational principle for the model leads to a relatively simple method for obtaining numerical solutions.

Rosenfeld, Y.↗

A rapid method for the approximate determination of nonlinear solutions Application to aerodynamic flows

A method for determining highly accurate approximations to families of strongly nonlinear solutions which are either continuous or discontinuous, and which represent variations in some arbitrary parameters, is developed and evaluated. The procedure consists of defining a unit perturbation by employing two or more nonlinear solutions which differ from one another by a nominal change in some geometric or flow parameter, and then using that unit perturbation to predict a family of related nonlinear solutions over a range of parameter variation. Coordinate straining is used in determining the unit perturbation to account for the movement of discontinuities and maxima of high-gradient regions due to the perturbation. Although the procedure is generally applicable, results are presented here for nonlinear aerodynamic applications. Attention is focused in particular on transonic flows which are strongly supercritical and exhibit large surface shock movement over the parametric range studied; and on subsonic flows which display large pressure variations in the stagnation and peak suction pressure regions. Flows past both isolated airfoils and compressor cascades involving a variety of flow and geometry parameter changes are considered. Comparisons with the corresponding 'exact' nonlinear solutions indicate a remarkable accuracy and range of validity of such a procedure. Computational time is trivial.

Stahara, S. S.↗

A model for unsteady transonic indicial responses

A method of characterizing the indicial curve for the response frequencies of flutter in the transonic regime by means of a time scale parameter is illustrated. Time linearization of the unsteady potential equation is assumed adequate and the perturbed potential about a steady flow is defined with appropriate boundary conditions. An indicial motion can then be obtained for a given mode, with the lift and drag coefficients which result from a step change in incidence or control surface deflection available to describe any kind of motion. The lift and drag coefficients are shown to result analytically by addition of the multiplicand involving the time scale parameter. If the coefficients are maximized, then at reduced frequencies only two parameters are necessary to form the indicial response curve.

Fung, K.-Y.↗

Macroscopic Lagrangian description of warm plasmas. II Nonlinear wave interactions

A macroscopic Lagrangian is simplified to the adiabatic limit and expanded about equilibrium, to third order in perturbation, for three illustrative cases: one-dimensional compression parallel to the static magnetic field, two-dimensional compression perpendicular to the static magnetic field, and three-dimensional compression. As examples of the averaged-Lagrangian method applied to nonlinear wave interactions, coupling coefficients are derived for interactions between two electron plasma waves and an ion acoustic wave, and between an ordinary wave, an electron plasma wave, and an ion acoustic wave.

Kim, H.↗

Overview of magnetic bearing control and linearization approaches for annular magnetically suspended devices

An overview of magnetic bearing control and linearization approaches which have been considered for annular magnetically suspended devices is presented. These devices include the Annular Momentum Control Device and the Annular Suspension and Pointing System. Two approaches were investigated for controlling the magnetic actuator. One approach involves controlling the upper and lower electromagnets differentially about a bias flux. The bias flux can either be supplied by permanent magnets in the magnetic circuit or by bias currents. In the other approach, either the upper electromagnet or the lower electromagnet is controlled depending on the direction of force required. One advantage of the bias flux is that for small gap perturbations about a fixed operating point, the force-current characteristic is linear. Linearization approaches investigated for individual element control include an analog solution of the nonlinear electromagnet force equation and a microprocessor-based table lookup method.

Groom, N. J.↗

Static and dynamic characteristics of parallel-grooved seals

Presented is an analytical method to determine static and dynamic characteristics of annular parallel-grooved seals. The governing equations were derived by using the turbulent lubrication theory based on the law of fluid friction. Linear zero- and first-order perturbation equations of the governing equations were developed, and these equations were analytically investigated to obtain the reaction force of the seals. An analysis is presented that calculates the leakage flow rate, the torque loss, and the rotordynamic coefficients for parallel-grooved seals. To demonstrate this analysis, we show the effect of changing number of stages, land and groove width, and inlet swirl on stability of the boiler feed water pump seals. Generally, as the number of stages increased or the grooves became wider, the leakage flow rate and rotor-dynamic coefficients decreased and the torque loss increased.

Iwatsubo, Takuzo↗

Frequency-shaped large-angle maneuvers

The paper considers the problem of maneuvering a flexible spacecraft through large angles in finite time. The basic control problem is divided into two parts. The first part consists of generating a frequency-shaped open-loop solution for the nonlinear rigid body as the nominal solution. The resulting two-point boundary-value problem is solved by introducing a continuation method for altering the mass distribution and boundary conditions for the spacecraft. For the second part, a feedback control is designed by linearizing the flexible body response about several points along the rigid body nominal solution. The perturbation gains are designed by using a frequency-shaped cost functional approach. The gains are linearly interpolated to produce smooth control time-histories as the linear piecewise constant plant models change during the maneuver.

Chun, Hon M.↗

Robust control with structured perturbations

Two important problems in the area of control systems design and analysis are discussed. The first is the robust stability using characteristic polynomial, which is treated first in characteristic polynomial coefficient space with respect to perturbations in the coefficients of the characteristic polynomial, and then for a control system containing perturbed parameters in the transfer function description of the plant. In coefficient space, a simple expression is first given for the l(sup 2) stability margin for both monic and non-monic cases. Following this, a method is extended to reveal much larger stability region. This result has been extended to the parameter space so that one can determine the stability margin, in terms of ranges of parameter variations, of the closed loop system when the nominal stabilizing controller is given. The stability margin can be enlarged by a choice of better stabilizing controller. The second problem describes the lower order stabilization problem, the motivation of the problem is as follows. Even though the wide range of stabilizing controller design methodologies is available in both the state space and transfer function domains, all of these methods produce unnecessarily high order controllers. In practice, the stabilization is only one of many requirements to be satisfied. Therefore, if the order of a stabilizing controller is excessively high, one can normally expect to have a even higher order controller on the completion of design such as inclusion of dynamic response requirements, etc. Therefore, it is reasonable to have a lowest possible order stabilizing controller first and then adjust the controller to meet additional requirements. The algorithm for designing a lower order stabilizing controller is given. The algorithm does not necessarily produce the minimum order controller; however, the algorithm is theoretically logical and some simulation results show that the algorithm works in general.

Keel, Leehyun↗

Galaxy interactions and strength of nuclear activity

Analysis of data in the literature for differential velocities and projected separations of nearby Seyfert galaxies with possible companions shows a clear difference in projected separations between type 1's and type 2's. This kinematic difference between the two activity classes reinforces other independent evidence that their different nuclear characteristics are related to a non-nuclear physical distinction between the two classes. The differential velocities and projected separations of the galaxy pairs in this sample yield mean galaxy masses, sizes, and mass to light ratios which are consistent with those found by the statistical methods of Karachentsev. Although the galaxy sample discussed here is too small and too poorly defined to provide robust support for these conclusions, the results strongly suggest that nuclear activity in Seyfert galaxies is associated with gravitational perturbations from companion galaxies, and that there are physical distinctions between the host companions of Seyfert 1 and Seyfert 2 nuclei which may depend both on the environment and the structure of the host galaxy itself.

Simkin, S. M.↗

A linear programming approach to characterizing norm bounded uncertainty from experimental data

The linear programming spectral overbounding and factorization (LPSOF) algorithm, an algorithm for finding a minimum phase transfer function of specified order whose magnitude tightly overbounds a specified nonparametric function of frequency, is introduced. This method has direct application to transforming nonparametric uncertainty bounds (available from system identification experiments) into parametric representations required for modern robust control design software (i.e., a minimum-phase transfer function multiplied by a norm-bounded perturbation).

Scheid, R. E.↗

Theoretical Study of the First Transition Row Oxides and Sulfides

The first transition row oxides and sulfides are studied using several different levels of theory. The calculations show the bonding mechanism in the sulfides and oxides to be very similar. For the oxides, accurate experimental data allow the theoretical methods to be calibrated. The same level of theory is used to study the sulfides where there is far less experimental information. For ScO through MnO and CuO the coupled cluster singles and doubles technique including a perturbational estimate of the unliked triple excitations [CCSD(T)] yields spectroscopic constants ((tau)e, (omega)e, and D0) in good agreement with experiment. The triple excitations are found to be very important in achieving this accuracy. For FeO to NiO, the single determinant self-consistent-field (SCF) approach yields pi orbitals that are localized on the metal or oxygen. This appears to cause problems for the single reference techniques; this is discussed in detail for NiO. The complete-active-space SCF/internally contracted averaged coupled pair functional approach (CASSCF/ICACPF) works well for FeO to NiO. The calculation of accurate dipole moments is found to be very difficult.

Bauschlicher, Charles W., Jt.↗

Theoretical Study of Infrared and Raman Spectra of Hydrated Magnesium Sulfate Salts

Harmonic and anharmonic vibrational frequencies, as well as infrared and Raman intensities, are calculated for MgSO4.nH20 (n=1-3). Electronic structure theory at the second order Moller-Plesset perturbation theory (MP2) level with a triple-zeta + polarization (TZP) basis set is used to determine the geometry, properties, and vibrational spectra of pure and hydrated MgSO4 salts. The direct vibrational self-consistent field (VSCF) method and its correlation corrected (CC-VSCF) extension are used to determine anharmonic corrections to vibrational frequencies and intensities for the pure MgSO4 and its complex with one water molecule. Very significant differences are found between vibrational of water molecules in complexes with MgSO4 and pure water. Some of the O-H stretching frequencies are shifted to the red very significantly (by up to 1500-2000/cm) upon complexation with magnesium sulfate. They should be observed between 1700 and 3000/cm in a region very different from the corresponding O-H stretch frequency region of pure water (3700-3800/cm). In addition, the SO2 stretching vibrations are found at lower frequency regions than the water vibrations. They can serve as unique identifiers for the presence of sulfate salts. The predicted infrared and Raman spectra should be of valuable help in the design of future missions and analysis of observed data from the ice surface of Jupiter's moon Europa that possibly contains hydrated MgSO4 salts.

Chaban, Galina M.↗

Inertial navigation systems analysis.

This volume offers the avionic systems engineer a fundamental exposition of the mechanization and error analysis of inertial navigation systems. While the material is applicable to spacecraft and undersea navigation, emphasis is placed upon terrestrial applications on or slightly above the earth's surface. As a result, practical considerations are geared toward those aircraft navigation systems of particular current interest. Extensive use is made of perturbation techniques to develop linearized system equations, whose solutions closely approximate those obtained by nonlinear differential equations. A unified error analysis technique is developed that is applicable to virtually all system configurations. The technique provides a greatly simplified method for comparing the performance of competing system configurations.

Britting, K. R.↗