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 703 records · Page 39

Investigation of the effects of external current systems on the MAGSAT data utilizing grid cell modeling techniques

The feasibility of modeling magnetic fields due to certain electrical currents flowing in the Earth's ionosphere and magnetosphere was investigated. A method was devised to carry out forward modeling of the magnetic perturbations that arise from space currents. The procedure utilizes a linear current element representation of the distributed electrical currents. The finite thickness elements are combined into loops which are in turn combined into cells having their base in the ionosphere. In addition to the extensive field modeling, additional software was developed for the reduction and analysis of the MAGSAT data in terms of the external current effects. Direct comparisons between the models and the MAGSAT data are possible.

Klumpar, D. M.↗

Fourier-Legendre spectral methods for incompressible channel flow

An iterative collocation technique is described for modeling implicit viscosity in three-dimensional incompressible wall bounded shear flow. The viscosity can vary temporally and in the vertical direction. Channel flow is modeled with a Fourier-Legendre approximation and the mean streamwise advection is treated implicitly. Explicit terms are handled with an Adams-Bashforth method to increase the allowable time-step for calculation of the implicit terms. The algorithm is applied to low amplitude unstable waves in a plane Poiseuille flow at an Re of 7500. Comparisons are made between results using the Legendre method and with Chebyshev polynomials. Comparable accuracy is obtained for the perturbation kinetic energy predicted using both discretizations.

Zang, T. A.↗

Investigation of dynamic noise affecting geodynamics information in a tethered subsatellite

Work performed as part of an investigation of noise affecting instrumentation in a tethered subsatellite, was studied. The following specific topics were addressed during the reporting period: a method for stabilizing the subsatellite against the rotational effects of atmospheric perturbation was developed; a variety of analytic studies of tether dynamics aimed at elucidating dynamic noise processes were performed; a novel mechanism for coupling longitudinal and latitudinal oscillations of the tether was discovered, and random vibration analysis for modeling the tethered subsatellite under atmospheric perturbation were studied.

Gullahorn, G. E.↗

Aeroelastic analysis of launch vehicles in transonic flight

A numerical study of the aeroelastic stability of launch vehicles in transonic flight is performed. The unsteady aerodynanic flowfield about the vehicle is simulated using the three dimensional Reynolds-Averaged Navier-Stokes equations with the account for turbulent mixing done by a two-layer algebraic eddy viscosity model. The structural-dynamic formulation considers free-free flexural vibration of an elongated beam with variable properties and the equations are cast in modal form. The two sets of equations are integrated simultaneously in time, after a initial perturbation, to ascertain the aeroelastic stability of a given configuration. The method propossed is described in detail, and results that illustrate its application are presented.

Azevedo, Joao Luiz F.↗

Kinetic effects of Alfven wave nonlinearity. I - Ponderomotive density fluctuations

The Vlasov theory is used to study kinetic corrections to fluid descriptions of Alfven wave nonlinearity. The method is to obtain an expression for the second-order perturbed distribution function produced by a nonlinear Alfven wave. From this distribution function a kinetically correct expression is obtained for the plasma density perturbation associated with an envelope-modulated Alfven wave. This kinetic theory result differs substantially from the fluid expression when the plasma beta is greater than about 1, and the electron and ion temperatures are approximately equal. This result is of interest because density fluctuations are an observationally accessible indicator of wave nonlinearity in solar system Alfven waves. It also will assist in the determination of properties of Alfven waves in the interstellar medium. Finally, this analysis also yields a kinetically correct expression for u, the magnetic field-aligned component of the plasma fluid velocity.

Spangler, Steven R.↗

A high level Ab initio study of the anionic hydrogen-bonded complexes FH-CN-, FH-NC-, H2O-CN- and H2O-NC-

HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.

Lee, Timothy J.↗

Signal processing applied to photothermal techniques for materials characterization

There is a need to make noncontact measurements of material characteristics in the microgravity environment. Photothermal and photoacoustics techniques offer one approach for attaining this capability since lasers can be used to generate the required thermal or acoustic signals. The perturbations in the materials that can be used for characterization can be detected by optical reflectance, infrared detection or laser detection of photoacoustics. However, some of these laser techniques have disadvantages of either high energy pulsed excitation or low signal-to-noise ratio. Alternative signal processing techniques that have been developed can be applied to photothermal or photoacoustic instrumentation. One fully coherent spread spectrum signal processing technique is called time delay spectrometry (TDS). With TDS the system is excited using a combined frequency-time domain by employing a linear frequency sweep excitation function. The processed received signal can provide either frequency, phase or improved time resolution. This signal processing technique was shown to outperform other time selective techniques with respect to noise rejection and was recently applied to photothermal instrumentation. The technique yields the mathematical equivalent of pulses yet the input irradiances are orders of magnitude less than pulses with the concomitant reduction in perturbation of the sample and can increase the capability of photothermal methods for materials characterization.

