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 109 records · Page 6

Free-Surface and Contact Line Motion of Liquid in Microgravity

This project involves fundamental studies of the role of nonlinearity in determining the motion of liquid masses under the principal influences of surface tension, viscosity and inertia. Issues to be explored are relevant to aspects of terrestrial processes, as well as being immediately applicable to fluid management in a low-gravity environment. Specific issues include: (1) the mechanic's of liquid masses in large-amplitude motions, (2) the influence of bounding surfaces on the motion, and (3) the ability of such surfaces to control liquid motion by wetting forces, especially when they are augmented by various surface treatments. Mathematical techniques include asymptotic analysis of the governing equations, for problem simplification, and numerical simulation, using both boundary-element and finite-difference methods. The flow problem is divided into an 'outer' or inviscid potential-flow region and one or more inner, or viscous dominated, regions. Relevant to one inner region, the vicinity of the contact line, we discuss time-dependent simulation of slow droplet motion, on a surface of variable wettability, using the lubrication approximation. The simulation uses a disjoining pressure model and reproduces realistic wetting-dewetting behavior.

Schwartz, Leonard W.↗

Wave propagation in a random medium

A simple technique is used to derive statistical characterizations of the perturbations imposed upon a wave (plane, spherical or beamed) propagating through a random medium. The method is essentially physical rather than mathematical, and is probably equivalent to the Rytov method. The limitations of the method are discussed in some detail; in general they are restrictive only for optical paths longer than a few hundred meters, and for paths at the lower microwave frequencies. Situations treated include arbitrary path geometries, finite transmitting and receiving apertures, and anisotropic media. Results include, in addition to the usual statistical quantities, time-lagged functions, mixed functions involving amplitude and phase fluctuations, angle-of-arrival covariances, frequency covariances, and other higher-order quantities.

Lee, R. W.↗

Combustion instability with partial length acoustic liners

The development of analytical techniques and numerical methods for the prediction of the stability behavior of liquid propellant rocket combustors is discussed. The studies involved nonlinear wave propagation effects and nonlinear combustion zone, liner, and nozzle responses. A second primary extension of the basic theory was directed at the inclusion of the effects of injector face baffles. Attention was restricted to the linear (small amplitude) problem, finite Mach numbers, three dimensional oscillations, and distributed combustion effects.

Mitchell, C. E.↗

Three-dimensional modelling in magnetotelluric and magnetic variational sounding

The Galerkin finite-element method is used to obtain approximate solutions for the three-dimensional induction problem. A rectangular conductive prism is considered as an example, and solutions are obtained for linear and circularly polarized incident plane-wave fields. Magnetotelluric tensor impedances and magnetic transfer functions are computed. Polar diagrams of the tensor impedances and magnetic transfer functions along with their amplitude contour maps are presented. The dimensionality parameter, skew, is contoured at the surface of the earth. It is shown that the relative amplitudes and shapes of the additional and principal impedance polar diagrams can be used to determine the dimensionality of geoelectrical structures. Stations with skew values greater than 0.2 are significantly influenced by the three-dimensionality of the geoelectric structure. The amplitudes of the magnetic transfer function and the orientations of its polar diagrams exhibit large anomalies in the vicinity of the intersection of the lateral contacts.

Reddy, I. K.↗

Computation of two-dimensional turbulent flow at subsonic Mach numbers over thick trailing edges

An implicit time marching finite difference method is used to predict two dimensional turbulent flow at a Reynolds number of 440,000 and a Mach number of 0.574 over a shortened NACA 0012 airfoil with a trailing edge of 4.5% thickness and semicircular shape. The flow is found to be unsteady but periodic in the trailing edge region. Thus, lift and drag fluctuate at small amplitudes around mean values and at distinct frequencies.

Drescher, P.↗

Flutter Analysis of the Transonic Truss-Braced Wing Aircraft Using Transonic Correction

