A least squares based finite volume method for the Cahn-Hilliard and Cahn-Hilliard-reaction equations

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:273

Article

Dargaville, S.¹  Farrell, T. W.² 

[1] Univ London Imperial Coll Sci Technol & Med, Dept Earth Sci & Engn, Appl Modelling & Computat Grp, London SW7 2AZ, England

[2] Queensland Univ Technol, Brisbane, Qld 4001, Australia

关键词: Cahn-Hilliard;    Cahn-Hilliard-reaction;    Finite volume;    Least-squares;    LiFePO4;    Lithium-ion;   

DOI  :  10.1016/j.cam.2014.06.020

来源: Elsevier

【 摘 要 】

A vertex-centred finite volume method (FVM) for the Cahn-Hilliard (CH) and recently proposed Singh et al. (2008) and Burch (2009) Cahn-Hilliard-reaction (CHR) equations is presented. Information at control volume faces is computed using a high-order least-squares approach based on Taylor series approximations. This least-squares problem explicitly includes the variational boundary condition (VBC) that ensures that the discrete equations satisfy all of the boundary conditions. We use this approach to solve the CH and CHR equations in one and two dimensions and show that our scheme satisfies the VBC to at least second order. For the CH equation we show evidence of conservative, gradient stable solutions, however for the CHR equation, strict gradient-stability is more challenging to achieve. (C) 2014 Elsevier B.V. All rights reserved.

【 授权许可】

Free