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 163 records · Page 9

On the Possibilities of Predicting Geomagnetic Secular Variation with Geodynamo Modeling

We use our MoSST core dynamics model and geomagnetic field at the core-mantle boundary (CMB) continued downward from surface observations to investigate possibilities of geomagnetic data assimilation, so that model results and current geomagnetic observations can be used to predict geomagnetic secular variation in future. As the first attempt, we apply data insertion technique to examine evolution of the model solution that is modified by geomagnetic input. Our study demonstrate that, with a single data insertion, large-scale poloidal magnetic field obtained from subsequent numerical simulation evolves similarly to the observed geomagnetic variation, regardless of the initial choice of the model solution (so long it is a well developed numerical solution). The model solution diverges on the time scales on the order of 60 years, similar to the time scales of the torsional oscillations in the Earth's core. Our numerical test shows that geomagnetic data assimilation is promising with our MoSST model.

Kuang, Wei-Jia↗

Galactic oscillations

Several oscillations have been identified in spherical galaxy models. These are normal mode oscillations in a stable galaxy. Each has its own distinct period and spatial form, and each rings without detectable damping through a Hubble time. The most important are: (1) a simple radial pulsation (fundamental mode), in which all parts of the galaxy move inward or outward with the same phase; and (2) a second spherically symmetrical radial mode with one node, so material inside the node moves outward when material outside moves inward. Numerical experiments suggest that normal mode oscillations may be present in nearly all galaxies at a considerably higher amplitude than has previously been thought. Amplitudes typically run a few percent of equilibrium values, and periods are around 50-300 Myrs in typical galaxies. These time scales are long enough that gas trapped near the center could cool during an oscillation cycle, allowing star formation activity. The second mode oscillations could cause bursts of star formation.

Smith, B. F.↗

High resolution absorption cross-sections and band oscillator strengths of the Schumann-Runge absorption bands of isotopic oxygen, (O-18)2, at 79 K

Cross-sections of (O-18)2 at 79 K have been obtained from photoabsorption measurements at various pressures throughout the wavelength region 177.8-197.8 nm with a 6.65 m photoelectric scanning spectrometer equipped with a 2400 lines/mm grating and having an instrumental width (FWHM) of 0.0013 nm. The measured absorption cross-sections of the Schumann-Runge bands (14,0) through (2,0) are, with the exception of the (12,0) band, independent of the instrumental width. The measured cross-sections are presented graphically here and are available at wavenumber intervals of about 0.1/cm as numerical compilations stored on magnetic tape. Band oscillator strengths of those bands have been determined by direct numerical integration of the measured absolute cross-sections and are in excellent agreement with these theoretically calculated values.

Yoshino, K.↗

Baseline Computational Fluid Dynamics Methodology for Longitudinal-Mode Liquid-Propellant Rocket Combustion Instability

A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Litchford, R. J.↗

Flow oscillations of spike-tipped bodies

The oscillating flow field around a spike-tipped body has been simulated by the numerical solution of the 3-D compressible Navier-Stokes equations at a Mach number of 3.0 and a nominal Reynolds number of 7.87 x 10 to the 6th/m. Computations were performed using a vectorized 3-D Navier-Stokes program on the STAR 100 computer. Numerical solutions confirmed the experimental result that the self-sustained oscillation occurs within a limited range of the protruded spike length to shoulder height ratio. The numerical result predicted correctly the discrete frequency range as well as the rms pressure intensity. The detailed flow structure is also presented and discussed.

Shang, J. S.↗

Numerical solution of the viscous hypersonic flow past blunted cones at angle of attack

