139 research outputs found
A Second Order Fully-discrete Linear Energy Stable Scheme for a Binary Compressible Viscous Fluid Model
We present a linear, second order fully discrete numerical scheme on a
staggered grid for a thermodynamically consistent hydrodynamic phase field
model of binary compressible fluid flow mixtures derived from the generalized
Onsager Principle. The hydrodynamic model not only possesses the variational
structure, but also warrants the mass, linear momentum conservation as well as
energy dissipation. We first reformulate the model in an equivalent form using
the energy quadratization method and then discretize the reformulated model to
obtain a semi-discrete partial differential equation system using the
Crank-Nicolson method in time. The numerical scheme so derived preserves the
mass conservation and energy dissipation law at the semi-discrete level. Then,
we discretize the semi-discrete PDE system on a staggered grid in space to
arrive at a fully discrete scheme using the 2nd order finite difference method,
which respects a discrete energy dissipation law. We prove the unique
solvability of the linear system resulting from the fully discrete scheme. Mesh
refinements and two numerical examples on phase separation due to the spinodal
decomposition in two polymeric fluids and interface evolution in the gas-liquid
mixture are presented to show the convergence property and the usefulness of
the new scheme in applications
Positivity-preserving, energy stable numerical schemes for the Cahn-Hilliard equation with logarithmic potential
We present and analyze finite difference numerical schemes for the Allen
Cahn/Cahn-Hilliard equation with a logarithmic Flory Huggins energy potential.
Both the first order and second order accurate temporal algorithms are
considered. In the first order scheme, we treat the nonlinear logarithmic terms
and the surface diffusion term implicitly, and update the linear expansive term
and the mobility explicitly. We provide a theoretical justification that, this
numerical algorithm has a unique solution such that the positivity is always
preserved for the logarithmic arguments. In particular, our analysis reveals a
subtle fact: the singular nature of the logarithmic term around the values of
and 1 prevents the numerical solution reaching these singular values, so
that the numerical scheme is always well-defined as long as the numerical
solution stays similarly bounded at the previous time step. Furthermore, an
unconditional energy stability of the numerical scheme is derived, without any
restriction for the time step size. The unique solvability and the
positivity-preserving property for the second order scheme are proved using
similar ideas, in which the singular nature of the logarithmic term plays an
essential role. For both the first and second order accurate schemes, we are
able to derive an optimal rate convergence analysis, which gives the full order
error estimate. The case with a non-constant mobility is analyzed as well. We
also describe a practical and efficient multigrid solver for the proposed
numerical schemes, and present some numerical results, which demonstrate the
robustness of the numerical schemes
- …