A monolithic finite element method for phase-field modeling of fully Eulerian fluid–structure interaction

authored by: Navid Valizadeh, Xiaoying Zhuang, Timon Rabczuk
Abstract: In this paper, we introduce a fully-monolithic, implicit finite element method designed for investigating fluid–structure interaction problems within a fully Eulerian framework. Our approach employs a coupled Navier–Stokes Cahn–Hilliard phase-field model, recently developed by Mokbel et al. (2018). This model adeptly addresses significant challenges such as large solid deformations, topology changes, stable incorporation of surface tensions, and eliminates the need for remeshing methods. While the original model was primarily tested for axisymmetric problems, our work extends its application to encompass a range of two- and three-dimensional verification tests. Additionally, we advance the model to handle multi-solid–fluid interaction scenarios through the integration of a multi-body contact algorithm. Assuming both the solid and fluid to be incompressible, we describe them using Navier–Stokes equations. For the solid, a hyperelastic neo-Hookean material is assumed, and the elastic solid stress is computed based on the left Cauchy–Green deformation tensor, which is governed by an Oldroyd-B like equation. We employ a residual-based variational multiscale method for solving the full Navier–Stokes equations, a stabilized Galerkin finite element method using Streamline-Upwind/Petrov–Galerkin (SUPG) stabilization for solving the Oldroyd-B equation, and a mixed finite element splitting scheme for the Cahn–Hilliard equation. The system of partial differential equations is solved using an implicit, monolithic scheme based on the generalized-α time integration method. Our approach is verified through two-dimensional numerical examples, including the deformation of an elastic wall by flow, the deformation and motion of a solid disk in a lid-driven cavity flow, and the bouncing of an elastic ball, showcasing the method's ability to handle solid-wall contact. Furthermore, we extend the application to multi-body contact problems and verify the model's accuracy by solving three-dimensional benchmark problems, such as the motion of an elastic solid sphere in lid-driven cavity flow and the falling of an elastic sphere onto an elastic block.
Organisation(s): Institute of Photonics
External Organisation(s): Bauhaus-Universität Weimar
Type: Article
Journal: Computer Methods in Applied Mechanics and Engineering
Volume: 435
No. of pages: 23
ISSN: 0045-7825
Publication date: 15.02.2025
Publication status: Published
Peer reviewed: Yes
ASJC Scopus subject areas: Computational Mechanics, Mechanics of Materials, Mechanical Engineering, General Physics and Astronomy, Computer Science Applications
Electronic version(s): https://doi.org/10.1016/j.cma.2024.117618 (Access: Open)
: Details in the research portal "Research@Leibniz University"

BibTeX

@article{dff4d62fb7d84cd0807e2f24e90e043b,
title = "A monolithic finite element method for phase-field modeling of fully Eulerian fluid–structure interaction",
abstract = "In this paper, we introduce a fully-monolithic, implicit finite element method designed for investigating fluid–structure interaction problems within a fully Eulerian framework. Our approach employs a coupled Navier–Stokes Cahn–Hilliard phase-field model, recently developed by Mokbel et al. (2018). This model adeptly addresses significant challenges such as large solid deformations, topology changes, stable incorporation of surface tensions, and eliminates the need for remeshing methods. While the original model was primarily tested for axisymmetric problems, our work extends its application to encompass a range of two- and three-dimensional verification tests. Additionally, we advance the model to handle multi-solid–fluid interaction scenarios through the integration of a multi-body contact algorithm. Assuming both the solid and fluid to be incompressible, we describe them using Navier–Stokes equations. For the solid, a hyperelastic neo-Hookean material is assumed, and the elastic solid stress is computed based on the left Cauchy–Green deformation tensor, which is governed by an Oldroyd-B like equation. We employ a residual-based variational multiscale method for solving the full Navier–Stokes equations, a stabilized Galerkin finite element method using Streamline-Upwind/Petrov–Galerkin (SUPG) stabilization for solving the Oldroyd-B equation, and a mixed finite element splitting scheme for the Cahn–Hilliard equation. The system of partial differential equations is solved using an implicit, monolithic scheme based on the generalized-α time integration method. Our approach is verified through two-dimensional numerical examples, including the deformation of an elastic wall by flow, the deformation and motion of a solid disk in a lid-driven cavity flow, and the bouncing of an elastic ball, showcasing the method's ability to handle solid-wall contact. Furthermore, we extend the application to multi-body contact problems and verify the model's accuracy by solving three-dimensional benchmark problems, such as the motion of an elastic solid sphere in lid-driven cavity flow and the falling of an elastic sphere onto an elastic block.",
keywords = "Cahn–Hilliard phase-field model, Finite element method, Frictionless contact, Fully Eulerian fluid–structure interaction, Residual-based variational multiscale method",
author = "Navid Valizadeh and Xiaoying Zhuang and Timon Rabczuk",
note = "Publisher Copyright: {\textcopyright} 2024 The Authors",
year = "2025",
month = feb,
day = "15",
doi = "10.1016/j.cma.2024.117618",
language = "English",
volume = "435",
journal = "Computer Methods in Applied Mechanics and Engineering",
issn = "0045-7825",
publisher = "Elsevier BV",
}

Back