Hypersonic viscous flow over spherically blunted cones of large half angle is computed at small angles of attack in the plane of symmetry of the flow field. Time-dependent viscous shock-layer equations in body-oriented coordinates are used to describe the flow field. The shock wave is treated as a discontinuity, across which the Rankine-Hugoniot relations are used to compute the flow conditions behind the shock. A time-marching second-order finite-difference method is used to solve the equations for a perfect gas. The local CFL (Courant-Friedrich-Lewy) time increment is used to advance the solution in time at each grid point. A fourth-order damping is used to damp the oscillations in the flow quantities. The numerical results of the present analysis for quantities such as shock standoff distance, surface-pressure distribution, and heating rates compare well with the existing theoretical and experimental results.

Kumar, A.↗

Optical antenna gain. II - Receiving antennas

Expressions are developed for the gain of a centrally obscured, circular optical antenna used as the collecting and focusing optics in a laser receiver, involving losses due to (1) incoming light blockage by central obscuration, (2) energy spillover at the detector, and (3) the effect of local oscillator distribution in the case of heterodyne or homodyne detection. Numerical results are presented for direct detection and for three types of local oscillator distribution (uniform, Gaussian, and matched).

Degnan, J. J.↗

Dynamics of super cloud clusters, westerly wind bursts, 30-60 day oscillations and ENSO - An unified view

Results are presented on numerical experiments that were carried out to investigate the mechanisms of the observed variabilities in wind and convection associated with supercloud clusters (SCCs), westerly wind bursts, and 30-60 day oscillations in the western Pacific region. It is shown that the generation of a 30-60 day eastward propagating precipitation pattern in the Lau and Peng (1987) model, which can be identified as SCC, is accompanied by convective clusters coming in opposite direction to that of the SCC itself. The results suggest that the westward propagating cloud clusters are produced at the initial stage of the 30-60 day disturbance due to mutual adjustment of the large-scale flow and heating.

Lau, K.-M.↗

Structure and Dynamics of the Quiet Solar Chromosphere

The grant supported research on the structure of the quiet, nonmagnetic chromosphere and on wave excitation and propagation in both the nonmagnetic chromosphere and the magnetic network. The work on the structure of the chromosphere culminated in the recognition that between two competing views of the solar chromosphere, older models by Avrett and collaborators (referred to as VAL) and the newer, dynamical model by Carlsson & Stein (referred to as CS), the clear decision is in favor of the older models, and this in spite of the evident lack of physics, which does not include wave motion and oscillations. The contrast between the static VAL models and the dynamical CS model can be stated most succinctly by comparing the temperature variation implied by the VAL models and the temperature fluctuations of the CS model, which are, respectively, of the order of 10% for the VAL model (at heights where hydrogen is 50% ionized) and a factor of 10 (at the upper boundary of their chromospheric model). The huge fluctuations of the CS model have never been observed, whereas the smaller temperature variations of the VAL models are consistent with ground-based and space-based observations. While it should be obvious which model describes the Sun and which one fails, the case is far from settled in the minds of solar physicists. Thus, much educational work remains to be done and, of course, more research to develop arguments that make the case more convincing. The research on waves and oscillations has been based on a unified theory of excitation of acoustic waves in the field-free atmosphere and of transverse and longitudinal waves in magnetic flux tubes located in the magnetic network by noting, first, that impulsive excitation of all these waves in gravitationally stratified media leads to oscillations at the respective cutoff frequencies and, second, that the observed oscillation frequencies in the nonmagnetic and magnetic parts of the chromosphere match corresponding cutoff frequencies in the upper solar photosphere. The dynamical simulations by Carlsson and Stein have been most instructive and of fundamental importance for understanding wave propagation in a stratified medium by their "flaws", the most important of which is an intensity excess at the H2v emission peak in the H line of Ca II, which is surprising since the observed intensity should have been an upper limit to the simulated intensity. The only plausible explanation for a predicted intensity that is higher than observed is that energy is spread horizontally in upward propagation in the Sun, but not in the plane-wave modeling of CS (as well as by almost everyone else). Investigations of the horizontal size of the region disturbed by the upward-propagating shock in the acoustic-wave propagation implies that the waves in H2v bright-point oscillations emanate from a point source with a diameter corresponding to the width of an intergranular lane, about 100 km, and reach a size of about 4000 km in the upper layers of the chromosphere. Linear, analytic modeling of waves emanating from a point source in a stratified atmosphere shows that the upward-expanding propagation channel does not have sharp boundaries and that the shape of the wave front depends on the order of the wave behind the initial pulse. Otherwise, the behavior of the linear waves resembles that of the nonlinear shock waves observed in the Sun. Research that needs to be done to firm up the conclusions reached above concerns the numerical simulations of nonlinear waves and oscillations in a three-dimensional stratified atmosphere with impulsive excitation, and observations linking directly the horizontal size of the disturbed area in upward propagation to individual waves.

