Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Wave wave interactions”

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 667 records · Page 37

Simulation of unsteady flows through stator and rotor blades of a gas turbine using the Chimera method

A two-dimensional model to solve compressible Navier-Stokes equations for the flow through stator and rotor blades of a turbine is developed. The flow domains for the stator and rotor blades are coupled by the Chimera method that makes grid generation easy and enhances accuracy because the area of the grid that have high turning of grid lines or high skewness can be eliminated from the computational domain after the grids are generated. The results of flow computations show various important features of unsteady flows including the acoustic waves interacting with boundary layers, Karman vortex shedding from the trailing edge of the stator blades, pulsating incoming flow to a rotor blade from passing stator blades, and flow separation from both suction and pressure sides of the rotor blades.

Nakamura, S.↗

Numerical simulation of receptivity and transition in a boundary layer on a flat plate with a suction hole

A Navier-Stokes solver has been modified to allow for unsteady boundary conditions to be imposed at the inflow boundary. Because the streamwise derivatives are discretized with three-point compact differences, there is no need for special treatment of the points adjacent to the inflow boundary, provided that all necessary derivatives are properly specified at the boundary points. The unsteady boundary conditions can be obtained from a previous computation, from a solution of the parabolized stability equations (PSE), or from an analytical exact or approximate solution of the Navier-Stokes equations, e.g. an asymptotic expansion solution to model a sound wave interacting with the boundary layer.

Fasel, Hermann F.↗

Numerical simulation of shock-induced combustion past blunt bodies using shock-fitting technique

Two-dimensional axisymmetric, reacting viscous flow over blunt projectiles is computed to study shock-induced combustion at Mach 5.11 and Mach 6.46 in hydrogen-air mixture. A finite-difference, shock-fitting method is used to solve the complete set of Navier-Stokes and species conservation equations. In this approach, the bow shock represents a boundary of the computational domain and is treated as a discontinuity across which Rankine-Hugoniot conditions are applied. All interior details of the flow such as compression waves, reaction front, and the wall boundary layer are captured automatically in the solution. Since shock-fitting approach reduces the amount of artificial dissipation, all the intricate details of the flow are captured much more clearly than has been possible with the shock-capturing approach. This has allowed an improved understanding of the physics of shock-induced combustion over blunt projectiles and the numerical results can now be explained more readily with one-dimensional wave-interaction model than before.

Ahuja, J. K.↗

Parametric study of shock-induced combustion in a hydrogen air system

A numerical parametric study is conducted to simulate shock-induced combustion under various free-stream conditions and varying blunt body diameter. A steady combustion front is established if the free-stream Mach number is above the Chapman-Jouguet speed of the mixture, whereas an unsteady reaction front is established if the free-stream Mach number is below or at the Chapman-Jouguet speed of the mixture. The above two cases have been simulated for Mach 5.11 and Mach 6.46 with a projectile diameter of 15 mm. Mach 5.11, which is an underdriven case, shows an unsteady reaction front, whereas Mach 6.46, which is an overdriven case, shows a steady reaction front. Next for Mach 5. 11 reducing the diameter to 2.5 mm causes the instabilities to disappear, whereas, for Mach 6.46 increasing the diameter of the projectile to 225 mm causes the instabilities to reappear, indicating that Chapman-Jouguet speed is not the only deciding factor for these instabilities to trigger. The other key parameters are the projectile diameter, induction time, activation energy and the heat release. The appearance and disappearance of the instabilities have been explained by the one-dimensional wave interaction model.

Ahuja, J. K.↗

Response of the thermosphere and ionosphere to geomagnetic storms

Four numerical simulations have been performed, at equinox, using a coupled thermosphere-ionosphere model, to illustrate the response of the upper atmosphere to geomagnetic storms. The storms are characterized by an increase in magnetospheric energy input at high latitude for a 12-hour period; each storm commences at a different universal time (UT). The initial response at high latitude is that Joule heating raises the temperature of the upper thermosphere and ion drag drives high-velocity neutral winds. The heat source drives a global wind surge, from both polar regions, which propagates to low latitudes and into the opposite hemisphere. The surge has the character of a large-scale gravity wave with a phase speed of about 600 m/s. Behind the surge a global circulation of magnitude 100 m/s is established at middle latitudes, indicating that the wave and the onset of global circulation are manifestations of the same phenomena. A dominant feature of the response is the penetration of the surge into the opposite hemisphere where it drives poleward winds for a few hours. The global wind surge has a preference for the night sector and for the longitude of the magnetic pole and therefore depends on the UT start time of the storm. A second phase of the meridional circulation develops after the wave interaction but is also restricted, in this case by the buildup of zonal winds via the Coriolis interaction. Conservation of angular momentum may limit the buildup of zonal wind in extreme cases. The divergent wind field drives upwelling and composition change on both height and pressure surfaces. The composition bulge responds to both the background and the storm-induced horizontal winds; it does not simply rotate with Earth. During the storm the disturbance wind modulates the location of the bulge; during the recovery the background winds induce a diurnal variation in its position. Equatorward winds in sunlight produce positive ionospheric changes during the main driving phase of the storm. Negative ionospheric phases are caused by increases of molecular nitrogen in regions of sunlight, the strength of which depends on longitude and the local time of the sector during the storm input. Regions of positive phase in the ionosphere persist in the recovery period due to decreases in mean molecular mass in regions of previous downwelling. Ion density changes, expressed as a ratio of disturbed to quiet values, exhibit a diurnal variation that is driven by the location of the composition bulge; this variation explains the ac component of the local time variation of the observed negative storm phase.

