Fully-Decoupled and Second-Order Time-Accurate Scheme for the Cahn-Hilliard Ohta-Kawaski Phase-Field Model of Diblock Copolymer Melt Confined in Hele-Shaw Cell

被引:2
作者
Cao, Junying [1 ]
Zhang, Jun [2 ]
Yang, Xiaofeng [3 ]
机构
[1] Guizhou Minzu Univ, Sch Data Sci & Informat Engn, Guiyang 550025, Peoples R China
[2] Guizhou Univ Finance & Econ, Computat Math Res Ctr, Guizhou Key Lab Big Data Stat Anal, Guiyang 550025, Peoples R China
[3] Univ South Carolina, Dept Math, Columbia, SC 29208 USA
基金
美国国家科学基金会;
关键词
Phase-field; Explicit-IEQ; Darcy; Decoupled; Energy Stability; Diblock copolymer melt; STABLE NUMERICAL SCHEME; CONTACT LINE MODEL; MICROPHASE SEPARATION; EFFICIENT; APPROXIMATIONS; ALGORITHMS; DENSITIES; FINGERS; GROWTH; FLOWS;
D O I
10.1007/s40304-022-00298-3
中图分类号
O1 [数学];
学科分类号
0701 ; 070101 ;
摘要
In this work, we consider numerical approximations of the phase-field model of diblock copolymer melt confined in Hele-Shaw cell, which is a very complicated coupled nonlinear system consisting of the Darcy equations and the Cahn-Hilliard type equations with the Ohta-Kawaski potential. Through the combination of a novel explicit-Invariant Energy Quadratization approach and the projection method, we develop the first full decoupling, energy stable, and second-order time-accurate numerical scheme. The introduction of two auxiliary variables and the design of two auxiliary ODEs play a vital role in obtaining the full decoupling structure while maintaining energy stability. The scheme is also linear and unconditional energy stable, and the practical implementation efficiency is also very high because it only needs to solve a few elliptic equations with constant coefficients at each time step. We strictly prove that the scheme satisfies the unconditional energy stability and give a detailed implementation process. Numerical experiments further verify the convergence rate, energy stability, and effectiveness of the developed algorithm.
引用
收藏
页码:479 / 504
页数:26
相关论文
共 63 条
[1]   Low viscosity contrast fingering in a rotating Hele-Shaw cell [J].
Alvarez-Lacalle, E ;
Ortín, J ;
Casademunt, J .
PHYSICS OF FLUIDS, 2004, 16 (04) :908-924
[2]  
Amaya A, 2017, IEEE ICCE
[3]   An island of stability in a sea of fingers: emergent global features of the viscous-flow instability [J].
Bischofberger, Irmgard ;
Ramachandran, Radha ;
Nagel, Sidney R. .
SOFT MATTER, 2015, 11 (37) :7428-7432
[4]  
Brazovskii S.A., 1975, Sov. Phys. JETP, V41, P85, DOI DOI 10.1142/9789814317344_0016
[5]   VISCOUS FINGERING IN LIQUID-CRYSTALS [J].
BUKA, A ;
PALFFYMUHORAY, P ;
RACZ, Z .
PHYSICAL REVIEW A, 1987, 36 (08) :3984-3989
[6]   Experiments in a rotating Hele-Shaw cell [J].
Carrillo, L ;
Magdaleno, FX ;
Casademunt, J ;
Ortin, J .
PHYSICAL REVIEW E, 1996, 54 (06) :6260-6267
[7]   Diffuse-interface approach to rotating Hele-Shaw flows [J].
Chen, Ching-Yao ;
Huang, Yu-Sheng ;
Miranda, Jose A. .
PHYSICAL REVIEW E, 2011, 84 (04)
[8]   Efficient numerical scheme for a new hydrodynamically-coupled conserved Allen-Cahn type Ohta-Kawaski phase-field model for diblock copolymer melt [J].
Chen, Chuanjun ;
Zhang, Jun ;
Yang, Xiaofeng .
COMPUTER PHYSICS COMMUNICATIONS, 2020, 256
[9]   Fast, provably unconditionally energy stable, and second-order accurate algorithms for the anisotropic Cahn-Hilliard Model [J].
Chen, Chuanjun ;
Yang, Xiaofeng .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2019, 351 :35-59
[10]   Efficient numerical scheme for a dendritic solidification phase field model with melt convection [J].
Chen, Chuanjun ;
Yang, Xiaofeng .
JOURNAL OF COMPUTATIONAL PHYSICS, 2019, 388 :41-62