Students: D.V. Le
Faculty: J. Peraire
Collaborators: B.C. Khoo (National University, Singapore)
We employ an immersed interface method (IIM) for solving low-to-moderate Reynolds number laminar incompressible flows involving moving interfaces and rigid boundaries. One of the challenges of these problems is that the fluid motion, the flexible interface motion and the interaction with the immersed rigid boundaries must be computed simultaneously. This is necessary to account for the complex interaction between the fluid and the immersed boundaries. An example of interface problems that we consider is shown in Fig. 1.
In the immersed interface method, the flow is assumed to exist everywhere in the domain. The effect of the immersed boundaries on the fluid is accounted for through a singular force F which acts along the interface. The presence of this force results in solutions to the Navier-Stokes equations, which may be non-smooth across the interface, i.e., there may be jumps in pressure and in the derivatives of both pressure and velocity at the interface. An essential ingredient of the immersed interface method is the relation between the jumps in the solutions and their derivatives, and the applied singular forces. The basic idea of the immersed interface method is to discretize the Navier-Stokes equations on a uniform Cartesian grid and to account for the singular forces by explicitly incorporating the jumps in the solutions and their derivatives into the difference equations. The main advantage of the IIM is that the solutions of the Navier-Stokes equations on a uniform mesh can be done very efficiently with the use of fast solvers, and at the same time, complex geometrical changes can be handled in a rather seamless manner. The drawback of this method is that a special discretization of the Navier-Stokes equations near the immersed boundaries needs to be performed to maintain both accuracy and stability.In our algorithm, the immersed boundaries are represented by a number of Lagrangian control points. In order to guarantee that the no-slip condition on the rigid boundary is satisfied, singular forces are applied on the fluid. The force strength at the rigid boundary is determined by solving a small system of equations at each time step. For flexible boundaries, the forces that the boundary exerts on the fluid are computed from the constitutive relation of the flexible boundary and are applied to the fluid through the jump conditions. The position of the flexible boundary is updated implicitly using a quasi-Newton method (BFGS) within each timestep. The Navier-Stokes equations are discretized on a staggered Cartesian grid by a second order accurate projection method for pressure and velocity.
Our contributions include the introduction of a novel numerical algorithm for solving the incompressible Navier-Stokes equations in the presence of rigid boundaries. In our approach, both flexible and rigid boundaries can be considered simultaneously. We recall that most of the current methods can only handle flexible boundaries and the rigid boundary is required to be aligned with the computational grid.
Selected Publications and Presentations:
Rigid particles in a fluid under the action of gravity