Fuller-Rowell, T. J.↗

Damping of orbital inclinations by bending waves

An inclined secondary orbiting in a disk will launch bending waves from resonance sites where the Doppler shifted forcing frequency matches the disk's natural frequency for vertical oscillations. These vertical resonances are of two types: external resonances falling interior and exterior to the perturber's semimajor axis that excite its inclination and coorbiting resonances that fall at the perturber's orbit and damp its inclination. We show that torques from coorbiting resonances dominate the bending wave interaction for a constant density disk. In this case the inclination ultimately decay and an estimate of the characteristic time scale for this process is made.

Ward, William R.↗

Efficient implementation of weighted ENO schemes

In this paper, we further analyze, test, modify and improve the high order WENO (weighted essentially non-oscillatory) finite difference schemes of Liu, Osher and Chan. It was shown by Liu et al. that WENO schemes constructed from the r-th order (in L1 norm) ENO schemes are (r+1)-th order accurate. We propose a new way of measuring the smoothness of a numerical solution, emulating the idea of minimizing the total variation of the approximation, which results in a 5-th order WENO scheme for the case r = 3, instead of the 4-th order with the original smoothness measurement by Liu et al. This 5-th order WENO scheme is as fast as the 4-th order WENO scheme of Liu et al., and both schemes are about twice as fast as the 4-th order ENO schemes on vector supercomputers and as fast on serial and parallel computers. For Euler systems of gas dynamics, we suggest computing the weights from pressure and entropy instead of the characteristic values to simplify the costly characteristic procedure. The resulting WENO schemes are about twice as fast as the WENO schemes using the characteristic decompositions to compute weights, and work well for problems which do not contain strong shocks or strong reflected waves. We also prove that, for conservation laws with smooth solutions, all WENO schemes are convergent. Many numerical tests, including the 1D steady state nozzle flow problem and 2D shock entropy wave interaction problem, are presented to demonstrate the remarkable capability of the WENO schemes, especially the WENO scheme using the new smoothness measurement, in resolving complicated shock and flow structures. We have also applied Yang's artificial compression method to the WENO schemes to sharpen contact discontinuities.

Jiang, Guang-Shan↗

Investigation of shock-induced combustion past blunt projectiles

A numerical study is conducted to simulate shock-induced combustion in premixed hydrogen-air mixtures at various free-stream conditions and parameters. Two-dimensional axisymmetric, reacting viscous flow over blunt projectiles is computed to study shock-induced combustion at Mach 5.11 and Mach 6.46 in hydrogen-air mixture. A seven-species, seven reactions finite rate hydrogen-air chemical reaction mechanism is used combined with a finite-difference, shock-fitting method to solve the complete set of Navier-Stokes and species conservation equations. The study has allowed an improved understanding of the physics of shock-induced combustion over blunt projectiles and the numerical results can now be explained more readily with one-dimensional wave-interaction model.

Ahuja, J. K.↗

Numerical Simulation of Shock-Induced Combustion Past Blunt Bodies Using Shock-Fitting Technique

Two-dimensional axisymmetric, reacting viscous flow over blunt projectiles is computed to study shock induced combustion at Mach 5.11 and Mach 6.46 in hydrogen-air mixture. A finite-difference, shock-fitting method is used to solve the complete set of Navier Stokes and species conservation equations. In this approach, the bow shock represents a boundary of the computational domain and is treated as a discontinuity across which Rankine-Hugoniot conditions are applied. All interior details of the flow such as compression waves, reaction front, and the wall boundary layer are captured automatically in the solution. Since shock-fitting approach reduces the amount of artificial dissipation, all the intricate details of the flow are captured much more clearly than has been possible with the shock-capturing approach. This has allowed an improved understanding of the physics of shock-induced combustion over blunt projectiles and the numerical results can now be explained more readily with one dimensional wave-interaction model than before.

