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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Micromechanics of crenulated fibers in carbon/carbon composites

The influence of crenulated noncircular fibers on the micromechanical stress states due to a transverse strain and to a temperature change in carbon/carbon composites is examined using the finite element method. Stresses at the interface of both fully bonded and fully disbonded fibers having two crenulation amplitudes and with two fiber volume fractions are presented. In each case, these interface stresses are compared to stresses at the interface of circular fibers which have the same degree of disbond and fiber volume fraction and are under the same loading conditions. For the disbonded cases, deformed meshes showing locations of fiber/matrix contact are also included. In addition to the interface stress states, selected composite properties are also computed and compared in each case examined. Interest in studying noncircular fibers stems from a desire to increase the transverse properties of carbon/carbon by introducing a mechanical interlocking between the fiber and the matrix. Results presented here indicate that this interlocking does in fact occur. Evidence from the interface stress data suggests, however, that any possible advantage of this interlocking may be outweighed by the disadvantage of stress concentrations which arise at the interface due to the crenulated geometry of the fibers.

Carapella, E. E.↗

NASA NTRS1991

A method for the geometrically nonlinear analysis of compressively loaded prismatic composite structures

A method was developed for the geometrically nonlinear analysis of the static response of thin-walled stiffened composite structures loaded in uniaxial or biaxial compression. The method is applicable to arbitrary prismatic configurations composed of linked plate strips, such as stiffened panels and thin-walled columns. The longitudinal ends of the structure are assumed to be simply supported, and geometric shape imperfections can be modeled. The method can predict the nonlinear phenomena of postbuckling strength and imperfection sensitivity which are exhibited by some buckling-dominated structures. The method is computer-based and is semi-analytic in nature, making it computationally economical in comparison to finite element methods. The method uses a perturbation approach based on the use of a series of buckling mode shapes to represent displacement contributions associated with nonlinear response. Displacement contributions which are of second order in the model amplitudes are incorported in addition to the buckling mode shapes. The principle of virtual work is applied using a finite basis of buckling modes, and terms through the third order in the model amplitudes are retained. A set of cubic nonlinear algebraic equations are obtained, from which approximate equilibrium solutions are determined. Buckling mode shapes for the general class of structure are obtained using the VIPASA analysis code within the PASCO stiffened-panel design code. Thus, subject to some additional restrictions in loading and plate anisotropy, structures which can be modeled with respect to buckling behavior by VIPASA can be analyzed with respect to nonlinear response using the new method. Results obtained using the method are compared with both experimental and analytical results in the literature. The configurations investigated include several different unstiffened and blade-stiffening panel configurations, featuring both homogeneous, isotropic materials, and laminated composite material.

Stoll, Frederick↗

NASA NTRS1985

Cellular interface morphologies in directional solidification. III - The effects of heat transfer and solid diffusivity

The shape and stability of two-dimensional finite-amplitude cellular interfaces arising during directional solidification are compared for several solidification models that account differently for latent heat released at the interface, unequal thermal conductivities of melt and solid, and solute diffusivity in the solid. Finite-element analysis and computer-implemented perturbation methods are used to analyze the families of steadily growing cellular forms that evolve from the planar state. In all models a secondary bifurcation between different families of finite-amplitude cells exists that halves the spatial wavelength of the stable interface. The quantitative location of this transition is very dependent on the details of the model. Large amounts of solute diffusion in the solid retard the growth of large-amplitude cells.

Ungar, Lyle H.↗

NASA NTRS1995

Category 5 problem solution using an unstructured finite volume algorithm

For the simulation of flows with complex geometries, unstructured finite volume methods have proven to be very popular, and simulations of a large number of flows have been done with good results using this approach. Since most of the simulations to date were done for steady flows, it is not clear that present unstructured finite volume algorithms can accurately track the unsteady propagation of acoustic waves in a computation. Therefore, there is a need to assess the accuracy of these methods for acoustic calculations. In this paper, we perform the numerical simulation of a very small amplitude acoustic wave incident on the non-uniform steady flow in a quasi- 1 D convergent-divergent nozzle using an unstructured finite volume algorithm with piece-wise linear, least square reconstruction, Roe flux difference splitting, and second-order MacCormack time marching. First, the spatial accuracy of the algorithm is evaluated for the steady flow by running the simulation with a sequence of successively finer meshes. Then the unsteady numerical solution with the acoustic perturbation is presented.

Bui, Trong T.↗

NASA NTRS1991

Special class of nonlinear damping models in flexible space structures

A special class of nonlinear damping models is investigated in which the damping force is proportional to the product of positive integer or the fractional power of the absolute values of displacement and velocity. For a one-degree-of-freedom system, the classical Krylov-Bogoliubov 'averaging' method is used, whereas for a distributed system, both an ad hoc perturbation technique and the finite difference method are employed to study the effects of nonlinear damping. The results are compared with linear viscous damping models. The amplitude decrement of free vibration for a single mode system with nonlinear models depends not only on the damping ratio but also on the initial amplitude, the time to measure the response, the frequency of the system, and the powers of displacement and velocity. For the distributed system, the action of nonlinear damping is found to reduce the energy of the system and to pass energy to lower modes.

