Published October 2019 | Version v1
Journal article

A comparative study of interface-conforming ALE-FE scheme and diffuse interface AMR-LB scheme for interfacial dynamics

  • 1. Department of Mathematics, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon (Hong Kong)
  • 2. Beijing Computational Science Research Center, Beijing 100193 (China)
  • 3. Department of Civil & Environmental Engineering and Earth Sciences, University of Notre Dame, South Bend, IN 46556 (United States)
  • 4. BP Institute and Department of Engineering, University of Cambridge, Cambridge CB3 0EZ (United Kingdom)
  • 5. Department of Mathematics & Statistics, Old Dominion University, Norfolk, VA 23529 (United States)

Description

Highlights: • Comparing the spurious velocities computed by the ALE-FEM and AMR-LBM for a static bubble. • Simulations for the dynamics of suspended bubbles in the absence of gravity. • Simulations for the dynamics of rising bubbles driven by buoyancy force. • The Ohnesorge number determines whether the bubble dynamics is in the underdamped or the overdamped regime. • The ALE-FEM and AMR-LBM are compared on their numerical capability and efficiency. -- Abstract: In this work, a comparative study for two simulation methods is conducted for interfacial flows in two dimensions: an arbitrary Lagrangian Eulerian (ALE) finite element method (FEM) on interface-conforming meshes and a phase field lattice Boltzmann method (LBM) on Cartesian meshes with quadtree adaptive mesh refinement (AMR). The methods are validated by simulations of a bubble without and with buoyancy force. In particular, a suspended bubble with initial nonequilibrium shape and a rising bubble driven by buoyancy force are simulated to validate the methods. Additional simulations of the breakup of a rising bubble and the bubble interaction with a horizontal wall are used to quantify the efficacy and efficiency of the two methods. It is observed that the phase field LBM is more dissipative due to the nature of the diffuse interface method used to capture the interfaces. Overall, the results obtained from both methods agree well with each other when the effects due to the numerical artifacts intrinsic to the diffuse interface method can be neglected. Also, the LBM is in general more efficient and easier to be parallelized.

Availability note (English)

Available from http://dx.doi.org/10.1016/j.jcp.2019.06.048

Additional details

Identifiers

DOI
10.1016/j.jcp.2019.06.048;
PII
S0021999119304565;

Publishing Information

Journal Title
Journal of Computational Physics (Print)
Journal Volume
395
Journal Page Range
p. 602-619
ISSN
0021-9991
CODEN
JCTPAH

INIS

Country of Publication
Netherlands
Country of Input or Organization
International Atomic Energy Agency (IAEA)
INIS RN
54127135
Subject category
S71: CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS;
Descriptors DEI
COMPUTERIZED SIMULATION; FINITE ELEMENT METHOD; GRAVITATION; LAGRANGIAN FUNCTION
Descriptors DEC
CALCULATION METHODS; FUNCTIONS; MATHEMATICAL SOLUTIONS; NUMERICAL SOLUTION; SIMULATION

Optional Information

Copyright
Copyright (c) 2019 Elsevier Inc. All rights reserved.