Search NASA⌕ Search

Engineering topics

Turkel, E.

Publications and source records attributed to Turkel, E..

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Large Eddy Simulation in the Computation of Jet Noise

Noise can be predicted by solving Full (time-dependent) Compressible Navier-Stokes Equation (FCNSE) with computational domain. The fluctuating near field of the jet produces propagating pressure waves that produce far-field sound. The fluctuating flow field as a function of time is needed in order to calculate sound from first principles. Noise can be predicted by solving the full, time-dependent, compressible Navier-Stokes equations with the computational domain extended to far field - but this is not feasible as indicated above. At high Reynolds number of technological interest turbulence has large range of scales. Direct numerical simulations (DNS) can not capture the small scales of turbulence. The large scales are more efficient than the small scales in radiating sound. The emphasize is thus on calculating sound radiated by large scales.

Mankbadi, R. R.↗
NASA NTRS1998

High-Accuracy Compact MacCormack-Type Schemes for Computational Aeroacoustics

Using MacCormack-type methods, a new class of highly accurate compact MacCormack-type schemes is derived which does not require a tridiagonal matrix inversion to obtain the spatial derivatives. Two examples are shown, and results of these schemes for three linear and nonlinear CAA Benchmark Problems are presented.

Hixon, R.↗
NASA NTRS1998

On Some Numerical Dissipation Schemes

Several schemes for introducing an artificial dissipation into a central difference approximation to the Euler and Navier Stokes equations are considered. The focus of the paper is on the convective upwind and split pressure (CUSP) scheme, which is designed to support single interior point discrete shock waves. This scheme is analyzed and compared in detail with scalar dissipation and matrix dissipation (MATD) schemes. Resolution capability is determined by solving subsonic, transonic, and hypersonic flow problems. A finite-volume discretization and a multistage time-stepping scheme with multigrid are used to compute solutions to the flow equations. Numerical solutions are also compared with either theoretical solutions or experimental data. For transonic airfoil flows the best accuracy on coarse meshes for aerodynamic coefficients is obtained with a simple MATD scheme. The coarse-grid accuracy for the original CUSP scheme is improved by modifying the limiter function used with the scheme, giving comparable accuracy to that obtained with the MATD scheme. The modifications reduce the background dissipation and provide control over the regions where the scheme can become first order.

Swanson, R. C.↗
NASA NTRS1997

Comparison of Several Dissipation Algorithms for Central Difference Schemes

Several algorithms for introducing artificial dissipation into a central difference approximation to the Euler and Navier Stokes equations are considered. The focus of the paper is on the convective upwind and split pressure (CUSP) scheme, which is designed to support single interior point discrete shock waves. This scheme is analyzed and compared in detail with scalar and matrix dissipation (MATD) schemes. Resolution capability is determined by solving subsonic, transonic, and hypersonic flow problems. A finite-volume discretization and a multistage time-stepping scheme with multigrid are used to compute solutions to the flow equations. Numerical results are also compared with either theoretical solutions or experimental data. For transonic airfoil flows the best accuracy on coarse meshes for aerodynamic coefficients is obtained with a simple MATD scheme.

Swanson, R. C.↗
NASA NTRS1995

Pressure updating methods for the steady-state fluid equations

We consider the steady state equations for a compressible fluid. Since we wish to solve for a range of speeds we must consider the equations in conservation form. For transonic speeds these equations are of mixed type. Hence, the usual approach is to add time derivatives to the steady state equations and then march these equations in time. One then adds a time derivative of the density to the continuity equation, a derivative of the momentum to the momentum equation and a derivative of the total energy to the energy equation. This choice is dictated by the time consistent equations. However, since we are only interested in the steady state this is not necessary. Thus we shall consider the possibility of adding a time derivative of the pressure to the continuity equation and similar modifications for the energy equation. This can then be generalized to adding combinations of time derivatives to each equation since these vanish in the steady state. When using acceleration techniques such as residual smoothing, multigrid, etc. these are applied to the pressure rather than the density. Hence, the code duplicates the behavior of the incompressible equations for low speeds.

