Nonconforming Schwarz-Spectral Element Methods For Incompressible Flow

Abstract

We present scalable implementations of spectral-element-based Schwarz overlapping (overset) methods for the incompressible Navier-Stokes (NS) equations. Our SEM-based overset grid method is implemented at the level of the NS equations, which are advanced independently within separate subdomains using interdomain velocity and pressure boundary-data exchanges at each timestep or sub-timestep. Central to this implementation is a general, robust, and scalable interpolation routine, {\em gslib-findpts}, that rapidly determines the computational coordinates (processor $p$, element number $e$, and local coordinates $(r,s,t)\in \hat{\Omega }:=[-1,1{]}^3$) for any arbitrary point ${\mathbf{x}}^{\ast }=(x^{\ast },y^{\ast },z^{\ast })\in \Omega \subset {\mathrm{I}\mathrm{R}}^3$. The communication kernels in $gslib$ execute with at most $\log P$ complexity for $P$ MPI ranks, have scaled to $P>10^6$, and obviate the need for development of any additional MPI-based code for the Schwarz implementation. The original interpolation routine has been extended to account for multiple overlapping domains. The new implementation discriminates the possessing subdomain by distance to the domain boundary, such that the interface boundary data is taken from the inner-most interior points. We present application of this approach to several heat transfer and fluid dynamic problems, discuss the computation/communication complexity and accuracy of the approach, and present performance measurements for $P>12,000$.

12 pages, 10 figures, submitted to Computers & Fluids for special issue for the 30th International Conference on Parallel Computational Fluid Dynamics

References in corpus (3)

• Direct numerical simulation of particulate flows with an overset grid method
• Particle transport in turbulent curved pipe flow
• High-order overset grid method for detecting particle impaction on a cylinder in a cross flow

Cited by in corpus (11)

• Efficient Exascale Discretizations: High-Order Finite Element Methods
• Simulation-Driven Optimization of High-Order Meshes in ALE Hydrodynamics
• Multirate Timestepping for the Incompressible Navier-Stokes Equations in Overlapping Grids
• hr-adaptivity for nonconforming high-order meshes with the target matrix optimization paradigm
• Accelerating High-Order Mesh Optimization Using Finite Element Partial Assembly on GPUs
• High-Order Mesh Morphing for Boundary and Interface Fitting to Implicit Geometries
• Combining high-order metric interpolation and geometry implicitization for curved $r$-adaption
• Direct Numerical Simulation of Rotating Ellipsoidal Particles using Moving Nonconforming Schwarz-Spectral Element Method
• Stability analysis of a singlerate and multirate predictor-corrector scheme for overlapping grids
• Adaptive Surface Fitting and Tangential Relaxation for High-Order Mesh Optimization
• Objective early identification of kinematic instabilities in shear flows