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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Scattering from dielectric structures above impedance surfaces and resistive sheets

Interest in understanding of electromagnetic interaction with rough surfaces has prompted the study of scattering from typical dielectric humps over impedance surfaces. It is shown that the Green's function of the problem for a resistive sheet resembles that of the impedance surface. Hence both problems are considered here. A numerical solution for the scattered field of a two-dimensional dielectric object, possibly inhomogeneous, with arbitrary cross section above the impedance surface or resistive sheet is sought. First the Green's function of the problem is derived based on the exact image theory. This form of the Green's function is amenable to numerical computation. Then the induced polarization currents are calculated by casting the integral equations into a matrix equation via the method of moments. Numerical problems in calculation of the Green's function when both source and observation points are close to the surface are discussed. Comparison of numerical results with a perturbation solution shows excellent agreement between the two methods.

Sarabandi, Kamal↗
NASA NTRS1977

A perturbation solution for interlaminar stresses in bidirectional laminates

The paper considers a thin elastic bidirectional symmetric composite laminate subjected to a uniform axial strain. All interlaminar stresses are determined by a zeroth-order perturbation analysis of the governing equations. The method of solution considers the laminate to be composed of an interior region (removed from the free edge) and a boundary layer region. The interior region is analyzed by dropping terms related to powers of the thickness-to-width ratio in the dimensionless elasticity equations. The boundary layer region is studied by introducing a stretching transformation to the governing equations. A uniform solution is then formed by satisfying Prandtl's matching principle of perturbation theory. Results for a four-layer bidirectional graphite-epoxy laminate are compared with existing finite-difference results. It is shown that the perturbation analysis provides more insight into the stress distribution near the free edge.

Hsu, P. W.↗
NASA NTRS1955

Stress analysis of circular semimonocoque cylinders with cutouts

A method is presented for analyzing the stresses about cutouts in circular semimonocoque cylinders with flexible rings. The method involves the use of so-called perturbation stress distributions which are superposed on the stress distribution that would exist in the structure with no cutout in such a way as to give the effects of a cutout. The method can be used for any loading case for which the structure without the cutout can be analyzed and is sufficiently versatile to account for stringer and shear reinforcement about the cutout.

Mccomb, Harvey G , Jr↗
NASA NTRS1982

A general method to determine the stability of compressible flows

Several problems were studied using two completely different approaches. The initial method was to use the standard linearized perturbation theory by finding the value of the individual small disturbance quantities based on the equations of motion. These were serially eliminated from the equations of motion to derive a single equation that governs the stability of fluid dynamic system. These equations could not be reduced unless the steady state variable depends only on one coordinate. The stability equation based on one dependent variable was found and was examined to determine the stability of a compressible swirling jet. The second method applied a Lagrangian approach to the problem. Since the equations developed were based on different assumptions, the condition of stability was compared only for the Rayleigh problem of a swirling flow, both examples reduce to the Rayleigh criterion. This technique allows including the viscous shear terms which is not possible in the first method. The same problem was then examined to see what effect shear has on stability.

Guenther, R. A.↗
NASA NTRS1982

Numerical comparisons of nonlinear convergence accelerators

As part of a continuing program of numerical tests of convergence accelerators, the iterated Aitken's Delta-squared method, Wynn's epsilon algorithm, Brezinski's theta algorithm, and Levin's u transform are compared on a broad range of test problems: linearly convergence alternating, monotone, and irregular-sign series, logarithmically convergent series, power method and Bernoulli method sequences, alternating and monotone asymptotic series, and some perturbation series arising in applications. In each category either the epsilon algorithm or the u transform gives the best results of the four methods tested. In some cases differences among methods are slight, and in others they are quite striking.

Smith, D. A.↗
NASA NTRS1990

The effects of forcing on a single stream shear layer and its parent boundary layer

