CN112380652B - 一种冷却微通道的设计方法 - Google Patents

一种冷却微通道的设计方法 Download PDF

Info

Publication number
CN112380652B
CN112380652B CN202011271733.0A CN202011271733A CN112380652B CN 112380652 B CN112380652 B CN 112380652B CN 202011271733 A CN202011271733 A CN 202011271733A CN 112380652 B CN112380652 B CN 112380652B
Authority
CN
China
Prior art keywords
design
channel
microchannel
flow
micro
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
CN202011271733.0A
Other languages
English (en)
Other versions
CN112380652A (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.)
Sichuan Changhong Electric Co Ltd
Original Assignee
Sichuan Changhong Electric Co Ltd
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 Sichuan Changhong Electric Co Ltd filed Critical Sichuan Changhong Electric Co Ltd
Priority to CN202011271733.0A priority Critical patent/CN112380652B/zh
Publication of CN112380652A publication Critical patent/CN112380652A/zh
Application granted granted Critical
Publication of CN112380652B publication Critical patent/CN112380652B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Cooling Or The Like Of Electrical Apparatus (AREA)
  • Cooling Or The Like Of Semiconductors Or Solid State Devices (AREA)

Abstract

本发明公开了一种冷却微通道的设计方法，属于电子产品散热技术领域。该设计方法，包括以下步骤：步骤1：根据热源尺寸选定微通道设计尺寸，根据热源功率确定冷却工质流量；步骤2：根据选定的微通道设计尺寸确立拓扑优化的设计域；步骤3：设计域内采用动量方程描述流体流动，采用能量方程描述固体和流体的温度分布；步骤4：利用目标函数验证优化的设计域内流体及固体的分布，完成微通道装置的流道的三维拓扑结构设计。本发明采用拓扑结构优化的方式来设计微通道的流动结构，该结构具有流动阻力小和传热热阻小的双重优势。

Description

