Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite amplitude method”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Three-dimensional elastic-plastic finite-element analysis of fatigue crack propagation

Fatigue cracks are a major problem in designing structures subjected to cyclic loading. Cracks frequently occur in structures such as aircraft and spacecraft. The inspection intervals of many aircraft structures are based on crack-propagation lives. Therefore, improved prediction of propagation lives under flight-load conditions (variable-amplitude loading) are needed to provide more realistic design criteria for these structures. The main thrust was to develop a three-dimensional, nonlinear, elastic-plastic, finite element program capable of extending a crack and changing boundary conditions for the model under consideration. The finite-element model is composed of 8-noded (linear-strain) isoparametric elements. In the analysis, the material is assumed to be elastic-perfectly plastic. The cycle stress-strain curve for the material is shown Zienkiewicz's initial-stress method, von Mises's yield criterion, and Drucker's normality condition under small-strain assumptions are used to account for plasticity. The three-dimensional analysis is capable of extending the crack and changing boundary conditions under cyclic loading.

Goglia, G. L.↗
NASA NTRS1986

Wave-interactions in three-dimensional boundary layers

A Fourier-Chebyshev spectral method is used to study the nonlinear stability and disturbance wave-interactions in a two-dimensional boundary layer on a curved wall and a three-dimensional boundary layer on a rotating disk. Goertler and crossflow instability, respectively, are the dominant mechanisms in these boundary layers. For the range of parameters considered, the present computational model does not support the Goertler/TS interaction of the type proposed by Nayfeh (1981). According to his model, oblique TS waves with spanwise wavelengths twice that of the finite-amplitude Goertler vortex are excited in the boundary layer. The present nonlinear calculations show that both the Goertler and crosssflow vortices develop strong instabilities in their harmonics with half the wavelength of the primary vortices.

Malik, M. R.↗
NASA NTRS1987

Numerical modeling of the vortex/airfoil interaction

A modeling of the vortex-airfoil interaction is presented in which the finite-area of the real vortices is taken into consideration. Two vortex models are used. In the first, a disturbed piece of vorticity layer is simulated by four rows of discrete vortices of small strength. In the second, a number of discrete vortices is arranged within a circle. The first model may simulate a shear layer or a wake, while the second, a well-formed vortex. The method was applied to the calculation of the pressure induced on the surface of the airfoil by the interacting vortex. Both models give similar results. It was found that for large distances of the vortex from the surface of the airfoil, the consideration or not of the finite-area of the vortex is not a significant factor in determining the induced pressure field. However, when the distance of the vortex from the surface is reduced, its shape is distorted and the induced pressure pulses have lower amplitude than the ones induced by an equivalent point vortex. In the limit, where the vortex impinges on the leading edge of the airfoil, it is split into two and the time dependent pressure coefficient takes even negative values at some time intervals.

Panaras, Argyris G.↗
NASA NTRS2010

Comparison of Methods to Predict Lower Bound Buckling Loads of Cylinders Under Axial Compression

Results from a numerical study of the buckling response of two different orthogrid stiffened circular cylindrical shells with initial imperfections and subjected to axial compression are used to compare three different lower bound buckling load prediction techniques. These lower bound prediction techniques assume different imperfection types and include an imperfection based on a mode shape from an eigenvalue analysis, an imperfection caused by a lateral perturbation load, and an imperfection in the shape of a single stress-free dimple. The STAGS finite element code is used for the analyses. Responses of the cylinders for ranges of imperfection amplitudes are considered, and the effect of each imperfection is compared to the response of a geometrically perfect cylinder. Similar behavior was observed for shells that include a lateral perturbation load and a single dimple imperfection, and the results indicate that the predicted lower bounds are much less conservative than the corresponding results for the cylinders with the mode shape imperfection considered herein. In addition, the lateral perturbation technique and the single dimple imperfection produce response characteristics that are physically meaningful and can be validated via testing.

Haynie, Waddy T.↗
NASA NTRS2002

Computational Study of Near-limit Propagation of Detonation in Hydrogen-air Mixtures

A computational investigation of the near-limit propagation of detonation in lean and rich hydrogen-air mixtures is presented. The calculations were carried out over an equivalence ratio range of 0.4 to 5.0, pressures ranging from 0.2 bar to 1.0 bar and ambient initial temperature. The computations involved solution of the one-dimensional Euler equations with detailed finite-rate chemistry. The numerical method is based on a second-order spatially accurate total-variation-diminishing (TVD) scheme, and a point implicit, first-order-accurate, time marching algorithm. The hydrogen-air combustion was modeled with a 9-species, 19-step reaction mechanism. A multi-level, dynamically adaptive grid was utilized in order to resolve the structure of the detonation. The results of the computations indicate that when hydrogen concentrations are reduced below certain levels, the detonation wave switches from a high-frequency, low amplitude oscillation mode to a low frequency mode exhibiting large fluctuations in the detonation wave speed; that is, a 'galloping' propagation mode is established.

