A DLM immersed boundary method based wave-structure interaction solver for high density ratio multiphase flows

Nishant Nangia, Neelesh A Patankar*, Amneet Pal Singh Bhalla

*Corresponding author for this work

Research output: Contribution to journalArticle

Abstract

In this paper we present a robust immersed boundary (IB) method for high density ratio multiphase flows that is capable of modeling complex wave-structure interaction (WSI) problems arising in marine and coastal engineering applications. The IB/WSI methodology is enabled by combining the distributed Lagrange multiplier (DLM) method of Sharma and Patankar (2005) [13] with a robust level set method based multiphase flow solver. The fluid solver integrates the conservative form of the variable-coefficient incompressible Navier-Stokes equations using a hybrid preconditioner and ensures consistent transport of mass and momentum at a discrete level. The consistent transport scheme preserves the numerical stability of the method in the presence of large density ratios found in problems involving air, water, and an immersed structure. The air-water interface is captured by the level set method on an Eulerian grid, whereas the free-surface piercing immersed structure is represented on a Lagrangian mesh. The Lagrangian structure is free to move on the background Cartesian grid without conforming to the grid lines. The fluid-structure interaction (FSI) coupling is mediated via Peskin's regularized delta functions in an implicit manner, which obviates the need to integrate the hydrodynamic stress tensor on the complex surface of the immersed structure. The IB/WSI numerical scheme is implemented within an adaptive mesh refinement (AMR) framework, in which the Lagrangian structure and the air-water interface are embedded on the finest mesh level to capture the thin boundary layers and the vortical structures arising from WSI. We use a well balanced gravitational force discretization that eliminates spurious velocity currents in the hydrostatic limit due to density variation in the three phases (air, water and solid). We also show that using a non-conservative and an inconsistent fluid solver can lead to catastrophic failure of the numerical scheme for large density ratio variations that are prevalent in WSI applications. An effective wave generation and absorption technique for a numerical wave tank is presented and used to simulate a benchmark case of water wave distortion due to a submerged structure. The numerical scheme is tested on several benchmark WSI problems from numerical and experimental literature in both two and three dimensions to demonstrate the applicability of the IB/WSI method to practical marine and coastal engineering problems.

Original languageEnglish (US)
Article number108804
JournalJournal of Computational Physics
Volume398
DOIs
StatePublished - Dec 1 2019

Fingerprint

Immersed Boundary Method
Multiplier Method
Lagrange multipliers
multiphase flow
Multiphase Flow
Multiphase flow
Interaction
interactions
Marine engineering
Coastal engineering
air
grids
Immersed Boundary
water
Air
mesh
fluids
Water
Numerical Scheme
engineering

Keywords

  • Adaptive mesh refinement
  • Distributed Lagrange multipliers
  • Fictitious domain method
  • Fluid-structure interaction
  • Numerical wave tank
  • Stokes wave

ASJC Scopus subject areas

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • Physics and Astronomy(all)
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Cite this

@article{d3fc511a9d454beaafcdb4775cd1b86b,
title = "A DLM immersed boundary method based wave-structure interaction solver for high density ratio multiphase flows",
abstract = "In this paper we present a robust immersed boundary (IB) method for high density ratio multiphase flows that is capable of modeling complex wave-structure interaction (WSI) problems arising in marine and coastal engineering applications. The IB/WSI methodology is enabled by combining the distributed Lagrange multiplier (DLM) method of Sharma and Patankar (2005) [13] with a robust level set method based multiphase flow solver. The fluid solver integrates the conservative form of the variable-coefficient incompressible Navier-Stokes equations using a hybrid preconditioner and ensures consistent transport of mass and momentum at a discrete level. The consistent transport scheme preserves the numerical stability of the method in the presence of large density ratios found in problems involving air, water, and an immersed structure. The air-water interface is captured by the level set method on an Eulerian grid, whereas the free-surface piercing immersed structure is represented on a Lagrangian mesh. The Lagrangian structure is free to move on the background Cartesian grid without conforming to the grid lines. The fluid-structure interaction (FSI) coupling is mediated via Peskin's regularized delta functions in an implicit manner, which obviates the need to integrate the hydrodynamic stress tensor on the complex surface of the immersed structure. The IB/WSI numerical scheme is implemented within an adaptive mesh refinement (AMR) framework, in which the Lagrangian structure and the air-water interface are embedded on the finest mesh level to capture the thin boundary layers and the vortical structures arising from WSI. We use a well balanced gravitational force discretization that eliminates spurious velocity currents in the hydrostatic limit due to density variation in the three phases (air, water and solid). We also show that using a non-conservative and an inconsistent fluid solver can lead to catastrophic failure of the numerical scheme for large density ratio variations that are prevalent in WSI applications. An effective wave generation and absorption technique for a numerical wave tank is presented and used to simulate a benchmark case of water wave distortion due to a submerged structure. The numerical scheme is tested on several benchmark WSI problems from numerical and experimental literature in both two and three dimensions to demonstrate the applicability of the IB/WSI method to practical marine and coastal engineering problems.",
keywords = "Adaptive mesh refinement, Distributed Lagrange multipliers, Fictitious domain method, Fluid-structure interaction, Numerical wave tank, Stokes wave",
author = "Nishant Nangia and Patankar, {Neelesh A} and Bhalla, {Amneet Pal Singh}",
year = "2019",
month = "12",
day = "1",
doi = "10.1016/j.jcp.2019.07.004",
language = "English (US)",
volume = "398",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",

}

