Energy stable multigrid method for local and non-local hydrodynamic models for freezing

Abstract In this paper we present a numerical method for hydrodynamic models that arise from time dependent density functional theories of freezing. The models take the form of compressible Navier–Stokes equations whose pressure is determined by the variational derivative of a free energy, which is a functional of the density field. We present unconditionally energy stable and mass conserving implicit finite difference methods for the models. The methods are based on a convex splitting of the free energy and that ensures that a discrete energy is non-increasing for any choice of time and space step. The methods are applicable to a large class of models, including both local and non-local free energy functionals. The theoretical basis for the numerical method is presented in a general context. The method is applied to problems using two specific free energy functionals: one local and one non-local functional. A nonlinear multigrid method is used to solve the numerical method, which is nonlinear at the implicit time step. The non-local functional, which is a convolution operator, is approximated using the Discrete Fourier Transform. Numerical simulations that confirm the stability and accuracy of the numerical method are presented.

[1]  T. Ala‐Nissila,et al.  Extended phase diagram of the three-dimensional phase field crystal model , 2010, Journal of physics. Condensed matter : an Institute of Physics journal.

[2]  Cheng Wang,et al.  Energy stable and efficient finite-difference nonlinear multigrid schemes for the modified phase field crystal equation , 2012, J. Comput. Phys..

[3]  A. Voigt,et al.  A Navier-Stokes phase-field crystal model for colloidal suspensions. , 2013, The Journal of chemical physics.

[4]  Cheng Wang,et al.  A convergent convex splitting scheme for the periodic nonlocal Cahn-Hilliard equation , 2014, Numerische Mathematik.

[5]  James F Lutsko,et al.  A dynamical theory of nucleation for colloids and macromolecules. , 2011, The Journal of chemical physics.

[6]  A. Archer,et al.  Dynamical density functional theory for molecular and colloidal fluids: a microscopic approach to fluid mechanics. , 2008, The Journal of chemical physics.

[7]  M. Yussouff,et al.  First-principles order-parameter theory of freezing , 1979 .

[8]  Hartmut Löwen,et al.  Derivation of the phase-field-crystal model for colloidal solidification. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[9]  Cheng Wang,et al.  Stable and efficient finite-difference nonlinear-multigrid schemes for the phase field crystal equation , 2009, J. Comput. Phys..

[10]  Nikolas Provatas,et al.  Phase-Field Methods in Materials Science and Engineering , 2010 .

[11]  Cheng Wang,et al.  Second order convex splitting schemes for periodic nonlocal Cahn-Hilliard and Allen-Cahn equations , 2014, J. Comput. Phys..

[12]  Steven M. Wise,et al.  An Energy Stable and Convergent Finite-Difference Scheme for the Modified Phase Field Crystal Equation , 2011, SIAM J. Numer. Anal..

[13]  M. Grant,et al.  Modeling elastic and plastic deformations in nonequilibrium processing using phase field crystals. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  Steven M. Wise,et al.  Convergence Analysis of a Second Order Convex Splitting Scheme for the Modified Phase Field Crystal Equation , 2012, SIAM J. Numer. Anal..

[15]  P. Español,et al.  Derivation of dynamical density functional theory using the projection operator technique. , 2009, The Journal of chemical physics.

[16]  A. Yoshimori.  Microscopic derivation of time-dependent density functional methods. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[17]  P. Chavanis.  Brownian particles with long- and short-range interactions , 2010, 1008.2898.

[18]  Junseok Kim.  Phase-Field Models for Multi-Component Fluid Flows , 2012 .

[19]  James F. Lutsko,et al.  Recent Developments in Classical Density Functional Theory , 2010 .

[20]  Cheng Wang,et al.  An Energy-Stable and Convergent Finite-Difference Scheme for the Phase Field Crystal Equation , 2009, SIAM J. Numer. Anal..

[21]  D. J. Eyre.  Unconditionally Gradient Stable Time Marching the Cahn-Hilliard Equation , 1998 .

[22]  Pedro Tarazona,et al.  Dynamic density functional theory of fluids , 1999 .

[23]  J. Lowengrub,et al.  Kinetic density functional theory of freezing. , 2013, The Journal of chemical physics.

[24]  Peter Smereka,et al.  Kinetic Monte Carlo simulation of strained heteroepitaxial growth with intermixing , 2009, 0903.4882.

[25]  Grigorios A. Pavliotis,et al.  The Overdamped Limit of Dynamic Density Functional Theory: Rigorous Results , 2012, Multiscale Model. Simul..

[26]  L. Gránásy,et al.  Nonlinear hydrodynamic theory of crystallization , 2013, Journal of physics. Condensed matter : an Institute of Physics journal.