Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Oscillation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Numerical solution of the incompressible Navier-Stokes equations

The current work is initiated in an effort to obtain an efficient, accurate, and robust algorithm for the numerical solution of the incompressible Navier-Stokes equations in two- and three-dimensional generalized curvilinear coordinates for both steady-state and time-dependent flow problems. This is accomplished with the use of the method of artificial compressibility and a high-order flux-difference splitting technique for the differencing of the convective terms. Time accuracy is obtained in the numerical solutions by subiterating the equations in psuedo-time for each physical time step. The system of equations is solved with a line-relaxation scheme which allows the use of very large pseudo-time steps leading to fast convergence for steady-state problems as well as for the subiterations of time-dependent problems. Numerous laminar test flow problems are computed and presented with a comparison against analytically known solutions or experimental results. These include the flow in a driven cavity, the flow over a backward-facing step, the steady and unsteady flow over a circular cylinder, flow over an oscillating plate, flow through a one-dimensional inviscid channel with oscillating back pressure, the steady-state flow through a square duct with a 90 degree bend, and the flow through an artificial heart configuration with moving boundaries. An adequate comparison with the analytical or experimental results is obtained in all cases. Numerical comparisons of the upwind differencing with central differencing plus artificial dissipation indicates that the upwind differencing provides a much more robust algorithm, which requires significantly less computing time. The time-dependent problems require on the order of 10 to 20 subiterations, indicating that the elliptical nature of the problem does require a substantial amount of computing effort.

Rogers, Stuart E.↗
NASA NTRS1986

Compton-heated winds and coronae above accretion disks. II - Instability and oscillations

The stability and evolution of windy accretion disks is investigated in detail. The basic disk evolution equations are briefly recapitulated, and an idealized analytic treatment of the wind and viscosity is used to show that steady disk flow is indeed unstable for sufficiently large ratio of the mass loss rate in the wind to the central accretion rate. Numerical solutions for a more realistic and appropriate expression for the mass loss rate and the standard ad hoc alpha model prescription for the viscosity are presented. The application of these results to real systems with Compton-heated winds is discussed, and a general formula for the oscillation period is given. The prediction is compared with observed periodic behavior of Galactic X-ray sources and AGNs.

Shields, G. A.↗
NASA NTRS1979

Computation of aerodynamic interference effects on oscillating airfoils with controls in ventilated subsonic wind tunnels

Lift interference effects are discussed based on Bland's (1968) integral equation. A mathematical existence theory is utilized for which convergence of the numerical method has been proved for general (square-integrable) downwashes. Airloads are computed using orthogonal airfoil polynomial pairs in conjunction with a collocation method which is numerically equivalent to Galerkin's method and complex least squares. Convergence exhibits exponentially decreasing error with the number n of collocation points for smooth downwashes, whereas errors are proportional to 1/n for discontinuous downwashes. The latter can be reduced to 1/n to the m+1 power with mth-order Richardson extrapolation (by using m = 2, hundredfold error reductions were obtained with only a 13% increase of computer time). Numerical results are presented showing acoustic resonance, as well as the effect of Mach number, ventilation, height-to-chord ratio, and mode shape on wind-tunnel interference. Excellent agreement with experiment is obtained in steady flow, and good agreement is obtained for unsteady flow.

Fromme, J. A.↗
NASA NTRS1994

Nonlinear dynamics of drops and bubbles and chaotic phenomena

Nonlinear phenomena associated with the dynamics of free drops and bubbles are investigated analytically, numerically and experimentally. Although newly developed levitation and measurement techniques have been implemented, the full experimental validation of theoretical predictions has been hindered by interfering artifacts associated with levitation in the Earth gravitational field. The low gravity environment of orbital space flight has been shown to provide a more quiescent environment which can be utilized to better match the idealized theoretical conditions. The research effort described in this paper is a closely coupled collaboration between predictive and guiding theoretical activities and a unique experimental program involving the ultrasonic and electrostatic levitation of single droplets and bubbles. The goal is to develop and to validate methods based on nonlinear dynamics for the understanding of the large amplitude oscillatory response of single drops and bubbles to both isotropic and asymmetric pressure stimuli. The first specific area on interest has been the resonant coupling between volume and shape oscillatory modes isolated gas or vapor bubbles in a liquid host. The result of multiple time-scale asymptotic treatment, combined with domain perturbation and bifurcation methods, has been the prediction of resonant and near-resonant coupling between volume and shape modes leading to stable as well as chaotic oscillations. Experimental investigations of the large amplitude shape oscillation modes of centimeter-size single bubbles trapped in water at 1 G and under reduced hydrostatic pressure, have suggested the possibility of a low gravity experiment to study the direct coupling between these low frequency shape modes and the volume pulsation, sound-radiating mode. The second subject of interest has involved numerical modeling, using the boundary integral method, of the large amplitude shape oscillations of charged and uncharged drops in the presence of a static or time-varying electric field. Theoretically predicted non linearity in the resonant frequency of the fundamental quadrupole mode has been verified by the accompanying experimental studies. Additional phenomena such as hysteresis in the frequency response of ultrasoncially levitated droplets in the presence of a time varying electric field, and mode coupling in the oscillations of ultrasonically modulated droplets, have also been uncovered. One of the results of this ground-based research has been the identification and characterization of phenomena strictly associated with the influence of the gravitational field. This has also allowed us to identify the specific requirements for potential microgravity investigations yielding new information not obtainable on Earth.

