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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

A finite difference method for the solution of the transonic flow around harmonically oscillating wings

A finite difference method for the solution of the transonic flow about a harmonically oscillating wing is presented. The partial differential equation for the unsteady transonic flow was linearized by dividing the flow into separate steady and unsteady perturbation velocity potentials and by assuming small amplitudes of harmonic oscillation. The resulting linear differential equation is of mixed type, being elliptic or hyperbolic whereever the steady flow equation is elliptic or hyperbolic. Central differences were used for all derivatives except at supersonic points where backward differencing was used for the streamwise direction. Detailed formulas and procedures are described in sufficient detail for programming on high speed computers. To test the method, the problem of the oscillating flap on a NACA 64A006 airfoil was programmed. The numerical procedure was found to be stable and convergent even in regions of local supersonic flow with shocks.

Ehlers, E. F.↗
NASA NTRS1983

Radial segregation induced by natural convection and melt/solid interface shape in vertical Bridgman growth

The roles of natural convection in the melt and the shape of the melt/solid interface on radial dopant segregation are analyzed for a prototype of vertical Bridgman crystal growth system by finite element methods that solve simultaneously for the velocity field in the melt, the shape of the solidification isotherm, and the temperature distribution in both phases. Results are presented for crystal and melt with thermophysical properties similar to those of gallium-doped germanium in Bridgman configurations with melt below (thermally destabilizing) and above (stabilizing) the crystal. Steady axisymmetric flow are classified according to Rayleigh number as either being nearly the growth velocity, having a weak cellular structure or having large amplitude cellular convention. The flows in the two Bridgman configurations are driven by different temperature gradients and are in opposite directions. Finite element calculations for the transport of a dilute dopant by these flow fields reveal radial segregation levels as large as sixty percent of the mean concentration. Segregation is found most severe at an intermediate value of Rayleigh number above which the dopant distribution along the interface levels as the intensity of the flow increases.

Chang, C. J.↗
NASA NTRS1981

A remark on the theory of measuring thermal diffusivity by the modified Angstrom's method

A theory of the measurement of the thermal diffusivity of a sample by the modified Angstrom method is developed for the case in which radiative heat loss from the end surface of the sample is not negligible, and applied to measurements performed on lunar samples. Formulas allowing sample thermal diffusivity to be determined from the amplitude decay and phase lag of a temperature wave traveling through the sample are derived for a flat disk sample for which only heat loss from the end surface is important, and a sample of finite diameter and length for which heat loss through the end and side surfaces must be considered. It is noted that in the case of a flat disk, measurements at a single angular frequency of the temperature wave are sufficient, while the sample of finite diameter and length requires measurements at two discrete angular frequencies. Comparison of the values of the thermal diffusivities of two lunar samples of dimensions approximately 1 x 1 x 2 cm derived by the present methods and by the Angstrom theory for a finite bar reveals them to differ by not more than 5%, and indicates that more refined data are required as the measurement theory becomes more complicated.

Horai, K.-I.↗
NASA NTRS1979

A simplified computer program for the prediction of the linear stability behavior of liquid propellant combustors

A program for predicting the linear stability of liquid propellant rocket engines is presented. The underlying model assumptions and analytical steps necessary for understanding the program and its input and output are also given. The rocket engine is modeled as a right circular cylinder with an injector with a concentrated combustion zone, a nozzle, finite mean flow, and an acoustic admittance, or the sensitive time lag theory. The resulting partial differential equations are combined into two governing integral equations by the use of the Green's function method. These equations are solved using a successive approximation technique for the small amplitude (linear) case. The computational method used as well as the various user options available are discussed. Finally, a flow diagram, sample input and output for a typical application and a complete program listing for program MODULE are presented.

Mitchell, C. E.↗
NASA NTRS1990

Shear-layer-edge interaction - Simulation by finite-area vortices

The impingement of a row of finite-area vortices on an edge is presently used to efficiently simulate the shear layer/edge interaction, yielding support for the hypothesis that the pressure waves emitted from an impingement edge are generated by the vortices/edge interaction. A parametric application of this method shows that pressure wave amplitude is a function of the length of the succession of vortices and that frequency of their release; this amplitude decreases with decreasing vortex spacing while succession length remains constant, or when succession length decreases while the number of vortices remains constant.

Panaras, Argyris G.↗
NASA NTRS1998

Nonlinear Interaction of Detuned Instability Waves in Boundary-Layer Transition: Amplitude Equations