Forcing and its effect on fluid flows has become an accepted tool in the study and control of flow systems. It has been used both as a diagnostic tool, to explore the development and interaction of coherent structures, and as a method of controlling the behavior of the flow. A number of forcing methods have been used in order to provide a perturbation to the flow; among these are the use of an oscillating trailing edge, acoustically driven slots, external acoustic forcing, and mechanical piston methods. The effect of a planar mechanical piston forcing on a single stream shear layer is presented; it can be noted that this is one of the lesser studied free shear layers. The single stream shear layer can be characterized by its primary flow velocity scale and the thickness of the separating boundary layer. The velocity scale is constant over the length of the flow field; theta (x) can be used as a width scale to characterize the unforced shear layer. In the case of the forced shear layer the velocity field is a function of phase time and definition of a width measure becomes somewhat problematic.

Haw, Richard C.↗
NASA NTRS1985

A Bayesian approach to nonlinear inversion

Powerful methods are now available for solving linear parametric inverse problems. However, many inverse problems which arise in geohysics are nonlinear. Fortunately, it is possible to treat most of these with the air of linear perturbation theory and liner inversion. But a convenient method is needed for assessing the importance of nonlinearity in these quasi-linear problems. The present paper provides such a method. Matsu'ura and Jackson (1984) have presented a simple algorithm for evaluating the asymptotic covariance matrix fo estimation errors. In the present investigation, aspects of linear inversion are discussed, taking into account linear parametric inverse problems, nonuniqueness, prior information, confidence limits, conditional and marginal statistics, the relative importance of the prior and observational data, and standardized variables. Attention is also given to nonlinear inversion, and the application of the considered approaches to a number of examples.

Jackson, D. D.↗
NASA NTRSDate not supplied

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

To understand the behavior of charged lunar regolith when perturbed by lunar landers, it is important to couple the grain dynamics with mechanical and electrical particle interactions. To accomplish this, improvements have been made to a discrete element method (DEM) software package to include both short- and long-range interactions between spherical particles. Short-range interactions rely on contact between the particles, such as electrical conduction and triboelectric charge transfer driven by work functions. Long-range interactions act at a distance between every pairing of particles, such as electrostatic forces and gravitational forces. Results from simulations between a few particles are compared with theory to verify these added behaviors prior to scaling up to more complex scenarios. The radii, initial charges, electrical conductivities, work functions, and separation of the particles are varied and the resultant charges as well as the time required to reach the final state are determined.

Electrostatics↗
NASA NTRS1985

A comparison of finite-difference and finite-element methods for calculating free edge stresses in composites

It is pointed out that edge delamination is a serious failure mechanism for laminated composite materials. Various numerical methods have been utilized in attempts to calculate the interlaminar stress components which precede delamination in a laminate. There are, however, discrepancies regarding the results provided by different methods, taking into account a finite-difference procedure, a perturbation procedure, and finite element approaches. The present investigation has the objective to assess the capacity of a finite difference method to predict the character and magnitude of the interlaminar stress distributions near an interface corner. A second purpose of the investigation is to determine if predictions by finite element method in-plane, interlaminar stress components near an interface corner represent actual laminate behavior.

Bauld, N. R., Jr.↗
NASA NTRS1990

Numerical simulation of spatially-evolving instability

A computational study of the spatial stability of plane Poiseuille flow is presented. The numerical scheme employs a time-splitting method to integrate the full Navier-Stokes equations using spectral collocation/finite-difference discretization on a non-staggered mesh. The eigenvalue decomposition procedure is applied for the solution of the Poisson equations using the capacitance matrix technique. The buffer domain method is incorporated for the outflow boundary conditions. The input perturbation velocities are obtained by solving the Orr-Sommerfeld equation for the nonlinear eigenvalue problem employing the companion matrix method. Computational results are compared with the linear theory for two-dimensional disturbances.

Danabasoglu, G.↗
NASA NTRS1994

Application of an improved Nelson-Nguyen analysis to eccentric, arbitrary profile liquid annular seals