Trinh, Eugene H.↗
NASA NTRS2004

High-Order Shock-Capturing Methods for Modeling Dynamics of the Solar Atmosphere

We use one-dimensional high-order central shock capturing numerical methods to study the response of various model solar atmospheres to forcing at the solar surface. The dynamics of the atmosphere is modeled with the Euler equations in a variable-sized flux tube in the presence of gravity. We study dynamics of the atmosphere suggestive of spicule formation and coronal oscillations. These studies are performed on observationally-derived model atmospheres above the quiet sun and above sunspots. To perform these simulations, we provide a new extension of existing second- and third- order shock-capturing methods to irregular grids. We also solve the problem of numerically maintaining initial hydrostatic balance via the introduction of new variables in the model equations and a careful initialization mechanism. We find several striking results: all model atmospheres respond to a single impulsive perturbation with several strong shock waves consistent with the rebound-shock model. These shock waves lift material and the transition region well into the initial corona, and the sensitivity of this lift to the initial impulse depends non-linearly on the details of the atmosphere model. We also reproduce an observed 3-minute coronal oscillation above sunspots compared to 5-minute oscillations above the quiet sun.

Bryson, Steve↗
NASA NTRS1996

A new methodology to determine kinetic parameters for one- and two-step chemical models

In this paper, a new methodology to determine kinetic parameters for simple chemical models and simple transport properties classically used in DNS of premixed combustion is presented. First, a one-dimensional code is utilized to performed steady unstrained laminar methane-air flame in order to verify intrinsic features of laminar flames such as burning velocity and temperature and concentration profiles. Second, the flame response to steady and unsteady strain in the opposed jet configuration is numerically investigated. It appears that for a well determined set of parameters, one- and two-step mechanisms reproduce the extinction limit of a laminar flame submitted to a steady strain. Computations with the GRI-mech mechanism (177 reactions, 39 species) and multicomponent transport properties are used to validate these simplified models. A sensitivity analysis of the preferential diffusion of heat and reactants when the Lewis number is close to unity indicates that the response of the flame to an oscillating strain is very sensitive to this number. As an application of this methodology, the interaction between a two-dimensional vortex pair and a premixed laminar flame is performed by Direct Numerical Simulation (DNS) using the one- and two-step mechanisms. Comparison with the experimental results of Samaniego et al. (1994) shows a significant improvement in the description of the interaction when the two-step model is used.

Mantel, T.↗
NASA NTRS2014

Numerical Simulation of a High-Lift Configuration with Embedded Fluidic Actuators

Numerical simulations have been performed for a vertical tail configuration with deflected rudder. The suction surface of the main element of this configuration is embedded with an array of 32 fluidic actuators that produce oscillating sweeping jets. Such oscillating jets have been found to be very effective for flow control applications in the past. In the current paper, a high-fidelity computational fluid dynamics (CFD) code known as the PowerFLOW(Registered TradeMark) code is used to simulate the entire flow field associated with this configuration, including the flow inside the actuators. The computed results for the surface pressure and integrated forces compare favorably with measured data. In addition, numerical solutions predict the correct trends in forces with active flow control compared to the no control case. Effect of varying yaw and rudder deflection angles are also presented. In addition, computations have been performed at a higher Reynolds number to assess the performance of fluidic actuators at flight conditions.

Vatsa, Veer N.↗
NASA NTRS1975

Further developments in the prediction of oscillatory aerodynamics in mixed transonic flow

This paper presents improvements made to the oscillatory transonic aerodynamic kernel function method given in AIAA Paper No. 74-359. The improved method is fully developed for flows with continuously varying local Mach number and oscillating normal shocks. The boundary conditions developed for oscillating shocks require that a line doublet be added at the foot of the shock which represents the loads induced by shock movement. Numerical results show good agreement with experiment and significant deviation from linear aerodynamic theory. The method is also used to predict wind-tunnel wall interference effects.-

Cunningham, A. M., Jr.↗
NASA NTRS1975

