Document Type

Article

Journal/Book Title/Conference

SIAM Journal on Scientific Computing

Volume

40

Issue

1

Publisher

Society for Industrial and Applied Mathematics

Publication Date

1-30-2018

First Page

B138

Last Page

B167

DOI

https://doi.org/10.1137/17M1111759

Abstract

We develop spatial-temporally second-order, energy stable numerical schemes for two classes of hydrodynamic phase field models of binary viscous fluid mixtures of different densities. One is quasi-incompressible while the other is incompressible. We introduce a novel energy quadratization technique to arrive at fully discrete linear schemes, where in each time step only a linear system needs to be solved. These schemes are then shown to be unconditionally energy stable rigorously subject to periodic boundary conditions so that a large time step is plausible. Both spatial and temporal mesh refinements are conducted to illustrate the second-order accuracy of the schemes. The linearization technique developed in this paper is so general that it can be applied to any thermodynamically consistent hydrodynamic theories so long as their energies are bounded below. Numerical examples on coarsening dynamics of two immiscible fluids and a heavy fluid drop settling in a lighter fluid matrix are presented to show the effectiveness of the proposed linear schemes. Predictions by the two fluid mixture models are compared and discussed, leading to our conclusion that the quasiincompressible model is more reliable than the incompressible one.

Included in

Mathematics Commons

Share

COinS