Yungster, S.↗
NASA NTRS2004

Numerical Simulations of Buoyancy Effects in low Density Gas Jets

This paper deals with the computational analysis of buoyancy effects in the near field of an isothermal helium jet injected into quiescent ambient air environment. The transport equations of helium mass fraction coupled with the conservation equations of mixture mass and momentum were solved using a staggered grid finite volume method. Laminar, axisymmetric, unsteady flow conditions were considered for the analysis. An orthogonal system with non-uniform grids was used to capture the instability phenomena. Computations were performed for Earth gravity and during transition from Earth to different gravitational levels. The flow physics was described by simultaneous visualizations of velocity and concentration fields at Earth and microgravity conditions. Computed results were validated by comparing with experimental data substantiating that buoyancy induced global flow oscillations present in Earth gravity are absent in microgravity. The dependence of oscillation frequency and amplitude on gravitational forcing was presented to further quantify the buoyancy effects.

Satti, R. P.↗
NASA NTRS1990

Numerical investigation of laminar-turbulent transition in a flat plate wake

Lamina-turbulent transition of high-deficit flat plate wakes is investigated by direct numerical simulations using the complete Navier Stokes equations. The simulations are based on a spatial model so that both the base flow and the disturbance flow can develop in the downstream direction. The Navier Stokes equations are used in a vorticity-velocity form and are solved using a combination of finite difference and spectral approximations. Fourier series are used in the spanwise direction. Second-order finite-differences are used to approximate the spatial derivatives in the streamwise and transverse directions. For the temporal discretion, a combination of ADI, Crank-Nicolson, and Adams-Bashforth methods is employed. The discretized velocity equations are solved using fast Helmholtz solvers. Code validation is accomplished by comparison of the numerical results to both linear stability and to experiments. Calculations of two- and/or three-dimensional sinuous and mode disturbances in the wake of flat plate are undertaken. For calculations of two-dimensional disturbances, the wake is forced at an amplitude level so that nonlinear disturbance development may be observed. In addition, the forcing amplitude is varied in order to determine its effect on the disturbance behavior. To investigate the onset of three-dimensionality, the wake is forced with a small-amplitude three-dimensional disturbance and a larger amplitude two-dimensional disturbance. The two-dimensional forcing amplitude is varied in order to determine its influence on the three-dimensional flow field.

Fasel, Hermann F.↗
NASA NTRS1993

Analysis of microstrip dipoles and slots transversely coupled to a microstrip line using the FDTD method

Printed dipole elements and their complement, linear slots, are elementary radiators that have found use in low-profile antenna arrays. Low-profile antenna arrays, in addition to their small size and low weight characteristics, offer the potential advantage of low-cost, high-volume production with easy integration with active integrated circuit components. The design of such arrays requires that the radiation and impedance characteristics of the radiating elements be known. The FDTD (Finite-Difference Time-Domain) method is a general, straight-forward implementation of Maxwell's equations and offers a relatively simple way of analyzing both printed dipole and slot elements. Investigated in this work is the application of the FDTD method to the analysis of printed dipole and slot elements transversely coupled to an infinite transmission line in a multilayered configuration. Such dipole and slot elements may be used in dipole and slot series-fed-type linear arrays, where element offsets and interelement line lengths are used to obtain the desired amplitude distribution and beam direction, respectively. The design of such arrays is achieved using transmission line theory with equivalent circuit models for the radiating elements. In an equivalent circuit model, the dipole represents a shunt impedance to the transmission line, where the impedance is a function of dipole offset, length, and width. Similarly, the slot represents a series impedance to the transmission line. The FDTD method is applied to single dipole and slot elements transversely coupled to an infinite microstrip line using a fixed rectangular grid with Mur's second order absorbing boundary conditions. Frequency-dependent circuit and scattering parameters are obtained by saving desired time-domain quantities and using the Fourier transform. A Gaussian pulse excitation is applied to the microstrip transmission line, where the resulting reflected signal due to the presence of the radiating element is used to determine the equivalent element impedance.

Tulintseff, A. N.↗
NASA NTRS1990

Instabilities and subharmonic resonances of subsonic heated round jets, volume 2

