@@ -22,12 +22,12 @@ The liquid rises until it reaches a steady state corresponding to a straight flu
Gravity is neglected in this simulation. 4,000 time increments are performed with a time step of 2,00x10^-4 s. The liquid viscosity is 1.0 mPa.s, the air viscosity is 1.0 µPa.s, the liquid/air surface tension is 37.07 mN/m, the liquid/fiber energy is 30.03 mN/m, and the air/fiber surface energy is 60.9 mN/m. The normal velocity is imposed to be vanishing using an augmented Lagrangian technique at each node of the fiber boundaries. Additionally, the normal velocity is set to vanish (Dirichlet condition) at the left and right boundaries of the computational domain. Finally, a normal stress of 1 Pa is considered as Neumann's condition at the bottom boundary.
The mesh consists of 141,000 triangles and 71,200 nodes. The simulation ran for 4.8 hours on 16 cores using MUMPS (direct solver) for solving the Stokes system.
The mesh consists of 141,000 triangles and 71,200 nodes. The simulation ran for 4.8 hours on 16 cores using MUMPS (direct solver) for solving the Stokes system, and 3.5 hours with an iterative solver (Generalized Conjugate Residual method preconditioned by the Schur complement (mass matrix)).
Note that the pressure scale range is set to fixed values of [-2, +2] to avoid large variations in scale between increments. Additionally, the 'volume' of liquid is plotted against time. It evolves as t^0.45:
Additionally, the 'volume' of liquid is plotted against time. It evolves as t^0.4:
