×

Peng-Robinson状态方程扩散界面模型的无条件能量稳定线性格式。 (英语) Zbl 1398.65305号

小结:本文研究了用Peng-Robinson状态方程描述石油工业中烃类流体真实状态的扩散界面模型的数值解。由于该模型中源项的强非线性,如何设计适当的时间离散化以在离散水平上保持系统的能量耗散规律是一个重大挑战。基于“不变能量求积”方法和惩罚公式,我们开发了有效的一阶和二阶时间步长格式来求解单组分两相流体问题。在这两种方案中,所得到的时间半离散化导致在每个时间步长具有对称正定空间算子的线性系统。我们严格证明了它们在时间离散意义下的无条件能量稳定性。为了验证所提线性方案的准确性和稳定性,并通过与实验室数据的比较,研究目标模型的物理可靠性,还进行了二维和三维空间的各种数值模拟。

MSC公司:

65N30型 含偏微分方程边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65牛顿50 涉及偏微分方程的边值问题的网格生成、细化和自适应方法
49S05号 物理学变分原理
PDF格式BibTeX公司 XML格式引用
全文: DOI程序

参考文献:

[1] Chen,R.,Ji,G.,Yang,X.,Zhang,H.:相场囊泡膜模型的解耦能量稳定方案。J.计算。物理学。302, 509-523 (2015) ·Zbl 1349.76842号 ·doi:10.1016/j.jcp.2015.09.025
[2] Copetti,M.,Elliott,C.:具有对数自由能的Cahn-Hilliard方程的数值分析。数值。数学。63, 39-65 (1992) ·Zbl 0762.65074号 ·doi:10.1007/BF01385847
[3] Davis,H.T.:相、界面和薄膜的统计力学。VCH,纽约(1996)·Zbl 0884.0006
[4] Du,Q.,Li,M.,Liu,C.:相场Navier-Stokes囊泡-流体相互作用模型的分析。数字化信息系统。Contin公司。动态。系统。B 8,539-556(2007)·Zbl 1172.74016号 ·doi:10.3934/dcdsb.2007.8.539
[5] Du,Q.,Liu,C.,Wang,X.:囊泡膜弹性弯曲能量数值研究中的相场方法。J.计算。物理学。198, 450-468 (2004) ·Zbl 1116.74384号 ·doi:10.1016/j.jcp.2004.01.029
[6] Du,Q.,Liu,C.,Wang,X.:模拟囊泡膜在三维弹性弯曲能量下的变形。J.计算。物理学。212, 757-777 (2005) ·Zbl 1086.74024号 ·doi:10.1016/j.jcp.2005.07.020
[7] Elliott,C.M.,Garcke,H.:关于简并迁移率的Cahn-Hilliard方程。SIAM J.数学。分析。27, 404-423 (1996) ·Zbl 0856.35071号 ·doi:10.1137/S0036141094267662
[8] Eyre,D.J.:无条件梯度稳定时间推进Cahn-Hilliard方程。在:微观结构演化的计算和数学模型(加利福尼亚州旧金山,1998年),材料研究学会专题讨论会论文集,第529卷,第39-46页。宾夕法尼亚州沃伦代尔(1998)·兹比尔1397.76070
[9] Firoozabadi,A.:油气藏热力学。McGraw-Hill,纽约(1999)
[10] Frenkel,D.,Smit,B.:理解分子模拟:从算法到应用。学术出版社,加利福尼亚州圣地亚哥(2001)·Zbl 0889.65132号
[11] Guillén-González,F.,Tierra,G.:关于Cahn-Hilliard扩散界面模型的线性方案。J.计算。物理学。234, 140-171 (2013) ·Zbl 1284.35025号 ·doi:10.1016/j.jp.2012.09.020
[12] Han,D.,Brylev,A.,Yang,X.,Tan,Z.:两相不可压缩流相场模型的二阶全离散能量稳定格式的数值分析。科学杂志。计算。70, 965-989 (2016) ·兹比尔1397.76070 ·doi:10.1007/s10915-016-0279-5
[13] Haugen,K.B.,Firoozabadi,A.:多组分非平衡流体相界面的成分。化学杂志。物理学。130, 064707 (2009) ·doi:10.1063/1.3072793
[14] Kim,J.:多组分流体流动的相场模型。Commun公司。计算。物理学。12, 613-661 (2012) ·Zbl 1373.76030号 ·doi:10.4208/cicp.301110.040811a
[15] Kou,J.,Sun,S.:利用流体界面梯度理论模拟表面张力的自适应有限元方法。J.计算。申请。数学。255, 593-604 (2014) ·Zbl 1291.76205号 ·doi:10.1016/j.cam.2013.06.017
[16] Kou,J.,Sun,S.:具有几何平均影响参数的多组分两相界面模型的数值方法。SIAM J.科学。计算。37,B543-B569(2015)·Zbl 1326.76062号 ·数字对象标识代码:10.1137/140969579
[17] Kou,J.,Sun,S.:用Peng-Robinson状态方程和各种边界条件模拟多组分两相界面模型的无条件稳定方法。J.计算。申请。数学。291, 158-182 (2016) ·Zbl 1329.76178号 ·doi:10.1016/j.cam.2015.02.037
[18] Kou,J.,Sun,S.,Wang,X.:利用流体界面梯度理论模拟多组分混合物表面张力的有效数值方法。计算。方法应用。机械。工程292、92-106(2015)·Zbl 1423.76250号 ·doi:10.1016/j.cma.2014.10.023
[19] Lin,H.,Duan,Y.Y.:丙烷(r-290)和异丁烷(r-600a)的表面张力测量(253至333)K.J.Chem。工程数据481360-1363(2003)·doi:10.1021/je034093m
[20] Miqueu,C.,Mendiboure,B.,Graciaa,A.,Lachaise,J.:用流体界面梯度理论模拟多组分混合物的表面张力。工业工程化学。第44号决议,3321-3329(2005年)·doi:10.1021/ie049086l
[21] Peng,D.,Robinson,D.:一个新的双常数状态方程。工业工程化学。芬丹。15, 59-64 (1976) ·doi:10.1021/i160057a011
[22] 乔,Z.,孙,S.:使用带有彭-罗宾逊状态方程的扩散界面模型进行双相流体模拟。SIAM J.科学。计算。36,B708-B728(2014)·Zbl 1305.76117号 ·数字对象标识代码:10.1137/130933745
[23] Rongy,L.,Haugen,K.B.,Firoozabadi,A.:二元非平衡流体相中费基扩散和自然对流的混合。AIChE J.58,1336-1345(2012)·doi:10.1002/aic.12685
[24] Shen,J.,Yang,X.:Allen-Cahn和Cahn-Hilliard方程的数值近似。数字化信息系统。Conti公司。动态。系统。A 28,1669-1691(2010)·Zbl 1201.65184号 ·doi:10.3934/dcds.2010.28.1669
[25] Sun,S.,Wheeler,M.F.:多孔介质中反应输运的对称和非对称间断Galerkin方法。SIAM J.数字。分析。43, 195-219 (2005) ·Zbl 1086.76043号 ·doi:10.1137/S003614290241708X
[26] Sun,S.,Wheeler,M.F.:基于局部问题的后验误差估计量,用于反应迁移的间断Galerkin近似。计算。地质科学。11, 87-101 (2007) ·Zbl 1117.76061号 ·doi:10.1007/s10596-007-9041-2
[27] van der Waals,J.:连续密度变化假设下的毛细现象热力学理论。《统计物理学杂志》。20, 197-244 (1893) ·Zbl 1245.82006年
[28] Wang,X.,Ju,L.,Du,Q.:相场弹性弯曲能量模型的高效稳定指数时间差分Runge-Kutta方法。J.计算。物理学。316, 21-38 (2016) ·Zbl 1349.92015年9月 ·doi:10.1016/j.jcp.2016.04.004
[29] Wang,C.,Wise,S.M.:修正相场晶体方程的能量稳定且收敛的有限差分格式。SIAM J.数字。分析。49, 945-969 (2011) ·Zbl 1230.82005年 ·doi:10.1137/090752675
[30] Wheeler,M.F.、Wick,T.、Wollner,W.:加压裂缝相场方法的增强图解法。计算。方法应用。机械。工程271、69-85(2014)·Zbl 1296.65170号 ·doi:10.1016/j.cma.2013.12.005
[31] Yang,X.:均聚物共混物相场模型的线性、一阶和二阶无条件能量稳定数值格式。J.计算。物理学。327, 294-316 (2016) ·Zbl 1373.82106号 ·doi:10.1016/j.jcp.2016.09.029
[32] Yang,X.:二元流体-表面活性剂系统Cahn-Hilliard相场模型的数值近似。科学杂志。计算。(2017). doi:10.1007/s10915-017-0508-6·Zbl 1456.65080号 ·doi:10.1007/s10915-017-0508-6
[33] Yang,X.,Han,D.:相场晶体方程的线性一阶和二阶无条件能量稳定格式。J.计算。物理学。330, 1116-1134 (2017) ·Zbl 1380.65209号 ·文件编号:10.1016/j.jcp.2016.10.020
[34] Yang,X.,Ju,L.:相场弹性弯曲能量模型的具有无条件能量稳定性的有效线性格式。计算。方法应用。机械。工程315691-712(2017)·Zbl 1439.74165号 ·doi:10.1016/j.cma.2016.10.041
[35] Yang,X.,Ju,L.:二元流体-表面活性剂相场模型的线性和无条件能量稳定格式。计算。方法应用。机械。工程3181005-1029(2017)·兹比尔1439.76029 ·doi:10.1016/j.cma.2017.02.011
[36] Yang,X.,Zhao,J.,Wang,Q.:基于不变能量求积方法的分子束外延生长模型的数值近似。J.计算。物理学。333, 104-127 (2017) ·Zbl 1375.82121号 ·doi:10.1016/j.jcp.2016.12.025
[37] Zhao,J.,Wang,Q.,Yang,X.:基于不变能量正交化方法的相场枝晶生长模型的数值近似。埋。J.数字。方法工程110、279-300(2017)·Zbl 1365.74138号 ·doi:10.1002/nme.5372
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。