Kinetics and thermodynamics during pyrolysis of scrap printed circuit board by TGA
-
摘要: 废旧电路板(SPCB)是一种典型的有机废弃物,可通过热解技术实现其资源化利用。采用热重分析技术(TGA)对其热解特性进行研究,揭示热解过程反应动力学和热力学。实验在氮气气氛下,考察了不同升温速率(5、10、15 ℃/min)对SPCB热失重特性的影响,结果表明热解过程主要发生在250 ~ 400 ℃温度区间,随着升温速率增大,SPCB热失重(TG)曲线逐渐向高温方向偏移,在对应的热失重速率(DTG)曲线中,存在一个明显的失重峰,且峰值温度不断增加,热滞后现象显著。采用Flynn-Wall-Ozawa(FWO)模型、Kissinger-Akahira-Sunose(KAS)模型和Friedman(FM)模型进行动力学分析,拟合得到平均表观活化能(Ea)分别为168.46、167.31、234.84 kJ/mol,活化能均随转化率增加而相应增大。利用FWO模型对热力学参数进行计算,在相同升温速率下,随着转化率的增大,吉布斯自由能变(ΔG)逐渐降低,对应的焓变(ΔH)和熵变(ΔS)不断增加;在相同转化率时,ΔH和ΔS随升温速率增加稍有降低,而ΔG逐渐增加。Abstract: Scrap printed circuit board (SPCB) is a typical organic waste, which could be utilized as a resource by pyrolysis technology. The pyrolysis characteristics of SPCB were studied by thermogravimetric analysis (TGA) to reveal the reaction kinetics and thermodynamics during the pyrolysis process. Under N2 atmosphere, the effects of different heating rates, e.g. 5 ℃/min, 10 ℃/min and 15 ℃/min, on the thermal decomposition behavior of SPCB were investigated in detail. The results observed showed that the pyrolysis process was mainly occurred in the temperature range of about 250 ℃ to 400 ℃. With the increase of heating rates, the thermogravimetry (TG) curves of SPCB was gradually deviated to the high temperature direction. For the derivative thermogravimetry (DTG) curves, there was an obvious mass loss peak, and the corresponding peak temperature kept increasing, indicating that there was a significant thermal hysteresis phenomenon. The Flynn-Wall-Ozawa (FWO) model, Kissinger-Akahira-Sunose (KAS) model and Friedman (FM) model were used to analyze the thermogravimetric data. The average apparent activation energies (Ea) were 168.46 kJ/mol, 167.31 kJ/mol and 234.84 kJ/mol, respectively, which increased with the increase of conversion rate. The thermodynamic parameters were calculated using FWO model. Under the same heating rate, the Gibbs free energy change (ΔG) was gradually decreased with the increase of conversion rate, while the corresponding enthalpy change (ΔH) and entropy change (ΔS) were all increased. Meanwhile, at the same conversion, both of ΔH and ΔS decreased slightly with the increase of heating rate, while ΔG was increased gradually.
-
Keywords:
- scrap printed circuit board /
- pyrolysis characteristics /
- kinetics /
- thermodynamics
-
电磁搅拌利用线圈激发的交变磁场渗透到金属熔体中, 在其中感生起感应电流, 感应电流与当地磁场互相作用产生电磁力。该搅拌方式不与坯体直接接触, 却能够有效的改善铸坯表面、皮下以及内部质量。目前在金属冶炼行业应用相对广泛的电磁搅拌器, 基本可以按照搅拌区域和激发磁场形态进行分类, 其中各模式磁场的电磁搅拌工艺针对不同的铸坯形状和材料特性有不同的效果。近年来, 许多钢铁企业、有色金属生产企业以及各大科研院所着力于开发新型的复合型材料, 这些材料在具有高性能的同时也对冶金工艺提出了新的要求[1-2]。单一磁场形态的电磁搅拌已经难以达到高性能复合材料的冶金要求, 由此湖南理工学院和湖南科美达电气股份有限公司联合设计开发了一款多模式磁场电磁搅拌器, 用来产生旋转磁场、行波磁场以及螺旋磁场, 可根据不同材料特性或铸坯形状, 在金属冶炼过程中选择不同的磁场形态进行电磁搅拌, 以求达到最佳的冶金效果。
电磁搅拌的过程受到磁流体力学和传热学等诸多领域的耦合作用[3-5], 是一个非常复杂的过程, 虽然多方面因素影响着电磁搅拌效果, 但整体而言搅拌磁场模式、励磁电流强度、励磁电流频率以及搅拌的时间是制约电磁搅拌效用的主要参数。为了研究分析各模式磁场的分布规律以及各主要工艺参数对磁场的影响, 本文利用有限元分析软件Ansoft Maxwell对多模式磁场电磁搅拌器建立了有限元数值分析模型, 主要对磁场中的磁感应强度和电磁力的分布特征以及数值变化规律进行研究, 以验证本设计的有效性, 并为后续设计提供参考。
1 多模式磁场电磁搅拌器的结构与原理
多模式磁场电磁搅拌器主体是由非磁性外壳、铁芯、旋转线圈绕组以及行波线圈绕组组成, 在铁芯和线圈之间加入了压圈和连接片等固定结构的零件。其装配图见图 1。共有12个线圈绕组, 其中6个线圈绕组垂直缠绕在铁芯上, 为线圈Ⅰ, 负责激发旋转磁场; 另外6个线圈绕组水平嵌于铁芯中, 为线圈Ⅱ, 负责激发行波磁场。2组线圈同时输入励磁电流, 则激发出螺旋磁场, 螺旋磁场的螺旋角度由旋转磁场线圈电流和行波磁场线圈电流的大小和方向共同确定。
2 基于Ansoft Maxwell的有限元分析
2.1 有限元仿真模型的建立
在保证计算精度的前提下, 为了提高模型的收敛性以及简化计算过程, 多模式磁场电磁搅拌器的计算模型采用如下假设条件:
1) 铁芯附近的磁感应强度分布均匀, 从截面穿出的磁力线垂直于截面;
2) 只计算铁芯表面部分对于搅拌器内腔磁场产生影响;
3) 结晶器内的金属为熔融状态;
4) 不考虑结晶器内流场变化的影响;
5) 结晶器与熔体之间有绝缘间隙。
以多模式磁场电磁搅拌器实验平台的结构参数为参考, 忽略细小结构件以及细小间隙构建三维几何模型, 以达到减小数值仿真模型的计算量和增加数值仿真模型收敛性的目的[6]。多模式磁场电磁搅拌器的结构参数如表 1所列, 结构模型如图 2所示。
表 1 电磁搅拌器结构尺寸与相关参数Table 1. Structure dimension and relevant parameters of electromagnetic stirrer2.2 磁场计算的数学模型
麦克斯韦方程组(Maxwell Equations)为描述电场、磁场与电流密度之间关系的方程组, 通过这个方程组可以描述出电磁场的规律。
(1) (2) (3) (4) 式(1)法拉第电磁感应定律, 说明了磁感应强度随时间的变化产生电磁场强度的环流; 式(2)安培定律, 表明传导电流密度可以产生磁场的环流, 其方向按右螺旋法则; 式(3)高斯定律, 表明磁力线皆为首尾相接的闭合曲线, 没有初始点, 也没有终止点[7-9]; 式(4)状态方程, 表明磁感应强度与磁场强度之间的关系; 结合总体方程组能表示出变化的磁场和电场之间相互作用形成稳态电磁场的基本原理[10-11]。利用这组方程理论上可以解决大部分宏观的电磁场问题。
2.3 材料属性定义与激励源设置
此分析为电磁场分析, 导磁区域为线圈、铁芯、空气、结晶器和金属熔体, 导电区域为线圈、结晶器和金属熔体, 材料属性定义见表 2。
表 2 有限元模型中材料的物理性质参数Table 2. Physical property parameters of materials in the finite element modelWD315-50硅钢片的相对磁导率的曲线如图 3所示。
图 3为WD315-50硅钢片在磁化过程中磁感强度(B)与磁场强度(H)之间关系的曲线, 当H逐渐增大时, B急剧增加, 几乎呈直线上升, 当H进一步增大时, B的增加开始变得缓慢, 到达一定值后, B几乎不再增加, 即达到了饱和。
对于多模式磁场电磁搅拌器, 在对旋转磁场线圈绕组和行波磁场线圈绕组同时通入三相交流电激励后, 在搅拌器内腔中激发出呈正弦变换的螺旋磁场, 如果我们减去旋转磁场型励磁线圈绕组的激励, 则得到行波磁场型的磁场形态。相反, 如果减去行波磁场型励磁线圈绕组的激励, 则可得到旋转磁场型的磁场形态[12]。若三相交流电之间的相位差为120°, 则施加至线圈截面上的电流激励可表示为:
(5) (6) (7) 式中:A1、A2、A3分别代表为第一相、第二相和第三相电流强度; A0为电流强度幅值; w为电源角频率, w=2πf; f为频率; t为时间[13-15]; 电磁场采用瞬态磁场计算方法, 时间步长设置为一个周期的1/20, 计算时长为两个周期。
3 多模式电磁搅拌数值模拟
在搭建完多模式磁场电磁搅拌器的有限元分析模型后, 通过改变电流幅值、电流频率和磁场类型等计算参数, 完成了多组不同计算参数的仿真计算。选取搅拌器内结晶器中心位置为原点, 以此选取结晶器内轴向和径向的磁感应强度数值和电磁力数值, 进行对比分析。
3.1 实测结果验证
对应多模式磁场电磁搅拌器中的旋转磁场和行波磁场, 利用高斯计对其结晶器中心轴线上以中心轴线坐标原点为中心等距离的5个点进行磁感应强度的测量, 之后利用数值仿真模型对相同位置的磁感应强度进行数值计算, 两者的结果对比如图 4、图 5所示。
从图 4、图 5中可以看出, 高斯计测量结果与数值模拟计算结果偏差较小, 两者基本一致, 说明所建数值分析模型是准确和可靠的[16]。由于螺旋磁场是由旋转磁场和行波磁场在空间中的叠加形成, 所以所建模型同样能对螺旋磁场进行准确可靠的数值模拟仿真。
3.2 各模式磁场矢量分析
选取电磁搅拌工艺参数为200 A、2 Hz, 分别对旋转磁场、行波磁场以及螺旋磁场进行数值仿真模拟, 截取结晶器中心位置横截面以及纵截面磁感应强度矢量图, 得出各模式磁场矢量分布特征, 见图 6、图 7和图 8。
由图 6可知, 旋转磁场磁感应强度呈现出的变换规律为成对变换, 电流随时间变化, 磁场波则按三相交流电相序沿圆周旋转, 从而导致金属溶液做旋转运动。
由图 7可知, 行波磁场磁感应强度呈现出的变化规律为环形纵向正弦变换且左右对称, 这种磁场分布使金属溶液的流速在纵向呈正弦变换, 从而形成一个左右对称的大环流, 流动中心与结晶器轴线重合。
由图 8可知, 螺旋磁场可以看成旋转磁场与行波磁场在空间上的叠加, 其方向由旋转磁场方向和行波磁场方向共同确定, 可以通过改变旋转磁场或者行波磁场的线圈电流的大小和方向来改变螺旋磁场的螺旋角[17-18]。
3.3 电流强度对磁场分布的影响
选取结晶器中轴线和结晶器中心横截面径向的磁感应强度数据, 电流频率固定为4 Hz, 分析各磁场模式下, 电流强度为50, 150, 200, 250 A时对应的磁感应强度分布情况。
图 9所示为搅拌频率为4 Hz时, 不同磁场模式下励磁电流强度对应结晶器中轴线以及径向各处磁感应强度值的影响。结果表明:随着励磁电流的增大, 不同磁场模式下的磁感应强度均有所增加, 不过磁场分布结构基本不改变; 相同的励磁电流下, 轴向分布呈现中心位置附近达到最大值, 向两侧逐渐减弱, 径向分布呈现出由中心向结晶器边缘逐渐增强的趋势[19-20]。
3.4 搅拌频率对磁场分布的影响
选取结晶器中轴线和结晶器中心横截面径向的实验数据, 分析各磁场模式下, 电流强度固定为200 A, 电流频率为0.5, 1, 2, 4, 6, 10 Hz时对应的磁感应强度分布如图 10所示。
图 10所示为励磁电流强度为200 A时, 不同磁场模式下结晶器中轴线以及径向各处磁感应强度随励磁电流频率的变化趋势。由图 10可见, 无论在哪种搅拌方式下, 在频率在0~2 Hz区间, 频率越大, 结晶器内的磁感应强度越大。而在2~10 Hz区间, 频率越大, 结晶器内的磁感应强度反而越小, 在2 Hz左右, 磁感应强度达到最大值。
3.5 电磁力计算结果分析
以前面定义的多模式电磁搅拌数值分析模型为例, 稀土铝熔体作为金属熔体, 选取200 A、2 Hz工艺参数进行仿真模拟, 计算结果分析如图 11和12所示。
由图 11可知, 铸坯的中心轴向电磁力在轴向上呈"两端小、中间大"的分布特征; 而铸坯在中心横截面上切向电磁力的径向分布, 呈现出从铸坯表面向中心缩减的特征, 铸坯中心受到的电磁力较小。
铸坯横截面电磁力随电流强度的变化曲线如图 12所示。可以看出, 电磁力的大小与电流强度大小呈正比关系; 电磁力的大小与电流频率的关系在0~3.5 Hz区间, 由于结晶器内磁感应强度有一个快速上升的过程, 且电流频率的增大导致感生电流渗入深度缩减至结晶器内有效工作区域, 导致电磁力大小随着频率增大而增大。而在3.5~10 Hz区间, 由于金属熔体受到集肤效应的影响[21], 熔体内的感生电流对金属熔体的渗入深度继续减小, 逐渐退出结晶器内有效工作区域, 搅拌器内磁感应强度也逐渐减弱, 随之产生的电磁力也随着频率增加而减小, 即电磁搅拌效率降低。
4 结束语
采用数值模拟仿真和实测相结合的方法, 对多模式磁场电磁搅拌器的磁场分布情况进行比较。对比分析了不同磁场模式下, 结晶器内的磁感应强度分布以及金属熔体所受到的电磁力情况, 得到如下结论:
1) 磁场分布特征:3种磁场类型作用下, 不同时刻的磁场分布形式相似, 只是磁场的相位发生了变化; 结晶器中心轴线磁感应强度数值的变化皆为"两头小, 中间大", 沿半径方向磁感应强度数值的变化则为结晶器边缘处最大, 向中心方向不断衰减, 与实际测量结果一致。
2) 电磁力分布规律:3种磁场类型作用下, 径向电磁力边缘最大, 逐渐向中心衰减; 轴向电磁力中下部较大, 沿轴线向两端逐渐减小。
3) 结晶器内磁感应强度在常规工况参数区间内, 与励磁电流强度成正比, 励磁电流强度越大则磁感应强度越大。磁感应强度与励磁电流频率的关系则表现为, 在0~2 Hz区间成正比关系, 在2~10 Hz区间成反比关系, 整个系统在2 Hz左右得到最大磁感应强度值。
4) 结晶器内金属熔体受到的电磁力大小与励磁电流强度成正比。在0~3.5 Hz区间, 频率越大, 结晶器内的磁感应强度越大; 而在3.5~10 Hz区间, 频率越大, 结晶器内的电磁力反而越小, 整个系统在3.5 Hz左右获得电磁力最大值。
王庆龙 -
表 1 样品物化特征分析
Table 1 Analysis of physical and chemical characteristics of samples
样品 工业分析/(ad, %) 元素分析/(daf, %) 高位热值/(kJ/g) 固定碳* 水分 挥发分 灰分 C H N S O* SPCB 6.15 0.56 34.27 59.02 55.25 5.55 1.59 0.48 37.13 9.59 注: *为差减法求算;ad为空气干燥基;daf为干燥无灰基。表 2 KAS模型、FWO模型和FM模型热解动力学参数结果
Table 2 Pyrolysis kinetics parameters of KAS, FWO and FM model
α KAS 模型 FWO模型 FM模型 Eα/(kJ/mol) A/s-1 R2 Eα/(kJ/mol) A/s-1 R2 Eα/(kJ/mol) A/s-1 R2 0.1 121.30 6.20×108 0.99 124.46 1.46×109 0.99 134.29 4.11×109 0.99 0.2 123.15 8.00×108 0.99 126.31 1.86×109 0.99 141.95 2.34×1010 0.99 0.3 127.64 1.98×109 0.99 130.63 4.33×109 0.99 157.86 7.85×1011 0.99 0.4 133.31 6.31×109 0.99 136.07 1.29×1010 0.99 177.15 4.23×1013 0.99 0.5 144.23 6.06×1010 0.99 146.52 1.07×1011 0.99 201.76 5.81×1015 0.99 0.6 163.36 2.97×1012 0.99 164.78 4.20×1012 0.99 235.74 4.20×1018 0.99 0.7 200.99 5.24×1015 0.99 200.68 5.09×1015 0.99 294.57 2.93×1023 0.99 0.8 324.48 1.25×1026 0.97 318.27 4.96×1025 0.98 535.40 2.09×1043 0.99 平均 167.31 168.46 234.84 表 3 各文献中SPCB热解过程活化能
Table 3 Eaof SPCB pyrolysis from various literature
表 4 不同升温速率下SPCB热解过程热力学参数
Table 4 Thermodynamic parameters of SPCB pyrolysis at different heating rates
升温速率 ) 5 ℃ /min 0.1 119.74 171.21 -89.24 7.01×108 0.2 121.54 171.14 -85.99 1.05×109 0.3 125.83 170.98 -78.29 2.66×109 0.4 131.24 170.79 -68.56 8.63×109 0.5 141.65 170.43 -49.91 8.21×1010 0.6 159.85 169.87 -17.37 4.16×1012 0.7 195.68 168.92 46.39 9.05×1015 0.8 313.14 166.71 253.87 6.40×1026 平均值 163.58 170.01 -11.14 8.00×1025 10 ℃ /min 0.1 119.63 172.76 -89.72 6.78×108 0.2 121.43 172.68 -86.56 1.00×109 0.3 125.71 172.52 -79.04 2.49×109 0.4 131.13 172.32 -69.56 7.84×109 0.5 141.54 171.95 -51.36 7.04×1010 0.6 159.76 171.38 -19.62 3.23×1012 0.7 195.60 170.41 42.55 5.79×1015 0.8 313.09 168.13 244.79 2.16×1026 平均值 163.49 171.52 -13.57 2.71×1025 15 ℃ /min 0.1 119.55 174.10 -90.50 6.27×108 0.2 121.35 174.03 -87.39 9.20×108 0.3 125.63 173.86 -80.00 2.25×109 0.4 131.05 173.66 -70.67 6.96×109 0.5 141.47 173.29 -52.78 6.02×1010 0.6 159.69 172.70 -21.57 2.59×1012 0.7 195.55 171.71 39.55 4.07×1015 0.8 313.05 169.40 238.30 9.99×1025 平均值 163.42 172.84 -15.63 1.25×1025 -
[1] HE W Z,LI G M,MA X F,et al. WEEE recovery strategies and the WEEE treatment status in China[J]. Journal of Hazardous Materials,2006,136(3):502-512.
[2] ANNAMALAI M,GURUMURTHY K. Characterization of end-of-life mobile phone printed circuit boards for its elemental composition and beneficiation analysis[J]. Journal of the Air & Waste Management Association,2021,71(3):315-327.
[3] HEACOCK M,KELLY C B,ASANTE K A,et al. E-waste and harm to vulnerable populations:a growing global problem[J]. Environmental Health Perspectives,2016,124(5):550-555.
[4] ZENG X L,GONG R Y,CHEN W Q,et al. Uncovering the recycling potential of “new” WEEE in China[J]. Environmental Science & Technology,2016,50(3):1347-1358.
[5] 王芳芳,赵跃民,张涛,等. 废弃线路板中金属资源的物理回收[J]. 矿产综合利用,2017(2):1-7. [6] 朱振华,沈天晓. 国内废线路板处理现状与思考[J]. 冶金管理,2022(17):12-14. [7] BIDINI G,FANTOZZI F,BARTOCCI P,et al. Recovery of precious metals from scrap printed circuit boards through pyrolysis[J]. Journal of Analytical and Applied Pyrolysis,2015,111:140-147.
[8] CHAUHAN G,JADHAO P R,PANT K K,et al. Novel technologies and conventional processes for recovery of metals from waste electrical and electronic equipment:challenges & opportunities-A review[J]. Journal of Environmental Chemical Engineering,2018,6(1):1288-1304.
[9] 刘书杰,申士富,叶力佳,等. 印刷电路板资源综合利用技术研究[J]. 矿冶,2011,20(2):101-104. [10] ZENG X L,MATHEWS J A,LI J H. Urban mining of E-waste is becoming more cost-effective than virgin mining[J]. Environmental Science & Technology,2018,52(8):4835-4841.
[11] 王芳,刘玉卿,海热提. 防溴型环氧树脂电路板热解动力学分析[J]. 化工学报,2011,62(10):2945-2950. [12] CAO R,ZHOU R S,LIU Y Q,et al. Research on the pyrolysis characteristics and mechanisms of waste printed circuit boards at fast and slow heating rates[J]. Waste Management,2022,149:134-145.
[13] WU D,DING G Z,CHI W F,et al. Research on the pyrolysis kinetics of resin powder on waste printed circuit board with different particle sizes at different heating rates:inspiration for the pyrolysis mechanism[J]. Journal of Thermal Analysis and Calorimetry,2022,147(14):8047-8059.
[14] ALENEZI R A,AL-FADHLI F M. Thermal degradation kinetics of waste printed circuit boards[J]. Chemical Engineering Research and Design,2018,130:87-94.
[15] 刘法球,周少基,向敏,等. 滤泥的热解特性及动力学分析[J]. 南方农业学报,2019,50(11):2552-2558. [16] 刘承飞,李江平,刘大方,等. 废旧电脑印刷电路板的热解特性及动力学分析[J]. 有色金属科学与工程,2022,13(1):38-43. [17] OZAWA T. A new method of analyzing thermogravimetric data[J]. Bulletin of the Chemical Society of Japan,1965,38(11):1881-1886.
[18] FLYNN J H,WALL L A. A quick,direct method for the determination of activation energy from thermogravimetric data[J]. Journal of Polymer Science Part B:Polymer Letters,1966,4(5):323-328.
[19] KISSINGER H E. Reaction kinetics in differential thermal analysis[J]. Analytical Chemistry,1957,29(11):1702-1706.
[20] FRIEDMAN H L. Kinetics of thermal degradation of char-forming plastics from thermogravimetry. Application to a phenolic plastic[J]. Journal of Polymer Science Part C:Polymer Symposia,2007,6(1):183-195.
[21] KIM Y S,KIM Y S,KIM S H. Investigation of thermodynamic parameters in the thermal decomposition of plastic waste-waste lube oil compounds[J]. Environmental Science & Technology,2010,44(13):5313-5317.
[22] OLSZAK-HUMIENIK M,MOZEJKO J. Thermodynamic functions of activated complexes created in thermal decomposition processes of sulphates[J]. Thermochimica Acta,2000,344(1/2):73-79.
[23] QUAN C,LI A M,GAO N B. Research on pyrolysis of PCB waste with TG-FTIR and py-GC/MS[J]. Journal of Thermal Analysis and Calorimetry,2012,110(3):1463-1470.
[24] SHEN Y F,CHEN X M,GE X L,et al. Thermochemical treatment of non-metallic residues from waste printed circuit board:pyrolysis vs. combustion[J]. Journal of Cleaner Production,2018,176:1045-1053.
[25] CHEN Y,YANG J K,ZHANG Y,et al. Kinetic simulation and prediction of pyrolysis process for non-metallic fraction of waste printed circuit boards by discrete distributed activation energy model compared with isoconversional method[J]. Environmental Science and Pollution Research,2018,25(4):3636-3646.
[26] CHEN S Y,LI R,SHEN Y Q,et al. Self-catalytic pyrolysis thermodynamics of waste printed circuit boards with co-existing metals[J]. Frontiers of Environmental Science & Engineering,2022,16(11):146.
[27] 陈斌. 废旧电路板金属资源分级分离回收及过程控制研究[D]. 合肥:中国科学技术大学,2020. [28] CHEN C X,MA X Q,HE Y. Co-pyrolysis characteristics of microalgae Chlorella vulgaris and coal through TGA[J]. Bioresource Technology,2012,117:264-273.
[29] SHARARA M, SADAKA S.Thermogravimetric analysis of swine manure solids obtained from farrowing, and growing-finishing farms[J].Journal of Sustainable Bioenergy Systems, 2014, 4(1): 75-86.
[30] 郭晓娟,张于峰,秦贯丰. 废弃印刷线路板热解表观动力学模型研究[J]. 化工进展,2012,31(11):2395-2399. [31] ABD-ELGHANY M,KLAPÖTKE T M,ELBEIH A. Investigation of 2,2,2-trinitroethyl-nitrocarbamate as a high energy dense oxidizer and its mixture with Nitrocellulose (thermal behavior and decomposition kinetics)[J]. Journal of Analytical and Applied Pyrolysis,2017,128:397-404.
[32] TURMANOVA S C,GENIEVA S D,DIMITROVA A S,et al. Non-isothermal degradation kinetics of filled with rise husk ash polypropene composites[J]. Express Polymer Letters,2008,2(2):133-146.
[33] VUPPALADADIYAM A K,LIU H,ZHAO M,et al. Thermogravimetric and kinetic analysis to discern synergy during the co-pyrolysis of microalgae and swine manure digestate[J]. Biotechnology for Biofuels,2019,12:170-187.
[34] XU Y L,CHEN B L. Investigation of thermodynamic parameters in the pyrolysis conversion of biomass and manure to biochars using thermogravimetric analysis[J]. Bioresource Technology,2013,146:485-493.