基于SIMULINK平台的吸附式制冷机模块化建模与仿真.pdf

发布时间：2022-05-30 发布人：admin 分类：说明书资料大小：0.73M 资料格式：pdf 举报版权申诉

weixin_38551938-12114479-16359647438751677562.pdf-第1页.png

第1页 / 共8页

weixin_38551938-12114479-16359647438751677562.pdf-第2页.png

第2页 / 共8页

weixin_38551938-12114479-16359647438751677562.pdf-第3页.png

第3页 / 共8页

weixin_38551938-12114479-16359647438751677562.pdf-第4页.png

第4页 / 共8页

weixin_38551938-12114479-16359647438751677562.pdf-第5页.png

第5页 / 共8页

weixin_38551938-12114479-16359647438751677562.pdf-第6页.png

第6页 / 共8页

weixin_38551938-12114479-16359647438751677562.pdf-第7页.png

第7页 / 共8页

weixin_38551938-12114479-16359647438751677562.pdf-第8页.png

第8页 / 共8页

文本预览

http://www.paper.edu.cn 基于 SIMULINK 平台的吸附式制冷机模块化建模与仿真* 林芃，王如竹，吴静怡上海交通大学制冷与低温工程研究所，上海（200240） E-mail：linpeng0478@sjtu.edu.cn 摘要：文章运用集总参数法建立了吸附式制冷机动态数学模型，以 Matlab-Simulink 为仿真平台实现了 10kW 级硅胶-水回热回质型吸附式制冷机的模块化仿真，考察了制冷机运行过程中吸附床温、冷却水、冷冻水出口水温的动态变化以及回热回质对于系统运行的作用；研究了不同运行工况对制冷量以及 COP 的影响。为冷热电联产系统的模块化全过程仿真奠定了良好的基础。关键词：吸附式制冷；模块化；仿真；动态性能；冷热电联产系统中图分类号：TB651 1. 引言冷热电联供（CCHP）系统是将发电、制冷、供暖设备进行集成与组合的先进用能方式，能够实现能源的梯级利用。上海交通大学采用 16kW 微型燃气内燃发电机和新颖的热驱动吸附式制冷机，构成微型冷热电联产系统。图 1 采用燃气内燃机与吸附式制冷机的微型冷热电联产系统由于 CCHP 系统设备类型多、过程耦合复杂，建立各子系统动态模型，实现系统全过程仿真非常必要，利用 Matlab-Simulink 平台建立热驱动制冷机动态模型，可方便的与发电、供热设备进行连接与扩展，研究 CCHP 系统动态性能。 Simulink 是用来对动态系统进行建模、仿真和分析的软件包，它提供了一种图形化的交互环境，非常直观，容易掌握，可利用 Matlab 中的丰富资源，建立分层的多级仿真模型[1]。若具备了 CCHP 系统的各种部件模型，只需简单连接就可实现系统的动态仿真。硅胶-水吸附式制冷机能够有效利用低品位热能，可以被 65 C 的热源驱动，是一种环境友好的制冷方式。上海交通大学制冷与低温工程研究所已研制出 10kW 级制冷量的余热型热水驱动吸附式制冷机 [2]，可将其作为冷热电联产系统的制冷部件，建立微型冷热电联供系统。本文基于 Simulink 平台对硅胶-水吸附式制冷机进行模块化建模与仿真，研究制冷机的动态变化规律以及变工况下的运行特性。 2. 吸附式制冷机系统描述图 2 为上海交通大学发明的硅胶-水吸附式冷水机组示意图。它由两个单床吸附式制冷系统复合而成的双床连续制冷系统。整个系统由三个真空腔组成，左右为两个由吸附床、冷凝器和蒸发器组成的吸附/解吸工作腔，底部为热管隔离蒸发器的工作腔。吸附床为管翅式紧凑换热器，冷 *本课题得到教育部博士点基金（项目编号：20040248055）的资助。 -1-