Kalkofen, Wolfgang↗

Thermal relaxation oscillations in horizontal-branch stars

In a recent numerical investigation of the effect of trapped 'cosmions' on horizontal-branch (HB) stars, thermal pulses were observed to occur if the energy transfer provided by these weakly interacting particles was chosen to be efficient enough to break convection. A simple analytic two-zone model of the core of an HB star is presented which allows the interpretation of these pulsations as oscillations of thermal relaxation between the small energy generating 'nuclear core' and the overlying 'gravothermal buffer'. It is shown that any novel form of energy transfer which breaks convection will lead to this behavior where the usual stellar 'conductive equilibrium' is replaced by a limit cycle. The star follows this cycle, oscillating with a period roughly given by the Kelvin-Helmholtz time scale of the core.

Salati, Pierre↗

Model of Semidiurnal Pseudo Tide in the High-Latitude Upper Mesosphere

We present numerical results for the m = 1 meridional winds of semi diurnal oscillations in the high-latitude upper mesosphere, which are generated in the Numerical Spectral Model (NSM) without solar excitations of the tides. Identified with heuristic computer runs, the pseudo tides attain amplitudes that are, at times, as large as the non-migrating tides produced with standard solar forcing. Under the influence of parameterized gravity waves, the nonlinear NSM generates internal oscillations like the quasi-biennial oscillation, that are produced with periods favored by the dynamical properties of the system. The Coriolis force would favor at polar latitudes the excitation of the 12-hour periodicity. This oscillation may help explain the large non-migrating semidiurnal tides that are observed in the region with ground-based and satellite measurements.

Talaat, E. R.↗

Nonlinear evolution of magnetic flux ropes. I - Low-beta limit

We study the nonlinear self-similar evolution of a cylindrical magnetic flux tube with two components of the magnetic field, axial and azimuthal. We restrict ourselves to the case of a plasma of low beta. Introducing a special class of configurations we call 'separable fields', we reduce the problem to an ordinary differential equation. Two cases are to be distinguished: (1) when the total field minimizes on the symmetry axis, the magnetic configuration inexorably collapses, and (2) when, on the other hand, the total field maximizes on the symmetry axis, the magnetic configuration behaves analogously to a nonlinear oscillator. Here we focus on the latter case. The effective potential of the motion contains two terms: a strong repulsive term and a weak restoring term associated with the pinch. We solve the nonlinear differential equation of motion numerically and find that the period of oscillations grows exponentially with the energy of the oscillator. Our treatment emphasizes the role of the force-free configuration as the lowest potential energy state about which the system oscillates.

Osherovich, V. A.↗

Computer study of convection of weakly ionized plasma in a nonuniform magnetic field.

A weakly ionized plasma in a strong and nonuniform magnetic field exhibits an instability analogous to the flute instability in a fully ionized plasma. The instability sets in at a critical magnetic field. To study the final state of the plasma after the onset of the instability, the plasma equations are integrated numerically assuming a certain initial spectrum of small disturbances. In the regime studied, numerical results indicate a final steadily oscillating state consisting of a single finite amplitude mode together with a time-independent modification of the original equilibrium. These results agree with the analytic results obtained by Simon in the slightly supercritical regime. As the magnetic field is increased further, the wavelength of the final oscillation becomes nonunique. There exists a subinterval in the unstable wave band. Final stable oscillation with a wavelength in this subinterval can be established if the initial disturbance has a sufficiently strong component at the particular wavelength.

