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 667 records · Page 37

Subsonic and supersonic indicial aerodynamics and aerodynamic transfer function for complex configurations

A general theory for indicial-potential-compressible aerodynamics around complex configurations is presented. The motion is assumed to consist of constant subsonic or supersonic speed (steady state) and small perturbations around the steady state. Using the finite-element method to discretize the space problem, a set of differential-difference equations in time relating the potential to its normal derivative on the surface of the body was obtained. The aerodynamics transfer function was derived by using standard method of operational calculus.

Morino, L.↗

Stability of a viscous fluid in a rectangular cavity in the presence of a magnetic field

The stability of an electrically conducting fluid subjected to two dimensional disturbance was investigated. A physical system consisting of two parallel infinite vertical plates which are thermally insulated was studied. An external magnetic field of constant strength was applied to normal plates. The fluid was heated from below so that a steady temperature gradient was maintained in the fluid. The governing equations were derived by perturbation technique, and solutions were obtained by a modified Galerkin method. It was found that the presence of the magnetic field increases the stability of the physical system and instability can occur in the form of neutral or oscillatory instability.

Liang, C. Y.↗

A computer program for wing subsonic aerodynamic performance estimates including attainable thrust and vortex lift effects

Numerical methods incorporated into a computer program to provide estimates of the subsonic aerodynamic performance of twisted and cambered wings of arbitrary planform with attainable thrust and vortex lift considerations are described. The computational system is based on a linearized theory lifting surface solution which provides a spanwise distribution of theoretical leading edge thrust in addition to the surface distribution of perturbation velocities. The approach used relies on a solution by iteration. The method also features a superposition of independent solutions for a cambered and twisted wing and a flat wing of the same planform to provide, at little additional expense, results for a large number of angles of attack or lift coefficients. A previously developed method is employed to assess the portion of the theoretical thrust actually attainable and the portion that is felt as a vortex normal force.

Carlson, H. W.↗

Three-dimensional baroclinic instability of a Hadley cell for small Richardson number

A three-dimensional, linear stability analysis of a baroclinic flow for Richardson number, Ri, of order unity is presented. The model considered is a thin horizontal, rotating fluid layer which is subjected to horizontal and vertical temperature gradients. The basic state is a Hadley cell which is a solution of the complete set of governing, nonlinear equations and contains both Ekman and thermal boundary layers adjacent to the rigid boundaries; it is given in a closed form. The stability analysis is also based on the complete set of equations; and perturbation possessing zonal, meridional, and vertical structures were considered. Numerical methods were developed for the stability problem which results in a stiff, eighth-order, ordinary differential eigenvalue problem. The previous work on three-dimensional baroclinic instability for small Ri was extended to a more realistic model involving the Prandtl number, sigma, and the Ekman number, E, and to finite growth rates and a wider range of the zonal wavenumber.

Antar, B. N.↗

Toward the routine aerodynamic analysis of complex configurations

An evaluation is made of the ability of the TranAir CFD code to routinely compute the aerodynamic characteristics of complex subsonic and supersonic aircraft configurations. TranAir solves the full-potential equation for transonic flow about completely arbitrary geometries, using the surface-paneling PanAir technique in geometry definition. The uniform global grid may be locally refined in regions where flow properties are rapidly changing, such as regions where shocks arise, and around wing leading edges. Unlike panel method codes, TranAir solutions are not undermined by small-perturbation assumptions. Illustrative results are presented for such configurations as the F-16A with wingtip-mounted missiles and underwing fuel tanks, a generic fighter configuration, and a model of NASA-Ames' 12-ft Pressure Wind Tunnel.

Madson, Michael D.↗

Matched asymptotic expansion of the Hamilton-Jacobi-Bellman equation for aeroassisted plane-change maneuvers

In this paper we develop a general procedure for constructing a matched asymptotic expansion of the Hamilton-Jacobi-Bellman equation based on the method of characteristics. The development is for a class of perturbation problems whose solution exhibits two-time-scale behavior. A regular expansion for problems of this type is inappropriate since it is not uniformly valid over a narrow range of the independent variable. 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 one where the left boundary condition lies to the left of the singular region. A simple example is used to illustrate the procedure, and its potential application to aeroassisted plane change is described.

Calise, Anthony J.↗

Numerical derivative techniques for trajectory optimization