Rooney, James A.↗

A high-level ab initio study of the anionic hydrogen-bonded complexes FH-CN(-), FH-NC(-), H2O-CN(-), and H2O-NC(-)

HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.

Lee, Timothy J.↗

On the formulation of a minimal uncertainty model for robust control with structured uncertainty

In the design and analysis of robust control systems for uncertain plants, representing the system transfer matrix in the form of what has come to be termed an M-delta model has become widely accepted and applied in the robust control literature. The M represents a transfer function matrix M(s) of the nominal closed loop system, and the delta represents an uncertainty matrix acting on M(s). The nominal closed loop system M(s) results from closing the feedback control system, K(s), around a nominal plant interconnection structure P(s). The uncertainty can arise from various sources, such as structured uncertainty from parameter variations or multiple unsaturated uncertainties from unmodeled dynamics and other neglected phenomena. In general, delta is a block diagonal matrix, but for real parameter variations delta is a diagonal matrix of real elements. Conceptually, the M-delta structure can always be formed for any linear interconnection of inputs, outputs, transfer functions, parameter variations, and perturbations. However, very little of the currently available literature addresses computational methods for obtaining this structure, and none of this literature addresses a general methodology for obtaining a minimal M-delta model for a wide class of uncertainty, where the term minimal refers to the dimension of the delta matrix. Since having a minimally dimensioned delta matrix would improve the efficiency of structured singular value (or multivariable stability margin) computations, a method of obtaining a minimal M-delta would be useful. Hence, a method of obtaining the interconnection system P(s) is required. A generalized procedure for obtaining a minimal P-delta structure for systems with real parameter variations is presented. Using this model, the minimal M-delta model can then be easily obtained by closing the feedback loop. The procedure involves representing the system in a cascade-form state-space realization, determining the minimal uncertainty matrix, delta, and constructing the state-space representation of P(s). Three examples are presented to illustrate the procedure.

Belcastro, Christine M.↗

Ab initio characterization of ClOOH - Implications for atmospheric chemistry

The equilibrium structure, dipole moment, harmonic vibrational frequencies, and infrared intensities of ClOOH are determined using the CCSD(T) (singles and doubles coupled-cluster theory plus a perturbational estimate of the effects of connected triple excitations) electronic structure method in conjunction with a TZ2P (triple xi plus double polarization) basis set. The heat of formation of CIOOH is determined (using two different isodesmic reactions) to be +1.5 +/- 1 kcal/mol at 0 K or +0.2 +/- 1 kcal/mol at 298.15 K. Using the computed heat of formation, we examined the stability of ClOOH with respect to the ClO + OH, ClOO + H, and HOO + Cl dissociation limits. Since ClOOH is found to be quite stable, it is argued that the chemistry of ClOOH should be included in any accurate modeling of the stratosphere.

Lee, Timothy J.↗

Guidance law development for aeroassisted transfer vehicles using matched asymptotic expansions

This report addresses and clarifies a number of issues related to the Matched Asymptotic Expansion (MAE) analysis of skip trajectories, or any class of problems that give rise to inner layers that are not associated directly with satisfying boundary conditions. The procedure for matching inner and outer solutions, and using the composite solution to satisfy boundary conditions is developed and rigorously followed to obtain a set of algebraic equations for the problem of inclination change with minimum energy loss. A detailed evaluation of the zeroth order guidance algorithm for aeroassisted orbit transfer is performed. It is shown that by exploiting the structure of the MAE solution procedure, the original problem, which requires the solution of a set of 20 implicit algebraic equations, can be reduced to a problem of 6 implicit equations in 6 unknowns. A solution that is near optimal, requires a minimum of computation, and thus can be implemented in real time and on-board the vehicle, has been obtained. Guidance law implementation entails treating the current state as a new initial state and repetitively solving the zeroth order MAE problem to obtain the feedback controls. Finally, a general procedure is developed for constructing a MAE solution up to first order, of the Hamilton-Jacobi-Bellman equation based on the method of characteristics. The development is valid for a class of perturbation problems whose solution exhibits two-time-scale behavior. A regular expansion for problems of this type is shown to be inappropriate since it is not valid over a narrow range of the independent variable. That is, it is not uniformly valid. Of particular interest here is the manner in which matching and boundary conditions are enforced when the expansion is carried out to first order. Two cases are distinguished-one where the left boundary condition coincides with, or lies to the right of, the singular region, and another one where the left boundary condition lies to the left of the singular region. A simple example is used to illustrate the procedure where the obtained solution is uniformly valid to O(Epsilon(exp 2)). The potential application of this procedure to aeroassisted plane change is also described and partially evaluated.

