CN114154262B - 一种交叉流印刷电路板式换热器正向设计方法 - Google Patents

一种交叉流印刷电路板式换热器正向设计方法 Download PDF

Info

Publication number
CN114154262B
CN114154262B CN202111454588.4A CN202111454588A CN114154262B CN 114154262 B CN114154262 B CN 114154262B CN 202111454588 A CN202111454588 A CN 202111454588A CN 114154262 B CN114154262 B CN 114154262B
Authority
CN
China
Prior art keywords
primary side
secondary side
flow
heat exchange
inlet
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.)
Active
Application number
CN202111454588.4A
Other languages
English (en)
Other versions
CN114154262A (zh
Inventor
李昕泽
张大林
李新宇
田文喜
秋穗正
苏光辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN202111454588.4A priority Critical patent/CN114154262B/zh
Publication of CN114154262A publication Critical patent/CN114154262A/zh
Application granted granted Critical
Publication of CN114154262B publication Critical patent/CN114154262B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F28—HEAT EXCHANGE IN GENERAL
    • F28D—HEAT-EXCHANGE APPARATUS, NOT PROVIDED FOR IN ANOTHER SUBCLASS, IN WHICH THE HEAT-EXCHANGE MEDIA DO NOT COME INTO DIRECT CONTACT
    • F28D9/00—Heat-exchange apparatus having stationary plate-like or laminated conduit assemblies for both heat-exchange media, the media being in contact with different sides of a conduit wall
    • F28D9/0031—Heat-exchange apparatus having stationary plate-like or laminated conduit assemblies for both heat-exchange media, the media being in contact with different sides of a conduit wall the conduits for one heat-exchange medium being formed by paired plates touching each other
    • F28D9/0037—Heat-exchange apparatus having stationary plate-like or laminated conduit assemblies for both heat-exchange media, the media being in contact with different sides of a conduit wall the conduits for one heat-exchange medium being formed by paired plates touching each other the conduits for the other heat-exchange medium also being formed by paired plates touching each other
    • 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
    • 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)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Thermal Sciences (AREA)
  • Mechanical Engineering (AREA)
  • Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)

Abstract

本发明公开了一种交叉流印刷电路板式换热器正向设计方法，该方法设计对象结构包括换热器芯体、冷热侧接口腔体，其中芯体采用一次侧、二次侧板交替布置；该方法基于规则结构采用标准单元法进行一、三维耦合数值计算，形成一次侧与二次侧任意单元流体温度的数值解，得到换热板组温度分布，并参照温度分布结果对换热板进行参数设计；采用有限元对换热板组最大热应力进行验证，并结合多孔介质模型确定一次侧与二次侧工质流动分布及各流道均匀性；由于整体结构层数叠加方向温度均匀，采用已知的换热板组温度分布进行有限元验证；该专利提供了完整的印刷电路板式换热器性能设计及安全验证流程方法，极大的降低了设计成本，提升了设计效率。

Description

