×

含有可溶性表面活性剂的相场移动接触线模型。 (英语) Zbl 1453.76146号

摘要:提出了含有可溶性表面活性剂的两相系统的相场移动接触线模型。通过引入标量辅助变量,将原自由能泛函转化为等价形式,得到了一个新的调节系统。所得模型由两个Cahn-Hilliard型方程和不可压缩的变密度Navier-Stokes方程以及运动接触线的广义Navier边界条件组成。我们证明了该模型满足总能量随时间的耗散。为了数值求解这样一个复杂系统,我们发展了一个具有无条件能量稳定性的非线性耦合格式。采用基于压力稳定的分裂方法求解Navier-Stokes方程。采用一些微妙的隐式-显式处理来离散对流和应力项。人为地添加了一个稳定项,以平衡与流固界面处表面能相关的显式非线性项。我们严格证明了该方案能够保持离散能量耗散。采用高效的交错网格有限差分方法进行空间离散。二维和三维数值结果表明了该方案的准确性和能量稳定性。利用我们的模型和数值格式,我们研究了液滴在固体壁上的润湿行为。数值结果表明,表面活性剂可以通过改变接触角的大小来影响液滴的润湿性能。

MSC公司:

76M20码 有限差分方法在流体力学问题中的应用
76D05型 不可压缩粘性流体的Navier-Stokes方程
76T06型 液-液双组分流动
76D45型 不可压缩粘性流体的毛细管(表面张力)
6500万06 含偏微分方程初值和初边值问题的有限差分方法
65Z05个 科学应用
65个M12 含偏微分方程初值和初边值问题数值方法的稳定性和收敛性

软件:

COMSOL公司
PDF格式BibTeX公司 XML格式引用
全文: 内政部 链接

参考文献:

