Numerical analysis of geometrical effect on drop deformation and breakup between rotating cylinders
Abstract
This study numerically investigates the effects of geometrical curvature on drop deformation and breakup in shear flow between counter-rotating concentric cylinders. Simulations were performed using a geometric Volume-of-Fluid (VOF) method incorporating plicRDF interface reconstruction and isoAdvector interface advection, combined with adaptive mesh refinement. Three circular geometries with different values of the non-dimensional geometrical curvature, κF, were considered and compared with a planar shear flow system. At fixed Reynolds and Capillary numbers (Re and Ca), the breakup morphology differed substantially between the planar and circular systems. In the planar system, the drop underwent relatively symmetric neck thinning, whereas in the circular systems, the drop exhibited asymmetric deformation and pronounced azimuthal elongation. At higher curvature, the drop formed a thin filament that subsequently fragmented, producing multiple satellite drops. Analysis of the internal velocity field and the magnitude of the rate-of-deformation tensor revealed distinct flow and deformation structures associated with these breakup behaviors. Furthermore, the critical Re for breakup, Recr, decreased with increasing κF at nvestigated Ca values. In contrast, the maximum dimensionless elongation length immediately before breakup at Re=Recr increased with κF. These findings suggest that geometrical curvature affects not only the breakup threshold but also the deformation process leading to breakup, providing insights into drop breakup in curved shear flows.
Related
- Critical deformation state of binary droplet breakup in turbulence
- Influence of nanoparticle shapes in nanofluid film boiling on vertical cylinders: A numerical study
- Metamorphosis of fractional instantons on a twisted T4 with a double-trace deformation: a numerical study
- Gravitational waves from core collapse of rotating very-massive stars: 3D numerical relativity computation
- Rotating Neutron Star Migrations as a Standardized Test for 3+1 Numerical Relativity
BibTeX
@misc{masuda2026numerical,
title = {{Numerical analysis of geometrical effect on drop deformation and breakup between rotating cylinders}},
author = {Hayato Masuda and Mitsuhiro Ohta and Tomislav Marić and Dieter Bothe},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.39154},
eprint = {2610.39154},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.39154},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - UNPB
TI - Numerical analysis of geometrical effect on drop deformation and breakup between rotating cylinders
AU - Hayato Masuda
AU - Mitsuhiro Ohta
AU - Tomislav Marić
AU - Dieter Bothe
PY - 2026
DA - 2026/10/29
PB - arXiv
DO - 10.48550/arxiv.2610.39154
UR - https://arxiv.org/abs/2610.39154
AB - This study numerically investigates the effects of geometrical curvature on drop deformation and breakup in shear flow between counter-rotating concentric cylinders. Simulations were performed using a geometric Volume-of-Fluid (VOF) method incorporating plicRDF interface reconstruction and isoAdvector interface advection, combined with adaptive mesh refinement. Three circular geometries with different values of the non-dimensional geometrical curvature, κ_F, were considered and compared with a planar shear flow system. At fixed Reynolds and Capillary numbers (Re and Ca), the breakup morphology differed substantially between the planar and circular systems. In the planar system, the drop underwent relatively symmetric neck thinning, whereas in the circular systems, the drop exhibited asymmetric deformation and pronounced azimuthal elongation. At higher curvature, the drop formed a thin filament that subsequently fragmented, producing multiple satellite drops. Analysis of the internal velocity field and the magnitude of the rate-of-deformation tensor revealed distinct flow and deformation structures associated with these breakup behaviors. Furthermore, the critical Re for breakup, Re_cr, decreased with increasing κ_F at nvestigated Ca values. In contrast, the maximum dimensionless elongation length immediately before breakup at Re=Re_cr increased with κ_F. These findings suggest that geometrical curvature affects not only the breakup threshold but also the deformation process leading to breakup, providing insights into drop breakup in curved shear flows.
ER -