Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spectral methods”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

A numerically exact full wave packet approach to molecule-surface scattering

A numerically exact spectral method for solving the time-dependent Schroedinger equation in spherical coordinates is described. The angular dependence of the wave function is represented on a two-dimensional grid of evenly spaced points. The fast Fourier transform algorithm is used to transform between the angle space representation of the wave function and its conjugate representation in momentum space. The time propagation of the wave function is evaluated using an expansion of the time evolution operator as a series of Chebyshev polynomials. Calculations performed for a model system representing H2 scattering from a rectangular corrugated surface yield transition probabilities that are in excellent agreement with those obtained using the close-coupling wave packet (CCWP) method. However, the new method is found to require substantially more computation time than the CCWP method because of the large number of grid points needed to represent the angular dependence of the wave function and the variation in the number of terms required in the Chebyshev representation of the time evolution operator.

Mowrey, R. C.↗
NASA NTRS1989

A numerical study of the interaction between unsteady free-stream disturbances and localized variations in surface geometry

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small free-stream disturbance and a small localized variation of the surface geometry has been carried out using both finite-difference and spectral methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T-S waves generated by the interaction between the free-stream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T-S waves.

Bodonyi, R. J.↗
NASA NTRS1990

Numerical simulation of interactions between Goertler vortices and Tollmien-Schlichting waves

The problem of nonlinear development of Goertler vortices and interaction with Tollmein-Schlichting (TS) waves is considered within the framework of incompressible Navier-Stokes equations which are solved by a Fourier-Chebyshev spectral method. It is shown that two-dimensional waves can be excited in the flow modulated by Goertler vortices. Due to nonlinear effects, this interaction further leads to the development of oblique waves with spanwise wavelength equal to the Goertler vortex wavelength. Interaction is also considered of oblique waves with spanwise wavelength twice that of Goertler vortices.

Malik, M. R.↗
NASA NTRS1989

Understanding transition and turbulence through direct simulations

Direct simulations consist in solving the full Navier-Stokes equations, without any turbulence model, and describing all the detailed features of the flow. Usually the flows are three-dimensional and time-dependent and contain both coarse and fine structures, which makes the numerical task very challenging in terms of both the algorithm and the computational effort. Most of the work until now has involved spectral methods, which are highly accurate but not very flexible in terms of geometry or complex equations. For that reason, future work will also rely on high-order finite-difference or other methods. Direct simulations complement experimental work, and both contribute to the theory and the empirical knowledge of turbulence. Once such a simulation has been shown to be accurate, the flow field is completely known in three dimensions and time, including the pressure, the vorticity and any other quantity. On the other hand, most simulations to date solved the incompressible equations in rather simple geometries, and direct simulations will always be limited to moderate Reynolds numbers. Extensive simulations have been conducted in homogeneous turbulence, channel flows, boundary layers, and mixing layers. Much effort is devoted to addressing flows with compressibility and chemical reactions, and to new geometries such as a backward-facing step.

Spalart, P. R.↗
NASA NTRS1990

Analysis of spectral projectors in one-dimensional domains

A class of projection operators with values in a subspace of polynomials is analyzed. These projection operators are related to the Hilbert spaces involved in the numerical analysis of spectral methods. They are, in the first part of the paper, the standard Sobolev spaces and, in the second part, some weighted Sobolev spaces, the weight of which is related to the orthoronality relation satisfied by the Chebyshev polynomials. These results are used to study the approximation of a model fourth-order problem.

Maday, Y.↗
NASA NTRS1991

A drag measurement technique for free piston shock tunnels

A new technique is described for measuring drag with 100-microsecond rise time on a nonlifting model in a free piston shock tunnel. The technique involves interpretation of the stress waves propagating within the model and its support. A finite element representation and spectral methods are used to obtain a mean square optimal estimate of the time history of the aerodynamic loading. Thus, drag is measured instantaneously and the previous restriction caused by the mechanical time constant of balances is overcome. The effectiveness of the balance is demonstrated by measuring the drag on cones with 5 and 15 deg semi-vertex angles in nominally Mach 5.6 flow with stagnation enthalpies from 2.6 to 33 MJ/kg.

Sanderson, S. R.↗
NASA NTRS1990

Enstrophy budget in decaying compressible turbulence