一种冷却微通道的设计方法
技术领域
本发明涉及电子产品散热技术领域，更具体的说是涉及一种冷却微通道的设计方法。
背景技术
随着社会的发展，高功率电子产品越来越多，集成的芯片热流密度越来越高。已有研究表明电子芯片的温度每升高10℃，故障率会增加一倍，温度过高甚至可能引起内部损坏，造成安全事故。因此电子产品散热愈加重要。
微通道是一种重要的冷却方式，具有散热能力强、体积小、噪音小等优点。微通道的截面多为矩形和圆形等形状，但传统微通道存在通道压力损失较大、温度分布不均匀等缺点，传统微通道一般不能达到最优的冷却效果。
发明内容
本发明的目的在于提供一种冷却微通道的设计方法，以期解决背景技术中存在的技术问题。
为了实现上述目的，本发明采用以下技术方案：
一种冷却微通道的设计方法，包括以下步骤：
步骤1：根据热源尺寸选定微通道设计尺寸，根据热源功率确定冷却工质流量；
步骤2：根据选定的微通道设计尺寸确立拓扑优化的设计域；
步骤3：设计域内采用动量方程描述流体流动，采用能量方程描述固体和流体的温度分布；
步骤4：利用目标函数验证优化的设计域内流体及固体的分布，完成微通道装置的流道的三维拓扑结构设计。
所述步骤2中：根据选定的微通道设计尺寸确立拓扑优化的设计域，包括：
(1)过滤设计域内设计变量；
(2)对设计域内相关参数进行插值。
所述过滤设计域内设计变量，包括：设计域内设计变量采用式1进行过滤：
Figure BDA0002777895410000021
其中r为过滤半径，
Figure BDA0002777895410000022
为过虑后的设计变量，γ为设计变量；
采用式2对过虑后的设计变量进一步过滤得到输出设计变量
Figure BDA0002777895410000023
Figure BDA0002777895410000024
其中，β和γβ分别为投影斜率和投影点。
所述对设计域内相关参数进行插值，包括：设计域内相关参数采用下式3进行插值：
Figure BDA0002777895410000025
αmin和αmax分别为反渗透率的最小值和最大值，其中，αmin取0，αmax取1×105，q为罚参数；
Figure BDA0002777895410000026
为输出设计变量，取值范围为0-1，0代表固体，1代表流体；/>
Figure BDA0002777895410000027
为有效导热系数，由下式4决定：
Figure BDA0002777895410000028
其中ks为固体导热系数，为kf流体导热系数，p为罚参数。
所述动量方程如下式5所示：
Figure BDA0002777895410000029
所述和能量方程如下式6表示：
Figure BDA00027778954100000210
其中，ρ为密度，u为速度场，p为压力，μ为流体动力粘度，α为多孔介质的反渗透率，cp为比热容，k为导热系数，T为温度。
所述利用目标函数验证优化的设计域内流体及固体的分布，完成微通道流道的三维拓扑结构设计，包括：优化设计目标函数考虑流动阻力和微通道热阻最小，如下式7所示：
Figure BDA0002777895410000031
其中ω1和ω2分别为
Figure BDA0002777895410000032
和/>
Figure BDA0002777895410000033
的权重，取值范围在0-1之间；/>
Figure BDA0002777895410000034
和
Figure BDA0002777895410000035
分别为归一化的流动耗散功和微通道热阻；/>
Figure BDA0002777895410000036
为微通道热阻最低时的流动耗散功，/>
Figure BDA0002777895410000037
为流动耗散功最小时的微通道热阻；
Figure BDA0002777895410000038
和/>
Figure BDA0002777895410000039
分别如下式8和9所示
Figure BDA00027778954100000310
Figure BDA00027778954100000311
其中Ω为设计域，Tave,source为热源平均温度，Tinlet为流体入口温度，Q为热源功率；
优化算法采用移动渐近线优化算法，通过上述算法获得设计域内流体及固体的分布，完成微通道流道的三维拓扑结构设计。
本发明的技术方案还包括将步骤4中得到的三维拓扑结构投影成二维平面，去除不规则边界，再将二维拉伸成三维拓扑结构。
进一步的，所述微通道设计域为矩形、圆形、椭圆形状。
进一步的，所述冷却工质包括去离子水、去离子水和乙二醇混合溶液、甲醇、丙酮中的一种。
进一步的，所述固体包括铝、铜、铁和PDMS中的一种。
进一步的，所述热源为芯片、激光器、LED等，热流密度范围为0-1000W/cm2。
本发明与现有技术相比具有的有益效果是：
本发明提供的设计方法能够提高微通道的换热能力，减小热阻，降低热源温度，同时降低流动阻力，减小功耗，效果远优于传统微通道。
附图说明
图1为利用本发明的一种冷却微通道的设计方法设计得到的微通道。
图2为一种普通微通道。
图3为压差随流量变化关系。
图4为温度随流量变化关系。
图5为温度随热流密度变化关系。
图6为温度随入口温度变化关系。
图7中7a、7b、7c为不同优化微通道流体域。
图8为利用本发明的一种冷却微通道的设计方法设计得的微通道的应用图。
具体实施方式
下面结合实施例对本发明作进一步的描述，所描述的实施例仅仅是本发明一部分实施例，并不是全部的实施例。基于本发明中的实施例，本领域的普通技术人员在没有做出创造性劳动前提下所获得的其他所用实施例，都属于本发明的保护范围。
实施例1:
一种冷却微通道的设计方法，包括以下步骤：
根据热源尺寸选定微通道设计尺寸，根据热源功率，确定冷却工质流量。根据选定的微通道设计尺寸确立拓扑优化的设计域，设计域内采用动量方程描述流体流动，采用能量方程描述固体和流体的温度分布。该设计域为三维几何结构。
动量方程和能量方程如下所示：
Figure BDA0002777895410000041
Figure BDA0002777895410000042
其中，ρ为密度，u为速度场，p为压力，μ为流体动力粘度，α为多孔介质的反渗透率，cp为比热容，k为导热系数，T为温度。
为了解决网格依赖问题，需要对设计变量进行过滤，设计域内设计变量γ采用下式进行过滤：
Figure BDA0002777895410000051
其中r为过滤半径，
Figure BDA0002777895410000052
为过虑后的设计变量。
为了减少流体于固体之间的灰色区域，采用下式进一步过滤
Figure BDA0002777895410000053
其中
Figure BDA0002777895410000054
为输出设计变量，β和γβ分别为投影斜率和投影点。
设计域内相关参数采用下式进行插值：
Figure BDA0002777895410000055
αmin和αmax分别为反渗透率的最小值和最大值，αmin一般取0，αmax一般取1×105，q为罚参数。
Figure BDA0002777895410000056
为输出设计变量，取值范围为0-1，0代表固体，1代表流体。/>
Figure BDA0002777895410000057
为有效导热系数，由下列方程式决定：
Figure BDA0002777895410000058
其中ks为固体导热系数，为kf流体导热系数，p为罚参数。
优化设计目标函数考虑流动阻力和微通道热阻最小，如下式所示
Figure BDA0002777895410000059
其中ω1和ω2分别为
Figure BDA00027778954100000510
和/>
Figure BDA00027778954100000511
的权重，取值范围在0-1之间。/>
Figure BDA00027778954100000512
和
Figure BDA00027778954100000513
分别为归一化的流动耗散功和微通道热阻。/>
Figure BDA00027778954100000514
为微通道热阻最低时的流动耗散功，/>
Figure BDA00027778954100000515
为流动耗散功最小时的微通道热阻。/>
Figure BDA00027778954100000516
和/>
Figure BDA00027778954100000517
分别如下式所示
Figure BDA00027778954100000518
Figure BDA00027778954100000519
其中Ω为设计域，Tave,source为热源平均温度，Tinlet为流体入口温度，Q为热源总功率。
通过上述算法获得设计域内流体及固体的分布，完成微通道流道的设计。
优化算法采用移动渐近线优化算法(MMA)，通过上述算法获得设计域内流体及固体的分布，完成微通道流道的三维拓扑结构设计。
考虑拓扑结构流固边界光滑性，可先将三维模型投影成二维平面，去除不规则边界，再将二维拉伸成三维拓扑结构。
所述的微通道设计域可以为矩形、圆形、椭圆等形状。所述的冷却工质包括去
具体参数入口流速为0.1m/s，微通道表面热流密度为7000W/m2，设计域长120mm，宽60mm。流体相关参数为kf＝0.61W/(m.K)，μ＝0.001Pa.s，ρf＝1000kg/m3，Cp,f＝4200J/(kg.K)，固体相关参数为ks＝237W/(m.K)，ρs＝2700kg/m3，Cp,s＝900J/(kg.K)，根据前述具体参数，利用本发明的拓扑优化方法获得如图1所示的优化微通道。同如图2所示的普通微通道相比较，压力最大可降低63Pa，最大温度可降低3℃，如图3、4、5、6所示。更换系数ω1和ω2的大小，可获得不同拓扑结构的微通道，不同的流体域如图7a、7b、7c所示。
图8是利用本发明的拓扑优化方法获得的微通道冷却热源装置示意图，在示意图中，微通道2紧贴热源1，两者之间一般涂敷有导热硅脂，并用螺钉固定。冷却工质从入口3进入，经过微通道后，从出口4流出，冷却工质一般包括去离子水、去离子水和乙二醇混合溶液、有机冷却工质、甲醇、丙酮等。固体结构可以为铝、铜、铁和PDMS等材料。热源可以为芯片、激光器和LED等。热流密度范围为0-1000W/cm2。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种冷却微通道的设计方法，其特征在于，包括以下步骤：
步骤1：根据热源尺寸选定微通道设计尺寸，根据热源功率确定冷却工质流量；微通道紧贴热源；
步骤2：根据选定的微通道设计尺寸确立拓扑优化的设计域；
步骤3：设计域内采用动量方程描述流体流动，采用能量方程描述固体和流体的温度分布；
步骤4：利用目标函数验证优化的设计域内流体及固体的分布，完成微通道装置的流道的三维拓扑结构设计；将三维拓扑结构投影成二维平面，去除不规则边界，再将二维拉伸成三维拓扑结构；
所述步骤2中：根据选定的微通道设计尺寸确立拓扑优化的设计域，包括：
(1)过滤设计域内设计变量；
对设计域内相关参数进行插值；
所述过滤设计域内设计变量，包括：设计域内设计变量采用式1进行过滤：
Figure QLYQS_1
其中r为过滤半径，
Figure QLYQS_2
为过虑后的设计变量，γ为设计变量；
采用式2对过虑后的设计变量进一步过滤得到输出设计变量
Figure QLYQS_3
Figure QLYQS_4
其中，β和γβ分别为投影斜率和投影点；
所述对设计域内相关参数进行插值，包括：设计域内相关参数采用下式3进行插值：
Figure QLYQS_5
αmin和αmax分别为反渗透率的最小值和最大值，其中，αmin取0，αmax取1×105，q为罚参数；
Figure QLYQS_6
为输出设计变量，取值范围为0-1，0代表固体，1代表流体；/>
Figure QLYQS_7
为有效导热系数，由下式4决定：
Figure QLYQS_8
其中ks为固体导热系数，为kf流体导热系数，p为罚参数。
2.根据权利要求1所述的一种冷却微通道的设计方法，其特征在于，所述动量方程如下式5所示：
Figure QLYQS_9
所述能量方程如下式6表示：
Figure QLYQS_10
其中，ρ为密度，u为速度场，p为压力，μ为流体动力粘度，α为多孔介质的反渗透率，cp为比热容，k为导热系数，T为温度。
3.根据权利要求1所述的一种冷却微通道的设计方法，其特征在于，所述利用目标函数验证优化的设计域内流体及固体的分布，完成微通道流道的三维拓扑结构设计，包括：优化设计目标函数考虑流动阻力和微通道热阻最小，如下式7所示：
Figure QLYQS_11
其中ω1和ω2分别为
Figure QLYQS_12
和/>
Figure QLYQS_13
的权重，取值范围在0-1之间；/>
Figure QLYQS_14
和/>
Figure QLYQS_15
Figure QLYQS_16
分别为归一化的流动耗散功和微通道热阻；/>
Figure QLYQS_17
为微通道热阻最低时的流动耗散功，/>
Figure QLYQS_18
为流动耗散功最小时的微通道热阻；
Figure QLYQS_19
和/>
Figure QLYQS_20
分别如下式8和9所示
Figure QLYQS_21
Figure QLYQS_22
其中Ω为设计域，Tave，source为热源平均温度，Tinlet为流体入口温度，Q为热源功率；
优化算法采用移动渐近线优化算法，通过上述算法获得设计域内流体及固体的分布，完成微通道流道的三维拓扑结构设计。
4.根据权利要求1所述的一种冷却微通道的设计方法，其特征在于，所述微通道设计域为矩形、圆形、椭圆形状。
5.根据权利要求1所述的一种冷却微通道的设计方法，其特征在于，所述冷却工质包括去离子水、去离子水和乙二醇混合溶液、甲醇、丙酮中的一种。
6.根据权利要求1所述的一种冷却微通道的设计方法，其特征在于，所述固体包括铝、铜、铁和PDMS中的一种。
7.根据权利要求1所述的一种冷却微通道的设计方法，其特征在于，所述热源为芯片、激光器、LED，热流密度范围为0-1000W/cm2。
CN202011271733.0A 2020-11-13 2020-11-13 一种冷却微通道的设计方法 Active CN112380652B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011271733.0A CN112380652B (zh) 2020-11-13 2020-11-13 一种冷却微通道的设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011271733.0A CN112380652B (zh) 2020-11-13 2020-11-13 一种冷却微通道的设计方法

