Engineering topics

Zang, T. A.

Publications and source records attributed to Zang, T. A..

At least 73 records · Page 4

NASA NTRS1983

Preconditioners for the spectral multigrid method

The systems of algebraic equations which arise from spectral discretizations of elliptic equations are full and direct solutions of them are rarely feasible. Iterative methods are an attractive alternative because Fourier transform techniques enable the discrete matrix-vector products to be computed with nearly the same efficiency as is possible for corresponding but sparse finite difference discretizations. For realistic Dirichlet problems preconditioning is essential for acceptable convergence rates. A brief description of Chebyshev spectral approximations and spectral multigrid methods for elliptic problems is given. A survey of preconditioners for Dirichlet problems based on second-order finite difference methods is made. New preconditioning techniques based on higher order finite differences and on the spectral matrix itself are presented. The preconditioners are analyzed in terms of their spectra and numerical examples are presented.

Phillips, T. N.↗

NASA NTRS1983

Pseudospectral solution of two-dimensional gas-dynamic problems

Chebyshev pseudospectral methods are used to compute two dimensional smooth compressible flows. Grid refinement tests show that spectral accuracy can be obtained. Filtering is not needed if resolution is sufficiently high and if boundary conditions are carefully prescribed.

Kopriva, D. A.↗

NASA NTRS1983

Viscous, resistive magnetohydrodynamic stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one-dimensional magnetohydrodynamic quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds-like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds-like numbers. Marginal stability curves, growth rates versus Reynolds-like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result that appears general is that instability has been found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three-dimensional instabilities may exist, similar to those in Poiseuille and Couette flow.

Dahlburg, R. B.↗

NASA NTRS1983

Spectral methods for partial differential equations

Origins of spectral methods, especially their relation to the Method of Weighted Residuals, are surveyed. Basic Fourier, Chebyshev, and Legendre spectral concepts are reviewed, and demonstrated through application to simple model problems. Both collocation and tau methods are considered. These techniques are then applied to a number of difficult, nonlinear problems of hyperbolic, parabolic, elliptic, and mixed type. Fluid dynamical applications are emphasized.

Hussaini, M. Y.↗

NASA NTRS1983

Spectral methods for inviscid, compressible flows

Report developments in the application of spectral methods to two dimensional compressible flows are reviewed. A brief introduction to spectral methods -- their history and especially their implementation -- is provided. The stress is on those techniques relevant to transonic flow computation. The spectral multigrid iterative methods are discussed with application to the transonic full potential equation. Discontinuous solutions of the Euler equations are considered. The key element is the shock fitting technique which is briefly explained.

Hussaini, M. Y.↗

NASA NTRS1983

Spectral methods for the Euler equations

Spectral methods for compressible flows are introduced in relation to finite difference and finite element techniques within the framework of the method of weighted residuals. Current spectral collocation methods are put in historical context. The basic concepts of both Fourier and Chebyshev spectral collocation methods are provided. Filtering strategies for both shock-fitting and shock-capturing approaches are also presented. Fourier shock capturing techniques are evaluated using a one-dimensional, periodic astrophysical 'nozzle' problem. Examples of shock-fitting approaches include a shock/acoustic wave interaction, shock/vortex interaction, and the classical blunt body problem. While the shock capturing spectral method does not yet show a clear advantage over second-order finite differences, equivalent accuracy can be obtained using shock fitting with far fewer grid points.

Hussaini, M. Y.↗

NASA NTRS1983

Psuedospectral 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.

Zang, T. A.↗

NASA NTRS1983

Spectral multigrid methods for elliptic equations 2

A detailed description of spectral multigrid methods is provided. This includes the interpolation and coarse-grid operators for both periodic and Dirichlet problems. The spectral methods for periodic problems use Fourier series and those for Dirichlet problems are based upon Chebyshev polynomials. An improved preconditioning for Dirichlet problems is given. Numerical examples and practical advice are included.

Zang, T. A.↗

NASA NTRS1983

Spectral multigrid methods with applications to transonic potential flow

Spectral multigrid methods are demonstrated to be a competitive technique for solving the transonic potential flow equation. The spectral discretization, the relaxation scheme, and the multigrid techniques are described in detail. Significant departures from current approaches are first illustrated on several linear problems. The principal applications and examples, however, are for compressible potential flow. These examples include the relatively challenging case of supercritical flow over a lifting airfoil.

Streett, C. L.↗

NASA NTRS1983

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.

Zang, T. A.↗

NASA NTRS1983

Two- and Three-Dimensional Galaxy Simulations

Dynamics of galaxies with up to 100,000 stars are simulated in set of computer programs available from COSMIC. In these models, large number of simulation stars move in combined gravitational field produced by stars themselves (self-consistent field) or by external matter (fixed field).

Hohl, F.↗

NASA NTRS1983

Viscous, resistive MHD stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one dimensional magnetohydrodynamic (MHD) quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds like numbers. Marginal stability curves, growth rates versus Reynolds like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result which appears general is that instability was found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three dimensional instabilities may exist, similar to those in Poiseuille and Couette flow.

Dahlburg, R. B.↗

NASA NTRS1983

Pseudospectral calculation of shock turbulence interactions

Zang, T. A.↗

NASA NTRS1982

Spectral multigrid methods for elliptic equations

Previously announced in STAR as N81-33888

Zang, T. A.↗

NASA NTRS1982

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.↗

NASA NTRS1982

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.↗

NASA NTRS1982

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.

Salas, M. D.↗

NASA NTRS1982

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.↗