An improved dynamic analysis for liquid annular seals with arbitrary profile based on a method, first proposed by Nelson and Nguyen, is presented. An improved first order solution that incorporates a continuous interpolation of perturbed quantities in the circumferential direction, is presented. The original method uses an approximation scheme for circumferential gradients, based on Fast Fourier Transforms (FFT). A simpler scheme based on cubic splines is found to be computationally more efficient with better convergence at higher eccentricities. A new approach of computing dynamic coefficients based on external specified load is introduced. This improved analysis is extended to account for arbitrarily varying seal profile in both axial and circumferential directions. An example case of an elliptical seal with varying degrees of axial curvature is analyzed. A case study based on actual operating clearances of an interstage seal of the Space Shuttle Main Engine High Pressure Oxygen Turbopump is presented.

Padavala, Satyasrinivas↗
NASA NTRS1994

Evaluating the design of an earth radiation budget instrument with system simulations. Part 2: Minimization of instantaneous sampling errors for CERES-I

Much of the new record of broadband earth radiation budget satellite measurements to be obtained during the late 1990s and early twenty-first century will come from the dual-radiometer Clouds and Earth's Radiant Energy System Instrument (CERES-I) flown aboard sun-synchronous polar orbiters. Simulation studies conducted in this work for an early afternoon satellite orbit indicate that spatial root-mean-square (rms) sampling errors of instantaneous CERES-I shortwave flux estimates will range from about 8.5 to 14.0 W/m on a 2.5 deg latitude and longitude grid resolution. Rms errors in longwave flux estimates are only about 20% as large and range from 1.5 to 3.5 W/sq m. These results are based on an optimal cross-track scanner design that includes 50% footprint overlap to eliminate gaps in the top-of-the-atmosphere coverage, and a 'smallest' footprint size to increase the ratio in the number of observations lying within to the number of observations lying on grid area boundaries. Total instantaneous measurement error also depends on the variability of anisotropic reflectance and emission patterns and on retrieval methods used to generate target area fluxes. Three retrieval procedures from both CERES-I scanners (cross-track and rotating azimuth plane) are used. (1) The baseline Earth Radiaton Budget Experiment (ERBE) procedure, which assumes that errors due to the use of mean angular dependence models (ADMs) in the radiance-to-flux inversion process nearly cancel when averaged over grid areas. (2) To estimate N, instantaneous ADMs are estimated from the multiangular, collocated observations of the two scanners. These observed models replace the mean models in computation of satellite flux estimates. (3) The scene flux approach, conducts separate target-area retrievals for each ERBE scene category and combines their results using area weighting by scene type. The ERBE retrieval performs best when the simulated radiance field departs from the ERBE mean models by less than 10%. For larger perturbations, both the scene flux and collocation methods produce less error than the ERBE retrieval. The scene flux technique is preferable, however, because it involves fewer restrictive assumptions.

Stowe, Larry↗
NASA NTRS1994

An initial investigation into methods of computing transonic aerodynamic sensitivity coefficients

The primary accomplishments of the project are as follows: (1) Using the transonic small perturbation equation as a flowfield model, the project demonstrated that the quasi-analytical method could be used to obtain aerodynamic sensitivity coefficients for airfoils at subsonic, transonic, and supersonic conditions for design variables such as Mach number, airfoil thickness, maximum camber, angle of attack, and location of maximum camber. It was established that the quasi-analytical approach was an accurate method for obtaining aerodynamic sensitivity derivatives for airfoils at transonic conditions and usually more efficient than the finite difference approach. (2) The usage of symbolic manipulation software to determine the appropriate expressions and computer coding associated with the quasi-analytical method for sensitivity derivatives was investigated. Using the three dimensional fully conservative full potential flowfield model, it was determined that symbolic manipulation along with a chain rule approach was extremely useful in developing a combined flowfield and quasi-analytical sensitivity derivative code capable of considering a large number of realistic design variables. (3) Using the three dimensional fully conservative full potential flowfield model, the quasi-analytical method was applied to swept wings (i.e. three dimensional) at transonic flow conditions. (4) The incremental iterative technique has been applied to the three dimensional transonic nonlinear small perturbation flowfield formulation, an equivalent plate deflection model, and the associated aerodynamic and structural discipline sensitivity equations; and coupled aeroelastic results for an aspect ratio three wing in transonic flow have been obtained.