Publications (2)

Publication Number Publication Date
CN112380652A CN112380652A (zh) 2021-02-19
CN112380652B true CN112380652B (zh) 2023-06-13

Family

ID=74582638

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011271733.0A Active CN112380652B (zh) 2020-11-13 2020-11-13 一种冷却微通道的设计方法

Country Status (1)

Country Link
CN (1) CN112380652B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115092896A (zh) * 2022-06-09 2022-09-23 蓬莱嘉信染料化工股份有限公司 一种微通道反应器制备亚硝酰硫酸的方法及其应用

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107122527A (zh) * 2017-04-10 2017-09-01 西安电子科技大学 一种基于拓扑优化的冷板流道设计方法
CN110944492A (zh) * 2019-12-03 2020-03-31 上海理工大学 液冷流道仿生优化设计方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11568098B2 (en) * 2018-05-16 2023-01-31 Autodesk, Inc Designing convective cooling channels
CN111832204B (zh) * 2020-07-02 2022-12-09 西安交通大学 一种热流耦合结构的非迭代式拓扑优化方法
CN112084591B (zh) * 2020-09-03 2022-09-13 西安电子科技大学 一种基于三维拓扑优化的散热器冷却通道设计方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107122527A (zh) * 2017-04-10 2017-09-01 西安电子科技大学 一种基于拓扑优化的冷板流道设计方法
CN110944492A (zh) * 2019-12-03 2020-03-31 上海理工大学 液冷流道仿生优化设计方法

