 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1815 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[求助]
關(guān)于二元混合溶劑溶解度數(shù)據(jù)利用wilson、NRTL模擬的1stOPT代碼 已有2人參與
|
|||
|
各位大神，我測得了不同比例的某物質(zhì)在二元混合溶劑中的溶解度隨溫度的變化，欲選用wilson、NRTL方程對其回歸求得參數(shù)，下面是我的數(shù)據(jù)，但是1stOp軟件t程序不大會寫，總是提示function有錯，請求高手指教 x2=0 283.3 3.02194E-05 288.31 0.000157704 293.3 0.0004041 298.25 0.000686379 303.15 0.001026437 308.18 0.001405125 313.05 0.001810204 318.1 0.002726428 323.05 0.003769643 x2=0.05 283.13 0.000349201 288.15 0.00046566 293.05 0.000981802 298.1 0.001500779 303.23 0.002474932 308.3 0.003373367 313.35 0.004337507 318.35 0.005850682 323.3 0.007716984 x2=0.1 283.2 0.000828269 288.23 0.001597956 293.16 0.002383646 298.24 0.003548518 303.2 0.004900833 308.15 0.006296129 313.15 0.008064994 318.2 0.010216362 323.2 0.013075324 x2=0.15 283.15 0.002970424 288.15 0.003905376 293.06 0.005280804 298.15 0.006616207 303.15 0.008243561 308.15 0.010303864 313.2 0.013119118 318.15 0.016480102 323.07 0.020616259 x2=0.2 283.1 0.005204799 288.15 0.006279544 293.05 0.007921776 298.05 0.00986987 303.15 0.012566941 308.15 0.015421326 313.18 0.019498661 318.25 0.024651302 323.15 0.030872327 x2=0.25 283.22 0.008048279 288.27 0.009556376 293.3 0.011981522 298.2 0.014520227 303.2 0.01789219 308.1 0.022389299 313.05 0.027783872 318.07 0.035092411 323.1 0.044652622 x2=0.3 283.05 0.011149853 288.05 0.013651819 293.06 0.016540268 298.07 0.020332495 303.1 0.024782096 308.23 0.030544405 313.25 0.037978064 318.35 0.048406476 323.2 0.061434365 x2=0.35 283.15 0.015106036 288.22 0.018196111 293.26 0.021838514 298.2 0.026852153 303.07 0.032794157 308.15 0.040592546 313 0.050666982 318.1 0.064596043 323.1 0.083519503 x2=0.4 283.1 0.020213298 288.15 0.024960264 293.25 0.029919577 298.3 0.035835999 303.2 0.042800335 308.25 0.052897495 313.17 0.065920884 318.21 0.082267938 323.2 0.106310795 x2=0.45 283.09 0.022163167 288.1 0.028017074 293.1 0.034843674 298.25 0.0420929 303.15 0.051962363 308.15 0.065005896 313.2 0.08215546 318.25 0.104250344 323.1 0.137575799 x2=0.5 283.15 0.027972969 288.13 0.03399462 293.05 0.040605327 298.01 0.048112851 303.1 0.058477267 308.05 0.072571638 313.2 0.09047995 318.25 0.114083364 323.15 0.147697283 下面是我的程序 VarConstant x2=; Constant R=8.314,V1=140.118,V2=97.9,V3=108.1053,H=23100,Tm=349.8; Parameters g12,g21,g31,g13,g32,g23; Variable T,x1; constStr a21=V1/V2*exp(-g21/R/T),a12=V2/V1*exp(-g12/R/T),a31=V1/V3*exp(-g31/R/T),a13=V3/V1*exp(-g13/R/T),a23=V3/V2*exp(-g23/R/T),a32=V2/V3*exp(-g32/R/T); constStr lnr=1-ln(x1+a12*x2+a13*(1-x1-x2))-x1/(x1+a12*x2+a13*(1-x1-x2))-a21*x2/(x2+a21*x1+a23(1-x1-x2))-a31*(1-x1-x2)/((1-x1-x2)+a31*x1+a32*x2); constStr x1+x2+x3=1; Function x1=exp(-H/R*(1/T-1/Tm)-lnr); Data; //T,x1 283.3 3.02194E-05 288.31 0.000157704 293.3 0.0004041 298.25 0.000686379 303.15 0.001026437 308.18 0.001405125 313.05 0.001810204 318.1 0.002726428 323.05 0.003769643 data; 283.13 0.000349201 288.15 0.00046566 293.05 0.000981802 298.1 0.001500779 303.23 0.002474932 308.3 0.003373367 313.35 0.004337507 318.35 0.005850682 323.3 0.007716984 data; 283.2 0.000828269 288.23 0.001597956 293.16 0.002383646 298.24 0.003548518 303.2 0.004900833 308.15 0.006296129 313.15 0.008064994 318.2 0.010216362 323.2 0.013075324 data; 283.15 0.002970424 288.15 0.003905376 293.06 0.005280804 298.15 0.006616207 303.15 0.008243561 308.15 0.010303864 313.2 0.013119118 318.15 0.016480102 323.07 0.020616259 data; 283.1 0.005204799 288.15 0.006279544 293.05 0.007921776 298.05 0.00986987 303.15 0.012566941 308.15 0.015421326 313.18 0.019498661 318.25 0.024651302 323.15 0.030872327 data; 283.22 0.008048279 288.27 0.009556376 293.3 0.011981522 298.2 0.014520227 303.2 0.01789219 308.1 0.022389299 313.05 0.027783872 318.07 0.035092411 323.1 0.044652622 data; 283.05 0.011149853 288.05 0.013651819 293.06 0.016540268 298.07 0.020332495 303.1 0.024782096 308.23 0.030544405 313.25 0.037978064 318.35 0.048406476 323.2 0.061434365 data; 283.15 0.015106036 288.22 0.018196111 293.26 0.021838514 298.2 0.026852153 303.07 0.032794157 308.15 0.040592546 313 0.050666982 318.1 0.064596043 323.1 0.083519503 data; 283.1 0.020213298 288.15 0.024960264 293.25 0.029919577 298.3 0.035835999 303.2 0.042800335 308.25 0.052897495 313.17 0.065920884 318.21 0.082267938 323.2 0.106310795 data; 283.09 0.022163167 288.1 0.028017074 293.1 0.034843674 298.25 0.0420929 303.15 0.051962363 308.15 0.065005896 313.2 0.08215546 318.25 0.104250344 323.1 0.137575799 data; 283.15 0.027972969 288.13 0.03399462 293.05 0.040605327 298.01 0.048112851 303.1 0.058477267 308.05 0.072571638 313.2 0.09047995 318.25 0.114083364 323.15 0.147697283 不知道哪里出現(xiàn)問題了，請求指教！ |
回復(fù)此樓» 猜你喜歡
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有182人回復(fù)
-
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有28人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有22人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有0人回復(fù)
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 自己這一段時間過柱子的總結(jié)
已經(jīng)有11人回復(fù)
月只藍
無蟲
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計算模擬區(qū)
3樓2014-07-14 15:33:47
2樓2014-07-14 15:08:53
dingd
主管區(qū)長
- 計算強帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
4樓2014-07-14 19:32:07
|
大神您好，跪求指導(dǎo)，是這樣的：我做的是固體物質(zhì)A在不同比例的兩種混合溶劑（B和C）中的溶解度數(shù)據(jù)回歸 x1是溶液中溶質(zhì)的摩爾分率，x2、x3是溶液中兩種溶劑B和C的摩爾分率，即有關(guān)系x1+x2+x3=1，其實我上邊列出的那一大列數(shù)據(jù)中的x2應(yīng)該改為xB，是溶劑B在混合溶劑B+C中所占比例，xB+xC=1，因為隨著溶質(zhì)A不斷加入溶液中，x1、x2、x3都在變化，不是固定值。 比如這一組數(shù)據(jù)， xB=0.05 T xA 283.13 0.000349201 288.15 0.00046566 293.05 0.000981802 298.1 0.001500779 303.23 0.002474932 308.3 0.003373367 313.35 0.004337507 318.35 0.005850682 323.3 0.007716984 是指在混合溶劑中B摩爾分率為0.05時，A在混合溶劑B+C中摩爾溶解度x1隨溫度變化。 單對此組數(shù)據(jù)寫程序還好，但是這樣的話每個組成下都會擬合出6個參數(shù)，上面的10組數(shù)據(jù)就會出現(xiàn)10組參數(shù)，沒法統(tǒng)一，我想整體擬合，最終只有6個參數(shù)，請求指導(dǎo) |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 溶解度模型.JPG
2014-07-15 11:18:03, 46.48 K
5樓2014-07-15 11:15:29
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0805求調(diào)劑 +5 | 是水分 2026-03-31 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 求調(diào)劑，一志愿北林食品與營養(yǎng)095500，301分，已過六級，有科研經(jīng)歷 +3 | 快樂儲蓄罐 2026-03-31 | 3/150 |
|
|
[考研] 070300化學(xué)354求調(diào)劑 +15 | 101次希望 2026-03-28 | 15/750 |
|
|
|
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +8 | 睿08 2026-03-27 | 8/400 |
|
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +7 | 朵朵- 2026-03-26 | 9/450 |
|
|
|
[考研] 266分，求材料冶金能源化工等調(diào)劑 +8 | 哇呼哼呼哼 2026-03-27 | 10/500 |
|
|
|
[考研] 313求調(diào)劑 +6 | 賣個關(guān)子吧 2026-03-31 | 6/300 |
|
|
|
[考研] 274求調(diào)劑 +6 | xiao愛同學(xué) 2026-03-30 | 6/300 |
|
|
|
[考研] 08工科，295，接受跨專業(yè)調(diào)劑 +6 | lmnlzy 2026-03-30 | 6/300 |
|
|
|
[考研] 一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑 +8 | 披星河 2026-03-30 | 8/400 |
|
|
|
[考研] 279求調(diào)劑 +12 | j的立方 2026-03-29 | 12/600 |
|
|
|
[考研] 材料化工340求調(diào)劑 +3 | jhx777 2026-03-30 | 3/150 |
|
|
|
[考研] 329求調(diào)劑 +8 | 星野? 2026-03-26 | 8/400 |
|
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +7 | 爭取九點睡 2026-03-28 | 8/400 |
|
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|