專利名稱：蒙特卡洛模擬預測微管式反應器內(nèi)自由基共聚的方法
技術(shù)領(lǐng)域：
本發(fā)明涉及高分子聚合技術(shù)領(lǐng)域，具體為一種蒙特卡洛模擬預測微管式反應器內(nèi)自由基共聚的方法。
背景技術(shù)：
圓底燒瓶是合成化學中用來合成毫克至克級別的經(jīng)典的化學反應容器。但其散熱效果不佳和存在“過熱點”，所以不能很好的控制劇烈放熱的反應。相對于圓底燒瓶，微管式反應器的小尺寸提供了獨特的性能，這些性能是傳統(tǒng)的反應容器不能做到的。比如，因為較短的擴散路徑，所以他們可以做到使溶液非常快的混合。而這就提高了傳質(zhì)效率、提供了快速動力學反應的理想條件。另外，和傳統(tǒng)反應器lOO~lOOOm2.πΓ3相比，這種設(shè)備有著非常大的比表面積，大約10000 50000 m2.πΓ3，熱交換性能更優(yōu)，這樣可以使得快速升降溫得以實現(xiàn)。因為高效的熱量控制和顯著的混合效率，可以被抑制副產(chǎn)物和過熱點的形成，從而提高對劇烈放熱反應的控制。另外，通過升高微反應器內(nèi)的溫度和壓力可以改善慢反應。而且，使用小尺寸的連續(xù)流體反應裝置，保證了那些需要處理少量有潛在放熱、毒性和易爆等不太適合常規(guī)反應器的反應的操作者的安全性。相對于傳統(tǒng)的過程，微管式裝置可以在反應過程中便利的控制諸如溫度、保留時間、反應物當量比等反應參數(shù)。Monte Carlo (蒙特卡洛)模擬方法是一種有效的隨機統(tǒng)計方法，非常適合研究包含多個基元反應的共聚體系，可以得到傳統(tǒng)的解析處理得不到的聚合物反應動力學、分子量分布，共聚組成和序列分布等信息。將微管式反應器應用于化學生產(chǎn)中，流體在微管中邊流動邊發(fā)生反應生成共聚物。共聚物的序列分布受原料的投入比，流體流動的的速率等參數(shù)影響。使用Monte Carlo模擬方法可以預測出最合理的工藝條件。