Solutions to Kuessner's integral equation in unsteady flow using local basis functions

The computational procedure and numerical results are presented for a new method to solve Kuessner's integral equation in the case of subsonic compressible flow about harmonically oscillating planar surfaces with controls. Kuessner's equation is a linear transformation from pressure to normalwash. The unknown pressure is expanded in terms of prescribed basis functions and the unknown basis function coefficients are determined in the usual manner by satisfying the given normalwash distribution either collocationally or in the complex least squares sense. The present method of solution differs from previous ones in that the basis functions are defined in a continuous fashion over a relatively small portion of the aerodynamic surface and are zero elsewhere. This method, termed the local basis function method, combines the smoothness and accuracy of distribution methods with the simplicity and versatility of panel methods. Predictions by the local basis function method for unsteady flow are shown to be in excellent agreement with other methods. Also, potential improvements to the present method and extensions to more general classes of solutions are discussed.

Fromme, J. A.↗
NASA NTRS1976

Natural modes of oscillation of rotating flexible structures about nontrivial equilibrium

This paper is concerned with the vibrational characteristics of rotating flexible structures in the neighborhood of steady rotation. It is assumed that the structure undergoes deformations during the steady rotation, so that the equilibrium is nontrivial. The object of the paper is to formulate the eigenvalue problem associated with small oscillations of the structure about the nontrivial equilibrium and to develop efficient methods to compute the system natural frequencies and modal vectors. A numerical example showing the procedure for the calculation of spacecraft modes is presented.

Meirovitch, L.↗
NASA NTRS1978

Some observations on the mechanism of aircraft wing rock

A pressure scale model of Northrop F-5A was tested in NASA Ames Research Center Eleven-Foot Transonic Tunnel to simulate the wing rock oscillations in a transonic maneuver. For this purpose, a flexible model support device was designed and fabricated which allowed the model to oscillate in roll at the scaled wing rock frequency. Two tunnel entries were performed to acquire the pressure (steady state and fluctuating) and response data when the model was held fixed and when it was excited by flow to oscillate in roll. Based on these data, a limit cycle mechanism was identified which supplied energy to the aircraft model and caused the Dutch roll type oscillations, commonly called wing rock. The major origin of the fluctuating pressures which contributed to the limit cycle was traced to the wing surface leading edge stall and the subsequent lift recovery. For typical wing rock oscillations, the energy balance between the pressure work input and the energy consumed by the model aerodynamic and mechanical damping was formulated and numerical data presented.

Hwang, C.↗
NASA NTRS1981

Computation of separated flows by a vortex-tracing algorithm

Numerical solutions for two-dimensional, time-dependent, separated flows around bodies are obtained, using a new version of the vortex method. This method provides an efficient representation of flows involving large regions of separation. The modifications incorporated in the new version, which improve its accuracy, versatility, and computing speed, are described. The computer cost is only of the order of the 3/2 power on N, instead of N-squared, for each step with N vortices. Arbitrary shapes can be treated; a conformal mapping is not required. Special attention is paid to the viscous character of the solution and to the accurate computation of the pressure distribution at the body surface. The vortex solution for the outer flow is coupled to an inner solution for the attached part of the boundary layer. Numerical results are presented for several bluff bodies exhibiting dependence on Reynolds number, for stationary airfoils under steady or transient conditions and for oscillating airfoils, including dynamic stalls. These results are compared with other available results, analytical or experimental, and demonstrate the enhanced reliability and accuracy of the improved method.

Spalart, P. R.↗
NASA NTRS1985

On the wake hazard alleviation associated with roll oscillations of wake-generating aircraft

An explanation is provided for the difference in wake vortex alleviation achieved by roll oscillations in flight tests with B-747 and L-1011 transport aircraft. Both aircraft had their landing flaps extended and several spoilers deployed. Numerical analysis shows that the growth in amplitude of the initial waves in the vortex filament is brought about by the sinusoidal instability. In the case of the B-747, growth is enhanced by a vortex whose strength is about the same as the tip vortex which is shed near the fuselage by the inboard end of the flaps. Conversely, the L-1011 is estimated to shed a negligible fuselage vortex and to have a relatively strong wingtip vortex. These characteristics bring about a rotation and amplification of the initial waves in the vortex filaments in the wake of the B-747 but not in the L-1011 vortex wake. An aircraft following the B-747 would then experience only intermittent encounters with the intense parts of the wake vortices so that the time-averaged wake-induced rolling moment is substantially reduced.

Rossow, V. J.↗
NASA NTRS1983

Advanced algorithm for orbit computation

