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

Accurate computation of surface stresses and forces with immersed boundary methods

Many immersed boundary methods solve for surface stresses that impose the velocity boundary conditions on an immersed body. These surface stresses may contain spurious oscillations that make them ill-suited for representing the physical surface stresses on the body. Moreover, these inaccurate stresses often lead to unphysical oscillations in the history of integrated surface forces such as the coefficient of lift. While the errors in the surface stresses and forces do not necessarily affect the convergence of the velocity field, it is desirable, especially in fluid-structure interaction problems, to obtain smooth and convergent stress distributions on the surface. To this end, we show that the equation for the surface stresses is an integral equation of the first kind whose ill-posedness is the source of spurious oscillations in the stresses. We also demonstrate that for sufficiently smooth delta functions, the oscillations may be filtered out to obtain physically accurate surface stresses. The filtering is applied as a post-processing procedure, so that the convergence of the velocity field is unaffected. We demonstrate the efficacy of the method by computing stresses and forces that converge to the physical stresses and forces for several test problems

Fish schooling as a basis for vertical axis wind turbine farm design

Most wind farms consist of horizontal axis wind turbines (HAWTs) due to the high power coefficient (mechanical power output divided by the power of the free-stream air through the turbine cross-sectional area) of an isolated turbine. However when in close proximity to neighbouring turbines, HAWTs suffer from a reduced power coefficient. In contrast, previous research on vertical axis wind turbines (VAWTs) suggests that closely-spaced VAWTs may experience only small decreases (or even increases) in an individual turbine's power coefficient when placed in close proximity to neighbours, thus yielding much higher power outputs for a given area of land. A potential flow model of inter-VAWT interactions is developed to investigate the effect of changes in VAWT spatial arrangement on the array performance coefficient, which compares the expected average power coefficient of turbines in an array to a spatially-isolated turbine. A geometric arrangement based on the configuration of shed vortices in the wake of schooling fish is shown to significantly increase the array performance coefficient based upon an array of 16x16 wind turbines. Results suggest increases in power output of over one order of magnitude for a given area of land as compared to HAWTs.Comment: Submitted for publication in BioInspiration and Biomimetics. Note: The technology described in this paper is protected under both US and international pending patents filed by the California Institute of Technolog

A fast lattice Green's function method for solving viscous incompressible flows on unbounded domains

A computationally efficient method for solving three-dimensional, viscous, incompressible flows on unbounded domains is presented. The method formally discretizes the incompressible Navier–Stokes equations on an unbounded staggered Cartesian grid. Operations are limited to a finite computational domain through a lattice Green's function technique. This technique obtains solutions to inhomogeneous difference equations through the discrete convolution of source terms with the fundamental solutions of the discrete operators. The differential algebraic equations describing the temporal evolution of the discrete momentum equation and incompressibility constraint are numerically solved by combining an integrating factor technique for the viscous term and a half-explicit Runge–Kutta scheme for the convective term. A projection method that exploits the mimetic and commutativity properties of the discrete operators is used to efficiently solve the system of equations that arises in each stage of the time integration scheme. Linear complexity, fast computation rates, and parallel scalability are achieved using recently developed fast multipole methods for difference equations. The accuracy and physical fidelity of solutions are verified through numerical simulations of vortex rings

A parallel fast multipole method for elliptic difference equations

A new fast multipole formulation for solving elliptic difference equations on unbounded domains and its parallel implementation are presented. These difference equations can arise directly in the description of physical systems, e.g. crystal structures, or indirectly through the discretization of PDEs. In the analog to solving continuous inhomogeneous differential equations using Green's functions, the proposed method uses the fundamental solution of the discrete operator on an infinite grid, or lattice Green's function. Fast solutions O(N)O(N) are achieved by using a kernel-independent interpolation-based fast multipole method. Unlike other fast multipole algorithms, our approach exploits the regularity of the underlying Cartesian grid and the efficiency of FFTs to reduce the computation time. Our parallel implementation allows communications and computations to be overlapped and requires minimal global synchronization. The accuracy, efficiency, and parallel performance of the method are demonstrated through numerical experiments on the discrete 3D Poisson equation

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 3700 are used to verify the accuracy and physical fidelity of the formulation

Fish schooling as a basis for vertical axis wind turbine farm design

Most wind farms consist of horizontal axis wind turbines (HAWTs) due to the high power coefficient (mechanical power output divided by the power of the free-stream air through the turbine cross-sectional area) of an isolated turbine. However when in close proximity to neighbouring turbines, HAWTs suffer from a reduced power coefficient. In contrast, previous research on vertical axis wind turbines (VAWTs) suggests that closely-spaced VAWTs may experience only small decreases (or even increases) in an individual turbine's power coefficient when placed in close proximity to neighbours, thus yielding much higher power outputs for a given area of land. A potential flow model of inter-VAWT interactions is developed to investigate the effect of changes in VAWT spatial arrangement on the array performance coefficient, which compares the expected average power coefficient of turbines in an array to a spatially-isolated turbine. A geometric arrangement based on the configuration of shed vortices in the wake of schooling fish is shown to significantly increase the array performance coefficient based upon an array of 16x16 wind turbines. Results suggest increases in power output of over one order of magnitude for a given area of land as compared to HAWTs.Comment: Submitted for publication in BioInspiration and Biomimetics. Note: The technology described in this paper is protected under both US and international pending patents filed by the California Institute of Technolog

Accurate computation of surface stresses and forces with immersed boundary methods

Many immersed boundary methods solve for surface stresses that impose the velocity boundary conditions on an immersed body. These surface stresses may contain spurious oscillations that make them ill-suited for representing the physical surface stresses on the body. Moreover, these inaccurate stresses often lead to unphysical oscillations in the history of integrated surface forces such as the coefficient of lift. While the errors in the surface stresses and forces do not necessarily affect the convergence of the velocity field, it is desirable, especially in fluid–structure interaction problems, to obtain smooth and convergent stress distributions on the surface. To this end, we show that the equation for the surface stresses is an integral equation of the first kind whose ill-posedness is the source of spurious oscillations in the stresses. We also demonstrate that for sufficiently smooth delta functions, the oscillations may be filtered out to obtain physically accurate surface stresses. The filtering is applied as a post-processing procedure, so that the convergence of the velocity field is unaffected. We demonstrate the efficacy of the method by computing stresses and forces that converge to the physical stresses and forces for several test problems

Immersed Boundary Lattice Green Function methods for External Aerodynamics

In this paper, we document the capabilities of a novel numerical approach - the immersed boundary lattice Green's function (IBLGF) method - to simulate external incompressible flows over complex geometries. This new approach is built upon the immersed boundary method and lattice Green's functions to solve the incompressible Navier-Stokes equations. We show that the combination of these two concepts allows the construction of an efficient and robust numerical framework for the direct numerical and large-eddy simulation of external aerodynamic problems at moderate to high-Reynolds numbers