發(fā)明內(nèi)容
本發(fā)明的目的在于提供一種可以預測出微管式反應器內(nèi)最合理的工藝條件的蒙特卡洛模擬方法。該方法簡便可行，更加接近真實過程。本發(fā)明提出的蒙特卡洛模擬預測微管式反應器內(nèi)傳統(tǒng)自由基共聚的方法，是從共聚反應的基元反應出發(fā)，模擬原料在微管式反應器中邊流動邊發(fā)生化學反應，最后得到的產(chǎn)物共聚組成和序列分布均勻，得到最合理的工藝條件。這種算法的優(yōu)點是不需要做各種假設(shè)，只要知道基元反應，讓隨機發(fā)生，并且根據(jù)流速使得原料流動，并加入物種的擴散，因而更加接近真實過程。而且可以將預設(shè)的工藝條件加入模擬中，通過觀察不同工藝條件下得到的共聚產(chǎn)物的序列分布來判斷工藝條件的合理性，避免大量的實驗消耗時間和成本。本發(fā)明提出的蒙特卡洛模擬預測微管式反應器內(nèi)傳統(tǒng)自由基共聚的方法，具體步驟如下:
(I)確定聚合體系中所有的基元反應和擴散反應，從聚合物手冊中查出速率常數(shù)的數(shù)值，確定各種基元反應和擴散反應的速率常數(shù)，確定各種基元反應的速率常數(shù)K.和擴散速率常數(shù)Kdiff—ma。，再將所述的所有的基元反應的宏觀速率常數(shù)K.轉(zhuǎn)化為蒙特卡洛模擬的微觀反應速率常數(shù)1^。 、擴散反應的宏觀擴散速率常數(shù)Kdiffma。轉(zhuǎn)化為蒙特卡洛模擬的微觀速率常數(shù)K—。，其中:對于一級反應，Kfflac=Kfflicro;對于不同物種之間的二級反應，Kfflicro=Kfflac/(VNa);對于同一物種之間的二級反應，Kmicro=2Kmac/(VNa), Na為阿伏伽德羅常數(shù)，大小為
6.02X IO23, V為反應體系的體積，Kmicro=Kmac*n2/L2, η為子體系的個數(shù)，L為體系的總長度。對三元共聚反應,其基元反應包括自由基引發(fā)反應、鏈增長反應、終止反應。鏈增長部分包括Ml自由基分別與單體Ml與Μ2、Μ3的加成反應，M2自由基分別與單體Ml與M2、M3的加成反應,一共九組基元反應反應；終止反應包括Ml自由基自終止反應；Μ2自由基自終止反應；Μ3自由基自終止反應。Ml自由基與M2自由基終止反應；M1自由基與M3自由基終止反應；M2自由基與M3自由基終止反應。對于三元共聚反應，三元自由基共聚體系的速率常數(shù)包括引發(fā)速率常數(shù)Ki，Ml自由基與單體Ml加成增長速率常數(shù)Kl I:M1自由基與單體M2增長速率常數(shù)K12 ;M1自由基與單體M3增長速率常數(shù)K13 ；M2自由基與單體Ml增長速率常數(shù)K21 ；M2自由基與單體M2增長速率常數(shù)K22 ；M2自由基與單體Ml增長速率常數(shù)K23 ；M3自由基與單體Ml增長速率常數(shù)K31 ；M3自由基與單體M2增長速率常數(shù)K32 ；M3自由基與單體M3增長速率常數(shù)K33 ；M1自由基與Ml自由基終止速率常數(shù)Ktl I ；M1自由基與M2自由基終止速率常數(shù)Ktl2 ;M1自由基與M3自由基終止速率常數(shù)Ktl3 ；M2自由基與M2自由基終止速率常數(shù)Kt22 ；M2自由基與M3自由基終止速率常數(shù)Kt23 ；M3自由基與M3自由基終止速率常數(shù)Kt33。對于擴散反應，大分子的擴散速率常數(shù)和小分子的擴散速率常數(shù)是不一樣的，不同物種在不同溶質(zhì)中的擴散速率常數(shù)也是不一樣的。(2)預設(shè)微管式反應器工藝條件:包括微管式反應器的容積、加入管中反應的不同原料的進料量和管內(nèi)流體的流速。

