Parallel computation of three-dimensional flows using overlapping grids with adaptive mesh refinement

This paper describes an approach for the numerical solution of time-dependent partial differential equations in complex three-dimensional domains. The domains are represented by overlapping structured grids, and block-structured adaptive mesh refinement (AMR) is employed to locally increase the grid resolution. In addition, the numerical method is implemented on parallel distributed-memory computers using a domain-decomposition approach. The implementation is flexible so that each base grid within the overlapping grid structure and its associated refinement grids can be independently partitioned over a chosen set of processors. A modified bin-packing algorithm is used to specify the partition for each grid so that the computational work is evenly distributed amongst the processors. All components of the AMR algorithm such as error estimation, regridding, and interpolation are performed in parallel. The parallel time-stepping algorithm is illustrated for initial-boundary-value problems involving a linear advection-diffusion equation and the (nonlinear) reactive Euler equations. Numerical results are presented for both equations to demonstrate the accuracy and correctness of the parallel approach. Exact solutions of the advection-diffusion equation are constructed, and these are used to check the corresponding numerical solutions for a variety of tests involving different overlapping grids, different numbers of refinement levels and refinement ratios, and different numbers of processors. The problem of planar shock diffraction by a sphere is considered as an illustration of the numerical approach for the Euler equations, and a problem involving the initiation of a detonation from a hot spot in a T-shaped pipe is considered to demonstrate the numerical approach for the reactive case. For both problems, the accuracy of the numerical solutions is assessed quantitatively through an estimation of the errors from a grid convergence study. The parallel performance of the approach is examined for the shock diffraction problem.

[1]  William D. Henshaw,et al.  A High-Order Accurate Parallel Solver for Maxwell's Equations on Overlapping Grids , 2005, SIAM J. Sci. Comput..

[2]  Frederick Stern,et al.  RANS simulation of a container ship using a single-phase level-set method with overset grids and the prognosis for extension to a self-propulsion simulator , 2006 .

[3]  W. Henshaw,et al.  Composite overlapping meshes for the solution of partial differential equations , 1990 .

[4]  R. Meakin.  Moving body overset grid methods for complete aircraft tiltrotor simulations , 1993 .

[5]  David I. Gottlieb,et al.  The Theoretical Accuracy of Runge-Kutta Time Discretizations for the Initial Boundary Value Problem: A Study of the Boundary Error , 1995, SIAM J. Sci. Comput..

[6]  Jacob Yström,et al.  Some properties of the Upper Convected Maxwell model for viscoelastic fluid flow , 1993 .

[7]  Michael Shelley,et al.  A moving overset grid method for interface dynamics applied to non-Newtonian Hele-Shaw flow , 2004 .

[8]  B. Kreiss,et al.  Construction of a Curvilinear Grid , 1983 .

[9]  W. Henshaw,et al.  An adaptive numerical scheme for high-speed reactive flow on overlapping grids , 2003 .

[10]  Henrik Johansson,et al.  A performance characterization of load balancing algorithms for parallel SAMR applications , 2006 .

[11]  G. Starius,et al.  On composite mesh difference methods for hyperbolic differential equations , 1980 .

[12]  L. Fuchs,et al.  Calculation of flows using three‐dimensional overlapping grids and multigrid methods , 1995 .

[13]  Bernice Leland.  Grace , 1920, The Psychological clinic.

[14]  Joel H. Ferziger,et al.  NUMERICAL COMPUTATION OF UNSTEADY INCOMPRESSIBLE FLOW IN COMPLEX GEOMETRY USING A COMPOSITE MULTIGRID TECHNIQUE , 1991 .

[15]  John B. Bell,et al.  Parallelization of structured, hierarchical adaptive mesh refinement algorithms , 2000 .

[16]  William D. Henshaw.  Mappings for Overture A Description of the Mapping Class and Doc-umentation for Many Useful Mappings , 2002 .

[17]  Abdellatif Agouzal,et al.  A New Approach to Domain Decomposition Methods with Non-matching Grids , 2007 .

[18]  Robert Meakin,et al.  Composite Overset Structured Grids , 1998 .

[19]  G. Starius,et al.  Composite mesh difference methods for elliptic boundary value problems , 1977 .

[20]  G. Starius,et al.  Constructing orthogonal curvilinear meshes by solving initial value problems , 1977 .

[21]  Bruce Hendrickson,et al.  Dynamic load balancing in computational mechanics , 2000 .

[22]  William D. Henshaw,et al.  A study of detonation diffraction in the ignition-and-growth model , 2006 .

[23]  Jeffrey W. Banks,et al.  A high-resolution Godunov method for compressible multi-material flow on overlapping grids , 2006, J. Comput. Phys..

[24]  K. Olson,et al.  PARAMESH: A Parallel, Adaptive Grid Tool , 2006 .

[25]  Joel H. Saltz,et al.  A manual for the multiblock PARTI runtime primitives revision 4 , 1993 .

[26]  Kristi D. Brislawn,et al.  Adaptively-refined overlapping grids for the numerical solution of systems of hyperbolic conservation laws , 1995 .

[27]  Pieter G. Buning,et al.  Development of a large scale Chimera grid system for the Space Shuttle Launch Vehicle , 1993 .

[28]  D Kwak,et al.  Computational approach for probing the flow through artificial heart devices. , 1997, Journal of biomechanical engineering.

[29]  N. Anders Petersson,et al.  A numerical method to calculate the two-dimensional flow around an underwater obstacle , 1992 .

[30]  Thomas H. Pulliam,et al.  Recent enhancements to OVERFLOW , 1997 .

[31]  W. D. Henshaw,et al.  A study of detonation propagation and diffraction with compliant confinement , 2008 .

[32]  Andrew M. Wissink,et al.  Parallel clustering algorithms for structured AMR , 2005, J. Parallel Distributed Comput..

[33]  William D. Henshaw,et al.  Moving overlapping grids with adaptive mesh refinement for high-speed reactive and non-reactive flow , 2006, J. Comput. Phys..

[34]  Robert B. Ross,et al.  Using MPI-2: Advanced Features of the Message Passing Interface , 2003, CLUSTER.

[35]  Bharat K. Soni,et al.  Handbook of Grid Generation , 1998 .

[36]  J. L. Steger,et al.  On the use of composite grid schemes in computational aerodynamics , 1987 .

[37]  D W Schwendeman,et al.  Mechanisms of detonation formation due to a temperature gradient , 2002 .

[38]  P. Buning,et al.  Numerical simulation of the integrated space shuttle vehicle in ascent , 1988 .

[39]  Isidore Rigoutsos,et al.  An algorithm for point clustering and grid generation , 1991, IEEE Trans. Syst. Man Cybern..

[40]  Manish Parashar,et al.  Dynamic structured partitioning for parallel scientific applications with pointwise varying workloads , 2006, Proceedings 20th IEEE International Parallel & Distributed Processing Symposium.