Skip to main content

CMG热采STARS模拟器粘度非线性混合的使用

关键字用法:

  *VSMIXCOMP comp_des
  *VSMIXENDP xlow xhigh
  *VSMIXFUNC f(1)...f(11)

  comp_des  在粘度非线性混合中关键组分的名称或序号。
  xlow     对应于第一个表格元素的横坐标。允许的范围是0~xhigh。
  xhigh    对应于最后一个表格元素的横坐标。允许的范围是xlow~1。
  f(1)...f(11)   表格元素定义的非线性混合规则函数,这个函数应该是单调递增的,并且应该是平滑的。

  缺省:
    如果*VSMIXCOMP不存在,则假设对所有组分进行线性混合。
    如果*VSMIXENDP不存在,则假设xlow=0,xhigh=1。
    如果*VSMIXFUNC不存在,则表格输入f(i)=(i-1)/10,i=1到11,对应于0到1之间的线性间隔。

算法:

  首先线性混合的规则为

              
  使用表征非线性的代替,其中a为关键组分。但是仍然需要满足和为1的条件,即
              
  由上式可以得出
              
  因此对于粘度的非线性混合规则可以表示为
              
  对于任意的就有三个范围去得到相应的的值
              
              的值从输入的表中得到
              
 可以用下图表示其含义:
 通过上述非线性混合的式子可以推导出的计算方法
              
              
  将和中间等距间隔的9个值分别代入上述式子中,得到*VSMIXFUNC关键字中需要输入的11个值。
  如果,则非线性混合规则退化为线性混合规则。

查看输出:

  通过*OUTPRN和*GRID的*VISCCMP子关键字,以及*OUTSRF *GRID的*VISWCMP和*VISOCMP子关键字和*SPECIAL检查粘度非线性混合关键组分的组成。

例子:

  以下数据和计算都是在同一温度下进行。一种沥青粘度为,一种纯溶解剂的浓度为,主要组分为溶解剂。其两者混合后的粘度表示为下图所示:

  相应的数据如下表所示:
0.000300000
0.015268143
0.030191506
0.045109375
0.06050090
0.07518580
0.0905816
0.1051808
0.120824
0.135633
0.150604
1.000600
  对于上述关系,溶解剂浓度从0.0到0.15之间可认为是非线性的,根据算法中计算的公式,就可以得到相应的值:
0.00012.610
0.01512.500.018064
0.03012.160.072227
0.04511.600.162360
0.06010.820.288023
0.0759.830.447603
0.0908.670.634493
0.1057.500.822474
0.1206.700.949045
0.1356.450.991380
0.1506.400.998959
  关键字设置为:
    *COMPNAME 'WATER' 'BITUMEN' 'SOLVENT'
    *VISCTABLE
    **TEMP
     75     0     300000   600
    *VSMIXCOMP 'SOLVENT'
    *VSMIXENDP 0.0 0.15
    *VSMIXFUNC
     **0.0  0.015    0.03   0.045   0.06    0.075
       0   0.0018064 0.072227 0.162360 0.288023  0.447603
     **0.09   0.105   0.12    0.135  0.15
      0.634493 0.822474 0.949045 0.991380 0.998959
 注:此例子为CMG帮助手册中的例子。

 一个模拟程序中只有一相可以使用非线性混合规则,即水相或者油相。如果是水相的非线性混合,需要注意的是在*MODEL关键字下的组分数和*COMPNAME关键字下的组分顺序,以及相应的各组分性质关键字的设置顺序。

说明:

  粘度非线性混合选项的计算方式与*DNMIXCOMP关键字说明的密度选项相同,唯一不同的是进行加权的量是粘度的对数。

  以上内容由本人参考帮助文档及模拟经验所得,如有错误还请您不吝赐教。

Comments

Popular posts from this blog

使用PHP Webhook方式打造Telegram Bot

一、找BotFather拿到bot token     在telegram中私聊BotFather建立自己的bot,给bot取名,名字必须要以bot结尾。建好后自己的bot就有一个唯一的token,类似下面的一串字符 164354723:AAEjT6-IyNoXjt7miD0dwa-P5VmDTtHQC8 二、确认bot响应文件的位置     在写好bot响应文件后,要把bot放在网络上的一个位置,并且这个位置必须要加密的,即以https开头的一串网址。比如响应文件的名称为telbot.php,把它放在下面这个网址的位置: https://my.webhost.com/ 164354723:AAEjT6-IyNoXjt7miD0dwa-P5VmDTtHQC8 /telbot.php 上面网址中的红色设置和bot的token一样是为了确定这个唯一的位置,当然也可以任意设置。 三、告诉Telegram响应文件的位置 Telegram用下面网址的形式来设定webhook响应方式 https://api.telegram.org/bot [myauthorization-token] /setwebhook?url= [myboturl] 按照上面的网址形式,把自己创建的bot的token以及响应文件的位置填入,然后在浏览器中运行一下即可设置成功。比如: https://api.telegram.org/bot164354723:AAEjT6-IyNoXjt7miD0dwa-P5VmDTtHQC8/setwebhook?url=https://my.webhost.com/164354723:AAEjT6-IyNoXjt7miD0dwa-P5VmDTtHQC8/telbot.php 设置成功后,页面会显示下面的内容: {"ok":true,"result":true,"description":"Webhook is already set"} 四、在Telegram中给自己的bot发消息进行验证 php响应文件例子 <?php  define('BOT_TOKEN', 'YOURBOT:TOK