The interaction of vorticity and shock waves in decaying compressible turbulence is investigated by direct numerical simulation of the Navier-Stokes equation using a spectral method with 64-cubed collocation points. The Taylor microscale Reynolds number decreases from 36 to 17 and the rms Mach number from 0.8 to 0.4 during the decay. It is observed that vorticity is created inside curved shocks through baroclinic interaction and intensified by the compression of fluid elements. Vortex stretching is, on the other hand, dominant outside shocks and is responsible for enhancing the strength of vorticity. A dipole of vorticity is typically created at the vertex of curved shocks, the structure of which is explained in terms of jump relations of velocity derivatives across a curved shock.

Kida, Shigeo↗
NASA NTRS1991

Studies on the dynamics of Venus ionosphere - Effects of varying the ionopause height and the onset of turbulence

The conservation equations of plasma dynamics in the upper ionosphere of Venus have been solved by using a spectral method in the horizontal and finite differencing in the vertical direction. The effect of varying the ionopause height on the computed nightside ion densities is investigated. These ion densities show a sharp decrease as the ionopause altitude is reduced to 300 km. The effect of viscous forces on the horizontal plasma flow is investigated for a wide range of values of the coefficient of viscosity. The Reynolds numbers characteristics of the flow are calculated and the conditions for the onset of turbulence discussed. It is found that the Reynolds number can be large (greater than 1000) in the subsolar region for a coefficient of viscosity of up to 1.6 x 10 to the -10th g/cm s. The influence of magnetic fields on viscosity is also discussed.

Singhal, R. P.↗
NASA NTRS1992

Spectral solution of inviscid supersonic flows over wedges and axisymmetric cones

A shock-fitted multidomain spectral collocation method is used to solve both steady and unsteady inviscid supersonic flows over bodies. New aspects of the method include two subdomain interface types and a zonal solution procedure to get efficient convergence to steady-state in supersonic regions. Steady-state examples include flow over a sharp cone, a hyperbolic cone and a hyperbolic wedge. For the cone, the exact flow solution is used to show that the method is spectrally accurate. As an example of an unsteady flow, a calculation of the interaction of a free-stream hot ring with the flow over a sharp cone is presented.

Kopriva, David A.↗
NASA NTRS1991

Towards large eddy and direct simulation of complex turbulent flows

Recent advances in the methodology for direct numerical simulation of turbulent flows and some of the current applications are reviewed. It is argued that high-order finite difference schemes yield solutions with comparable accuracy to the spectral methods with the same number of degrees of freedom. The effects of random inflow conditions on the downstream evolution of turbulence are discussed.

Moin, Parviz↗
NASA NTRS1992

The eigenvalue spectrum of the Rayleigh equation for a plane shear layer

The eigenvalue spectrum of the Rayleigh equation is examined using three different solution techniques. In particular, a simple second-order finite difference scheme and two spectral methods, the Chebyshev tau and Chebyshev collocation methods, are used to discretize the equation. All of the approximation methods are shown to be capable of predicting the discrete spectrum as well as the continuous spectrum associated with the critical point singularity for the Rayleigh equation. The global eigenvalue methods considered here provide an efficient way of obtaining either an approximation to the complete eigenvalue spectrum or initial guesses for a local shooting procedure for the discrete part of the spectrum.

Liou, William W.-W.↗
NASA NTRS1993

Spectral solution of the viscous blunt-body problem

The viscous blunt-body problem is solved with a shock-fitted Chebyshev spectral method. No explicit artificial viscosity or filtering is needed to obtain smooth, converged solutions. The method is applied to two problems. First, results for the flow over a right circular cylinder in the Mach number range of 5.5-6.0 are compared with experimental data. Second, a solution for a Mach 25 flow over a hyperbolic cone is compared with a viscous shock-layer calculation.

Kopriva, David A.↗
NASA NTRS1993

Numerical simulation of laminar-turbulent transition in a spatially-developing flat plate wake

Laminar-turbulent transition of an incompressible flat-plate wake is investigated by direct numerical integration of the Navier-Stokes equations. For the numerical integration, a combination of finite-difference and spectral methods along with an ADI/Crank-Nicolson/Adams-Bashforth time integration scheme is employed. Subject to 2D forcing, the wake exhibited a rapidly-growing fundamental disturbance that quickly saturated. This saturation was due partly to the stabilizing effect of the mean flow distortion. Downstream of the saturation point, disturbance energy was concentrated in the fundamental disturbance, the second harmonic, and the mean flow distortion component. At large amplitude levels, a Karman vortex street formed. Variations in the 2D forcing level did not alter the qualitative behavior of the disturbances. Simulations of 3D breakdown indicates that the presence of large-amplitude, 2D disturbances tends to initially suppress small-amplitude 3D disturbance growth. Following this initial suppression, a resumption of 3D growth is observed that may have been due to a secondary instability mechanism. For high levels of 3D disturbance energy, lambda-shaped vortical structures formed between adjacent Karman vortices.