The adoption of robust numerical optimization techniques in trajectory simulation programs has resulted in powerful design and analysis tools. These trajectory simulation/optimization programs are widely used, and a representative list includes the GTS system, the POST program, and newer collocation methods such as OTIS and FONPAC. All of these programs rely on optimization algorithms which require objective function and constraint gradient data during the iteration process. However, most trajectory optimization problems lack simple analytical expressions for these derivatives. In the general case a function evaluation involves integrating aerodynamic, propulsive, and gravity forces over multiple trajectory phases with complex control models. With the newer collocation methods, the integration is replaced by defect constraints and cubic approximations for the state. While analytic gradient expressions can sometimes be derived for trajectory optimization problems, the derivation is cumbersome, time consuming, and prone to mistakes. Fortunately, an alternate method exists for the gradient evaluation, namely finite difference approximations. In this paper some finite difference gradient techniques developed for use with the GTS system are presented. These techniques include methods for computing first and second partial derivatives of single and multiple sets of functions. A key feature of these methods is an error control mechanism which automatically adjusts the perturbation size to obtain accurate derivative values.

Hallman, Wayne P.↗

Accurate Quartic Force Fields and Vibrational Frequencies for HCN and HNC

The quartic force fields of HCN and HNC are determined using atomic natural orbital one-particle basis sets of spdf/spd and spdfg/spdf quality in conjunction with the CCSD(T) electron correlation method (singles and doubles coupled-cluster theory plus a perturbational estimate of the effects of connected triple excitations). The HCN force field is in good agreement with a recent experimentally derived force field and also with the force field recently computed by Wong and Bacskay. On the basis of the good agreement obtained for HCN, it is argued that the ab initio quartic force field for HNC is superior to a prior force field derived from experiment. The harmonic frequencies of HNC are predicted to be 3822 +/- 10,472 +/- 5, and 2051 +/-10/cm for omega1, omega2, and omega3, respectively; the experimentally derived values are above these values and fall outside the estimated uncertainties. Using the quartic force field, spectroscopic constants are predicted for HNC based on a vibrational second-order perturbation theory analysis. It is also asserted that the gas-phase fundamental v(sub 3) for HNC is slightly lower than the matrix isolation value. The range of validity of the quartic force fields is investigated by comparison of variational vibrational energies computed with the quartic force fields to those obtained from our recently reported global HCN/HNC potential energy surface and also to experimental data.

Lee, Timothy J.↗

Accurate Quartic Force Fields and Vibrational Frequencies for HCN and HNC

The quartic force fields of HCN and HNC are determined using atomic natural orbital one-particle basis sets of spdf/spd and spdfg/spdf quality in conjunction with the CCSD(T) electron correlation method (singles and doubles coupled-cluster theory plus a perturbation estimate of the effects of connected triple excitations). The HCN force field is in good agreement with a recent experimentally derived force field and also with the force field recently computed by Wong and Bacskay. On the basis of the good agreement obtained for HCN, it is argued that the ab initio quartic force field for HNC is superior to a prior force field derived from experiment. The harmonic frequencies of HNC are predicted to be 3822 +/- 10, 472 +/- 5, and 2051 +/- 10 cm(exp -1) for omega(sub 1), omega(sub 2), and omega(sub 3), respectively; the experimentally derived values are above these values and fall outside the estimated uncertainties. Using the quartic force field, spectroscopic constants are predicted for HNC based on a vibrational second-order perturbation theory analysis. It is also asserted that the gas-phase fundamental nu(sub 3) for HNC is slightly lower than the matrix isolation value. The range of validity of the quartic force fields is investigated by comparison of variational vibrational energies computed with the quartic force fields to those obtained from our recently reported global HCN/HNC potential energy surface and also to experimental data.

Lee, Timothy J.↗

Critical Analyses of Data Differences Between FNMOC and AFGWC Spawned SSM/I Datasets

Antenna temperatures and the corresponding geolocation data from the five sources of the Special Sensor Microwave/Imager data from the Defense Meteorological Satellite Program F11 satellite have been characterized. Data from the Fleet Numerical Meteorology and Oceanography Center (FNMOC) have been compared with data from other sources to define and document the differences resulting from different processing systems. While all sources used similar methods to calculate antenna temperatures, different calibration averaging techniques and other processing methods yielded temperature differences. Analyses of the geolocation data identified perturbations in the FNMOC and National Environmental Satellite, Data and Information Service data. The effects of the temperature differences were examined by generating rain rates using the Goddard Scattering Algorithm. Differences in the geophysical precipitation products are directly attributable to antenna temperature differences.

Ritchie, Adrian A., Jr.↗

Dynamic Stability of Uncertain Laminated Beams Under Subtangential Loads

