求助：Sorption 计算吸附等温线

霍亚金

Sorption 计算吸附等温线时，高压下，逸度和压力怎么转换，有没有相关的脚本或者软件

ABetaCarw

1. ################################################################################
   
2. #!perl                                                                         #
   
3. #                                                                              #
   
4. # PENG-ROBINSON EQUATION OF STATE                                              #
   
5. #                                                                              #
   
6. # Author: Reinier Akkermans (Biovia)                                           #
   
7. # Version: 1.0                                                                 #
   
8. # Tested on: Materials Studio 7.0                                              #
   
9. # Required modules: Materials Visualizer                                       #
   
10. #                                                                              #
    
11. # Description: This script calculates the fugacity and pressure for a range of #
    
12. # volume using the Peng-Robinson Equation of state. The data is output in a    #
    
13. # study table                                                                  #
    
14. #                                                                              #
    
15. ################################################################################
    
16. 
    
17. 
    
18. use strict;
    
19. use Getopt::Long;
    
20. use MaterialsScript qw(:all);
    
21. 
    
22. 
    
23. use constant
    
24. {
    
25.     GAS_CONSTANT => 8.314, # J/mol/K
    
26.     # constants in the Peng-Robinson EoS
    
27.     PR_a => 0.457235,
    
28.     PR_b => 0.077796,
    
29.     PR_kappa0 => 0.37464,
    
30.     PR_kappa1 => 1.54226,
    
31.     PR_kappa2 => -0.26992,
    
32. };
    
33. 
    
34. 
    
35. ################################################################################
    
36. #  BEGIN USER INPUT                                                            #
    
37. ################################################################################
    
38. 
    
39. 
    
40. # example: CO2
    
41. my $Tc = 304.25; # temperature at critical point in K
    
42. my $Pc = 7380000; # pressure at critical point in Pa
    
43. my $omega = 0.228; # acentric factor
    
44. 
    
45. 
    
46. my $temperature = 298; # K   
    
47. my $volume_start = 0.0004; # m^3/mol
    
48. my $volume_end = 0.1; # m^3/mol
    
49. my $volume_steps = 100;
    
50. 
    
51. 
    
52. ################################################################################
    
53. #  END USER  INPUT                                                             #
    
54. ################################################################################
    
55. 
    
56. 
    
57. my $Tr = $temperature/$Tc;
    
58. my $RTc = GAS_CONSTANT*$Tc; # J/mol
    
59. my $b = PR_b*$RTc/$Pc; # m^3/mol
    
60. my $kappa = PR_kappa0+PR_kappa1*$omega+PR_kappa2*$omega**2;
    
61. my $alpha = (1+$kappa*(1-sqrt($Tr)))**2;
    
62. my $aTc = PR_a*$RTc**2/$Pc; # Pa*(m^3/mol)^2
    
63. my $aT = $aTc*$alpha; # Pa*(m^3/mol)^2
    
64. my $RT = GAS_CONSTANT*$temperature;
    
65. 
    
66. 
    
67. my $std = Documents->New("PengRobinson.std");
    
68. $std->ColumnHeading(0) = "Volume (m^3/mol)";
    
69. $std->ColumnHeading(1) = "Pressure (kPa)";
    
70. $std->ColumnHeading(2) = "Fugacity (kPa)";
    
71. $std->ColumnHeading(3) = "A";
    
72. $std->ColumnHeading(4) = "B";
    
73. $std->ColumnHeading(5) = "Z";
    
74. $std->ColumnHeading(6) = "check";
    
75. 
    
76. 
    
77. my $scaleFactor = exp(log($volume_end/$volume_start)/($volume_steps-1));
    
78. my $volume = $volume_start;
    
79. for(my $i = 0; $i < $volume_steps; $i++)
    
80. {
    
81.     $volume *= $scaleFactor if $i > 0;
    
82.     my $pressure = $RT/($volume-$b)
    
83.         -$aT/($volume*($b+$volume)+$b*($volume-$b));
    
84.    
    
85.     my $A = $aT*$pressure/$RT**2;
    
86.     my $B = $b*$pressure/$RT;
    
87.     my $Z = $pressure*$volume/$RT;
    
88.    
    
89.     # sanity check; should evaluate to zero
    
90.     my $check = $Z**3+($B-1)*$Z**2+($A-2*$B-3*$B**2)*$Z+$B**3+$B**2-$A*$B;
    
91. 
    
92. 
    
93.     my $lnfP = $Z-1
    
94.         - log($Z-$B)
    
95.         - $A/(sqrt(8)*$B)*log(($Z+(1+sqrt(2))*$B)/($Z-(1-sqrt(2))*$B));
    
96.         
    
97.     my $fugacity = exp($lnfP)*$pressure;
    
98.    
    
99.     $std->Cell($i,0) = $volume;
    
100.     $std->Cell($i,1) = $pressure/1000;
     
101.     $std->Cell($i,2) = $fugacity/1000;
     
102.     $std->Cell($i,3) = $A;
     
103.     $std->Cell($i,4) = $B;
     
104.     $std->Cell($i,5) = $Z;
     
105.     $std->Cell($i,6) = $check;
     
106. }
     

复制代码

脚本来自Reinier Akkermans (Biovia)

霍亚金

ABetaCarw 发表于 2021-6-8 11:33
脚本来自Reinier Akkermans (Biovia)

这个脚本怎么用,是在ms里吗？

ABetaCarw

霍亚金 发表于 2021-6-8 14:50
这个脚本怎么用，是用matlab吗

哎……本想多说两句的，还是不说了。是MS的perl脚本。至于怎么用，你百度或者看手册吧。

霍亚金

ABetaCarw 发表于 2021-6-8 14:53
哎……本想多说两句的，还是不说了。是MS的perl脚本。至于怎么用，你百度或者看手册吧。

好的，感谢感谢

霍亚金

ABetaCarw 发表于 2021-6-8 14:53
哎……本想多说两句的，还是不说了。是MS的perl脚本。至于怎么用，你百度或者看手册吧。

老师，您有计算不同温度压力下的气体密度的脚本吗
可以分享一下吗

ABetaCarw

霍亚金 发表于 2021-6-8 15:30
老师，您有计算不同温度压力下的气体密度的脚本吗
可以分享一下吗

要不还是看看物化书吧……

霍亚金

ABetaCarw 发表于 2021-6-8 15:43
要不还是看看物化书吧……

老师，我用这个脚本算氢气，怎么逸度系数大于1，这正常吗