Fiterman, A.↗
NASA NTRS1994

Boundary Conditions for Jet Flow Computations

Ongoing activities are focused on capturing the sound source in a supersonic jet through careful large eddy simulation (LES). One issue that is addressed is the effect of the boundary conditions, both inflow and outflow, on the predicted flow fluctuations, which represent the sound source. In this study, we examine the accuracy of several boundary conditions to determine their suitability for computations of time-dependent flows. Various boundary conditions are used to compute the flow field of a laminar axisymmetric jet excited at the inflow by a disturbance given by the corresponding eigenfunction of the linearized stability equations. We solve the full time dependent Navier-Stokes equations by a high order numerical scheme. For very small excitations, the computed growth of the modes closely corresponds to that predicted by the linear theory. We then vary the excitation level to see the effect of the boundary conditions in the nonlinear flow regime.

Hayder, M. E.↗
NASA NTRS1993

Central difference TVD schemes for time dependent and steady state problems

A Runge-Kutta formula in time is presently used to advance schemes in which central differences are used to solve the time-dependent Euler equations; a second difference is added near shocks as an artificial viscosity to reduce the given scheme to a first-order upwind one at shocks. A matrix-valued dissipation is introduced and compared with the scalar viscosity; a connection is shown between this artificial viscosity and flux limiters. The use of various flux limiters for this central difference scheme is compared.

Jorgenson, P.↗
NASA NTRS1993

Aspects of a high-resolution scheme for the Navier-Stokes equations

In this paper we emphasize the importance of the form of the numerical dissipation model in computing accurate viscous flow solutions. A high-resolution scheme for viscous flows based on three-point central differencing and a matrix dissipation is considered. The various components of this scheme, including 'entropy fix', limiter function, and boundary-point dissipation are discussed. By analyzing boundary-point dissipation stencils, we confirm that with the matrix dissipation model the normal numerical dissipation terms in the streamwise momentum equation are independent of the Reynolds number. Such independence is not achieved with a scalar dissipation form. The accuracy of the central-difference scheme, with and without matrix dissipation, and the flux-difference split scheme of Roe, which is classified as a high-resolution scheme, is compared. For this comparison, three high Reynolds number laminar flows are considered. Solutions of the Navier-Stokes equations are obtained for low-speed flow over a flat plate, transonic flow over an airfoil with transition near the leading edge, and hypersonic flow over a compression ramp. The emphasis of the comparison is primarily on the details of the viscous flows. The necessity of the high-resolution property is revealed.

Swanson, R. C.↗
NASA NTRS1993

Preconditioning and the limit to the incompressible flow equations

The use of preconditioning methods to accelerate the convergence to a steady state for both the incompressible and compressible fluid dynamic equations are considered. The relation between them for both the continuous problem and the finite difference approximation is also considered. The analysis relies on the inviscid equations. The preconditioning consists of a matrix multiplying the time derivatives. Hence, the steady state of the preconditioned system is the same as the steady state of the original system. For finite difference methods the preconditioning can change and improve the steady state solutions. An application to flow around an airfoil is presented.

Turkel, E.↗
NASA NTRS1992

An effective multigrid method for high-speed flows

The use is considered of a multigrid method with central differencing to solve the Navier-Stokes equations for high speed flows. The time dependent form of the equations is integrated with a Runge-Kutta scheme accelerated by local time stepping and variable coefficient implicit residual smoothing. Of particular importance are the details of the numerical dissipation formulation, especially the switch between the second and fourth difference terms. Solutions are given for 2-D laminar flow over a circular cylinder and a 15 deg compression ramp.

Swanson, R. C.↗
NASA NTRS1992

Central difference TVD and TVB schemes for time dependent and steady state problems