凝器采用管壳式换热器，蒸发器采用了热管隔离技术使发生器解吸侧与冷媒水传热管之间实现热隔离。了完整的吸附式制冷循环，制冷机可以连续输出冷冻水。 http://www.paper.edu.cn 该机组可自动以回热回质吸附式制冷 hot water outlet 循环方式运行，两个吸附床交替进行吸附制冷，主要包括以下过程：（1）左床解吸、右床吸附过程。驱动热源热水通入左侧吸附床中，使得左侧吸附床内吸附剂升温，左真空腔内蒸汽压力升高，当压力超过左冷凝器温度对应的饱和蒸汽压时，左冷凝器开始冷凝制冷剂；此时，左隔离器蒸发面不工作，其温度升高至冷凝温度，该温度高于热管工作腔内的蒸发器温度，从而实现了左隔离器与蒸发器之间的热隔离。与此同时，制冷剂进入右吸附床，右吸附床被冷却水冷却，开始降温吸附，右真空腔内的制冷剂蒸汽压力随之下降，当压力低于右隔离器温度对应的饱和蒸汽压时，右隔离器的蒸发面开始（1）蒸发制冷，其温度迅速降低，热管工作腔底部的蒸发器蒸发出来的热管工质蒸汽在（2）右隔离器传热表面上凝结，从而输出制冷量。（2）从左到右的回质过程。当左床解吸/ 右床吸附过程临近结束时，回质真空阀打（3）（4）开，左腔内的制冷剂蒸汽就会在较大的压差作用下迅速流到右腔，导致左隔离器温度降低而右隔离器温度升高，左右两腔体（5）内的压力迅速趋于平衡。同时，左吸附床解吸出的制冷剂蒸汽通过回质阀流入右吸附床内被吸附，实现二次解吸和吸附过程。（3）从左到右的回热过程。当左右腔体内的压力接近平衡时，关闭回质真空阀，打开相应阀门进行两个吸附床之间的回热。左吸附床内的热水进入右吸附床中，将其中的冷水排出，同时右吸附床内的冷水进入左吸附床中，将其中的热水排除，完成回热过程[3]。相应的，还应包括右床解吸、左床吸附过程；从右到左的回质过程；从右到左的回热过程，在此不再赘述。这样才组成 -2- hot water inlet valve components cooling water inlet cooling water inlet A remained water methanol chilled water inlet adsorber adsorber condenser condenser isolator isolator vacuum valve for mass recovery B remained water water evaporator C chilled water outlet remained water 图 2 硅胶-水吸附制冷机系统结构示意图 3. 吸附式制冷系统数学模型 3.1 系统基本假设系统仿真数学模型采用集总参数法。为了简化计算，本文做出了如下假设：整个吸附床内部的温度和蒸汽压力都是均匀的；制冷剂被吸附剂均匀地吸附，并且在吸附剂内部凝聚为液体；忽略吸附床与冷凝器或蒸发器之间的压差；除了热水、冷却水和冷冻水与外界有换热外，忽略系统散失到环境中的热/ 冷量。冷凝器中制冷剂蒸汽冷凝为饱和水；蒸发器中的制冷剂蒸发成饱和制冷剂蒸汽。 3.2 制冷机动态数学模型 3.2.1 吸附方程采用非平衡吸附率方程[4]： dx )] d τ （1） E exp[ − R 2 p D s /( RT b a 15 = 0 ⋅ * ( x − x ) 式中表面扩散系数 Ds0 ＝ 2.54×10-4m2/s，硅胶表明的活化能 Ea＝ 4.2×104J/mol，硅胶颗粒的平均半径 Rp＝ 7.1×10-4m，R 为水蒸气的气体常数。