Also Published As

Publication number Publication date
CN112380652A (zh) 2021-02-19

Similar Documents

Publication Publication Date Title
CN109860947B (zh) 一种动力电池包可变流道主动热管理控制方法及系统
Collins et al. A permeable-membrane microchannel heat sink made by additive manufacturing
CN107677152B (zh) 一种微通道液冷冷板
JP5754031B2 (ja) 高熱伝導率体を使用する電気エネルギー蓄積装置または電気化学エネルギー発生装置の温度を変更するシステム
CN113258097B (zh) 一种船用氢燃料电池冷却系统的控制方法
CN112380652B (zh) 一种冷却微通道的设计方法
CN103997880B (zh) 一种微通道散热器及由其组成的微机电产品散热系统装置
CN108847497B (zh) 一种车用燃料电池热管理系统
CN209626374U (zh) 液冷板、电池包及车辆
CN206775904U (zh) 一种采用复合相变材料和液态金属散热的热控装置
CN113094944B (zh) 一种微流道散热器及其细观尺度拓扑优化设计方法
CN113873849B (zh) 一种自适应调节半浸没式液冷散热腔体、循环系统及应用
CN114626260A (zh) 一种逆流换热器流道的拓扑优化设计方法
CN108375130A (zh) 一种变频空调驱动模块的散热装置及其控制方法
CN111463179A (zh) 基于弹性湍流的超低流阻微通道液冷换热器及其制造方法
CN115963903A (zh) 一种浸没液冷服务器装置及控制方法
Wang et al. A review: The development of crucial solar systems and corresponding cooling technologies
CN108418545A (zh) 一种加入多孔传热表面的微喷射流冷却板及其制造方法
CN106197081B (zh) 用于气井开采的驱动装置的空冷机
CN108627044A (zh) 一种用于超临界二氧化碳回热器变截面机翼型高效换热通道设计方法
CN112930082B (zh) 吸油烟机控制板散热方法
CN208780012U (zh) 一种用于超临界二氧化碳回热器变截面机翼型高效换热通道
CN204830957U (zh) 大齿距大波浪型流道抗污染换热器
CN203859726U (zh) 一种用于太阳能电池的纳米流体石墨微通道散热器
CN117317463B (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