Calise, Anthony J.↗

Multi-Sensors Observations of Pre-Earthquake Signals. What We Learned from the Great Tohoku Earthquake?

The lessons learned from the Great Tohoku EQ (Japan, 2011) will affect our future observations and an analysis is the main focus of this presentation. Multi-sensors observations and multidisciplinary research is presented in our study of the phenomena preceding major earthquakes Our approach is based on a systematic analysis of several physical and environmental parameters, which been reported by others in connections with earthquake processes: thermal infrared radiation; temperature; concentration of electrons in the ionosphere; radon/ion activities; and atmospheric temperature/humidity [Ouzounov et al, 2011]. We used the Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model, one of several possible paradigms [Pulinets and Ouzounov, 2011] to interpret our observations. We retrospectively analyzed the temporal and spatial variations of three different physical parameters characterizing the state of the atmosphere, ionosphere the ground surface several days before the March 11, 2011 M9 Tohoku earthquake Namely: (i) Outgoing Long wave Radiation (OLR) measured at the top of the atmosphere; (ii) Anomalous variations of ionospheric parameters revealed by multi-sensors observations; and (iii) The change in the foreshock sequence (rate, space and time); Our results show that on March 8th, 2011 a rapid increase of emitted infrared radiation was observed and an anomaly developed near the epicenter with largest value occurring on March 11 at 07.30 LT. The GPS/TEC data indicate an increase and variation in electron density reaching a maximum value on March 8. Starting from this day in the lower ionosphere there was also observed an abnormal TEC variation over the epicenter. From March 3 to 11 a large increase in electron concentration was recorded at all four Japanese ground-based ionosondes, which returned to normal after the main earthquake. We use the Japanese GPS network stations and method of Radio Tomography to study the spatiotemporal structure of ionospheric perturbations, and to distinguish ionospheric responses to processes of EQ preparation against the effects of other factors. The 2-D snapshots of the electron density over Japan showed abnormal increase over the maximum stress during the night, a few hours before the main shock. Our results from recording atmospheric and ionospheric conditions during the earthquake indicate the presence of anomalies in the atmosphere and ionospheres occurring consistently over regions of maximum stress near the epicenter. Due to their long duration (hours and days) and spatial appearance (only over the Sendai region) these results do not appear to be caused by meteorological or magnetic activity. They reveal the existence of atmospheric and ionospheric phenomena occurring prior to the earthquake, which indicates new evidence of a distinct coupling between the lithosphere and atmosphere/ionosphere. Similar results have been reported before the catastrophic events in Chile (M8.8, 2010), Italy (M6.3, 2009) and Sumatra (M9.3, 2004).

Ouzonounov, D.↗

A semi-analytic theory for the motion of a close-earth artificial satellite with drag

A semi-analytic method is used to estimate the decay history/lifetime and to generate orbital ephemeris for close-earth satellites perturbed by the atmospheric drag and earth oblateness due to the spherical harmonics J2, J3, and J4. The theory maintains efficiency through the application of the theory of a method of averaging and employs sufficient numerical emphasis to include a rather sophisticated atmospheric density model. The averaged drag effects with respect to mean anomaly are evaluated by a Gauss-Legendre quadrature while the averaged variational equations of motion are integrated numerically with automatic step size and error control.

Liu, J. J. F.↗

Studies on the interference of wings and propeller slipstreams

The small disturbance potential flow theory is applied to determine the lift of an airfoil in a nonuniform parallel stream. The given stream is replaced by an equivalent stream with a certain number of velocity discontinuities, and the influence of these discontinuities is obtained by the method of images. Next, this method is extended to the problem of an airfoil in a nonuniform stream of smooth velocity profile. This model allows perturbation velocity potential in a rotational undisturbed stream. A comparison of these results with numerical solutions of Euler equations indicates that, although approximate, the present method provides useful information about the interaction problem while avoiding the need to solve the Euler equations.