一种交叉流印刷电路板式换热器正向设计方法
技术领域
本发明专利涉及印刷电路板式换热器领域，具体涉及一种交叉流印刷电路板式换热器正向设计方法。
背景技术
印刷电路板式换热器是一种高效、紧凑式换热器，这种换热器具有换热效率高、耐高温、耐高压等优点，是深海浮式液化设备中主换热器的首选，也在核电、火电等领域有广泛的应用。
在此背景下，近些年来印刷电路板式换热器受到了广泛关注。但针对这种类型的换热器并没有优秀的流程化、模块化设计方法。目前在换热性能设计时没有形成准确的经验公式，且换热器芯体通道数量巨大，难以进行完整的有限元分析，因此通常参考板式换热器经验公式结合局部有限元法进行计算，但局部模型无法对流道均匀性进行评估，进而对由边缘流道冷却流量不足和局部过热等安全隐患无法预测。并且由于结构形式和加工工艺的特点，该方法原理上只能粗略获得如换热面积等换热器通用参数，无法直接得到换热器芯体结构如流道尺寸、间距等直接参数，结果上时常导致出口温度不准确，温度分布存在偏差。同时温度分布的不确定性也会导致结构设计保守，并使设计方通过使用更大尺寸的芯体和更厚的壁厚以保证换热性能和结构安全，最终使换热器失去了其高效、紧凑的优势，同时也造成了原材料和空间的浪费。
发明内容
为解决上述问题，在实现换热性能准确、结构安全的前提下，保证设计的经济性和高效性，本发明提供了一种交叉流印刷电路板式换热器正向设计方法，从芯体结构出发，创新性的提出了标准单元法，对换热板组进行单元化建模，并对其单元内的流动换热关系给出了数值解，实现温度分布计算的实时性、流程化、模块化。基于数值结果，该方法可以方便地从设计需求出发，根据输入参数对换热器性能及结构尺寸进行逐一对应设计，并验证其安全性。
为了达到上述目的，本发明的具体技术方案如下：
一种交叉流印刷电路板式换热器正向设计方法，采用标准单元法对换热器的温度分布进行数值计算；采用多孔介质法对换热器流动分布与压降进行计算；给出了换热器热应力最大值的具体位置及评估方法；
该方法设计对象换热器结构包括一次侧换热板1、二次侧换热板2、上盖板3、下盖板4、一次侧入口5、一次侧出口6、二次侧入口7及二次侧出口8；其中一次侧换热板1与二次侧换热板2交替布置，流动方向呈90°夹角，形成叠加方向周期性变化的温度分布，上盖板3和下盖板4作为边缘承压件，一次侧入口5及一次侧出口6结构相同，采用对称设计作为热侧压力边界，二次侧入口7及二次侧出口8结构相同，采用对称设计作为冷侧压力边界；一次侧换热板1与二次侧换热板2采用半圆形流道；
进行设计时，已知或能够初步评估以下参数：换热器所选结构材料，一次侧与二次侧流体工质设计压力P1和P2、一次侧与二次侧流体工质流量U1和U2、一次侧与二次侧流体工质进口温度T10和T20、一次侧与二次侧流体工质出口温度T1∞和T2∞，一次侧与二次侧换热板厚度D1与D2，换热器所处空间对一次侧换热板长度的限制尺寸L10，二次侧换热板长度的限制尺寸L20和高度的限制尺寸h0；
采用设计方法后得到的计算结果包括：一次侧与二次侧换热板温度分布数值结果，一次侧与二次侧工质在各自流道内的分布均匀性，一次侧与二次侧的工质压降，换热器整体热应力分布，一次侧与二次侧换热板流道长度L1和L2，一次侧与二次侧流道数N1和N2，一次侧与二次侧流道半径r1和r2，一次侧与二次侧流道间距d1和d2，上盖板、下盖板、一次侧与二次侧进口和一次侧与二次侧出口作为压力边界所需的结构厚度E；
设计方法在描述过程中使用的英文字母及后方数字编号具有以下规则：英文字母与其表达的定义为一一对应关系，对于流量U，英文字母后的1，2分别表示一次侧与二次侧，U1表示一次侧流量，U2表示二次侧流量；下角标0表示入口或初始值，则U10表示一次侧入口流量，下角标∞表示最终值或出口值，则U1∞表示一次侧出口流量；
对于下述步骤中的标准单元，其命名规则如下：编号No.ij表示一次侧第i条管道和二次侧第j条管道形成的标准单元，与进出口无关的参数：对流换热系数表示为Hij，与进出口有关的参数：一次侧入口流量表示为U1ij-1，一次侧出口流量为U1ij，二次侧入口流量为U2i-1j，二次侧出口流量为U2ij，编号U1ij-1也表示No.i j-1单元一次侧出口流量，但由于标准单元具备连续性，No.i j-1单元一次侧出口与No.ij单元一次侧入口为同一位置，因此两者使用同一编号；
所述正向设计方法具体包括以下步骤：
步骤1，选取标准单元及其尺寸：标准单元为具有一对交叉流动的一次侧流道和二次侧流道上下叠加所形成的单元结构，单个一次侧热板和二次侧热板形成的换热板组结构为多个标准单元正交组合而成；具体尺寸为一次侧流道半径r1＝α1*D1，二次侧流道半径r2＝α2*D2，取l1＝β1*r2作为标准单元一次侧流道截取长度，取l2＝β2*r1作为单元二次侧流道截取长度，纵向高度h＝D1+D2；其中α1、α2和β1、β2分别为标准单元尺寸系数，取值范围为1/2<α1，α2<1，2<β1，β2<5；
步骤2，建立标准单元流动换热模型：取No.11单元，根据空间限制的一次换热板长度L10，二次侧换热板长度L20，高度h0，且每个标准单元包括一个一次侧流道，一个二次侧流道，则计算初始一次侧流道数N10＝L10/l1，初始一次侧流道数N20＝L20/l2，初始换热板组层数为C0＝H0/h，由此得一次侧流量U110＝U1/(N10*C0)，二次侧流量U201＝U2/(N20*C0)，并且对于No.11单元，一次侧与二次侧入口工质温度、压力和物性即为换热器的入口工质温度、压力和物性，以此分别求出一次侧和二次侧入口位置雷诺数Re1和Re2，由于工质物性在单个标准单元内变化小，此处做均匀性假设，即入口位置雷诺数为标准单元内雷诺数，对后续无量纲参数做相同处理；同理求得一次侧和二次侧的普朗特数Pr1和Pr2，依据雷诺数与普朗特数数值范围，选取合适的对流换热系数公式，求得一次侧与二次侧的努塞尔数Nu1和Nu2，进而依据努塞尔数与对流换热系数公式求得一次侧与二次侧对流换热系数H111与H211；总体来看对流换热系数为关于标准单元入口温度Tij、标准单元流道半径r，标准单元流道截取长度l，入口流量Uij的函数式，即Hij＝f(Tij，r，l，Uij)；将对流换热系数带入热阻公式计算得到一次侧与二次侧热阻R111与R211，此处同样采用均匀性假设，并根据选取结构材料导热系数计算得到一次侧与二次侧流道间结构热阻Rm11；根据换热平衡关系得到一次侧与二次侧的壁面温度Tw111＝T110-(T110-T201)*R111/(R111+R211+Rm11)(结构热阻在上文提及),Tw211＝T201-(T201-T110)*R211/(R111+R211+Rm11)；根据对流换热公式Q＝hS(T-Tw)计算得到标准单元内一次侧与二次侧换热量，并根据均匀性假设采用入口温度下的定压比热容cp和能量公式Q＝UcpΔT得到工质流经标准单元后温度变化量，即可得到T111＝T110+ΔT111，T211＝T201+ΔT211；
步骤3，建立标准单元组成的换热板组平面流动换热模型：对于No.12单元，一次侧入口温度为No.11单元一次侧出口温度，二次侧入口温度为T201，也就是说，对于No.ij单元，一次侧入口温度为No.ij-1的一次侧出口温度，即T1ij-1，二次侧入口温度为No.i-1j的二次侧出口温度，即T2i-1j，对于整个标准单元组成的换热板组，任意已知入口温度情况即能求出出口温度；由于T110和T201已知，采用递归法得T1ij和T2ij，即任意流道数目下整个平面上每个标准单元一次侧与二次侧的出口温度，也就是换热板组完整的温度分布数值结果；接下来根据温度分布数值结果计算换热器参数流道长度L，流道数N，流道半径r，流道间距d，假设经过i个一次侧流道，j个二次侧流道后，两侧出口温度符合设计要求，那么得到一次侧流道数N1＝i，流道尺寸L1＝l1*i，二次侧流道数目N2＝j，流道尺寸L2＝l2*j,流道间距d1＝l2，d2＝l1；由于所有温度都为已知量的数值迭代算式，因此实时调整所有已知量，并根据当前流量下的温度分布实现对流道数、流道尺寸和层数的设计；
步骤4：芯体热应力迭代优化：根据已得到的换热板组温度分布，建立热固耦合模型对换热板组进行热应力验证，由于No.11为一次侧与二次侧温差最大处，因此需计算No.11至No.ij，共i*j个标准单元组成的结构的热应力，如出现热应力超标，在不改变标准单元尺寸的条件下单独改变一次侧与二次侧流道半径r1和r2，调整流道间壁厚，由此带来的流量变化通过调节换热板组层数C抵消；由于温度分布为已知量的代数式，因此能够即时完成计算并进行迭代优化，直至得到符合应力标准的一次侧与二次侧流道半径r1和r2；
步骤5：工质流动分布均匀性与压降计算：根据多孔介质原理，压降Δp＝av2+bv，式中a＝c2*0.5*ρ*Δn，其中c2为惯性阻力系数，ρ为流体密度，Δn为多孔介质区域厚度即流道尺寸L；式中b＝Δn*μ/α，其中1/α为粘性阻力系数，μ为动力粘度；由于流道分布均匀，根据计算已知的换热器结构，选取部分入口面积进行阻力验证，计算中需先假设多个速度v进行试算，得到Δp与速度v之间的多项式关系，通过线性回归求出惯性阻力系数和粘性阻力系数；其中速度v在取值时应采用与芯体入口速度相同数量级并将入口速度包含在取值范围内，以减少流动特性改变引起的误差；对堆芯一次侧与二次侧建立完整尺寸三维流体模型，三维流体模型阻力系数采用计算得到的惯性阻力系数和粘性阻力系数，并通过计算流体动力学方法CFD计算得到流体工质在流道内的分流均匀性和堆芯压降；
步骤6：换热器整体热应力计算：采用有限元法基于换热板组温度分布对换热器整体结构进行热应力计算，确定上盖板(3)和下盖板(4)以及一次侧入口(5)、一次侧出口(6)、二次侧入口(7)和二次侧出口(8)壁厚，确保整体结构热应力符合标准。
步骤1中α1＝2/3，α2＝2/3，β1＝3，β2＝3。
步骤5中部分入口面积为10％-50％入口面积。
和现有技术相比较，本发明具备如下优点：
1、本发明从流体传热原理出发，不借助传统换热器经验公式，以换热器结构为基础，通过数值计算得到换热器芯体温度场，避免了原理性误差，计算精度高，结果准确，减少了工程材料与空间的浪费，具备更强的经济型。
2、本发明采用多孔介质模型，解决了经验公式与局部有限元计算无法得到完整流场的问题，避免的流量分配不均和局部温度过高导致的安全隐患。
3、本发明设计流程短、效率高、计算流程化可操作性强，很大程度缩短了分析与设计的时间。
附图说明
图1为一种交叉流印刷电路板式换热器正向设计方法流程示意图；
图2为设计对象整体结构示意图；
图3为设计对象芯体结构示意图；
图4a为标准单元结构示意图，图4b为图4a的局部放大图；
图5为局部热应力示意图；
图6为流场均匀性示意图；
图7为整体应力示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
一种交叉流印刷电路板式换热器正向设计方法，采用标准单元法对换热器的温度分布进行数值计算；采用多孔介质法对换热器流动分布与压降进行计算；给出了换热器热应力最大值的具体位置及评估方法；
该方法设计对象换热器结构如图2所示，包括一次侧换热板1、二次侧换热板2、上盖板3、下盖板4、一次侧入口5、一次侧出口6、二次侧入口7及二次侧出口8；其中换热器热板排布如图3所示，一次侧换热板1与二次侧换热板2交替布置，流动方向呈90°夹角，形成叠加方向周期性变化的温度分布，上盖板3和下盖板4作为边缘承压件，一次侧入口5及一次侧出口6结构相同，采用对称设计作为热侧压力边界，二次侧入口7及二次侧出口8结构相同，采用对称设计作为冷侧压力边界；一次侧换热板1与二次侧换热板2采用半圆形流道；
进行设计时，已知或能够初步评估以下参数：换热器所选结构材料，一次侧与二次侧流体工质设计压力P1和P2、一次侧与二次侧流体工质流量U1和U2、一次侧与二次侧流体工质进口温度T10和T20、一次侧与二次侧流体工质出口温度T1∞和T2∞，一次侧与二次侧换热板厚度D1与D2，换热器所处空间对一次侧换热板长度的限制尺寸L10，二次侧换热板长度的限制尺寸L20和高度的限制尺寸h0；
采用设计方法后得到的计算结果包括：一次侧与二次侧换热板温度分布数值结果，一次侧与二次侧工质在各自流道内的分布均匀性，一次侧与二次侧的工质压降，换热器整体热应力分布，一次侧与二次侧换热板流道长度L1和L2，一次侧与二次侧流道数N1和N2，一次侧与二次侧流道半径r1和r2，一次侧与二次侧流道间距d1和d2，上盖板、下盖板、一次侧与二次侧进口和一次侧与二次侧出口作为压力边界所需的结构厚度E；
设计方法在描述过程中使用的英文字母及后方数字编号具有以下规则：英文字母与其表达的定义为一一对应关系，对于流量U，英文字母后的1，2分别表示一次侧与二次侧，U1表示一次侧流量，U2表示二次侧流量；下角标0表示入口或初始值，则U10表示一次侧入口流量，下角标∞表示最终值或出口值，则U1∞表示一次侧出口流量；
对于下述步骤中的标准单元，其命名规则如下：编号No.ij表示一次侧第i条管道和二次侧第j条管道形成的标准单元，与进出口无关的参数：对流换热系数表示为Hij，与进出口有关的参数：一次侧入口流量表示为U1ij-1，一次侧出口流量为U1ij，二次侧入口流量为U2i-1j，二次侧出口流量为U2ij，编号U1ij-1也表示No.i j-1单元一次侧出口流量，但由于标准单元具备连续性，No.i j-1单元一次侧出口与No.ij单元一次侧入口为同一位置，因此两者使用同一编号；
如图1所示，所述正向设计方法具体包括以下步骤：
步骤1，选取标准单元及其尺寸，标准单元结构如图4a和图4b所示：标准单元为具有一对交叉流动的一次侧流道和二次侧流道上下叠加所形成的单元结构，单个一次侧热板和二次侧热板形成的换热板组结构为多个标准单元正交组合而成；具体尺寸为一次侧流道半径r1＝α1*D1，二次侧流道半径r2＝α2*D2，取l1＝β1*r2作为标准单元一次侧流道截取长度，取l2＝β2*r1作为单元二次侧流道截取长度，纵向高度h＝D1+D2；其中α1、α2和β1、β2分别为标准单元尺寸系数，取值范围为1/2<α1，α2<1，2<β1，β2<5；推荐采用α1＝2/3，α2＝2/3，β1＝3，β2＝3。
步骤2，建立标准单元流动换热模型：取No.11单元，根据空间限制的一次换热板长度L10，二次侧换热板长度L20，高度h0，且每个标准单元包括一个一次侧流道，一个二次侧流道，则计算初始一次侧流道数N10＝L10/l1，初始一次侧流道数N20＝L20/l2，初始换热板组层数为C0＝H0/h，由此得一次侧流量U110＝U1/(N10*C0)，二次侧流量U201＝U2/(N20*C0)，并且对于No.11单元，一次侧与二次侧入口工质温度、压力和物性即为换热器的入口工质温度、压力和物性，以此分别求出一次侧和二次侧入口位置雷诺数Re1和Re2，由于工质物性在单个标准单元内变化小，此处做均匀性假设，即入口位置雷诺数为标准单元内雷诺数，对后续无量纲参数做相同处理；同理求得一次侧与二次侧的普朗特数Pr1和Pr2，依据雷诺数与普朗特数数值范围，选取合适的对流换热系数公式，求得一次侧与二次侧的努塞尔数Nu1和Nu2，进而依据努塞尔数与对流换热系数公式求得一次侧与二次侧对流换热系数H111与H211；总体来看对流换热系数为关于标准单元入口温度Tij、标准单元流道半径r，标准单元流道截取长度l，入口流量Uij的函数式，即Hij＝f(Tij，r，l，Uij)；将对流换热系数带入热阻公式计算得到一次侧与二次侧热阻R111与R211，此处同样采用均匀性假设，并根据选取结构材料导热系数计算得到一次侧与二次侧流道间结构热阻Rm11；根据换热平衡关系得到一次侧与二次侧的壁面温度Tw111＝T110-(T110-T201)*R111/(R111+R211+Rm11)(结构热阻在上文提及),Tw211＝T201-(T201-T110)*R211/(R111+R211+Rm11)；根据对流换热公式Q＝hS(T-Tw)计算得到标准单元内一次侧与二次侧换热量，并根据均匀性假设采用入口温度下的定压比热容cp和能量公式Q＝UcpΔT得到工质流经标准单元后温度变化量，即可得到T111＝T110+ΔT111，T211＝T201+ΔT211；
步骤3，建立标准单元组成的换热板组平面流动换热模型：对于No.12单元，一次侧入口温度为No.11单元一次侧出口温度，二次侧入口温度为T201，也就是说，对于No.ij单元，一次侧入口温度为No.ij-1的一次侧出口温度，即T1ij-1，二次侧入口温度为No.i-1j的二次侧出口温度，即T2i-1j，对于整个标准单元组成的换热板组，任意已知入口温度情况即能求出出口温度；由于T110和T201已知，采用递归法得T1ij和T2ij，即任意流道数目下整个平面上每个标准单元一次侧与二次侧的出口温度，也就是换热板组完整的温度分布数值结果；接下来根据温度分布数值结果计算换热器参数流道长度L，流道数N，流道半径r，流道间距d，假设经过i个一次侧流道，j个二次侧流道后，两侧出口温度符合设计要求，那么得到一次侧流道数N1＝i，流道尺寸L1＝l1*i，二次侧流道数目N2＝j，流道尺寸L2＝l2*j,流道间距d1＝l2，d2＝l1；由于所有温度都为已知量的数值迭代算式，因此实时调整所有已知量，并根据当前流量下的温度分布实现对流道数、流道尺寸和层数的设计，计算得到局部一次侧温度分布如下表1所示，与其对应标准单元编号的局部二次侧温度分布如下表2所示，左上角第一个数据为标准单元No.11温度，表格边框为管道方向示意；
表1局部一次侧温度分布
表2局部二次侧温度分布
步骤4：芯体热应力迭代优化：根据已得到的换热板组温度分布，建立热固耦合模型对换热板组进行热应力验证，由于No.11为一次侧与二次侧温差最大处，因此需计算No.11至No.ij，共i*j个标准单元组成的结构的热应力，此处选取i，j＝6，结果如图5所示，如出现热应力超标，在不改变标准单元尺寸的条件下单独改变一次侧与二次侧流道半径r1和r2，调整流道间壁厚，由此带来的流量变化通过调节换热板组层数C抵消；由于温度分布为已知量的代数式，因此能够即时完成计算并进行迭代优化，直至得到符合应力标准的一次侧与二次侧流道半径r1和r2；
步骤5：工质流动分布均匀性与压降计算：根据多孔介质原理，压降Δp＝av2+bv，式中a＝c2*0.5*ρ*Δn，其中c2为惯性阻力系数，ρ为流体密度，Δn为多孔介质区域厚度即流道尺寸L；式中b＝Δn*μ/α，其中1/α为粘性阻力系数，μ为动力粘度；由于流道分布均匀，根据计算已知的换热器结构，选取部分入口面积通常为10％-50％，进行阻力验证，计算中需先假设多个速度v进行试算，得到Δp与速度v之间的多项式关系，通过线性回归求出惯性阻力系数和粘性阻力系数；其中速度v在取值时应采用与芯体入口速度相同数量级并将入口速度包含在取值范围内，以减少流动特性改变引起的误差；对堆芯一次侧与二次侧建立完整尺寸三维流体模型，三维流体模型阻力系数采用计算得到的惯性阻力系数和粘性阻力系数，并通过计算流体动力学方法CFD计算得到流体工质在流道内的分流均匀性和堆芯压降，结果如图6所示；
步骤6：换热器整体热应力计算：采用有限元法基于换热板组温度分布对换热器整体结构进行热应力计算，结果如图7所示，并确定上盖板(3)和下盖板(4)以及一次侧入口(5)、一次侧出口(6)、二次侧入口(7)和二次侧出口(8)壁厚，确保整体结构热应力符合标准。

