79 research outputs found
A numerical method for the quasi-incompressible Cahn-Hilliard-Navier-Stokes equations for variable density flows with a discrete energy law
In this paper, we investigate numerically a diffuse interface model for the
Navier-Stokes equation with fluid-fluid interface when the fluids have
different densities \cite{Lowengrub1998}. Under minor reformulation of the
system, we show that there is a continuous energy law underlying the system,
assuming that all variables have reasonable regularities. It is shown in the
literature that an energy law preserving method will perform better for
multiphase problems. Thus for the reformulated system, we design a finite
element method and a special temporal scheme where the energy law is preserved
at the discrete level. Such a discrete energy law (almost the same as the
continuous energy law) for this variable density two-phase flow model has never
been established before with finite element. A Newton's method is
introduced to linearise the highly non-linear system of our discretization
scheme. Some numerical experiments are carried out using the adaptive mesh to
investigate the scenario of coalescing and rising drops with differing density
ratio. The snapshots for the evolution of the interface together with the
adaptive mesh at different times are presented to show that the evolution,
including the break-up/pinch-off of the drop, can be handled smoothly by our
numerical scheme. The discrete energy functional for the system is examined to
show that the energy law at the discrete level is preserved by our scheme
Numerical Complete Solution for Random Genetic Drift by Energetic Variational Approach
In this paper, we focus on numerical solutions for random genetic drift
problem, which is governed by a degenerated convection-dominated parabolic
equation. Due to the fixation phenomenon of genes, Dirac delta singularities
will develop at boundary points as time evolves. Based on an energetic
variational approach (EnVarA), a balance between the maximal dissipation
principle (MDP) and least action principle (LAP), we obtain the trajectory
equation. In turn, a numerical scheme is proposed using a convex splitting
technique, with the unique solvability (on a convex set) and the energy decay
property (in time) justified at a theoretical level. Numerical examples are
presented for cases of pure drift and drift with semi-selection. The remarkable
advantage of this method is its ability to catch the Dirac delta singularity
close to machine precision over any equidistant grid.Comment: 22 pages, 11 figures, 2 table
- …