Engineering PapersSearch

Engineering topics

Mcdonald, Jeffrey D.

Publications and source records attributed to Mcdonald, Jeffrey D..

Data-parallel lower-upper relaxation method for reacting flows

The implicit lower-upper symmetric Gauss-Seidel (LU-SGS) method of Yoon and Jameson is modified for use on massively parallel computers. The method has been implemented on the Thinking Machines CM-5 and the MasPar MP-1 and MP-2, where large percentages of the theoretical peak floating point performance are obtained. It is shown that the new data-parallel LU relaxation method has better convergence properties than the original method for two different inviscid compressible flow simulations. The convergence is also improved for five-species reacting air computations. The performance of the method on various partitions of the CM-5 and on the MasPar computers is discussed. The new method shows promise for the efficient simulation of very large perfect gas and reacting flows.

Candler, Graham V.

Models of thermal relaxation mechanics for particle simulation methods

A comparative study is conducted of the simulation of molecular vibrational relaxation in rarefied flows by means of existing models and a novel phenomenological microscopic model. This improved technique iterates between translation-rotation and rotation-vibration exchanges; equilibrium is in this way promoted if the model for the latter process is compatible with quantized oscillators. The iteration-equipartition model retains computational simplicity and promotes thermal equilibrium, even in the case of nondegenerate anharmonic quantized oscillators and multispecies gas mixtures.

Haas, Brian L.

Statistical discrete particle simulation studies

A particle simulation code was developed that is suited to Multiple Instruction Multiple Data multiprocessor computers. The resulting code is designed to be portable between a large class of multiprocessor computer architectures, including hypercubes (e.g., Intel iPSC/860), shared memory machines (e.g., Cray, SGI, DASH), and simple uniprocessors (e.g., SUN, VAX, IBM PC). This code provides a powerful general framework that can be adapted to specific applications, such as thermomechanical modeling, general 3-D geometry support, or visualization support via a codeveloped CPlot data reduction and visualization system. Chemistry models were implemented and validated with comparison to continuum solutions of thermochemically relaxing gas mixtures and with experimental results from the high speed flow about a circular cylinder. Furthermore, the exchange of energy between translational and internal modes was modeled. New models were developed, combining greater efficiency than earlier phenomenological models. They offer greater physical detail by addressing vibrational energy as a discretely distributed quantity.

Mcdonald, Jeffrey D.

Particle simulation in a multiprocessor environment

A parallel implementation of a particle simulation method that is portable between a wide class of multiprocessor computers is presented. A fine grain spatial decomposition is utilized where several subdomains having a regular structure are computed at each processing node. This leads directly to an efficient and straightforward load balancing scheme if the number of subdomains at each processor is permitted to vary in an appropriate manner. Three dimensional simulations incorporating full thermochemical nonequilibrium are possible using the resulting code. Vectorizable algorithms are retained from earlier work allowing efficient use of deeply pipelined node processors where available. Performance results are presented from three different machine architectures demonstrating the portability of the code. On a 128-node Intel iPSC/860, performance is twice that of a single Cray-Y/MP CPU running a highly vectorized simulation code. Speedup is linear over the full range of number of processors on all target machines, indicating scalability of the method to higher degrees of parallelism.

Mcdonald, Jeffrey D.

Validation of chemistry models employed in a particle simulation method

The chemistry models employed in a statistical particle simulation method, as implemented in the Intel iPSC/860 multiprocessor computer, are validated and applied. Chemical relaxation of five-species air in these reservoirs involves 34 simultaneous dissociation, recombination, and atomic-exchange reactions. The reaction rates employed in the analytic solutions are obtained from Arrhenius experimental correlations as functions of temperature for adiabatic gas reservoirs in thermal equilibrium. Favorable agreement with the analytic solutions validates the simulation when applied to relaxation of O2 toward equilibrium in reservoirs dominated by dissociation and recombination, respectively, and when applied to relaxation of air in the temperature range 5000 to 30,000 K. A flow of O2 over a circular cylinder at high Mach number is simulated to demonstrate application of the method to multidimensional reactive flows.

Haas, Brian L.

Three dimensional Discrete Particle Simulation about the AFE geometry

The Discrete Particle Simulation method, due to Baganoff, has recently been extended to allow representation of gases composed of multiple species, to general power-law molecular interactions and to permit flows in thermal non-equilibrium. Particular attention has been paid to the implementation of this physics while retaining the efficiency of the original algorithm. Here, the enhanced algorithm is applied to the simulation of the flow field about the Aeroassisted Flight Experiment (AFE) vehicle with the same flight parameters as in a previous paper. The enhancements to the algorithm are introduced and comparisons are made to the previous calculation.

Feiereisen, William J.

Application of a vectorized particle simulation in high-speed near-continuum flow

An efficient particle simulation technique, developed for use on vector architecture based supercomputers for studying hypersonic rarefied gas flows is employed to simulate the complex wake generated by Mach six flow over a 10 deg half-angle wedge for freestream Reynolds numbers of 1780 and 3560. Data obtained are compared against higher Reynolds number experimental results. Simulations utilized as many as 10 to the 5th computational cells and 10 to the 7th simulated particles having power-law interaction potentials. A code performance of 1.8-2.4 microsec of Cray-2 CPU time to process a single particle per timestep is achieved. Diffuse adiabatic and isothermal wedge surface models are used in this investigation. Although the wedge geometry is two-dimensional, the simulation incorporates a width-wise direction, resulting in a three-dimensional computation.

Woronowicz, Michael S.

Three dimensional discrete particle simulation of an AOTV

A discrete Particle Simulation method, recently formulated by Baganoff, is discussed in the context of its application to the simulation of the flow field about the Aeroassisted Flight Experiment (AFE). As a basis for discussion the current algorithm is first described. Because of the use of a cubic Cartesian mesh, the representation of the geometry is different than that of other particle methods and an algorithm for its generation is discussed. The method is applied to test problems and then to the AFE calculation with the use of 9.52 million particles and 432,000 cells.

Feiereisen, William J.

Vectorization of a particle simulation method for hypersonic rarefied flow

An efficient particle simulation technique for hypersonic rarefied flows is presented at an algorithmic and implementation level. The implementation is for a vector computer architecture, specifically the Cray-2. The method models an ideal diatomic Maxwell molecule with three translational and two rotational degrees of freedom. Algorithms are designed specifically for compatibility with fine grain parallelism by reducing the number of data dependencies in the computation. By insisting on this compatibility, the method is capable of performing simulation on a much larger scale than previously possible. A two-dimensional simulation of supersonic flow over a wedge is carried out for the near-continuum limit where the gas is in equilibrium and the ideal solution can be used as a check on the accuracy of the gas model employed in the method. Also, a three-dimensional, Mach 8, rarefied flow about a finite-span flat plate at a 45 degree angle of attack was simulated. It utilized over 10 to the 7th particles carried through 400 discrete time steps in less than one hour of Cray-2 CPU time. This problem was chosen to exhibit the capability of the method in handling a large number of particles and a true three-dimensional geometry.

Mcdonald, Jeffrey D.