Computational and analytical techniques which simplify the solution of complex problems in orbit mechanics, Astrodynamics and Celestial Mechanics were developed. The major tool of the simplification is the substitution of transformations in place of numerical or analytical integrations. In this way the rather complicated equations of orbit mechanics might sometimes be reduced to linear equations representing harmonic oscillators with constant coefficients.

Szenbehely, V.↗
NASA NTRS1988

Stability of swirling gas flows

The stability of inviscid swirling gas flows to small nonaxisymmetric perturbations is considered. For small Brunt-Vaisala frequencies, the problem reduces to the classical Sturm-Liouville form and the oscillation theorem can be applied. The resulting necessary and sufficient stability condition is compared to various criteria in the literature and a limited numerical study of isothermal rigidly rotating Poiseuille flow. For given azimuthal and axial wavenumbers, it is found numerically that the higher inertial modes become unstable for successively lower Rossby numbers and that this sequence of critical values approaches the theoretical value from above.

Hultgren, Lennart S.↗
NASA NTRS1988

On the pulsational overstability in narrowly confined viscous rings

Viscous accretion disks and rings are subject to pulsational overstability if the effective viscous stress is a rapidly increasing function of the disks' (or rings') surface density. Such a phenomenon is induced by the transfer of energy from shear to normal mode oscillations through viscous stress. The paper examines the global response of viscous accretion disks or rings to pulsational overstabilities. The analytic and preliminary numerical results indicate that all modes with wavelength longer than the scale height of the disk or ring are overstable. From the results, some examples of the complex time-dependent structure which may be induced by this overstability are illustrated. Finally, results are applied to narrow planetary rings to address the issues of maintenance of uniform precession and origin of eccentricity.

Papaloizou, J. C. B.↗
NASA NTRS1991

Anomalous temperature and zonal wind in the tropical upper stratosphere, 1982/1983

Observed temperature and zonal wind anomalies during 1982-1983 are described in relation to the 20-yr rocketsonde data record. Temperature data from the Nimbus 7 SAMS indicate that in mid-1982, cooling occurred in the upper and middle tropical stratosphere, starting shortly after the eruption of El Chichon (April 1982). Rocketsonde data from Kwajalein gave additional evidence of cooling and revealed that zonal winds in the upper and middle stratosphere were anomalous in 1982 and also in 1983. Rocketsonde observations were consistent with balance winds derived from SAMS. Upper-level cooling was linked to anomalous easterlies near the equator and therefore could not be interpreted as a radiative (nondynamical) response to El Chichon. The structure of an expected 2D dynamical response to aerosol heating is examined with a numerical model and compared to the data. It is suggested that internal sources of atmospheric variability (eg., extratropical forcing and quasi-biennial oscillation) contributed to the equatorial cooling observed in 1982.

Dunkerton, Timothy J.↗
NASA NTRS2000

Downward Link of Solar Activity Variations Through Wave Driven Equatorial Oscillations (QBO and SAO)

Signatures of the 11-year solar activity/irradiance cycle are observed in the Quasi Biennial Oscillation (QBO) of the lower stratosphere. At these altitudes, the QBO is understood to be the result of "downward control" exerted by the wave mean flow interactions that drive the phenomenon. It is reasonable then to speculate that the QBO is a natural conduit to lower altitudes of solar activity variations in radiance (SAV). To test this hypothesis, we conducted experiments with a 2D version of our Numerical Spectral Model that incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GW). To account for the SAV, we change the solar heating rate on a logarithmic scale from 0.1% at the surface to 1% at 50 kin to 10% at 100 km. With the same GW flux, we then conduct numerical experiments to evaluate the magnitude of the solar activity irradiance effect (SAE) on the zonal circulation at low latitudes. The numerical results obtained show that, under certain conditions, the SAE is significant in the zonal circulation and does extend to lower altitudes where the SAV is small. The differences in the wind velocities can be as large as 5 m/s at 20 kin. We carried out two numerical experiments with integrations over more than 20 years: 1) With the QBO period "tuned" to be 30 months, of academic interest but instructive, the seasonal cycle in the solar forcing [through the Semi-annual Oscillation (SAO)] acts as a strong pacemaker to produce a firm lock on the period and phase of the QBO. The SAE then shows up primarily as a distinct but relatively weak amplitude modulation. 2) With the QBO period between 30 and 34 (or less than 30, presumably) months, the seasonal phase lock is weak compared with (1). The SAV in the seasonal cycle then causes variations in the QBO period and phase, and this amplifies the SAE to produce relatively large variations in the wind field. We conclude that, under realistic conditions as in (2), the solar seasonal forcing, with activity variations in radiance, causes the phase and period of the QBO to change to produce a relatively large solar activity effect in the zonal circulation at low latitudes and low altitudes.

Mengel, J. G.↗