Saved in:
| Main Authors: | , , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | https://arxiv.org/abs/2508.19700 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866912556224872448 |
|---|---|
| author | Benkley, Tyler Deparis, Simone Ricci, Paolo Mortensen, Andreas |
| author_facet | Benkley, Tyler Deparis, Simone Ricci, Paolo Mortensen, Andreas |
| contents | An iterative Finite Element method predicated on a linearisation of the weak form around a reference configuration is derived for general, three-dimensional, free-surface flows, including systems with moving contact lines. The method is a rigorous generalisation of the Total Linearisation Method that was proposed by Kruyt et al. (1988) for two-dimensional flows with contact angles limited to $90^\circ$. In contrast to existing numerical methods for free-surface-flow problems, the present linearisation produces a weak form that is devoid of displacement degrees of freedom in the bulk, thus nearly halving the size of the linear system when compared to standard linearised methods. A novel preconditioner, whose implementation is made possible by the size reduction, is employed to solve the large resulting monolithic Jacobian systems with the Generalised Minimum Residual Method. The proposed method and the preconditioner are shown to be effective on two numerical examples of capillary flows, namely (i) the cylindrical die swell problem, solved both in 3D Cartesian coordinates and under the Ansatz of axisymmetry; and (ii) an enclosed 2D thermo-capillary problem. For the die swell, numerical results are validated by existing experimental results and prior simulations, and confirm both the extension to 3D and theoretical convergence rates. For the thermo-capillary problem, simulations verify earlier calculations. Additional simulations are also carried out for a new range of contact angles made possible by the extension. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2508_19700 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Generalisation of the Total Linearisation Method to Three-dimensional Free-Surface Flows Benkley, Tyler Deparis, Simone Ricci, Paolo Mortensen, Andreas Fluid Dynamics An iterative Finite Element method predicated on a linearisation of the weak form around a reference configuration is derived for general, three-dimensional, free-surface flows, including systems with moving contact lines. The method is a rigorous generalisation of the Total Linearisation Method that was proposed by Kruyt et al. (1988) for two-dimensional flows with contact angles limited to $90^\circ$. In contrast to existing numerical methods for free-surface-flow problems, the present linearisation produces a weak form that is devoid of displacement degrees of freedom in the bulk, thus nearly halving the size of the linear system when compared to standard linearised methods. A novel preconditioner, whose implementation is made possible by the size reduction, is employed to solve the large resulting monolithic Jacobian systems with the Generalised Minimum Residual Method. The proposed method and the preconditioner are shown to be effective on two numerical examples of capillary flows, namely (i) the cylindrical die swell problem, solved both in 3D Cartesian coordinates and under the Ansatz of axisymmetry; and (ii) an enclosed 2D thermo-capillary problem. For the die swell, numerical results are validated by existing experimental results and prior simulations, and confirm both the extension to 3D and theoretical convergence rates. For the thermo-capillary problem, simulations verify earlier calculations. Additional simulations are also carried out for a new range of contact angles made possible by the extension. |
| title | Generalisation of the Total Linearisation Method to Three-dimensional Free-Surface Flows |
| topic | Fluid Dynamics |
| url | https://arxiv.org/abs/2508.19700 |