Search NASA⌕ Search

SEARCH · Search NASA

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Numerical simulation of boundary layers. Part 3: Turbulence and relaminarization in sink flows

Direct numerical simulations of sink-flow boundary layers, with acceleration parameters K between 1.5 and 3.0 x 10 to the -6 power, are presented. The three-dimensional, time-dependent Navier-Stokes equations are solved numerically using a spectral method, with about 10 to the -6 degrees of freedom. The flow is assumed to be statistically steady, and self-similar. A multiple-scale approximation and periodic conditions are applied to the fluctuations. The turbulence is studied using instantaneous and statistical results. Good agreement with the experiments of Jones and Launder is observed. Two effects of the favorable pressure gradient are to extend the logarithmic layer, and to alter the energy balance of the turbulence near the edge of the boundary layer. At very low Reynolds number the logarithmic layer is shortened and slightly displaced, but wall-layer streaks are present even at the lowest values of the Reynolds number for which turbulence can be sustained. Large quiescent patches appear in the flow. Relaminarization occurs at K = 3.0 x 10 to the -6 power, corresponding to a Reynolds number, based on momentum thickness, of about 330.

Spalart, Philippe R.↗
NASA NTRS1988

Semi-discrete approximations to nonlinear systems of conservation laws; consistency and L(infinity)-stability imply convergence

A convergence theory for semi-discrete approximations to nonlinear systems of conservation laws is developed. It is shown, by a series of scalar counter-examples, that consistency with the conservation law alone does not guarantee convergence. Instead, a notion of consistency which takes into account both the conservation law and its augmenting entropy condition is introduced. In this context it is concluded that consistency and L(infinity)-stability guarantee for a relevant class of admissible entropy functions, that their entropy production rate belongs to a compact subset of H(loc)sup -1 (x,t). One can now use compensated compactness arguments in order to turn this conclusion into a convergence proof. The current state of the art for these arguments includes the scalar and a wide class of 2 x 2 systems of conservation laws. The general framework of the vanishing viscosity method is studied as an effective way to meet the consistency and L(infinity)-stability requirements. How this method is utilized to enforce consistency and stability for scalar conservation laws is shown. In this context we prove, under the appropriate assumptions, the convergence of finite difference approximations (e.g., the high resolution TVD and UNO methods), finite element approximations (e.g., the Streamline-Diffusion methods) and spectral and pseudospectral approximations (e.g., the Spectral Viscosity methods).

Tadmor, Eitan↗
NASA NTRS1988

Navier-Stokes simulations of axisymmetric vortex rings

The solution of the incompressible Navier-Stokes equations by means of a spectral method is used in the present numerical investigation of viscous axisymmetric vortex rings, in the cases of a single ring over the 0.001-1000 Re number range, and a pair of interacting rings. The propagation speeds of vortex rings of varying Re are computed and compared with the results of Saffman's (1970) theory; the present error estimate is found to be smaller than Saffman's, and to decrease with increasing Re at fixed core/ring radius ratio. The 'leapfrogging' of two vortex rings of equal sign and Re=1000 each is observed, indicating severe strain of the vortex core through the first pass and merging during the second pass.

Stanaway, S. K.↗
NASA NTRS1988

Direct simulation of a turbulent boundary layer up to R sub theta = 1410

The turbulent boundary layer on a flat plate, with zero pressure gradient, is simulated numerically at four stations between R sub theta = 225 and R sub theta = 1410. The three-dimensional time-dependent Navier-Stokes equations are solved using a spectral method with up to about 10 to the 7th grid points. Periodic spanwise and streamwise conditions are applied, and a multiple-scale procedure is applied to approximate the slow streamwise growth of the boundary layer. The flow is studied, primarily, from a statistical point of view. The solutions are compared with experimental results. The scaling of the mean and turbulent quantities with Reynolds number is compared with accepted laws, and the significant deviations are documented. The turbulence at the highest Reynolds number is studied in detail. The spectra are compared with various theoretical models. Reynolds-stress budget data are provided for turbulence-model testing.

Spalart, Philippe R.↗
NASA NTRS1988

Turbulent magnetohydrodynamic density fluctuations

A spectral-method numerical code is used to compute mass-density fluctuation spectra in turbulent magnetofluids. The computations are used to test and extend the analytical theory of density variations in slightly compressible magnetofluids given by Montgomery, et al. (1987) and used to infer inertial-range density-fluctuation spectra for the nearby interstellar medium and solar wind. A local equation of state is assumed, relating density to pressure. Constant, scalar resistivities and viscosities are used. In the limit of low Mach numbers and high mechanical-to-magnetic pressure ratios, the fit of the computations to the analytical theory is seen to be close.