Hu, Anren↗

NASA NTRS1990

Numerical simulation of the control of the three-dimensional transition process in boundary layers

Surface heating techniques to control the three-dimensional laminar-turbulent transition process are numerically investigated for a water boundary layer. The Navier-Stokes and energy equations are solved using a fully implicit finite difference/spectral method. The spatially evolving boundary layer is simulated. Results of both passive and active methods of control are shown for small amplitude two-dimensional and three-dimensional disturbance waves. Control is also applied to the early stages of the secondary instability process using passive or active control techniques.

Kral, L. D.↗

NASA NTRS1992

A probabilistic approach to the evaluation of fatigue damage in a space propulsion system injector element

The fatigue damage of a space propulsion system component is computed using probabilistic structural analysis methods. The analysis takes into account the variations in static and dynamic loads, the uncertainty in structural damping, and the scatter in material fatigue resistance. The key elements of the probabilistic approach include: (1) a numerical engine model for describing the global component interface loads consistent with engine balance, (2) models for computing the local mode boundary conditions on the component, (3) a static structural analysis model that captures the strains at a damage critical location as a function of engine performance variables, (4) a finite element model for assessment of the random amplitude stress response due to random base and pressure excitations with uncertain power and correlation length, and (5) advanced first order reliability methods for computing the probabilities associated with the fatigue damage.

Rajagopal, K. R.↗

NASA NTRS1995

Application of essentially nonoscillatory methods to aeroacoustic flow problems

A finite-difference essentially nonoscillatory (ENO) method has been applied to several of the problems prescribed for the workshop sponsored jointly by the Institute for Computer Applications in Science and Engineering and by NASA Langley Research Center entitled 'Benchmark Problems in Computational Aeroacoustics'. The workshop focused on computational challenges specific to aeroacoustics. Among these are long-distance propagation of a short-wavelength disturbance, propagation of small-amplitude disturbances, and nonreflective boundary conditions. The shock capturing-capability inherent to the ENO method effectively eliminates oscillations near shock waves without the need to add and tune dissipation or filter terms. The method-of-lines approach allows the temporal and spatial operators to be chosen separately in accordance with the demands of a particular problem. The ENO method was robust and accurate for all problems in which the propagating wave was resolved with 8 or more points per wavelength. The finite-wave-model boundary condition, a local nonlinear acoustic boundary condition, performed well for the one-dimensional problems. The buffer-domain approach performed well for the two-dimensional test problem. The amplitude of nonphysical reflections were less than 1 percent of the exiting wave's amplitude.

Atkins, Harold L.↗

NASA NTRS1981

Linear/nonlinear behavior in unsteady transonic aerodynamics

The accurate calculation of the aerodynamic forces in unsteady transonic flow requires the solution of the nonlinear flow equations. The aeroelastician, on the other hand, seeks to treat his problems (flutter, for example) by means of linear equations whenever possible. He may do this, even when the underlying flow is nonlinear, if the perturbation forces are linear over some (perhaps small) range of unsteady amplitude of motion. This paper assesses the range of parameters over which linear behavior occurs. In particular calculations are made for an NACA 64A006 airfoil oscillating in pitch over a range of amplitudes, frequencies, and Mach numbers. The primary aerodynamic method used is the well known LTRAN2 code of Ballhaus and Goorjian that provides a finite-difference solution to the low frequency, small disturbance, two-dimensional potential flow equation. Comparisons are made with linear subsonic theory, local linearization, and, for steady flow, with the full potential equation code of Bauer, Garabedian, and Korn.

Dowell, E. H.↗

NASA NTRS1989

Finite element large-amplitude free and forced vibrations of rectangular thin composite plates

A finite element formulation is presented for determining the large-amplitude free and steady-state forced vibration response of arbitrarily laminated anisotropic composite rectangular thin plates. The nonlinear stiffness and harmonic force matrices of an arbitrarily laminated composite rectangular plate element are developed for nonlinear free and forced vibration analyses. The linearized updated-mode method with nonlinear time function approximation is employed for the solution of the system nonlinear eigenvalue equations. The amplitude-frequency relations for convergence with gridwork refinement, different boundary conditions, aspect ratios, lamination angles and number of plies are presented. The finite element results are compared with available approximate continuum solutions.

Chiang, C. K.↗

NASA NTRS1991

Singularities in water waves and Rayleigh-Taylor instability

Singularities in inviscid two-dimensional finite-amplitude water waves and inviscid Rayleigh-Taylor instability are discussed. For the deep water gravity waves of permanent form, through a combination of analytical and numerical methods, results describing the precise form, number, and location of singularities in the unphysical domain as the wave height is increased are presented. It is shown how the information on the singularity in the unphysical region has the same form as for deep water waves. However, associated with such a singularity is a series of image singularities at increasing distances from the physical plane with possibly different behavior. Furthermore, for the Rayleigh-Taylor problem of motion of fluid over a vacuum and for the unsteady water wave problem, integro-differential equations valid in the unphysical region are derived, and how these equations can give information on the nature of singularities for arbitrary initial conditions is shown.