Claims (3)

1.一种交叉流印刷电路板式换热器正向设计方法，其特征在于：采用标准单元法对换热器的温度分布进行数值计算；采用多孔介质法对换热器流动分布与压降进行计算；给出了换热器热应力最大值的具体位置及评估方法；
该方法设计对象换热器结构包括一次侧换热板(1)、二次侧换热板(2)、上盖板(3)、下盖板(4)、一次侧入口(5)、一次侧出口(6)、二次侧入口(7)及二次侧出口(8)；其中一次侧换热板(1)与二次侧换热板(2)交替布置，流动方向呈90°夹角，形成叠加方向周期性变化的温度分布，上盖板(3)和下盖板(4)作为边缘承压件，一次侧入口(5)及一次侧出口(6)结构相同，采用对称设计作为热侧压力边界，二次侧入口(7)及二次侧出口(8)结构相同，采用对称设计作为冷侧压力边界；一次侧换热板(1)与二次侧换热板(2)采用半圆形流道；
进行设计时，已知或能够初步评估以下参数：换热器所选结构材料，一次侧与二次侧流体工质设计压力P1和P2、一次侧与二次侧流体工质流量U1和U2、一次侧与二次侧流体工质进口温度T10和T20、一次侧与二次侧流体工质出口温度T1∞和T2∞，一次侧与二次侧换热板厚度D1与D2，换热器所处空间对一次侧换热板长度的限制尺寸L10，二次侧换热板长度的限制尺寸L20和高度的限制尺寸h0；
采用设计方法后得到的计算结果包括：一次侧与二次侧换热板温度分布数值结果，一次侧与二次侧工质在各自流道内的分布均匀性，一次侧与二次侧的工质压降，换热器整体热应力分布，一次侧与二次侧换热板流道长度L1和L2，一次侧与二次侧流道数N1和N2，一次侧与二次侧流道半径r1和r2，一次侧与二次侧流道间距d1和d2，上盖板、下盖板、一次侧与二次侧进口和一次侧与二次侧出口作为压力边界所需的结构厚度E；
设计方法在描述过程中使用的英文字母及后方数字编号具有以下规则：英文字母与其表达的定义为一一对应关系，对于流量U，英文字母后的1，2分别表示一次侧与二次侧，U1表示一次侧流量，U2表示二次侧流量；下角标0表示入口或初始值，则U10表示一次侧入口流量，下角标∞表示最终值或出口值，则U1∞表示一次侧出口流量；
对于下述步骤中的标准单元，其命名规则如下：编号No.ij表示一次侧第i条管道和二次侧第j条管道形成的标准单元，与进出口无关的参数：对流换热系数表示为Hij，与进出口有关的参数：一次侧入口流量表示为U1ij-1，一次侧出口流量为U1ij，二次侧入口流量为U2i-1j，二次侧出口流量为U2ij，编号U1ij-1也表示No.i j-1单元一次侧出口流量，但由于标准单元具备连续性，No.i j-1单元一次侧出口与No.ij单元一次侧入口为同一位置，因此两者使用同一编号；
所述正向设计方法具体包括以下步骤：
步骤1，选取标准单元及其尺寸：标准单元为具有一对交叉流动的一次侧流道和二次侧流道上下叠加所形成的单元结构，单个一次侧热板和二次侧热板形成的换热板组结构为多个标准单元正交组合而成；具体尺寸为一次侧流道半径r1＝α1*D1，二次侧流道半径r2＝α2*D2，取l1＝β1*r2作为标准单元一次侧流道截取长度，取l2＝β2*r1作为单元二次侧流道截取长度，纵向高度h＝D1+D2；其中α1、α2和β1、β2分别为标准单元尺寸系数，取值范围为1/2<α1，α2<1，2<β1，β2<5；
步骤2，建立标准单元流动换热模型：取No.11单元，根据空间限制的一次换热板长度L10，二次侧换热板长度L20，高度h0，且每个标准单元包括一个一次侧流道，一个二次侧流道，则计算初始一次侧流道数N10＝L10/l1，初始一次侧流道数N20＝L20/l2，初始换热板组层数为C0＝H0/h，由此得一次侧流量U110＝U1/(N10*C0)，二次侧流量U201＝U2/(N20*C0)，并且对于No.11单元，一次侧与二次侧入口工质温度、压力和物性即为换热器的入口工质温度、压力和物性，以此分别求出一次侧和二次侧入口位置雷诺数Re1和Re2，由于工质物性在单个标准单元内变化小，此处做均匀性假设，即入口位置雷诺数为标准单元内雷诺数，对后续无量纲参数做相同处理；同理求得一次侧和二次侧的普朗特数Pr1和Pr2，依据雷诺数与普朗特数数值范围，选取合适的对流换热系数公式，求得一次侧与二次侧的努塞尔数Nu1和Nu2，进而依据努塞尔数与对流换热系数公式求得一次侧与二次侧对流换热系数H111与H211；总体来看对流换热系数为关于标准单元入口温度Tij、标准单元流道半径r，标准单元流道截取长度l，入口流量Uij的函数式，即Hij＝f(Tij，r，l，Uij)；将对流换热系数带入热阻公式计算得到一次侧与二次侧热阻R111与R211，此处同样采用均匀性假设，并根据选取结构材料导热系数计算得到一次侧与二次侧流道间结构热阻Rm11；根据换热平衡关系得到一次侧与二次侧的壁面温度Tw111＝T110-(T110-T201)*R111/(R111+R211+Rm11),Tw211＝T201-(T201-T110)*R211/(R111+R211+Rm11)；根据对流换热公式Q＝hS(T-Tw)计算得到标准单元内一次侧与二次侧换热量，并根据均匀性假设采用入口温度下的定压比热容cp和能量公式Q＝UcpΔT得到工质流经标准单元后温度变化量，即可得到T111＝T110+ΔT111，T211＝T201+ΔT211；
步骤3，建立标准单元组成的换热板组平面流动换热模型：对于No.12单元，一次侧入口温度为No.11单元一次侧出口温度，二次侧入口温度为T201，也就是说，对于No.ij单元，一次侧入口温度为No.ij-1的一次侧出口温度，即T1ij-1，二次侧入口温度为No.i-1j的二次侧出口温度，即T2i-1j，对于整个标准单元组成的换热板组，任意已知入口温度情况即能求出出口温度；由于T110和T201已知，采用递归法得T1ij和T2ij，即任意流道数目下整个平面上每个标准单元一次侧与二次侧的出口温度，也就是换热板组完整的温度分布数值结果；接下来根据温度分布数值结果计算换热器参数流道长度L，流道数N，流道半径r，流道间距d，假设经过i个一次侧流道，j个二次侧流道后，两侧出口温度符合设计要求，那么得到一次侧流道数N1＝i，流道尺寸L1＝l1*i，二次侧流道数目N2＝j，流道尺寸L2＝l2*j,流道间距d1＝l2，d2＝l1；由于所有温度都为已知量的数值迭代算式，因此实时调整所有已知量，并根据当前流量下的温度分布实现对流道数、流道尺寸和层数的设计；
步骤4：芯体热应力迭代优化：根据已得到的换热板组温度分布，建立热固耦合模型对换热板组进行热应力验证，由于No.11为一次侧与二次侧温差最大处，因此需计算No.11至No.ij，共i*j个标准单元组成的结构的热应力，如出现热应力超标，在不改变标准单元尺寸的条件下单独改变一次侧与二次侧流道半径r1和r2，调整流道间壁厚，由此带来的流量变化通过调节换热板组层数C抵消；由于温度分布为已知量的代数式，因此能够即时完成计算并进行迭代优化，直至得到符合应力标准的一次侧与二次侧流道半径r1和r2；
步骤5：工质流动分布均匀性与压降计算：根据多孔介质原理，压降Δp＝av2+bv，式中a＝c2*0.5*ρ*Δn，其中c2为惯性阻力系数，ρ为流体密度，Δn为多孔介质区域厚度即流道尺寸L；式中b＝Δn*μ/α，其中1/α为粘性阻力系数，μ为动力粘度；由于流道分布均匀，根据计算已知的换热器结构，选取部分入口面积进行阻力验证，计算中需先假设多个速度v进行试算，得到Δp与速度v之间的多项式关系，通过线性回归求出惯性阻力系数和粘性阻力系数；其中速度v在取值时应采用与芯体入口速度相同数量级并将入口速度包含在取值范围内，以减少流动特性改变引起的误差；对堆芯一次侧与二次侧建立完整尺寸三维流体模型，三维流体模型阻力系数采用计算得到的惯性阻力系数和粘性阻力系数，并通过计算流体动力学方法CFD计算得到流体工质在流道内的分流均匀性和堆芯压降；
步骤6：换热器整体热应力计算：采用有限元法基于换热板组温度分布对换热器整体结构进行热应力计算，确定上盖板(3)和下盖板(4)以及一次侧入口(5)、一次侧出口(6)、二次侧入口(7)和二次侧出口(8)壁厚，确保整体结构热应力符合标准。
2.根据权利要求1所述的一种交叉流印刷电路板式换热器正向设计方法，其特征在于：步骤1中α1＝2/3，α2＝2/3，β1＝3，β2＝3。
3.根据权利要求1所述的一种交叉流印刷电路板式换热器正向设计方法，其特征在于：步骤5中部分入口面积为10％-50％入口面积。
CN202111454588.4A 2021-12-01 2021-12-01 一种交叉流印刷电路板式换热器正向设计方法 Active CN114154262B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111454588.4A CN114154262B (zh) 2021-12-01 2021-12-01 一种交叉流印刷电路板式换热器正向设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111454588.4A CN114154262B (zh) 2021-12-01 2021-12-01 一种交叉流印刷电路板式换热器正向设计方法