We use central differences to solve the time dependent Euler equations. The schemes are all advanced using a Runge-Kutta formula in time. Near shocks, a second difference is added as an artificial viscosity. This reduces the scheme to a first order upwind scheme at shocks. The switch that is used guarantees that the scheme is locally total variation diminishing (TVD). For steady state problems it is usually advantageous to relax this condition. Then small oscillations do not activate the switches and the convergence to a steady state is improved. To sharpen the shocks, different coefficients are needed for different equations and so a matrix valued dissipation is introduced and compared with the scalar viscosity. The connection between this artificial viscosity and flux limiters is shown. Any flux limiter can be used as the basis of a shock detector for an artificial viscosity. We compare the use of the van Leer, van Albada, mimmod, superbee, and the 'average' flux limiters for this central difference scheme. For time dependent problems, we need to use a small enough time step so that the CFL was less than one even though the scheme was linearly stable for larger time steps. Using a total variation bounded (TVB) Runge-Kutta scheme yields minor improvements in the accuracy.

Jorgenson, P.↗
NASA NTRS1992

Central difference TVD and TVB schemes for time dependent and steady state problems

We use central differences to solve the time dependent Euler equations. The schemes are all advanced using a Runge-Kutta formula in time. Near shocks, a second difference is added as an artificial viscosity. This reduces the scheme to a first order upwind scheme at shocks. The switch that is used guarantees that the scheme is locally total variation diminishing (TVD). For steady state problems it is usually advantageous to relax this condition. Then small oscillations do not activate the switches and the convergence to a steady state is improved. To sharpen the shocks, different coefficients are needed for different equations and so a matrix valued dissipation is introduced and compared with the scalar viscosity. The connection between this artificial viscosity and flux limiters is shown. Any flux limiter can be used as the basis of a shock detector for an artificial viscosity. We compare the use of the van Leer, van Albada, mimmod, superbee, and the 'average' flux limiters for this central difference scheme. For time dependent problems, we need to use a small enough time step so that the CFL was less than one even though the scheme was linearly stable for larger time steps. Using a total variation bounded (TVB) Runge-Kutta scheme yields minor improvements in the accuracy.

Jorgenson, P.↗
NASA NTRS1991

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that removes the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional flow over a blunt biconic.

Turkel, E.↗
NASA NTRS1991

An effective multigrid method for high-speed flows

The use is considered of a multigrid method with central differencing to solve the Navier-Stokes equations for high speed flows. The time dependent form of the equations is integrated with a Runge-Kutta scheme accelerated by local time stepping and variable coefficient implicit residual smoothing. Of particular importance are the details of the numerical dissipation formulation, especially the switch between the second and fourth difference terms. Solutions are given for 2-D laminar flow over a circular cylinder and a 15 deg compression ramp.

Swanson, R. C.↗
NASA NTRS1991

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time-dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that remove the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock-capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional viscous flow over a blunt biconic.

Turkel, E.↗
NASA NTRS1989

Global properties of pseudospectral methods

The present application of polynomial interpolation methods to function-approximation and numerical solutions for hyperbolic and elliptic PDEs allows the explicit construction of the derivative matrix for a general sequence of collocation points. An evaluation of the effect of several factors on the performance of these methods indicates an inability to interpret global methods in terms of local ones; the accuracy of the approximation will differ when the function's large gradients occur near the center of the region or near the boundary, irrespective of the boundary vicinity's collocation-point density.

Solomonoff, A.↗
NASA NTRS1987

Wave phenomena in a high Reynolds number compressible boundary layer

The behavior of spatially unstable waves in a high Reynolds number compressible laminar boundary layer is investigated by solution of the laminar two-dimensional compressible Navier-Stokes equations (solved to fourth-order accuracy) over a flat plate with a fluctuating disturbance generated at the inflow. A significant nonlinear distortion is produced, in qualitative agreement with experimental data. It is shown that increasing compressibility can significantly stabilize the flow over a flat plate, and that the mechanism of phase cancellation is a viable mechanism for the control of growing disturbances.

Bayliss, A.↗