Shebalin, John V.↗
NASA NTRS1988

Single-grid spectral collocation for the Navier-Stokes equations

The aim of the paper is to study a collocation spectral method to approximate the Navier-Stokes equations: only one grid is used, which is built from the nodes of a Gauss-Lobatto quadrature formula, either of Legendre or of Chebyshev type. The convergence is proven for the Stokes problem provided with inhomogeneous Dirichlet conditions, then thoroughly analyzed for the Navier-Stokes equations. The practical implementation algorithm is presented, together with numerical results.

Bernardi, Christine↗
NASA NTRS1989

Numerical simulation of Goertler/Tollmien-Schlichting wave-interaction

The problem of nonlinear development of Goertler vortices and interaction with Tollmien-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 NTRS1988

Cut-and-connect of two antiparallel vortex tubes

Motivated by an early conjecture that vortex cut-and-connect plays a key role in mixing and production of turbulence, helicity and aerodynamic noise, the cross-linking of two antiparallel viscous vortex tubes via direct numerical simulation is studied. The Navier-Stokes equations are solved by a dealiased pseudo-spectral method with 64 cubed grid points in a periodic domain for initial Reynolds numbers Re up to 1000. The vortex tubes are given an initial sinusoidal perturbation to induce a collision and keep the two tubes pressed against each other as annihilation continues. Cross-sectional and wire plots of various properties depict three stages of evolution: (1) Inviscid induction causing vortex cores to first approach and form a contact zone with a dipole cross-section, and then to flatten and stretch; (2) Vorticity annihilation in the contact zone accompanied by bridging between the two vortices at both ends of the contact zone due to a collection of cross-linked vortex lines, now orthogonal to the initial vortex tubes. The direction of dipole advection in the contact zone reverses; and (3) Threading of the remnants of the original vortices in between the bridges as they pull apart. The crucial stage 2 is shown to be a simple consequence of vorticity annihilation in the contact zone, link-up of the un-annihilated parts of vortex lines, and stretching and advection by the vortex tube swirl of the cross-linked lines, which accumulate at stagnation points in front of the annihilating vortex dipole. It is claimed that bridging is the essence of any vorticity cross-linking and that annihilation is sustained by stretching of the dipole by the bridges. Vortex reconnection details are found to be insensitive to asymmetry. Modeling of the reconnection process is briefly examined. The 3D spatial details of scalar transport (at unity Schmidt number), enstrophy production, dissipation and helicity are also examined.

Melander, Mogens V.↗
NASA NTRS1989

Absolute and convective instability of a cylinder wake

The stability behavior of a circular cross-section cylinder's wake at Reynolds number values of up to 45 is presently investigated by means of local linear stability theory. The steady-wake profiles computed are Navier-Stokes solutions of a uniform, incompressible viscous flow around a cylinder obtained by a spectral method. An absolutely unstable region is found to begin to form at a Reynolds number of about 20, and grows with incresing Reynolds number. The onset of global instability response must be characterized by a critical length of an absolutely unstable region; a critical Reynolds number criterion and preferred frequency are proposed based on linear stability analysis.

Yang, X.↗
NASA NTRS1990

Unsteady three-dimensional marginal separation, including breakdown

A situation involving a three-dimensional marginal separation is considered, where a (steady) boundary layer flow is on the verge of separating at a point (located along a line of symmetry/centerline). At this point, a triple-deck is included, thereby permitting a small amount of interaction to occur. Unsteadiness is included within this interaction region through some external means. It is shown that the problem reduces to the solution of a nonlinear, unsteady, partial-integro system, which is solved numerically by means of time-marching together with a pseudo-spectral method spatially. A number of solutions to this system are presented which strongly suggest a breakdown of this system may occur, at a finite spatial position, at a finite time. The structure and details of this breakdown are then described.

Duck, Peter W.↗
NASA NTRS1989

Multidomain spectral solution of shock-turbulence interactions

The use of a fitted-shock multidomain spectral method for solving the time-dependent Euler equations of gasdynamics is described. The multidomain method allows short spatial scale features near the shock to be resolved throughout the calculation. Examples presented are of a shock-plane wave, shock-hot spot and shock-vortex street interaction.

Kopriva, David A.↗
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↗