*x 由经典平衡吸附方程确定[5]： * x = TA ( b [) ⋅ TP ( w s TP ( s b ) ) （2） TB ( b ) ] 式中 Ps( wT )和 Ps( bT )分别为水蒸气温 0 = = + + + + + + A 0 B （4）（3） TA b1 TB b1 TA 3 b3 TB 3 b3 TA 2 b2 TB 2 b2 度和硅胶温度对应下水的饱和压力；其中 TA ( b ) TB b ) ( 上式中的 A0…B3 为常数，取值如下： A0 ＝ -6.5314 ， A1 ＝ 0.72452×10-1 ， A2 ＝ -0.23951×10-3，A3＝0.25493×10-6； B0＝-15.587，B1＝0.15915，B2＝-50612×10-3， B3＝0.53290×10-6。水的饱和压力(kPa)和饱和温度(K)之间的关系： TPs )( )15. .0 273 0013802 − − T ( .0 0857427 )15. + − 273 3 T ( .0 0000888 = T ( )15.273 2 − .0 4709375 + （5） http://www.paper.edu.cn , out = T − a T − a T ad T ad KAad 为吸附床的换热性能系数，kW/K。 KA − ad cm （7） exp( wpw in . ) , 3.2.3 冷凝器能量平衡方程冷凝器中没有流体滞留，输入热量为制冷剂带来的显热和潜热，输出热量为被冷却水带走的热量。 dTMc [ c δ ⋅ ⋅−= cu 1 d τ cmTT ) ( ⋅ +− cool c dxMc dxML des des ⋅+ wv d d τ τ T T ( )] ⋅ − cool in cool out , , wp . ⋅ ⋅ ⋅ ⋅ a a a c （8） ) （9） out , T cool T cool T − c T − c = exp( − KA c c ⋅ m wp , cool in . Mc 冷凝器金属的质量，15.2kg；L 为水的汽化潜热，2500kJ/kg；Tc 冷凝器温度，K； Tcool,in 和 Tcool,out 分别为冷凝器进出口水温， K； coolm 为加热/冷却流体的质量流量， 1.11kg/s；KAc 为冷凝器的换热性能系数， kW/K。 a ⋅ + c wp , 3.2.2 吸附床能量平衡方程[2] d {[ d τ qM ⋅ st McMc x ) ⋅+⋅ cu dxMc ⋅ d τ T ) ad out , cM ( a dx d τ T ( ⋅ inad cm ⋅ w ⋅+ al 1( −+ ) δ 1 adt , − + = wp . wv ⋅ ⋅ ⋅ ⋅ a a , ] ⋅ T } b fin , ad TT ( − e a ) （6） Ma 为单个吸附床内吸附剂的质量， 50kg；ca 为吸附剂的比热，kJ/(kg·K)；cp,w 和 cwv 分别为水和水蒸汽的比热，kJ/(kg·K)； ccu 和 cal 分别为传热管材料铜和翅片材料铝的比热，kJ/(kg·K)；Mt,ad 为传热管质量， 21.8kg；Mfin,ad 为翅片质量，10.9kg；qst 为等压吸附/解吸热，2800kJ/kg；Te 和 Ta 分别为蒸发器和吸附床温度，K；Tad,in 和 Tad,out 分别为吸附床进出口水温，K； wm 为加热/冷却流体的质量流量，1.11kg/s。 1δ = 1, ⎧ ⎨ 0, ⎩ process desorption adsorption process 3.2.4 蒸发器能量平衡方程蒸发器中有液体制冷剂滞留，输入蒸发器的热量包括冷媒水带入的热量和冷凝器输送制冷剂所带入的热量；输出蒸发器的热量包括制冷剂蒸发吸热所带走的热量。在系统仿真过程中，不考虑蒸发器内 ⋅ ⋅+ cu 部复杂的传热过程，把热管蒸发器作为一个整体考虑，能量平衡方程为： d TMcMc [( wp e , d τ T cm [ )] ⋅+ − + ⋅ δδ out chill chill 1 2 , dxM dxML ) des des ⋅ δ 2 d d τ τ dxML )[ ads ⋅− δ a 1 d τ TT ) ( ⋅−⋅ e c we , T ( in chill , 1(] −=⋅ 1( −− c wp , wp . ) ⋅ ⋅ ⋅ ⋅ ⋅ a a e （10） ≤ > T e T e （11） exp( − KA e c ⋅ m chill wp , ) （12） = T 1, ⎧ c ⎨ T 0, ⎩ c T − e T − e = 其中 δ 2 , T chill T chill out in . -3-

Me 为蒸发器金属的质量，65kg；Me,w 为蒸发器内液体制冷剂(水)的质量，kg； Tchill,in 和 Tchill,out 分别为蒸发器冷冻水进出口水温，K； chillm 为冷冻水质量流量，0.5kg/s， KAe 为蒸发器的换热性能系数，kW/K。 3.2.5 回质过程 http://www.paper.edu.cn a wv 2 )( v wv − b ) = RT wv + ( P wv v （18） 3.2.6 回热过程回热过程中，阀门切换，两床传热介质连通，通过流动介质传热，吸附床驻留的冷在回质过程中，冷凝器处于空闲状态，水进入解吸床进行冷却，而解吸中驻留的热解吸腔内的蒸发器将开始蒸发，吸附腔内的水进入吸附床进行加热。吸附床、冷凝器和蒸发器开始出现冷凝现象。吸附床的吸附将蒸发器各换热器的能量方程在形式上与正由吸附床与解吸腔之间的压差驱动，解吸床常情况下形式一致，此时，吸附床的进口水的解吸由解吸床与吸附腔之间的压差驱动，对应吸附方程中的 Ps 分别取解吸腔和吸附腔的压力，除此之外，描述吸附/解吸床的传热模型方程与其他过程中的形式完全一致。最大的不同是蒸发器的能量平衡方程： d Mc wp , d τ δ m + 4 其中： m ⎧ ⎨ m ⎩ chamber desorbing adsorbing chamber TMc ⋅ cu e T ( chill for , , for （13） L ⋅−= T chill δ 3 cond e, evap e, δ 3 )] [( = chill wp . we , − + c , out , in ] ) e ⋅ ⋅ ⋅ ⋅ ⋅ （14）和 e,condm 分别为回质过程中解吸腔 e,evapm 内水蒸发器的蒸发流量和吸附腔内水蒸发器的冷凝流量，kg/s。温是解吸床的出口水温，解吸床的进口水温则是吸附床的出口水温，因此，在回热过程中有以下特征关系式： T ads , = in T des , out ； Tdes = in, T ads , out 3.2.7 系统性能指标的求取制冷量： τ∫ = 0 Q ref cycle c wp , ⋅ m chill − T chill out , ) d τ in , T ( ⋅ chill τ cycle （19） cycleτ 为循环周期。加热量： τ∫ = 0 wp , c cycle ⋅ Q h inh , Tm ( ⋅ h τ cycle − T h , out ) d τ （20） δ 4 = T 1, ⎧ e ⎨ T 0, ⎩ e ≤ > T chill T chill , in , in （15） COP = Q ref Q h （21）制冷剂质量平衡方程： dxM des d τ m e m e evap − = + a , , cond + dxM ads d τ ⋅ a = m mr （16） mrm 为回质流量，kg/s。解吸腔与吸附腔内压力差： P wv , des − P wv , abs = mv 2 wv mr 2A 2 （17） vwv 为水蒸气的比容，m3/kg；A 为回质通道的截面积，m2。添加范德瓦尔斯状态方程使方程组封闭： -4- 4. 基于 Simulink 平台系统建模 4.1 利用 Simulink 建立吸附式制冷机动态模型 Simulink 是一个交互式动态系统建模、仿真和分析工具，用户可采用交互或批处理的方式，方便得更换参数来进行“what-if”式的分析仿真，在控制系统、DSP 系统、通信系统中已经有了很广泛的应用[1]。在热力系统中有非常多的集中参数模型，其数学形式表现为耦合的微分方程组，通常我们将微分

http://www.paper.edu.cn 方程进行离散，利用 C、Fortran 或 Matlab 进行编程，其实利用 Simulink 来求解这些方程组是非常直观和方便的。另外，Simulink 模型的模块化使它很容易进行扩展，若已具图 3 所示。将部件模型进行子系统封装，即可完成吸附式制冷机的模块化仿真，将其热水入口与烟气换热器模型相连。烟气换热器的烟气有燃气发电机，蓄热水箱，吸附式制冷机以入口可与今后将要完成的燃气发电机模型及换热器模块，只需用鼠标进行简单连接就可以实现 CCHP 系统的模块化仿真。利用 Simulink 建立的吸附式制冷机部件模型如进行连接，这样就完成了冷热电联产系统的模块化仿真，如图 4。图 3 吸附式制冷机部件模型图 4 CCHP 系统模块化仿真 5. 系统仿真结果与分析 5.1 吸附式制冷机全过程动态仿真系统采用变步长 ode45 求解器（龙格库塔 4、5 阶算法）求解耦合微分方程组，仿真额定工况：热源温度 85℃，冷却水入 -5- 口温度 30℃，冷冻水入口温度 20℃，循环周期 2400s，其中回质时间 180s，回热时间 20s。仿真结果：Qref=9.116kW，COP=0.412 实验结果：Qref=9.38kW，COP=0.43[6] 吸附床温度动态变化如图 5 所示：

http://www.paper.edu.cn 图 5 吸附床温度动态变化过程 0-1200s 为制冷机启动时运行状态，由于是初启动，无回热回质过程；1200s-1380s 为回质过程，1380s-1400s 为回热过程。1400s 后为将外界热源冷源重新接入后的吸附床的加热、冷却过程。1200s-2400s 为正常运行情况下的半个循环周期（区别于初始开机过程）。由于传热传质效率所限，吸附床吸附/解析过程时间较长，每个吸附床需要 40 分钟完成加热解析/冷却吸附的制冷循环。用于加热/冷却吸附床的热水/冷水出口温度动态变化，如图 6：图 6 加热/冷却吸附床的热水/冷水出口温度如图 6 所示，吸附床热/冷水出口温度变化与吸附床温度变化趋势大致相同。冷凝器冷却水出水以及蒸发器冷冻水出水温度动态变化如图 7 所示：图 7 冷却水出水以及冷冻水出水温度从图 7 中 1200s-1300s 可清楚地看到回质回热过程对于冷却水出水以及冷冻水出水温度产生的波动。图 8 制冷机运行过程中吸附率的变化图 8 表明了制冷机工作过程中吸附床吸附率的变化，为了提高单位循环制冷量以及 COP，采用了回质提高吸附床的吸附量；采用回热节约了外界的热量输入。这两种方式显著提高了系统的制冷性能。从图中可看到吸附床吸附量在 0.23-0.07kg/kg 之间循环，回质过程实现了二次解析和吸附。 5.2 变工况对系统制冷量及 COP 的影响冷热电联产系统是一个多种设备协调运行的复杂系统，运行过程中工况经常会发生变化，因此，研究变工况对于部件运行性能的影响是非常有必要的。改变系统工况时，变化热源、冷却水或冷冻水其中一个量，其余量保持额定工况不变。 5.2.1 改变热源温度 -6-

http://www.paper.edu.cn 表 3 冷冻水进口水温对制冷机性能的影响 COP 冷冻水进口 refQ kW 水温℃ 10 12 14 16 18 20 22 5.117 5.953 6.761 7.557 8.351 9.116 9.876 0.345 0.362 0.376 0.388 0.4 0.412 0.424 冷冻水入口水温对 COP 的影响也比较大，随着冷媒水入口水温的提高 refQ 和 COP 都显著增大。 6. 结论本文通过吸附式制冷机动态数学模型，基于 Simulink 平台实现了制冷机模块化仿真，模拟了吸附式制冷机动态运行过程，以及变工况对于系统性能的影响，为冷热电联产系统模块化全过程仿真奠定了基础。具体的结论如下：（1）通过 Simulink 平台可以实现吸附式制冷的动态仿真，并且可与原动机连接组合成冷热电联产系统进行仿真研究。（2）分析仿真计算结果，本文研究的吸附式制冷机标准工况下制冷量为 9.116kW， COP 为 0.412，循环周期 40 分钟，如何提高系统的传热传质效率是吸附式制冷机取得突破的关键。（3）相比加热热源温度，冷却水与冷冻水入口水温对系统 COP 影响更大，当系统运行在变工矿状态下，应注意对冷却水与冷冻水入口水温的控制。表 1 热源温度对制冷机性能影响热水温度℃ 65 70 75 80 85 refQ kW 5.886 6.847 7.61 8.398 9.116 COP 0.376 0.398 0.404 0.409 0.412 随着热源温度的升高，制冷量迅速增大，COP 也呈增大趋势，但增大并不明显，表明硅胶-水吸附式制冷机具有能够较好利用低品位热源的特性。 refQ 和 COP 的增加幅度随着热源温度的提高而逐渐减小。在冷热电联产系统中，若发电机烟气排气温度降低，根据仿真结果表明吸附机的 COP 下降不大，变工况时比较稳定，因此将吸附机作为冷热电联产的制冷部件是合适的。 5.2.2 改变冷却水进口水温表 2 冷却水进口水温对制冷机性能的影响 COP 冷却水进口 refQ kW 水温℃ 26 28 30 32 34 11.675 10.345 9.116 7.849 6.61 0.439 0.426 0.412 0.398 0.385 随着冷却水进口水温提高， refQ 下降， COP 也迅速下降，相比热源，冷却水水温的波动对 COP 的影响更大，因此夏季工况时冷却水温的控制尤为重要。 5.2.3 改变冷冻水进口水温 -7-

http://www.paper.edu.cn 参考文献 [1] 姚俊,马松辉.Simulink 建模与仿真[M].西安：西安电子科技大学出版社，2002：9-12. [2] D.C. Wang, Z.Z. Xia, J.Y. Wu, R.Z. Wang, H. Zhai, W.D. Dou. Study of a novel silica-gel water adsorption chiller. Part I. Design and performance prediction[J]. International Journal of Refrigeration,2005,28(7)：1073-1083. [3] 王如竹,代彦军.太阳能制冷[M].北京：化学工业出版社,2007：220-222. [4] A. Sokoda, M. Suzuki. Fundamental study on solar powered adsorption cooling system[J]. Journal of Chemical Engineering Japan,1984,17(1)：52-57. [5] B.B. Saha, S. Koyama, J.B. Lee, K. Kuwahara, K.C.A. Alam, Y. Hamamoto, A. Akisawa, T. Kashiwagi. Performance evaluation of a low-temperature waste heat driven multi-bed adsorption chiller[J]. International Journal of Multiphase Flow 2003,29：1249-1263. [6] J. Di, J.Y. Wu, Z.Z. Xia, R.Z. Wang. Theoretical and experimental study on characteristics of a novel silica gel-water chiller under the conditions of variable heat source temperature[J]. International Journal of Refrigeration, 2007, 30：515-526. [7] R.Z. Wang, J.Y. Wu, Y.J. Dai, W. Wang, et al. Adsorption Refrigeration[J], China Machine Press, 2002. [8] 王如竹,吴静怡,代彦军,王文,姜周曙.吸附式制冷[M].北京：机械工业出版社,2001 [9] 吴静怡,王如竹,许煜雄.连续回热型吸附式空调/热泵循环特性与动态性能(I)模型与仿真[J].化工学报,2002,53(2)：144-149. Modularized Modeling And Simulation Of An Adsorption Chiller Based On Matlab-Simulink Institute of Refrigeration and Cryogenics, Shanghai Jiao Tong University, Shanghai, Lin Peng, Wang Ruzhu, Wu Jingyi China (200240) Abstract In this article, a dynamic model for the 10kW silica gel-water adsorption chiller is developed using lumped parameters method. Matlab-Simulink is used as the tool to realize the modularized simulation of this chiller based on its transient model. From the simulation results, the temperature of the adsorption bed, the outlet temperature of cooling and cold water and the impact of the heat and mass recovery are investigated. Also the influence on refrigeration capacity and COP under different working conditions is discussed in the article. The work has laid a foundation for further analysis and simulation of the CCHP (combined cooling, heating and powering) system. Keywords: adsorption chiller; modularization; simulation; dynamic performances; CCHP system -8-

分享到：

赞收藏

资料库

基于SIMULINK平台的吸附式制冷机模块化建模与仿真.pdf

相关推荐

开发技术

热门标签

最新资料