Simulation Of Three-dimensional Flow Field Research Articles

Parallel dynamic load-balancing algorithm for three-dimensional adaptive unstructured grids

Adaptive local grid refinement/coarsening results in unequal distribution of work load among the processors of a parallel system. A novel method for balancing the load in cases of dynamically changing tetrahedral grids is developed. The approach employs local exchange of cells among processors to redistribute the load equally. An important part of the load-balancing algorithm is the method employed by a processor to determine which cells within its subdomain are to be exchanged. Two such methods are presented and compared. The strategy for load balancing is based on the divide-and-conquer approach that leads to an efficient parallel algorithm. This method is implemented on a distributed-memory multiple instruction multiple data system. I. Introduction C OMPUTATIONAL fluid dynamics (CFD) has advanced rapidly over the last two decades, and it is recognized as a valuable tool for engineering design. However, numerical simulation of three-dimensi onal flowfields remains very expensive even with the use of current vector supercomputers. Vector computers have accelerated computations only by one or two orders of magnitude compared with scalar machines. An evolving approach in computer architectures is the design of scalable massively parallel computers, wherein a number of processors are involved in executing different portions of the job. Parallel computing appears to be a promising approach for future design applications of CFD. State-of-the-art parallel architectures can be broadly classified into single instruction multiple data (SIMD), shared memory multiple instruction multiple data (MIMD), and partitioned memory MIMD architectures. Shared memory MIMD architectures such as the Cray Y-MP are extensions of pipelined vector processors with additional facilities to enable effective utilization of the available multiple processors. However, the scalability of such architectures is severely limited due to the inherent bottleneck of the common, shared memory. On the other hand, SIMD architectures such as the CM-2 are based on the lockstep paradigm of parallel computation wherein a large number of processors execute the same instructions on local data. Although this enhances scalability, the overhead associated with communication among these processors can sometimes prove to be a major bottleneck.1 Partitioned memory MIMD architectures provide a good compromise between the other two types of systems. The user has the flexibility to allocate data as well as tasks to each individual processor, thereby enabling fine tuning of the application to the underlying architecture. There is a price to be paid, however, in terms of additional user responsibility to coordinate the processors by exchanging the relevant information via message passing. In the case of an adaptive grid algo

Evaluation of a bounded high-resolution scheme for combustor flow computations

The present paper is concerned with the application of a high-resolution scheme for a finite volume discretization method. The capability of the chosen scheme is evaluated for the numerical simulation of three-dimensional flowfields. Special emphasis lies in the calculation of the processes in the mixing zone of gas turbine combustors. Thus, numerical results are compared with detailed measurements of velocity and temperature in a model mixing zone. The velocity field as well as the temperature field strongly depend on the scheme used for the discretization of the momentum and energy equations. The proposed high-resolution scheme and the well-known QUICK scheme yield similar results, which are superior to those obtained with the simple UPWIND scheme. For the simulation of turbulent transport, the standard fc, € model is employed. The effect of numerical diffusion caused by the usual UPWIND discretization of the convective transport of A and e is assessed by the application of the high-resolution scheme of the discretization of the A and e equations. As a result, it is found that the velocity and temperature fields are insensitive to the scheme that is applied for the discretization of the A and e equations.

A time-split finite-volume algorithm for three-dimensional flowfieldsimulation

A general finite-volume algorithm is developed for solving three-dimensional, time-dependent, compressible Navier-Stokes equations for high Reynolds number flows over an arbitrary geometry. This algorithm adapts MacCormack's (1982) explicit-implicit scheme to a time-split, three-dimensional finite-volume concept in a general coordinate system. It is shown that the thin-layer approximation in all three spatial directions significantly reduces the evaluation of viscous terms and allows the algorithm to solve more complicated geometries with all boundaries in two or all three directions. The calculated results using this method are found to be in good agreement with the experimental measurements of a blunt-fin induced shock wave and boundary-layer interaction problems. Observations of the existence of peak pressure, primary horseshoe and secondary vortices, and reversed supersonic zones show that computational fluid dynamics can effectively supplement the wind tunnel tests for aerodynamic design as well as for understanding basic fluid dynamics.