This paper describes a flutter analysis method for the Transonic Truss-Braced Wing aircraft using a vortex lattice method coupled to an unsteady transonic correction method to account for unsteady aerodynamics in transonic flow. A steady-state vortex-lattice model of the Truss-Braced Wing aircraft is developed using vortex lattice code VSPAERO. A transonic and viscous flow correction method is implemented in the VSPAERO model to account for steady-state transonic and viscous flow effects using transonic small disturbance code TSFOIL coupled to an in-house integral boundary layer code. In addition, a wing-strut interference correction method is developed to account for the transonic interference aerodynamics in the strut juncture region using high-fidelity CFD code FUN3D. A structural dynamic finite-element model of the Truss-Braced Wing aircraft is developed using BEAM3D in-house finite-element code and is coupled to the VSPAERO. The BEAM3D model includes a geometric nonlinearity due to the tension in the strut which causes a deflection-dependent nonlinear stiffness. An unsteady transonic correction method is developed to better capture the unsteady aerodynamics in transonic flow. The unsteady transonic correction method makes use of the Theodorsen’s theory to account for the amplitude and phase shift of the unsteady lift coefficient in transonic flow. A preliminary flutter analysis of the Truss-Braced Wing aircraft is conduct to illustrate the unsteady transonic correction approach.

Nhan Nguyen↗

Nonlinear vibrations of thin arbitrarily laminated composite plates subjected to harmonic excitations using DKT elements

A finite element formulation is presented for determining the large-amplitude free and steady-state forced vibration response of arbitrarily laminated anisotropic composite thin plates using the Discrete Kirchhoff Theory (DKT) triangular elements. The nonlinear stiffness and harmonic force matrices of an arbitrarily laminated composite triangular 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, triangular plates, different boundary conditions, lamination angles, number of plies, and uniform versus concentrated loads are presented.

Chiang, C. K.↗

Prediction of finite pressure oscillations in stable rocket motors

A methodology is outlined for predicting the amplitude of forced vibrations in the acoustic cavity of a solid rocket motor. The equation for forced vibration of the motor cavity acoustic system is written by parallel with the acoustic mechanical analogy. Acoustic and aeroacoustic theory are used to predict the frequency and intensity of vortex systems or turbulence created by passage of the mean flow over geometric discontinuities in the motor port. Approximate methods are presented for coupling the acoustic field due to the flow noise with the chamber acoustics and for summing the effect upon the multiple acoustic modes.

Hessler, R. O.↗

Acoustic Treatment Design Scaling Methods: Numerical Simulation of the Nonlinear Acoustic Impedance of a Perforated Plate Single-Degree-of-Freedom Resonator Using a Time-Domain Finite Difference Method - Volume 4

Single-degree-of-freedom resonators consisting of honeycomb cells covered by perforated facesheets are widely used as acoustic noise suppression liners in aircraft engine ducts. The acoustic resistance and mass reactance of such liners are known to vary with the intensity of the sound incident upon the panel. Since the pressure drop across a perforated liner facesheet increases quadratically with the flow velocity through the facesheet, this is known as the nonlinear resistance effect. In the past, two different empirical frequency domain models have been used to predict the Sound Pressure Level effect of the incident wave on the perforated liner impedance, one that uses the incident particle velocity in isolated narrowbands, and one that models the particle velocity as the overall velocity. In the absence of grazing flow, neither frequency domain model is entirely accurate in predicting the nonlinear effect that is measured for typical perforated sheets. The time domain model is developed in an attempt to understand and improve the model for the effect of spectral shape and amplitude of multi-frequency incident sound pressure on the liner impedance. A computer code for the time-domain finite difference model is developed and predictions using the models are compared to current frequency-domain models.

Kraft, R. E.↗

Finite-amplitude thermal convection in a spherical shell

The properties of finite-amplitude thermal convection for a Boussinesq fluid contained in a spherical shell are investigated. All nonlinear terms are retained in the equations, and both axisymmetric and nonaxisymmetric solutions are studied. The velocity is expanded in terms of poloidal and toroidal vectors. Spherical surface harmonics resolve the horizontal structure of the flow, but finite differences are used in the vertical. With a few modifications, the transform method developed by Orszag (1970) is used to calculate the nonlinear terms, while Green's function techniques are applied to the poloidal equation and diffusion terms.

