CN116415449A - 顺酐反应器的设计方法、设计系统及信息数据处理终端 - Google Patents
顺酐反应器的设计方法、设计系统及信息数据处理终端 Download PDFInfo
- Publication number
- CN116415449A CN116415449A CN202310686938.2A CN202310686938A CN116415449A CN 116415449 A CN116415449 A CN 116415449A CN 202310686938 A CN202310686938 A CN 202310686938A CN 116415449 A CN116415449 A CN 116415449A
- Authority
- CN
- China
- Prior art keywords
- flow
- maleic anhydride
- temperature
- anhydride reactor
- reactor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- FPYJFEHAWHCUMM-UHFFFAOYSA-N maleic anhydride Chemical compound O=C1OC(=O)C=C1 FPYJFEHAWHCUMM-UHFFFAOYSA-N 0.000 title claims abstract description 140
- 238000000034 method Methods 0.000 title claims abstract description 90
- 238000013461 design Methods 0.000 title claims abstract description 89
- 238000012545 processing Methods 0.000 title claims abstract description 8
- 238000009826 distribution Methods 0.000 claims abstract description 116
- 150000003839 salts Chemical class 0.000 claims abstract description 45
- 238000004458 analytical method Methods 0.000 claims abstract description 26
- 238000004088 simulation Methods 0.000 claims abstract description 23
- 230000008878 coupling Effects 0.000 claims abstract description 20
- 238000010168 coupling process Methods 0.000 claims abstract description 20
- 238000005859 coupling reaction Methods 0.000 claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims description 93
- 239000012530 fluid Substances 0.000 claims description 49
- 230000000149 penetrating effect Effects 0.000 claims description 36
- 230000035945 sensitivity Effects 0.000 claims description 25
- 238000006243 chemical reaction Methods 0.000 claims description 18
- 238000012544 monitoring process Methods 0.000 claims description 12
- 238000005457 optimization Methods 0.000 claims description 12
- 238000010206 sensitivity analysis Methods 0.000 claims description 12
- 238000012360 testing method Methods 0.000 claims description 11
- 238000010521 absorption reaction Methods 0.000 claims description 8
- 238000012546 transfer Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 230000001105 regulatory effect Effects 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- 230000005484 gravity Effects 0.000 claims description 5
- 238000005520 cutting process Methods 0.000 claims description 3
- 239000003054 catalyst Substances 0.000 description 7
- UHOVQNZJYSORNB-UHFFFAOYSA-N Benzene Chemical compound C1=CC=CC=C1 UHOVQNZJYSORNB-UHFFFAOYSA-N 0.000 description 6
- 230000003647 oxidation Effects 0.000 description 6
- 238000007254 oxidation reaction Methods 0.000 description 6
- 238000003860 storage Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- WFDIJRYMOXRFFG-UHFFFAOYSA-N Acetic anhydride Chemical compound CC(=O)OC(C)=O WFDIJRYMOXRFFG-UHFFFAOYSA-N 0.000 description 3
- LGRFSURHDFAFJT-UHFFFAOYSA-N Phthalic anhydride Natural products C1=CC=C2C(=O)OC(=O)C2=C1 LGRFSURHDFAFJT-UHFFFAOYSA-N 0.000 description 3
- JHIWVOJDXOSYLW-UHFFFAOYSA-N butyl 2,2-difluorocyclopropane-1-carboxylate Chemical compound CCCCOC(=O)C1CC1(F)F JHIWVOJDXOSYLW-UHFFFAOYSA-N 0.000 description 3
- 239000006227 byproduct Substances 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 150000001336 alkenes Chemical class 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 125000004432 carbon atom Chemical group C* 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- IJDNQMDRQITEOD-UHFFFAOYSA-N n-butane Chemical compound CCCC IJDNQMDRQITEOD-UHFFFAOYSA-N 0.000 description 2
- JRZJOMJEPLMPRA-UHFFFAOYSA-N olefin Natural products CCCCCCCC=C JRZJOMJEPLMPRA-UHFFFAOYSA-N 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 239000002994 raw material Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 150000008064 anhydrides Chemical class 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 239000002826 coolant Substances 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 239000002778 food additive Substances 0.000 description 1
- 235000013373 food additive Nutrition 0.000 description 1
- 235000003086 food stabiliser Nutrition 0.000 description 1
- 230000020169 heat generation Effects 0.000 description 1
- 231100000086 high toxicity Toxicity 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000010687 lubricating oil Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 239000003973 paint Substances 0.000 description 1
- 239000000575 pesticide Substances 0.000 description 1
- 239000004014 plasticizer Substances 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000003381 stabilizer Substances 0.000 description 1
- 229920003002 synthetic resin Polymers 0.000 description 1
- 239000000057 synthetic resin Substances 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Physical Or Chemical Processes And Apparatus (AREA)
Abstract
本发明公开了一种顺酐反应器的设计方法、设计系统及信息数据处理终端,包括:S1、计算顺酐反应器模型的流场和温度场;S2、优化反应器内部流道;S3、设计换热介质的出口和入口的流量分配板;S4、在熔盐泵正常流量的60%~100%范围内,更改熔盐流量,重复计算得到若干个温度场;根据得到的各个温度场,选取温度最高点截面、径向温差最高点截面设置测温点;S5、完成顺酐反应器流道设计。本发明通过对顺酐反应器的流场‑温度场耦合仿真分析,获得顺酐反应器内的流场、温度场分布,并根据这些流场和温度场信息,对顺酐反应器的结构进行优化,确保顺酐反应器的管束区流场和温度场的均匀性。
Description
技术领域
本发明属于顺酐反应器技术领域,具体涉及一种顺酐反应器的设计方法、设计系统及信息数据处理终端。
背景技术
顺丁烯二酸酐(简称顺酐,MA)又名马来酸酐,是目前世界上产量仅次于苯酐和醋酐的第三大有机酸酐,被广泛应用于合成树脂、涂料、农药、医药、增塑剂、润滑油添加剂、食品添加剂和稳定剂等方面。顺酐的生产工艺路线按照原料的不同可分为正丁烷氧化法、苯氧化法、C4烯烃氧化法和苯酐副产法。其中苯酐副产法顺酐的产量有限,C4烯烃氧化法因副产物较多已被淘汰,苯氧化法因其毒性大、环境污染和碳原子利用率低等特点被逐步淘汰。正丁烷氧化法因其原料价廉、污染相对较轻和碳原子利用率高等优点,已在顺酐生产中居于主导地位。
众所周知,顺酐反应器是实现顺酐生产工艺的重要设备之一;针对不同结构的顺酐反应器,其内部的流速分布和温度分布不尽相同,即使是同一顺酐反应器,其内部不同位置的流速分布和温度分布也不尽相同,而流速分布和温度分布的差异化直接影响反应器催化剂的使用寿命,进而影响整个装置的运行周期和运行成本,为此,设计开发一种顺酐反应器的设计方法及设计系统具有重要的意义。
发明内容
本发明为解决公知技术中存在的技术问题,提供一种顺酐反应器的设计方法、设计系统及信息数据处理终端,运用计算流体力学的方法,通过调整换热介质的出口/入口的流量分配板的结构;对顺酐反应器内部流道进行优化设计,包括调整折流板边上穿流孔的大小与数量、调整折流板板间距、调整中心管直径,实现换热介质在反应器内部均匀流动;进行流体流场和温度场的耦合,获得反应器内部温度场数值,选定顺酐反应器空间内熔盐流动性最差的点进行熔盐温度监测。控制反应器任意横向截面上温度监测点的径向温差最大值不大于±2℃;轴向温差最大值不大于5℃;进而达到避免催化剂出现超温现象,保证顺酐反应器安全高效生产的目的。
本发明通过对顺酐反应器的流场-温度场耦合仿真分析,获得顺酐反应器内的流场分布,特别是关键位置(进出口流量分配板、折流板、管束区)的流速分布和管束区的温度场分布,从而根据这些流场和温度场信息,对顺酐反应器的流道结构进行优化,确保顺酐反应器的管束区流场和温度场的均匀性。
本发明的第一目的是提供一种顺酐反应器的设计方法,包括:
S1、计算顺酐反应器模型的流场和温度场;具体包括:
S101、根据顺酐反应器结构,建立计算流体力学模型;
S102、设置计算流体力学模型的参数:进口熔盐的质量流量和温度、管束区的阻力、松弛因子;
S103、根据反应管热负荷试验数据,设置温度场计算的前置关联条件;
S104、对流场和温度场进行耦合分析计算:根据顺酐反应器入口的质量流量、熔盐的密度、运动粘度和管束区的流通面积、传热管的直径,计算管束区的雷诺数;根据雷诺数选择流动模型,若雷诺数小于104选择层流模型,若雷诺数不小于104选择湍流模型,通过模拟计算得到速度场分布;根据温度场计算的前置关联条件,得到温度场的分布;
S105、基于边界条件完成流场和温度场的计算,并获得初始管束区轴向和径向的温度分布,通过选取不同的横向界面和竖向界面,获得轴向温差最大值和径向温差最大值,并判断径向温差最大值和轴向温差最大值是否满足设计要求,所述设计要求为:径向温差最大值是否不大于±2℃,轴向温差最大值是否不大于5℃;当初始设计结构不满足上述设计要求时,执行S2,当初始设计结构满足上述设计要求时,执行S3;
S2、优化反应器内部流道;具体包括:
首先对流道参数进行敏感性分析,流道参数包括折流板上穿流孔的开孔范围、折流板上穿流孔尺寸、折流板内径、折流板外径、折流板间距和中心管尺寸;
然后对流道参数敏感性进行大小排序;
最后根据敏感性从大到小的排序,调整流道参数;调整完成后返回至S101;
S3、设计换热介质的出口和入口的流量分配板;具体为:
通过计算流体力学的方法,确定外循环错流模型下不等面积的流量分配板尺寸的分布规律;在选取开孔直径为定值的前提下,获得不等面积流量分配板上的开孔率,对进入反应器的流量进行调节;最终得到开孔率S的规律为:
S=0.0024x2+0.0156x+0.3161;
其中,S为开孔率,x为流量分配板中心线和入口管中心线之间的夹角或流量分配板中心线和出口管中心线之间的夹角;
S4、在熔盐泵正常流量的60%~100%范围内,更改熔盐流量,重复计算得到若干个温度场;
根据得到的各个温度场,选取温度最高点截面、径向温差最高点截面设置测温点;
S5、完成顺酐反应器流道设计;具体包括:
将折流板上穿流孔的开孔范围、折流板上穿流孔尺寸、折流板内径、折流板外径、折流板间距、中心管尺寸、出口和入口的流量分配板开孔率的数值作为顺酐反应器施工设计的数据。
优选地,所述顺酐反应器包括圆柱形的空腔壳体,在空腔壳体的内腔设置有环形结构的折流板和中心管,在空腔壳体的侧壁开设有左右对称的入口管和出口管;S101具体为:
S1011、建立顺酐反应器的三维流场计算模型;
S1012、通过布尔运算,提取三维流场计算模型中的流体域模型;
S1013、将流体域模型沿X平面和Y平面进行切割,得到1/4的流体域模型作为计算流体力学分析模型;所述空腔壳体的中心轴位于X平面和Y平面内,且X平面和Y平面相互垂直。
优选地,S102具体为:
基于顺酐反应器的运行工况,设置进口熔盐的质量流量和温度;
根据顺酐反应器内部的管束排布方式,设置管束区的阻力;
S103具体为:根据试验得到的顺酐反应器管侧的化学反应特性,设置管束区的吸热功率、放热功率。
优选地,S104中,壁面采用无滑移边界,壁面函数采用two-layer all y+treatment;动量方程添加重力源项;边界条件设定为质量流量入口、压力出口;将动量方程中速度参数的松弛因子调节至0.15;将动量方程中压力参数的松弛因子调节至0.03;将湍动能参数的松弛因子调节至0.15;将湍流耗散率参数的松弛因子调节至0.2;将能量参数的松弛因子调节至0.3;计算过程中监测质量、速度、湍动能和湍流耗散率的残差、进出口的压降和温差、计算域的最大速度和最大温度,当监测量稳定且波动幅度不大于1%时,认为计算收敛。
本发明的第二目的是提供一种顺酐反应器的设计系统,包括:
计算模块:计算顺酐反应器模型的流场和温度场;具体包括:
S101、获取计算流体力学分析模型;
S102、设置下列参数:进口熔盐的质量流量和温度、管束区的阻力;
S103、根据反应管热负荷试验数据,设置温度场计算的前置关联条件;
S104、对流场和温度场进行耦合分析计算:根据顺酐反应器入口的质量流量、熔盐的密度、运动粘度和管束区的流通面积、传热管的直径,计算管束区的雷诺数;根据雷诺数选择流动模型,若雷诺数小于104选择层流模型,若雷诺数不小于104选择湍流模型,通过模拟计算得到速度场分布;根据温度场计算的前置关联条件,得到温度场的分布;
S105、基于边界条件完成流场和温度场的计算,并获得初始管束区轴向和径向的温度分布,通过选取不同的横向界面和竖向界面,获得轴向温差最大值和径向温差最大值,并判断径向温差最大值和轴向温差最大值是否满足设计要求,所述设计要求为:径向温差最大值是否不大于±2℃,轴向温差最大值是否不大于5℃;当初始设计结构不满足上述设计要求时,执行优化模块,当初始设计结构满足上述设计要求时,执行进出口设计模块和测温点选取模块;
优化模块:优化反应器内部流道;具体包括:
首先将折流板上穿流孔的开孔范围、折流板上穿流孔的面积、折流板内径、折流板外径、折流板间距和中心管尺寸逐一作为敏感性因素进行敏感性分析,然后根据分析结果,优化每个敏感性因素的大小;
进出口设计模块:设计换热介质的出口和入口的流量分配板开孔率;具体为:
通过计算流体力学的方法,确定外循环错流模型下不等面积的流量分配板尺寸的分布规律;在选取开孔直径为定值的前提下,获得不等面积流量分配板上的开孔率,对进入反应器的流量进行调节;最终得到开孔率S的规律为:
S=0.0024x2+0.0156x+0.3161;
其中,S为开孔率,x为流量分配板中心线和入口管中心线之间的夹角或流量分配板中心线和出口管中心线之间的夹角;
测温点选取模块:在熔盐泵正常流量的60%~100%范围内,更改S102中的熔盐流量,重复计算得到若干个温度场;
根据得到的各个温度场,选取温度最高点截面、径向温差最高点截面设置测温点;
流道设计模块:完成顺酐反应器流道设计;具体包括:
将折流板上穿流孔的开孔范围、折流板上穿流孔尺寸、折流板内径、折流板外径、折流板间距、中心管尺寸、出口和入口的流量分配板开孔率的数值作为顺酐反应器施工设计的数据。
本发明的第三目的是提供一种信息数据处理终端,用于实现上述顺酐反应器的设计方法。
本发明具有的优点和积极效果是:
本发明采用计算流体力学方法,建立顺酐反应器的三维流场计算模型,采用多孔介质建模的方法进行流场分析,并进行流场和温度场的耦合计算,通过调节动量方程和能量方程的松弛因子,确保计算的收敛性,获得可视化的顺酐反应器流场和温度场分布,并对顺酐反应器结构进行迭代优化,提高顺酐反应器管束区的流场和温度场的均匀性;
本发明能够获得可视化的顺酐反应器的流场和温度场分布,并对顺酐反应器结构进行迭代优化,通过分析顺酐反应器三维流场计算模型的计算结果,获取顺酐反应器关键位置的流场信息(如进出口流量分配板、折流板),以提高顺酐反应器管束区流场和温度场均匀性为优化目标,对进出口流量分配板的开孔尺寸和数量、折流板的开孔尺寸和数量进行迭代,最终获得管束区流场和温度场分布更加均匀的顺酐反应器几何结构;
本发明中的realizable k-e湍流模型采用雷诺平均方法模拟湍流,而realizablek-e湍流模型属于雷诺平均方法的其中一类。将流体的瞬态运动分解为平均运动和脉动运动,其中脉动项采用雷诺应力体现,再根据涡粘模型假设,引入湍动能和湍流耗散率模拟雷诺应力对流动的影响;
本发明中的two-layer all y+ treatment壁面函数在受到流体粘性影响较大的壁面附近,采用半经验公式计算壁面与充分发展湍流区域之间的粘性影响区域;
通过采用上述技术方案,特别是动量方程中松弛因子的调节方法和流场-温度场耦合分析方法,确保顺酐反应器流场仿真计算的收敛性,并获得顺酐器内可视化的流场分布和管束区的温度场分布,随后根据计算获得的流场分布和温度场分布对顺酐反应器的结构进行优化,使得顺酐反应器管束区的流场和温度场更加均匀,提高顺酐反应器的性能。
附图说明
图1为顺酐反应器的结构图;
图2为本申请优选实施例中构建的计算模型;
图3为本申请优选实施例中网格模型的示意图;
图4为本申请优选实施例中顺酐反应器的流线图;
图5为本申请优选实施例中出口截面的速度云图;
图6为本申请优选实施例中流量分配板处速度分布云图;
图7为本申请优选实施例中反应管催化剂填装图;
图8为本申请优选实施例的流程图;
图9为本申请优选实施例中吸热-放热曲线。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效,兹列举以下实施例,并配合附图详细说明如下:
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的技术方案,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
Velocity:Magnitude(m/s)为速度大小;
请参阅图1至图9。
一种顺酐反应器的设计方法,如图1所示:顺酐反应器包括圆柱形的空腔壳体,在空腔壳体的内腔设置有环形结构的折流板和中心管,在空腔壳体的侧壁开设有左右对称的入口管和出口管;上述顺酐反应器的设计方法包括如下步骤:
S1、计算顺酐反应器模型的流场和温度场;具体包括:
S101、获取计算流体力学分析模型;具体为:
首先,根据顺酐反应器结构,建立计算流体力学分析模型;采用三维建模软件Autodesk Inventor构建顺酐反应器的三维流场计算模型,顺酐反应器的三维流场计算模型与实体模型结构和尺寸一致,通过布尔运算,取出三维流场计算模型中的流场区域,构建流体域模型;
然后,由于顺酐反应器总体呈现沿中心圆筒的X平面和Y平面的对称结构,所述空腔壳体的中心轴位于X平面和Y平面内,且X平面和Y平面相互垂直,因此将流体域模型沿X平面和Y平面进行切割,得到其的1/4的流体域模,即顺酐反应器的计算流体力学分析模型,这样在不影响计算精度的情况下可以减小计算规模,从而提高计算效率。构建的计算模型如图2所示;
最后,将计算模型导出为.stp格式用于网格划分;具体为:
对流体区域进行网格划分;在对计算流体力学分析模型的流场进行分析之前,需要对计算区域在空间上进行网格划分,网格划分采用STAR-CCM+软件进行,采用非结构化的多面体网格。按照流通区域可以将计算流体力学分析模型划分为:进口接管区、进口流量分配板区、管束区、折流板区、出口流量分配板区和出口管区。为保证获得更准确的顺酐反应器进出口位置和折流板处流动工质的流速分布,在进口流量分配板区、折流板区和出口流量分配板区对网格进行加密。网格模型如图3所示;
仿真计算:在完成对计算流体力学分析模型的网格划分后,需要对其进行流场仿真,计算采用STAR-CCM+软件进行,计算为三维稳态不可压缩流动,工质的密度、运动粘度、定压比热容和导热系数均为随温度变化的函数;
S102、设置计算流体力学模型的参数:进口熔盐的质量流量和温度、管束区的阻力、松弛因子;具体为:
基于顺酐反应器的运行工况,设置进口熔盐的质量流量和温度;
根据顺酐反应器特有的内部管束排布方式,设置管束区的阻力;设置管束区阻力的方法为:根据fluent软件设置管束区多孔介质阻力本构关系,通过同类型设备现场压降数据对压降计算结果进行修正。具体的,根据1万吨、2万吨、2.5万吨产能相同布管规格的压降实测数据,推导出布管区最小单位管束组合的压降数值的实际数值推导公式,将该公式用于修正流体力学模拟计算中管束区阻力的矫正,得到符合实际工况的管束区的阻力,确保模拟计算关键参数的正确性;
根据顺酐反应器管侧的化学反应特性设置管束区的吸热功率和放热功率;吸热功率和放热功率设置方法如下:采用流场和温度场耦合的计算方法,管束区的热源设置采用分布的体积热源。由于顺酐反应器管侧的入口温度为150℃左右,而壳侧的熔盐温度通常超过400℃,因此在顺酐反应器中心轴向1/3以下位置,管侧气体吸收壳侧熔盐和自身化学反应的热量而升温,导致壳侧熔盐温度降低。在顺酐反应器中心轴向1/3以上位置,管侧气体在催化剂的作用下发生放热的化学反应,壳侧熔盐吸收该反应热而升温。壳侧熔盐沿轴向的吸热和放热的体积热源分布根据顺酐催化剂单管试验,将图7各分段测量得到的温度数据进行拟合,获得特有的吸热-放热曲线。在图7中:A表示反应管外径,B表示反应管内径,C表示反应管长度,D表示管板厚度,E表示支撑弹簧高度,F表示顶部空间,G表示底部惰性段高度,H表示催化剂装填高度,I表示顶部惰性段高度,J表示底部空间。
S103、根据反应管热负荷试验数据,设置温度场计算的前置关联条件(管束区的吸热功率、放热功率),将温度场与流场进行耦合关联,根据流场的计算结果,模拟出温度场的分布;具体为:
流热耦合计算,主要解决流动与结构的换热计算,例如换热器的温度场计算,通过该方法可以获得较为准确的结构表面的对流换热边界条件。使用Fluent模块完成流热耦合计算,也可以使用Fluent+结构热模块和系统耦合器完成计算,基于ANSYS Workbench平台可以完成单向和双向流热耦合技术。通过将温度场与流场的耦合关联,根据流场的计算结果,模拟出温度场的分布,为设置催化剂测温点提供依据。
S104、对流场和温度场进行耦合分析计算:根据顺酐反应器入口的质量流量、熔盐的密度、运动粘度和管束区的流通面积、传热管的直径,计算管束区的雷诺数;根据雷诺数选择流动模型,若雷诺数小于104选择层流模型,若雷诺数不小于104选择湍流模型,通过模拟计算得到速度场分布;根据温度场计算的前置关联条件,得到温度场的分布;具体为:
首先对流场进行初步分析,根据顺酐反应器入口的质量流量、熔盐的密度、运动粘度和管束区的流通面积、传热管的直径等参数,计算管束区的雷诺数,计算结果表明,当雷诺数大于104,属于湍流,此时的流动模型采用realizable k-e湍流模型;当雷诺数不大于104,属于层流,采用层流模型;壁面采用无滑移边界,壁面函数采用 two-layer all y+treatment方法。动量方程添加重力源项,值为9.81 m/s2。边界条件设定为质量流量入口,压力出口。为了确保计算的收敛性,针对动量方程中的速度参数,其松弛因子调节至0.15;针对动量方程中的压力参数,其松弛因子调节至0.03;针对湍动能参数,其松弛因子调节至0.15;针对湍流耗散率参数,其松弛因子调节至0.2;针对能量参数,其松弛因子调节至0.3。计算过程中监测质量、速度、湍动能和湍流耗散率的残差、进出口的压降和温差、计算域的最大速度和最大温度,当上述监测量稳定且波动幅度不大于1%时,认为计算收敛。
S105、基于边界条件完成流场和温度场的计算,并获得初始管束区轴向和径向的温度分布,通过选取不同的横向界面和竖向界面,获得轴向温差最大值和径向温差最大值,并判断径向温差最大值和轴向温差最大值是否满足设计要求,所述设计要求为:径向温差最大值是否不大于±2℃,轴向温差最大值是否不大于5℃;当初始设计结构不满足上述设计要求时,执行S2,当初始设计结构满足上述设计要求时,执行S3;
S2、优化反应器内部流道;具体包括:
首先对流道参数进行敏感性分析,流道参数包括折流板上穿流孔的开孔范围、折流板上穿流孔尺寸、折流板内径、折流板外径、折流板间距和中心管尺寸;
然后对流道参数敏感性进行大小排序;
最后根据敏感性从大到小的排序,调整流道参数;调整完成后返回至S101;
S2.1折流板上穿流孔的开孔范围和面积;
开孔区域的确定:根据初始顺酐反应器的流速矢量图和温度场可知:在折流板内侧边缘和外侧边缘会产生漩涡,进而在这些区域出现熔盐流动的滞死区,最终导致管束区局部换热不均匀,导致局部温度过热。因此需要在折流板的内侧和外侧设计穿流孔使得部分熔盐从穿流孔流出,避免出现漩涡而造成局部换热不均。通过对穿流孔的直径和数量的敏感性分析和迭代计算,获得最佳的穿流孔设计参数,从而避免流动滞死区和传热恶化的出现。
开孔范围和面积的确定:穿流孔的具体设计方法如下:首先通过对初始顺酐反应器结构的流场分析获得漩涡影响的区域,该区域作为需要布置穿流孔的区域,穿流孔的区域以多孔介质的方式加入到顺酐反应器的计算流体力学分析模型中。随后初步选定穿流孔的孔径和数量,通过下式计算穿流孔的多孔介质阻力系数。
其中:k1为多孔介质阻力系数,k0为局部阻力系数,a为沿程阻力系数,Re为雷诺数,l为折流板厚度,D为穿流孔直径,F0为总的孔内流通面积(即单个孔的面积乘以孔的数量),F1为折流板开孔区域的面积。其中:初步设定的穿流孔的数量和范围是根据以往工程经验,参考顺酐反应器折流板上的开孔直径和范围进行设置的。由于产能及反应器直径、冷却介质流量等参数的不同,以往的工程经验无法直接应用在新设计的反应器中;上述三个计算公式中: Re、l、D为已知或已设置的参数,k1、k0、a为需要计算求取的未知数,上述数据是“步骤S102设置管束区的阻力”的关键参数。根据流场、温度场模拟结果,进行判断,径向温差最大值是否不大于±2℃,轴向温差最大值是否不大于5℃,如果不满足上述要求,就需要根据计算结果进行评判,选择优化的方向,是增加或是减少孔径(孔径的范围为0~12mm);增加或减少开孔的范围(调整范围为0~1800mm);或两个同时调整。当模拟计算结果满足“径向温差最大值是否不大于±2℃,轴向温差最大值是否不大于5℃”时,所选取的孔径及开孔范围就是优化解。
S2.2调整折流板内径、折流板外径尺寸;
S2.3调整折流板间距;
S2.4调整中心管尺寸;
S2.2-S2.4 工作内容相同,统一阐述:
由于上述参数的调整均为体型参数,需先对上述参数(折流板内径、外径尺寸、折流板间距、中心管尺寸)进行其与温度场数值的敏感性分析,根据敏感系数及其计算过程,按范围和比例调整参数进行反复计算,选取最优解。综合上述影响因素的最优解所获得的管束区轴向和径向的温度分布,通过选取不同的横向和竖向界面,获得轴向温差的最大值和径向温差的最大值,判断是否满足径向温差最大值≤±2℃;轴向温差最大值≤5℃的设计要求。若选取值无法满足上述要求时,选取敏感系数高的影响因素进行调整,直至满足上述要求为止。
敏感性分析是识别特定参数对最终结果影响程度的有效辅助手段,可以提高模型计算的可靠性和准确性。这里指的特定参数就是指:折流板上穿流孔的开孔范围、折流板上穿流孔的面积、折流板内径、折流板外径、折流板间距和中心管尺寸;
以出口和入口的流量分配板开孔率(开孔数量)为例;
敏感性分析的方法有很多,本方法是以敏感系数作为评判敏感性的依据,敏感系数越高,则该参数与最终结果(温度场数值)的敏感性越高,对计算结果影响越大。
敏感系数的计算方法:
取倍数关系的5-7个“特定参数”,并根据这些参数计算得到对应的温度场数值。计算温度场数值/对应的“特定参数”,并绘制对应关系的折线图。
如果折线斜率的最大值与平均值的比值小于2,则取斜率的平均值作为敏感系数。
如果折线斜率的最大值与平均值的比值大于2,则在斜率最大值的数值范围内重复步骤(1)至(3),直至满足步骤(2)。
敏感系数的作用:
根据敏感系数计算过程,能够确定该“特定参数”在哪个范围和比例内进行调整最为合理。
敏感系数的计算结果,能够确定不同的“特定参数”哪个对温度场计算结果影响更大;
S3、设计换热介质的出口和入口的流量分配板开孔率(开孔数量);具体为:
通过计算流体力学的方法,确定外循环错流模型下不等面积的流量分配板尺寸的分布规律;在选取开孔直径为定值的前提下,获得不等面积流量分配板上的开孔率,对进入反应器的流量进行调节;最终得到开孔率S的规律为:
S=0.0024x2+0.0156x+0.3161;
其中,S为开孔率,x为流量分配板中心线和入口管中心线之间的夹角或流量分配板中心线和出口管中心线之间的夹角;单位为度。
按照上述规律开孔使得入口和出口处的流量分配均匀,进而消除由于环道处开孔位置不同而造成的熔盐流动的不均匀性。
S4、在熔盐泵正常流量的60%~100%范围内,更改熔盐流量,重复计算得到若干个温度场;
根据得到的各个温度场,选取温度最高点截面、径向温差最高点截面设置测温点;
S5、根据S2中敏感性分析方法,按照熔盐100%流量的工况设计顺酐反应器流道;具体包括:
首先将折流板上穿流孔的开孔范围、折流板上穿流孔的面积、折流板内径、折流板外径、折流板间距、中心管尺寸、出口和入口的流量分配板开孔率逐一作为敏感性因素进行敏感性分析,然后根据分析结果,优化得到每个敏感性因素的大小。
本发明通过实践数据修正模拟计算核心参数,提高模拟计算的精度和准确性;
通常的流体力学仿真模拟计算都是以大型商业软件自带的模型、边界条件为基础的,往往缺乏实践验证,用于大型反应器模拟计算时,容易“失真”,无法直接用于指导设计。本发明在计算的关键参数和节点上,通过采用实践数据“矫正”的方法,使模拟结果更加接近于实际工况,能够更好的指导工程设计。具体包括对①多孔介质阻力的修正②流场-温度场耦合的修正。
本发明利用计算流体动力学模拟的方式,优化了流量分配板设计计算方法。
传统的设计方法是根据动量距定理进行公式推导,得出流量分配板尺寸的理论推导数值。这种方法仅从入口管道几何尺寸、流量分配板所对应的角度等因素去考虑,通过控制流量分配板上开孔率(开孔面积)来使流体通过流量分配板之后的压降值相同的方法,达到使流体分配均匀的目的。不同直径、不同布管规律、折流板上穿流孔的分布范围等参数均会造成内部区域阻力的变化,进而影响上述计算方法的准确性。本设计方法通过模拟计算迭代之后,获得优化后流场的阻力等参数,并使用这一参数对上述公式进行计算修正,再对出入口的流量分配板进行开孔尺寸计算,相对于传统的计算方法,更加合理和精准。
本发明通过调整下列参数:①折流板上穿流孔的开孔范围和尺寸;②折流板内径、外径尺寸;③折流板间距;④中心管尺寸;⑤熔盐出、入口流量分配板尺寸;
对顺酐反应器内部流场、温度场进行优化设计;
本发明通过改变折流板开穿流孔区域局部多孔介质阻力、调整松弛因子的方式,实现降低原有超大规模网格模型计算量,保证计算收敛的目的。
本申请对顺酐反应器的计算模型进行后处理中,采用STAR-CCM+软件进行,获得顺酐反应器的流线图、各截面的速度分布、进出口流量分配板的流速分布等流场信息以及管束区的温度分布等信息。顺酐反应器的流线如图4所示,出口截面的速度云图如图5所示,这些流场和温度场计算结果对顺酐反应器的结构优化具有重要的参考价值。
以上针对顺酐反应器的流场-温度场耦合仿真计算可以获得可视化的流场和温度场分析结果,更重要的是获得折流板区和进出口流量分配板区域的精细化的流场信息以及管束区的温度场信息,为顺酐反应器的结构优化提供理论依据。
本申请根据顺酐反应器的压力、温度、流量等边界条件,使用通用的计算流体力学软件进行计算,并通过后处理获得反应器内部熔盐流动的流场和温度场,包括进出口流量分配板位置的流速分布、折流板区域流速分布、管束区温度场分布等信息,用于顺酐反应器的评估。
流场和温度场耦合计算:顺酐反应器管束区的化学反应放热通过在能量方程中添加热源项实现,该热源项的形式为单位体积单位时间的发热量,该热源随空间位置的变化以拟合函数的形式添加到能量方程中。
一种顺酐反应器流场设计系统,用于执行上述的设计方法,包括:
计算模块:计算顺酐反应器模型的流场和温度场;具体包括:
S101、根据顺酐反应器结构,建立计算流体力学模型;
S102、设置计算流体力学模型的参数:进口熔盐的质量流量和温度、管束区的阻力、松弛因子;具体为:
基于顺酐反应器的运行工况,设置进口熔盐的质量流量和温度;
根据顺酐反应器内部的管束排布方式,设置管束区的阻力;
S103、根据反应管热负荷试验数据,设置温度场计算的前置关联条件;根据试验得到的顺酐反应器管侧的化学反应特性,设置管束区的吸热功率、放热功率。
S104、对流场和温度场进行耦合分析计算:根据顺酐反应器入口的质量流量、熔盐的密度、运动粘度和管束区的流通面积、传热管的直径,计算管束区的雷诺数;根据雷诺数选择流动模型,若雷诺数小于104选择层流模型,若雷诺数不小于104选择湍流模型,通过模拟计算得到速度场分布;根据温度场计算的前置关联条件,得到温度场的分布;本实施例中,壁面采用无滑移边界,壁面函数采用two-layer all y+ treatment;动量方程添加重力源项;边界条件设定为质量流量入口、压力出口;将动量方程中速度参数的松弛因子调节至0.15;将动量方程中压力参数的松弛因子调节至0.03;将湍动能参数的松弛因子调节至0.15;将湍流耗散率参数的松弛因子调节至0.2;将能量参数的松弛因子调节至0.3;计算过程中监测质量、速度、湍动能和湍流耗散率的残差、进出口的压降和温差、计算域的最大速度和最大温度,当监测量稳定且波动幅度不大于1%时,认为计算收敛;
S105、基于边界条件完成流场和温度场的计算,并获得初始管束区轴向和径向的温度分布,通过选取不同的横向界面和竖向界面,获得轴向温差最大值和径向温差最大值,并判断径向温差最大值和轴向温差最大值是否满足设计要求,所述设计要求为:径向温差最大值是否不大于±2℃,轴向温差最大值是否不大于5℃;当初始设计结构不满足上述设计要求时,执行优化模块,当初始设计结构满足上述设计要求时,执行进出口设计模块;
优化模块:优化反应器内部流道;具体包括:
首先对流道参数进行敏感性分析,流道参数包括折流板上穿流孔的开孔范围、折流板上穿流孔尺寸、折流板内径、折流板外径、折流板间距和中心管尺寸;
然后对流道参数敏感性进行大小排序;
最后根据敏感性从大到小的排序,调整流道参数;调整完成后返回至S101;
进出口设计模块:设计换热介质的出口和入口的流量分配板;具体为:
通过计算流体力学的方法,确定外循环错流模型下不等面积的流量分配板尺寸的分布规律;在选取开孔直径为定值的前提下,获得不等面积流量分配板上的开孔率,对进入反应器的流量进行调节;最终得到开孔率S的规律为:
S=0.0024x2+0.0156x+0.3161;
其中,S为开孔率,x为流量分配板中心线和入口管中心线之间的夹角或流量分配板中心线和出口管中心线之间的夹角;
测温点选取模块:在熔盐泵正常流量的60%~100%范围内,更改熔盐流量,重复计算得到若干个温度场;
根据得到的各个温度场,选取温度最高点截面、径向温差最高点截面设置测温点;
流道设计模块:完成顺酐反应器流道设计;具体包括:
将折流板上穿流孔的开孔范围、折流板上穿流孔尺寸、折流板内径、折流板外径、折流板间距、中心管尺寸、出口和入口的流量分配板开孔率的数值作为顺酐反应器施工设计的数据。
一种信息数据处理终端,用于实现上述顺酐反应器的设计方法。
一种计算机可读存储介质,包括指令,当其在计算机上运行时,使得计算机执行上述顺酐反应器的设计方法。
在上述实施例中,可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用全部或部分地以计算机程序产品的形式实现,所述计算机程序产品包括一个或多个计算机指令。在计算机上加载或执行所述计算机程序指令时,全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中,或者从一个计算机可读存储介质向另一个计算机可读存储介质传输,例如,所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL)或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输)。所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质,(例如,软盘、硬盘、磁带)、光介质(例如,DVD)、或者半导体介质(例如固态硬盘Solid StateDisk(SSD))等。
以上所述仅是对本发明的较佳实施例而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施例所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。
Claims (9)
1.一种顺酐反应器的设计方法,其特征在于,包括:
S1、计算顺酐反应器模型的流场和温度场;具体包括:
S101、根据顺酐反应器结构,建立计算流体力学模型;
S102、设置计算流体力学模型的参数:进口熔盐的质量流量和温度、管束区的阻力、松弛因子;
S103、根据反应管热负荷试验数据,设置温度场计算的前置关联条件;
S104、对流场和温度场进行耦合分析计算:根据顺酐反应器入口的质量流量、熔盐的密度、运动粘度和管束区的流通面积、传热管的直径,计算管束区的雷诺数;根据雷诺数选择流动模型,若雷诺数小于104选择层流模型,若雷诺数不小于104选择湍流模型,通过模拟计算得到速度场分布;根据温度场计算的前置关联条件,得到温度场的分布;
S105、基于边界条件完成流场和温度场的计算,并获得初始管束区轴向和径向的温度分布,通过选取不同的横向界面和竖向界面,获得轴向温差最大值和径向温差最大值,并判断径向温差最大值和轴向温差最大值是否满足设计要求,所述设计要求为:径向温差最大值是否不大于±2℃,轴向温差最大值是否不大于5℃;当初始设计结构不满足上述设计要求时,执行S2,当初始设计结构满足上述设计要求时,执行S3;
S2、优化反应器内部流道;具体包括:
首先对流道参数进行敏感性分析,流道参数包括折流板上穿流孔的开孔范围、折流板上穿流孔尺寸、折流板内径、折流板外径、折流板间距和中心管尺寸;
然后对流道参数敏感性进行大小排序;
最后根据敏感性从大到小的排序,调整流道参数;调整完成后返回至S101;
S3、设计换热介质的出口和入口的流量分配板;具体为:
通过计算流体力学的方法,确定外循环错流模型下不等面积的流量分配板尺寸的分布规律;在选取开孔直径为定值的前提下,获得不等面积流量分配板上的开孔率,对进入反应器的流量进行调节;最终得到开孔率S的规律为:
S=0.0024x2+0.0156x+0.3161;
其中,S为开孔率,x为流量分配板中心线和入口管中心线之间的夹角或流量分配板中心线和出口管中心线之间的夹角;
S4、在熔盐泵正常流量的60%~100%范围内,更改熔盐流量,重复计算得到若干个温度场;
根据得到的各个温度场,选取温度最高点截面、径向温差最高点截面设置测温点;
S5、完成顺酐反应器流道设计;具体包括:
将折流板上穿流孔的开孔范围、折流板上穿流孔尺寸、折流板内径、折流板外径、折流板间距、中心管尺寸、出口和入口的流量分配板开孔率的数值作为顺酐反应器施工设计的数据。
2.根据权利要求1所述的顺酐反应器的设计方法,其特征在于,所述顺酐反应器包括圆柱形的空腔壳体,在空腔壳体的内腔设置有环形结构的折流板和中心管,在空腔壳体的侧壁开设有左右对称的入口管和出口管;S101具体为:
S1011、建立顺酐反应器的三维流场计算模型;
S1012、通过布尔运算,提取三维流场计算模型中的流体域模型;
S1013、将流体域模型沿X平面和Y平面进行切割,得到1/4的流体域模型作为计算流体力学分析模型;所述空腔壳体的中心轴位于X平面和Y平面内,且X平面和Y平面相互垂直。
3.根据权利要求1所述的顺酐反应器的设计方法,其特征在于,S102具体为:
基于顺酐反应器的运行工况,设置进口熔盐的质量流量和温度;
根据顺酐反应器内部的管束排布方式,设置管束区的阻力;
S103具体为:根据反应管热负荷试验数据,设置管束区的吸热功率、放热功率。
4.根据权利要求1所述的顺酐反应器的设计方法,其特征在于,S104中,壁面采用无滑移边界,壁面函数采用two-layer all y+ treatment;动量方程添加重力源项;边界条件设定为质量流量入口、压力出口;将动量方程中速度参数的松弛因子调节至0.15;将动量方程中压力参数的松弛因子调节至0.03;将湍动能参数的松弛因子调节至0.15;将湍流耗散率参数的松弛因子调节至0.2;将能量参数的松弛因子调节至0.3;计算过程中监测质量、速度、湍动能和湍流耗散率的残差、进出口的压降和温差、计算域的最大速度和最大温度,当监测量稳定且波动幅度不大于1%时,认为计算收敛。
5.一种顺酐反应器的设计系统,其特征在于,包括:
计算模块:计算顺酐反应器模型的流场和温度场;具体包括:
S101、根据顺酐反应器结构,建立计算流体力学模型;
S102、设置计算流体力学模型的参数:进口熔盐的质量流量和温度、管束区的阻力、松弛因子;
S103、根据反应管热负荷试验数据,设置温度场计算的前置关联条件;
S104、对流场和温度场进行耦合分析计算:根据顺酐反应器入口的质量流量、熔盐的密度、运动粘度和管束区的流通面积、传热管的直径,计算管束区的雷诺数;根据雷诺数选择流动模型,若雷诺数小于104选择层流模型,若雷诺数不小于104选择湍流模型,通过模拟计算得到速度场分布;根据温度场计算的前置关联条件,得到温度场的分布;
S105、基于边界条件完成流场和温度场的计算,并获得初始管束区轴向和径向的温度分布,通过选取不同的横向界面和竖向界面,获得轴向温差最大值和径向温差最大值,并判断径向温差最大值和轴向温差最大值是否满足设计要求,所述设计要求为:径向温差最大值是否不大于±2℃,轴向温差最大值是否不大于5℃;当初始设计结构不满足上述设计要求时,执行优化模块,当初始设计结构满足上述设计要求时,执行进出口设计模块;
优化模块:优化反应器内部流道;具体包括:
首先对流道参数进行敏感性分析,流道参数包括折流板上穿流孔的开孔范围、折流板上穿流孔尺寸、折流板内径、折流板外径、折流板间距和中心管尺寸;
然后对流道参数敏感性进行大小排序;
最后根据敏感性从大到小的排序,调整流道参数;调整完成后返回至S101;
进出口设计模块:设计换热介质的出口和入口的流量分配板;具体为:
通过计算流体力学的方法,确定外循环错流模型下不等面积的流量分配板尺寸的分布规律;在选取开孔直径为定值的前提下,获得不等面积流量分配板上的开孔率,对进入反应器的流量进行调节;最终得到开孔率S的规律为:
S=0.0024x2+0.0156x+0.3161;
其中,S为开孔率,x为流量分配板中心线和入口管中心线之间的夹角或流量分配板中心线和出口管中心线之间的夹角;
测温点选取模块:在熔盐泵正常流量的60%~100%范围内,更改熔盐流量,重复计算得到若干个温度场;
根据得到的各个温度场,选取温度最高点截面、径向温差最高点截面设置测温点;
流道设计模块:完成顺酐反应器流道设计;具体包括:
将折流板上穿流孔的开孔范围、折流板上穿流孔尺寸、折流板内径、折流板外径、折流板间距、中心管尺寸、出口和入口的流量分配板开孔率的数值作为顺酐反应器施工设计的数据。
6.根据权利要求5所述的顺酐反应器的设计系统,其特征在于,所述顺酐反应器包括圆柱形的空腔壳体,在空腔壳体的内腔设置有环形结构的折流板和中心管,在空腔壳体的侧壁开设有左右对称的入口管和出口管;S101具体为:
S1011、建立顺酐反应器的三维流场计算模型;
S1012、通过布尔运算,提取三维流场计算模型中的流体域模型;
S1013、将流体域模型沿X平面和Y平面进行切割,得到1/4的流体域模型作为计算流体力学分析模型;所述空腔壳体的中心轴位于X平面和Y平面内,且X平面和Y平面相互垂直。
7.根据权利要求5所述的顺酐反应器的设计系统,其特征在于,S102具体为:
基于顺酐反应器的运行工况,设置进口熔盐的质量流量和温度;
根据顺酐反应器内部的管束排布方式,设置管束区的阻力;
S103具体为:根据反应管热负荷试验数据,设置管束区的吸热功率、放热功率。
8.根据权利要求5所述的顺酐反应器的设计系统,其特征在于,S104中,壁面采用无滑移边界,壁面函数采用two-layer all y+ treatment;动量方程添加重力源项;边界条件设定为质量流量入口、压力出口;将动量方程中速度参数的松弛因子调节至0.15;将动量方程中压力参数的松弛因子调节至0.03;将湍动能参数的松弛因子调节至0.15;将湍流耗散率参数的松弛因子调节至0.2;将能量参数的松弛因子调节至0.3;计算过程中监测质量、速度、湍动能和湍流耗散率的残差、进出口的压降和温差、计算域的最大速度和最大温度,当监测量稳定且波动幅度不大于1%时,认为计算收敛。
9.一种信息数据处理终端,其特征在于,用于实现权利要求1-4任一项所述的顺酐反应器的设计方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310686938.2A CN116415449B (zh) | 2023-06-12 | 2023-06-12 | 顺酐反应器的设计方法、设计系统及信息数据处理终端 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310686938.2A CN116415449B (zh) | 2023-06-12 | 2023-06-12 | 顺酐反应器的设计方法、设计系统及信息数据处理终端 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116415449A true CN116415449A (zh) | 2023-07-11 |
CN116415449B CN116415449B (zh) | 2023-08-22 |
Family
ID=87049606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310686938.2A Active CN116415449B (zh) | 2023-06-12 | 2023-06-12 | 顺酐反应器的设计方法、设计系统及信息数据处理终端 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116415449B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101898103A (zh) * | 2010-02-24 | 2010-12-01 | 南京钟腾化工有限公司 | 年产2万吨顺酐的熔盐内循环式反应装置和反应工艺 |
CN106279078A (zh) * | 2015-06-29 | 2017-01-04 | 常州瑞华化工工程技术有限公司 | 制备顺酐的方法和低压降径向等温反应器 |
CN108280300A (zh) * | 2018-01-24 | 2018-07-13 | 南京罕华流体技术有限公司 | 基于计算流体力学的大数据流量计开发方法 |
CN108345714A (zh) * | 2018-01-11 | 2018-07-31 | 武汉科技大学 | 一种内环向射流稳压腔参数设计的数值模拟方法 |
CN109433116A (zh) * | 2018-12-29 | 2019-03-08 | 常州瑞凯化工装备有限公司 | 用于强放热化学反应过程的管壳式轴向列管反应器 |
CN113280669A (zh) * | 2021-06-08 | 2021-08-20 | 中国科学院理化技术研究所 | 折流板的设计方法及内置折流板的蓄冷/热器 |
CN114444413A (zh) * | 2022-01-21 | 2022-05-06 | 西安交通大学 | 一种板状燃料堆芯亚通道级三维热工水力分析方法 |
CN115421460A (zh) * | 2022-09-15 | 2022-12-02 | 重庆大学 | 基于计算机数值模拟的铸造残余应力控制优化方法及应用 |
CN115659760A (zh) * | 2022-11-16 | 2023-01-31 | 河北省机电一体化中试基地有限公司 | 一种基于Fluent的热熔胶胶枪分析和结构优化方法 |
CN115738921A (zh) * | 2022-11-29 | 2023-03-07 | 东方电气集团东方锅炉股份有限公司 | 一种反应器进出口多腔室均匀冷却的列管式顺酐反应器系统 |
-
2023
- 2023-06-12 CN CN202310686938.2A patent/CN116415449B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101898103A (zh) * | 2010-02-24 | 2010-12-01 | 南京钟腾化工有限公司 | 年产2万吨顺酐的熔盐内循环式反应装置和反应工艺 |
CN106279078A (zh) * | 2015-06-29 | 2017-01-04 | 常州瑞华化工工程技术有限公司 | 制备顺酐的方法和低压降径向等温反应器 |
CN108345714A (zh) * | 2018-01-11 | 2018-07-31 | 武汉科技大学 | 一种内环向射流稳压腔参数设计的数值模拟方法 |
CN108280300A (zh) * | 2018-01-24 | 2018-07-13 | 南京罕华流体技术有限公司 | 基于计算流体力学的大数据流量计开发方法 |
CN109433116A (zh) * | 2018-12-29 | 2019-03-08 | 常州瑞凯化工装备有限公司 | 用于强放热化学反应过程的管壳式轴向列管反应器 |
CN113280669A (zh) * | 2021-06-08 | 2021-08-20 | 中国科学院理化技术研究所 | 折流板的设计方法及内置折流板的蓄冷/热器 |
CN114444413A (zh) * | 2022-01-21 | 2022-05-06 | 西安交通大学 | 一种板状燃料堆芯亚通道级三维热工水力分析方法 |
CN115421460A (zh) * | 2022-09-15 | 2022-12-02 | 重庆大学 | 基于计算机数值模拟的铸造残余应力控制优化方法及应用 |
CN115659760A (zh) * | 2022-11-16 | 2023-01-31 | 河北省机电一体化中试基地有限公司 | 一种基于Fluent的热熔胶胶枪分析和结构优化方法 |
CN115738921A (zh) * | 2022-11-29 | 2023-03-07 | 东方电气集团东方锅炉股份有限公司 | 一种反应器进出口多腔室均匀冷却的列管式顺酐反应器系统 |
Non-Patent Citations (1)
Title |
---|
冯宏 等: ""正丁烷法制顺酐氧化反应器温度调节控制装置探究"", 《天津化工》, vol. 33, no. 3, pages 48 - 49 * |
Also Published As
Publication number | Publication date |
---|---|
CN116415449B (zh) | 2023-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cheng et al. | Surrogate based multi-objective design optimization of lithium-ion battery air-cooled system in electric vehicles | |
CN113255229B (zh) | 一种基于联合仿真的燃料组件多学科结构设计优化方法 | |
WO2022099713A1 (zh) | 一种基于 overset 模型的低温碳化炉丝束加热过程三维模拟方法 | |
Balaton et al. | Operator training simulator process model implementation of a batch processing unit in a packaged simulation software | |
CN112699620A (zh) | 基于计算流体力学的反应堆堆芯热工水力特性分析方法 | |
Qu et al. | A combined scheme of generalized finite difference method and Krylov deferred correction technique for highly accurate solution of transient heat conduction problems | |
CN111027112B (zh) | 一种针对快堆棒束组件耦合传热模型的多孔介质模拟方法 | |
CN106897520B (zh) | 一种含有模糊参数的传热系统可靠性分析方法 | |
CN115659908B (zh) | 一种印刷电路板换热器的三单元非平衡多孔介质方法 | |
CN106844853B (zh) | 结合阻力和能量分布包含格架搅混效应的子通道分析方法 | |
CN114266171B (zh) | 一种u型管蒸汽发生器全耦合共轭传热计算方法 | |
CN112528572B (zh) | 一种基于overset模型的低温碳化炉丝束加热过程三维模拟方法 | |
CN104657589B (zh) | 一种管壳式换热器多孔介质系数计算方法 | |
WO2024016621A1 (zh) | 反应堆试验模型的规模确定方法、装置和计算机设备 | |
CN111611698B (zh) | 一种超薄体传热仿真方法 | |
CN114444413A (zh) | 一种板状燃料堆芯亚通道级三维热工水力分析方法 | |
CN116415449B (zh) | 顺酐反应器的设计方法、设计系统及信息数据处理终端 | |
Shahri et al. | Multi-objective optimization to minimize pumping power and flow non-uniformity at the outlets of a distributor manifold using CFD simulations and ANN rapid predictions | |
Lin et al. | CFD investigation for a 7-pin wrapped-wire fuel assembly with different wires | |
CN115017674A (zh) | 一种基于Modelica的管壳式换热器建模仿真方法 | |
CN110728072B (zh) | 一种确定数字反应堆计算流体力学分析网格尺寸的方法 | |
CN110728030B (zh) | 基于表面式凝汽器数值模拟汽水两侧耦合的方法 | |
CN116757108B (zh) | 换热器的模拟方法、装置及设备 | |
CN114462336B (zh) | 一种核反应堆主管道冷却剂平均温度计算方法 | |
CN116227060A (zh) | 一种基于正交试验的管壳式换热器结构优化分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |