基于自适应变步长同伦法的循环流程收敛算法

doi:10.11949/0438-1157.20231400

化工学报 ›› 2024, Vol. 75 ›› Issue (7): 2604-2612.DOI: 10.11949/0438-1157.20231400

基于自适应变步长同伦法的循环流程收敛算法

李洪瑞¹(), 黄纯西¹, 洪小东²(), 廖祖维¹(), 王靖岱¹, 阳永荣¹

^1.化学工程联合国家重点实验室，浙江大学，浙江杭州 310027
^2.功能材料智能设计与制造浙江省工程研究中心，浙江大学杭州国际科创中心，浙江杭州 311215

收稿日期:2023-12-30 修回日期:2024-03-19 出版日期:2024-07-25 发布日期:2024-08-09
通讯作者: 洪小东,廖祖维
作者简介:李洪瑞（2001—），男，硕士研究生，22328110@zju.edu.cn
基金资助:
国家自然科学基金项目(U22A20415);浙江省尖兵领雁计划项目(2022C01SA442617);浙江省自然科学基金项目(LQ24B060001)

An adaptive variable-step homotopy-based algorithm for process simulation with cyclic streams

Hongrui LI¹(), Chunxi HUANG¹, Xiaodong HONG²(), Zuwei LIAO¹(), Jingdai WANG¹, Yongrong YANG¹

^1.State Key Laboratory of Chemical Engineering, Zhejiang University, Hangzhou 310027, Zhejiang, China
^2.Engineering Research Center of Functional Materials Intelligent Manufacturing of Zhejiang Province, ZJU-Hangzhou Global Scientific and Technological Innovation Center, Zhejiang University, Hangzhou 311215, Zhejiang, China

Received:2023-12-30 Revised:2024-03-19 Online:2024-07-25 Published:2024-08-09
Contact: Xiaodong HONG, Zuwei LIAO

摘要/Abstract

摘要：

包含循环流股的流程模拟面临初值难给定和传统算法收敛困难的问题。针对该问题提出了自适应变步长同伦算法，通过自适应调整同伦参数步长处理中间过程收敛失败的情形，提高收敛效率，并采用回溯线搜索法将迭代变量限制在定义域内，以提升算法的鲁棒性。构建了基于直接迭代法、Wegstein法、Broyden法和牛顿法的同伦算法，分析了同伦参数和辅助函数的影响。以淤浆法工艺生产高密度聚乙烯作为案例，模拟结果表明同伦算法可以提高模型在不同工况下的收敛性。商业流程模型软件Aspen Plus内置的求解器最高仅在21%测试工况下可以获得可行解，而同伦算法收敛案例占比可达88%。

关键词: 过程系统, 算法, 化学过程, 模拟, 同伦法

Abstract:

The simulation of processes with cyclic streams faces challenges in providing initial values and encounters difficulties in convergence with traditional algorithms. To address this issue, an adaptive variable-step homotopy-based algorithm is proposed. The algorithm adjusts the step size in the homotopy parameter during the process, handling situations where convergence fails in intermediate stages, which improve the convergence efficiency of the algorithm. A backtracking line search method is employed to constrain the iterative variables within the defined domain, enhancing the algorithm’s robustness. Homotopy algorithms are constructed based on direct iteration, Wegstein, Broyden, and Newton methods. The impact of homotopy parameters and auxiliary functions is analyzed. Taking the production of high-density polyethylene using the slurry method as a case study, simulation results indicate that the homotopy algorithm improves the convergence of the model under different operating conditions. The built-in solver of the commercial process model software Aspen Plus can only obtain feasible solutions in 21% of the test conditions, while the homotopy algorithm converges in 88% of cases.

Key words: process system, algorithm, chemical processes, simulation, homotopy method

中图分类号:

TQ 021.8

李洪瑞, 黄纯西, 洪小东, 廖祖维, 王靖岱, 阳永荣. 基于自适应变步长同伦法的循环流程收敛算法[J]. 化工学报, 2024, 75(7): 2604-2612.