(3)使用蒙特卡洛算法，基于各種基元反應和分批加料的工藝條件建立自由基共聚反應的算法，其動力學的蒙特卡洛模擬步驟為:
(a)確定基元反應和分批加料的工藝條件，輸入所有基元反應速率常數(shù)和所有種類化學物種的分子個數(shù)，并將宏觀反應速率常數(shù)變換為微觀反應速率常數(shù)，將時間t設(shè)定為O ;
(b)將微管分為η個單元，計算每個單元的容積、并計算每個單元內(nèi)的所有種類的化學物種的分子個數(shù)，和各種基元反應、擴散反應的速率以及速率加和，計算各反應概率和各單元反應的總概率；
(c)產(chǎn)生兩個單位區(qū)間內(nèi)均勻分布的隨機數(shù)，判斷下一步反應將位于哪一單元；
(d)產(chǎn)生兩個單位區(qū)間內(nèi)均勻分布的隨機數(shù)，判斷下一步將發(fā)生何種反應，相應的調(diào)整反應所涉及的化學物種的分子個數(shù)，計算時間間隔，累加到時間t上去；
(e)根據(jù)確定的微管式反應器的工藝條件，在設(shè)定的時間點上加入一定比例的原料改變相應的化學物種的分子個數(shù)，并根據(jù)流速計算因為流動而改變的化學物種的分子個數(shù)；多次重復步驟(b) (c) (d) (e)，即得到各化學物種隨時間的演化；
(f)將物種的變化分別記錄在不同數(shù)組中，并進行統(tǒng)計；
共聚聚合反應的模擬中，所述化學物種為引發(fā)劑、單體、進入共聚物的單體和共聚物；其中各個鏈自由基都用下述四個參數(shù):標號、鏈長、自由基種類和鏈段分布來標記，并分別儲存在兩個關(guān)聯(lián)的數(shù)據(jù)中，便于統(tǒng)計。其應用程序的運算過程可以做這樣形象化的描述:首先依據(jù)步驟(I)中確定的基元反應的速率參數(shù)來計算研究體系中各基元反應的概率，然后由隨機數(shù)來決定下一步將在哪個單元內(nèi)發(fā)生何種基元反應。在某一反應發(fā)生時，統(tǒng)計所涉及到的所有相應化學物種的變化，其中包括各種單體的個數(shù)，鏈自由基的長度，共聚物中不同單體的鏈段的長度，共聚物的鏈長，并將各種單體分子個數(shù)，引發(fā)劑分子個數(shù)，各種鏈自由基個數(shù)，各種鏈自由基個數(shù)的長度，共聚物中不同單體的鏈段的長度，聚合物分子量分別儲存到相應的數(shù)組里。同時每根鏈自由基有四個參數(shù)用來標記標號，鏈長，自由基種類，鏈段分布。例如Ml長鏈自由基與單體M2發(fā)生增長反應時，體系中單體M2的數(shù)量需相應地減去一個單體，Ml長鏈自由基數(shù)量也減少一個，并在表征Ml長鏈自由基的數(shù)組里，隨機的挑選一條長度為L的分子鏈，將其長度賦值為L+1，然后Ml長鏈自由基數(shù)量減少一個，同時將此長鏈自由基中記錄序列分布的數(shù)列數(shù)組中的Ml的鏈段長度S記錄到統(tǒng)計序列分布的數(shù)組中。產(chǎn)生一個長度為S的Ml單體的連段。將M2長鏈自由基數(shù)量相應增加一個，其長度為L+l，M2鏈段長度S賦值為I。例如Ml長鏈自由基與單體Ml發(fā)生增長反應時，體系中單體Ml的數(shù)量需相應地減去一個單體，并在表征Ml長鏈自由基的數(shù)組里，隨機的挑選一條長度為L的分子鏈，將其長度賦值為L+1,然后將長度為L的鏈的數(shù)目減去I，而將長度為L+1的鏈的數(shù)目上加上I,同時將此長鏈自由基中記錄序列分布的數(shù)組中的Ml鏈段長度S賦值S+1.例如Ml長鏈自由基與M2長鏈自由基發(fā)生終止反應時，需在表征Ml長鏈自由基的數(shù)組里，隨機的挑選一條長度為m的分子鏈，在表征M2長鏈自由基的數(shù)組里，隨機的挑選一個長度為η的分子鏈，如果是發(fā)生歧化終止反應，產(chǎn)生分子量分別是m和η兩根聚合物，并在表征聚合物的數(shù)組中記錄。如果是發(fā)生偶合終止反應，產(chǎn)生一根分子量是m+n聚合物，并在表征聚合物的數(shù)組中記錄。例如Ml長鏈自由基從第η單元擴散至η+1單元，則第η單元內(nèi)在表征Ml長鏈自由基的數(shù)組里，隨機的挑選一條長度為L的分子鏈，將其長度、序列長度賦值給第η+1單元內(nèi)一條新的鏈，η單元內(nèi)的Ml長鏈自由基的數(shù)目減去1，自由基總個數(shù)減去1，η+1單元內(nèi)的Ml長鏈自由基個數(shù)加上1，自由基總的個數(shù)加上I。每一步反應之后由隨機數(shù)來決定其和上一步反應之間的時間間隔并累加到時間t上，即可得到各參量隨反應時間變化的信息。計算流動可以做這樣形象化的描述:除了第一單元，每個單元流出的流體體積等于流入該單元的流體體積，而第一單元流出的流體體積等于加入該單元的原料體積之和。流出的流體組分組成和該單元中的組成比例相同。如第三單元流入第四單元。每隔一段時間，用流速乘以時間計算從第三、四單元流出的流體體積，該流體體積內(nèi)的組分比例與第三、四單元內(nèi)組分比例相同，例如第三單元中流出的nl個Ml單體，n2個鏈長為L的Ml長鏈自由基，第四單元中流出的ml個Ml單體，m2個鏈長為L的Ml長鏈自由基，第四單元中原有k2個Ml單體，k3個鏈長為L的Ml長鏈，則第四單元中在該時間點上的組成為nl-ml+K2個Ml單體，n2+K3_m2個鏈長為L的Ml長鏈自由基。如第一單元原有nl個Ml單體，n2個鏈長為L的Ml長鏈自由基，投入kl個Ml單體，第一單元中流出ml個Ml單體，m2個鏈長為L的Ml長鏈自由基，則第一單元中在該時間點上的組成為nl-ml+K個Ml單體，n2~m2個鏈長為L的Ml長鏈自由基。(4)對步驟(3)模擬結(jié)果進行統(tǒng)計，得到不同的轉(zhuǎn)化率階段的序列分布。在模擬程序運算過程中，已經(jīng)對體系中的每一物種，包括單體，引發(fā)劑，自由基，聚合物的變化進行精確的跟蹤記錄，并分別記錄在不同數(shù)組中，隨著反應的進行，程序會對各種參數(shù)進行統(tǒng)計，并記錄在數(shù)組中，我們可以得到所需的結(jié)果，具體結(jié)果如下:
(a)通過跟蹤統(tǒng)計所有單體的變化可以推導出微管式反應器的工藝條件對反應體系的動力學的影響；
(b)通過跟蹤統(tǒng)計已經(jīng)進入共聚物的單體Ml和M2的變化可以推導導出微管式反應器的工藝條件對共聚物組成的影響；
(C)通過跟蹤統(tǒng)計每個聚合物聯(lián)眾所有的鏈段長度可以推導導出微管式反應器的工藝條件對共聚物鏈段序列分布的影響。(5)根據(jù)步驟(3)、(4)得到序列分布結(jié)果考察預設(shè)的微管式反應器的工藝條件的合理性，重新設(shè)定工藝條件。直到得到滿意的序列分布結(jié)果和相對應的分批加料的工藝條件，將這個工藝條件應用于共聚反應的試驗中。在化學反應動力學中，化學反應的微觀本質(zhì)是分子之間的反應概率問題，是一個隨機過程。根據(jù)對隨機事件的抽樣原理，和某一時間間隔內(nèi)所發(fā)生的化學反應的種類均由單元區(qū)間內(nèi)均勻分布的隨機數(shù)來決定，若將化學物種的濃度定義為每單元體積中的摩爾數(shù)，在實際模擬時，相繼發(fā)生的兩次基元反應之間的時間間隔t，不是一個常數(shù)，而是一個隨機變量，它由單元區(qū)間內(nèi)均勻分布的隨機數(shù)來決定。本發(fā)明中，化學物種包括引發(fā)劑，單體，進入共聚物的單體，共聚物。本發(fā)明中，所述的共聚產(chǎn)物的鏈段序列分布，指不同單體在聚合物中的鏈段長度的分布趨勢。本發(fā)明的優(yōu)點在于:本方法簡單可行，可以預測出微管式反應器內(nèi)最合理的工藝條件。


