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 217 records · Page 12

Candle Flames in Non-Buoyant Atmospheres

This paper addresses the behavior of a candle flame in a long-duration, quiescent microgravity environment both on the space Shuttle and the Mir Orbiting Station (OS). On the Shuttle, the flames became dim blue after an initial transient where there was significant yellow (presumably soot) in the flame. The flame lifetimes were typically less than 60 seconds. The safety-mandated candlebox that contained the candle flame inhibited oxygen transport to the flame and thus limited the flame lifetime. 'Me flames on the Mir OS were similar, except that the yellow luminosity persisted longer into the flame lifetime because of a higher initial oxygen concentration. The Mir flames burned for as long as 45 minutes. The difference in the flame lifetime between the Shuttle and Mir flames was primarily the redesigned candlebox that did not inhibit oxygen transport to the flame. In both environments, the flame intensity and the height-to-width ratio gradually decreased as the ambient oxygen content in the sealed chamber slowly decreased. Both sets of experiments showed spontaneous, axisymmetric flame oscillations just prior to extinction. The paper also presents a numerical model of candle flame. The model is detailed in the gas-phase, but uses a simplified liquid/wick phase. 'Me model predicts a steady flame with a shape and size quantitatively similar to the Shuttle and Mir flames. ne model also predicts pre-extinction flame oscillations if the decrease in ambient oxygen is small enough.

Dietrich, D. L.↗

Dynamic Characteristics of the Partially Filled Rotating Dewar of the Gravity Probe-B Spacecraft

The dynamical behavior of fluids, in particular the effect of surface tension on partially-filled rotating fluids (cryogenic liquid helium and helium vapor) in a full scale Gravity Probe-B Spacecraft propellant dewar tank imposed by various frequencies of gravity jitters have been investigated. The study of the liquid-vapor interface oscillations due to various frequencies of gravity jitter under different dewar rotating speeds and different levels of background gravity have been numerically simulated. Results disclose the time sequence evolution for the excitation of large amplitude oscillation waves at the liquid-vapor interface. The study shows that slosh waves excited inside the spacecraft propellant tank are characterized by the lowest frequency of the waves initiated, frequencies of the gravity jitters imposed on the propellant system, the levels of background gravity environment, and dewar rotating speeds. Conditions for suppression and amplification of the slosh waves are discussed. It also shows that fluid stress distribution exerted on the walls of the rotating dewar are closely related to the characteristics of slosh waves excited on the liquid-vapor interface in the rotating dewar tank. This can provide a set of data leading toward the control of spacecraft imbalance caused by the uneven fluid stress distribution from slosh waves.

Hung, R. J.↗

The varieties of symmetric stellar rings and radial caustics in galaxy disks

Numerical, restricted three-body and analytic calculations are used to study the formation and propagation of cylindrically symmetric stellar ring waves in galaxy disks. It is shown that such waves can evolve in a variety of ways, depending on the amplitude of the perturbation and the potential of the target galaxy. Rings can thicken as they propagate outward, remain at a nearly constant width, or be pinched off at large radii. Multiple, closely spaced rings can result from a low-amplitude collision, while an outer ring can appear well-separated from overlapping inner rings or an apparent lens structure in halo-dominated potentials. All the single-encounter rings consist of paired fold caustics. The simple, impulsive, kinematic oscillation equations appear to provide a remarkably accurate model of the numerical simulations. Simple analytic approximations to these equations permit very good estimates of oscillation periods and amplitudes, the evolution of ring widths, and ring birth and propagation characteristics.

Struck-Marcell, Curtis↗

Theoretical studies of Resonance Enhanced Stimulated Raman Scattering (RESRS) of frequency doubled Alexandrite laser wavelength in cesium vapor

The third phase of research will focus on the propagation and energy extraction of the pump and SERS beams in a variety of configurations including oscillator structures. In order to address these questions a numerical code capable of allowing for saturation and full transverse beam evolution is required. The method proposed is based on a discretized propagation energy extraction model which uses a Kirchoff integral propagator coupled to the three level Raman model already developed. The model will have the resolution required by diffraction limits and will use the previous density matrix results in the adiabatic following limit. Owing to its large computational requirements, such a code must be implemented on a vector array processor. One code on the Cyber is being tested by using previously understood two-level laser models as guidelines for interpreting the results. Two tests were implemented: the evolution of modes in a passive resonator and the evolution of a stable state of the adiabatically eliminated laser equations. These results show mode shapes and diffraction losses for the first case and relaxation oscillations for the second one. Finally, in order to clarify the computing methodology used to exploit the speed of the Cyber's computational speed, the time it takes to perform both of the computations previously mentioned to run on the Cyber and VAX 730 must be measured. Also included is a short description of the current laser model (CAVITY.FOR) and a flow chart of the test computations.

