Search NASA⌕ Search

Engineering topics

Hussaini, M. Y.

Publications and source records attributed to Hussaini, M. Y..

At least 163 records · Page 9

Navier-Stokes solution on the CYBER-203 by a pseudospectral technique

A three-level, time-split, mixed spectral/finite difference method for the numerical solution of the three-dimensional, compressible Navier-Stokes equations has been developed and implemented on the Control Data Corporation (CDC) CYBER-203. This method uses a spectral representation for the flow variables in the streamwise and spanwise coordinates, and central differences in the normal direction. The five dependent variables are interleaved one horizontal plane at a time and the array of their values at the grid points of each horizontal plane is a typical vector in the computation. The code is organized so as to require, per time step, a single forward-backward pass through the entire data base. The one-and two-dimensional Fast Fourier Transforms are performed using software especially developed for the CYBER-203.

Lambiotte, J. J.↗

Second-order explicit finite-difference methods for transient-flow analysis

Three second-order accurate numerical methods - MacCormack's method, Lambda scheme and Gabutti scheme - are introduced to solve the quasi-linear, hyperbolic partial differential equations describing transient flows in closed conduits. The details of these methods and the treatment of boundary conditions are presented and the results computed by using these methods for a typical piping system are compared. It is shown that for the same accuracy, second-order methods require considerably lesser number of computational nodes and computer time as compared to those required by the first-order methods.

Chaudhry, M. H.↗

Pseudospectral calculation of shock turbulence interactions

A Chebyshev-Fourier discretization with shock fitting is used to solve the unsteady Euler equations. The method is applied to shock interactions with plane waves and with a simple model of homogeneous isotropic turbulence. The plane wave solutions are compared to linear theory. Previously announced in STAR as N83-27952

Zang, T. A.↗

Iterative spectral methods and spectral solutions to compressible flows

A spectral multigrid scheme is described which can solve pseudospectral discretizations of self-adjoint elliptic problems in O(N log N) operations. An iterative technique for efficiently implementing semi-implicit time-stepping for pseudospectral discretizations of Navier-Stokes equations is discussed. This approach can handle variable coefficient terms in an effective manner. Pseudospectral solutions of compressible flow problems are presented. These include one dimensional problems and two dimensional Euler solutions. Results are given both for shock-capturing approaches and for shock-fitting ones.

Hussaini, M. Y.↗

Shock-fitted Euler solutions to shock-vortex interactions

The interaction of a shock wave with a hot spot, a single vortex and a vortex street is studied within the framework of the two dimensional compressible Euler equations. The numerical results obtained by the pseudospectral method and the finite difference MacCormack method are compared. In both the methods the shock wave is fitted as a boundary of the computational domain.

Salas, M. E.↗

Accurate numerical solution of compressible, linear stability equations

The present investigation is concerned with a fourth order accurate finite difference method and its application to the study of the temporal and spatial stability of the three-dimensional compressible boundary layer flow on a swept wing. This method belongs to the class of compact two-point difference schemes discussed by White (1974) and Keller (1974). The method was apparently first used for solving the two-dimensional boundary layer equations. Attention is given to the governing equations, the solution technique, and the search for eigenvalues. A general purpose subroutine is employed for solving a block tridiagonal system of equations. The computer time can be reduced significantly by exploiting the special structure of two matrices.

Malik, M. R.↗

Numerical computations of turbulence amplification in shock wave interactions

Numerical computations are presented which illustrate and test various effects pertinent to the amplification and generation of turbulence in shock wave-turbulent boundary layer interactions. Several fundamental physical mechanisms are identified. Idealizations of these processes are examined by nonlinear numerical calculations. The results enable some limits to be placed on the range of validity of existing linear theories. Additional results are given which are of a fundamentally nonlinear nature.

Zang, T. A.↗

Shock-fitted Euler solutions to shock-vortex interactions

The interaction of a planar shock wave with one or more vortexes is computed using a pseudospectral method and a finite difference method. The development of the spectral method is emphasized. In both methods the shock wave is fitted as a boundary of the computational domain. The results show good agreement between both computational methods. The spectral method is, however, restricted to smaller time steps and requires use of filtering techniques. Previously announced in STAR as N82-28061

Salas, M. D.↗

Multigrid method with weighted mean scheme

Multigrid algorithms based on the weighted mean scheme are developed for the solution of the two dimensional incompressible Navier-Stokes equations. They are applied to two typical problems encountered in engineering applications, namely, the convection diffusion problem of the Benard convection cell, and the driven cavity problem. An analysis of the smoothing rates and stability is given. The efficiency of the multigrid method is investigated.