圖1為實施例1的裝置示意圖。圖2為實施例1中轉(zhuǎn)化率隨時間的變化。圖3為實施例1中的聚合度與聚合度分布隨時間的變化關(guān)系。圖4為實施例1中的序列長度與序列長度分布隨時間的變化關(guān)系。圖5為實施例1的序列長度分布曲線。
具體實施例方式下面結(jié)合實施例對本發(fā)明作進一步詳細說明。實施例1
圖1為實施例1的裝置示意圖。物料從I處進，從N處出。物料在管中流動的長度等于微管總長。模擬54°C溫度下，丙烯腈(AN)、衣康酸(Itaconic acid)、丙烯酸甲酯(MA)在微管式反應器內(nèi)自由基三元共聚，溶劑為40%的硫氰酸鈉水溶液。具體實施步驟如下:
(I)確定聚合體系中所有的基元反應和擴散反應，
A.氧化還原反應產(chǎn)生兩個初級自由基；
B.初級自由基和丙烯腈、衣康酸、丙烯酸甲酯加成生成丙烯腈增長自由基、衣康酸增長自由基、丙烯酸甲酯增長自由基，共3個基元反應；
C.丙烯腈增長自由基、衣康酸增長自由基、丙烯酸甲酯增長自由基分別和丙烯腈單體、衣康酸單體、丙烯酸甲酯單體加成，共9個基元反應；
D.初級自由基、丙烯腈增長自由基、衣康酸增長自由基、丙烯酸甲酯增長自由基相互終止反應(此時只考慮雙基終止)，則共10個基元反應；
E.氧化劑、還原劑、丙烯腈單體、衣康酸單體、丙烯酸甲酯單體、初級自由基、丙烯腈增長自由基、衣康酸增長自由基、丙烯酸甲酯增長自由基、終止分子的擴散反應，共11個擴散反應。(2)從聚合物手冊中查出速率常數(shù)的數(shù)值，確定各種基元反應的速率常數(shù)Kma。，再將所述的所有的基元反應的宏觀速率常數(shù)K.轉(zhuǎn)化為蒙特卡洛模擬的微觀反應速率常數(shù)Kmicro;對于一級反應Kmac=Kmira。。對于不同物種之間的二級反應^^^洲仏對于同一物種之間的二級反應Kmiero=2Kmae/(VNa) ,Na為阿伏伽德羅常數(shù)，大小為6.02X 1023，V為反應體系的體積。對于擴散反應，Kfflicro=Kfflac*n2/L2, η為子體系的個數(shù)，L為體系的總長度。(3)預設(shè)微管式反應器的工藝條件:
Α.微管式反應器的容積為:4.71*10_5 L ；
B.進料量:溶劑為3.3*10_6 L/h，丙烯腈的進料量為6.719*10_5 L/h，衣康酸的進料量為5.9*10_5L/h，丙烯酸甲酯的進料量為2.24*10_6L/h，氧化劑的進料量為2.455*10_5L/h，還原劑的進料量為8.457*10_5L/h。(4)基于各種基元反應和預設(shè)微管式反應器工藝條件，使用蒙特卡洛算法，建立自由基共聚反應的算法，其動力學的蒙特卡洛模擬步驟為:
Ca)輸入步驟(I)中所有基元反應的微觀反應速率常數(shù)及所有種類化學物種的分子個數(shù)，將時間t設(shè)定為O;
(b)將微管分為30個單元，計算每個單元的容積、并計算每個單元內(nèi)的所有種類的化學物種的分子個數(shù)，和各種基元反應、擴散反應的速率以及速率加和，計算各反應概率和各單元反應的總概率；
(C)產(chǎn)生兩個單位區(qū)間內(nèi)均勻分布隨機數(shù)，判斷下一步反應位于哪一單元，
(d)產(chǎn)生兩個單位區(qū)間內(nèi)均勻分布的隨機數(shù)，判斷下一步將發(fā)生何種反應，相應的調(diào)整反應所涉及的化學物種的分子個數(shù)，計算時間間隔，累加到時間t上去；
(e)根據(jù)步驟(2)中的所述的預設(shè)的微管式反應器的工藝條件，在設(shè)定的時間點上加入一定比例的原料改變相應的化學物種的分子個數(shù)，并根據(jù)流速計算因為流動而改變的化學物種的分子個數(shù)；
(f)多次重復步驟(b)(c) (d) (e)，即可得到各化學物種隨時間的變化情況；將所述化學物種的變化分別記錄在不同數(shù)組中，并進行統(tǒng)計，得到模擬結(jié)果； 其中步驟(3)中所述化學物種包括引發(fā)劑、單體、進入共聚物的單體和共聚物；步驟(f)中所述統(tǒng)計所述化學物種的變化記錄在不同數(shù)組，即將各種單體分子的個數(shù)，引發(fā)劑分子的個數(shù)，各種鏈自由基的個數(shù)，各種鏈自由基的長度，共聚物中不同單體的鏈段的長度，聚合物的分子量分別儲存到相應的數(shù)組里。特別指出，每一個鏈自由基都需用四個參數(shù)用來標記；標號，鏈長，自由基種類，鏈段分布。(5)對步驟(4)得到的模擬結(jié)果進行統(tǒng)計，得到序列分布隨時間變化的曲線。(6)根據(jù)得到的序列分布結(jié)果考察預設(shè)的微管式反應器的工藝條件的合理性，重新設(shè)定工藝條件，重復步驟(3)、(4)直到得到滿意的序列分布結(jié)果和相對應的分批加料的工藝條件，將這個工藝條件應用于共聚反應的試驗中。本發(fā)明模擬得到的結(jié)果如圖2、圖3、圖4、圖5所示。本實施例中，物料從進口進入微管，在微管中歷經(jīng)和微管長度相等的路程，排出微管外。排出微管外的產(chǎn)物組成，分子量，分子量分布，平均序列長度，序列分布會達到一個定態(tài)，定態(tài)由投料量，投料比，微管中的流速、微管管長決定。任意調(diào)節(jié)一個或幾個參量，可以改變產(chǎn)物的定態(tài)。從圖可知，當100°C時，微管:L:0.6m ；D:0.01m;
流速-.Van= 0.000001* (67.19/0.8)/3600 L,
Vita= 0.000001*59.0/3600 L,
Vma= 0.000001* (2.24/0.96)/3600 L,
Voxide= 0.000001* (24.55*0.72)/3600 L,·Vreduce= 0.000001* (84.57*0.72)/3600 L,
Vwater= 0.000001*33.0/3600 L,
Vtotal=7.136Ε-8 L，
(AN:1TA:MA:Oxide:Reduce=800:4:17:1:3)
且各種物種的分散系數(shù)都為O時，產(chǎn)物的定態(tài)為，轉(zhuǎn)化率為37%，鏈長達到430，鏈長分布系數(shù)為1.93，平均序列長度為90，序列長度分布為2.3。
權(quán)利要求
1.一種蒙特卡洛模擬預測微管式反應器內(nèi)自由基共聚的方法，其特征在于具體步驟如下: (1)確定聚合體系中所有的基元反應和擴散反應，從聚合物手冊中查出各速率常數(shù)的數(shù)值，確定各種基元反應的速率常數(shù)K.和擴散速率常數(shù)Kdiff—ma。，再將所述的所有的基元反應的宏觀速率常數(shù)K.轉(zhuǎn)化為蒙特卡洛模擬的微觀反應速率常數(shù)Kmta。、擴散反應的宏觀擴散速率常數(shù)Kdiff 。轉(zhuǎn)化為蒙特卡洛模擬的微觀速率常數(shù)Kdiff mi_，其中:對于一級反應，Kniac=Kniic^ ;對于不同物種之間的二級反應，Kmicro=Kmac/(VNa);對于同一物種之間的二級反應，Kmicro=2Kmac/(VNa)，Na為阿伏伽德羅常數(shù)，大小為6.02X IO23, V為反應體系的體積，Kmicro=Kmac*n2/L2, η為子體系的個數(shù)，L為體系的總長度； (2)預設(shè)微管式反應器的工藝條件:包括微管式反應器的容積、加入管中反應的不同原料的進料量和管內(nèi)流體的流速； (3)基于各種基元反應和預設(shè)微管式反應器工藝條件，使用蒙特卡洛算法，建立自由基共聚反應的算法，其動力學的蒙特卡洛模擬步驟為: Ca)輸入步驟(I)中所有基元反應的微觀反應速率常數(shù)及所有種類化學物種的擴散微管擴散系數(shù)，將時間t設(shè)定為O; (b)將微管分為η個單元，計算每個單元的容積、并計算每個單元內(nèi)的所有種類的化學物種的分子個數(shù)，和各種基元反應、擴散反應的速率以及速率加和，計算各反應概率和各單元反應的總概率； (c)產(chǎn)生兩個單位區(qū)間內(nèi)均勻分布的隨機數(shù)，判斷下一步反應位于哪一單元； (d)產(chǎn)生兩個單位區(qū)間內(nèi)均勻分布的隨機數(shù)，判斷下一步將發(fā)生何種基元反應或者何種擴散反應，相應的調(diào)整反應所涉及的化學物種的分子個數(shù)，計算時間間隔，累加到時間t上去; (e)根據(jù)步驟(2)中的所述的預設(shè)的微觀式反應器的工藝條件，在設(shè)定的時間點上加入一定比例的原料改變相應的化學物種的分子個數(shù)，并根據(jù)流速計算因為流動而改變的化學物種的分子個數(shù)； (f)多次重復步驟(b)(c) (d) (e)，即得到各化學物種隨時間的變化情況；將所述化學物種的變化分別記錄在不同數(shù)組中，并進行統(tǒng)計，得到模擬結(jié)果； 其中:步驟(3)中所述化學物種為引發(fā)劑、單體、進入共聚物的單體和共聚物；步驟(f)中所述將所述化學物種的變化記錄在不同數(shù)組中，即將各種單體分子的個數(shù)，引發(fā)劑分子的個數(shù)，各種鏈自由基的個數(shù)，各種鏈自由基的長度，共聚物中不同單體的鏈段的長度，聚合物的分子量分別儲存到相應的數(shù)組里，其中各個鏈自由基都用下述四個參數(shù):標號、鏈長、自由基種類和鏈段分布來標記； (4)對步驟(3)得到的模擬結(jié)果進行統(tǒng)計，得到序列分布隨時間變化的曲線； (5)根據(jù)得到的序列分布結(jié)果考察預設(shè)的微管式反應器的工藝條件的合理性，重新設(shè)定工藝條件，重復步驟(3)、(4)直到得到滿意的序列分布結(jié)果和相對應的分批加料的工藝條件，將這個工藝條件應用于共聚反應的試驗中。
2.根據(jù)權(quán)利要求1所述的方法，其特征在于:步驟(4)中，對模擬結(jié)果進行統(tǒng)計，其具體方法如下: (a)通過跟蹤統(tǒng)計所有單體的變化， 推導出微管式反應器的工藝條件對反應體系的動力學的影響； (b)通過跟蹤統(tǒng)計已經(jīng)進入共聚物的單體的變化，推導出微管式反應器的工藝條件對共聚物組成的影響； Ce)通過跟蹤統(tǒng)計每個聚合物鏈長所有的鏈段長度，推導出微管式反應器的工藝條件對共聚物鏈段序列分布的影響。
3.根據(jù) 權(quán)利要求1所述的方法，其特征在于:所述方法應用于三元或以上聚合反應中。
全文摘要
本發(fā)明涉及高分子聚合技術(shù)領(lǐng)域，具體為一種蒙特卡洛模擬預測微管式反應器內(nèi)自由基共聚的方法。本發(fā)明從共聚反應的基元反應出發(fā)，描述單體、引發(fā)劑等原料進入微管中反應一段時間后流出的全過程，用模擬的方法得到流出的產(chǎn)物組分的組成，高分子鏈的聚合度和序列分布等隨時間的變化?？疾煳⒐苁椒磻鞯墓に嚄l件對聚合反應的影響，得到最佳的投料比。該方法簡便可行。
文檔編號G06F19/00GK103150477SQ20131008074
公開日2013年6月12日 申請日期2013年3月14日 優(yōu)先權(quán)日2013年3月14日
發(fā)明者郜澤慧, 何軍坡, 張紅東, 楊玉良 申請人:復旦大學