A fast immersed boundary method for external incompressible viscous flows using lattice Green's functions

TL;DR

This paper introduces a novel parallel immersed boundary method utilizing lattice Green's functions for efficient simulation of three-dimensional viscous incompressible flows on unbounded domains, significantly reducing computational cost.

Contribution

The paper presents a new immersed boundary method that combines lattice Green's functions, adaptive grids, and specialized solvers for efficient, accurate flow simulations on unbounded domains.

Findings

Accurately simulates flows around flat plates and spheres at Reynolds numbers up to 3,700.

Reduces computational time by limiting operations to regions near the immersed surface.

Demonstrates high fidelity and efficiency in complex flow simulations.

Abstract

A new parallel, computationally efficient immersed boundary method for solving three-dimensional, viscous, incompressible flows on unbounded domains is presented. Immersed surfaces with prescribed motions are generated using the interpolation and regularization operators obtained from the discrete delta function approach of the original (Peskin's) immersed boundary method. Unlike Peskin's method, boundary forces are regarded as Lagrange multipliers that are used to satisfy the no-slip condition. The incompressible Navier-Stokes equations are discretized on an unbounded staggered Cartesian grid and are solved in a finite number of operations using lattice Green's function techniques. These techniques are used to automatically enforce the natural free-space boundary conditions and to implement a novel block-wise adaptive grid that significantly reduces the run-time cost of solutions by limiting operations to grid cells in the immediate vicinity and near-wake region of the immersed surface. These techniques also enable the construction of practical discrete viscous integrating factors that are used in combination with specialized half-explicit Runge-Kutta schemes to accurately and efficiently solve the differential algebraic equations describing the discrete momentum equation, incompressibility constraint, and no-slip constraint. Linear systems of equations resulting from the time integration scheme are efficiently solved using an approximation-free nested projection technique. The algebraic properties of the discrete operators are used to reduce projection steps to simple discrete elliptic problems, e.g. discrete Poisson problems, that are compatible with recent parallel fast multipole methods for difference equations. Numerical experiments on low-aspect-ratio flat plates and spheres at Reynolds numbers up to 3,700 are used to verify the accuracy and physical fidelity of the formulation.