*Postgraduate, School of Naval Architecture and Maritime, Zhejiang Ocean University*

*Postgraduate, School of Naval Architecture and Maritime, Zhejiang Ocean University*

**3. Dynamic modelling of a flexible deep-water riser subjected to the internal slug-flow using 7-DOF ANCF combined with Arbitrary Lagrangian-Eulerian description (ALE-ANCF) (2020.11~2021.03)**

The two researches presented above is to understand the underlying principles of the finite element method and the implementation of FEM.

On the basis of Fortran based 6-DOF ANCF program, the **7-DOF model — ANCF combined with Arbitrary Lagrangian-Eulerian description (ALE-ANCF)** is proposed to investigate the dynamic response of the flexible deep-water riser subjected to the internal slug-flow.

Two cases are carried out to validate the implementation of the present model. A third case is carried to study the dynamic response of the steel catenary riser considering the internal fluid moving inside the riser.

First, a free-falling elastic beam under gravity in the air is set up to validate and demonstrate the implementation of the ALE-ANCF model. The initial setting and configuration of the elastic beam are the same as “2. Dynamic modelling of an elastic beam using 6-DOF Absolute Nodal Coordinate Formulation (2020.07~2020.10)“. Note that in this case, all the nodes’ nodal material coordinates have been **fixed**, which is done to ensure the results obtained from this ALE-ANCF model will keep consistent with conventional ANCF when mesh nodes don’t move in the simulation.

**Fig. 1** shows the impact of number of elements on the dynamic response of the cable. It can be viewed that the position of last node display no significant difference when the number of elements is more than 40. Thus, for the purpose of improving the computational efficiency and the accuracy, **40** elements are used to conduct the comparison with Berzeri and Shabana (2000)** **whose resluts are displayed in **Fig. 2**. The comparison results demonstrate our code results keep consistent with the results from the conventional ANCF method.

To validate the proposed ALE-ANCF model experiencing a complex external loads, the hydrodynamic drag force and buoyancy are considered. In this case, the free-falling of a rubber tether submerged in the still water is investigated. The evaluation of hydrodynamic drag force can refer to Takehara (2011)** **and Choo Y et al. (1971).

In conclusion, the implemented model ALE-ANCF model is acceptable through these two validation cases.

In this section, the dynamic behavior of a steel catenary riser (SCR) subjected to the internal upward slug-flow is investigated. For simplicity, only the partial riser is assumed to contain the inertial fluid — represented by the variable-length elements controlled by the moving material coordinates. The physical parameters of this scale model of SCR are taken from Morooka and Tsukada (2013). The schematic of SCR’s initial configuration is plotted in **Fig. 5** using 3 elements. The time step is set to be 0.0001 s, and the simulation is set as 60.0 s.

The stable dynamic configuration of SCR without the movement of internal slug-flow is firstly obtained to observe the riser’s variety of configurations subjected to the internal slug-flow. In this case, the first 30.0 s is set to obtain the stable dynamic configuration. After that, the internal slug-flow begin to move upward.

The comparisons of riser configuration are presented in **Fig. 6**. I can be viewed that the node 2 moves gradually to the upper left with node 3 shifts to the lower right with the movement of internal slug-flow.

To further reveal the impact of upward slug-flow on the structure of SCR, **Fig. 7** presents the displacement time history in the y-direction and z-direction of nodes 2 and 3.

Note that the average amplitude of SCR’s nodes with internal slug-flow is more significant than when there is no internal slug-flow movement inside the SCR.

Furthermore, node 2 gradually moves to the upper left under the influence of moving internal slug-flow in the global frame as shown in **(a)** and **(b)**, at the same time, node 3 shifts to the lower right can be seen in (c) and (d).

The reason for this phenomenon is that the movement of internal slug-flow controlled by the moving material coordinates corresponds to the introduction of external force to push the variable-length fluid element moving upward.

This phenomenon indicates that the internal slug-flow generates a superposition of the structure deformation, which aggravates the SCR oscillates.

- 79,880
- 38,951