
Citation: | JIA Demin, WANG Xiaoyan, WANG Peilun, YU Binbin, XU Yu, CHEN Xiumin, ZHAO Jianhui. Simulation Study of Influence of Structural Parameters on Entrainment Coefficient of Ejector Based on Response Surface Method[J]. Journal of Southwest Jiaotong University, 2024, 59(5): 1167-1175. doi: 10.3969/j.issn.0258-2724.20220232 |
To analyze the significant influencing factors of the entrainment coefficient of the ejector, a two-dimensional numerical model of the compressible flow of the ejector with air as the working medium was established, and its calculation accuracy was validated by the experimental data. Meanwhile, the calculation matrix was designed by utilizing the D-optimal experimental design method. Based on the least-squares method, the response surface prediction model of the entrainment coefficient with a second-order form was constructed, and the significant parameters of the ejection coefficient and their interaction were simulated based on the constructed model. The research results show that the coincidence between the predicted and calculated values of the entrainment coefficients proves the accuracy of the response surface prediction model; the interaction between the length of the diffuser section, the mixing section length, the diameter of the mixing section, and the distance from the nozzle outlet to the inlet of the mixing section, the interaction between the mixing section diameter and mixing section length, and the interaction between the mixing section length and the diffusion angle of the diffuser section are the key factors affecting the entrainment coefficient because their
随着我国交通行业和互联网信息技术的快速发展,物流服务信息化水平和专业能力得到了大幅度提高. 物流服务行业的迅速发展,有助于提升产品的运输效率,降低企业物流成本,提升企业自身的核心竞争力. 在激烈的竞争中,物流服务需求方的目标是用最低的成本获得最满意的物流服务,而物流服务供给方要在物流需求方满意的同时追求较高的利润,供需双方的交易是典型的不交叉关系,两者的匹配决策为自主自愿的市场化运作,其合作特征符合双边匹配理论基本要求,适用双边匹配理论和匹配决策方法进行研究分析.
在物流服务供需双边匹配决策问题中,通常物流服务供需匹配的供方或需方需要对对方做出一定的偏好排序,即供给方或需求方需要具备对对方的“敏感性偏好”,在该问题的研究中,根据供需双方各自的评价指标给出完全偏好序信息是常见的研究方法. 因此,基于完全偏好序的双边匹配问题研究备受关注,樊治平等[1]在基于完全偏好序信息的一对一双边匹配问题中,以双边主体的匹配满意度最大构建了多目标优化模型,利用双边主体的最高可接受偏好序算法进行求解获得匹配结果. 孔德财等[2]针对一对一双边匹配问题,考虑匹配的稳定性、公平性和满意性,以匹配主体的满意度最大及匹配主体间的满意度绝对差异最小建立了多目标优化模型. 姜艳萍等[3]针对基于偏好序的双边匹配问题,提出了具有抗操作和抗自亏性的匹配方法. 张笛等[4]在考虑匹配特征基础上,提出了多重偏好下的一对一多阶段双边匹配方法. Dalzell等[5]认为当用于匹配的分类变量有误差时,可以同时估计回归模型和匹配记录,提出了一种贝叶斯匹配方法. Li等[6]提出了一种基于双重犹豫模糊偏好信息的一对一双边匹配方法,在稳定匹配约束下,以双方的满意度最大、差异度最小建立多目标优化模型. Yu等[7]针对人事岗位匹配问题提出了新的直觉模糊Choquet积分聚2集算子,构建了一个直观的模糊一对一双边匹配模型. 李铭洋等[8]针对基于序值偏好信息的一对多双边匹配问題构建了以匹配主体序值之和最小为目标的多目标优化模型. Zhang等[9]对多对多双边匹配的匹配因素进行了研究,构建了一种顺序联结反应模型. Chen等[10]研究了弱偏好序的多对多双边匹配问题,提出了一种新的Pareto稳定匹配算法. Klaus等[11]利用两阶段非揭示机制研究了多对多匹配问题. Zhang等[12]提出了一种基于失望理论的不完全模糊偏好序一对一双边匹配决策理论,构建了确定最佳匹配解的双目标优化模型. Jiao等[13]考虑了多对多双边匹配中的激励相容性问题. Fan等[14]考虑心理因素将预期的各代理对相对代理的偏好序数基于不确定偏好序数计算,根据失望理论得到可能的匹配结果. Chen等[15]提出了一种根据参与者的主观偏好修正一些关键目标,有效地提高了双边参与者的满意度. 林杨等[16-19]也进行了相应的研究.
既有文献对一对一、多对多双边匹配问题进行了一定的研究,主要从整体满意度角度考虑匹配主体进行匹配,而对匹配主体中个体的满意度及其均衡性较少涉及,个体满意度未得到充分满足,很难保证匹配的公平性. 在保证整体满意度前提下,考虑双边匹配主体中个体在匹配中的满意度,主体中个体的匹配满意度较均衡时才能达到公平的双边匹配. 鉴于此,本文充分考虑整体满意度和个体满意度,构建基于完全偏好序物流服务供需多对多双边匹配两阶段优化模型,并基于NSGA-Ⅱ(non-dominated sorting genetic algorithm-Ⅱ)和线性规划方法设计模型求解算法.
假设所有物流服务供需双方都要接受决策者的指令,决策者给出多种匹配方案后,物流服务供需方根据各自不同的偏好选择匹配方案. 假设在物流服务供需双方数量充足的情况下,物流服务需求方其集合为
基于完全偏好序的物流服务供需双边匹配模型中的相关参数与变量定义如下:
αij=l(rij), |
βij=g(tij), |
式中:
按照完全偏好序的物流服务供需方双边匹配问题的描述,本文规定物流服务供需双方偏好序的倒数即为物流服务供需双方对对方的满意度,则满意度
αij=1rij, | (1) |
βij=1tij. | (2) |
基于完全偏好序向量
本文构建多对多物流服务供需双边匹配两阶段优化模型,第1阶段以物流服务需求方的整体满意度最大及物流服务需求方个体满意度的方差最小构建匹配模型(同理建立物流服务供给方匹配模型),利用遗传算法进行模型求解得到物流服务供需方理想点的双边匹配Pareto解集;第2阶段构建以供需双方满意度与理想点满意度差最小的多目标优化模型,利用遗传算法进行模型求解得到双边匹配Pareto解集.
1) 针对物流服务需求方
构建第1阶段物流服务需求方优化模型1:
maxZ1=m∑i=1n∑j=1αijxij, | (3) |
minZ2=∑i(∑jαijxij−fA)2∑i∑jxij; | (4) |
s.t.
n∑j=1xij⩽pi, | (5) |
m∑i=1xij⩽qj, | (6) |
n∑j=1αijxij>hi∑k=11n−(k−1), | (7) |
m∑i=1βijxij>tj∑λ=11m−(λ−1). | (8) |
第1阶段物流服务需求方优化模型1,只考虑物流服务需求方主体A总体满意度和物流服务需求方
2) 针对物流服务供给方B,依据满意度值矩阵
构建第1阶段物流服务供给方优化模型2,模型2只考虑物流服务供给方主体B总体满意度和物流服务需求方
第2阶段模型中,考虑物流服务供需主体的整体满意度与供需双方在匹配时理想值的差值最小. 构建第2阶段优化模型3:
minZ5=Z∗1−∑i∑jαijxij, | (10) |
minZ6=Z∗3−∑i∑jβijxij; | (11) |
s.t.
n∑j=1xij⩽pi, | (12) |
m∑i=1xij⩽qj. | (13) |
优化模型3中:
模型1为双目标非线性0-1规划问题,优化目标分别为整体满意度
1) 个体编码
如图1所示,个体采用多层0-1编码,总长度为
2) 适应度函数
针对物流服务需求方,
3) 种群进化策略
一般情况下,通过交叉获得新个体的方法来保证在进化过程中群体多样性. 由于本模型中个体的行和列均有约束,个体的合法性采用惯常的交叉方法无法得到保证,基于上述考虑,设计一种能够满足模型中约束条件的行(列)交叉算子,计算步骤如下:
步骤1 随机选中两个父个体,交换其奇数行(列)生成两个子个体;
步骤2 针对每一个子个体,从0到
步骤3 对于非法个体,选择不满足步骤2中条件的列(行),定位其奇数行中第1 次出现元素1的行列值,尝试与本行(列)第1个元素0交换位置,交换后分别检查交换列(行)是否满足步骤2中条件,如果满足,则修正后的子个体合法,如果不满足,重新尝试与本行(列)其他元素0交换位置,再次检查,直到条件满足.
以某物流供需服务为例,共有4个物流服务需求方和6个物流服务供给方,4个物流服务方可匹配的供给方个数分别不超过2、3、2、2,6个物流服务供给方可匹配的需求方个数分别不超过2、2、2、2、2、1. 以如图2所示两个父个体的行交叉为例,首先将父个体1和父个体2的第1、3行交换,生成子个体1和子个体2,下一步分别检查子个体1和子个体2的合法性.
对于子个体1,第6列的元素和为2,不满足可匹配的需求方个数约束,需要进行修正,第6列第1个元素1出现在第1行,首先与第1行第2列的元素0交换,交换后分别检查第2列和第6列,第6列元素和为1,满足要求,但第2列元素和为3,不满足可匹配的需求方个数约束,再次尝试第1行第3列的元素0交换,交换后分别检查第3列和第6列,均满足可匹配的需求方个数约束,子个体1合法. 同理,对于子个体2,行交叉后也不满足可匹配的需求方个数约束,将第3列第3行的元素1和第1列第3行的元素0交换后,个体合法.
为防止早熟收敛,在交叉的基础上引入适度的变异. 变异操作采用行(列)交换变异法得到新个体,核心思想是随机选择个体,同时基于随机方法选择两行(列),将被选个体的选中行进行交换,生成一个新个体,并检查个体的合法性. 如果个体非法,从导致个体非法的行中找到第1个元素1,与被交换行中的相同位置元素交换位置,交换后分别检查交换列(行)是否满足步骤2中条件,如果满足,则修正后的子个体合法,如果不满足,继续在导致个体非法的行中查找元素1,与被交换行中的相同位置元素交换位置,再次检查,直到条件满足.
以图3所示个体为例,选中个体的2、3行,进行交换,交换后第2行满足可匹配的需求方个数约束,但第3行元素之和为3,不满足可匹配的需求方个数约束,在第3行中查找第1个元素1,其位于第3行第1列,与第2行第1列元素进行交换,交换后2、3行同时满足可匹配的需求方个数约束,个体合法.
针对第1阶段优化得到的Pareto最优解中,根据对整体满意度和个体满意度的偏好,供需双方决策者分别选择合理的整体满意度作为第2阶段模型的理想点
加权将多目标优化模型转化为单目标优化模型4:
minF=ε(Z∗1−m∑i=1m∑j=1αijxij)+(1−ε)(Z∗3−m∑i=1m∑j=1βijxij); | (15) |
s.t.
n∑j=1xij⩽pi, | (16) |
m∑i=1xij⩽qj. | (17) |
在现实的物流服务双边匹配问题中,
某地区的中介公司收到(
需求方 | B1 | B2 | B3 | B4 | B5 |
A1 | 2 | 1 | 3 | 4 | 5 |
A2 | 1 | 5 | 4 | 3 | 2 |
A3 | 4 | 2 | 5 | 1 | 3 |
A4 | 2 | 4 | 1 | 5 | 3 |
A5 | 4 | 1 | 5 | 2 | 3 |
A6 | 3 | 1 | 5 | 2 | 4 |
供给方 | A1 | A2 | A3 | A4 | A5 | A6 |
B1 | 3 | 1 | 2 | 5 | 4 | 6 |
B2 | 6 | 4 | 1 | 5 | 2 | 3 |
B3 | 4 | 3 | 2 | 5 | 6 | 1 |
B4 | 5 | 4 | 1 | 3 | 6 | 2 |
B5 | 1 | 3 | 4 | 6 | 5 | 2 |
针对表1和表2的排序信息,假设每个物流服务需求方A可选择服务供应方最大数量分别为2、3、2、4、2、1,而物流服务供应方B可向需求方提供匹配的最大数量为3、2、2、2、1. 根据求解方法,在第1阶段考虑个体满意度,分别针对需求方A和供给方B进行优化,可分别获得供需双方的整体满意度和个体满意度Pareto前沿(如图4、5所示).
假设需求方A和供给方B可接受的个体满意度阈值均为0.22,则供需双方可选择的方案如图4、5所示,供需双方整体满意度的理想点为5.98和5.23,在此理想点下,根据基于理想点的模型求解方法,可获得如图6所示的整体满意度Pareto最优前沿.
对于图6中整体满意度Pareto最优解,可获得需求方A对供给方B间的匹配关系,假设供需双方可接受的整体满意度值为(5.93,5.20),则匹配关系如图7所示.
当不考虑个体满意度时,可以构建多对多物流服务供需双边匹配整体满意度优化模型,以供需双方匹配满意度最大为目标函数的多目标优化模型,求解中先求出物流服务供需双方不考虑对方满意度的情况下自身的理想点满意度,再通过物流服务供需双方实际最大满意度与理想点满意度的差距最小建立多目标优化模型,利用线性规划方法进行模型求解得到双边匹配Pareto解集. 基于NSGA-Ⅱ的多目标优化算法获得整体满意度和个体满意度的Pareto解,为决策者提供了体现整体和个体满意度不同偏好的决策方案,指导决策者获得供需双方满意的匹配关系.
根据基于理想点法对算例求解,其理想点为(7.33,6.75),可获得如图8所示的整体满意度Pareto最优前沿.
从图8中可以看出:需求方A对供给方B的整体满意度和B对需求方A的整体满意度是此消彼长的关系,即需求方A对供给方B取最大整体满意度时,供给方B对需求方A的整体满意度最低,因此,在决策过程中,需要供需双方谈判以获得彼此可接受的满意度.
对于图8中任意解,可获得需求方A对供给方B间的匹配关系,如供需双方可接受的满意度值为(5.10,6.75),则匹配关系如图9所示. 在此匹配关系下,需求方A和供给方B的个体满意度方差分别0.57和0.22,个体满意度越小,表示个体之间的差异越小,对于本方案,需求方A的个体满意度差异较大,选择此方案导致需求方之间满意度不公平,容易导致矛盾产生.
由图6与图8的Pareto前沿对比可以看出:考虑个体满意度差异的情况下得到的整体满意度明显劣于未考虑个体满意度的情况,相对未考虑个体满意度情况,引入个体满意度时会导致双方整体满意度下降,即为使每个个体能够得到满意、公平的匹配结果,需要牺牲双方的整体的满意度. 表明考虑个体满意度最优时,整体满意度会受到一定的影响,但其结果更贴近现实情况,并能得到公平的匹配结果,因为个体满意度离理想点越近,匹配结果越接近个体对匹配结果的期望.
本文中同时考虑整体满意度和个体满意度均衡性的供需双方多对多双边匹配两阶段优化模型,能够刻画物流服务中供需双方的利益关系,既能考虑物流服务供需方整体满意度,又兼顾供需方匹配个体满意度的均衡性. 后续研究将考虑物流服务供需方个体差异化和优先级的稳定匹配问题,以提高实用性和适应性.
[1] |
FU W N, LIU Z L, LI Y X, et al. Numerical study for the influences of primary steam nozzle distance and mixing chamber throat diameter on steam ejector performance[J]. International Journal of Thermal Sciences, 2018, 132: 509-516. doi: 10.1016/j.ijthermalsci.2018.06.033
|
[2] |
SAMAKÉ O, GALANIS N, SORIN M. On the design and corresponding performance of steam jet ejectors[J]. Desalination, 2016, 381: 15-25. doi: 10.1016/j.desal.2015.11.027
|
[3] |
CHEN W X, CHONG D T, YAN J J, et al. The numerical analysis of the effect of geometrical factors on natural gas ejector performance[J]. Applied Thermal Engineering, 2013, 59(1/2): 21-29.
|
[4] |
BANASIAK K, PALACZ M, HAFNER A, et al. A CFD-based investigation of the energy performance of two-phase R744 ejectors to recover the expansion work in refrigeration systems: an irreversibility analysis[J]. International Journal of Refrigeration, 2014, 40: 328-337. doi: 10.1016/j.ijrefrig.2013.12.002
|
[5] |
HAKKAKI-FARD A, AIDOUN Z, OUZZANE M. A computational methodology for ejector design and performance maximisation[J]. Energy Conversion and Management, 2015, 105: 1291-1302. doi: 10.1016/j.enconman.2015.08.070
|
[6] |
ZHU Y H, CAI W J, WEN C Y, et al. Numerical investigation of geometry parameters for design of high performance ejectors[J]. Applied Thermal Engineering, 2009, 29(5/6): 898-905.
|
[7] |
CHONG D T, YAN J J, WU G S, et al. Structural optimization and experimental investigation of supersonic ejectors for boosting low pressure natural gas[J]. Applied Thermal Engineering, 2009, 29(14/15): 2799-2807.
|
[8] |
NIKIFOROW K, KOSKI P, KARIMÄKI H, et al. Designing a hydrogen gas ejector for 5 kW stationary PEMFC system-CFD-modeling and experimental validation[J]. International Journal of Hydrogen Energy, 2016, 41(33): 14952-14970. doi: 10.1016/j.ijhydene.2016.06.122
|
[9] |
王子瑞. 基于响应面法的贯流式水轮机多目标优化设计[D]. 西安: 西安理工大学, 2012.
|
[10] |
OMIDVAR A, GHAZIKHANI M, MODARRES RAZAVI S M R. Entropy analysis of a solar-driven variable geometry ejector using computational fluid dynamics[J]. Energy Conversion and Management, 2016, 119: 435-443. doi: 10.1016/j.enconman.2016.03.090
|
[11] |
LI S Y, YAN J, LIU Z, et al. Optimization on crucial ejector geometries in a multi-evaporator refrigeration system for tropical region refrigerated trucks[J]. Energy, 2019, 189: 116347.1-116347.14.
|
[12] |
BANASIAK K, HAFNER A, ANDRESEN T. Experimental and numerical investigation of the influence of the two-phase ejector geometry on the performance of the R744 heat pump[J]. International Journal of Refrigeration, 2012, 35(6): 1617-1625. doi: 10.1016/j.ijrefrig.2012.04.012
|
[13] |
SOLMAZ H, ARDEBILI S M S, CALAM A, et al. Prediction of performance and exhaust emissions of a CI engine fueled with multi-wall carbon nanotube doped biodiesel-diesel blends using response surface method[J]. Energy, 2021, 227: 1205181.1-1205181.13.
|
[1] | WANG Xiaoyan, KANG Tianqin, WANG Peilun, JIA Demin, XU Yu, ZHAO Jianhui. Simulation Research on Ejector for Natural Gas Engine with High Exhaust Gas Recirculation Rate[J]. Journal of Southwest Jiaotong University, 2024, 59(2): 467-476. doi: 10.3969/j.issn.0258-2724.20210950 |
[2] | CAI Huikun, SU Lijun, LIAO Yidai, WENG Zeju, XU Chen. Heat Dissipation Characteristics Analysis of Surface Air-Oil Heat Exchanger Applied in Aircraft Engine[J]. Journal of Southwest Jiaotong University, 2021, 56(1): 214-220. doi: 10.3969/j.issn.0258-2724.20180964 |
[3] | ZHANG Yafeng, ZHANG Shaohua, ZHOU Gang, ZHANG Jiyang, QING Tao, ZHOU Ningning. Strain Rate Sensitivity of Porous Cotton-Phenolic Bearing Retainer[J]. Journal of Southwest Jiaotong University, 2021, 56(4): 847-852, 863. doi: 10.3969/j.issn.0258-2724.20190457 |
[4] | TIAN Guanghui, SHEN Mingrong, YANG Ruifang, ZHANG Qingzhao, WANG Zhen. Shear Rheological Characteristics and Nonlinear Constitutive Model of Serrate Structure Surface[J]. Journal of Southwest Jiaotong University, 2018, 53(5): 1017-1025. doi: 10.3969/j.issn.0258-2724.2018.05.019 |
[5] | DENG Shasha, SHEN Huoming, LIU Lang, TANG Huaiping. Research of Calculation Method for Wave Forces Acting on Large-Scale Bridge Piers Based on Diffraction Theory[J]. Journal of Southwest Jiaotong University, 2018, 53(2): 266-271. doi: 10.3969/j.issn.0258-2724.2018.02.007 |
[6] | ZHANG Xianmin, DONG Qian, LÜ, Yaozhi. Mechanical Responses of Pavement under Aircrafts with Different Main Landing Gears[J]. Journal of Southwest Jiaotong University, 2014, 27(4): 675-681. doi: 10.3969/j.issn.0258-2724.2014.04.018 |
[7] | HUO Zhiqin, RU Yi, HAN Songchen. Analysis Model of Transport Aircraft Veering off Runway during Landing Phase[J]. Journal of Southwest Jiaotong University, 2012, 25(5): 895-900. doi: 10.3969/j.issn.0258-2724.2012.05.026 |
[8] | WANG Zhen, LI Guhua, ZHANG Zejiang, HUANG Tao. Effect of Fire Fighting Water on Compressive Strength of Concrete in Tunnel Fire[J]. Journal of Southwest Jiaotong University, 2011, 24(2): 223-227. doi: 10.3969/j.issn.0258-2724.2011.02.008 |
[9] | TAN Daiming, QI Taiyue. Sensitivity Analysis of Geoelectric Parameters in Transient Electromagnetic Method for Tunnels[J]. Journal of Southwest Jiaotong University, 2008, 21(5): 679-684. |
[10] | REN Qiliang, PENG Qiyuan. Sensitivity Analysis of Traffic Impact of Large Buildings on Traffic of Road Network[J]. Journal of Southwest Jiaotong University, 2006, 19(5): 663-668. |
[11] | WANG Feng-yan, XU Jian-ping, WU Song-rong, XU Jun-feng. Improvement ofTransientResponse Performances of M icroprocessor Voltage Regulator[J]. Journal of Southwest Jiaotong University, 2005, 18(2): 153-157. |
[12] | JIANG Gao, SONGXue-jun, CHENDe-xin. Comprehensive EvaluationM odel for Difficulty Coefficient ofRoad Adm inistration[J]. Journal of Southwest Jiaotong University, 2005, 18(3): 412-416. |
[13] | LIUXiu-feng, ZHANG Ai-li. H igh-DimensionalAffine Codes[J]. Journal of Southwest Jiaotong University, 2005, 18(2): 179-183. |
[14] | FAN Li-li, WUNa, MA Yue. Parameter Analyses in the Model of Profit Sharing and Valuation of Technical Assets[J]. Journal of Southwest Jiaotong University, 2002, 15(6): 693-697. |
[15] | ZHAO Can-hui. The Asynchronous Excitation Model for the Seismic Response Analysis of Long-Span Bridges[J]. Journal of Southwest Jiaotong University, 2002, 15(3): 236-240. |
[16] | LIUHai-yan, SUNHong, DUWen. Sensitivity Analysis of Equilibrium Network Flow[J]. Journal of Southwest Jiaotong University, 2002, 15(4): 439-442. |
[17] | ZhangXiaoxia, Pan Wei, Luo Bin. Output Characteristics of Vertical Cavity Surface Emitting Lasers[J]. Journal of Southwest Jiaotong University, 1999, 12(5): 480-484. |
[18] | ZHANG Xiao-Xia- , BO Wei, Luo- Bin. Output Characteristics of Vertical Cavity Surface Emitting Lasers[J]. Journal of Southwest Jiaotong University, 1999, 12(5): 480-484. |
需求方 | B1 | B2 | B3 | B4 | B5 |
A1 | 2 | 1 | 3 | 4 | 5 |
A2 | 1 | 5 | 4 | 3 | 2 |
A3 | 4 | 2 | 5 | 1 | 3 |
A4 | 2 | 4 | 1 | 5 | 3 |
A5 | 4 | 1 | 5 | 2 | 3 |
A6 | 3 | 1 | 5 | 2 | 4 |
供给方 | A1 | A2 | A3 | A4 | A5 | A6 |
B1 | 3 | 1 | 2 | 5 | 4 | 6 |
B2 | 6 | 4 | 1 | 5 | 2 | 3 |
B3 | 4 | 3 | 2 | 5 | 6 | 1 |
B4 | 5 | 4 | 1 | 3 | 6 | 2 |
B5 | 1 | 3 | 4 | 6 | 5 | 2 |