130 research outputs found
A fast immersed boundary method for external incompressible viscous flows using lattice Green's functions
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.Comment: 32 pages, 9 figures; preprint submitted to Journal of Computational
Physic
Stable Fourth Order Stream-Function Methods for Incompressible Flows with Boundaries
Fourth-order stream-function methods are proposed for the time dependent, incompressible Navier-Stokes and Boussinesq equations. Wide difference stencils are used instead of compact ones and the boundary terms are handled by extrapolating the stream-function values inside the computational domain to grid points outside, up to fourth-order in the noslip condition. Formal error analysis is done for a simple model problem, showing that this extrapolation introduces numerical boundary layers at fifth-order in the stream-function. The fourth-order convergence in velocity of the proposed method for the full problem is shown numerically
- …