Finite element methods for particulate flows

Abstract

In Chapter 4 we described finite difference methods for the direct numerical simulations of fluid–solid systems in which the local flow field around moving particles is resolved numerically without modeling. Here we describe methods of the finite element type for the same purpose. The first method is termed the ALE (arbitrary Lagrangian-Eulerian) particle mover. The ALE particle mover uses a technique based on a combined formulation of the fluid and particle momentum equations, together with an arbitrary Lagrangian–Eulerian moving unstructured finite element mesh technique to deal with the movement of the particles. It was first developed by Hu et al . (1992, 1996, 2001). The method has been used to solve particle motions in both Newtonian and viscoelastic fluids in two- and three-dimensional flow geometries. It also handles particles of different sizes, shapes, and materials. The second method is based on a stabilized space-time finite element technique to solve problems involving moving boundaries and interfaces. In this method, the temporal coordinate is also discretized using finite elements. The deformation of the spatial domain with time is reflected simply in the deformation of the mesh in the temporal coordinate. Using this technique, Johnson and Tezduyar (1999) were able to simulate the sedimentation of 1000 spheres in a Newtonian fluid at a Reynolds number of 10. The third numerical method is termed the DLM (distributed Lagrange multiplier) particle mover. The basic idea of the DLM particle mover is to extend the problem from a time-dependent geometrically complex domain to a stationary fictitious domain, which is larger but simpler.