MatLab中patch函数的基本用法

patch是用来构建多边形的一个基本函数。 用法一 patch(X,Y,C) patch(X,Y,Z,C) patch( 'XData' ,X, 'YData' ,Y) patch( 'XData' ,X, 'YData' ,Y, 'ZData' ,Z) 1.1 说明 patch(X,Y,C)用来构建一个或者多个可填充的多边形,其使用X和Y作为每个点的坐标值,patch将会按顺序连接每个点。如果要得到一个多边形,将X和Y设置为向量;如果要得到多个多边形,将X和Y设置为矩阵,没一列对应一个多边形。C决定多边形的颜色,可以是系统认定的字符,也可以是一个数值,也可以是RGB向量。 patch(X,Y,Z,C)用来构建三维坐标下的多边形。 patch(‘XData’,X,’YData’,Y)和patch(‘XData’,X,’YData’,Y,’ZData’,Z)的用法与patch(X,Y,C)和patch(X,Y,Z,C)的用法类似,只是不设定颜色。 1.2 例子 1.2.1 x = [ 0 1 1 0 ] ; y = [ 0 0 1 1 ] ; patch(x,y, 'red' ) x和y都是1*4的向量,表示将四个点(0,0)、(1,0)、(1,1)和(0,1)依次连接,最后闭合形成一个四边形,设定颜色为红色。 1.2.2 x2 = [ 2 5 ; 2 5 ; 8 8 ] ; y2 = [ 4 0 ; 8 2 ; 4 0 ] ; patch(x2,y2, 'green' ) x2和y2都是3*2的向量,两列表示画两个多边形。第一个多边形连接的点依次是(2,4)、(2,8)和(8,4),第二个多边形连接的点依次是(5,0)、(5,2)和(8,0),颜色设定为绿色。 1.2.3 如果上例的三角形第一个是红色,第二个是绿色,那么patch代码修改为 x2 = [ 2 5 ; 2 5 ; 8 8 ] ; y2 = [ 4 0 ; 8 2 ; 4 0 ] ; patch(x2(:, 1 ),y2(:, 1 ), 'red' ) pat

telegram中的Sci-Hub机器人,又一文献下载利器

或许你看到标题会问什么是telegram,什么是Sci-Hub?请听我一一道来。 什么是Sci-Hub Sci-Hub是一个线上 数据库 ,其上提供48,000,000篇科学学术论文和文章。网站透过“.edu”代理服务器访问相关页面,每天会上传新的论文文章。2011年,哈萨克研究生亚历珊卓·艾尔巴金(Alexandra  Elbakyan)因为研究论文成本过高,每篇论文在付费墙机制下通常需要花费30美元,而决定成立Sci-Hub。2014年,学术界开始预测网站将会发展为类似Napster的服务。不过到了2015年,学术出版社爱思唯尔向纽约地方法院提交诉讼,指控Sci-Hub已经侵犯版权。纽约地方法院在2015年10月28日仍下令Sci-Hub原本使用的网域名称“Sci-Hub.org”必须终止。爱思唯尔在法院上获得胜诉后,一群研究人员、作家和艺术家则连署一封表态支持Sci-Hub和创世纪图书馆的公开信,声称这次诉讼对于世界各地的研究人员是“重大打击”,并指出:“它同样贬低我们、作者、编辑和读者。它寄生于我们的劳动,它阻挠我们为大众服务,它阻拦我们进入。”而该计划于11月因法院命令中止后,在同一个月内便改用网域名称“.io”重新上线,并开放使用Tor浏览。2016年1月时,Sci-Hub平均每天约有200,000人访问,Sci-Hub则声称网站服务每天平均有数十万次档案请求。  Sci-Hub是目前已知第一个提供大量自动且免费的付费学术论文的网站,使用者不需要事前订阅或付款,就能够使用原本存放在付费数据库的论文文章,并提供搜寻原先出版社网站内的文件档案服务。 以上介绍来源于维基百科词条 Sci-Hub Sci-Hub网站被屡次下线,但是又通过更换域名重新上线。以下三个网址经测试可以使用:  http://www.sci-hub.bz/   http://www.sci-hub.ac/   http://www.sci-hub.cc/   广大学者将自己的文章发表至学术期刊(免费或者支付版面费),然而当需要查看其他学者的文章时还需要向出版商付费,你是不是也觉得这完全阻碍了科学文化的传播。艾尔巴金在为自己辩护时援引联合国《世界人权宣言》第二十七条所提的:“人人有权自由参加社会之文化生活,欣赏艺术,并共同襄享科学进步及其利益。”关于这个问题,另一位被称