Молекулярно-динамічне моделювання розподілу газоподібних продуктів піролізи поліетилену високої густини з використанням силового поля ReaxFF
983–997 (2025)
PACS numbers: 02.70.Ns, 07.05.Tp, 31.15.-p, 34.20.Gj, 36.20.-r, 82.30.Lp, 82.35.Np
Отримано 4 квітня 2025 р.; після доопрацювання — 16 травня 2025 р.
Об'єктом дослідження є процес піролізи поліетилену високої густини (ПВГ). Метою роботи є дослідження якісного та кількісного розподілів цінних і шкідливих газоподібних продуктів піролізи ПВГ та визначення раціональних температурних режимів перебігу процесу з максимальним виходом цінних продуктів і мінімізацією шкідливого впливу на навколишнє середовище. Дослідження ґрунтуються на методах молекулярно-динамічного (МД) моделювання у нанорозмірній обчислювальній комірці з використанням реакційного силового поля ReaxFF і спеціялізованих програмних продуктів Materials Studio i LAMMPS. Методика досліджень включає: створення в Materials Studio атомістичного моделю полімерного ланцюга з додаванням молекул Оксигену, проведення його геометричної оптимізації та мінімізації енергії й врівноваження у канонічному ізохорно-ізотермічному ансамблі (NVT) за температури у 300 К із кроком інтегрування по часу в 1 фс упродовж 5 нс; моделювання процесу піролізи в LAMMPS за температур у 600°C, 800°C, 1000°C, 1400°C, 1800°C, 2000°C. За результатами МД-моделювання піролізи побудовано графіки виходу кількости молекул основних газоподібних продуктів залежно від часу, що включають O2, H2, CO, C2H2, C2H4, C2H6, CH2O, CH4, CO2. Показано, що з наближенням концентрації Оксигену до мінімуму більш значущими стають зміни у концентраціях проміжних і кінцевих газоподібних продуктів. Серед основних газоподібних продуктів піролізи найбільший вихід (≥10 молекул) мають такі компоненти: водень (H2), метан (CH4), етилен (C2H4), ацетилен (C2H2) та монооксид Карбону (СО). Проведено оптимізацію температурного режиму піролізи ПВГ за критерієм максимального виходу основних газоподібних продуктів, за результатами якої визначено, що за температури у 1000°C відбувається достатньо інтенсивне виділення H2, CH4, C2H4 та СО, а подальше підвищення температури до 1400°C, 1800°C та 2000°C приводить до збільшення виходу H2 та зменшення виділення СО й деяких інших вуглеводнів (наприклад, CH4). Тому температуру у 1000°C взято за оптимальну. Одержані результати сприятимуть більш глибшому розумінню процесів піролізного розкладання ПВГ на нанорівні та розробці ефективних методів утилізації полімерних відходів.
КЛЮЧОВІ СЛОВА: піроліза, поліетилен високої густини, молекулярна динаміка, силове поле ReaxFF, моделювання, газоподібні продукти
ЛІТЕРАТУРА
- BlueWeaveConsulting, High Density Polyethylene HDPE Market; https://www.blueweaveconsulting.com/report/high-density-polyethylene-hdpe-market
- S. D. Anuar Sharuddin, F. Abnisa, W. M. A. Wan Daud, and M. K. Aroua, Energ. Convers. Manage., 115: 308 (2016); https://doi.org/10.1016/j.enconman.2016.02.037
- N. Jasminská, T. Brestovič, and M. Čarnogurská, Acta Mech. Autom., 7, No. 1: 20 (2013); https://doi.org/10.2478/ama-2013-0004
- Y.-F. Huang, P.-T. Chiueh, W.-H. Kuan, and S.-L. Lo, Energy, 89: 974 (2015); https://doi.org/10.1016/j.energy.2015.06.035
- Y. K. Choi, S.-J. Park, S. Park, S. Kim, N. R. Kern, J. Lee, and W. Im, J. Chem. Theory Comput., 17, No. 4: 2431 (2021); https://doi.org/10.1021/acs.jctc.1c00169
- K. Chenoweth, A. C. T. van Duin, and W. A. Goddard, J. Phys. Chem. A, 112, No. 5: 1040 (2008); https://doi.org/10.1021/jp709896w
- A. C. T. van Duin, S. Dasgupta, F. Lorant, and W. A. Goddard, J. Phys. Chem. A, 105, No. 41: 9396 (2001); https://doi.org/10.1021/jp004368u
- T. B. Y. Chen, A. C. Y. Yuen, B. Lin, L. Liu, A. L. P. Lo, Q. N. Chan, J. Zhang, S. C. P. Cheung, and G. H. Yeoh, J. Anal. Appl. Pyrol., 153: 104931 (2021); https://doi.org/10.1016/j.jaap.2020.104931
- BIOVIA Materials Studio; https://www.3ds.com/products/biovia/materials-studio
- LAMMPS Molecular Dynamics Simulator; https://www.lammps.org/
- I. V. Omelchuk and A. Y. Karvatskii, Visnyk of Vinnytsia Polytechnical Institute, 172, No. 1: 128 (2024) (in Ukrainian); https://doi.org/10.31649/1997-9266-2024-172-1-128-137
- J. Kong, K. Zhou, X. Ren, Y. Chen, Y. Li, and P. Meng, J. Anal. Appl. Pyrol., 169: Article ID 105847 (2023); https://doi.org/10.1016/j.jaap.2022.105847
- D. Hong, P. Li, T. Si, and X. Guo, Energy, 218: Article ID 119553 (2021); https://doi.org/10.1016/j.energy.2020.119553
- T. P. Senftle, S. Hong, Md. M. Islam, S. B. Kylasa, Y. Zheng, Y. K. Shin, C. Junkermeier, R. Engel-Herbert, M. J. Janik, H. M. Aktulga, T. Verstraelen, A. Grama, and A. C. T. van Duin, npj Comput. Mater., 2, Iss. 1: Article No. 15011 (2016); https://doi.org/10.1038/npjcompumats.2015.11
- Experimental Data for O₂ (Oxygen Diatomic) (2022); https://cccbdb.nist.gov/exp2x.asp?casno=7782447&charge=0
- msi2lmp Program Convert Tool LAMMPS Module Github Website; https://github.com/lammps/lammps/blob/develop/tools/msi2lmp/README
- R. L. C. Akkermans, N. A. Spenley, and S. H. Robertson, Mol. Simulat., 39, Nos. 14-15: 1153 (2013); https://doi.org/10.1080/08927022.2013.843775
- H. A. Posch, W. G. Hoover, and F. J. Vesely, Phys. Rev. A, 33, No. 6: 4253 (1986); https://doi.org/10.1103/physreva.33.4253
- H. M. Aktulga, J. C. Fogarty, S. A. Pandit, and A. Y. Grama, Parallel Comput., 38, Iss. 4-5: 245 (2012); https://doi.org/10.1016/j.parco.2011.08.005
- J. Lange, F. Gomes de Souza Jr., M. Nele, F. W. Tavares, I. S. V. Segtovich, G. C. Q. da Silva, and J. C. Pinto, Macromol. Theor. Simul., 25, No. 1: 45 (2016); https://doi.org/10.1002/mats.201500049