Dratler, D. I.↗
NASA NTRS1994

Modeling of the coupled magnetospheric and neutral wind dynamos

This report summarizes the progress made in the first year of NASA Grant No. NAGW-3508 entitled 'Modeling of the Coupled Magnetospheric and Neutral Wind Dynamos.' The approach taken has been to impose magnetospheric boundary conditions with either pure voltage or current characteristics and solve the neutral wind dynamo equation under these conditions. The imposed boundary conditions determine whether the neutral wind dynamo will contribute to the high-latitude current system or the electric potential. The semi-annual technical report, dated December 15, 1993, provides further detail describing the scientific and numerical approach of the project. The numerical development has progressed and the dynamo solution for the case when the magnetosphere acts as a voltage source has been evaluated completely using spectral techniques. The simulation provides the field-aligned current distribution at high latitudes due to the neutral wind dynamo. A number of geophysical conditions can be simulated to evaluate the importance of the neutral wind dynamo contribution to the field-aligned current system. On average, field-aligned currents generated by the neutral wind dynamo contributed as much as 30 percent to the large-scale field-aligned current system driven by the magnetosphere. A term analysis of the high-latitude neutral wind dynamo equation describing the field aligned current distribution has also been developed to illustrate the important contributing factors involved in the process. The case describing the neutral dynamo response for a magnetosphere acting as a pure current generator requires the existing spectral code to be extended to a pseudo-spectral method and is currently under development.

Thayer, Jeffrey P.↗
NASA NTRS1994

Calculations of Arctic ozone chemistry using objectively analyzed data in a 3-D CTM

A three-dimensional chemical transport model (CTM) (Kaminski, 1992) has been used to study the evolution of the Arctic ozone during the winter of 1992. The continuity equation has been solved using a spectral method with Rhomboidal 15 (R15) truncation and leap-frog time stepping. Six-hourly meteorological fields from the Canadian Meteorological Center global objective analysis routines run at T79 were degraded to the model resolution. In addition, they were interpolated to the model time grid and were used to drive the model from the surface to 10 mb. In the model, processing of Cl(x) occurred over Arctic latitudes but some of the initial products were still present by mid-January. Also, the large amounts of ClO formed in the model in early January were converted to ClNO3. The results suggest that the model resolution may be insufficient to resolve the details of the Arctic transport during this time period. In particular, the wind field does not move the ClO(x) 'cloud' to the south over Europe as seen in the MLS measurements.

Kaminski, J. W.↗
NASA NTRS1994

Direct numerical simulation of sheared turbulent flow

The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.

Harris, Vascar G.↗
NASA NTRS1994

Nonlinear decay of magnetic helicity in magnetohydrodynamic turbulence with a mean magnetic field

We show that the magnetic helicity associated with fluctuations in homogeneous incompressible magnetohydrodynamic (MHD) turbulence with a mean magnetic field decays in time because of nonlinear processes. Evidence is obtained numerically, by use of both dissipative and nondissipative spectral method simulations. The described effect stands in contrast to expectations based on studies of MHD turbulence without an applied mean field, in which magnetic helicity is transferred nonlinearly to long wavelengths and is preserved in time because of selective decay when dissipation is present. The process of nonlinear decay is described in terms of a generalized ideal helicity invariant and is characterized by the transient production of a mean induced electric field aligned with the applied magnetic field, an effect reminiscent of the alpha effect in dynamo theory. A simple phenomenological model for the decay process is proposed.

Stribling, Troy↗
NASA NTRS1995

Test particle study of minor ions in solar wind turbulence

We perform a parameter study of the temporal evolution of a test particle distribution function in MHD turbulence. The turbulent fields are calculated using a pseudo-spectral method and periodic boundary conditions on a regular grid of 180(exp 3) points, appropriate for incompressible, homogeneous and isotropic turbulence. Initially, the kinetic and the magnetic energy are equal on the average. Both, deterministic and random initial conditions are used, in the former case with zeros of the magnetic field located at grid points, in the latter case located by interpolation between grid points. The evolution of the minor ion distribution function is studied in detail as these turbulent fields evolve, developing strong current and vorticity sheets. Using the full collisionless equation of motion for the test particles, the efficiency of nonlinear interactions can be studied. The results are compared to theoretical predictions and are then discussed in connection with the observations of the dynamical properties of solar wind minor ions derived from in situ observations.

Zurbuchen, Th.↗