Analytical solutions of layered Poiseuille flows in the diffuse interface model

TL;DR

This paper derives analytical solutions for layered Poiseuille flows in the diffuse interface model, compares them with sharp interface models, and proposes new viscosity models to improve simulation accuracy.

Contribution

It provides analytical solutions for steady layered Poiseuille flows in the diffuse interface model and introduces new viscosity models to enhance simulation fidelity.

Findings

Analytical solutions differ significantly from sharp interface models at high viscosity ratios.

Numerical results agree well with analytical solutions, validating the theoretical approach.

Proposed symmetrical viscosity models improve the accuracy of diffuse interface simulations.

Abstract

Based on the two-phase macroscopic governing equations in the phase field model, the governing equations and analytical solutions for the steady-state layered Poiseuille flows in the diffuse interface (DI) model are derived and analyzed. Then, based on three dynamic viscosity models commonly used in the literature, the corresponding analytical solutions of the velocity profile are obtained. Under the condition of high dynamic viscosity ratio, the analytical solution of DI model may be significantly different from that of the sharp interface (SI) model, and the degree of deviation depends on the dynamic viscosity model and the interface thickness. Therefore, the numerical simulation of layered Poiseuille flow with DI model should be compared with the analytical solution of DI model with the same dynamic viscosity model. A direct comparison of the numerical solution results with the SI analytical solution could misinterpret the model error with the numerical error. In addition, the direct numerical simulation data and the DI analytical solutions agree well, which validates the theoretical results. Finally, a new set of symmetrical dynamic viscosity models is proposed and recommended for the simulation of two-phase flows in the DI model, which makes both the viscosity profiles and velocity profiles close to the SI model.