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 397 records · Page 22

Middle Atmosphere Transport Properties of Assimilated Datasets

One of the most compelling reasons for performing data assimilation in the middle atmosphere is to obtain global, balanced datasets for studies of trace gas transport and chemistry. This is a major motivation behind the Goddard Earth observation System-Data Assimilation System (GEOS-DAS). Previous studies have shown that while this and other data assimilation systems can generally obtain good estimates of the extratropical rotational velocity field, the divergent part of the dynamical field is deficient; this impacts the "residual circulation" and leads to spurious trace gas transport on seasonal and interannual timescales. These problems are impacted by the quality and the method of use of the observational data and by deficiencies in the atmospheric general circulation model. Whichever the cause at any place and time, the "solution" is to introduce non-physical forcing terms into the system (the so-called incremental analysis updates); these can directly (thermal) or indirectly (mechanical) affect the residual circulation. This paper will illustrate how the divergent circulation is affected by deficiencies in both observations and models. Theoretical considerations will be illustrated with examples from the GEOS-DAS and from simplified numerical experiments. These are designed to isolate known problems, such as the inability of models to sustain a quasi-biennial oscillation and sparse observational constraints on tropical dynamics, or radiative inconsistencies in the presence of volcanic aerosols.

Pawson, Steven↗

Middle Atmospheric Transport Properties of Assimilated Datasets

One of the most compelling reasons for performing data assimilation in the middle atmosphere is to obtain global, balanced datasets for studies of trace gas transport and chemistry. This is a major motivation behind the Goddard Earth observation System-Data Assimilation System (GEOS-DAS). Previous studies have shown that while this and other data assimilation systems can generally obtain good estimates of the extratropical rotational velocity field, the divergent part of the dynamical field is deficient; this impacts the "residual circulation" and leads to spurious trace gas transport on seasonal and interannual timescales. These problems are impacted by the quality and the method of use of the observational data and by deficiencies in the atmospheric general circulation model. Whichever the cause at any place and time, the "solution" is to introduce non-physical forcing terms into the system (the so-called incremental analysis updates); these can directly (thermal) or indirectly (mechanical) affect the residual circulation. This paper will illustrate how the divergent circulation is affected by deficiencies in both observations and models. Theoretical considerations will be illustrated with examples from the GEOS-DAS and from simplified numerical experiments. These are designed to isolate known problems, such as the inability of models to sustain a quasi-biennial oscillation and sparse observational constraints on tropical dynamics, or radiative inconsistencies in the presence of volcanic aerosols.

Pawson, Steven↗

Numerical Investigation of 'Transonic Resonance' with a Convergent-Divergent Nozzle

At pressure ratios lower than the design value, convergent-divergent (C-D) nozzles often undergo a flow resonance accompanied by the emission of acoustic tones. The phenomenon, driven by the unsteady shock within the divergent section of the nozzle, has been studied experimentally by Zaman et al. In this paper, the space-time conservation element solution element (CE/SE) method is employed to numerically investigate the phenomenon. The computations are performed for a given nozzle geometry for several different pressure ratios. Sustained 'limit cycle' oscillations are encountered in all cases. The oscillation frequencies, their variation with pressure ratio including a 'stage jump', agree well with the experimental results. The unsteady flow data confirm that stage 1 of the resonance (fundamental) involves a one-quarter standing wave while stage 2 (third harmonic) involves a three-quarter standing wave within the divergent section of the nozzle. Details of the shock motion, and the flow and near acoustic field, are documented for one case each of stages 1 and 2.

Loh, Ching Y.↗

Large-Scale Alfvenic Impulses on the Sun: How They Are Generated and What We Learn From Them

NASA GSFC The Sun's atmosphere hosts a wide variety of magnetosonic disturbances. These wave modes are detected, almost exclusively, by examining images of the Sun's magnetic atmosphere and looking for propagating distortions. Although none of the Sun's plasma parameters are measured directly, we derive a great deal of information from these observations. In fact, by modeling these propagating disturbances, we may be able to derive the most accurate estimates plasma parameters. From observations absorption, refraction, reflection, and coupling of numerous wave modes, we advance our knowledge of the Sun's magnetic field, temperature, density, and current. The Sun's continuous oscillation, coronal mass ejections, flares, and other dynamic phenomena can produce wave disturbances which are observable from near-Earth space. Several of these disturbances have been traced from the inner corona out into the heliosphere. From the generation of these disturbances, we are able to learn about the phenomena which create them as well as the media through which they re-propagating. The presentation will include a discussion of the generation of Alfvenic disturbances on the Sun, ways we observe these disturbances, and how recent advances in modeling and analysis have brought us closer to determining solar in situ parameters.