Ahuja, J. K.↗

Coupling Processes Between Atmospheric Chemistry and Climate

This is the first semi-annual report for NAS5-97039 summarizing work performed for January 1997 through June 1997. Work in this project is related to NAS1-20666, also funded by NASA ACMAP. The work funded in this project also benefits from work at AER associated with the AER three-dimensional isentropic transport model funded by NASA AEAP and the AER two-dimensional climate-chemistry model (co-funded by Department of Energy). The overall objective of this project is to improve the understanding of coupling processes between atmospheric chemistry and climate. Model predictions of the future distributions of trace gases in the atmosphere constitute an important component of the input necessary for quantitative assessments of global change. We will concentrate on the changes in ozone and stratospheric sulfate aerosol, with emphasis on how ozone in the lower stratosphere would respond to natural or anthropogenic changes. The key modeling tools for this work are the AER two-dimensional chemistry-transport model, the AER two-dimensional stratospheric sulfate model, and the AER three-wave interactive model with full chemistry.

Ko, M. K. W.↗

Design of Three-Dimensional Hypersonic Inlets with Rectangular to Elliptical Shape Transition

A methodology has been devised for the design of three-dimensional hypersonic inlets which include a rectangular to elliptical shape transition. This methodology makes extensive use of inviscid streamtracing techniques to generate a smooth shape transition from a rectangular-like capture to an elliptical throat. Highly swept leading edges and a significantly notched cowl enable use of these inlets in fixed geometry configurations. The design procedure includes a three dimensional displacement thickness calculation and uses established correlations to check for boundary layer separation due to shock wave interactions. Complete details of the design procedure are presented and the characteristics of a modular inlet with rectangular to elliptical shape transition and a design point of Mach 7.1 are examined. Comparison with a classical two-dimensional inlet optimized for maximum total pressure recovery indicates that this three-dimensional inlet demonstrates good performance even well below its design point.

Smart, M. K.↗

Oscillator Strengths for Fine-Structure Transitions in S III

Oscillator strengths and transition probabilities for transitions among the fine-structure levels of the terms belonging to the 3s(sup 2)3p(sup 2), 3s3p(sup 3), 3s(sup 2)3p3d, 3s(sup 2)3p4s, 3s(sup 2)3p4p, and 3s(sup 2)3p4d configurations of S III are calculated using extensive configuration-interaction wave functions. The relativistic effects in intermediate coupling are incorporated by means of the Breit-Pauli Hamiltonian. Small adjustments to the diagonal elements of the Hamiltonian matrices have been made so that the energy splittings are as close as possible to the experimental values. The present results are compared with other available calculations and experiments.

Tayal, S. S.↗

Oscillator Strengths of Allowed and Intercombination Transitions in Neutral Sulfur

We have calculated oscillator strengths and transition probabilities of electric-dipole allowed and intercombination transitions from fine-structure levels of the ground 3s(sup 2)3p(sup 4) configuration to the levels belonging to configurations 3s(sup 2)3p(sup 3)4s, 3s(sup 2) 3p(sup 3)5s, 3(sup 2)3p(sup 3)3d, 3s(sup 2)3p(sup 3)4d of neutral sulfur. Extensive configuration-interaction wave functions are used to represent these levels. The relativistic corrections have been included through the Breit-Pauli Hamiltonian. The results are compared with previous theoretical calculations and with measurements.

Tayal, S. S.↗

Collision Strengths for Electron Collisional Excitation of S II

Electron collisional excitation strengths for inelastic transitions in S II are calculated using the R-matrix method in a 19-state (3s(sup 2)3p(sup 3)(sup 4)S(sup o), (sup 2)D(sup o), (sup 2)p(sup o), 3s3p(sup 4)(sup 4)P, (sup 2)D, (sup 2)S, 3S(sup 2)3p(sup 2)3d(sup 2)P, (sup 4)F, (sup 4)D, (sup 2)F, (sup 4)P, 3s(sup 2)3p(sup 2)4s(sup 4)P, (sup 2)P, 3s(sup 2)3p(sup 2)4p(sup 2)S(s o), (sup 4)D(sup o), (sup 4)P(sup o), (sup 2)D(sup o), (sup 4)S(sup o), (sup 2)P(sup o)) close-coupling approximation. These target states are represented by extensive configuration-interaction wave functions that give excitation energies and oscillator strengths that are usually in good agreement with the experimental values and the available accurate calculations. The present results for collision strengths are in very good agreement with the recent merged beams energy loss measurement of Liao et al. and agree reasonably well with the 18-state R-matrix calculation of Ramsbottom, Bell, & Stafford, but show significant differences from the 12-state R-matrix calculation of Cai & Pradhan.