Publications (2)

Publication Number Publication Date
CN114154262A CN114154262A (zh) 2022-03-08
CN114154262B true CN114154262B (zh) 2023-08-15

Family

ID=80455810

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111454588.4A Active CN114154262B (zh) 2021-12-01 2021-12-01 一种交叉流印刷电路板式换热器正向设计方法

Country Status (1)

Country Link
CN (1) CN114154262B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114595590B (zh) * 2022-03-22 2023-09-12 中国核动力研究设计院 再生式换热器偏离设计工况下换热功率分析方法及系统
CN115326005B (zh) * 2022-10-14 2022-12-13 中国核动力研究设计院 换热芯块微流道形变值的测量方法、装置、设备及介质
CN115791244B (zh) * 2023-02-06 2023-04-28 中国核动力研究设计院 一种模块式微通道紧凑换热实验本体、方法、设备及介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109443056A (zh) * 2018-09-14 2019-03-08 中国科学院工程热物理研究所 双面交错印刷电路板式换热板和换热器
CN110160380A (zh) * 2019-06-03 2019-08-23 中国矿业大学 一种宽通道板式换热器及换热器粒子群优化结构设计方法
JP2021173428A (ja) * 2020-04-20 2021-11-01 三菱アルミニウム株式会社 空調用熱交換器コアの評価方法と評価装置および評価プログラム

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109443056A (zh) * 2018-09-14 2019-03-08 中国科学院工程热物理研究所 双面交错印刷电路板式换热板和换热器
CN110160380A (zh) * 2019-06-03 2019-08-23 中国矿业大学 一种宽通道板式换热器及换热器粒子群优化结构设计方法
JP2021173428A (ja) * 2020-04-20 2021-11-01 三菱アルミニウム株式会社 空調用熱交換器コアの評価方法と評価装置および評価プログラム

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
胡彧;袁宏;王海松;成翔;刘晓波.一种超临界CO_2印刷电路板式换热器热工水力特性数值研究.科技视界.2020,(19),全文. *

