Daozhi Han

Daozhi Han, Xiaoming Wang

We propose a novel second order in time numerical scheme for Cahn-Hilliard-Navier- Stokes phase field model with matched density. The scheme is based on second order convex-splitting for the Cahn-Hilliard equation and pressure-projection for the Navier-Stokes equation. We show that the scheme is mass-conservative, satisfies a modified energy law and is therefore unconditionally stable. Moreover, we prove that the scheme is uncondition- ally uniquely solvable at each time step by exploring the monotonicity associated with the scheme. Thanks to the weak coupling of the scheme, we design an efficient Picard iteration procedure to further decouple the computation of Cahn-Hilliard equation and Navier-Stokes equation. We implement the scheme by the mixed finite element method. Ample numerical experiments are performed to validate the accuracy and efficiency of the numerical scheme.

Daozhi Han, Yifeng Hou, Roger Temam

In this article,we study theoretically and numerically the interaction of a vortex induced by a rotating cylinder with a perpendicular plane. We show the existence of weak solutions to the swirling vortex models by using the Hopf extension method, and by an elegant contradiction argument, respectively. We demonstrate numerically that the model could produce phenomena of swirling vortex including boundary layer pumping and two-celled vortex that are observed in potential line vortex interacting with a plane and in a tornado.

Wenbin Chen, Daozhi Han, Xiaoming Wang

We propose and analyze two novel decoupled numerical schemes for solving the Cahn-Hilliard-Stokes-Darcy (CHSD) model for two-phase flows in karstic geometry. In the first numerical scheme, we explore a fractional step method (operator splitting) to decouple the phase-field (Cahn-Hilliard equation) from the velocity field (Stokes-Darcy fluid equations). To further decouple the Stokes-Darcy system, we introduce a first order pressure stabilization term in the Darcy solver in the second numerical scheme so that the Stokes system is decoupled from the Darcy system and hence the CHSD system can be solved in a fully decoupled manner. We show that both decoupled numerical schemes are uniquely solvable, energy stable, and mass conservative. Ample numerical results are presented to demonstrate the accuracy and efficiency of our schemes.

Gang Chen, Daozhi Han, Jiaxuan Liu et al.

We propose a hybridizable discontinuous Galerkin (HDG) method combined with convex-concave splitting for the temporal discretization of the convective Cahn-Hilliard equation. The convection term is discretized explicitly without stabilization, yielding three key advantages: (1) unconditional stability, (2) preservation of the optimal convergence rate for piecewise constant approximations, and (3) a symmetric system after local elimination, enabling efficient solver via minimal residual methods. We establish optimal convergence rates in the $L^2$ norm for both the scalar and flux variables for any polynomial degree $k \geq 0$. To achieve optimal $L^2$-norm estimates, we introduce a specialized HDG elliptic projection operator and analyze its approximation properties. Within the HDG framework, local elimination is employed to reduce the degrees of freedom associated with the globally coupled unknowns, and the scalar variables exhibit superconvergence. Finally, numerical experiments validate the theoretical convergence rates and demonstrate the effectiveness of the proposed method.

Daozhi Han

We propose a novel decoupled unconditionally stable numerical scheme for the simulation of two-phase flow in a Hele-Shaw cell which is governed by the Cahn-Hilliard-Hele-Shaw system (CHHS) with variable viscosity. The temporal discretization of the Cahn-Hilliard equation is based on a convex-splitting of the associated energy functional. Moreover, the capillary forcing term in the Darcy equation is separated from the pressure gradient at the time discrete level by using an operator-splitting strategy. Thus the computation of the nonlinear Cahn-Hilliard equation is completely decoupled from the update of pressure. Finally, a pressure-stabilization technique is used in the update of pressure so that at each time step one only needs to solve a Poisson equation with constant coefficient. We show that the scheme is unconditionally stable. Numerical results are presented to demonstrate the accuracy and efficiency of our scheme.