Thompson, Barbara↗

Observations of Intermediate-Scale Diurnal Waves in the Equatorial Mesosphere and Lower Thermosphere

The purpose of this study is to seek observational evidence for diurnal, vertically propagating inertia-gravity waves (IGWs) in the mesosphere and lower thermosphere (MLT). Numerical modeling studies indicate that vertically propagating IGWs excited by tropical heating provide a significant source of momentum for the semiannual oscillation (SAO) in the equatorial zonal mean winds. The power spectrum of these waves has a strong diurnal component. We analyze global patterns of ascending and descending node differences in MLT satellite temperatures, which are assumed in this study to be proxies for waves of diurnal period. Juxtaposed upon the planetary-scale features are localized variations with longitudinal wavelengths ranging from 25 deg to 50 deg, and with vertical wavelengths between 13 and 25 km. These intermediate-scale variations are spatially coherent, and persist for several weeks. Their amplitudes generally increase with altitude, while their phase structures suggest both eastward and westward propagation. The variance of wave numbers 9-17 in the MLT is examined in relation to the underlying stratospheric zonal mean zonal gradient winds. Stronger variances generally coincide with periods where the underlying zonal mean winds are relatively light, or unidirectional. The weak inverse association between variance strength and zonal wind magnitude is suggestive of a wave filtering mechanism.

Lieberman, R. S.↗

Weather Research and Forecasting Model Wind Sensitivity Study at Edwards Air Force Base, CA

NASA prefers to land the space shuttle at Kennedy Space Center (KSC). When weather conditions violate Flight Rules at KSC, NASA will usually divert the shuttle landing to Edwards Air Force Base (EAFB) in Southern California. But forecasting surface winds at EAFB is a challenge for the Spaceflight Meteorology Group (SMG) forecasters due to the complex terrain that surrounds EAFB, One particular phenomena identified by SMG is that makes it difficult to forecast the EAFB surface winds is called "wind cycling". This occurs when wind speeds and directions oscillate among towers near the EAFB runway leading to a challenging deorbit bum forecast for shuttle landings. The large-scale numerical weather prediction models cannot properly resolve the wind field due to their coarse horizontal resolutions, so a properly tuned high-resolution mesoscale model is needed. The Weather Research and Forecasting (WRF) model meets this requirement. The AMU assessed the different WRF model options to determine which configuration best predicted surface wind speed and direction at EAFB, To do so, the AMU compared the WRF model performance using two hot start initializations with the Advanced Research WRF and Non-hydrostatic Mesoscale Model dynamical cores and compared model performance while varying the physics options.

Watson, Leela R.↗

Genesis of Twin Tropical Cyclones as Revealed by a Global Mesoscale Model: The Role of Mixed Rossby Gravity Waves

In this study, it is proposed that twin tropical cyclones (TCs), Kesiny and 01A, in May 2002 formed in association with the scale interactions of three gyres that appeared as a convectively coupled mixed Rossby gravity (ccMRG) wave during an active phase of the Madden-Julian Oscillation (MJO). This is shown by analyzing observational data, including NCEP reanalysis data and METEOSAT 7 IR satellite imagery, and performing numerical simulations using a global mesoscale model. A 10-day control run is initialized at 0000 UTC 1 May 2002 with grid-scale condensation but no sub-grid cumulus parameterizations. The ccMRG wave was identified as encompassing two developing and one non-developing gyres, the first two of which intensified and evolved into the twin TCs. The control run is able to reproduce the evolution of the ccMRG wave and thus the formation of the twin TCs about two and five days in advance as well as their subsequent intensity evolution and movement within an 8-10 day period. Five additional 10-day sensitivity experiments with different model configurations are conducted to help understand the interaction of the three gyres, leading to the formation of the TCs. These experiments suggest the improved lead time in the control run may be attributed to the realistic simulation of the ccMRG wave with the following processes: (1) wave deepening (intensification) associated with a reduction in wavelength and/or the intensification of individual gyres, (2) poleward movement of gyres that may be associated with boundary layer processes, (3) realistic simulation of moist processes at regional scales in association with each of the gyres, and (4) the vertical phasing of low- and mid-level cyclonic circulations associated with a specific gyre.

forecast↗

The QBO and weak external forcing by solar activity: A three dimensional model study

