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

Research output: Contribution to journalArticleResearchpeer review

Authors

  • Navid Valizadeh
  • Xiaoying Zhuang
  • Timon Rabczuk

Research Organisations

External Research Organisations

  • Bauhaus-Universität Weimar
View graph of relations

Details

Original languageEnglish
Article number117618
Number of pages23
JournalComputer Methods in Applied Mechanics and Engineering
Volume435
Early online date10 Dec 2024
Publication statusE-pub ahead of print - 10 Dec 2024

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

ASJC Scopus subject areas

Cite this

A monolithic finite element method for phase-field modeling of fully Eulerian fluid–structure interaction. / Valizadeh, Navid; Zhuang, Xiaoying; Rabczuk, Timon.
In: Computer Methods in Applied Mechanics and Engineering, Vol. 435, 117618, 15.02.2025.

Research output: Contribution to journalArticleResearchpeer review

Valizadeh N, Zhuang X, Rabczuk T. A monolithic finite element method for phase-field modeling of fully Eulerian fluid–structure interaction. Computer Methods in Applied Mechanics and Engineering. 2025 Feb 15;435:117618. Epub 2024 Dec 10. doi: 10.1016/j.cma.2024.117618
Download
@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 = "2024",
month = dec,
day = "10",
doi = "10.1016/j.cma.2024.117618",
language = "English",
volume = "435",
journal = "Computer Methods in Applied Mechanics and Engineering",
issn = "0045-7825",
publisher = "Elsevier",

}

Download

TY - JOUR

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

AU - Valizadeh, Navid

AU - Zhuang, Xiaoying

AU - Rabczuk, Timon

N1 - Publisher Copyright: © 2024 The Authors

PY - 2024/12/10

Y1 - 2024/12/10

N2 - 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.

AB - 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.

KW - Cahn–Hilliard phase-field model

KW - Finite element method

KW - Frictionless contact

KW - Fully Eulerian fluid–structure interaction

KW - Residual-based variational multiscale method

UR - http://www.scopus.com/inward/record.url?scp=85211218484&partnerID=8YFLogxK

U2 - 10.1016/j.cma.2024.117618

DO - 10.1016/j.cma.2024.117618

M3 - Article

AN - SCOPUS:85211218484

VL - 435

JO - Computer Methods in Applied Mechanics and Engineering

JF - Computer Methods in Applied Mechanics and Engineering

SN - 0045-7825

M1 - 117618

ER -