Hongrui LI, Chunxi HUANG, Xiaodong HONG, Zuwei LIAO, Jingdai WANG, Yongrong YANG. An adaptive variable-step homotopy-based algorithm for process simulation with cyclic streams[J]. CIESC Journal, 2024, 75(7): 2604-2612.

图/表 7

参考文献 33

1	Pistikopoulos E N, Barbosa-Povoa A, Lee J H, et al. Process systems engineering—the generation next?[J]. Computers & Chemical Engineering, 2021, 147: 107252.
2	Sinnott R, Towler G. Chemical Engineering Design[M]. Amsterdam: Elsevier, 2020.
3	Al-Malah K I M. Aspen Plus^®: Chemical Engineering Applications[M]. Hoboken, NJ, USA: John Wiley & Sons Inc., 2016.
4	Haydary J. Chemical Process Design and Simulation: Aspen Plus and Aspen Hysys Applications[M]. Hoboken, NJ, USA:John Wiley & Sons Inc., 2018.
5	Song J Y, Jiang Z Y, Ding Y L. Analysis and evaluation of material flow in different steel production processes by gPROMS-based simulation[J]. Energy Procedia, 2019, 158: 4218-4223.
6	Georgakis C, Chin S T, Wang Z, et al. Data-driven optimization of an industrial batch polymerization process using the design of dynamic experiments methodology[J]. Industrial & Engineering Chemistry Research, 2020, 59(33): 14868-14880.
7	Burre J, Kabatnik C, Al-Khatib M, et al. Global flowsheet optimization for reductive dimethoxymethane production using data-driven thermodynamic models[J]. Computers & Chemical Engineering, 2022, 162: 107806.
8	Schweidtmann A M, Zhang D D, von Stosch M. A review and perspective on hybrid modeling methodologies[J]. Digital Chemical Engineering, 2024, 10: 100136.
9	邱彤, 戴一阳, 赵永臣. 化工过程模拟——理论与实践[M]. 北京: 化学工业出版社, 2020.
	Qiu T, Dai Y Y, Zhao Y C. Chemical Process Simulation: Theory and Practice[M]. Beijing: Chemical Industry Press, 2020.
10	Biegler L T, Hughes R R. Feasible path optimization with sequential modular simulators[J]. Computers & Chemical Engineering, 1985, 9(4): 379-394.
11	Kang Y Y, Luo Y Q, Yuan X G. Recent progress on equation-oriented optimization of complex chemical processes[J]. Chinese Journal of Chemical Engineering, 2022, 41: 162-169.
12	Jiménez-Islas H, Martínez-González G M, Navarrete-Bolaños J L, et al. Nonlinear homotopic continuation methods: a chemical engineering perspective review[J]. Industrial & Engineering Chemistry Research, 2013, 52(42): 14729-14742.
13	Bonilla J, Diehl M, Logist F, et al. A convexity-based homotopy method for nonlinear optimization in model predictive control[J]. Optimal Control Applications and Methods, 2010, 31(5): 393-414.
14	Asadi J, Jalali Farahani F. Optimization of dimethyl ether production process based on sustainability criteria using a homotopy continuation method[J]. Computers & Chemical Engineering, 2018, 115: 161-178.
15	胡晖. 同伦新算法在精馏模拟中的应用[D]. 天津: 天津大学, 2004.
	Hu H. Application of new homotopy algorithm in the field of distillation simulation[D]. Tianjin: Tianjin University, 2004.
