An Immersed Interface Method for the Incompressible Navier-Stokes Equations

We present an immersed interface algorithm for the incompressible Navier Stokes equations. The interface is represented by cubic splines which are interpolated through a set of Lagrangian control points. The position of the control points is implicitly updated using the fluid velocity. The forces that the interface exerts on the fluid are computed from the constitutive relation of the interface and are applied to the fluid through jumps in the pressure and jumps in the derivatives of pressure and velocity. A projection method is used to time advance the Navier-Stokes equations on a uniform cartesian mesh. The Poisson-like equations required for the implicit solution of the diffusive and pressure terms are solved using a fast Fourier transform algorithm. The position of the interface is updated implicitly using a quasi-Newton method (BFGS) within each timestep. Several examples are presented to illustrate the flexibility of the presented approach.Singapore-MIT Alliance (SMA

An Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains

We present an immersed interface method for the incompressible Navier Stokes equations capable of handling rigid immersed boundaries. The immersed boundary is represented by a set of Lagrangian control points. In order to guarantee that the no-slip condition on the boundary is satisfied, singular forces are applied on the fluid at the immersed boundary. The forces are related to the jumps in pressure and the jumps in the derivatives of both pressure and velocity, and are interpolated using cubic splines. The strength of singular forces is determined by solving a small system of equations at each time step. The Navier-Stokes equations are discretized on a staggered Cartesian grid by a second order accurate projection method for pressure and velocity.Singapore-MIT Alliance (SMA

Brownian dynamics of rigid particles in an incompressible fluctuating fluid by a meshfree method

A meshfree Lagrangian method for the fluctuating hydrodynamic equations (FHEs) with fluid-structure interactions is presented. Brownian motion of the particle is investigated by direct numerical simulation of the fluctuating hydrodynamic equations. In this framework a bidirectional coupling has been introduced between the fluctuating fluid and the solid object. The force governing the motion of the solid object is solely due to the surrounding fluid particles. Since a meshfree formulation is used, the method can be extended to many real applications involving complex fluid flows. A three-dimensional implementation is presented. In particular, we observe the short and long-time behaviour of the velocity autocorrelation function (VACF) of Brownian particles and compare it with the analytical expression. Moreover, the Stokes-Einstein relation is reproduced to ensure the correct long-time behaviour of Brownian dynamics.Comment: 24 pages, 2 figure

A comparative study of immersed-boundary interpolation methods for a flow around a stationary cylinder at low Reynolds number

The accuracy and computational efficiency of various interpolation methods for the implementation of non grid-confirming boundaries is assessed. The aim of the research is to select an interpolation method that is both efficient and sufficiently accurate to be used in the simulation of vortex induced vibration of the flow around a deformable cylinder. Results are presented of an immersed boundary implementation in which the velocities near nonconfirming boundaries were interpolated in the normal direction to the walls. The flow field is solved on a Cartesian grid using a finite volume method with a staggered variable arrangement. The Strouhal number and Drag coefficient for various cases are reported. The results show a good agreement with the literature. Also, the drag coefficient and Strouhal number results for five different interpolation methods were compared it was shown that for a stationary cylinder at low Reynolds number, the interpolation method could affect the drag coefficient by a maximum 2% and the Strouhal number by maximum of 3%. In addition, the bi-liner interpolation method took about 2% more computational time per vortex shedding cycle in companion to the other methods

An efficient parallel immersed boundary algorithm using a pseudo-compressible fluid solver

We propose an efficient algorithm for the immersed boundary method on distributed-memory architectures, with the computational complexity of a completely explicit method and excellent parallel scaling. The algorithm utilizes the pseudo-compressibility method recently proposed by Guermond and Minev [Comptes Rendus Mathematique, 348:581-585, 2010] that uses a directional splitting strategy to discretize the incompressible Navier-Stokes equations, thereby reducing the linear systems to a series of one-dimensional tridiagonal systems. We perform numerical simulations of several fluid-structure interaction problems in two and three dimensions and study the accuracy and convergence rates of the proposed algorithm. For these problems, we compare the proposed algorithm against other second-order projection-based fluid solvers. Lastly, the strong and weak scaling properties of the proposed algorithm are investigated