The non-equilibrium critical-layer analysis of a system of frequency-detuned resonant-triads is presented. In this part of the analysis, the system of partial differential critical-layer equations derived in Part I is solved analytically to yield the amplitude equations which are analyzed using a combination of asymptotic and numerical methods. Numerical solutions of the inviscid non-equilibrium oblique-mode amplitude equations show that the frequency-detuned self-interaction enhances the growth of the lower-frequency oblique modes more than the higher-frequency ones. All amplitudes become singular at the same finite downstream position. The frequency detuning delays the occurrence of the singularity. The spanwise-periodic mean-flow distortion and low-frequency nonlinear modes are generated by the critical-layer interaction between frequency-detuned oblique modes. The nonlinear mean flow and higher harmonics as well as the primary instabilities become as large as the base mean flow in the inviscid wall layer in the downstream region where the distance from the singularity is of the order of the wavelength scale.

Lee, Sang Soo↗
NASA NTRS1978

Non-linear effects on sound propagation in a variable area duct carrying high subsonic Mach number flows

The general one-dimensional non-linear equation for the acoustic velocity potential in a variable area duct carrying high subsonic Mach number flows is presented and solved numerically using an implicit finite difference scheme. For the linearized equation with no flow, the present scheme is compared with the exact solution and the fourth-order Runge-Kutta method with excellent agreement for dimensionless time periods. Non-linear solutions are more sensitive to Mach numbers and exciting amplitudes and less sensitive to exciting frequencies. In general, non-linear effects can be safely neglected for low Mach number flows but must be accounted for when high Mach number flows are encountered.

Lumsdaine, E.↗
NASA NTRS1984

Transition to turbulence in plane channel flows

Results obtained from a numerical simulation of the final stages of transition to turbulence in plane channel flow are described. Three dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time evolution of two and three dimensional finite amplitude disturbances. Computations are performed on CYBER-203 vector processor for a 32x51x32 grid. Results are presented for no-slip boundary conditions at the solid walls as well as for periodic suction blowing to simulate active control of transition by mass transfer. Solutions indicate that the method is capable of simulating the complex character of vorticity dynamics during the various stages of transition and final breakdown. In particular, evidence points to the formation of a lambda-shape vortex and the subsequent system of horseshoe vortices inclined to the main flow direction as the main elements of transition. Calculations involving periodic suction-blowing indicate that interference with a wave of suitable phase and amplitude reduces the disturbance growth rates.

Biringen, S.↗
NASA NTRS1973

The acoustics of aircraft engine-duct systems

Noise generated in aircraft engines is usually suppressed by acoustically treating the engine ducts. The optimization of this treatment requires an understanding of the transmission and attenuation of the acoustic waves. A critical review is presented of the state of the art regarding methods of determining the transmission and attenuation parameters and the effect on these parameters of (1) acoustic properties of liners, (2) the mean velocity, including uniform and shear profiles and nonparallel flow, (3) axial and transverse temperature gradients, (4) slowly and abruptly varying cross sections, and (5) finite-amplitude waves and nonlinear duct liners.

Nayfeh, A. H.↗
NASA NTRS1983

Transition to turbulence in plane channel flow

A numerical simulation of the final stages of transition to turbulence in plane channel flow is reported. Three dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time-evolution of two and three dimensional finite amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32x51x32 grid. Results are presented for no-slip boundary conditions at the solid walls as well as for periodic suction-blowing to simulate active control of transition by mass transfer. Solutions indicate that the method is capable of simulating the complex character of vorticity dynamics during the various stages of transition and final breakdown. In particular, evidence points to the formation of a lambda-shape vortex and the subsequent system of horseshoe vortices inclined to the main flow direction as the main elements of transition. Calculations involving suction-blowing indicate that interference with a wave of suitable phase and amplitude reduces the disturbance growth rates.

Biringen, S.↗
NASA NTRS1987

A linearized Euler analysis of unsteady flows in turbomachinery

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗
NASA NTRS1989

Calculation of unsteady flows in turbomachinery using the linearized Euler equations

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗
NASA NTRS1977

Flutter of articulated pipes at finite amplitude

The plane motion of an articulated pipe made of two segments is examined and the flow velocity at which flutter manifests itself is sought. The pressure in the reservoir feeding the pipe is kept constant. In contrast to previous works, the flow velocity is not taken as a prescribed parameter of the system but is left to follow the laws of motion. This approach requires a nonlinear formulation of the problem and the equations of motion are solved using Krylov-Bogoliubov's method. A graph of the amplitude of the limit cycles, as a function of the fluid-system mass ratio, is presented and conclusions are drawn as to the necessity of considering nonlinearities in the analysis.

Rousselet, J.↗
NASA NTRS1975

Effect of Nozzle Nonlinearities upon Nonlinear Stability of Liquid Propellant Rocket Motors

