For the 3D-axisymmetric performance studies included in Section , a self-implemented staggered grid finite difference (SGFD) methodology [1]} is employed in the Lagrangian coordinate system (where \(\mathbf {s} = s \in \mathbb {R}\)  is the arc length), where only the tension force given by
\(\mathbf {T}=Eh\left(1-\sqrt{\frac{\partial \mathbf {X}}{\partial s}\cdot \frac{\partial \mathbf {X}}{\partial s}}\right)\) 