Shiau, J. N.↗

Unsteady subsonic compressible flow around finite thickness wings.

A general formulation for the unsteady subsonic compressible potential flow around aircraft having arbitrary configurations is presented. An integral representation of the velocity potential is obtained. From this a linear integral equation relating the perturbation potential and its normal derivative (which is known from the boundary conditions) is derived. For the numerical solution of the integral equation, the surface of the aircraft is divided into small elements and the potential is assumed to be constant within each element. Numerical results are obtained for an oscillating finite-thickness wing and indicate good convergence and excellent agreement with existing lifting surface solutions.

Morino, L.↗

Vibration-translation energy transfer in anharmonic diatomic molecules. 2: The vibrational quantum number dependence

A semiclassical model of the inelastic collision between a vibrationally excited anharmonic oscillator and a structureless atom was used to predict the variation of thermally averaged vibration-translation rate coefficients with temperature and initial-state quantum number. Multiple oscillator states were included in a numerical solution for collinear encounters. The results are compared with CO-He experimental values for both ground and excited initial states using several simplified forms of the interaction potential. The numerical model was also used as a basis for evaluating several less complete but analytic models. Two computationally simple analytic approximations were found that successfully reproduced the numerical rate coefficients for a wide range of molecular properties and collision partners. Their limitations were also identified. The relative rates of multiple-quantum transitions from excited states were evaluated for several molecular types.

Mckenzie, R. L.↗

Vibration-translation energy transfer in anharmonic diatomic molecules. II - The vibrational quantum-number dependence

A semiclassical model of the inelastic collision between a vibrationally excited anharmonic oscillator and a structureless atom is used to predict the variation of thermally averaged vibrational-translational rate coefficients with temperature and initial-state quantum number. Multiple oscillator states are included in a numerical solution for collinear encounters. The results are compared with CO-He experimental values for both ground and excited initial states using several simplified forms of the interaction potential. The numerical model is also used as a basis for evaluating several less complete, but analytic, models. Two computationally simple analytic approximations are found that successfully reproduce the numerical rate coefficients for a wide range of molecular properties and collision partners. Their limitations are identified, and the relative rates of multiple-quantum transitions from excited states are evaluated for several molecular types.

Mckenzie, R. L.↗

Vibration of a hydrostatic gas bearing due to supply pressure oscillations

The vibration of a statically loaded, inherently compensated hydrostatic journal bearing due to oscillating supply pressure is investigated. Both angular and radial vibration modes are analyzed. The time-dependent Reynolds equation governing the pressure distribution between the oscillating journal and the sleeve is solved numerically together with the journal equation of motion to obtain the response characteristics of the bearing. The Reynolds equation and the equation of motion are simplified by applying regular perturbation theory for small displacements. The results presented include Bode plots of bearing oscillation gain and phase for a particular bearing configuration for various combinations of parameters over a range of frequencies, including the resonant frequency. The results are compared with the results of an earlier study involving the response of a similar bearing to oscillating exhaust pressure.

Branch, H. D.↗

Finite element or Galerkin type semidiscrete schemes

A finite element of Galerkin type semidiscrete method is proposed for numerical solution of a linear hyperbolic partial differential equation. The question of stability is reduced to the stability of a system of ordinary differential equations for which Dahlquist theory applied. Results of separating the part of numerical solution which causes the spurious oscillation near shock-like response of semidiscrete scheme to a step function initial condition are presented. In general all methods produce such oscillatory overshoots on either side of shocks. This overshoot pathology, which displays a behavior similar to Gibb's phenomena of Fourier series, is explained on the basis of dispersion of separated Fourier components which relies on linearized theory to be satisfactory. Expository results represented.

Durgun, K.↗