Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II
Abstract
Here, in this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes;more »
- Authors:
-
- Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
- Publication Date:
- Research Org.:
- Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); USDOE Laboratory Directed Research and Development (LDRD) Program
- OSTI Identifier:
- 1379922
- Alternate Identifier(s):
- OSTI ID: 1396594
- Grant/Contract Number:
- AC02-05CH11231
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Journal of Computational Physics
- Additional Journal Information:
- Journal Volume: 344; Journal Issue: C; Journal ID: ISSN 0021-9991
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; interface dynamics; discontinuous Galerkin methods; implicitly defined meshes; interfacial gauge methods; level set methods; high-order
Citation Formats
Saye, Robert. Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II. United States: N. p., 2017.
Web. doi:10.1016/j.jcp.2017.05.003.
Saye, Robert. Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II. United States. https://doi.org/10.1016/j.jcp.2017.05.003
Saye, Robert. Thu .
"Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II". United States. https://doi.org/10.1016/j.jcp.2017.05.003. https://www.osti.gov/servlets/purl/1379922.
@article{osti_1379922,
title = {Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II},
author = {Saye, Robert},
abstractNote = {Here, in this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.},
doi = {10.1016/j.jcp.2017.05.003},
journal = {Journal of Computational Physics},
number = C,
volume = 344,
place = {United States},
year = {Thu May 11 00:00:00 EDT 2017},
month = {Thu May 11 00:00:00 EDT 2017}
}
Web of Science
Works referenced in this record:
Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part I
journal, September 2017
- Saye, Robert
- Journal of Computational Physics, Vol. 344
Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems
journal, January 2002
- Arnold, Douglas N.; Brezzi, Franco; Cockburn, Bernardo
- SIAM Journal on Numerical Analysis, Vol. 39, Issue 5
Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi formulations
journal, November 1988
- Osher, Stanley; Sethian, James A.
- Journal of Computational Physics, Vol. 79, Issue 1
The Voronoi Implicit Interface Method for computing multiphase physics
journal, November 2011
- Saye, R. I.; Sethian, J. A.
- Proceedings of the National Academy of Sciences, Vol. 108, Issue 49
Analysis and applications of the Voronoi Implicit Interface Method
journal, July 2012
- Saye, R. I.; Sethian, J. A.
- Journal of Computational Physics, Vol. 231, Issue 18
Interfacial gauge methods for incompressible fluid dynamics
journal, June 2016
- Saye, Robert
- Science Advances, Vol. 2, Issue 6
Numerical solution of the Navier-Stokes equations
journal, January 1968
- Chorin, Alexandre Joel
- Mathematics of Computation, Vol. 22, Issue 104
On a new way of writing the Navier-Stokes equation. The Hamiltonian formalism
journal, June 1989
- Oseledets, V. I.
- Russian Mathematical Surveys, Vol. 44, Issue 3
A Hamiltonian theory for weakly interacting vortices
journal, December 1972
- Roberts, P. H.
- Mathematika, Vol. 19, Issue 2
On the Use of Higher-Order Projection Methods for Incompressible Turbulent Flow
journal, January 2013
- Almgren, A. S.; Aspden, A. J.; Bell, J. B.
- SIAM Journal on Scientific Computing, Vol. 35, Issue 1
A fourth-order auxiliary variable projection method for zero-Mach number gas dynamics
journal, January 2008
- Kadioglu, Samet Y.; Klein, Rupert; Minion, Michael L.
- Journal of Computational Physics, Vol. 227, Issue 3
Semi-implicit projection methods for incompressible flow based on spectral deferred corrections
journal, March 2004
- Minion, Michael L.
- Applied Numerical Mathematics, Vol. 48, Issue 3-4
An overview of projection methods for incompressible flows
journal, September 2006
- Guermond, J. L.; Minev, P.; Shen, Jie
- Computer Methods in Applied Mechanics and Engineering, Vol. 195, Issue 44-47
Stable and accurate pressure approximation for unsteady incompressible viscous flow
journal, May 2010
- Liu, Jian-Guo; Liu, Jie; Pego, Robert L.
- Journal of Computational Physics, Vol. 229, Issue 9
Application of a fractional-step method to incompressible Navier-Stokes equations
journal, June 1985
- Kim, J.; Moin, P.
- Journal of Computational Physics, Vol. 59, Issue 2
Accurate Projection Methods for the Incompressible Navier–Stokes Equations
journal, April 2001
- Brown, David L.; Cortez, Ricardo; Minion, Michael L.
- Journal of Computational Physics, Vol. 168, Issue 2
High-order methods for computing distances to implicitly defined surfaces
journal, January 2014
- Saye, Robert
- Communications in Applied Mathematics and Computational Science, Vol. 9, Issue 1
A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework
journal, March 2006
- Francois, Marianne M.; Cummins, Sharen J.; Dendy, Edward D.
- Journal of Computational Physics, Vol. 213, Issue 1
A continuum method for modeling surface tension
journal, June 1992
- Brackbill, J. U.; Kothe, D. B.; Zemach, C.
- Journal of Computational Physics, Vol. 100, Issue 2
Oscillations of soap bubbles
journal, July 2010
- Kornek, U.; Müller, F.; Harth, K.
- New Journal of Physics, Vol. 12, Issue 7
Fluid pipes
journal, September 2002
- Hancock, Matthew J.; Bush, John W. M.
- Journal of Fluid Mechanics, Vol. 466
A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs
journal, January 1998
- Karypis, George; Kumar, Vipin
- SIAM Journal on Scientific Computing, Vol. 20, Issue 1
High-Order Quadrature Methods for Implicitly Defined Surfaces and Volumes in Hyperrectangles
journal, January 2015
- Saye, R. I.
- SIAM Journal on Scientific Computing, Vol. 37, Issue 2
Patch-Recovery Filters for Curvature in Discontinuous Galerkin-Based Level-Set Methods
journal, February 2016
- Kummer, Florian; Warburton, Tim
- Communications in Computational Physics, Vol. 19, Issue 2
Works referencing / citing this record:
Multimesh finite element methods: Solving PDEs on multiple intersecting meshes
journal, January 2019
- Johansson, August; Kehlet, Benjamin; Larson, Mats G.
- Computer Methods in Applied Mechanics and Engineering, Vol. 343
Efficient operator-coarsening multigrid schemes for local discontinuous Galerkin methods
preprint, January 2018
- Fortunato, Daniel; Rycroft, Chris H.; Saye, Robert
- arXiv