Prabhu, R. K.↗

Eigenvalue perturbation models for flexible structures

It is pointed out that real parametric modal frequency and damping variation for lightly damped systems does not resemble disks on the complex plane. Using complex uncertainty (and therefore disklike) models can introduce conservativeness in the design method. An alternative means of developing suitable complex uncertainty models is presented. It involves treating the uncertainty as perturbations to the system eigenvalues and using a particular linear fractional transformation to cover this uncertainty. This approach is applicable to modes within the bandwidth of control. It can be used in conjunction with the standard perturbation modeling approaches. A simple SISO (single-input single-output) example, motivated by a flexible truss experiment at the Jet Propulsion Laboratory, is discussed in order to illustrate the proposed approach.

Smith, Roy S.↗

System Being Developed to Measure the Rotordynamic Characteristics of Air Foil Bearings

Because of the many possible advantages of oil-free engine operation, interest in using air lubricated foil-bearing technology in advanced oil-free engine concepts has recently increased. The Oil-Free Turbomachinery Program at the NASA Glenn Research Center at Lewis Field has partially driven this recent push for oil-free technology. The program's goal of developing an innovative, practical, oil-free gas turbine engine for aeropropulsion began with the development of NASA's high-temperature solid-lubricant coating, PS304. This coating virtually eliminates the life-limiting wear that occurs during the startup and shutdown of the bearings. With practically unlimited life, foil air bearings are now very attractive to rotating machinery designers for use in turbomachinery. Unfortunately, the current knowledge base of these types of bearings is limited. In particular, the understanding of how these types of bearings contribute to the rotordynamic stability of turbomachinery is insufficient for designers to design with confidence. Recent work in oil-free turbomachinery has concentrated on advancing the understanding of foil bearings. A high-temperature fiber-optic displacement probe system and measurement method were developed to study the effects of speed, load, temperature, and other environmental issues on the stiffness characteristics of air foil bearings. Since high temperature data are to be collected in future testing, the testing method was intentionally simplified to minimize the need for expensive test hardware. The method measures the displacement induced upon a bearing in response to an applied perturbation load. The early results of these studies, which are shown in the accompanying figure, indicate trends in steady state stiffness that suggest stiffness increases with load and decreases with speed. It can be seen, even from these data, that stiffness is not expected to change by orders of magnitude over the normal operating range of most turbomachinery; a promising sign for their eventual integration into oil-free turbomachines. Planned future testing will generate similar plots for stiffness changes with temperature and geometry, as well as damping data. The data collected by this method represent a critical step toward understanding how to successfully apply foil air bearings to future oil-free turbomachinery systems.

Howard, Samuel A.↗

Application of a one-strip integral method to the unsteady supersonic aerodynamics of an inclined flat surface

The method of integral relations is applied in a one-strip approximation to the perturbation equations governing small motions of an inclined, sharp-edged, flat surface about the mean supersonic steady flow. Algebraic expressions for low reduced-frequency aerodynamics are obtained and a set of ordinary differential equations are obtained for general oscillatory motion. Results are presented for low reduced-frequency aerodynamics and for the variation of the unsteady forces with frequency. The method gives accurate results for the aerodynamic forces at low reduced frequency which are in good agreement with available experimental data. However, for cases in which the aerodynamic forces vary rapidly with frequency, the results are qualitatively correct, but of limited accuracy. Calculations indicate that for a range of inclination angles near shock detachment such that the flow in the shock layer is low supersonic, the aerodynamic forces vary rapidly both with inclination angle and with reduced frequency.

Bennett, R. M.↗

Optimal perturbation models for averaged orbit generation

Averaging techniques applied to the variation of parameters (VOP) formulation of the equations of motion are being investigated as methods for long-term prediction of artificial satellite orbits. Analytically averaged equations were compared with numerically averaged equations with respect to accuracy and efficiency for computation of zonal and nonresonant third-body perturbations. Numerically averaged equations were also evaluated for computation of long-period effects from resonant third-body, tesseral harmonic, and atmospheric drag perturbations. Guidelines will be presented for application of averaged VOP equations to a broad class of orbits.

Long, A. C.↗