Consequences of submarine landslides include both their direct impact on offshore infrastructure, such as subsea electric cables and gas/oil pipelines, and their indirect impact via the generated tsunami. The simulation of submarine landslides and their consequences has been a long-standing challenge majorly due to the strong coupling among sliding sediments, seawater and infrastructure as well as the induced extreme material deformation during the complete process. In this paper, we propose a unified finite element formulation for solid and fluid dynamics based on a generalised Hellinger–Reissner variational principle so that the coupling of fluid and solid can be achieved naturally in a monolithic fashion. In order to tackle extreme deformation problems, the resulting formulation is implemented within the framework of the particle finite element method. The correctness and robustness of the proposed unified formulation for single-phase problems (e.g. fluid dynamics problems involving Newtonian/Non-Newtonian flows and solid dynamics problems) as well as for multi-phase problems (e.g. two-phase flows) are verified against benchmarks. Comparisons are carried out against numerical and analytical solutions or experimental data that are available in the literature. Last but not least, the possibility of the proposed approach for modelling submarine landslides and their consequences is demonstrated via a numerical experiment of an underwater slope stability problem. It is shown that the failure and post-failure processes of the underwater slope can be predicted in a single simulation with its direct threat to a nearby pipeline and indirect threat by generating tsunami being estimated as well.
Abstract Consequences of submarine landslides include both their direct impact on offshore infrastructure, such as subsea electric cables and gas/oil pipelines, and their indirect [...]
Papers Repository of the International Centre for Numerical Methods in Engineering (CIMNE) (2019). 12
Abstract
This paper presents a purely Lagrangian approach for the 3D simulation of Bingham free-surface uids and their interaction with deformable solid structures.
In the proposed numerical strategy, the fluid is handled using the Particle Finite Element Method (PFEM) to tackle the issues resulting from extreme changes of geometry, such as mesh distortion and free-surface evolution. Additionally, the Papanastasiou model is employed as a regularization technique to overcome the
computational difficulties associated with the classical Bingham model. The solid structure, on the other hand, is represented by the hypoelastic constitutive model and simulated using the conventional Finite Element Method (FEM). The coupling between the fluid and the structure is achieved via a monolithic approach, called Unified formulation. Several numerical examples are presented to illustrate the correctness and the robustness of the proposed formulation, in 2D and in 3D. Special attention is devoted to the analysis of the convergence behavior of the proposed computational framework, the effect of the regularization on the numerical results and the 3D effects. Moreover, detailed comparisons between the simulated results and experimental data are performed so that the concerned problems and results can serve as benchmarks.
Abstract This paper presents a purely Lagrangian approach for the 3D simulation of Bingham free-surface uids and their interaction with deformable solid structures.
In the proposed [...]
Int. Journal for Numerical and Analytical Methods in Geomechanics (2018). Vol. 42, pp. 1806-1822 (preprint)
Abstract
This paper presents a detailed numerical study of the retrogressive failure of
landslides in sensitive clays. The dynamic modelling of the landslides is carried
out using a novel continuum approach, the particle finite element method,
complemented with an elastoviscoplastic constitutive model. The multiwedge
failure mode in the collapse is captured successfully, and the multiple retrogressive
failures that have been widely observed in landslides in sensitive clays
are reproduced with the failure mechanism, the kinematics, and the deposition
being discussed in detail. Special attention has been paid to the role of the clay
sensitivity on each retrogressive failure as well as on the final retrogression distance
and the final run‐ out distance via parametric studies. Moreover, the
effects of the viscosity of sensitive clays on the failure are also investigated
for different clay sensitivities.
Abstract This paper presents a detailed numerical study of the retrogressive failure of
landslides in sensitive clays. The dynamic modelling of the landslides is carried
out [...]