聚合

聚合是小分子结合形成大分子(包含原始分子的重复结构单元)的化学过程。

通常情况下,在溶剂中,如果化学反应中同时存在引发剂,单体分子会结合在一起,形成聚合物链或 3D 网络。 聚合涉及许多反应步骤,如引发、传播、分支、终止。 将形成许多不同链长的自由基。 最终产物是包含不同长度和结构的聚合物分子的混合物。

此特征基于聚合的矩量法模型。 可通过标准矩量法获得如聚合物分子链长分布等主要属性。

Simcenter STAR-CCM+ 求解三个活性聚合物(即 λ 0 λ 1 λ 2 )和三个非活性聚合物(即 μ 0 μ 1 μ 2 )的标量传输。 除了下述源项以外,Simcenter STAR-CCM+ 还可用于添加定义的矩源。 Simcenter STAR-CCM+ 还可以求解引发剂、单体、自由基、溶剂和改性剂的标量传输。

矩传输方程的源项取决于聚合的子过程:

  • 引发剂分解
    1. EQUATION_DISPLAY
    Iikdi2R,i=1,2,...,Ni
    (3781)
  • 链引发
    2. EQUATION_DISPLAY
    R.+MkIR1
    (3782)

    其中,

    可使用以下公式获得引发剂消耗率:

    3. EQUATION_DISPLAY
    rIi=-kdi[Ii],i=1
    (3783)

    并且使用以下公式获得主自由基生成率:

    4. EQUATION_DISPLAY
    rR.=i=1Nd2fikdi[Ii]-kI[R.][M]
    (3784)

    其中,

    f i = 引发剂效率

  • 传播
    5. EQUATION_DISPLAY
    Rx+MkpRx+1
    (3785)

    其中,可使用以下公式获得单体消耗率:

    6. EQUATION_DISPLAY
    rM=-kp[M]λ0
    (3786)
  • 链转移至单体
    7. EQUATION_DISPLAY
    Rx+MktmPx+R1
    (3787)
  • 链转移至溶剂
    8. EQUATION_DISPLAY
    Rx+SktsPx+R1
    (3788)

    其中,可使用以下公式获得溶剂消耗率:

    9. EQUATION_DISPLAY
    rS=-kts[S]λ0
    (3789)
  • 链转移至聚合物
    10. EQUATION_DISPLAY
    Rx+PyktsRy+Px
    (3790)
  • β 裂解
    11. EQUATION_DISPLAY
    Rx+PykβPx+Rz+Py-z=
    (3791)
  • 与端子双键的反应
    12. EQUATION_DISPLAY
    Rx+Py=kdbRx+y
    (3792)
  • 通过歧化反应终止
    13. EQUATION_DISPLAY
    Rx+RyktdPx+Py
    (3793)
  • 通过合成终止
    14. EQUATION_DISPLAY
    Rx+RyktcPx+y.
    (3794)

Simcenter STAR-CCM+ 可求解矩的标量传输方程,如下所示:

  • 对于活性聚合物

    可使用 r λ 0 作为源项求解 λ 0

    15. EQUATION_DISPLAY
    λ0=n=1Rn
    (3795)

    其中, r λ 0 由以下公式定义:

16. EQUATION_DISPLAY
rλ0=i=1Nd2fkdi[Ii]-(ktc+ktd)λ20
(3796)

可使用 r λ 1 作为源项求解 λ 1

17. EQUATION_DISPLAY
λ1=n=1nRn
(3797)

其中, r λ 1 由以下公式定义:

18. EQUATION_DISPLAY
rλ1=i=1Nd2fkdi[Ii]-(ktm[M]+kts[S])λ1+kp[M]λ0-(ktc+ktd)λ0λ1+ktp(λ0μ2-λ1μ1)+kdbλ0μ1+kβ(12λ0μ2-λ1μ1)
(3798)

可使用 r λ 2 作为源项求解 λ 2

19. EQUATION_DISPLAY
λ2=n=1n2Rn
(3799)

其中, r λ 2 由以下公式定义:

20. EQUATION_DISPLAY
rλ2=i=1Nd2fkdi[Ii]-(ktm[M]+kts[S])λ2+2kp[M]λ1-(ktc+ktd)λ0λ2+ktp(λ0μ3-λ2μ1)+kdb(2λ1μ1+λ0μ2)+kβ(13λ0μ3-12λ0μ2+16λ0μ1-λ2μ1)
(3800)
  • 对于非活性聚合物

    可使用 r μ 0 作为源项求解 μ 0

    21. EQUATION_DISPLAY
    μ0=n=1Pn
    (3801)

    其中, r μ 0 由以下公式定义:

    22. EQUATION_DISPLAY
    rμ0=(ktm[M]+kts[S])λ0+(12ktc+ktd)λ20-kdbλ0μ0+kβλ0μ1
    (3802)

可使用 r μ 1 作为源项求解 μ 1

23. EQUATION_DISPLAY
μ1=n=1nPn
(3803)

其中, r μ 1 由以下公式定义:

24. EQUATION_DISPLAY
rμ1=(ktm[M]+kts[S])λ1+(ktc+ktd)λ0λ1-ktp(λ0μ2-λ1μ1)-kdbλ0μ1+kβ(λ1μ1-12λ0μ2)
(3804)

可使用 r μ 2 作为源项求解 μ 2

25. EQUATION_DISPLAY
μ2=n=1n2Pn
(3805)

其中, r μ 2 由以下公式定义:

26. EQUATION_DISPLAY
rμ2=(ktm[M]+kts[S])λ2+(ktc+ktd)λ0λ2+ktcλ21-ktp(λ0μ3-λ2μ1)-kdbλ0μ2+kβ[λ2μ1-λ0(23μ3+12μ2-16μ1)]
(3806)

Pladis 和 Kiparissides [779] 采用较高矩 μ 3 的封闭。

除了用于上述聚合物标量和聚合物矩的源项及其雅可比以外,还可以包括用于聚合物标量和矩源的用户自定义的速率和雅可比。

根据指定的传播热和引发热,相应源项将添加到能量方程。

通过求解以上六个矩方程的传输,可以获得有关聚合物分子结构(如聚合物的多分散性和分子量分布)的信息。 Simcenter STAR-CCM+ 为多分散性指数、聚合物平均数链长和聚合物平均质量链长提供后处理场函数。 所有活性和非活性矩还可用作后处理场函数。

多分散性指数由以下公式给定:

27. EQUATION_DISPLAY
μ 2 μ 0 μ 2 1
(3807)

聚合物平均数链长 (NACL) 由以下公式给定:

28. EQUATION_DISPLAY
μ 1 μ 0
(3808)

聚合物平均质量链长 (MACL) 由以下公式给定:

29. EQUATION_DISPLAY
μ 2 μ 1
(3809)