Numerical simulation of sedimentation of rectangular particle in Newtonian fluid

The sedimentation of a rectangular particle falling in a two-dimensional channel filled with Newtonian fluid was simulated with finite element arbitrary Lagrangian–Eulerian domain method. The numerical procedure was validated by comparison of the simulation results with existing numerical work. Moreover, good agreement was obtained between the simulation results and experimental measurements performed in the current study. The equilibrium position, stable orientation and drag coefficient of a rectangular particle for different particle Reynolds numbers (Re_p) were studied. The results show that there is a critical particle Reynolds number for the preferred orientation of a rectangular particle falling in a Newtonian fluid. When Re_p is smaller than the critical value, the particle falls with its long side parallel to gravity; otherwise the particle falls with its long side perpendicular to gravity. The critical particle Reynolds number is a decreasing function of the blockage ratio and aspect ratio. The distributions of pressure and shear stress on rectangular particle surface were analyzed. Moreover, the drag coefficient of the rectangular particle decreases as Re_p or the blockage ratio increases; however, it appears to be independent of aspect ratio.