A better understanding is attempted of the physical mechanisms leading to significant correlations between oscillations in the lower and middle stratosphere and solar variability associated with the sun's rotation. A global 3-d mechanistic model of the middle atmosphere is employed to investigate the effects of minor artificially induced perturbations. The aim is to explore the physical mechanisms of the dynamical response especially of the stratosphere to weak external forcing as it may result from UV flux changes due to solar rotation. First results of numerical experiments dealing about the external forcing of the middle atmosphere by solar activity were presented elsewhere. Different numerical studies regarding the excitation and propagation of weak perturbations have been continued since then. The model calculations presented are made to investigate the influence of the quasi-biennial oscillation (QBO) on the dynamical response of the middle atmosphere to weak perturbations by employing different initial wind fields which represent the west and east phase of the QBO.

Dameris, M.↗

Study of unsteady flow field over a forward-looking endoatmospheric hit-to-kill interceptor

Forward-looking recessed aperture interceptor has significant aero-optical and aero-thermal advantages. Previous experimental studies have shown that the flow field in front of a forward-looking cavity is unsteady and the bow shock oscillates at the cavity fundamental resonant frequency. In this study, an advanced CFD code is applied to study the above unsteady phenomena. The code is first validated against the experiments and good comparisons are found. The numerical parametric study shows that the existence of oscillatory bow shock is very sensitive to the cavity geometry. At a FOV of 60 deg, the initial transient quickly dampens out to a steady state. With a decrease of FOV, an unsteady oscillatory flow field is sustained after initial transient and the amplitude of oscillation is a function of FOV. For FOV of 20 deg, the amplitude of pressure oscillation is 25 percent of the mean value in the cavity. For a FOV of 10 deg, it can be as high as 50 percent.

Yang, H. Q.↗

Numerical solution of the two-dimensional Euler equations by second-order upwind difference schemes

Two time-level, five-point explicit and implicit upwind difference schemes based on the characteristic flux difference splitting concept have been developed for the two-dimensional Euler equations. The method is conservative, second-order accurate in time and space, and general coordinate systems are used to treat complex geometries. Nonlinear flux limiters are employed to yield oscillation free sharp shock profiles. Upstream interpolation is used to yield a class of higher-order upwind schemes which closely mimic the locally one-dimensional method of characteristics (with fixed time intervals) through operator splitting. Numerical results have been obtained for a plane shock reflection and for flow over a circular arc in a channel. Characteristics of upwind TVD schemes, as applied to two-dimensional flows with embedded shocks are discussed.

Yang, J. Y.↗

El Nino

An attempt is made to provide the background for a coupled model of ENSO (El Nino-Southern Oscillation) with emphasis placed on the oceanography (i.e. on El Nino). Observations of the normal annual cycle in the Pacific and of the evolution of a typical El Nino event are reviewed, and a theory for the oceanography of El Nino is proposed. The influence of SST anomalies on the tropical atmosphere is assessed, and results from a numerical model for the coupled system able to generate El Nino events are presented. Implications for the real ENSO cycle are discussed. In both the model and nature, ENSO has the character of a relaxation oscillation of the coupled system, and its cycle is aperiodic. Results on the predictability of dynamical systems show the impossibility of predicting ahead several events.

Cane, M. A.↗

Analysis of three-dimensional unsteady flow around oscillating wings

A method based on the Navier-Stokes equations was developed for determining analytically the three-dimensional unsteady flow patterns around oscillating wings. The Helmholz vorticity transport equations were discretized in three-dimensional finite element form from a variational formulation and integrated numerically. At each time step of the numerical integration the velocity field was calculated from the representation of the three-dimensional wing by a system of optimized distribution of vortices in space. During the numerical integration of the vorticity transport equations the time-dependent boundary conditions on the wing were specified as external constraint conditions. Examples of obtained results describing the three-dimensional unsteady flow around a wing were presented.

Bratanow, T.↗

A numerical study of the stability of radiative shocks

Attention is given to the oscillatory instability of optically thin radiative shocks in time-dependent numerical calculations of accretion flows onto degenerate dwarfs. The present nonlinear calculations yield good quantitative agreement with the linear results obtained for oscillation frequencies, damping rates, and critical alpha-values. The fundamental mode and the first overtone in the shock radius and luminosity variations can be clearly identified, and evidence is sometimes seen for the second overtone. Time-dependent calculations are also performed which include additional physics relevant to degenerate dwarf accretion, such as electron thermal conduction, unequal electron and ion temperatures, Compton cooling, and relativistic corrections to the bremsstrahlung cooling law. All oscillatory modes are found to be damped, and hence stable, in the case of a 1-solar mass white dwarf accreting in spherical symmetry.

Imamura, J. N.↗

Flutter and oscillating air-force calculations for an airfoil in a two-dimensional supersonic flow

