Computation of acoustic waves in a jet

A numerical treatment of acoustic waves in a jet is described. The full time dependent Euler equations are used in both linear and nonlinear formulations. The computational region of integration is artificially bounded and boundary conditions are developed to simulate outgoing waves and to enable the computational domain to be substantially restricted. Higher order methods and coordinate transformations are introduced to further reduce the number of grid points as well as to increase the efficiency of the program. Numerical results are presented for time harmonic sources as well as for sources with more complicated time dependence

Global collocation methods for approximation and the solution of partial differential equations

Polynomial interpolation methods are applied both to the approximation of functions and to the numerical solutions of hyperbolic and elliptic partial differential equations. The derivative matrix for a general sequence of the collocation points is constructed. The approximate derivative is then found by a matrix times vector multiply. The effects of several factors on the performance of these methods including the effect of different collocation points are then explored. The resolution of the schemes for both smooth functions and functions with steep gradients or discontinuities in some derivative are also studied. The accuracy when the gradients occur both near the center of the region and in the vicinity of the boundary is investigated. The importance of the aliasing limit on the resolution of the approximation is investigated in detail. Also examined is the effect of boundary treatment on the stability and accuracy of the scheme

Simulation of the fluctuating field of a forced jet

The fluctuating field of a jet excited by transient mass injection is simulated numerically. The model is developed by expanding the state vector as a mean state plus a fluctuating state. Nonlinear terms are not neglected and the effect of nonlinearity is studied. The results show a significant spectral broadening in the flow field due to the nonlinearity. In addition, large scale structures are broken down into smaller scales

Boundary conditions for the numerical solution of elliptic equations in exterior regions

Elliptic equations in exterior regions frequently require a boundary condition at infinity to ensure the well-posedness of the problem. Examples of practical applications include the Helmholtz equation and Laplace's equation. Computational procedures based on a direct discretization of the elliptic problem require the replacement of the condition on a finite artificial surface. Direct imposition of the condition at infinity along the finite boundary results in large errors. A sequence of boundary conditions is developed which provides increasingly accurate approximations to the problem in the infinite domain. Estimates of the error due to the finite boundary are obtained for several cases. Computations are presented which demonstrate the increased accuracy that can be obtained by the use of the higher order boundary conditions. The examples are based on a finite element formulation but finite difference methods can also be used

The numerical solution of the Helmholtz equation for wave propagation problems in underwater acoustics

The Helmholtz Equation (-delta-K(2)n(2))u=0 with a variable index of refraction, n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. A numerical algorithm was developed and a computer code implemented that can effectively solve this equation in the intermediate frequency range. The equation is discretized using the finite element method, thus allowing for the modeling of complicated geometrices (including interfaces) and complicated boundary conditions. A global radiation boundary condition is imposed at the far field boundary that is exact for an arbitrary number of propagating modes. The resulting large, non-selfadjoint system of linear equations with indefinite symmetric part is solved using the preconditioned conjugate gradient method applied to the normal equations. A new preconditioner is developed based on the multigrid method. This preconditioner is vectorizable and is extremely effective over a wide range of frequencies provided the number of grid levels is reduced for large frequencies. A heuristic argument is given that indicates the superior convergence properties of this preconditioner

On accuracy conditions for the numerical computation of waves

The Helmholtz equation (Delta + K(2)n(2))u = f with a variable index of refraction n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. Such problems can be solved numerically by first truncating the given unbounded domain and imposing a suitable outgoing radiation condition on an artificial boundary and then solving the resulting problem on the bounded domain by direct discretization (for example, using a finite element method). In practical applications, the mesh size h and the wave number K, are not independent but are constrained by the accuracy of the desired computation. It will be shown that the number of points per wavelength, measured by (Kh)(-1), is not sufficient to determine the accuracy of a given discretization. For example, the quantity K(3)h(2) is shown to determine the accuracy in the L(2) norm for a second-order discretization method applied to several propagation models

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

Feasibility study of a procedure to detect and warn of low level wind shear

A Doppler radar system which provides an aircraft with advanced warning of longitudinal wind shear is described. This system uses a Doppler radar beamed along the glide slope linked with an on line microprocessor containing a two dimensional, three degree of freedom model of the motion of an aircraft including pilot/autopilot control. The Doppler measured longitudinal glide slope winds are entered into the aircraft motion model, and a simulated controlled aircraft trajectory is calculated. Several flight path deterioration parameters are calculated from the computed aircraft trajectory information. The aircraft trajectory program, pilot control models, and the flight path deterioration parameters are discussed. The performance of the computer model and a test pilot in a flight simulator through longitudinal and vertical wind fields characteristic of a thunderstorm wind field are compared

Stability and control of compressible flows over a surface with concave-conves curvature

The active control of spatially unstable disturbances in a laminar, two-dimensional, compressible boundary layer over a curved surface is numerically simulated. The control is effected by localized time-periodic surface heating. We consider two similar surfaces of different heights with concave-convex curvature. In one, the height is sufficiently large so that the favorable pressure gradient is sufficient to stabilize a particular disturbance. In the other case the pressure gradient induced by the curvature is destabilizing. It is shown that by using active control that the disturbance can be stabilized. The results demonstrate that the curvature induced mean pressure gradient significantly enhances the receptivity of the flow localized time-periodic surface heating and that this is a potentially viable mechanism in air

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