A three dimensional, nonlinear nozzle admittance relation is developed by solving the wave equation describing finite amplitude oscillatory flow inside the subsonic portion of a choked, slowly convergent axisymmetric nozzle. This nonlinear nozzle admittance relation is then used as a boundary condition in the analysis of nonlinear combustion instability in a cylindrical liquid rocket combustor. In both nozzle and chamber analyses solutions are obtained using the Galerkin method with a series expansion consisting of the first tangential, second tangential, and first radial modes. Using Crocco's time lag model to describe the distributed unsteady combustion process, combustion instability calculations are presented for different values of the following parameters: (1) time lag, (2) interaction index, (3) steady-state Mach number at the nozzle entrance, and (4) chamber length-to-diameter ratio. In each case, limit cycle pressure amplitudes and waveforms are shown for both linear and nonlinear nozzle admittance conditions. These results show that when the amplitudes of the second tangential and first radial modes are considerably smaller than the amplitude of the first tangential mode the inclusion of nozzle nonlinearities has no significant effect on the limiting amplitude and pressure waveforms.

Padmanabhan, M. S.↗
NASA NTRS1994

Radiation boundary condition and anisotropy correction for finite difference solutions of the Helmholtz equation

In this paper finite-difference solutions of the Helmholtz equation in an open domain are considered. By using a second-order central difference scheme and the Bayliss-Turkel radiation boundary condition, reasonably accurate solutions can be obtained when the number of grid points per acoustic wavelength used is large. However, when a smaller number of grid points per wavelength is used excessive reflections occur which tend to overwhelm the computed solutions. Excessive reflections are due to the incompability between the governing finite difference equation and the Bayliss-Turkel radiation boundary condition. The Bayliss-Turkel radiation boundary condition was developed from the asymptotic solution of the partial differential equation. To obtain compatibility, the radiation boundary condition should be constructed from the asymptotic solution of the finite difference equation instead. Examples are provided using the improved radiation boundary condition based on the asymptotic solution of the governing finite difference equation. The computed results are free of reflections even when only five grid points per wavelength are used. The improved radiation boundary condition has also been tested for problems with complex acoustic sources and sources embedded in a uniform mean flow. The present method of developing a radiation boundary condition is also applicable to higher order finite difference schemes. In all these cases no reflected waves could be detected. The use of finite difference approximation inevita bly introduces anisotropy into the governing field equation. The effect of anisotropy is to distort the directional distribution of the amplitude and phase of the computed solution. It can be quite large when the number of grid points per wavelength used in the computation is small. A way to correct this effect is proposed. The correction factor developed from the asymptotic solutions is source independent and, hence, can be determined once and for all. The effectiveness of the correction factor in providing improvements to the computed solution is demonstrated in this paper.

Tam, Christopher K. W.↗
NASA NTRS1994

A method of measuring the amplitude-modulated vacuum field near a conducting mirror

Electromagnetic fields of the vacuum mode near a conducting mirror are modified with respect to those in free space, with their amplitudes having a sinusoidal spatial dependence from the mirror. Therefore if we combine this spatially amplitude-modulated vacuum field mode and intense coherent light with a beam splitter, we may detect this fluctuation of the vacuum mode in a homodyne detection scheme. It will give a new method to produce squeezed states of light with a single mirror placed close to an unused port of a beam splitter. We show that the amplitude fluctuation of the combined light can be reduced by a factor of 2 below that of the coherent light. We also discuss the limitations due to the finite line width of the laser and the effective absorption length of the photodiodes.

Youn, Sun-Hyun↗
NASA NTRS1992

The NASA/Industry Design Analysis Methods for Vibrations (DAMVIBS) Program - Boeing helicopters airframe finite element modeling

Finite-element modeling of the airframe vibration of the Army/Boeing CH-47D helicopter is conducted with comparisons to experimental data in an effort to improve the design process. A NASTRAN FEM is developed that is fully representative of the test configuration and includes the support fixture, shakers, and the aircraft/shaker suspension system. The analysis is conducted with specific attention given to the prediction of reasonable forced amplitudes throughout the airframe. Reasonable correlation is noted between the FEM and experimental results, although improved correlation can be obtained by including more accurate damping values and secondary effects such as stringer shear loading. It is shown that the general stress model does not provide an adequate dynamic analysis on which to base design improvements. A more detailed model is required that emphasizes highly detailed helicopter elements and employs a finer mesh particularly in the description of the mass distribution.

Gabel, R.↗
NASA NTRS2005

Convection in Ice I with Composite Newtonian/Non-Newtonian Rheology: Application to the Icy Galilean Satellites

Ice I exhibits a complex rheology at temperature and pressure conditions appropriate for the interiors of the outer ice I shells of Europa, Ganymede, and Callisto. We use numerical methods to determine the conditions required to trigger convection in an ice I shell with a stress-, temperature-, and grain-size-dependent rheology measured in laboratory experiments by Goldsby and Kohlstedt [2001] (henceforth GK2001). Triggering convection from an initially conductive ice shell with a non-Newtonian rheology for ice I requires that a finite-amplitude temperature perturbation be issued to the ice shell [2]. Here, we characterize the amplitude and wavelength of temperature perturbation required to initiate convection in the outer ice I shells of Europa, Ganymede, and Callisto using the GK2001 rheology for a range of ice grain sizes.

Barr, Amy C.↗