An account is given of the Possio theory of non stationary flow for small disturbances in a two-dimensional supersonic flow and of its application to determination of the aerodynamic forces on an oscillating airfoil. Further application is made to the problem of wing flutter in the degrees of freedom - torsion, bending, and aileron torsion. Numerical tables for flutter calculations are provided for numerous values of the Mach number greater than unity. Results for bending-torsion wing flutter are discussed. The static instabilities of divergence and aileron reversal are examined as is a one degree of freedom case of torsional oscillatory instability.

Garrick, I E↗

Flutter and oscillating air-force calculations for an airfoil in two-dimensional supersonic flow

A connected account is given of the Possio theory of non-stationary flow for small disturbances in a two-dimensional supersonic flow and of its application to the determination of the aerodynamic forces on an oscillating airfoil. Further application is made to the problem of wing flutter in the degrees of freedom - torsion, bending, and aileron rotations. Numerical tables for flutter calculations are provided for various values of the Mach number greater than unity. Results for bending-torsion wing flutter are shown in figures and are discussed. The static instabilities of divergence and aileron reversal are examined as is a one-degree-of-freedom case of torsional oscillatory instability.

Garrick, I E↗

A massively parallel computational approach to coupled thermoelastic/porous gas flow problems

A new computational scheme for coupled thermoelastic/porous gas flow problems is presented. Heat transfer, gas flow, and dynamic thermoelastic governing equations are expressed in fully explicit form, and solved on a massively parallel computer. The transpiration cooling problem is used as an example problem. The numerical solutions have been verified by comparison to available analytical solutions. Transient temperature, pressure, and stress distributions have been obtained. Small spatial oscillations in pressure and stress have been observed, which would be impractical to predict with previously available schemes. Comparisons between serial and massively parallel versions of the scheme have also been made. The results indicate that for small scale problems the serial and parallel versions use practically the same amount of CPU time. However, as the problem size increases the parallel version becomes more efficient than the serial version.

Shia, David↗

Some Numerical Experiments on Detonation Wave Propagation

In this paper we present the results of a series of numerical experiments done on the propagation and initiation of a detonation wave. The calculations are performed in one-dimension, with considerable grid resolution. Of particular interest are the following questions: (1) the nature of periodic and chaotic instabilities generated by the wave; (2) the influence of the grid resolution on these instabilities; (3) the influence of the 'quality' of the numerical scheme; and (4) the influence of 'noise'. In the calculations, we use a second-order Total Variation Diminishing (TVD) scheme as the basic numerical method, with grid spacings as low as a fraction of a micron. Detonations waves are generated at the closed end of a tube, and allowed to propagate for approximately 20 cm. The required energy for successful initiation of the detonation will be measured for different cases of grid resolution and numerical schemes. A modified version of the TVD scheme has also been devised, which allows for much lower numerical diffusion of the radical species in the exponentially growing region behind the shock. The effect of this modification will be demonstrated. Oscillations in peak pressure and induction length are seen to develop in some cases: the oscillations can go through a sequence of modes, from a regular, high frequency mode to a low frequency mode with period doubling. A chaotic regime can also be obtained. General conclusions on the quality of algorithms will be presented. We will also discuss the performance of a version of the code developed on the IBM SP2 parallel computer.

Cambier, Jean-Luc↗

Extension and Validation of the Pendulum Model for Longitudinal Solar Prominence Oscillations

Context. Longitudinal oscillations in prominences are common phenomena on the Sun. These oscillations can be used to infer the geometry and intensity of the filament magnetic field. Previous theoretical studies of longitudinal oscillations made two simplifying assumptions: uniform gravity and semicircular dips on the supporting flux tubes. However, the gravity is not uniform and realistic dips are not semicircular. Aims. Our aim is to understand the effects of including the nonuniform solar gravity on longitudinal oscillations and explore the validity of the pendulum model with different flux-tube geometries. Methods. We first derived the equation describing the motion of the plasma along the flux tube including the effects of nonuniform gravity, yielding corrections to the original pendulum model. We also computed the full numerical solutions for the normal modes and compared them with the new pendulum approximation. Results. We find that the nonuniform gravity introduces a significant modification in the pendulum model. We also found a cut-off period; i.e., the longitudinal oscillations cannot have a period longer than 167 min. In addition, considering different tube geometries, the period depends almost exclusively on the radius of curvature at the bottom of the dip. Conclusions. We conclude that nonuniform gravity significantly modifies the pendulum model. These corrections are important for prominence seismology, because the inferred values of the radius of curvature and minimum magnetic-field strength differ substantially from those of the old model. However, we find that the corrected pendulum model is quite robust and is still valid for noncircular dips.

Sun: corona↗