When a jet is perturbed by a periodic excitation of suitable frequency, a large-scale coherent structure develops and grows in amplitude as it propagates downstream. The structure eventually rolls up into vortices at some downstream location. The wavy flow associated with the roll-up of a coherent structure is approximated by a parallel mean flow and a small, spatially periodic, axisymmetric wave whose phase velocity and mode shape are given by classical (primary) stability theory. The periodic wave acts as a parametric excitation in the differential equations governing the secondary instability of a subharmonic disturbance. The (resonant) conditions for which the periodic flow can strongly destabilize a subharmonic disturbance are derived. When the resonant conditions are met, the periodic wave plays a catalytic role to enhance the growth rate of the subharmonic. The stability characteristics of the subharmonic disturbance, as a function of jet Mach number, jet heating, mode number and the amplitude of the periodic wave, are studied via a secondary instability analysis using two independent but complementary methods: (1) method of multiple scales, and (2) normal mode analysis. It is found that the growth rates of the subharmonic waves with azimuthal numbers beta = 0 and beta = 1 are enhanced strongly, but comparably, when the amplitude of the periodic wave is increased. Furthermore, compressibility at subsonic Mach numbers has a moderate stabilizing influence on the subharmonic instability modes. Heating suppresses moderately the subharmonic growth rate of an axisymmetric mode, and it reduces more significantly the corresponding growth rate for the first spinning mode. Calculations also indicate that while the presence of a finite-amplitude periodic wave enhances the growth rates of subharmonic instability modes, it minimally distorts the mode shapes of the subharmonic waves.

Ng, Lian Lai↗
NASA NTRS1988

A new technique for measuring extragalactic distances

A relatively direct technique of determining extragalactic distances is described. The method relies on measuring the luminosity fluctuations that arise from the counting statistics of the stars contributing the flux in each pixel of a high signal-to-noise CCD image of a galaxy. The amplitude of these fluctuations is inversely proportional to the distance of the galaxy. This approach bypasses most of the successive stages of calibration required in the traditional extragalactic distance ladder; the only serious drawback to this method is that it requires an accurate knowledge of the bright end of the luminosity function. Potentially, this method can produce accurate distances of elliptical galaxies and spiral bulges at distances out to about 20 Mpc. In this paper, it is explained how to calculate the value of the fluctuations, taking into account various sources of contamination and the effects of finite spatial resolution, and the feasibility and limitations of this technique are demonstrated via simulations and CCD images of M 32 and N 3379.

Tonry, John↗
NASA NTRS1997

Computational Considerations for the Simulation of Discontinuous Flows

The numerical study of aeroacoustic problems places stringent demands on the choice of a computational algorithm, because it requires the ability to propagate disturbances of small amplitude and short wavelength. The demands are particularly high when shock waves are involved, because the chosen algorithm must also resolve discontinuities in the solution. In a previous work the capabilities and deficiencies of shock-capturing methods for aeroacoustic problems were demonstrated using a high-order essentially nonoscillatory (ENO) numerical method. It was shown that first-order results are obtained when simulating time-dependent flows with discontinuities. The present study reaffirms this conclusion by comparing the ENO results with those obtained using a conventional linear scheme. A sixth-order-accurate compact implicit finite difference scheme is used to investigate various discontinuous flows. The design order of accuracy is achieved in the smooth regions of a steady-state, quasi-one-dimensional Euler test case, as well as in the time-dependent Burgers' equation. However, in the unsteady Euler sound-shock interaction, first-order results are obtained downstream of the shock. A comparison is made between the linear and nonlinear results, noting the advantages of each method. A discontinuous linear model problem is then used to identify the cause of the first-order results. Here, the nature of the solution error is quantified as being predominantly a numerical phase shift, and a post-processing procedure is demonstrated which increases the solution accuracy downstream of the discontinuity to second-order.

Carpenter Mark H.↗
NASA NTRS1997

Ultrasonic Nondestructive Characterization of Adhesive Bonds

Qualitative measurements of adhesion or binding forces can be accomplished, for example, by using the reflection coefficient of an ultrasound or by using thermal waves (Light and Kwun, 1989, Achenbach and Parikh, 1991, and Bostrom and wickham, 1991). However, a quantitative determination of binding forces is rather difficult. It has been observed that higher harmonics of the fundamental frequency are generated when an ultrasound passes through a nonlinear material. It seems that such non-linearity can be effectively used to characterize the bond strength. Several theories have been developed to model this nonlinear effect (Adler and Nagy, 1991; Achenbach and Parikh, 1991; Parikh and Achenbach, 1992; and Hirose and Kitahara, 1992; Anastasi and Roberts, 1992). Based on a microscopic description of the nonlinear interface binding force, a quantitative method was presented by Pangraz and Arnold (1994). Recently, Tang, Cheng and Achenbach (1997) made a comparison between the experimental and simulated results based on this theoretical model. A water immersion mode-converted shear wave through-transmission setup was used by Berndt and Green (1997) to analyze the nonlinear acoustic behavior of the adhesive bond. In this project, the nonlinear responses of an adhesive joint was investigated through transmission tests of ultrasonic wave and analyzed by the finite element simulations. The higher order harmonics were obtained in the tests. It is found that the amplitude of higher harmonics increases as the aging increases, especially the 3dorder harmonics. Results from the numerical simulation show that the material nonlinearity does indeed generate higher order harmonics. In particular, the elastic-perfect plastic behavior generates significant 3rd and 5th order harmonics.