Carlson, Leland A.↗
NASA NTRS1991

Experimental verification of modified command shaping using a flexible manipulator

A brief discussion introduces the original input shaping method applied to a system with varying parameters. A change in parameters causes a vibration in the system and a modified command shaping technique is created to eliminate this unwanted motion. Using a two degree of freedom (DOF) flexible manipulator, experiments are conducted by varying perturbations to input trajectories to compare the modified method to a previously used adaptive proportional plus derivative (P.D.) control method. The control scheme that produces the smaller magnitude of resonant vibration at the first natural frequency of the robot is considered the more effective control method.

Magee, David P.↗
NASA NTRS1988

Intermolecular potential for thermal H2O-He collisions

Theoretical potentials for rotational excitation of H2O by He were constructed via several methods, all of which start with a large basis set SCF interaction. The semiempirical Hartree-Fock with damped dispersion model adds a damped long-range attraction with parameters adjusted to fit experimental total differential cross sections. Purely ab initio potentials add correlation energies obtained via perturbation theory (MP2 and MP4) or a variational method (ICF1). Scattering calculations were performed on all surfaces to compare wih available beam scattering and pressure broadening data and to assess sensitivity of state-to-state rates to uncertainties in the potential. From comparison with the limited experimental data, the ICF1 surface appears to be marginally better than the MP4 surface. Thermal rates calculated from this surface should be accurate to better than 50 percent, at least for the larger, more important rates.

Palma, Amedeo↗
NASA NTRS1996

Summary of work on shock wave feature extraction in 3-D datasets

A method for extracting and visualizing shock waves from three dimensional data-sets is discussed. Issues concerning computation time, robustness to numerical perturbations, and noise introduction are considered and compared with other methods. Finally, results using this method are discussed.

Hesselink, Lambertus↗
NASA NTRS1974

Techniques for the analysis of geodynamic effects using laser data

A technique is described by which laser measurements may be used in an explicit definition of geodetic parameters at the one-meter level of resolution. Observations made by a tracking station at the Goddard Space Flight Center of the Beacon Explorer C satellite are analyzed with this technique to yield highly precise measures of perturbations in the satellite's inclination, including the effect of earth and ocean tides and variations in the station's latitude due to polar motion. The tracking configuration is described, and the basic technique is outlined (analysis of six-hour observations with a weighted least squares orbit determination method). Techniques for analyzing orbital errors and inclination perturbations are described, and it is shown that quarter-day spans of laser data can be employed to monitor the inclination of the satellite to the order of 0.01 arcseconds precision over a period of 17 months.

Dunn, P. J.↗
NASA NTRS1990

Differential equation based method for accurate approximations in optimization

A method to efficiently and accurately approximate the effect of design changes on structural response is described. The key to this method is to interpret sensitivity equations as differential equations that may be solved explicitly for closed form approximations, hence, the method is denoted the Differential Equation Based (DEB) method. Approximations were developed for vibration frequencies, mode shapes and static displacements. The DEB approximation method was applied to a cantilever beam and results compared with the commonly-used linear Taylor series approximations and exact solutions. The test calculations involved perturbing the height, width, cross-sectional area, tip mass, and bending inertia of the beam. The DEB method proved to be very accurate, and in most cases, was more accurate than the linear Taylor series approximation. The method is applicable to simultaneous perturbation of several design variables. Also, the approximations may be used to calculate other system response quantities. For example, the approximations for displacements are used to approximate bending stresses.

Pritchard, Jocelyn I.↗