Because of the inherent complexity of fiber-reinforced laminated composites, it can be challenging to manufacture composite structures according to their exact design specifications, resulting in unwanted material and geometric uncertainties. In this research, we focus on the deterministic and probabilistic stability analysis of laminated structures subject to subtangential loading, a combination of conservative and nonconservative tangential loads, using the dynamic criterion. Thus a shear-deformable laminated beam element, including warping effects, is derived to study the deterministic and probabilistic response of laminated beams. This twenty-one degrees of freedom element can be used for solving both static and dynamic problems. In the first-order shear deformable model used here we have employed a more accurate method to obtain the transverse shear correction factor. The dynamic version of the principle of virtual work for laminated composites is expressed in its nondimensional form and the element tangent stiffness and mass matrices are obtained using analytical integration The stability is studied by giving the structure a small disturbance about an equilibrium configuration, and observing if the resulting response remains small. In order to study the dynamic behavior by including uncertainties into the problem, three models were developed: Exact Monte Carlo Simulation, Sensitivity Based Monte Carlo Simulation, and Probabilistic FEA. These methods were integrated into the developed finite element analysis. Also, perturbation and sensitivity analysis have been used to study nonconservative problems, as well as to study the stability analysis, using the dynamic criterion.

Goyal, Vijay K.↗

Transition Detection at Cryogenic Temperatures Using a Carbon-Based Resistive Heating Layer Coupled with Temperature Sensitive Paint

This paper will highlight the development and application of a carbon-based resistive heating layer for use in transition detection at cryogenic temperatures at the National Transonic Facility (NTF) for full-flight Reynolds number testing. This study builds upon previous work that was successfully demonstrated at the 0.3-m Transonic Cryogenic Tunnel on a smaller-scale airfoil shape of regular geometry. However, the test performed at the NTF involved a semispan wing with complex geometry and significantly larger than previous tests. This required the development of new coatings to provide suitable resistances to provide adequate heating rates for transition detection. Successful implementation of this technology has the ability to greatly enhance transition detection experiments at cryogenic temperatures as well as reducing perturbation in the tunnel caused by more traditional transition detection methods.

Watkins, A. Neal↗

Applications of direct numerical simulation of turbulence in second order closures

This paper discusses two methods of developing models for the rapid pressure-strain correlation term in the Reynolds stress transport equation using direct numerical simulation (DNS) data. One is a perturbation about isotropic turbulence, the other is a perturbation about two-component turbulence -- an extremely anisotropic turbulence. A model based on the latter method is proposed and is found to be very promising when compared with DNS data and other models.

Shih, Tsan-Hsing↗

Perturbation analysis of internal balancing for lightly damped mechanical systems with gyroscopic and circulatory forces

Approximate expressions are developed for internally balanced singular values corresponding to the modes of mechanical systems with gyroscopic forces, light damping, and small circulatory forces. A brief overview is first given of the balanced realization model reduction method, including a discussion of recent work. The models considered are defined, and a perturbation analysis is used to show that the modal representation becomes asymptotically balanced as damping reduces to zero. The approximate balanced singular values are calculated, and a simple example of a flexible, dual-spin spacecraft is given as an illustration of the results.

Blelloch, P. A.↗

An algorithm for the empirical optimization of antenna arrays

A numerical technique is presented to optimize the performance of arbitrary antenna arrays under realistic conditions. An experimental-computational algorithm is formulated in which n-dimensional minimization methods are applied to measured data obtained from the antenna array. A numerical update formula is used to induce partial derivative information without requiring special perturbations of the array parameters. The algorithm provides a new design for the antenna array, and the method proceeds in an iterative fashion. Test case results are presented showing the effectiveness of the algorithm.

Blank, S.↗

Steady subsonic flow around finite-thickness wings

The general method for analyzing steady subsonic potential aerodynamic flow around a lifting body having arbitrary shape is presented. By using the Green function method, an integral representation for the potential is obtained. Under small perturbation assumption, the potential at any point, P, in the field depends only upon the values of the potential and its normal derivative on the surface of the body. Hence if the point P approaches the surface of the body, the representation reduces to an integral equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface. The question of uniqueness is examined and it is shown that, for thin wings, the operator becomes singular as the thickness approaches zero. This fact may yield numerical problems for very thin wings. However, numerical results obtained for a rectangular wing in subsonic flow show that these problems do not appear even for thickness ratio tau = .001. Comparison with existing results shows that the proposed method is at least as fast and accurate as the lifting surface theories.

Kuo, C. C.↗

An analytical evaluation of mascon effects in a semi-analytic theory

A semianalytic method to predict the orbit of a spacecraft orbiting a planet and perturbed by oblateness and a gravity anomaly (mascon) is presented. The Hamiltonian is first expressed in closed form in terms of the nonsingular canonical set of Poincare elements where the mascon is represented by a point-mass. The planet global potential is limited to the second and third harmonics. The mascon potential in harmonics is limited to second degree and first order. The original Lie-Hori method is then applied to transform the canonical equations of motion from the osculating element space to the mean element space. Numerical integration of the new equations is used.

Salama, Ahmed H.↗