Abstract
In this work, we implemented and compared two different methods to impose the rigidβbody motion constraint on a solid particle moving inside a fluid. We consider a fictitious domain method to easily manage the particle motion. As the solid as well as the fluid inertia are neglected, the particle can be discretized through its boundary only. The rigidβbody motion is imposed via Lagrange multipliers on the boundary. In the first method, such constraints are imposed in discrete points on the boundary (collocation), whereas in the second the constraint is imposed in a weak way on elements dividing the particle surface.
Two test problems, that is, a spherical and an ellipsoidal particle in a sheared Newtonian fluid, are chosen to compare the methods. In both cases, the analysis is carried out in 2D as well as in 3D.
The results show that for the collocation method an optimal number of collocation points exist leading to the smallest error. However, small variations in the optimal value can generate large deviations. In the weak implementation, the error is only mildly affected by the number of elements used to discretize the particle boundary and by the Lagrange multiplier's interpolation space. A further analysis is carried out to study the effect of an approximated integration of weak constraints.
A comparison between the two methods showed that the same accuracy can be achieved by using less constraints if the weak discretization is used. Finally, the rigidβbody motion imposed via weak constraints leads to better conditioned linear systems. Copyright Β© 2009 John Wiley & Sons, Ltd.