Lawandy, Nabil M.↗

Extension of a semi-implicit shock-capturing algorithm for 3-D fully coupled, chemically reacting flows in generalized coordinates

A semiimplicit high-resolution shock-capturing method for multidimensional systems of hyperbolic conservation laws with stiff source terms has been developed by Yee and Shinn (1987). The goal of this work is to extend this method to solve the three-dimensional fully coupled Navier-Stokes equations for a hypersonic chemically reacting flow in generalized coordinates. In this formulation, the global continuity equation was replaced by all the species continuity equations. The shock-capturing technique is a second-order-accurate, symmetric total-variation-diminishing method which accounts fully and directly for the coupling among the fluid and all the species. To verify the current approach, it was implemented into an existing computer code which contained the MacCormack method. Test results for a five-species reacting flow are shown to be oscillation-free around the shock, and the time spent per iteration only doubles when compared to the result using classical way of supplying numerical dissipation. The extra computation is more than justified by the elimination of spurious oscillation and nonlinear instability associated with the classical shock-capturing schemes in computing hypersonic reacting flows.

Shinn, J. L.↗

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗

Present obliquity oscillations of Mars - Fourth-order accuracy in orbital e and I

A long period analysis of solar system orbital evolution, correct to fourth order in orbital eccentricities and inclinations (Bretagnon, 1974), and an improved value of the planet's moment of inertia (Reasenberg, 1977) have been incorporated in a recalculation of the obliquity oscillations of Mars. A linearized solution predicts a maximum oscillation amplitude of 13.6 deg centered on a long-term average value of 24.4 deg. A numerical integration of the obliquity for the past 10,000,000 years is also presented. Epochs of minimal oscillation like the present occur at intervals of the order of 4 m.y.

Ward, W. R.↗

Evaluation of Turbulence Models for Unsteady Flows of an Oscillating Airfoil

Unsteady flowfields of a two-dimensional oscillating airfoil are calculated using an implicit, finite-difference, Navier Stokes numerical scheme. Five widely used turbulence models are used with the numerical scheme to assess the accuracy and suitability of the models for simulating the retreating blade stall of helicopter rotor in forward flight. Three unsteady flow conditions corresponding to an essentially attached flow, light-stall, and deep-stall cases of an oscillating NACA 0015 wing experiment were chosen as test cases for computations. Results of unsteady airloads hysteresis curves, harmonics of unsteady pressures, and instantaneous flowfield patterns are presented. Some effects of grid density, time-step size, and numerical dissipation on the unsteady solutions relevant to the evaluation of turbulence models are examined. Comparison of unsteady airloads with experimental data show that all models tested are deficient in some sense and no single model predicts airloads consistently and in agreement with experiment for the three flow regimes. The chief findings are that the simple algebraic model based on the renormalization group theory (RNG) offers some improvement over the Baldwin Lomax model in all flow regimes with nearly same computational cost. The one-equation models provide significant improvement over the algebraic and the half-equation models but have their own limitations. The Baldwin-Barth model overpredicts separation and underpredicts reattachment. In contrast, the Spalart-Allmaras model underpredicts separation and overpredicts reattachment.

Srinivasan, G. R.↗

Thermal acoustic oscillations, volume 2

A number of thermal acoustic oscillation phenomena and their effects on cryogenic systems were studied. The conditions which cause or suppress oscillations, the frequency, amplitude and intensity of oscillations when they exist, and the heat loss they induce are discussed. Methods of numerical analysis utilizing the digital computer were developed for use in cryogenic systems design. In addition, an experimental verification program was conducted to study oscillation wave characteristics and boiloff rate. The data were then reduced and compared with the analytical predictions.

Spradley, L. W.↗

Finite element analysis of unsteady incompressible flow around an oscillating obstacle of arbitrary shape.

An analytical procedure based on Navier-Stokes equations was developed for representing unsteady flow patterns around oscillating obstacles. A variational formulation of the Helmholtz vorticity equation was discretized in finite element form and integrated numerically. At each step of the numerical integration the velocity field around the obstacle was determined from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints at each time step of the numerical integration. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗

Computational aspects of unsteady flows

The calculation of unsteady flows and the development of numerical methods for solving unsteady boundary layer equations and their application to the flows around important configurations such as oscillating airfoils are presented. A brief review of recent work is provided with emphasis on the need for numerical methods which can overcome possible problems associated with flow reversal and separation. The zig-zag and characteristic box schemes are described in this context, and when embodied in a method which permits interaction between solutions of inviscid and viscous equations, the characteristic box scheme is shown to avoid the singularity associated with boundary layer equations and prescribed pressure gradient. Calculations were performed for a cylinder started impulsively from rest and oscillating airfoils. The results are presented and discussed. It is conlcuded that turbulence models based on an algebraic specification of eddy viscosity can be adequate, that location of translation is important to the calculation of the location of flow separation and, therefore, to the overall lift of an oscillating airfoil.