Tayal, S. S.↗

Coupling Processes between Atmospheric Chemistry and Climate

This is the fourth semi-annual report for NAS5-97039, covering the time period July through December 1998. The overall objective of this project is to improve the understanding of coupling processes between atmospheric chemistry and climate. Model predictions of the future distributions of trace gases in the atmosphere constitute an important component of the input necessary for quantitative assessments of global change. We will concentrate on the changes in ozone and stratospheric sulfate aerosol, with emphasis on how ozone in the lower stratosphere would respond to natural or anthropogenic changes. The key modeling tools for this work are the Atmospheric and Environmental Research (AER) two-dimensional chemistry-transport model, the AER two-dimensional stratospheric sulfate model, and the AER three-wave interactive model with full chemistry. For this six month period, we report on a modeling study of new rate constant which modify the NOx/NOy ratio in the lower stratosphere; sensitivity to changes in stratospheric water vapor in the future atmosphere; a study of N2O and CH4 observations which has allowed us to adjust diffusion in the 2-D CTM in order to obtain appropriate polar vortex isolation; a study of SF6 and age of air with comparisons of models and measurements; and a report on the Models and Measurements II effort.

Ko, Malcolm K. W.↗

Effective Collision Strengths for Electron Impact Excitation of Inelastic Transitions in S III

We have calculated electron collisional excitation strengths for all electric dipole forbidden, semi-forbidden, and allowed transitions among the lowest 17 LS states 3s(exp 2)3p(exp 2) P-3, D-1, S-1, 3s3p(exp 3)S-5(exp 0), D-3(exp 0), P-3(exp 0), P-1(exp 0), S-3(exp 0), D-1(exp 0), 3S(exp 2)3p3d D-1(exp 0), F-3(exp 0), P-3(exp 0), D-3(exp 0), F-3(exp 0), P-1(exp 0), and 3S(exp 2)3p4S P-3(exp 0), P-l(exp 0) of S III using the R-matrix method. These S m states are represented by fairly extensive configuration-interaction wave functions that yield excited state energies in close agreement with recent laboratory measurements. Rydberg series of resonances converging to the excited state thresholds are explicitly included in the scattering calculation. The effective collision strengths are determined assuming Maxwellian distribution of electron energies. These are listed over a wide temperature range ([0.5-10] x 10(exp 4) K) and compared, where possible, with other available calculations. Subject headings: atomic data - atomic processes

Tayal, S. S.↗

Spurious Numerical Oscillations in Numerical Simulation of Supersonic Flows Using Shock Capturing Schemes

The numerical simulation of transitional and turbulent processes in hypersonic boundary layers often involves a physical process of a shock-disturbance wave interaction in complex two-dimensional and three-dimensional flow fields. For such simulations, it is required that there be a high order of accuracy in capturing both the shock and the small disturbances. The purpose of this paper is to evaluate the viability of using high order shock capturing schemes to track small disturbances in a multi-dimensional steady hypersonic flow. The numerical methods that are to be studied are the Total Variation Diminishing (TVD) scheme, and Essentially Non-Oscillatory (ENO) scheme. This paper shows that the presence of numerical oscillations in the flow field solution may drastically hinder any attempt at tracking the propagation of any physical disturbances. It has been found that the numerical oscillations that exist for shock capturing methods may be significant enough to pollute a flow field containing small physical disturbances. The effects of the refinement of the grid do not reduce the oscillations, but rather they decrease the wavelength of the oscillations. It is shown that by aligning the shock with the grid, the amplitude of these spurious oscillations may be greatly reduced.

Lee, Theodore K.↗

Kinetic Modeling of Plasmaspheric Refilling

We have a new model of ion transport that we have applied to the problem of plasmaspheric flux tube refilling after a geomagnetic disturbance. Effects included are a time-varying ionospheric source, self-consistent Coulomb collisions, field-aligned electric field, and ion cyclotron wave heating. We see refilling rates similar to those of earlier observations and models, except when the electric field is included. In this case, the refilling rates are larger than previously predicted, yet still within the experimental limits. Wave interactions increase the trapped zone population and thus the density, especially during the few few hours of the refilling process. Results are compared with observations from DE/RIMS and Polar/TIDE, as well as hydrodynamic and kinetic particle model simulations.

Liemohn, M. W.↗