16	Chen W, Zhu L, Chen X, et al. Sensitivity embedded homotopy-based backtracking method for chemical process simulation[C]//2013 10th IEEE International Conference on Control and Automation (ICCA). Hangzhou, China: IEEE, 2013: 1274-1277.
17	Steffen V, Silva E A. Numerical methods and initial estimates for the simulation of steady-state reactive distillation columns with an algorithm based on tearing equations methodology[J]. Thermal Science and Engineering Progress, 2018, 6: 1-13.
18	Zhu L Y, Chen Z Q, Chen X, et al. Simulation and optimization of cryogenic air separation units using a homotopy-based backtracking method[J]. Separation and Purification Technology, 2009, 67(3): 262-270.
19	祝铃钰, 陈智强, 陈曦, 等. 大规模变工况流程模拟的回溯同伦法[J]. 高校化学工程学报, 2009, 23(4): 690-695.
	Zhu L Y, Chen Z Q, Chen X, et al. Homotopy based backtracking method for large-scale process simulations with load variation[J]. Journal of Chemical Engineering of Chinese Universities, 2009, 23(4): 690-695.
20	Chen W F, Shao Z J, Zhu L Y, et al. Homotopy with second-order correction based backtracking method for chemical process simulation[J]. Industrial & Engineering Chemistry Research, 2014, 53(39): 15080-15088.
21	Jalali F, Seader J D, Khaleghi S. Global solution approaches in equilibrium and stability analysis using homotopy continuation in the complex domain[J]. Computers & Chemical Engineering, 2008, 32(10): 2333-2345.
22	Malinen I, Kangas J, Tanskanen J. A new Newton homotopy based method for the robust determination of all the stationary points of the tangent plane distance function[J]. Chemical Engineering Science, 2012, 84: 266-275.
23	Malinen I, Kangas J, Ahola J, et al. A new homotopy-based strategy for the robust determination of all the feasible solutions for CSTR systems[J]. Periodica Polytechnica Chemical Engineering, 2016, 60(1): 8-23.
24	Malinen I, Tanskanen J. Modified bounded homotopies to enable a narrow bounding zone[J]. Chemical Engineering Science, 2008, 63(13): 3419-3430.
25	Paloschi J R. Bounded homotopies to solve systems of algebraic nonlinear equation[J].An International Journal of Computer Application in Chemical Engineering, 1995, 19(12): 1243-1254.
26	Wang Y, Topputo F. A homotopy method based on theory of functional connections[EB/OL]. 2019..
27	de la Luz López-González M, Quemada-Villagómez M L, Martinez-González G M, et al. A novel predictive homotopic path tracking algorithm to solve non-linear algebraic equations[J]. The Canadian Journal of Chemical Engineering, 2023, 101(6): 3382-3408.
28	Rahimian S K, Jalali F, Seader J D, et al. A new homotopy for seeking all real roots of a nonlinear equation[J]. Computers & Chemical Engineering, 2011, 35(3): 403-411.
29	Rahimian S K, Jalali F, Seader J D, et al. A robust homotopy continuation method for seeking all real roots of unconstrained systems of nonlinear algebraic and transcendental equations[J]. Industrial & Engineering Chemistry Research, 2011, 50(15): 8892-8900.
30	Kang X Q, Cheng H Y, Tong L W, et al. Development of a bifurcation analysis approach based on gPROMS platform[J]. Chinese Journal of Chemical Engineering, 2016, 24(12): 1742-1749.
31	Fan X Q, Sun J Y, Wang J D, et al. Stability analysis of ethylene polymerization in a liquid-containing gas-solid fluidized bed reactor[J]. Industrial & Engineering Chemistry Research, 2018, 57(16): 5616-5629.
32	Fan X Q, Sun J Y, Yang Y, et al. Thermal-stability analysis of ethylene-polymerization fluidized-bed reactors under condensed-mode operation through a TPM-PBM integrated model[J]. Industrial & Engineering Chemistry Research, 2019, 58(22): 9486-9499.
33	Khare N P, Seavey K C, Liu Y A, et al. Steady-state and dynamic modeling of commercial slurry high-density polyethylene (HDPE) processes[J]. Industrial and Engineering Chemistry Research, 2002, 41(23): 5601-5618.

编辑推荐

Metrics

阅读次数

全文

296

HTML			PDF

最新录用	在线预览	正式出版	最新录用	在线预览	正式出版
3	0	12	9	0	272

