Gorkem Simsek

Mahnaz Shokrpour Roudbari, Gorkem Simsek, E. Harald van Brummelen et al.

While various phase-field models have recently appeared for two-phase fluids with different densities, only some are known to be thermodynamically consistent, and practical stable schemes for their numerical simulation are lacking. In this paper, we derive a new form of thermodynamically-consistent quasi-incompressible diffuse-interface Navier-Stokes Cahn-Hilliard model for a two-phase flow of incompressible fluids with different densities. The derivation is based on mixture theory by invoking the second law of thermodynamics and Coleman-Noll procedure. We also demonstrate that our model and some of the existing models are equivalent and we provide a unification between them. In addition, we develop a linear and energy-stable time-integration scheme for the derived model. Such a linearly-implicit scheme is nontrivial, because it has to suitably deal with all nonlinear terms, in particular those involving the density. Our proposed scheme is the first linear method for quasi-incompressible two-phase flows with nonsolenoidal velocity that satisfies discrete energy dissipation independent of the time-step size, provided that the mixture density remains positive. The scheme also preserves mass. Numerical experiments verify the suitability of the scheme for two-phase flow applications with high density ratios using large time steps by considering the coalescence and break-up dynamics of droplets including pinching due to gravity.

Xunxun Wu, Kristoffer van der Zee, Gorkem Simsek et al.

While many methods exist to discretize nonlinear time-dependent partial differential equations (PDEs), the rigorous estimation and adaptive control of their discretization errors remains challenging. In this paper, we present a methodology for duality-based a posteriori error estimation for nonlinear parabolic PDEs, where the full discretization of the PDE relies on the use of an implicit-explicit (IMEX) time-stepping scheme and the finite element method in space. The main result in our work is a decomposition of the error estimate that allows to separate the effects of spatial and temporal discretization error, and which can be used to drive adaptive mesh refinement and adaptive time-step selection. The decomposition hinges on a specially-tailored IMEX discretization of the dual problem. The performance of the error estimates and the proposed adaptive algorithm is demonstrated on two canonical applications: the elementary heat equation and the nonlinear Allen-Cahn phase-field model.