From e7b5b1f981278c8e67766b8357f116cf70ceca9b Mon Sep 17 00:00:00 2001 From: Even Solbraa Date: Thu, 16 May 2019 17:03:44 +0200 Subject: [PATCH] added freezingCO2 script --- derby.log | 13 +++++++++++++ example/freezingCO2.asv | 10 ++++++++++ example/freezingCO2.m | 8 ++++++++ example/livetest.mlx | Bin 3695 -> 11338 bytes example/testUMRPhaseEnvelope.m | 2 +- 5 files changed, 32 insertions(+), 1 deletion(-) create mode 100644 example/freezingCO2.asv create mode 100644 example/freezingCO2.m diff --git a/derby.log b/derby.log index e69de29..303c389 100644 --- a/derby.log +++ b/derby.log @@ -0,0 +1,13 @@ +---------------------------------------------------------------- +Thu May 16 10:57:58 CEST 2019: +Booting Derby (version The Apache Software Foundation - Apache Derby - 10.12.1.1 - (1704137)) instance a816c00e-016a-bfdc-b7e3-00000448ab50 +on database directory classpath:data/neqsimthermodatabase in READ ONLY mode with class loader com.mathworks.jmi.CustomURLClassLoader@5b6813df. +Loaded from file:/C:/Users/esol/OneDrive%20-%20Equinor/programming/neqsimmatlab/ext/NeqSim.jar. +java.vendor=Oracle Corporation +java.runtime.version=1.8.0_121-b13 +user.dir=C:\Users\esol\Documents\MATLAB +os.name=Windows 10 +os.arch=amd64 +os.version=10.0 +derby.system.home=null +Database Class Loader started - derby.database.classpath='' diff --git a/example/freezingCO2.asv b/example/freezingCO2.asv new file mode 100644 index 0000000..308d4e6 --- /dev/null +++ b/example/freezingCO2.asv @@ -0,0 +1,10 @@ +system1 = neqsim.thermo.system.SystemSrkEos(250.0, 20.0); % Kelvin / bara +system1.addComponent('methane', 90.0); +system1.addComponent('CO2', 10.0); +system1.setMixingRule(2); +system1.setSolidPhaseCheck('CO2'); +freezt(system1); +system1 + +disp(['freezing temperature ',num2str(system1.getTemperature()-273.15)]) + diff --git a/example/freezingCO2.m b/example/freezingCO2.m new file mode 100644 index 0000000..8c0a8d5 --- /dev/null +++ b/example/freezingCO2.m @@ -0,0 +1,8 @@ +system1 = thermo('srk', 250.0, 20.0); % EoS / Temperature [Kelvin] / Pressure [bara] +system1.addComponent('methane', 90.0); +system1.addComponent('CO2', 10.0); +system1.setMixingRule('classic'); +system1.setSolidPhaseCheck('CO2'); +freezt(system1); +system1 +disp(['freezing temperature ',num2str(system1.getTemperature()-273.15)]) \ No newline at end of file diff --git a/example/livetest.mlx b/example/livetest.mlx index 48edbafcbaa0b90c78461d4c62bd43c13fb2a081..4c1933ab6448865dbceb97faefc8893127ba519f 100644 GIT binary patch delta 9068 zcmV-yBa__k9LiX*kOBirg|M@e0{a3EN`U z6-EJ#${JQC_iZZ64Xz?Zikm2r(qKMt7x#S;fBgiX(^Z+>?(_hIX=FdxepD!ydr#cu z++g_Uf7|Nn(~=Bh-Zp>#6Z#YKr%y=}wQ=5LfB)0Xm49UZ^lzVkPU>~?C7>E%zn}*EPHs{q*>{I%Ogs*;cZ~Xo z!pN*K^|sJPZf}3TF6aFzO|HRFm1<{?z~|J3yrDtY2HN;WIptZu`>+#l*g>ijpC--798w1y}?Psr;bwX0DkU42hrw1*HgS@4m+@4oxmVHU^d z^)>*dU7 z(frm;me()t1CuibPF*)j-4+Q8it6WvKj&-$rzbs3BWX#~6lSZ1Vv~VMz2UwT)jdw~ zfSQN?b<+l!h2W6DA{LTgNUYqj2;j8dE)BMdJC1)VmV-}3;_fs_p8T|NtNmQGv`HD) zK607{M=&xRaU5Uwd(S4y9gTfQa*&u#lmV>W%nTbpL)R&BPUnfYN4%^wd^}jiLtmj% zY*5NHHS4#mWi4gdG)0~VK(0cA&>G#f_IX{ zxo=FO%a9e|%Z4f92AK=I>kF|_i%&mk3>^Yu*_c!s;Hp^`%P<=erQ@#}2px2758)Ctpb!ATZON?d=8VDUG`kz(+J;t@b7j*{WlKXTR6tM@Z@ z2liUSPkdu6IT){EsB<^5lLeYft$o{IrSB=$gp%ELcgX00mRhdnXivr`2!dgly3y2f z59<9n_IF^9*UaaFQ$I<#c^>c!*H;mDs_*?IhZwu0R`Ggm$_72d>j{GxDTXO@h;o0m z)(5mlz)qU%1xl@kfZQ_zpTaPWs}bTN*i}^mF0F@CJ??a;B}23eOpOBpO0 z>MYBQE#_G_F-k0j0P@8%-1yoB>OJ3YyQXQ%dRYQ=VD2uxEo z%0~ypMLWL)&HO{{+hxC6U&bkDr)Qw~9B@=Q?Scy?mhxU_kmguV?Q6(5g8;S)dTaGp zaXqGZ^6d4NYOo3doKOb6bIoB~dCivGILoe-L)Txx5Kc8j*()9l(xcBP`qHMg%k1d9 z{ov$MM9+8TquHg}1-hF&_uzlY(+l9*VR3*uX_UK%oQwg6i5SfW6F|?sf$%BTj}+P- z00X(_l;d-t0Xc}%D_GnG?{@qqb=I*C*}t4Lsal`(f7Ql2snaAVY2i06yPr7WQKHDw z5KWTo`vVhN>1f*`-ZpMQKl!hgqI>4WsQwYV>x&PyjM?jN^ptTY?i7FI&?G#-SHk=_ zKJwR-cG6w5PNj{clJia(pV4%tAO1)iPYkT38ozDooJN8`_*JsJM2S!3hV4dV~ zoTg*YdROihV;@FxtMh+IL=J3azE2-!q4%g>o~#>pcD<{1^Ij|w$6gg@$2vkbRcugQ z*V`$R)!?M{)FhBE4RD9*LaJR7K4Mcn3`2VY;7%{)g(9M^43mbMp;fm!TptK zKl*hZyBgvUEYN<8MHu^RpIXw8iAV!^v>KD8y1+`NCu3B<3GRQu>DntB)?8t{i-nsi z@vX%L8x^oIYDY&}eZ>)B^yoJ2uTiySJCLbR&mUlKIc-Rl1n>%d42nbNhe~<@OAy*M z13F#wfKAmWy&?CC<3wAK(4s_Xn#zw2<%BWWs6+*O!nMjJsdI#J^;M9G>4lgd^Kora zVt#h)gL>O8vIBp@P+NX1tX_dJpxBbW`5`c`(w8={7jo#l@42wn9I2wWpix>p27hUU zq$$Nd3+>7Y*Ck4f*@5SNQDdQ8v!iVaaq(7!V|50OQwkftbXOlxI&PeH@>Wqn!pUk` zOzFEsWWE*4L|Yxs4mi@XjiL;sZ5#BRe_T1K+Q@Zuj(vYTpF;nQrTKbayg4;1@&H_V zod|~V;bz215xGVEL*NkZ`ZYH7&B@OXfTTi#YfA)~9C|)5+EHx~*SPCN);chktlhy} zr&p(vez2WP1399XAzt}5=p{!=0?gJtB;zJ+KRpQi7;-&inhxC>uKM!q;7h@|hge=o z78H}mI{JS;A9}yBYE;5ze_iM8WjoxJmeEaq0{}()Qm^sl9vv%AEBkOF?4}D3!F!mq z3f%>syYl$dQV$)(wanO`X01V?An~Ozf&F<)0)TuD!F&>+ubfyuK&%^aU&=BNj!h?j zXjvdRCkxxVoJ)Z}89Y)1SOHI1$ZUj`-!+~KgzkS_GjSYVmv=J~XQ}6owvfj*U^>o2 zqCkRCp3g0sef0(MX$%*%(d-CKCs;@;tx&P_G^hPtL^=Jr>BuYnp15 z+OAlZ<)jmx;2U{1EO%*`-A&WN095dp5RD~Azv2^6n0zKjJoid`IyXP>M<9F2GHOc^(!DOg>;ZL4;YUovdn84vZfs%!t9{_(d zG+35}Q>JTBuSF-ol`eg{aQ7@w$$J1_N*L~wNhO~$Ve1#PM?>|k)=6!d&CaA}1WF2H znr93(@zq=n>ye%af;L$Yz3cN+KGA8XVUzwO0Upg`#^WUh#CS^@w}Y&%BO!CZ_e#Sp zN@+!z+Qz7OfLxjF3&m8l#=EzV8VY~pbVpkTfTFi??X;04klCwHypLS&BE}Py$P?f@IeZv++9Q7WzpZexTlrY4&P{ERf0aavs zI#R4?1u3=?^t?sV%1-_>#czKnK~P>92Wp@t67AEg?-QXHtAhYtGz{bk98ZlRF?wJL!^u@^*i*F)=zg`s3JN{%Sv0wKW01ZBHB{A>8xWG(J8P=UL=E_Ss>* zZy)n!T^9L;JKY*v5G9Ux6&f(-{3UYsBfHq>+AamsT|vGv?fLSI^*21@=K{ZJe1~%v zu@rfHx%1wUt8(X@XyOG-%kz5Q7xc#*F7Ofh^pPh+mL7?v*93p$So#%CwY;wYOBHZb z*~N*GGC*6!9^*LYvh*)A;Uu01wqjh{CC0CFbeH1vAQRr})nqwG9d(0%_0pZZ`7{6K zvW*kAzC(o|0Yoez6fKc5zO{W!Y;l?C`GFPbxfT0I0R|!`iIG-w1H9>kl+u%6ZZ@|| zZ!mx6lq`cs<2Zl+@|oul<4$i2ZylsdY-m_75S}IQtSJxQeWQ9YbRay7PeRRNQ64oy-LLI5)ws z#-0Tb!5G3DTKUuZ;#`MU={G^Ho(Eh@94~8atxUFD#Or^k-lS1|OEQ`zT()zUmz+GT z@^DVwR&t&^Xp^vRgzb4pu`{YyNyaj>BY}~ZS+Q`x#73DXj`3;>eE+Cvg`Tl_>uG1_ z2X9I;lY<)iviEl_pLgq!7)wu3O zad#H*Zb6{!)?G`^y_iK|h`}Mp_4{t8fo9=wdhKu|<+op=vxD(@HyB_(mw45zF0sFB}a*MY9 zTRVTfx&ChVOe`x3f|IoR2*dU@3^anNxjOrLW_39R9Tm%oGSB{4$9HlVK0&;(|pwQxA>=+L8{=PO28sEI?);1GH12XhN z=(aZtm4taN8?lRgHT)|JtNO;h9T$@g?RmkHe55bC*4?$#3IWG7(*=75!;qjb9Pfs0Vbd6Q@ok3zm!qNav8FY-fc>YHfwZAKh2t9d4a8Xcz5{>m zV^T9A*x~(=fK{?yCcYt+p{d5d%%r?I`jqeH^kqaHS}Uy8*#nHw2hWlc*@an2Lxb ztlvV|)HiNePeTp^a`xO#0`pcX>dJrT>e+jIITusJjfx=aj-n{y$AF~8g_doRl}^wkCM%hYvev+sLMpRmE- zVDOC#*R$8-rz%Iox~QZFPXza8V)gV5*m(qEjVOzvfO`4yE@3Iwm4OSjuqJ=fJ$Np_ zCcr`B4=0Q=z7p6kBYs!mR1G-p1}y~a!v z&qBR5AbW#ca(Cj)Y5g&?0hxhbC=>zuyS**G4k|9Q%6ET9R!KVwVn8zUe(r43<3WDV z#O_?hu|K`ZOG1zS*}_Xzln#>O`hIs?MG*TqTPt|(hG!Zlzcem;uQY$(-7A7X``zui zaysIuR|Ks-%Ag25@lWGS2nX-G%NTb)&ex=AW~vJZ2kpFyfS;;dMW8Cx%4O89f#(=- zMqS(ixl$EJ2B84f&;;-dYZ;ea##&msO#i4yce2TlQjB=M%l&0k986(Hlc$*H9l6lZ za!&RacLlk8J#Uv~isgS{ooxe%RU5L}IU!xFxU$OPeYuj{d*ie0RK?gF=kxk%&wzJc zR6pZ-hXWidCz_Ne)0b_y$2rSd`+-miuE=edIxh-uXqPk*oes&Kk7kvGvw1J24%DWZ zbAFH6OKogUJO^i9ego|6r)Dc0pO*ob*;=T``%KGw@#JOu1bBa%EW5>1mZiYpWT6D< z_MtsfU26|c%*3J>18y&PVAoyVmd+_QyRRqDU!QMU^_QB#Ad<+c;Gk}aL2 z;yCa3xtr(%UOZuAmbLoPKI?r-h?d|AHYX%Hr#-F@O^^@yFCs{@Rd1Kb4N!g~YtSa=eUpRmEY-Yi!Kl(&(!0;p@*eH;v06)9x_qtzSnF8v8KtJpK8)I3FTot~{ zjv680^t%iq-8nkWV$Tz(-+;2~cmaVqGjGj$oT+*UAMgav0~IIxmMLIuUFCZLIo?U- z50wRdjmHc09q;nUc`d_ntYL??+=gHeT}=fHwr3I+m-<5l1{=In3O88F$co2(q*W#@ff& zHyEpWb+6lEb-oC;V%$KLU_Z)C?J-fRzUfg&~UX1je zNp!PtYfFDp1qRqy88+w4bZz348RBDAbxyooAE%r;YB#4J;M00Zi_zcOq`$3}o}rEJ zz_18vvq-ho(S67-tOCAOCPO1XEeF~Hm5f3XcjzITdFed4Orod-6T2V zgT@xcXP8D#3%FRx>k`ImfpX5~^y+!U(d8)cZa7B!49L5j=A1KO3{aR8Zka8(rjfom zls|Fhmiw2{;<)qWjbA2f>mm{M$0&|yl${9u(4yBx3Ztmj@q~cT# z1Vew$^Um{Zd)#r(56qF{GVwT$A8*r9`)$!+mm_EJ5-8Pk%cqyu@IJ_Bg^(eBywD}f z`TzgBPT9!!9fnhb6g@tyt0$H3BgP<=eoGwx{*QA7p7;j0ni|+>W_l)w1b&Vs5_*;y z+^b$UHA*9*-CU;g%lFFenYq=Ql!_78%e{Yi5$tK0aGD@M9ut(rX3pF>S7Z=pTBGS% ze?ZtJJZTifNEI6U_fvO2=`^5b+k;6=$}HNqG7&)Cr~*l1vA)BJhAz+a21To>!M&$e z`E}I=#h&NIj9=D4zj+b*g*)NW`6*Mje)%WQxj-`%1OcQ+vv2sYh`>-(SyeMVj(2|> zl}0FCMJ*1xgy^Zpr_h%n>qWTvMc~X}h&j0?qy7#ky5Hd|W)(!Dp}}o_#wo?@LZBfI4iF6Qwi*xi6CWg*qfQp&$6& zDzaq1^9~GPee^AE*#npcKXvgW0K0#@oWFY9h1L>=G+gRd@U`__$^wCtz$_2ydHz1hh)!lLQ zZn+41xpp@|+&6-VnFOR(=s^%e2?_$9gRQ->2xju=P?~+KA?AJ<2*UUAZF@B6)ej&Q zOA?RjT6D1hU%?fEF(o2ObP0dx<@$rEY78^5sTb(h4V|m2gVw>j?Q!?%s&ykyxW3GV z1^2}u_tUegy=MxuI1I&M{3b5xVuU!cGcM2#rty)4Ay$h?(UIp+Z`dnc6xk!IBYepgV4y}Of z!VU=Vtd28I;CWWs@qh!{9Y#MV4D-{D?&=GJ3?$IU?K`kHSUyB(>D>Ekk|d>rSb-&$em+wvYfU$_#;QaE z{uayUxP{7F=pnVns+4<+tX*)X(8G^n?0#G#@EaUO{i5rhX`Zk!?N=L)`!|P4OZwWe z3{#d26WQ)9sBVvX8+DTc#SlRp_Oa(?ih29)jTx|gG{eIV9X)?gCbktI==-*@N{Vj~ z9Yc_rE3U6$IyMI$^?XTAvE!JhUr``#bS|pryB;xw-mZSn>Rq^F5z|T@R1y0;P&YVH<2Y=j8u_LQl335}c z{pr{Dhpp&2n=F5PW?89YTh_j=w7d~aaNf~#G7w=0jNx5Z2t|+871;5U?EH69AwB7W zz3&v@oLWA_26BC&nfO`>={({gL4WwW=_)Ia!z&Mibhg6M$Gr&V+TVjEM>JWKoYR&r zIAtou`q04RPULD2a2qTiZ&hx(w_JD2aYKyFc5gfKnRke)yk`D6yGqtA(ynyJSC zf9eciyl;p>x@rlYJ9Xax+9r8zDi9kqBz|@tUg6S^xF8onZ^K|g;?8#+o}KXd2OM6YG{FDo|SJ=@fHIORxq z0z4P8N*JCm@ZQiYg~E+72P9-(U6($CJ$_ceKgpdoEDJ^vqJux`QSTX4GfF;r}eHuA~h%`ALLk+;%7VN+O&I)>``eah=QM(byd+?7}5j@9dUdHs&vz` zFLZy`z-0GcA{b68YecZnsBe~K@!V85Za43W1(Tbbn{fj)U;z_l9QJwKh+;=h6fq~9 zK*aGIm}Q^a$!ua^^Y>dkUJL`>u^pya=p94n$w%w)ko&#AgOJX5T1S8muY55y7<`xS zs1bIzGJe)$uT2PbkDg*;FsTBlx!sQ*5UYP~^e_P#9z6{p0N(rRwi`2r)>BiZh4L_T?JKD}fWq-=k#V4l zqkuY*D|0|mL8*7VAIQvu+wxqL{TGc{h!oetJPk!PV(o*nf&%JpIvOcFRQ2y)?^J&t zHyeMse}hoi=e-I2lcq-`uXQEiH>Wq{E!6Y=V_RepW)v${hTzhnO?s%eNKaF1U-J!_RCt= zX~>)Gb2qe6GLB}K3?KJv$Dg7q2?l?zs=nSe{iWk??y7!9etQ!YsPl_>Q!q?A(D*Y3 z1pCK_L#nT>Umt(DL8~@#vL(NT68`*+{rMIEzv%b>hyc(ZZ1tzVf1;6}kiY$|-spGr z#=ogI{>RgMO6reI2dbZD00GedNAdvR)2Dy^6ulDgC;nY5upGNz#Q}w;Hm`obX#6LXg#JB>fBp2;p8oRtmV|#7 zJ^GvI(SKOFSzlR3-T{PHlqQU-mK?8qeo$eoSpx^Nn>bq#5$T#6Qn*V#?A(+@du=st`puZ0W z`fV`K|1Myx!ITUD4)iA^`4PpwV;#)*0r@r+LH~eU&;EchF>QxWInW#ZpD6h+pTN&|vIa;R)L;Jg(eq!w`q58WGF{s$kL7l#Kk8if zpXT(#pI`47JSM=QKK*~|JSLx-pavTM()tOk|1tT;lyS07U!C!5_vw|7pVGPckp9!u zeu{!BnyX`MmG^TVe;I#_?b^JVJ|_B=J^zsZ<#6<0)B>XZTDZUb_@&bLqS047_zy#@ z*O4fYzZ}p10hN!|?T5eo4>Y0v_AdbI-#4G{lX2|ZRW9~l9ydw8oPhkE<5cbR!>_RZ z>pvj;UiADv3CF)3Kk)y%`0-EGe{)>)n*qN$F#6RQ(b4|z{ulNC0Z>Z=1QY-O00;m? zcZg16wF|ozlVAxme`RdJ2H1sWSDKiJ(RkOCGuB8;(^lF3y&Vp=%w2eS@;>Lu`6$+V z)kWW6o$PgqXhKi~TJ(+7ZHZp%O`N01d9EAY^%_ca0EeP=RS3gGZ{gYY2CSFhq9D?a z2~(n-_l9BY#10g9Ng#9xx4l)|htjs#aPi4oz$qastibbze|wIP!r074bW%2AR+=x{ zje@lfVrcv$+aH0s~5gu#ycv4oZcvPAY^;wsau? z0DV1^Obs<0Vzmpq76AYN5CQ-I8vpJ?H${vMq}zkf56?&A#}~cfXwN zbNm;ZG$B7hs#x&j9(Bh(0z4E^%rB37pU+=SYE4MOdrMWgQ*=Fg3BnR6Yd7}K?=cE6l3i8*}85u<-p_jpWExpsck7#=j` zjlct6!X@JEq+mb?zCx#QIv4$Q}vTfWN4dg=b}P(8}Uq_gJW{LL3wkaQLC()YWT{sc=iVyDa`K<@eK%Cm4_x zapHLMiA$%?Zd2SRbv_NcGihR_H*cjk-&uMti=tDJW`cj?WKJJXfi>gM-yge+gTWE$ zq;X)Sao9@Z@UhUSVGvwlQyYz`mByl##^SNiz=J8f!a*%7WoD&8$E`Sw_;K(!J@c+u zbf?DL4uiJCppON{FRU{q;=m#ohSalU$Dtiad$mFzR<8UPAR%gNFcs+qQ{dQE$imKG z+Rk9Qy9|F6==bpkC7<#HoP`AewHK{)*3MwIdkkXJ0G47wqKbx7QQ7ZnP`Z@$ZAk3J zDs4YU=nftvjY12le4*2?mgw!}sn@Dq86hwyVJ*BjC<|3T<@Rn8|6om zu7xv`2<mCo_KqS+D_iq1lxtCSo++HRX&o($ch5wtsJjgDtuXFHhd*JUJibdawHE3#^la zt`JQKia?8@mAb3YOS6d!6gkgz%lkn?g%03Qw605Gm>4WP*};JI5?mBS+A(1&wDaCD zjGfql;w}k<4&ipNiu+L74jV2$cn3HoWQ7%Y-g1A>@lhC?#fVPIRxC>MY5Pgg3Jg6c z(B5I1P<&>U^alf1JLd3$_F$E}3Qx&G?wwqyaU7E|o47DP#cz)f&(o^0)W=f^Ky`Ui zGhu;yXrmxtKZDK5Yj$_vY*3vNdK;6Qm=q07Sza(wyqEY_|5>QQZ?^Iu#ylpghEkSg zjN}d%##vC8_|FsTZ%|7E1e2R98ncurL;?dMY^IajDn1S(Y^F{jN+5EY0{{U36O$(^ yH3ffKSoal^VJj5^O$n2DD@OvG3X{PrNdlV-lff$+ldUNYlQb+426qeq0000XP;t)y diff --git a/example/testUMRPhaseEnvelope.m b/example/testUMRPhaseEnvelope.m index 63c4924..00a8cbd 100644 --- a/example/testUMRPhaseEnvelope.m +++ b/example/testUMRPhaseEnvelope.m @@ -1,5 +1,5 @@ fluid1 = thermo('UMR-PRU-EoS',273.15+55, 25.0); -fluid1.addComponent('methane', 6.0); +fluid1.addComponent('ethane', 26.0); fluid1.addComponent('n-hexane', 44.5); fluid1.createDatabase(1); fluid1.setMixingRule('UNIFAC_UMRPRU');