来源	本网站	其他网站

次数	79	217
比例	27%	73%

摘要

217

最新录用	在线预览	正式出版

34	0	183

来源	本网站	其他网站

次数	40	177
比例	18%	82%

模块	参数		初始值
反应器1	催化剂进料流率/(kg/h)		0.428
	乙烯进料流率/(kg/h)		5558
	氢气进料流率/(kg/h)		1.159
	正己烷进料流率/(kg/h)		17263
	反应器温度/℃		87
	反应器压力/bar		8.19
反应器2	乙烯进料流率/(kg/h)		4453
	氢气进料流率/(kg/h)		0.033
	正己烷进料流率/(kg/h)		9172
	反应器温度/℃		87
	反应器压力/bar		5.17
流股	氢气流率/(kmol/h)	乙烯流率/(kmol/h)	正己烷流率/(kmol/h)
R-Vap1	100	100	100
R-Vap2	100	100	100

模块	参数		初始值
反应器1	催化剂进料流率/(kg/h)		0.428
	乙烯进料流率/(kg/h)		5558
	氢气进料流率/(kg/h)		1.159
	正己烷进料流率/(kg/h)		17263
	反应器温度/℃		87
	反应器压力/bar		8.19
反应器2	乙烯进料流率/(kg/h)		4453
	氢气进料流率/(kg/h)		0.033
	正己烷进料流率/(kg/h)		9172
	反应器温度/℃		87
	反应器压力/bar		5.17
流股	氢气流率/(kmol/h)	乙烯流率/(kmol/h)	正己烷流率/(kmol/h)
R-Vap1	100	100	100
R-Vap2	100	100	100

组分	案例A		案例B
组分	R-Vap1	R-Vap2	R-Vap1	R-Vap2
氢气/(kmol/h)	41.450	19.726	69.114	52.680
乙烯/(kmol/h)	592.344	540.521	987.703	1445.184
正己烷/(kmol/h)	199.157	321.754	332.086	860.249

组分	案例A		案例B
组分	R-Vap1	R-Vap2	R-Vap1	R-Vap2
氢气/(kmol/h)	41.450	19.726	69.114	52.680
乙烯/(kmol/h)	592.344	540.521	987.703	1445.184
正己烷/(kmol/h)	199.157	321.754	332.086	860.249

[1]	吴哲明, 张碧云, 郑仁朝. 腈水解酶立体选择性改造及其合成布瓦西坦[J]. 化工学报, 2024, 75(7): 2633-2643.
[2]	韩志敏, 李江, 陈则齐, 刘威, 徐志明. 脉动流通道内不同纵向涡发生器的颗粒污垢特性[J]. 化工学报, 2024, 75(7): 2486-2496.
[3]	马君霞, 李林涛, 熊伟丽. 基于Tri-training GPR的半监督软测量建模方法[J]. 化工学报, 2024, 75(7): 2613-2623.
[4]	方立昌, 李梓龙, 陈博, 苏政, 贾莉斯, 王智彬, 陈颖. 基于相变微胶囊悬浮液的芯片阵列冷却特性研究[J]. 化工学报, 2024, 75(7): 2455-2464.
[5]	李子扬, 郑楠, 方嘉宾, 魏进家. 再压缩S-CO₂布雷顿循环性能分析及多目标优化[J]. 化工学报, 2024, 75(6): 2143-2156.
[6]	卢飞, 鲁波娜, 许光文. 气固微型流化床反应分析仪的理想流型判据分析[J]. 化工学报, 2024, 75(6): 2201-2213.
[7]	黄斌, 丰生杰, 傅程, 张威. 液滴撞击单丝铺展特性的数值研究[J]. 化工学报, 2024, 75(6): 2233-2242.
[8]	霍宗伟, 牛亚宾, 潘艳秋. 油水膜分离中高黏度油滴行为研究和影响因素分析[J]. 化工学报, 2024, 75(6): 2262-2273.
[9]	张晗, 张淑宁, 刘珂, 邓冠龙. 基于慢特征分析与最小二乘支持向量回归集成的草酸钴合成过程粒度预报[J]. 化工学报, 2024, 75(6): 2313-2321.
[10]	赵光耀, 杨明磊, 钱锋. 基于降方差采样策略的随机重构法[J]. 化工学报, 2024, 75(5): 1939-1950.
[11]	王金山, 王世学, 朱禹. 冷却表面温差对高温质子交换膜燃料电池性能的影响[J]. 化工学报, 2024, 75(5): 2026-2035.
[12]	李怡菲, 董新宇, 王为术, 刘璐, 赵一璠. 微肋板表面干冰升华喷雾冷却传热数值模拟[J]. 化工学报, 2024, 75(5): 1830-1842.
[13]	刘帆, 张芫通, 陶成, 胡成玉, 杨小平, 魏进家. 歧管式射流微通道液冷散热性能[J]. 化工学报, 2024, 75(5): 1777-1786.
[14]	李娟, 曹耀文, 朱章钰, 石雷, 李佳. 仿生正形尾鳍结构微通道流动与传热特性数值研究及结构优化[J]. 化工学报, 2024, 75(5): 1802-1815.
[15]	黄志鸿, 周利, 柴士阳, 吉旭. 耦合加氢装置优化的多周期氢网络集成[J]. 化工学报, 2024, 75(5): 1951-1965.

