We consider suspensions of rigid bodies in a two‐dimensional viscous fluid. Even with high‐fidelity numerical methods, unphysical contact between particles occurs because of spatial and temporal discretization errors. We extend a time stepping met...
We consider suspensions of rigid bodies in a two‐dimensional viscous fluid. Even with high‐fidelity numerical methods, unphysical contact between particles occurs because of spatial and temporal discretization errors. We extend a time stepping method that avoids overlap by imposing a minimum separation distance between all pairs of bodies. In its original form, the method discretizes interactions between different particles explicitly. Therefore, to avoid stiffness, a large minimum separation distance is used. In this paper, we introduce a new implicit time stepping method that is able to simulate dense suspensions with large time step sizes and a small minimum separation distance. The method is tested on various unbounded and bounded flows, and rheological properties of the resulting suspensions are computed.
We develop a time stepping method to simulate rigid body suspensions in a viscous fluid. By combining a contact algorithm with implicit interactions, suspensions with nearly‐touching bodies are simulated. We analyze the convergence of the method and use the simulations to investigate rheological properties of the flow.