A Barrier Method for Contact Avoiding Particles in Stokes Flow
Fuente:
arXiv
Guardado en:
| Autores principales: | , |
|---|---|
| Formato: | Preprint |
| Publicado: |
2023
|
| Materias: | |
| Acceso en línea: | |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866910478250278912 |
|---|---|
| author | Broms, Anna Tornberg, Anna-Karin |
| author_facet | Broms, Anna Tornberg, Anna-Karin |
| contents | Rigid particles in a Stokesian fluid can physically not overlap, as a thin layer of fluid always separates a particle pair, exerting increasingly strong repulsive forces on the bodies for decreasing separations. Numerically, resolving these lubrication forces comes at an intractably large cost even for moderate system sizes. Hence, it can typically not be guaranteed that particle collisions and overlaps do not occur in a dynamic simulation, independently of the choice of method to solve the Stokes equations. In this work, non-overlap constraints, in terms of the Euclidean distance between boundary points on the particles, are represented via a barrier energy. We solve for the minimum magnitudes of repelling contact forces between any particle pair in contact to correct for overlaps by enforcing a zero barrier energy at the next time level, given a contact-free configuration at a previous instance in time. The method is tested using a multiblob method to solve the mobility problem in Stokes flow applied to suspensions of spheres, rods and boomerang shaped particles. Collision free configurations are obtained at all instances in time. The effect of the contact forces on the collective order of a set of rods in a background flow that naturally promote particle interactions is also illustrated. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2301_01666 |
| institution | arXiv |
| publishDate | 2023 |
| record_format | arxiv |
| spellingShingle | A Barrier Method for Contact Avoiding Particles in Stokes Flow Broms, Anna Tornberg, Anna-Karin Fluid Dynamics Numerical Analysis 76D07 (Primary) 70F35, 90C33, 65Kxx (Secondary) Rigid particles in a Stokesian fluid can physically not overlap, as a thin layer of fluid always separates a particle pair, exerting increasingly strong repulsive forces on the bodies for decreasing separations. Numerically, resolving these lubrication forces comes at an intractably large cost even for moderate system sizes. Hence, it can typically not be guaranteed that particle collisions and overlaps do not occur in a dynamic simulation, independently of the choice of method to solve the Stokes equations. In this work, non-overlap constraints, in terms of the Euclidean distance between boundary points on the particles, are represented via a barrier energy. We solve for the minimum magnitudes of repelling contact forces between any particle pair in contact to correct for overlaps by enforcing a zero barrier energy at the next time level, given a contact-free configuration at a previous instance in time. The method is tested using a multiblob method to solve the mobility problem in Stokes flow applied to suspensions of spheres, rods and boomerang shaped particles. Collision free configurations are obtained at all instances in time. The effect of the contact forces on the collective order of a set of rods in a background flow that naturally promote particle interactions is also illustrated. |
| title | A Barrier Method for Contact Avoiding Particles in Stokes Flow |
| topic | Fluid Dynamics Numerical Analysis 76D07 (Primary) 70F35, 90C33, 65Kxx (Secondary) |
| url | https://arxiv.org/abs/2301.01666 |