Also Published As

Publication number Publication date
CN114154262A (zh) 2022-03-08

Similar Documents

Publication Publication Date Title
CN114154262B (zh) 一种交叉流印刷电路板式换热器正向设计方法
Prithiviraj et al. Three-dimensional numerical simulation of shell-and-tube heat exchangers. Part II: Heat transfer
Manglik et al. Heat transfer and pressure drop correlations for the rectangular offset strip fin compact heat exchanger
Matos et al. Optimally staggered finned circular and elliptic tubes in forced convection
Li et al. An approach based on the porous media model for numerical simulation of 3D finned-tubes heat exchanger
CN110610037B (zh) 逆流式回热器热力仿真方法
CN115659908B (zh) 一种印刷电路板换热器的三单元非平衡多孔介质方法
Gunnasegaran et al. The effect of geometrical parameters on heat transfer characteristics of compact heat exchanger with louvered fins
CN114444413A (zh) 一种板状燃料堆芯亚通道级三维热工水力分析方法
Huang et al. Multi-scale thermal analysis approach for the typical heat exchanger in automotive cooling systems
Du et al. Experimental study and numerical simulation of flow and heat transfer performance on an offset plate-fin heat exchanger
CN111027205B (zh) 一种液态金属冷却快堆堆芯盒间流数值模拟方法
CN102332041A (zh) 一种管带式散热器散热性能分析与结构设计系统
Xi et al. Numerical investigation on fluid-solid coupled heat transfer with variable properties in cross-wavy channels using half-wall thickness multi-periodic boundary conditions
Shahdad et al. Numerical analysis of the surface and geometry of plate fin heat exchangers for increasing heat transfer rate
Mustafa Maximization of heat transfer density rate from a single row of rhombic tubes cooled by forced convection based on constructal design
Yao et al. Numerical and experimental investigations of the heat transfer and fatigue life of a new W-type mini microchannel heat exchanger
Ge et al. Optimal arrangement design of a tube bundle in cross-flow using computational fluid dynamics and multi-objective genetic algorithm
Das et al. Multi-objective optimization of a novel butterfly-wing vortex generator fabricated in a rectangular microchannel based on CFD and NSGA-II genetic algorithm
Gao et al. Numerical and experimental investigation of additively manufactured shell-lattice copper heat exchanger
De Losier et al. The parametric study of an innovative offset strip-fin heat exchanger
Li et al. Study on iterative algorithm of full plate cross-flow and counter-flow printed circuit heat exchanger for fluoride-salt-cooled high-temperature advanced reactor
Duan et al. Topology optimization design of scramjet structures with forced convective heat transfer on unstructured meshes
Sheui et al. Numerical study of heat transfer in two-row heat exchangers having extended fin surfaces
CN112440416A (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