Multi-grid domain decomposition approach for solution of Navier–Stokes equations in primitive variable form

A new multi-grid (two-grid) pseudospectral element method has been carried out for solution of incompressible flow in terms of primitive variable formulation. The main objective of the proposed method is to apply the multi-grid technique solving the incompressible flow problems associated with three commonly encountered multi-grid environments. In domain decomposition terminology, it includes (I) partially overlapped subdomains, each of which has same types of grids; (II) partially overlapped subdomains, each of which has different types of grids; (III) local adaptive subdomains fully overlapped with the original computational domain (composite grids). The approach for flow problems, complex geometry or not, is to first divide the computational domain into a number of subdomains with the inter-overlapping area (partially or fully overlapped). In categories I and II, the fine-grid or coarse-grid subdomains can be defined by their representation, while in category III the fine-grid or coarse-grid subdomains are defined as usual. Next, implement the Schwarz Alternating Procedure (SAP) to exchange the data among subdomains, where the coarse-grid correction is used to remove the high frequency error that occurs when the data interpolation from the fine-grid subdomain to the coarse-grid subdomain is conducted. The strategy behind the coarse-grid correction is to adopt the operator of the divergence of velocity field, which intrinsically links the pressure equation, into this process. The solution of each subdomain can be efficiently solved by the direct (or iterative) eigenfunction expansion technique or preconditioned method with the least storage requirement, i.e. O(N²) in 2-D. Numerical results of (i) driven cavity flow (Re = 100,400) with Cartesian grids (category I) in each subdomain, (ii) driven cavity flow (Re = 3200) with local adaptive grids (category III) in each subdomain, and (iii) flow over a cylinder (Re = 250) with ‘O’ grids in one subdomain and Cartesian grids in another (category II) will be presented in the paper to account for the versatility of the proposed multi-grid method.