基于自适应变步长同伦法的循环流程收敛算法

An adaptive variable-step homotopy-based algorithm for process simulation with cyclic streams

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 7

参考文献 33

相关文章 15

编辑推荐

Metrics

本文评价

算法	同伦辅助函数	案例A			案例B
算法	同伦辅助函数	结果	流程运行次数	计算时间/s	结果	流程运行次数	计算时间/s
直接迭代法	—	失败	—	—	失败	—	—
Wegstein法	—	收敛	34	6	失败	—	—
Broyden法	—	收敛	28	5	失败	—	—
牛顿法	—	失败	—	—	失败	—	—
同伦-直接迭代法	牛顿同伦	失败	—	—	失败	—	—
同伦-直接迭代法	定点同伦	失败	—	—	失败	—	—
同伦-Wegstein法	牛顿同伦	收敛	74	35	收敛	792	374
同伦-Wegstein法	定点同伦	收敛	89	42	收敛	841	496
同伦-Broyden法	牛顿同伦	收敛	54	26	收敛	1577	874
同伦-Broyden法	定点同伦	收敛	58	28	收敛	1366	712
同伦-牛顿法	牛顿同伦	收敛	279	143	收敛	2711	1349
同伦-牛顿法	定点同伦	收敛	341	166	失败	—	—

算法	同伦辅助函数	案例C				案例D
算法	同伦辅助函数	收敛工况个数	收敛案例百分比/%	平均流程运行次数	平均计算时间/s	收敛工况个数	收敛案例百分比/%	平均流程运行次数	平均计算时间/s
直接迭代法	—	0	0	—	—	0	0	—	—
Wegstein法	—	7	13	46	8	10	18	42	7
Broyden法	—	7	13	41	7	12	21	39	7
牛顿法	—	4	7	235	43	6	11	192	36
同伦-直接迭代法	牛顿同伦	0	0	—	—	0	0	—	—
同伦-直接迭代法	定点同伦	0	0	—	—	0	0	—	—
同伦-Wegstein法	牛顿同伦	37	66	915	424	41	73	762	369
同伦-Wegstein法	定点同伦	35	63	976	488	45	80	797	391
同伦-Broyden法	牛顿同伦	45	80	1469	746	49	88	1225	641
同伦-Broyden法	定点同伦	42	75	1539	796	47	84	1376	721
同伦-牛顿法	牛顿同伦	34	61	3022	1579	38	68	2451	1365
同伦-牛顿法	定点同伦	32	57	2844	1464	33	59	2274	1217