Cebeci, T.↗

Diatomic gasdynamic lasers

Predictions from a numerical model of the vibrational relaxation of anharmonic diatomic oscillators in supersonic expansions are used to show the extent to which the small anharmonicity of gases like CO can cause significant overpopulations of upper vibrational states. When mixtures of CO and N2 are considered, radiative gain on many of the vibration-rotation transitions of CO is predicted. Experiments are described that qualitatively verify the predictions by demonstrating laser oscillation in CO-N2 expansions. The resulting CO-N2 gasdynamic laser displays performance characteristics that equal or exceed those of similar CO2 lasers.

Mckenzie, R. L.↗

Diatomic gasdynamic lasers.

Predictions from a numerical model of the vibrational relaxation of anharmonic diatomic oscillators in supersonic expansions are used to show the extent to which the small anharmonicity of gases like CO can cause significant overpopulations of upper vibrational states. When mixtures of CO and N2 are considered, radiative gain on many of the vibration-rotation transitions of CO is predicted. Experiments are described that qualitatively verify the predictions by demonstrating laser oscillation in CO-N2 expansions. The resulting CO-N2 gasdynamic laser displays performance characteristics that equal or exceed those of similar CO2 lasers.

Mckenzie, R. L.↗

The linear two-dimensional stability of inviscid vortex streets of finite-cored vortices

The stability of two-dimensional infinitesimal disturbances of the steady inviscid incompressible flow produced by infinite rows of finite-cored particles is studied. Numerical results are obtained for the growth rate and oscillation frequencies of disturbances of arbitrary subharmonic wavenumber, and the stability boundaries are calculated. The stabilization of the pairing instability by finite area demonstrated by Saffman and Schatzman (1982) is confirmed, as is Kida's (1982) result that this is not the most unstable disturbance when the area is finite. However, the finite area does not stabilize the inviscid Karman vortex street of finite-area vortices to infinitesimal two-dimensional disturbances of arbitrary wavelength. The finite area is always unstable except for one isolated value of the aspect ratio which depends upon the size of the vortices.

Meiron, D. I.↗

Traveling waves and chaos in thermosolutal convection

Numerical experiments on two-dimensional thermosolutal convection reveal oscillations in the form of traveling, standing, modulated, and chaotic waves. Transitions between these wave forms and steady convection are investigated and compared with theory. Such rich nonlinear behavior is possible in fluid layers of wide horizontal extent, and provides an explanation for waves observed in recent laboratory experiments with binary fluid mixtures.

Deane, A. E.↗

The Dynamics of Vapor Bubbles in Acoustic Pressure Fields

In spite of a superficial similarity with gas bubbles, the intimate coupling between dynamical and thermal processes confers to oscillating vapor bubbles some unique characteristics. This paper examines numerically the validity of some asymptotic-theory predictions such as the existence of two resonant radii and a limit size for a given sound amplitude and frequency. It is found that a small vapor bubble in a sound field of sufficient amplitude grows quickly through resonance and continues to grow thereafter at a very slow rate, seemingly indefinitely. Resonance phenomena therefore play a role for a few cycles at most, and reaching a limit size-if one exists at all-is found to require far more than several tens of thousands of cycles. It is also found that some small bubbles may grow or collapse depending on the phase of the sound field. The model accounts in detail for the thermo-fluid-mechanic processes in the vapor. In the second part of the paper, an approximate formulation valid for bubbles small with respect to the thermal penetration length in the vapor is derived and its accuracy examined, The present findings have implications for acoustically enhanced boiling heat transfer and other special applications such as boiling in microgravity.

Hao, Y.↗

Optimized Shapes of Ocsillating Resonators for Generating High-Amplitude Pressure Waves

It is well known that the resonator geometry strongly influences the resonant frequencies of an acoustical resonator and the generated nonlinear standing pressure waveform. Maximizing the ratio of maximum to minimum gas pressure at an end of an oscillating resonator by optimizing the cavity contour is investigated numerically. A quasi-Newton type scheme is used to find optimized axisymmetric resonator shapes to achieve the maximum pressure compression ratio. The acoustical field is solved using a one-dimensional model, and the resonance frequency shift and hysteresis effects are obtained through an automation scheme based on continuation methods. Results are presented from optimizing cone, horn-cone, and cosine resonator geometries. Significant performance improvement is found in the optimized shapes over others previously published. Different optimized shapes are found when starting with different initial guesses, indicating multiple local extrema. The numerical model is validated by comparing with the experimental results of a horn-cone shaped resonator.

Li, Xiao-Fan↗