Young, R. E.↗

Small Crack Growth and Fatigue Life Predictions for High-Strength Aluminium Alloys: Experimental and Fracture Mechanics Analysis - Part 1

The small crack effect was investigated in two high-strength aluminium alloys: 7075-T6 bare and LC9cs clad alloy. Both experimental and analytical investigations were conducted to study crack initiation and growth of small cracks. In the experimental program, fatigue tests, small crack and large crack tests A,ere conducted under constant amplitude and Mini-TWIST spectrum loading conditions. A pronounced small crack effect was observed in both materials, especially for the negative stress ratios. For all loading conditions, most of the fatigue life of the SENT specimens was shown to be crack propagation from initial material defects or from the cladding layer. In the analysis program, three-dimensional finite element and A weight function methods were used to determine stress intensity factors and to develop SIF equations for surface and corner cracks at the notch in the SENT specimens. A plastisity-induced crack-closure model was used to correlate small and large crack data, and to make fatigue life predictions, Predicted crack-growth rates and fatigue lives agreed well with experiments. A total fatigue life prediction method for the aluminum alloys was developed and demonstrated using the crack-closure model.

Wu, X. R.↗

An introduction to wavelet analysis in oceanography and meteorology - With application to the dispersion of Yanai waves

Wavelet analysis is a relatively new technique that is an important addition to standard signal analysis methods. Unlike Fourier analysis that yields an average amplitude and phase for each harmonic in a dataset, the wavelet transform produces an instantaneous estimate or local value for the amplitude and phase of each harmonic. This allows detailed study of nonstationary spatial or time-dependent signal characteristics. The wavelet transform is discussed, examples are given, and some methods for preprocessing data for wavelet analysis are compared. By studying the dispersion of Yanai waves in a reduced gravity equatorial model, the usefulness of the transform is demonstrated. The group velocity is measured directly over a finite range of wavenumbers by examining the time evolution of the transform. The results agree well with linear theory at higher wavenumber but the measured group velocity is reduced at lower wavenumbers, possibly due to interaction with the basin boundaries.

Meyers, Steven D.↗

Advanced stability analysis for laminar flow control

Five classes of problems are addressed: (1) the extension of the SALLY stability analysis code to the full eighth order compressible stability equations for three dimensional boundary layer; (2) a comparison of methods for prediction of transition using SALLY for incompressible flows; (3) a study of instability and transition in rotating disk flows in which the effects of Coriolis forces and streamline curvature are included; (4) a new linear three dimensional instability mechanism that predicts Reynolds numbers for transition to turbulence in planar shear flows in good agreement with experiment; and (5) a study of the stability of finite amplitude disturbances in axisymmetric pipe flow showing the stability of this flow to all nonlinear axisymmetric disturbances.

Orszag, S. A.↗

Numerical simulation of laminar-turbulent transition in a spatially-developing flat plate wake

Laminar-turbulent transition of an incompressible flat-plate wake is investigated by direct numerical integration of the Navier-Stokes equations. For the numerical integration, a combination of finite-difference and spectral methods along with an ADI/Crank-Nicolson/Adams-Bashforth time integration scheme is employed. Subject to 2D forcing, the wake exhibited a rapidly-growing fundamental disturbance that quickly saturated. This saturation was due partly to the stabilizing effect of the mean flow distortion. Downstream of the saturation point, disturbance energy was concentrated in the fundamental disturbance, the second harmonic, and the mean flow distortion component. At large amplitude levels, a Karman vortex street formed. Variations in the 2D forcing level did not alter the qualitative behavior of the disturbances. Simulations of 3D breakdown indicates that the presence of large-amplitude, 2D disturbances tends to initially suppress small-amplitude 3D disturbance growth. Following this initial suppression, a resumption of 3D growth is observed that may have been due to a secondary instability mechanism. For high levels of 3D disturbance energy, lambda-shaped vortical structures formed between adjacent Karman vortices.