Qu, Jianmin↗
NASA NTRS1991

Use of finite volume schemes for transition simulation

The use of finite-volume methods in the study of spatially and temporally evolving transitional flows over a flat plate is investigated. Schemes are developed with both central and upwind differencing. The compressible Navier-Stokes equations are solved with a Runge-Kutta time stepping scheme. Disturbances are determined using linear theory and superimposed at the inflow boundary. Time accurate integration is then used to allow temporal and spatial disturbance evolution. Characteristic-based boundary conditions are employed. The requirements of using finite-volume algorithms are studied in detail. Special emphasis is placed on difference schemes, grid resolution, and disturbance amplitudes. Moreover, comparisons are made with linear theory for small amplitude disturbances. Both subsonic and supersonic flows are considered, and it is shown that the locations of branch 1 and branch 2 of the neutral stability curve are well predicted, given sufficient resolution.

Fenno, Charles C., Jr.↗
NASA NTRS1991

Waveform effects of a metastable olivine tongue in subducting slabs

Velocity models of subducting slabs with a kinetically-depressed olivine to beta- and gamma-spinel transition are constructed, and the effect that such structures would have on teleseismic P waveforms are examined using a full-wave finite-difference method. These 2D calculations yielded waveforms at a range of distances in the downdip direction. The slab models included a wedge-shaped, low-velocity metastable olivine tongue (MOTO) to a depth of 670 km, as well as a plausible thermal anomaly; one model further included a 10-km-thick fast layer on the surface of the slab. The principal effect of MOTO is to produce grazing reflections at wide angles off the phase boundary, generating a secondary arrival 0 to 4 seconds after the initial arrival depending on the take-off angle. The amplitude and timing of this feature vary with the lateral location of the seismic source within the slab cross-section.

Vidale, John E.↗
NASA NTRS1983

Unsteady laminar boundary-layer calculations on oscillating configurations including backflow. Part 1: Flat plate, oscillating in its own plane

A finite difference method has been developed to calculate the unsteady boundary layer over an oscillating flat plate. Low- and high frequency approximations were used for comparison with numerical results. Special emphasis was placed on the behavior of the flow and on the numerical calculation procedure as soon as reversed flow has occurred over part of the oscillation cycle. The numerical method displayed neither problems nor singular behavior at the beginning of or within the reversed flow region. Calculations, however, came to a limit where the back-flow region reached the plate's leading edge in the case of high oscillation amplitudes. It is assumed that this limit is caused by the special behavior of the flow at the plate's leading edge where the boundary layer equations are not valid.

Geissler, W.↗
NASA NTRS1992

Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

Joslin, Ronald D.↗
NASA NTRS1986

Dynamic analysis and experiment methods for a generic space station model

Modal vibration tests in conjunction with finite element analysis were used to characterize a generic dynamic model. The model consists of five substructures to simulate the multi-body, low frequency nature of large space structures. Static tests were used to refine the substructure analytical models prior to full assemblage analysis. The effects of a cable suspension system are analyzed using prestressed vibration analysis. Coupling between a cable suspension mode and model bending mode is shown to be influenced by the distance from the model center of gravity to the cable-to-model attachment location. A damping characterization method using noncontacting exciters was used to measure amplitude dependent damping. Frequency and damping measurements in ambient air and at near-vacuum conditions are presented.

Belvin, W. K.↗
NASA NTRS1987

Prediction of high-speed rotor noise with a Kirchhoff formula

A new methodology has been developed to predict the impulsive noise generated by a transonic rotor blade. The formulation uses a full-potential finite-difference method to obtain the pressure field close to the blade. A Kirchhoff integral formulation is then used to extend these finite-difference results into the far-field. This Kirchhoff formula is written in a blade-fixed coordinate system. It requires initial data across a plane at the sonic radius. This data is provided by the finite-difference solution. Acoustic pressure predictions show excellent agreement with hover experimental data for two hover cases of 0.88 and 0.90 tip Mach number, the latter of which has delocalized transonic flow. These results represent the first successful prediction technique for peak pressure amplitudes using a computational code.

Purcell, Timothy W.↗