[1] Engblom,S。;Do-Quang,M。;Amberg,G。;Tornberg,A.-K.,《关于两相流体中表面活性剂的扩散界面建模和模拟》,Commun。计算。物理。,14, 879-915 (2013) ·Zbl 1373.76334号
[2] 刘,H。;Zhang,Y.,用可溶性表面活性剂对液滴动力学进行相场建模,J.Comput。物理。,229, 9166-9187 (2010) ·Zbl 1427.76187号
[3] 赖,M.-C。;Tseng,Y.-H。;Huang,H.,用浸没边界法对含表面活性剂的移动接触线进行数值模拟,Commun。计算。物理。,8, 735 (2010) ·Zbl 1364.76125号
[4] 刘,H。;Ba,Y。;吴,L。;Li,Z。;Xi,G。;Zhang,Y.,含不溶性表面活性剂液滴动力学的混合格子Boltzmann和有限差分方法,J.流体力学。,837, 381-412 (2018) ·Zbl 1419.76521号
[5] 徐,J.-J。;Yang,Y。;Lowengrub,J.,具有不溶性表面活性剂的两相流的水平集连续体方法,J.计算。物理。,231, 5897-5909 (2012) ·Zbl 1522.76105号
[6] 徐,J.-J。;Ren,W.,《带有移动接触线和不溶性表面活性剂的两相流的水平集方法》,J.Compute。物理。,263, 71-90 (2014) ·Zbl 1349.76554号
[7] 徐,J.-J。;Li,Z。;Lowengrub,J。;Zhao,H.,表面活性剂界面流动的水平集方法,J.Compute。物理。,212, 590-616 (2006) ·Zbl 1161.76548号
[8] 阿尔克,A。;Bothe,D.,基于流体体积法的流体界面可溶性表面活性剂三维数值模拟,fluid Dyn。马特。程序。,5, 345-372 (2009)
[9] 詹姆斯·A·J。;Lowengrub,J.,《不溶性表面活性剂界面流动的表面活性剂-流体体积守恒法》,J.Compute。物理。,201, 685-722 (2004) ·兹比尔1061.76062
[10] 穆拉多格鲁,M。;Tryggvason,G.,《用可溶性表面活性剂计算界面流动的前跟踪方法》,J.Compute。物理。,227, 2238-2262 (2008) ·Zbl 1329.76238号
[11] 范德斯曼,R。;Van der Graaf,S.,表面活性剂在平面和液滴界面上吸附的扩散界面模型,Rheol。《学报》,46,3-11(2006)
[12] 泰森,O。;Gompper,G.,《自发乳化的晶格-玻尔兹曼研究》,《欧洲物理学》。J.B,康登斯。物质复杂系统。,11, 91-100 (1999)
[13] 张,L。;康(Kang,Q.)。;姚,J。;高,Y。;孙,Z。;刘,H。;Valocchi,A.J.,基于数字核心技术的液气两相流孔隙尺度模拟,科学。中国,Technol。科学。,58, 1375-1384 (2015)
[14] 丰塞卡,I。;莫里尼,M。;Slastikov,V.,《泡沫稳定性中的表面活性剂:相场模型》,Arch。定额。机械。分析。,183, 411-456 (2007) ·Zbl 1107.76076号
[15] 顾S。;张,H。;Zhang,Z.,二元流体-表面活性剂体系的能量稳定有限差分格式,J.Compute。物理。,270, 416-431 (2014) ·Zbl 1349.76469号
[16] Teigen,K.E。;宋,P。;Lowengrub,J。;Voigt,A.,《含可溶性表面活性剂的两相流扩散界面法》,J.Compute。物理。,230, 375-393 (2011) ·Zbl 1428.76210号
[17] 寇,J。;Sun,S.,用Peng-Robinson状态方程对非等温扩散界面两相流进行热力学一致性模拟,J.Compute。物理。,371, 581-605 (2018) ·Zbl 1415.76644号
[18] Garcke,H。;Lam,K.F。;Stinner,B.,两相流中可溶性表面活性剂的扩散界面模拟,Commun。数学。科学。,12, 1475-1522 (2014) ·Zbl 1319.35309号
[19] Alpak,F。;萨马尔季奇,A。;Frank,F.,使用相场方法的有限差分实现的数字岩石图像上孔隙尺度两相流的分布式并行直接模拟器,J.Pet。科学。工程,166806-824(2018)
[20] Alpak,F.O。;里维埃尔,B。;Frank,F.,使用非平衡润湿边界条件直接模拟多孔介质中两相流的相场方法,计算。地质科学。,20, 881-908 (2016) ·Zbl 1391.76797号
[21] 沈杰。;Yang,X.,两相复杂流体相场模型的解耦能量稳定方案,SIAM J.Sci。计算。,36,B122-B145(2014)·Zbl 1288.76057号
[22] 沈杰。;Yang,X.,两相不可压缩流相场模型的解耦能量稳定格式,SIAM J.Numer。分析。,53, 279-296 (2015) ·Zbl 1327.65178号
[23] Jackmin,D.,使用相场建模计算两相Navier-Stokes流,J.Compute。物理。,155, 96-127 (1999) ·Zbl 0966.76060号
[24] 岳,P。;Feng,J.J。;刘,C。;沈,J.,《模拟复杂流体两相流的扩散界面法》,J.流体力学。,515, 293-317 (2004) ·Zbl 1130.76437号
[25] Sun,S.,《利用重力和毛细作用进行达西尺度相平衡建模》,J.Compute。物理。,399,第108908条pp.(2019)·Zbl 1453.76207号
[26] 丁·H。;拼写,P.D。;Shu,C.,具有大密度比的不可压缩两相流的扩散界面模型,J.Comput。物理。,226, 2078-2095 (2007) ·Zbl 1388.76403号
[27] 沈杰。;Yang,X.,不同密度和粘度的两相不可压缩流动的相场模型及其数值近似,SIAM J.Sci。计算。,32, 1159-1179 (2010) ·Zbl 1410.76464号
[28] 刘,C。;Riviere,B.,Cahn-Hilliard-Navier-Stokes方程的间断Galerkin方法的数值分析(2018),arXiv预印本
[29] 朱,G。;陈,H。;姚,J。;Sun,S.,流体动力学耦合相场模型的高效能量稳定方案,应用。数学。型号。,70, 82-108 (2019) ·Zbl 1462.76133号
[30] 弗兰克·F。;刘,C。;Alpak,F.O。;Berg,S。;Riviere,B.,《使用相场法在孔隙尺度图像上对流动进行直接数值模拟》,SPE J.,231833-831850(2018)
[31] 拉拉吉,M。;郭,H。;格兰特,M。;Zuckermann,M.J.,《表面活性剂对相分离动力学的影响》,J.Phys。康登斯。Matter,46715(1992)
[32] Komura,S。;Kodama,H.,油水表面活性剂系统的二阶参数模型,Phys。修订版E,551722(1997)
[33] Pätzold,G。;Dawson,K.,《表面活性剂和流体动力学存在下相分离的数值模拟》,Phys。E版,526908(1995)
[34] 托特,G.I。;Kvamme,B.,《表面活性剂辅助液相分离的Ginzburg-Landau型模型分析》,Phys。E版,91,第032404条pp.(2015)
[35] 李毅。;Kim,J.,与表面活性剂不混溶的二元混合物相场模型的比较研究,Eur.Phys。J.B,85,340(2012)
[36] 杨,X。;Ju,L.,相场弹性弯曲能量模型的无条件能量稳定性高效线性格式,计算。方法应用。数学。,315, 691-712 (2017) ·Zbl 1439.74165号
[37] 朱,G。;寇,J。;Sun,S。;姚,J。;Li,A.,相场表面活性剂模型的解耦能量稳定方案,计算。物理学。社区。,233, 67-77 (2018) ·Zbl 07694813号
[38] Yun,A。;李毅。;Kim,J.,水-油表面活性剂系统的新相场模型,应用。数学。计算。,229, 422-432 (2014) ·Zbl 1364.76038号
[39] Aland,S。;哈恩,A。;卡勒,C。;Nürnberg,R.,基于尖锐和扩散界面方法的Taylor流动与表面活性剂的比较模拟,(流体界面的传输过程(2017),Springer),639-661·Zbl 1391.76769号
[40] 朱,G。;寇,J。;Sun,S。;姚,J。;Li,A.,带流体流动的相场表面活性剂模型的数值近似,J.Sci。计算。,80, 223-247 (2019) ·Zbl 1448.76171号
[41] 高,M。;Wang,X.-P.,变密度和粘度运动接触线问题相场模型的有效方案,J.Compute。物理。,272, 704-718 (2014) ·Zbl 1349.76464号
[42] 王晓平。;钱,T。;Sheng,P.,化学图形表面上的移动接触线,J.流体力学。,605, 59-78 (2008) ·Zbl 1145.76045号
[43] 高,M。;Wang,X.-P.,移动接触线问题相场模型的梯度稳定格式,J.Compute。物理。,231, 1372-1386 (2012) ·Zbl 1408.76138号
[44] 朱,G。;姚,J。;李,A。;Sun,H。;Zhang,L.,二氧化碳强化采油的孔隙尺度调查,能源燃料,315324-5332(2017)
[45] 张,Z。;徐,S。;Ren,W.,《含不溶性表面活性剂的连续介质模型和运动接触线能量定律的推导》,Phys。流体,26,第062103条,第(2014)页
[46] 风扇,W。;Sun,H。;姚,J。;风扇,D。;Yang,Y.,考虑滑移效应的页岩系统内液体流动均匀化方法,J.Clean。生产,235146-157(2019)
[47] 钱,T。;王晓平。;Sheng,P.,不混溶流动的分子尺度接触线流体动力学,物理学。E版,68,第016306条,pp.(2003)
[48] 钱,T。;王晓平。;Sheng,P.,两相非混相流中移动接触线的幂律滑移剖面,Phys。修订稿。,93,第094501条pp.(2004)
[49] 钱,T。;Wang,X.-P。;Sheng,P.,《动接触线流体力学的变分方法》,J.流体力学。,564, 333-360 (2006) ·Zbl 1178.76296号
[50] Bonart,H。;卡勒,C。;Repke,J.-U.,使用热力学一致的Cahn-Hilliard Navier-Stokes模型对移动接触线问题进行能量稳定模拟的比较(2018),arXiv预印本
[51] 陈,J。;Sun,S。;Wang,X.-P.,自由流和多孔介质耦合系统中两相流模型的数值方法,J.Compute。物理。,268, 1-16 (2014) ·Zbl 1349.76187号
[52] 朱,G。;寇,J。;姚明,B。;Wu,Y.-s。;姚,J。;Sun,S.,带移动接触线和可溶性表面活性剂的两相流热力学一致性建模,J.流体力学。,879, 327-359 (2019) ·Zbl 1430.76475号
[53] Yu,H。;Yang,X.,具有可变密度和粘度的相场移动接触线模型的数值近似,J.Compute。物理。,334, 665-686 (2017) ·Zbl 1375.76201号
[54] 沈杰。;杨,X。;Yu,H.,相场移动接触线模型的高效能量稳定数值格式,J.Compute。物理。,284, 617-630 (2015) ·Zbl 1351.76184号
[55] 杨,X。;Ju,L.,二元流体-表面活性剂相场模型的线性和无条件能量稳定格式,计算。方法应用。数学。,318, 1005-1029 (2017) ·兹比尔1439.76029
[56] 沈杰。;徐,J。;Yang,J.,梯度流的标量辅助变量(SAV)方法,J.Compute。物理。,353, 407-416 (2018) ·Zbl 1380.65181号
[57] 沈杰。;徐,J。;Yang,J.,梯度流的一类新型高效稳健能量稳定方案,SIAM Rev.,61,474-506(2019)·Zbl 1422.65080号
[58] 美国Fjordholm。;米什拉,S。;Tadmor,E.,《不连续地形浅水方程的井平衡和能量稳定格式》,J.Compute。物理。,230, 5587-5609 (2011) ·Zbl 1452.35149号
[59] 杨,X。;Yu,H.,使用不变能量求积方法的相场移动接触线模型的有效二阶无条件稳定格式,SIAM J.Sci。计算。,40,B889-B914(2018)·Zbl 1395.65095号
[60] 科佩蒂,M。;Elliott,C.M.,具有对数自由能的Cahn-Hilliard方程的数值分析,Numer。数学。,63, 39-65 (1992) ·兹比尔0762.65074
[61] 吉蒙德,J.-L。;Salgado,A.,基于压力泊松方程的变密度不可压缩流的分裂方法,J.Compute。物理。,228, 2834-2846 (2009) ·Zbl 1159.76028号
[62] 寇,J。;Sun,S.,具有部分混溶性的多组分两相流热力学一致性建模与模拟,计算。方法应用。数学。,331, 623-649 (2018) ·Zbl 1439.76079号
[63] 寇,J。;Sun,S。;Wang,X.,多组分两相可压缩流动的线性解耦能量稳定数值方法,SIAM J.Numer。分析。,56, 3219-3248 (2018) ·Zbl 1404.65226号
[64] 陈,H。;Sun,S。;Zhang,T.,交错网格上不可压Navier-Stokes方程的一些全离散数值格式的能量稳定性分析,J.Sci。计算。,1-30 (2017)
[65] 穆卡勒德,F。;Mangani,L。;Darwish,M.,计算流体动力学中的有限体积法,(OpenFOAM和Matlab的高级介绍(2016)),3-8·兹比尔1329.76001
[66] 刘,P。;严,X。;姚,J。;Sun,S.,使用两相热-水文-化学耦合模型对碳酸盐岩酸化过程进行建模和分析,Chem。工程科学。(2019)
[67] 李,J。;Yu,B。;王,Y。;Tang,Y。;Wang,H.,使用单网格和多重网格方法研究对流扩散方程的复合格式的计算效率,J.Therm。科学。技术。,10,条款JTST0009 pp.(2015)
[68] 严,X。;Huang,Z。;姚,J。;张,Z。;刘,P。;李毅。;Fan,D.,使用等效连续介质模型和嵌入离散裂缝模型对裂隙-溶洞多孔介质中的水力-机械耦合进行数值模拟,Adv.Water Resour。,126, 137-154 (2019)
[69] 曾庆东。;姚,J。;Shao,J.,通过EDFM-XFEM混合方法研究各向异性多孔弹性介质中的水力压裂,计算。岩土工程。,105, 51-68 (2019)
[70] 王,X。;姚,J。;龚,L。;Sun,H。;Yang,Y。;张,L。;李毅。;Liu,W.,水力裂缝和裂缝网络中支撑剂沉积和运移特性的数值模拟,J.Pet。科学。工程,183,第106401条pp.(2019)
[71] 张,L。;Jing,W。;Yang,Y。;Yang,H。;郭毅。;Sun,H。;赵,J。;姚,J.,《利用数字岩心模拟技术进行渗透率计算的研究》,能源,123273(2019)
[72] 蔡,X。;马沙尔,H。;沃纳,M。;Deutschmann,O.,用于OpenFOAM®三维润湿过程数值模拟的自适应网格细化相场方法,(第二届多尺度多相过程工程国际研讨会(MMPE))。第二届多尺度多相过程工程国际研讨会,德国汉堡(2014)
[73] A.COMSOL,COMSOL Multiphysics User's Guide,2012年第4.3版。
[74] 风扇,W。;Sun,H。;姚,J。;风扇,D。;Yang,Y.,基于均匀化理论的考虑多种机制和非均质性的页岩气高标度输运模型,J.Pet。科学。工程,183,第106392条pp.(2019)
[75] 科斯拉,P。;Rubin,S.,一个对角占优的二阶精确隐式格式,计算。流体,2207-209(1974)·Zbl 0335.76009号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。