• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    高溫高壓CO2狀態(tài)方程研究

    2016-05-09 06:37:34郭向利龍新平
    含能材料 2016年5期
    關(guān)鍵詞:狀態(tài)方程熱力學(xué)偏差

    韓 勇, 郭向利, 龍新平

    (1. 中國工程物理研究院化工材料研究所, 四川 綿陽 621999; 2. 中國工程物理研究院, 四川 綿陽 621999)

    1 引 言

    炸藥爆轟時,爆轟產(chǎn)物處于高溫高壓狀態(tài),爆轟產(chǎn)物壓力由數(shù)十吉帕逐漸衰減至兆帕量級,溫度也由數(shù)千開爾文逐漸衰減至1000K甚至更低,建立能夠準(zhǔn)確描述爆轟產(chǎn)物氣體高溫、高壓的狀態(tài)方程對有效表征爆轟產(chǎn)物狀態(tài)變化過程具有十分重要的意義。

    目前,實現(xiàn)爆轟產(chǎn)物氣體所處的高溫、高壓試驗條件十分困難,關(guān)于高溫、高壓氣體狀態(tài)的實驗數(shù)據(jù)很少。分子動力學(xué)(MD)/蒙特卡洛方法(MC)是研究爆轟產(chǎn)物氣體高溫高壓熱力學(xué)性質(zhì)的有力工具,它們通過有限的分子數(shù)目,結(jié)合邊界的約束處理,能夠準(zhǔn)確計算氣體的宏觀熱力學(xué)性質(zhì)。然而MD方法和MC方法計算耗時,且無法直接應(yīng)用。解析形式的狀態(tài)方程則可以直接應(yīng)用于實際問題,形成子程序直接納入相關(guān)程序中應(yīng)用。在描述爆轟產(chǎn)物氣體方面,已有BKW[1]、VLW[2]、LJD[3]、JCZ[4]等解析狀態(tài)方程形式。其中,VLW狀態(tài)方程為我國吳雄教授借鑒相似理論,認(rèn)為各階維里系數(shù)在高溫下是相似的,高階維里系數(shù)可以通過二階維里系數(shù)求得,進而將維里物態(tài)方程以一種簡化形式寫出,形成了VLW爆轟產(chǎn)物狀態(tài)方程[5-6],其成功應(yīng)用于炸藥爆轟性能的理論計算。與BKW狀態(tài)方程相比,VLW狀態(tài)方程在描述爆轟環(huán)境下氣體高溫高壓熱力學(xué)狀態(tài)時有了一定改善。然而VLW狀態(tài)方程形式仍然存在一些不足,VLW狀態(tài)方程不能有效描述爆轟產(chǎn)物氣體組分的高溫高壓熱力學(xué)狀態(tài),從而使其所采用的爆轟產(chǎn)物勢參數(shù)缺乏必要的物理基礎(chǔ)支撐。作者[7]前期應(yīng)用VLW狀態(tài)方程計算H2O沖擊Hugoniot曲線的結(jié)果表明,即使在優(yōu)化后的勢參數(shù)條件下,VLW狀態(tài)方程仍然不能在寬的壓力范圍有效描述H2O的高溫高壓狀態(tài)。

    針對VLW狀態(tài)方程存在的上述不足,本研究提出了一種描述爆轟環(huán)境下高溫、高壓氣體的對比態(tài)維里型狀態(tài)方程VHL(Viral-Han-Long),該狀態(tài)方程基于LJ勢能函數(shù),通過對炸藥爆轟產(chǎn)物中的重要氣體組分CO2高溫高壓熱力學(xué)狀態(tài)的描述,表明對比態(tài)VHL狀態(tài)方程能夠有效地描述爆轟環(huán)境下CO2氣體的壓力、體積和溫度(pVT)熱力學(xué)關(guān)系。

    2 對比態(tài)VHL狀態(tài)方程

    理論上,任何氣體的狀態(tài)方程,都可以用維里形式描述。然而,實踐中隨著維里系數(shù)階數(shù)的提高,計算的復(fù)雜性迅速增大。Barker等[8]基于LJ勢能函數(shù),針對無量綱第三階、第四階、第五階系數(shù),進行了精確的理論計算,其以表格形式列出了特定溫度時的系數(shù)值,不適宜實際使用。在前期研究中,作者等[9]提出用一種簡化維里型狀態(tài)方程形式描述高溫甲烷氣體的熱力學(xué)pVT狀態(tài)。為有效描述爆轟環(huán)境下其它氣體產(chǎn)物的pVT熱力學(xué)關(guān)系,本研究提出了一種基于對比態(tài)原理的維里型狀態(tài)方程形式VHL,表達(dá)式見(1)。

    (1)

    C*=c1T*c2+c3T*c4+c5T*c6+c7T*c8+c9T*c10

    D*=d1T*d2+d3T*d4+d5T*d6+d7T*d8+d9T*d10

    E*=e1T*e2+e3T*e4

    b0=2πNσ3/3

    式中,p為氣體實際壓力,GPa;V為摩爾體積,cm3·mol-1;T為溫度,K;R為摩爾氣體常數(shù)8.3145,J·mol-1·K-1;N為阿佛加德羅常數(shù)6.022045×1023mol-1,第二階無量綱維里系數(shù)B*采用變步長辛卜生求積法近似計算獲得[10];C*、D*分別為第三階和第四階無量綱維里系數(shù),c1~c10、d1~d10為常數(shù),通過該表達(dá)式所得第三階和第四階無量綱維里系數(shù)值與理論值十分吻合; 對于第五階以上的無量綱維里系數(shù),則采用組合函數(shù)表示,其系數(shù)值e1~e4、a、b、f由無極性分子CH4的熱力學(xué)數(shù)據(jù)[11-12](溫度1000 K以上的112組pVT數(shù)據(jù))確定; 各階維里系數(shù)的擬合常數(shù)獲得方式見文獻[9]所述,具體數(shù)值見表1。T*為無量綱溫度,bCH4為67.21,w為對比態(tài)參量,VHL狀態(tài)方程通過參量w實現(xiàn)其它氣體與甲烷狀態(tài)方程的對比。ε、σ為LJ勢參數(shù)。

    表1VHL狀態(tài)方程中各階維里系數(shù)擬合常數(shù)值

    Table1Fitting constant values of each order viral coefficient in VHL Equation of State(EOS)

    coefficientvaluecoefficientvaluecoefficientvaluec1-0.96665 d1-2.31154 e12.15701c2-4.51039 d2-7.49608 e2-2.52736c3-2.58733 d3-2.56864 e30.27080c4-0.85487 d4-1.14947 e4-1.00190c52.08033 d52.09534 a-0.00243c6-0.51631 d6-0.84757 b0.00018c72.02825 d72.81486 f3.60800c8-2.15543 d8-4.44480c9-0.12489 d9-0.30669c10-7.65161 d10-11.36945

    3 CO2氣體高溫高壓熱力學(xué)狀態(tài)的計算

    由于液態(tài)CO2樣品制備條件苛刻、難度大,其高溫、高壓的基礎(chǔ)實驗數(shù)據(jù)較少。1990年初,Schott等[13]才發(fā)表了用化爆技術(shù)在5~30 GPa區(qū)域獲得的一組試驗點。Nellis等[14]用二級輕氣炮技術(shù)在25~70 GPa區(qū)域獲得了一組試驗點。劉福生等[15]利用二級輕氣炮作沖擊加載手段,獲得了CO2在20~60 GPa區(qū)域六個Hugoniot數(shù)據(jù)點。實驗數(shù)據(jù)點偏少,且缺乏溫度的直接測量數(shù)據(jù)。因此,本研究以Belonoshko等[12]的分子動力學(xué)計算數(shù)據(jù)作為CO2的基礎(chǔ)數(shù)據(jù),溫度范圍為718~4978 K,壓力范圍為0.5116~111.078 GPa,基于VHL狀態(tài)方程,應(yīng)用復(fù)形調(diào)優(yōu)法[16]優(yōu)化了CO2的LJ勢參數(shù),勢參數(shù)值與文獻值比較見表2所示。采用VHL、VLW狀態(tài)方程計算結(jié)果和分子動力學(xué)計算值比較見表3所示,不同溫度或壓力下CO2體積計算偏差如圖1、圖2所示。由表3、圖1和圖2可得,采用VHL狀態(tài)方程計算得CO2體積平均絕對偏差為0.971%,最大偏差為4.04%, VHL狀態(tài)方程計算所得體積偏差與壓力和溫度參量無明顯的相關(guān)性。VLW狀態(tài)方程計算所得體積偏差則與溫度具有明顯的相關(guān)性,在較低溫度下,計算所得體積偏差較大,最大偏差87.149%,隨著溫度的升高,計算體積偏差逐漸減小,但仍普遍高于VHL狀態(tài)方程計算結(jié)果,采用VLW狀態(tài)方程計算所得平均絕對偏差20.2%。

    表2本研究所采用CO2勢參數(shù)值與文獻值比較

    Table2Comparison of the potential parameter values of CO2used in this paper and literature ones

    (ε/k)/Kb0/mL·mol-1reference247.063.37[17]205.085.05[5]181.867.30thispaper

    圖1不同壓力下VHL狀態(tài)方程、VLW狀態(tài)方程預(yù)測CO2的體積誤差百分比

    Fig.1Error percentage of volume predicted by VHL EOS and VLW EOS at different pressure

    圖2不同溫度下VHL狀態(tài)方程、VLW狀態(tài)方程預(yù)測CO2的體積誤差百分比

    Fig.2Error percentage of volume predicted by VHL EOS and VLW EOS at different temperature

    為驗證對比態(tài)VHL狀態(tài)方程在更低壓力下的有效性,同時與VLW狀態(tài)方程計算結(jié)果比較,本研究引

    用NIST數(shù)據(jù)庫的數(shù)據(jù)[18],對CO2在1000 K,20~800 MPa的熱力學(xué)狀態(tài)進行了計算,結(jié)果見表4所示。采用對比態(tài)VHL狀態(tài)方程計算得CO2體積絕對平均偏差為0.698%,采用VLW狀態(tài)方程計算所得體積絕對平均偏差為11.988%。在固定溫度1000 K條件下,隨著壓力的增加,VLW計算所得體積偏差逐漸增大,在壓力為800 MPa時,體積偏差最大達(dá)-19.771%。其原因可能與VLW狀態(tài)方程形式的高階維里系數(shù)過度簡化有關(guān),隨著壓力增大,描述多個氣體分子同時相互作用的高階維里系數(shù)的準(zhǔn)確性要求提高,而VLW狀態(tài)方程中,除第二階維里系數(shù)與理論值相符合外,隨著維里系數(shù)階級的增大,其與理論值的差距也逐漸增大,故CO2的體積計算結(jié)果偏差隨壓力增大而增大。而本文所提出的對比態(tài)VHL狀態(tài)方程的第三階、第四階維里系數(shù)均與理論值吻合,高階維里系數(shù)則通過甲烷高溫高壓熱力學(xué)狀態(tài)數(shù)據(jù)優(yōu)化獲得,其具有扎實物理基礎(chǔ),因此,該狀態(tài)方程能夠很好描述CO2高溫狀態(tài)下較低壓力范圍內(nèi)的熱力學(xué)狀態(tài)。

    表3VHL、VLW狀態(tài)方程計算CO2高溫pVT關(guān)系

    Table3ThepVTrelation of CO2at high temperature calculated by the means of VHL EOS and VLW EOS

    No. T/K p/GPaV/cm3·mol-1MD[12]VHLVLWerror/% VHL VLW1718.80.51164040.9118.132.275-54.6752798.80.58294040.3125.830.775-35.4253920.20.63964040.7331.431.825-21.4254963.90.68354040.3632.260.900-19.350511130.77324040.4335.221.075-11.950612300.84604040.4436.721.100-8.200712960.92084039.8536.74-0.375-8.1508720.60.77453535.9514.892.714-57.4579774.70.84463535.5320.111.514-42.54310902.00.94643535.4625.661.314-26.68611805.71.4803030.3617.631.200-41.23312999.41.7433029.9922.55-0.033-24.8331311901.9913029.8024.87-0.667-17.1001415092.2903030.0427.480.133-8.4001516752.4393030.1728.400.567-5.3331618072.6683029.8128.49-0.633-5.0331719452.8323029.7528.86-0.833-3.80018773.82.1682727.2613.790.368-49.22719887.32.3622727.0217.49-0.515-35.6042010162.4882727.1120.00-0.184-26.3622111683.6992524.7119.42-1.160-22.3202219834.9522524.7723.39-0.920-6.4402327196.0672524.7224.72-1.120-1.1202436377.0652525.0426.030.1604.1202544197.9142525.1826.720.7206.8802650608.0822525.8727.813.48011.240

    Table3continued

    No. T/K p/GPaV/cm3·mol-1MD[12]VHLVLWerror/% VHL VLW27700.84.4822222.662.7994.040-87.14928812.94.7212222.4511.923.076-45.27129913.44.9552222.3114.032.433-35.5833010245.2902222.1215.431.561-29.1553112187.9302020.1915.350.950-23.25032186410.522019.5617.35-2.200-13.25033262812.182019.6318.86-1.850-5.70034432115.052019.9220.70-0.4003.50035497816.362019.8821.00-0.6005.00036803.211.351818.238.881.278-50.68937122812.471818.0013.310.000-26.05638164113.561818.0015.070.000-16.27839197714.521817.9915.94-0.056-11.44440241415.641818.0216.770.111-6.83341275616.741817.9617.17-0.222-4.61142316117.561818.0417.680.222-1.77843346918.241818.0617.990.333-0.05644783.019.471616.147.190.875-55.08745118620.621615.8811.14-0.750-30.37546164521.921615.8812.91-0.750-19.31347200622.941615.9213.79-0.500-13.81348231824.111615.9014.31-0.625-10.56249277025.071616.0015.010.000-6.18750307226.321615.9415.27-0.375-4.56351338427.041615.9915.59-0.063-2.56252395627.991616.1216.160.7501.00053803.726.221515.097.010.600-53.28054119527.391514.8510.27-1.000-31.53355160828.631514.8311.77-1.133-21.53356206929.971514.8712.81-0.867-14.60057240030.921514.9113.36-0.600-10.93358271131.951514.9213.77-0.533-8.20059307532.891514.9714.19-0.200-5.40060340133.671515.0114.520.067-3.20061393834.741515.1015.010.6670.06762787.235.671414.126.170.857-55.90063121937.021413.849.51-1.143-32.07164165938.361413.8110.91-1.357-22.07165205039.581413.8311.71-1.214-16.35766234040.501413.8512.17-1.071-13.07167279341.971413.8912.74-0.786-9.00068319042.681413.9613.19-0.286-5.78669340144.181413.9013.29-0.714-5.07170789.649.451313.145.701.077-56.12371119250.771312.868.57-1.077-34.05472159452.071312.819.80-1.462-24.64673198453.321312.8210.57-1.385-18.69274231354.381312.8411.07-1.231-14.84675288656.211312.8711.74-1.000-9.69276307957.431312.8511.89-1.154-8.53877355158.001312.9312.34-0.538-5.07778395658.711312.9912.67-0.077-2.53879782.069.771212.215.111.750-57.425

    Table3continued

    No. T/K p/GPaV/cm3·mol-1MD[12]VHLVLWerror/% VHL VLW80120871.241211.917.86-0.750-34.52581161372.641211.848.94-1.333-25.50882199173.981211.839.61-1.417-19.91783240775.371211.8410.17-1.333-15.25084285476.901211.8610.64-1.167-11.33385323878.241211.8710.97-1.083-8.58386349879.161211.8811.16-1.000-7.00087807.0100.71111.264.962.364-54.882881209102.21110.987.12-0.182-35.264891599103.61110.908.06-0.909-26.755901971104.91110.888.66-1.091-21.273912416106.61110.879.20-1.182-16.400922866108.11110.889.63-1.091-12.500933012109.21110.879.73-1.182-11.545943326109.61110.899.98-1.000-9.236953959111.11110.9210.41-0.727-5.364

    Note: error=100×(VEOS-VMD)/VMD

    表41000 K時CO2pVT關(guān)系的VHL、VLW狀態(tài)方程計算值及與NIST數(shù)據(jù)庫數(shù)據(jù)的比較

    Table4Comparison of the calculated values of CO2pVTrelation at 1000 K by VHL EOS and VLW EOS and the data of CO2in NIST database

    No.p/GPa V/cm3·mol-1 NIST[18] VLW VHLerror/% VLW VHL10.02434.60436.9437.400.5290.64820.04228.77230.1231.500.5811.18730.06161.37161.6163.700.1431.47240.08128.35127.5130.30-0.6621.52750.10108.89107.0110.50-1.7361.43860.1296.07393.3597.31-2.8341.28970.1486.96983.5787.96-3.9081.13580.1680.15376.1980.94-4.9440.98590.1874.84270.4175.48-5.9220.847100.2070.57865.7671.08-6.8260.716110.2267.06961.9167.47-7.6920.595120.2464.12458.6864.43-8.4900.482130.2661.61355.9261.84-9.2400.374140.2859.44253.5359.60-9.9460.270150.3057.54251.4457.64-10.6040.172160.3255.86249.5855.91-11.2460.079170.3454.36447.9354.36-11.835-0.009180.3653.01746.4552.97-12.387-0.092190.3851.79845.1151.71-12.912-0.170200.4050.68743.8950.56-13.410-0.242210.4249.66942.7749.52-13.890-0.309220.4448.73241.7548.55-14.327-0.371230.4647.86640.8047.66-14.762-0.430240.4847.06239.9246.83-15.176-0.484250.5046.31239.1146.07-15.551-0.533260.5245.61138.3545.35-15.919-0.578270.5444.95437.6444.67-16.270-0.621280.5644.33636.9744.04-16.614-0.660290.5843.75436.3543.45-16.922-0.697300.6043.20435.7642.89-17.230-0.732310.6242.68335.2042.36-17.532-0.764

    Table4Continued

    No.p/GPaV/cm3·mol-1NIST[18]VLWVHLerror/%VLWVHL320.6442.18934.6741.85-17.822-0.794330.6641.71934.1741.38-18.095-0.822340.6841.27233.7040.92-18.347-0.848350.7040.84533.2440.49-18.619-0.872360.7240.43832.8140.08-18.863-0.896370.7440.04932.4039.68-19.099-0.920380.7639.67632.0139.30-19.322-0.941390.7839.31931.6338.94-19.555-0.963400.8038.97631.2738.59-19.771-0.983

    Note: error=100×(VEOS-VNIST)/VNIST.

    4 結(jié) 論

    本研究提出了一種基于LJ勢能函數(shù)的對比態(tài)維里型狀態(tài)方程VHL用于描述爆轟環(huán)境下CO2的高溫高壓熱力學(xué)狀態(tài)。采用VHL狀態(tài)方程計算得CO2體積平均絕對偏差為0.971%,最大偏差為4.04%,采用VLW狀態(tài)方程計算所得平均絕對偏差20.2%,最大偏差87.149%。因此,在計算CO2高溫、中高壓熱力學(xué)狀態(tài)時,VHL狀態(tài)方程的計算準(zhǔn)確性得到了大幅度提高。VHL狀態(tài)方程計算所得體積偏差與壓力和溫度參量無明顯的相關(guān)性; VLW狀態(tài)方程計算所得體積偏差則與溫度具有明顯的相關(guān)性,隨著溫度的升高,計算體積偏差逐漸減小。

    參考文獻:

    [1] Mader C L. Numerical modeling of explosives and propellants[M]. CRC press, 2007: 377-380.

    [2] 吳雄, 龍新平, 何碧, 等. VLW 爆轟產(chǎn)物狀態(tài)方程[J]. 中國科學(xué): B 輯, 2008, 38(12): 1129-1132.

    WU Xiong,LONG Xin-ping, HE Bi, et al. VLW equation of state of detonation products[J].ScienceinChinaSeriesB:Chemistry, 2008, 38(12): 1129-1132.

    [3] Fickett W, Wood W W. Tables of the Lennard-Jones and Devonshire equation of state at high temperatures and densities[J].TheJournalofChemicalPhysics, 1952, 20(10): 1624-1626.

    [4] Fried L E, Howard W M. An accurate equation of state for the exponential-6 fluid applied to dense supercritical nitrogen[J].TheJournalofChemicalPhysics, 1998, 109(17): 7338-7348.

    [5] 吳雄, 龍新平, 何碧, 等. VLW狀態(tài)方程的回顧與展望[J]. 高壓物理學(xué)報, 1999, 13(1): 55-58.

    WU Xiong, LONG Xin-ping, HE Bi, et al. Review and look forward to the progress of VLW equation of state[J].ChineseJournalofHighPressurePhysics, 1999, 13(1): 55-58.

    [6] 龍新平, 何碧, 蔣曉華, 等. 論VLW狀態(tài)方程[J]. 高壓物理學(xué)報, 2003, 17(4): 247-254.

    LONG Xin-ping, HE Bi, JIANG Xiao-hua, et al. Discussions on the VLW equation of state[J].ChineseJournalofHighPressurePhysics, 2003, 17(4): 247-254.

    [7] 韓勇, 龍新平, 蔣治海, 等. 用VLW狀態(tài)方程計算水的沖擊Hugoniot曲線[J]. 爆炸與沖擊, 2010, 30(1): 17-20.

    HAN Yong, LONG Xin-ping, JIANG Zhi-hai, et al. A theoretical calculation for the hugoniot of water using the VLW equation of state[J].ExplosionandShockWaves, 2010, 30(1): 17-20.

    [8] Barker J A, Leonard P J, Pompe A. Fifth virial coefficients[J].TheJournalofChemicalPhysics, 1966, 44(11): 4206-4211.

    [9] 韓勇, 龍新平, 郭向利. 一種簡化維里型狀態(tài)方程預(yù)測高溫甲烷pVT關(guān)系[J]. 物理學(xué)報, 2014, 63(15): 150505.

    HAN Yong, LONG Xin-ping, GUO Xiang-li. Prediction of methane pVT relations at high temperatures by a simplified virial equation of state[J].ActaPhysSin, 2014, 63(15): 150505.

    [10] 韓勇, 龍新平, 黃毅民, 等. L-J, Exp-6 兩種形式勢能函數(shù)對計算無量綱第二維里系數(shù)的影響[J]. 含能材料, 2009, 17(5): 574-577.

    HAN Yong, LONG Xin-ping, HUANG Yi-min, et al. Effect of L-J or Exp-6 potential function on calculation of reduced second Viral coefficient[J].ChineseJournalofEnergeticMaterials(HannengCailliao), 2009, 17(5): 574-577.

    [11] Duan Z, M?ller N, Weare J H. Molecular dynamics simulation of PVT properties of geological fluids and a general equation of state of nonpolar and weakly polar gases up to 2000 K and 20,000 bar[J].GeochimicaetCosmochimicaActa, 1992, 56(10): 3839-3845.

    [12] Belonoshko A, Saxena S K. A molecular dynamics study of the pressure-volume-temperature properties of supercritical fluids: II. CO2, CH4, CO, O2and H2[J].GeochimicaetCosmochimicaActa, 1991, 55(11): 3191-3208.

    [13] Schott G L. Shock-compressed carbon dioxide: Liquid measurements and comparisons with selected models[J].HighPressureResearch, 1991, 6(3): 187-200.

    [14] Nellis W J, Mitchell A C, Ree F H, et al. Equation of state of shock-compressed liquds: Carbon dioxide and air[J].JChemPhys, 1991, 95: 5268-5272.

    [15] 劉福生, 陳先猛, 陳攀森, 等. 液態(tài)CO2高溫高密度狀態(tài)方程研究[J]. 高壓物理學(xué)報, 1998, 12(1): 28-33.

    LIU Fu-sheng, CHEN Xian-meng, CHEN Pan-sen, et al. Equation of state of liquod CO2at highe temperatures and high densities[J].ChineseJournalofHighPressurePhysics, 1998, 12(1): 28-33.

    [16] 徐士良. 常用算法程序集[M]. 北京: 清華大學(xué)出版社, 1995: 414.

    [17] Ben-Amotz D, Herschbach D R. Estimation of effective diameters for molecular fluids[J].JournalofPhysicalChemistry,1990, 94(3): 1038-1047.

    [18] Span R, Wagner W. A new equation of state for carbon dioxide covering the fluid region from the triple-point temperature to 1100 K at pressures up to 800 MPa[J].JournalofPhysicalandChemicalReferenceData, 1996, 25(6): 1509-1596.

    猜你喜歡
    狀態(tài)方程熱力學(xué)偏差
    LKP狀態(tài)方程在天然氣熱物性參數(shù)計算的應(yīng)用
    煤氣與熱力(2021年6期)2021-07-28 07:21:30
    如何走出文章立意偏差的誤區(qū)
    兩矩形上的全偏差
    基于隨機與區(qū)間分析的狀態(tài)方程不確定性比較
    Fe-C-Mn-Si-Cr的馬氏體開始轉(zhuǎn)變點的熱力學(xué)計算
    上海金屬(2016年1期)2016-11-23 05:17:24
    活塞的靜力學(xué)與熱力學(xué)仿真分析
    電子制作(2016年19期)2016-08-24 07:49:54
    關(guān)于均數(shù)與偏差
    用狀態(tài)方程模擬氨基酸水溶液的熱力學(xué)性質(zhì)
    用狀態(tài)方程模擬氨基酸水溶液的熱力學(xué)性質(zhì)
    一類非奇異黑洞的熱力學(xué)穩(wěn)定性
    日本免费a在线| 国产激情偷乱视频一区二区| 看十八女毛片水多多多| 最好的美女福利视频网| 深夜精品福利| 久久午夜亚洲精品久久| 精品久久久久久久久久久久久| 久久久欧美国产精品| 插逼视频在线观看| 久久久久久伊人网av| 亚洲欧美日韩卡通动漫| 亚洲av电影不卡..在线观看| 欧美绝顶高潮抽搐喷水| 亚洲成人中文字幕在线播放| 日日干狠狠操夜夜爽| 麻豆乱淫一区二区| 久久久久久久久中文| 不卡一级毛片| 精品午夜福利在线看| 卡戴珊不雅视频在线播放| 亚洲一区高清亚洲精品| 国产在线男女| 午夜福利在线在线| 欧美性猛交黑人性爽| 亚洲av免费高清在线观看| 亚洲图色成人| 人人妻人人澡人人爽人人夜夜 | 春色校园在线视频观看| 免费看日本二区| 欧美成人免费av一区二区三区| 可以在线观看毛片的网站| 国产激情偷乱视频一区二区| 自拍偷自拍亚洲精品老妇| 色综合站精品国产| 国产视频内射| 成人高潮视频无遮挡免费网站| 国产精品福利在线免费观看| 亚洲最大成人av| 熟女人妻精品中文字幕| 免费av观看视频| 国产爱豆传媒在线观看| 国产一区二区三区在线臀色熟女| 男人狂女人下面高潮的视频| 老熟妇乱子伦视频在线观看| 99热只有精品国产| 麻豆国产av国片精品| 亚洲人成网站高清观看| 久久精品91蜜桃| 高清日韩中文字幕在线| 精品欧美国产一区二区三| 丝袜喷水一区| 国产高清视频在线播放一区| 淫秽高清视频在线观看| 久久精品夜色国产| 有码 亚洲区| 国产亚洲av嫩草精品影院| 丝袜喷水一区| 国产一级毛片七仙女欲春2| 久久精品国产清高在天天线| 禁无遮挡网站| 可以在线观看毛片的网站| 老熟妇仑乱视频hdxx| 久久精品影院6| 最后的刺客免费高清国语| 一卡2卡三卡四卡精品乱码亚洲| 天天一区二区日本电影三级| 久久99热这里只有精品18| 夜夜爽天天搞| 一卡2卡三卡四卡精品乱码亚洲| 99久久九九国产精品国产免费| 婷婷色综合大香蕉| 最近中文字幕高清免费大全6| 免费看av在线观看网站| 不卡视频在线观看欧美| 日本 av在线| 夜夜爽天天搞| 中文字幕久久专区| 亚洲精品在线观看二区| 偷拍熟女少妇极品色| 欧美中文日本在线观看视频| 国产黄片美女视频| 亚洲精品色激情综合| 99久久精品国产国产毛片| av在线老鸭窝| 国产极品精品免费视频能看的| 久久久久久久久大av| 一边摸一边抽搐一进一小说| 天堂网av新在线| 国产精品一区二区性色av| 欧美精品国产亚洲| 欧美色欧美亚洲另类二区| 久久亚洲国产成人精品v| 日产精品乱码卡一卡2卡三| 中文字幕熟女人妻在线| 日本黄色片子视频| 国产精品1区2区在线观看.| 色综合亚洲欧美另类图片| 久久天躁狠狠躁夜夜2o2o| 国产 一区 欧美 日韩| 欧美三级亚洲精品| 国产精品爽爽va在线观看网站| 三级毛片av免费| 日韩欧美在线乱码| 久久这里只有精品中国| 欧美zozozo另类| 少妇熟女欧美另类| 午夜久久久久精精品| 国产精品一区二区免费欧美| 亚洲丝袜综合中文字幕| 国产精品亚洲美女久久久| 全区人妻精品视频| 少妇人妻一区二区三区视频| 国产成年人精品一区二区| 亚洲七黄色美女视频| 国产精品一二三区在线看| 国产乱人偷精品视频| 丝袜喷水一区| 婷婷亚洲欧美| 久久久午夜欧美精品| 国产精品国产三级国产av玫瑰| 岛国在线免费视频观看| 美女 人体艺术 gogo| 国产 一区 欧美 日韩| 国产探花极品一区二区| 热99在线观看视频| 乱系列少妇在线播放| 在线国产一区二区在线| 精华霜和精华液先用哪个| videossex国产| 99热这里只有是精品50| 大型黄色视频在线免费观看| 亚洲最大成人av| 中出人妻视频一区二区| 少妇高潮的动态图| 国产精品国产三级国产av玫瑰| 国产人妻一区二区三区在| 日韩欧美在线乱码| 一a级毛片在线观看| videossex国产| 青春草视频在线免费观看| 亚洲av中文av极速乱| 观看美女的网站| 一卡2卡三卡四卡精品乱码亚洲| 熟女电影av网| a级毛片a级免费在线| 成人鲁丝片一二三区免费| 国语自产精品视频在线第100页| 免费电影在线观看免费观看| a级毛色黄片| 白带黄色成豆腐渣| 亚洲熟妇熟女久久| 桃色一区二区三区在线观看| 欧美区成人在线视频| 久久人人爽人人爽人人片va| 亚洲欧美精品综合久久99| 亚洲精品一卡2卡三卡4卡5卡| 免费不卡的大黄色大毛片视频在线观看 | 午夜久久久久精精品| 99久久久亚洲精品蜜臀av| or卡值多少钱| 老司机福利观看| 99视频精品全部免费 在线| 99在线人妻在线中文字幕| 少妇熟女欧美另类| 精品熟女少妇av免费看| 少妇裸体淫交视频免费看高清| 久久精品国产亚洲av天美| 1000部很黄的大片| 欧美成人a在线观看| 黄色欧美视频在线观看| 精品久久久久久久久亚洲| 日韩高清综合在线| 成人永久免费在线观看视频| 深夜精品福利| 亚洲人成网站在线播| 欧美日韩综合久久久久久| 亚洲内射少妇av| 极品教师在线视频| 国产精品综合久久久久久久免费| 免费人成视频x8x8入口观看| 欧美高清成人免费视频www| 黄片wwwwww| 美女cb高潮喷水在线观看| 欧美高清成人免费视频www| 国产高清激情床上av| 午夜精品在线福利| 日韩精品中文字幕看吧| 哪里可以看免费的av片| 内地一区二区视频在线| 午夜福利在线在线| 老司机影院成人| 搡老熟女国产l中国老女人| 搡女人真爽免费视频火全软件 | www日本黄色视频网| 亚洲最大成人手机在线| 久久精品国产亚洲av涩爱 | 午夜福利视频1000在线观看| 欧美日韩国产亚洲二区| 波多野结衣高清作品| 国产亚洲欧美98| 全区人妻精品视频| 女人十人毛片免费观看3o分钟| 国产在线精品亚洲第一网站| 看非洲黑人一级黄片| 亚洲国产精品成人久久小说 | 日韩精品有码人妻一区| 非洲黑人性xxxx精品又粗又长| 婷婷亚洲欧美| 深爱激情五月婷婷| 大又大粗又爽又黄少妇毛片口| 亚洲自拍偷在线| 俺也久久电影网| 在线国产一区二区在线| 国产一区亚洲一区在线观看| 精品午夜福利视频在线观看一区| 日韩av不卡免费在线播放| 亚洲aⅴ乱码一区二区在线播放| 欧美潮喷喷水| 成人国产麻豆网| 亚洲欧美中文字幕日韩二区| 精品福利观看| 黄色视频,在线免费观看| 亚洲国产精品国产精品| www日本黄色视频网| 在线免费十八禁| 免费无遮挡裸体视频| 精品久久久久久久久av| av黄色大香蕉| 美女免费视频网站| 最近2019中文字幕mv第一页| 国产亚洲精品久久久久久毛片| 狂野欧美激情性xxxx在线观看| 别揉我奶头 嗯啊视频| 日韩av不卡免费在线播放| 久久午夜福利片| 午夜日韩欧美国产| 国产蜜桃级精品一区二区三区| 日韩 亚洲 欧美在线| 网址你懂的国产日韩在线| 成年女人永久免费观看视频| 久久人妻av系列| 亚洲美女黄片视频| 亚洲成人久久爱视频| 婷婷精品国产亚洲av| 欧美性猛交黑人性爽| 亚洲天堂国产精品一区在线| 久久精品夜夜夜夜夜久久蜜豆| 最近在线观看免费完整版| 少妇的逼水好多| 日日啪夜夜撸| 中文字幕熟女人妻在线| 久久久久国内视频| 欧洲精品卡2卡3卡4卡5卡区| 一个人看的www免费观看视频| 亚洲av电影不卡..在线观看| 国产日本99.免费观看| 精品久久久久久久久久免费视频| 亚洲第一区二区三区不卡| 久久中文看片网| 美女cb高潮喷水在线观看| 99热精品在线国产| 在线观看av片永久免费下载| 在线免费观看的www视频| 性色avwww在线观看| 啦啦啦观看免费观看视频高清| 国产爱豆传媒在线观看| 色综合站精品国产| 国产精品伦人一区二区| 国产毛片a区久久久久| 久久久久久久久大av| 日韩亚洲欧美综合| 国产av麻豆久久久久久久| 国产 一区精品| 精品人妻熟女av久视频| 悠悠久久av| 男女视频在线观看网站免费| 精品少妇黑人巨大在线播放 | 此物有八面人人有两片| 99久国产av精品国产电影| 亚洲国产精品合色在线| 免费大片18禁| 日韩在线高清观看一区二区三区| 国产黄a三级三级三级人| 亚洲图色成人| 搡女人真爽免费视频火全软件 | 中文字幕av成人在线电影| 春色校园在线视频观看| 国产熟女欧美一区二区| 午夜老司机福利剧场| 国产探花极品一区二区| 国产黄色小视频在线观看| av专区在线播放| 久久精品国产清高在天天线| 久久久久精品国产欧美久久久| 午夜福利18| 三级毛片av免费| 麻豆久久精品国产亚洲av| 97热精品久久久久久| 国产精品国产高清国产av| 激情 狠狠 欧美| 国产色婷婷99| 深夜精品福利| 淫妇啪啪啪对白视频| 亚洲欧美成人综合另类久久久 | 熟女人妻精品中文字幕| 午夜福利高清视频| 国产片特级美女逼逼视频| 国产高清不卡午夜福利| av卡一久久| 欧美+日韩+精品| 久久久久久国产a免费观看| 十八禁国产超污无遮挡网站| 日韩三级伦理在线观看| 老女人水多毛片| 草草在线视频免费看| 久久久国产成人精品二区| 天天躁夜夜躁狠狠久久av| 久久九九热精品免费| 波野结衣二区三区在线| 夜夜爽天天搞| 三级男女做爰猛烈吃奶摸视频| 亚洲七黄色美女视频| av在线亚洲专区| 激情 狠狠 欧美| 国产在线男女| 国产精品三级大全| 国产欧美日韩精品一区二区| 97在线视频观看| 99热精品在线国产| 久久久久久国产a免费观看| 免费观看精品视频网站| 国产黄a三级三级三级人| 天堂网av新在线| 免费不卡的大黄色大毛片视频在线观看 | 三级国产精品欧美在线观看| 99久久精品国产国产毛片| 欧美区成人在线视频| 亚洲五月天丁香| 成年版毛片免费区| av卡一久久| 免费无遮挡裸体视频| 亚洲av熟女| 久久久久国产网址| 亚洲欧美日韩卡通动漫| 国产白丝娇喘喷水9色精品| 精品免费久久久久久久清纯| 久久久久国产网址| 毛片一级片免费看久久久久| 小说图片视频综合网站| 男女做爰动态图高潮gif福利片| 亚洲av电影不卡..在线观看| 国产高清激情床上av| 99热这里只有精品一区| 久久久成人免费电影| 国产91av在线免费观看| 国产蜜桃级精品一区二区三区| 精品久久久久久成人av| 国产高清激情床上av| 亚洲久久久久久中文字幕| 久久精品国产自在天天线| 精品日产1卡2卡| 深夜a级毛片| 夜夜夜夜夜久久久久| 国产精品久久久久久亚洲av鲁大| 少妇熟女欧美另类| 国产真实伦视频高清在线观看| a级毛色黄片| 国产男靠女视频免费网站| 狂野欧美白嫩少妇大欣赏| 又粗又爽又猛毛片免费看| 成人av在线播放网站| 久久久久久九九精品二区国产| 成人av在线播放网站| 国产爱豆传媒在线观看| 麻豆乱淫一区二区| 91av网一区二区| 极品教师在线视频| 久久精品国产亚洲av香蕉五月| 欧美极品一区二区三区四区| 亚洲欧美成人综合另类久久久 | 蜜桃亚洲精品一区二区三区| 一级a爱片免费观看的视频| 女人被狂操c到高潮| 国产 一区精品| 女的被弄到高潮叫床怎么办| 久久婷婷人人爽人人干人人爱| 亚洲久久久久久中文字幕| 婷婷六月久久综合丁香| 国产在视频线在精品| 久久久久久久久大av| 久久久a久久爽久久v久久| 99久久中文字幕三级久久日本| 久久久a久久爽久久v久久| 亚洲18禁久久av| 精品久久久噜噜| 国产成人福利小说| 国产精品不卡视频一区二区| 天美传媒精品一区二区| 人人妻人人澡人人爽人人夜夜 | 老熟妇乱子伦视频在线观看| 俺也久久电影网| 永久网站在线| 在线免费观看不下载黄p国产| 日韩精品青青久久久久久| 精品无人区乱码1区二区| 欧美日本亚洲视频在线播放| 国产v大片淫在线免费观看| 最新中文字幕久久久久| 国产乱人偷精品视频| 国产精品国产三级国产av玫瑰| 黄色欧美视频在线观看| or卡值多少钱| 黄色配什么色好看| 中国美女看黄片| 欧美区成人在线视频| 男人和女人高潮做爰伦理| 精品人妻偷拍中文字幕| 桃色一区二区三区在线观看| 亚洲精品亚洲一区二区| 男人的好看免费观看在线视频| 小蜜桃在线观看免费完整版高清| 国产精品电影一区二区三区| 久久欧美精品欧美久久欧美| 久久天躁狠狠躁夜夜2o2o| 国产色婷婷99| 九九久久精品国产亚洲av麻豆| 久久人人爽人人片av| 身体一侧抽搐| 一个人观看的视频www高清免费观看| 中国美白少妇内射xxxbb| 国产黄a三级三级三级人| 国产色爽女视频免费观看| 国内揄拍国产精品人妻在线| 男人的好看免费观看在线视频| 亚洲经典国产精华液单| 超碰av人人做人人爽久久| 亚洲av中文字字幕乱码综合| 我的女老师完整版在线观看| 啦啦啦啦在线视频资源| 91午夜精品亚洲一区二区三区| 午夜福利视频1000在线观看| 人人妻人人澡人人爽人人夜夜 | 国产精品免费一区二区三区在线| 亚洲人成网站在线播| 蜜臀久久99精品久久宅男| 极品教师在线视频| 五月伊人婷婷丁香| 丝袜喷水一区| 国产精品爽爽va在线观看网站| 亚洲无线观看免费| 卡戴珊不雅视频在线播放| 欧美日韩一区二区视频在线观看视频在线 | 如何舔出高潮| 成人毛片a级毛片在线播放| 精品人妻熟女av久视频| 色在线成人网| 国产在视频线在精品| 国产黄a三级三级三级人| 亚洲激情五月婷婷啪啪| 精品人妻熟女av久视频| av中文乱码字幕在线| 如何舔出高潮| 日韩三级伦理在线观看| 中国美白少妇内射xxxbb| 色播亚洲综合网| 日本-黄色视频高清免费观看| 午夜爱爱视频在线播放| 日本色播在线视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品粉嫩美女一区| 亚洲,欧美,日韩| 久久天躁狠狠躁夜夜2o2o| 人妻夜夜爽99麻豆av| 97超级碰碰碰精品色视频在线观看| 美女免费视频网站| 亚洲av成人精品一区久久| 亚洲最大成人av| 亚洲色图av天堂| 好男人在线观看高清免费视频| 亚洲精品一区av在线观看| 亚洲av电影不卡..在线观看| 波多野结衣高清作品| 久久亚洲国产成人精品v| 夜夜爽天天搞| 午夜福利高清视频| 成熟少妇高潮喷水视频| 欧美精品国产亚洲| 一个人看视频在线观看www免费| 亚洲精品久久国产高清桃花| 国产高清视频在线播放一区| 亚洲欧美精品综合久久99| 99热只有精品国产| 欧美日韩国产亚洲二区| 亚洲va在线va天堂va国产| 嫩草影视91久久| 99热这里只有精品一区| 少妇的逼好多水| 国产成人一区二区在线| 久久久久九九精品影院| 国产精品久久久久久亚洲av鲁大| 亚洲图色成人| 丰满乱子伦码专区| 99久久九九国产精品国产免费| 老司机影院成人| 国产精品无大码| 久久人妻av系列| 亚洲性夜色夜夜综合| av天堂中文字幕网| 一级黄片播放器| avwww免费| 欧美一级a爱片免费观看看| 成人漫画全彩无遮挡| 一个人看的www免费观看视频| 简卡轻食公司| 久久精品夜夜夜夜夜久久蜜豆| 啦啦啦韩国在线观看视频| 亚洲七黄色美女视频| 不卡一级毛片| 如何舔出高潮| 亚洲熟妇熟女久久| 国产精品一区二区三区四区免费观看 | 亚洲最大成人中文| 最近在线观看免费完整版| 18禁黄网站禁片免费观看直播| 观看免费一级毛片| 搡老岳熟女国产| 欧美一级a爱片免费观看看| 校园人妻丝袜中文字幕| 亚洲一区二区三区色噜噜| 久久99热6这里只有精品| 我的女老师完整版在线观看| 亚洲无线观看免费| 99久久精品国产国产毛片| 村上凉子中文字幕在线| 乱码一卡2卡4卡精品| 俄罗斯特黄特色一大片| 免费观看在线日韩| 成人鲁丝片一二三区免费| 亚州av有码| 久久精品国产亚洲av涩爱 | 中文字幕久久专区| 亚洲av熟女| 一进一出抽搐gif免费好疼| 中文字幕免费在线视频6| 国产一区二区激情短视频| 看十八女毛片水多多多| 中文字幕精品亚洲无线码一区| 成人国产麻豆网| 中文字幕av在线有码专区| 欧美xxxx黑人xx丫x性爽| 日韩精品青青久久久久久| 又爽又黄a免费视频| av在线天堂中文字幕| 亚洲电影在线观看av| 老司机午夜福利在线观看视频| 天堂√8在线中文| 亚洲最大成人中文| 可以在线观看的亚洲视频| 男女下面进入的视频免费午夜| 精品国内亚洲2022精品成人| 一个人观看的视频www高清免费观看| 男女啪啪激烈高潮av片| 日本三级黄在线观看| 最近视频中文字幕2019在线8| 小说图片视频综合网站| 成人鲁丝片一二三区免费| 天堂动漫精品| 欧美不卡视频在线免费观看| 午夜激情欧美在线| 亚洲一级一片aⅴ在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 给我免费播放毛片高清在线观看| 亚洲精品成人久久久久久| 免费一级毛片在线播放高清视频| 少妇猛男粗大的猛烈进出视频 | 天美传媒精品一区二区| 99久久久亚洲精品蜜臀av| 高清毛片免费看| av在线天堂中文字幕| 国内精品美女久久久久久| 亚洲av电影不卡..在线观看| 综合色av麻豆| 亚洲人成网站在线播放欧美日韩| 男人狂女人下面高潮的视频| 精品久久久久久久久av| 一级黄色大片毛片| aaaaa片日本免费| 国产精品一区二区三区四区久久| 白带黄色成豆腐渣| 欧美色欧美亚洲另类二区| 天堂动漫精品| 婷婷精品国产亚洲av| 最新在线观看一区二区三区| 99在线人妻在线中文字幕| av天堂中文字幕网| 欧美中文日本在线观看视频| 麻豆国产97在线/欧美| 日本成人三级电影网站| 久久草成人影院| 久久精品综合一区二区三区| 在线天堂最新版资源| 亚洲av不卡在线观看| 噜噜噜噜噜久久久久久91| 你懂的网址亚洲精品在线观看 | aaaaa片日本免费| 波多野结衣高清作品| 人人妻人人澡人人爽人人夜夜 | 神马国产精品三级电影在线观看| 国产综合懂色| а√天堂www在线а√下载| 午夜福利高清视频| 国产又黄又爽又无遮挡在线| 两个人视频免费观看高清| 免费在线观看影片大全网站| 久久精品国产亚洲网站|