Roy S. Baty

A two-dimensional fast solver for arbitrary vortex distributions

A method which is capable of an efficient calculation of the two-dimensional stream function and velocity field produced by a large system of vortices is presented in this report. This work is based on the adaptive scheme of Carrier, Greengard, and Rokhlin with the added feature that the evaluation or target points do not have to coincide with the location of the source or vortex positions. A simple algorithm based on numerical experiments has been developed to optimize the method for cases where the number of vortices N{sub V} differs significantly from the number of target points N{sub T}. The ability to specify separate source and target fields provides an efficient means for calculating boundary conditions, trajectories of passive scalar quantities, and stream-function plots, etc. Test cases have been run to benchmark the truncation errors and CPU time savings associated with the method. For six terms in the series expansions, non-dimensional truncation errors for the magnitudes of the complex potential and velocity fields are on the order of 10{sup {minus}5} and 10{sup {minus}3} respectively. The authors found that the CPU time scales as {radical}(N{sub V}N{sub T}) for N{sub V}/N{sub T} in the range of 0.1 to 10. For {radical}(N{sub V}N{sub T}) less than 200, there is virtually no CPU time savings while for {radical}N{sub V}N{sub T} roughly equal to 20,000, the fast solver obtains solutions in about 1% of the time required for the direct solution technique depending somewhat upon the configuration of the vortex field and the target field.

Instability of jets of arbitrary geometry

This paper describes a calculation technique for determining the stability of jets of arbitrary cross section. In particular, elliptic and rectangular jets are considered. The numerical procedure involves both a conformal transformation between the computational domain and the physical plane and a solution of the transformed stability equation in the computational domain. Modern, efficient, conformal mappings are used for both simply and doubly connected domains. The numerical solution is based on a finite difference/pseudospectral discretization of the stability equation. The technique is verified by comparison with previous calculations for circular and elliptic jets. Calculations are performed for the stability of elliptic and rectangular jets of aspect ratio 2. Growth rates, phase velocities, and pressure eigenfunctions are calculated.

Global numerical prediction of bursting frequency in turbulent boundary layers

The frequencies of the bursting events associated with the streamwise coherent structures of spatially developing incompressible turbulent boundary layers were predicted. The structures were modeled as wavelike disturbances associated with the turbulent mean flow using a direct‐resonance theory. Global numerical solutions for the resonant eigenmodes of the Orr‐Sommerfeld and the vertical vorticity equations were developed. The global method involves the use of second and fourth order accurate finite difference formulae for the differential equations as well as the boundary conditions. The predicted resonance frequencies were found to agree very well with previous results using a local shooting technique and measured data.

A linear shock cell model for non-circular jets using conformal mapping with a pseudo-spectral hybrid scheme

The shock structure in non-circular supersonic jets is predicted using a linear model. This model includes the effects of the finite thickness of the mixing layer and the turbulence in the jet shear layer. A numerical solution is obtained using a conformal mapping grid generation scheme with a hybrid pseudo-spectral discretization method. The uniform pressure perturbation at the jet exit is approximated by a Fourier-Mathieu series. The pressure at downstream locations is obtained from an eigenfunction expansion that is matched to the pressure perturbation at the jet exit. Results are presented for a circular jet and for an elliptic jet of aspect ratio 2.0. Comparisons are made with experimental data.