Lustman, L. R.↗

Spectral multigrid methods for elliptic equations

An alternative approach which employs multigrid concepts in the iterative solution of spectral equations was examined. Spectral multigrid methods are described for self adjoint elliptic equations with either periodic or Dirichlet boundary conditions. For realistic fluid calculations the relevant boundary conditions are periodic in at least one (angular) coordinate and Dirichlet (or Neumann) in the remaining coordinates. Spectral methods are always effective for flows in strictly rectangular geometries since corners generally introduce singularities into the solution. If the boundary is smooth, then mapping techniques are used to transform the problem into one with a combination of periodic and Dirichlet boundary conditions. It is suggested that spectral multigrid methods in these geometries can be devised by combining the techniques.

Zang, T. A.↗

Vectored injection into compressible laminar and turbulent boundary layers

A fundamental analysis of two-dimensional supersonic boundary layer flow, both laminar and turbulent, is presented for a wide range of normal and nonnormal mass-transfer velocities. The analysis is based on the numerical solution of the Navier-Stokes equations, and results are compared with available theoretical and experimental data. Certain cases of practical importance, for which results are not presently available, are referred to.

Degani, D.↗

Mixed spectral/finite difference approximations for slightly viscous flows

A numerical scheme is presented, which employs a standard finite difference approximation for the viscous terms in high Reynolds number flows, and resorts to pseudo-spectral methods for the advection terms to greatly reduce the major source of numerical error without dramatically increasing computational cost. The spectral/finite difference (SFD) method evaluates the advection term and second-order differences to evaluate the diffusion term. The fully finite difference (FFD) method with second-order central differences on both terms is also used. The SFD method can handle strong shocks and can outperform the FFD method at moderate viscosity.

Zang, T. A.↗

Asymptotic features of shock-wave boundary-layer interaction

A semi-implicit method is applied to solve the Navier-Stokes equations numerically and to evaluate the features of the free-interaction phenomenon that occurs when a shock wave impinges on a Blasius boundary layer. Comparisons are made with predictions of the triple-deck theory and experiment. Results include pressure and skin-friction distribution in the free-interaction region for various values of Reynolds number.

Hussaini, M. Y.↗

Numerical experiments in boundary-layer stability

Numerical solution of the three-dimensional incompressible Navier-Stokes equations is used to study the instability of a flat-plate boundary layer in a manner analogous to the vibrating-ribbon experiments. Flow-field structures are observed which are very similar to those found in the vibrating-ribbon experiment to which computational initial conditions have been matched. Streamwise periodicity is assumed in the simulation so that the evolution occurs in time, but the events which constitute the instability are so similar to the spatially occurring ones of the laboratory that it seems clear the physical processes involved are the same. A spectral and finite difference numerical algorithm is employed in the simulation.

Wray, A.↗

Numerical simulation of supersonic cone flow at high angle of attack

A conical symmetry assumption is applied to the full Navier-Stokes equations resulting in an equation set containing time and two coordinate directions as independent variables. The set is integrated by use of a finite difference technique for the particular case of sharp cones at incidence. Solutions are obtained and compared with experiment for auxiliary conditions corresponding to both laminar and turbulent flows. Closure for the turbulent flow case is provided by use of a scalor eddy viscosity model based on the mixing length hypothesis. Modifications to the eddy viscosity model were found which led to excellent surface pressure and surface flow direction agreement for turbulent flow at low supersonic Mach numbers. It is apparent that further work on the turbulence model is necessary for agreement at higher supersonic Mach numbers.

Mcrae, D. S.↗

Supersonic viscous flow over cones at incidence

The conically symmetric Navier-Stokes equations are solved by MacCormack's method for the supersonic flow past sharp cones at incidence. To provide closure for the case of turbulent flows, a scalar eddy-viscosity model based on mixing length hypotheses is used. The results are compared with available experimental data.

Mcrae, D. S.↗

The general solution to the classical problem of the finite Euler-Bernoulli beam

An analytical solution is obtained for the problem of free and forced vibrations of a finite Euler-Bernoulli beam with arbitrary (partially fixed) boundary conditions. The effects of linear viscous damping, Winkler foundation, constant axial tension, a concentrated mass, and an arbitrary forcing function are included in the analysis. No restriction is placed on the values of the parameters involved, and the solution presented here contains all cited previous solutions as special cases.

Hussaini, M. Y.↗