An efficient phase-field method for turbulent multiphase flows
arXiv:2105.01865 · doi:10.1016/j.jcp.2021.110659
Abstract
With the aim of efficiently simulating three-dimensional multiphase turbulent flows with a phase-field method, we propose a new discretization scheme for the biharmonic term (the 4th-order derivative term) of the Cahn-Hilliard equation. This novel scheme can significantly reduce the computational cost while retaining the same accuracy as the original procedure. Our phase-field method is built on top of a direct numerical simulation solver, named AFiD (www.afid.eu) and open-sourced by our research group. It relies on a pencil distributed parallel strategy and a FFT-based Poisson solver. To deal with large density ratios between the two phases, a pressure split method [1] has been applied to the Poisson solver. To further reduce computational costs, we implement a multiple-resolution algorithm which decouples the discretizations for the Navier-Stokes equations and the scalar equation: while a stretched wall-resolving grid is used for the Navier-Stokes equations, for the Cahn-Hilliard equation we use a fine uniform mesh. The present method shows excellent computational performance for large-scale computation: on meshes up to 8 billion nodes and 3072 CPU cores, a multiphase flow needs only slightly less than 1.5 times the CPU time of the single-phase flow solver on the same grid. The present method is validated by comparing the results to previous studies for the cases of drop deformation in shear flow, including the convergence test with mesh refinement, and breakup of a rising buoyant bubble with density ratio up to 1000. Finally, we simulate the breakup of a big drop and the coalescence of O(10^3) drops in turbulent Rayleigh-Bénard convection at a Rayleigh number of , observing good agreement with theoretical results.
arXiv comment32 pages, 14 figures
References in corpus (16)
- Large scale dynamics in turbulent Rayleigh-Benard convection
- Two-phase binary fluids and immiscible fluids described by an order parameter
- Boundary layer structure in turbulent thermal convection and its consequences for the required numerical resolution
- Extended lifetime of respiratory droplets in a turbulent vapour puff and its implications on airborne disease transmission
- Transition to the ultimate regime in two-dimensional Rayleigh-Bénard convection
- Rayleigh-Bénard convection with a melting boundary
- AFiD-GPU: a versatile Navier-Stokes Solver for Wall-Bounded Turbulent Flows on GPU Clusters
- Transitions between turbulent states in rotating Rayleigh-Benard convection
- A multiple--resolution strategy for Direct Numerical Simulation of scalar turbulence
- Droplet-turbulence interactions and quasi-equilibrium dynamics in turbulent emulsions
- Droplets in homogeneous shear turbulence
- Flow organization and heat transfer in turbulent wall sheared thermal convection
- A parallel interaction potential approach coupled with the immersed boundary method for fully resolved simulations of deformable interfaces and membranes
- Wall roughness induces asymptotic ultimate turbulence
- Regime transitions in thermally driven high-Rayleigh number vertical convection
- Improved phase-field models of melting and dissolution in multi-component flows
Cited by in corpus (8)
- Heat transfer in turbulent Rayleigh-Bénard convection within two immiscible fluid layers
- Enhancing Heat Transport in Multiphase Rayleigh-Bénard Turbulence by Changing the Plate-Liquid Contact Angles
- Evaporating Rayleigh-Bénard convection: prediction of interface temperature and global heat transfer modulation
- Buoyancy-driven flow regimes for a melting vertical ice cylinder in saline water
- Turbulent Rayleigh-Bénard convection with bubbles attached to the plate
- A Graph Neural Network Simulation of Dispersed Systems
- Collective effects of neighbouring melting ice objects
- A FFT-accelerated multi-block finite-difference solver for massively parallel simulations of incompressible flows