Tanveer, S.↗

NASA NTRS1991

Singularities in water waves and Rayleigh-Taylor instability

Tanveer, S.↗

NASA NTRS2007

Investigation of a Parabolic Iterative Solver for Three-dimensional Configurations

A parabolic iterative solution procedure is investigated that seeks to extend the parabolic approximation used within the internal propagation module of the duct noise propagation and radiation code CDUCT-LaRC. The governing convected Helmholtz equation is split into a set of coupled equations governing propagation in the positive and negative directions. The proposed method utilizes an iterative procedure to solve the coupled equations in an attempt to account for possible reflections from internal bifurcations, impedance discontinuities, and duct terminations. A geometry consistent with the NASA Langley Curved Duct Test Rig is considered and the effects of acoustic treatment and non-anechoic termination are included. Two numerical implementations are studied and preliminary results indicate that improved accuracy in predicted amplitude and phase can be obtained for modes at a cut-off ratio of 1.7. Further predictions for modes at a cut-off ratio of 1.1 show improvement in predicted phase at the expense of increased amplitude error. Possible methods of improvement are suggested based on analytic and numerical analysis. It is hoped that coupling the parabolic iterative approach with less efficient, high fidelity finite element approaches will ultimately provide the capability to perform efficient, higher fidelity acoustic calculations within complex 3-D geometries for impedance eduction and noise propagation and radiation predictions.

Nark, Douglas M.↗

NASA NTRS1989

Advanced methods for 3-D inelastic structural analysis for hot engine structures

Three-dimensional Inelastic Analysis Methods are described. These methods were incorporated into a series of new computer codes embodying a progression of mathematical models (mechanics of materials, specialty finite element, boundary element) for streamlined analysis of hot engine structures such as: (1) combustor liners, (2) turbine blades, and (3) turbine vanes. These models address the effects of high temperatures and thermal/mechanical loadings on the local (stress/strain) and global (displacements, frequencies, amplitudes, buckling) structural behavior of the three respective components. The methods and the three computer codes, referred to as MOMM (Mechanics Of Materials Model), MHOST (MARC-Hot Section Technology), and BEST3D (Boundary Element Stress Technology), have been developed and are briefly described.

Chamis, C. C.↗

NASA NTRS1982

Toward a unified theory of atmospheric convective instability

A nonlinear three-dimensional truncated spectral model of shallow and moist Boussinesq convection indicates that parallel instability and thermal forcing are linked, in view of the fact that only one convective mode exists in which either or both mechanisms are operating to generate convection in the planetary boundary layer. It is also established that the wind field causes two-dimensional roll convection formation, an alignment of the convection with the wind in a preferred manner, and a propagation speed that is related to the wind component perpendicular to the roll axis. Latent heating is responsible for the decrease of the critical value of the environmental lapse rate in accordance with the slice method stability criterion. When only the upper part of the upward branch is moist and all of the downward branch is dry, latent heating also causes a finite-amplitude convective solution for Rayleigh number values lower than the critical value of linear analysis.

Shirer, H. N.↗

NASA NTRS1987

Wave envelope technique for multimode wave guide problems

A numerical solution technique for two-dimensional EM waveguide problems is developed analytically and demonstrated. The frequency-domain waveguide problem is transformed to a time-dependent problem and treated following the limiting-amplitude principle; the wave-envelope technique is applied; and the resulting problem is solved numerically using a finite-difference scheme based on a variation of the finite-area method (Baumeister, 1977 and 1984). Results for several sample problems, obtained in 10 min or less of execution time on PCs, are presented in tables.

Hariharan, S. I.↗

NASA NTRS2000

Generation and Radiation of Acoustic Waves from a 2-D Shear Layer

A parallel numerical simulation of the radiation of sound from an acoustic source inside a 2-D jet is presented in this paper. This basic benchmark problem is used as a test case for scattering problems that are presently being solved by using the Impedance Mismatch Method (IMM). In this technique, a solid body in the domain is represented by setting the acoustic impedance of each medium, encountered by a wave, to a different value. This impedance discrepancy results in reflected and scattered waves with appropriate amplitudes. The great advantage of the use of this method is that no modifications to a simple Cartesian grid need to be made for complicated geometry bodies. Thus, high order finite difference schemes may be applied simply to all parts of the domain. In the IMM, the total perturbation field is split into incident and scattered fields. The incident pressure is assumed to be known and the equivalent sources for the scattered field are associated with the presence of the scattering body (through the impedance mismatch) and the propagation of the incident field through a non-uniform flow. An earlier version of the technique could only handle uniform flow in the vicinity of the source and at the outflow boundary. Scattering problems in non-uniform mean flow are of great practical importance (for example, scattering from a high lift device in a non-uniform mean flow or the effects of a fuselage boundary layer). The solution to this benchmark problem, which has an acoustic wave propagating through a non-uniform mean flow, serves as a test case for the extensions of the IMM technique.

Agarwal, Anurag↗

NASA NTRS1985

Mixed models and reduction method for dynamic analysis of anisotropic shells

A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.

Noor, A. K.↗