Dratler, D. I.↗

On eigensolutions for discontinuous liners in a duct containing uniform mean flow

Sound attenuation in a rectangular acoustically lined duct containing uniform mean flow is analytically investigated using the generalized Wiener-Hopf technique. Uniqueness of the solution is enforced for lined sections of the finite axial extent by imposing edge conditions at the liner interface. Possible edge conditions are considered, including the Kutta condition, and the causal solution corresponding to edge conditions is considered the best choice. Solution methods such as the mode matching and singularity methods imply differing edge conditions, and results show that power attenuation is insensitive to the imposed edge conditions, although significant differences are observed for the reflection coefficient. The amplitude of the exponentially increasing instability mode in the lined section must be set to zero as a first approximation to the nonlinear situation, and results indicate that measurements of the reflection factor can be used to make a more definite decision about physically appropriate edge conditions.

Koch, W.↗

An optimized finite-difference scheme for wave propagation problems

Two fully-discrete finite-difference schemes for wave propagation problems are presented, a maximum-order scheme and an optimized (or spectral-like) scheme. Both combine a seven-point spatial operator and an explicit six-stage time-march method. The maximum-order operator is fifth-order in space and is sixth-order in time for a linear problem with periodic boundary conditions. The phase and amplitude errors of the schemes obtained using Fourier analysis are given and compared with a second-order and a fourth-order method. Numerical experiments are presented which demonstrate the usefulness of the schemes for a range of problems. For some problems, the optimized scheme leads to a reduction in global error compared to the maximum-order scheme with no additional computational expense.

Zingg, D. W.↗

A Method for Optimizing Non-Axisymmetric Liners for Multimodal Sound Sources

Central processor unit times and memory requirements for a commonly used solver are compared to that of a state-of-the-art, parallel, sparse solver. The sparse solver is then used in conjunction with three constrained optimization methodologies to assess the relative merits of non-axisymmetric versus axisymmetric liner concepts for improving liner acoustic suppression. This assessment is performed with a multimodal noise source (with equal mode amplitudes and phases) in a finite-length rectangular duct without flow. The sparse solver is found to reduce memory requirements by a factor of five and central processing time by a factor of eleven when compared with the commonly used solver. Results show that the optimum impedance of the uniform liner is dominated by the least attenuated mode, whose attenuation is maximized by the Cremer optimum impedance. An optimized, four-segmented liner with impedance segments in a checkerboard arrangement is found to be inferior to an optimized spanwise segmented liner. This optimized spanwise segmented liner is shown to attenuate substantially more sound than the optimized uniform liner and tends to be more effective at the higher frequencies. The most important result of this study is the discovery that when optimized, a spanwise segmented liner with two segments gives attenuations equal to or substantially greater than an optimized axially segmented liner with the same number of segments.

Watson, W. R.↗

Bow-tie antennas on a dielectric half-space - Theory and experiment

A new formulation is discussed for the rigous calculation of the radiation pattern of a bow-tie antenna of finite length and infinitesimal thickness, placed on a lossless dielectric substrate. The analysis is based on a representation of the current density on the metal surface of the antenna as a sum of an imposed (quasistatic) term and a set of current modes with unknown amplitudes. Free-space fields that are expressed in terms of continuous spectra of symmetrized plane waves are matched to the current modes using the method of moments. The resulting set of equations are solved for the unknown current amplitudes. The calculations show that for increasing bow length the antenna impedance spirals rapidly to a value predicted by transmission line theory. The theory also shows that the E-plane pattern of a two-wavelength, 60-deg bow-tie antenna is dominated by low-loss current modes propagating at the dielectric wavenumber. As the bow tie narrows, the loss of the modes increases, and the dominant wavenumber tends to the quasi-static value. Pattern measurements made at 94 GHz are shown to agree well with theoretical predictions. Measurements for a long-wire antenna, a linear array of bow-tie elements, and a log-periodic antenna are also presented.

Compton, Richard C.↗