A DLM immersed boundary method based wave-structure interaction solver for high density ratio multiphase flows. / Nangia, Nishant; Patankar, Neelesh A; Bhalla, Amneet Pal Singh.

In: Journal of Computational Physics, Vol. 398, 108804, 01.12.2019.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A DLM immersed boundary method based wave-structure interaction solver for high density ratio multiphase flows

AU - Nangia, Nishant

AU - Patankar, Neelesh A

AU - Bhalla, Amneet Pal Singh

PY - 2019/12/1

Y1 - 2019/12/1

N2 - In this paper we present a robust immersed boundary (IB) method for high density ratio multiphase flows that is capable of modeling complex wave-structure interaction (WSI) problems arising in marine and coastal engineering applications. The IB/WSI methodology is enabled by combining the distributed Lagrange multiplier (DLM) method of Sharma and Patankar (2005) [13] with a robust level set method based multiphase flow solver. The fluid solver integrates the conservative form of the variable-coefficient incompressible Navier-Stokes equations using a hybrid preconditioner and ensures consistent transport of mass and momentum at a discrete level. The consistent transport scheme preserves the numerical stability of the method in the presence of large density ratios found in problems involving air, water, and an immersed structure. The air-water interface is captured by the level set method on an Eulerian grid, whereas the free-surface piercing immersed structure is represented on a Lagrangian mesh. The Lagrangian structure is free to move on the background Cartesian grid without conforming to the grid lines. The fluid-structure interaction (FSI) coupling is mediated via Peskin's regularized delta functions in an implicit manner, which obviates the need to integrate the hydrodynamic stress tensor on the complex surface of the immersed structure. The IB/WSI numerical scheme is implemented within an adaptive mesh refinement (AMR) framework, in which the Lagrangian structure and the air-water interface are embedded on the finest mesh level to capture the thin boundary layers and the vortical structures arising from WSI. We use a well balanced gravitational force discretization that eliminates spurious velocity currents in the hydrostatic limit due to density variation in the three phases (air, water and solid). We also show that using a non-conservative and an inconsistent fluid solver can lead to catastrophic failure of the numerical scheme for large density ratio variations that are prevalent in WSI applications. An effective wave generation and absorption technique for a numerical wave tank is presented and used to simulate a benchmark case of water wave distortion due to a submerged structure. The numerical scheme is tested on several benchmark WSI problems from numerical and experimental literature in both two and three dimensions to demonstrate the applicability of the IB/WSI method to practical marine and coastal engineering problems.

AB - In this paper we present a robust immersed boundary (IB) method for high density ratio multiphase flows that is capable of modeling complex wave-structure interaction (WSI) problems arising in marine and coastal engineering applications. The IB/WSI methodology is enabled by combining the distributed Lagrange multiplier (DLM) method of Sharma and Patankar (2005) [13] with a robust level set method based multiphase flow solver. The fluid solver integrates the conservative form of the variable-coefficient incompressible Navier-Stokes equations using a hybrid preconditioner and ensures consistent transport of mass and momentum at a discrete level. The consistent transport scheme preserves the numerical stability of the method in the presence of large density ratios found in problems involving air, water, and an immersed structure. The air-water interface is captured by the level set method on an Eulerian grid, whereas the free-surface piercing immersed structure is represented on a Lagrangian mesh. The Lagrangian structure is free to move on the background Cartesian grid without conforming to the grid lines. The fluid-structure interaction (FSI) coupling is mediated via Peskin's regularized delta functions in an implicit manner, which obviates the need to integrate the hydrodynamic stress tensor on the complex surface of the immersed structure. The IB/WSI numerical scheme is implemented within an adaptive mesh refinement (AMR) framework, in which the Lagrangian structure and the air-water interface are embedded on the finest mesh level to capture the thin boundary layers and the vortical structures arising from WSI. We use a well balanced gravitational force discretization that eliminates spurious velocity currents in the hydrostatic limit due to density variation in the three phases (air, water and solid). We also show that using a non-conservative and an inconsistent fluid solver can lead to catastrophic failure of the numerical scheme for large density ratio variations that are prevalent in WSI applications. An effective wave generation and absorption technique for a numerical wave tank is presented and used to simulate a benchmark case of water wave distortion due to a submerged structure. The numerical scheme is tested on several benchmark WSI problems from numerical and experimental literature in both two and three dimensions to demonstrate the applicability of the IB/WSI method to practical marine and coastal engineering problems.

KW - Adaptive mesh refinement

KW - Distributed Lagrange multipliers

KW - Fictitious domain method

KW - Fluid-structure interaction

KW - Numerical wave tank

KW - Stokes wave

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

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

U2 - 10.1016/j.jcp.2019.07.004

DO - 10.1016/j.jcp.2019.07.004

M3 - Article

AN - SCOPUS:85071291683

VL - 398

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

M1 - 108804

ER -