CN113094944A - 一种微流道散热器及其细观尺度拓扑优化设计方法 - Google Patents
一种微流道散热器及其细观尺度拓扑优化设计方法 Download PDFInfo
- Publication number
- CN113094944A CN113094944A CN202110301573.8A CN202110301573A CN113094944A CN 113094944 A CN113094944 A CN 113094944A CN 202110301573 A CN202110301573 A CN 202110301573A CN 113094944 A CN113094944 A CN 113094944A
- Authority
- CN
- China
- Prior art keywords
- microscopic
- micro
- channel
- unit
- radiator
- 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
- 238000013461 design Methods 0.000 title claims abstract description 98
- 238000000034 method Methods 0.000 title claims abstract description 66
- 238000005457 optimization Methods 0.000 title claims abstract description 49
- 239000012530 fluid Substances 0.000 claims description 65
- 239000007787 solid Substances 0.000 claims description 18
- 239000000463 material Substances 0.000 claims description 16
- 238000010206 sensitivity analysis Methods 0.000 claims description 14
- 230000017525 heat dissipation Effects 0.000 claims description 11
- 238000012546 transfer Methods 0.000 claims description 9
- 239000012071 phase Substances 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 4
- 230000020169 heat generation Effects 0.000 claims description 4
- 230000005855 radiation Effects 0.000 claims description 3
- 239000007790 solid phase Substances 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 8
- 230000008878 coupling Effects 0.000 abstract description 4
- 238000010168 coupling process Methods 0.000 abstract description 4
- 238000005859 coupling reaction Methods 0.000 abstract description 4
- 230000008569 process Effects 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 3
- 238000001816 cooling Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 101100356682 Caenorhabditis elegans rho-1 gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000011160 research 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
- 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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于热流耦合拓扑优化设计相关技术领域，其公开了一种微流道散热器及其细观尺度拓扑优化设计方法，该方法包括以下步骤：首先，将微流道散热器进行局部特征提取以获得细观单元，并将所述细观单元作为设计对象，进而对所述细观单元进行有限元网格划分；接着，基于所述细观单元的设计域中的有限元网格，通过参数化水平集拓扑优化方法得到细观尺度下的微流道拓扑结构，并将细观尺度下的微流道拓扑结构进行阵列以得到微流道散热器。本发明将微流道散热器划分为多个细观尺度的单元，并以该细观单元作为对象进行拓扑优化设计，避免了宏观设计时为描述细小流道而划分大量有限元网格单元，节省了计算开销，提高了计算效率。
Description
技术领域
本发明属于热流耦合拓扑优化设计相关技术领域，更具体地，涉及一种微流道散热器及其细观尺度拓扑优化设计方法。
背景技术
液冷散热器在诸多工程领域中均有广泛应用，例如大功率电子设备、发动机、热交换器、航空航天飞行器等。例如，在航空航天方面，发动机再生冷却流道的设计本质上属于液冷散热器设计的范畴；在电子硬件领域，器件向着高热流密度、小型化、阵列化方向发展，这也对与之配套的散热系统设计提出了更高的要求。微流道散热器由于其比表面积高，可以实现数倍于传统流道散热器的换热系数，而且其体积小，可以进行多尺度、阵列化设计。基于经验的传统结构设计方法难以进一步提升散热效率，而借助结构优化技术来寻求结构某项或多项性能所对应的目标函数最优的设计方案已成为更为可行的设计手段。
结构优化技术中的拓扑优化设计方法由于其设计自由度高的特点，已被应用于散热器流道设计，流道散热拓扑优化目前主要包括两种方法，即变密度法和水平集法。变密度法求解稳定性较好，计算简单，但可能产生较多中间密度单元，导致其对结构边界描述不清晰。水平集方法可清晰描述结构边界，但传统水平集方法需要求解H-J PDE，计算较为复杂。参数化水平集方法可避免求解H-J PDE，但目前尚未应用到流道散热器设计领域；而且其研究多采用规则网格，难以适应带有复杂几何外形零件的设计需求，而且其多孔初始设计为人为选择，易导致初始设计依赖性。微流道散热器流道尺寸小，在宏观尺度进行整体设计需要划分大量网格单元，同时需要采用流道尺寸约束手段，计算开销较大。相应地，本领域存在着发展一种高效的微流道散热器拓扑优化方法的技术需求。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供了一种微流道散热器及其细观尺度拓扑优化设计方法，所述设计方法是一种热流耦合流道优化设计方法，其能够设计出微流道的微观构型，使得散热器的传热效率最大化，可以满足用于航空航天、电子工程中的散热器性能设计需求。
为实现上述目的，按照本发明的一个方面，提供了一种微流道散热器的细观尺度拓扑优化设计方法，所述方法包括以下步骤：
首先，将微流道散热器进行局部特征提取以获得细观单元，并将所述细观单元作为设计对象，进而对所述细观单元进行有限元网格划分；接着，基于所述细观单元的设计域中的有限网格单元得到细观尺度下的微流道拓扑结构，并将细观尺度下的微流道拓扑结构进行阵列以得到微流道散热器。
进一步地，该方法包括以下子步骤：
(1)将微流道散热器进行局部特征提取以获得细观单元，并将所述细观单元作为设计对象，进而对所述细观单元进行有限元网格划分；
(2)将紧支撑径向基函数分布于所述细观单元的设计域中的有限网格节点上，进而得到以水平集函数表示的散热器细观流道拓扑构型；
(3)基于该散热器细观流道拓扑构型得到有限元网格单元密度，继而得到细观流道拓扑构型对应散热器细观单元的物理模型；
(4)基于新得到的物理模型求解得到初始水平集函数后，逐步更新水平集函数的形状以确定最终的固体-流体分布，继而得到细观尺度下的微流道拓扑结构，并进行阵列以得到微流道散热器。
进一步地，步骤(1)中，将微流道散热器进行局部特征提取以获得细观单元，并将所述细观单元作为设计对象，根据所述细观单元中原有的微流道形式确定设计域及边界条件，继而根据设计域对微流道散热器的细观单元进行有限网格划分。
进一步地，步骤(2)中，将带有扩展系数的紧支撑径向基函数分布于所述细观单元的设计域中的有限网格节点上，并根据网格自适应地调整所述紧支撑径向基函数的支撑半径。
进一步地，基于径向基函数差值原理将所有基函数数值相加以得到参数化水平集函数；之后，基于所述参数化水平集函数及水平集理论描述所述设计域的固体及流体材料相，通过改变扩展系数使得所述参数化水平集函数具有不同形状，继而改变所述设计域中固体及流体材料分布，以得到以水平集函数表示的散热器细观流道拓扑构型。
进一步地，单个径向基函数的函数值表示为：
其中，
式中，x代表空间坐标；xI代表径向基函数中心点即控制点坐标；dmI用于控制支撑半径；dl为空间坐标为x的点与径向基函数中心点之间的欧氏距离；
支撑半径根据网格单元大小进行自适应调整，具体按照以下公式进行：
式中，xI代表径向基函数中心点坐标；xn代表相邻的基函数中心点坐标；N代表相邻基函数的数量；dmax为常数，通常在2-4范围内取值；n为相邻基函数的标号。
进一步地，步骤(3)中，基于所述散热器细观流道拓扑构型，采用Ersatz Material模型将水平集函数映射为连续变化的有限元网格单元密度，并采用带有温度方程的Navier-Stokes方程组构造设计域的不可压缩流体非等温强迫对流的物理模型，以表达散热流道中的流体传热行为；接着，根据得到的网格单元密度在各个网格单元中对Navier-Stokes方程组中的流体阻力、导热率和传热率进行固体相与流体相之间的线性插值，使得设计域中不同位置的网格单元分别具有固体或者流体的材料属性，继而得到流道拓扑构型对应散热器细观单元的物理模型。
进一步地，带有温度方程的Navier-Stokes方程组形式为：
式中，u为流体速度场；p为流体压力场；μ为流体动力粘度；ρ为流体密度；F为阻力项；c为流体热容；T为温度场；kf为流体导热率；
采用的插值公式为：
F＝-κu＝-[κmax+(κmin-κmax)·H(Φ)]u
式中，u为流体速度场；κmax为固体区域的阻力值，取一个较大的数，通常为105；κmin为流体区域阻力值，取一个较小的数，通常为0.01；H(Φ)为单元密度值，与水平集函数Ф相关；
得到的温度方程的插值公式为：
其中，β为产热系数。
进一步地，基于新得到的物理模型，采用非线性有限元方法对带有温度方程的Navier-Stokes方程组进行求解，以最大化散热效率作为目标函数，结合体积分数约束条件及有限元求解结果采用伴随方法进行敏度分析，并根据敏度分析结果和变密度拓扑优化理论获得初始水平集函数；接着，以该初始水平集函数为起点，根据有限元求解结果及敏度分析结果采用优化准则法逐步更新水平集函数的形状，以确定最终的固体-流体材料分布，继而得到细观尺度下的微流道拓扑结构，进而将细观尺度下的微流道拓扑结构进行阵列以得到完整的微流道散热器。按照本发明的另一个方面，提供了一种微流道散热器，所述微流道散热器是采用如上所述的微流道散热器的细观尺度拓扑优化设计方法设计而成的。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，本发明提供的微流道散热器及其细观尺度拓扑优化设计方法主要具有以下有益效果：
1.本发明将微流道散热器划分为多个细观尺度的单元，并以该细观单元作为对象进行拓扑优化设计，避免了宏观设计时为描述细小流道而划分大量有限元网格单元，节省了计算开销，提高了计算效率。
2.本发明采用变密度法进行进一步迭代以确定初始水平集函数，即拓扑优化中散热流道的初始拓扑构型，该初始设计为自动生成，避免了人工定义初始水平集函数对设计结果造成的人为影响，消除了设计的主观依懒性。
3.本发明采用参数化水平集方法描述散热流道拓扑构型，可将传统水平集中的Hamilton-Jacobi偏微分方程转化为常微分方程，不需要采用有限差分格式求解，因此不依赖结构化网络，同时采用了可以根据网格形状自适应调整支撑半径的紧支撑径向基函数，可以应用于具有非规则有限元网格的设计域。
4.本发明采用带有扩展系数的紧支撑径向函数对水平集函数进行拟合，在保证散热流道边界光滑性的同时，将标准水平集函数中的时间和空间两个耦合变量进行解耦，无需对水平集函数进行重新初始化，减少了计算时间消耗。
5.本发明的散热流道拓扑优化过程采用紧支撑径向基函数的扩展系数作为设计变量，可以采用优化准则法(OC)或者移动渐近线法(MMA)等成熟的梯度方法进行优化问题的求解，优化收敛过程稳定。
附图说明
图1是本发明提供的微流道散热器的细观尺度拓扑优化设计方法的流程示意图；
图2是本发明实施例的微流道散热器初始模型的示意图；
图3是根据图2中的微流道散热器初始模型的局部特征得到的设计域及边界条件示意图；
图4中的(a)、(b)分别是微流道散热器的初始拓扑构型及最终拓扑构型示意图；
图5是本发明涉及的迭代曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
请参阅图1，本发明提供的微流道散热器的细观尺度拓扑优化设计方法，所述方法主要包括以下步骤：
步骤一，将微流道散热器进行局部特征提取以获得细观单元，并将所述细观单元作为设计对象，根据所述细观单元中原有的微流道形式确定设计域及边界条件，继而根据设计域对微流道散热器的细观单元进行有限网格划分。
具体地，微流道散热器中的局部特征是指对人工设计的具有规则流道形状的微流道散热器进行局部可阵列特征的提取，继而确定细观单元形状、设计域及边界条件，该局部特征是指：对该特征形状执行平铺、圆周阵列等阵列化操作，可以重新得到原始的微流道散热器构型的特征。
步骤二，将带有扩展系数的紧支撑径向基函数(CS-RBFs)分布于所述细观单元的设计域中的有限网格节点上，并根据网格自适应地调整所述紧支撑径向基函数的支撑半径；接着，基于径向基函数差值原理将所有基函数数值相加以得到参数化水平集函数；之后，基于所述参数化水平集函数及水平集理论描述所述设计域的固体及流体材料相，通过改变扩展系数使得所述参数化水平集函数具有不同形状，继而改变所述设计域中固体及流体材料分布，以得到以水平集表示的散热器细观流道拓扑构型。
具体地，水平集描述形式是指构造一种标量函数场，成为水平集函数，函数值为零的点构成的集合即为零水平集，本实施方式以描述固-流边界，水平集函数值大于零的区域为流体域，小于零的则为固体域，形式为：
式中，Φ(x,t)为水平集函数值；x代表空间坐标；t代表时间，因为水平集函数是动态演化的；D代表设计域；Ω代表含边界的流体域。
紧支撑径向基函数(CS-RBFs)是指具有C2连续的紧支撑径向基函数，单个径向基函数的函数值表示为：
其中，
式中，x代表空间坐标；xI代表径向基函数中心点，即控制点坐标；dmI用于控制支撑半径。
由于基函数分布于有限元节点上，其中心坐标与有限元节点重合，为了适应非规则有限元网格，支撑半径可以根据网格单元大小进行自适应变化，具体实现方式为：
式中，xI代表径向基函数中心点坐标；xn代表相邻的基函数中心点坐标；N代表相邻基函数的数量；dmax为常数，通常在2-4范围内取值。其中，通过将每个基函数乘对应的扩展系数，并将所有带有扩展系数的基函数的函数值相加，以得到光滑的水平集函数数值，至此完成水平集函数的参数化。
步骤三，基于所述散热器细观流道拓扑构型，采用Ersatz Material模型将水平集函数映射为连续变化的有限元网格单元密度，并采用带有温度方程的Navier-Stokes方程组构造设计域的不可压缩流体非等温强迫对流的物理模型，以表达散热流道中的流体传热行为；接着，根据得到的网格单元密度在各个网格单元中对Navier-Stokes方程组中的流体阻力、导热率和传热率有关系数进行固体相与流体相之间的线性插值，使得细观设计域中不同位置的网格单元分别具有固体或者流体的材料属性，继而得到流道拓扑构型对应散热器细观单元的物理模型。
具体地，带有温度方程的Navier-Stokes方程组形式为：
式中，u为流体速度场；p为流体压力场；μ为流体动力粘度；ρ为流体密度；F为阻力项；c为流体热容；T为温度场；kf为流体导热率。其中，该方程组共含有三个方程，第一个方程为流体动量方程，第二个方程为流体不可压缩性方程，上述两个方程为Navier-Stokes方程组，描述了不可压缩流体的流动；第三个方程为温度方程，描述了流体的传热行为。
Ersatz Material模型是指在单个有限元网格单元中，将水平集函数竖直为正的区域所占单元面积(或者体积)的比例作为单元密度，由步骤二可以明显得出，流体区域内的单元密度为1，固体区域中的单元密度为0，边界处的单元密度介于0和1之间。此外，将该单元密度对流体和固体相关参数进行插值，即可得到流道拓扑结构对应的物理方程具体参数。
对于流体阻力，理论上固体中为无穷大，流体中为0，故插值形式表达为：
F＝-κu＝-[κmax+(κmin-κmax)·H(Φ)]u
式中，u为流体速度场；κmax为固体区域的阻力值，取一个较大的数，通常为105；κmin为流体区域阻力值，取一个较小的数，通常为0.01；H(Φ)为单元密度值，与水平集函数Ф相关。
对于导热和传热相关系数，固体和流体具有不同的热传导性质：
式中，ks表示固体导热率；Q表示热源强度，其他参数见步骤一中温度方程的描述，由此可以得到温度方程的插值形式：
其中，β为产热系数。
步骤四，基于新得到的物理模型，采用非线性有限元方法对带有温度方程的Navier-Stokes方程组进行求解，以最大化散热效率作为目标函数，结合约束条件及有限元求解结果采用伴随方法进行敏度分析，并根据敏度分析结果和变密度拓扑优化理论获得初始水平集函数；接着，以该初始水平集函数为起点，根据有限元求解结果及敏度分析结果采用优化准则法(OC)逐步更新水平集函数的形状，以确定最终的固体-流体材料分布，继而得到细观尺度下的微流道拓扑结构，进而将细观尺度下的微流道拓扑结构进行阵列以得到完整的微流道散热器。
具体地，最大化散热的目标函数为：
J＝∫D[1-H(Φ)]β(1-T)dΩ
约束通常为体积分数，完整的优化形式为：
Max J＝∫D[1-H(Φ)]β(1-T)dΩ
∫DH(Φ)dΩ≤Vmax
式中，Vmax是指最大体积分数。
敏度分析是基于上面的偏微分方程组约束，通过该约束方程组和目标函数构成拉格朗日函数，进而应用伴随变量方法进行敏度分析。由于温度方程和流体Navier-Stokes方程之间为弱耦合，即流动速度影响温度场进而温度场不影响流体速度，通过伴随方法得到的敏度形式为：
其中：
式中，
为Navier-Stokes方程在非线性有限元求解时自变量U的残差；
为温度方程在非线性有限元求解时自变量T的残差。
其中，通过变密度法获得初始拓扑的实施过程可分为三个子步骤，具体为：
4.1.1将设计域中的所有单元密度初始化为最大体积分数Vmax，并采用线性插值方式对方程系数进行插值得到具体的物理方程参数。
4.1.2将单元密度直接作为设计变量，采用非线性有限元物理进行方程的求解，利用得到的敏度分析公式进行敏度分析，然后采用优化准则法进一步迭代，从而更新密度场。
4.1.3将单元密度插值到节点，即参数化水平集函数的控制点处，将所有节点上的密度值同时减去所有节点密度值的中位数，并将每个节点上的密度值作为该节点上径向基函数的扩展系数，然后将基函数数值相加，即得到初始水平集函数，进而得到体积分数约为0.5的初始流道拓扑。
水平集拓扑优化的实施可分为三个子步骤，分别为：
4.2.1基于流道拓扑结构对应的物理模型，采用非线性有限元方法进行物理方程的求解，然后将有限元求解结果代入敏度分析公式进行敏度分析。
4.2.2利用敏度信息，根据参数化水平集拓扑优化原理将每个径向基函数的系数作为设计变量，采用优化准则法更新设计变量，继而更新水平集函数。
4.2.3判断迭代是否收敛，若未收敛则返回步骤4.2.1；若已经收敛，则根据结构的水平集函数得到细观流道拓扑结构，然后根据提取特征的方式将细观流道进行重新阵列化，从而得到完整的微流道散热器设计。
以下以一个实施例来对本发明进行进一步的详细说明。
实施例
本发明实施例提供了一种细观尺度散热微流道拓扑优化设计方法，且本实施例中以圆柱壁面微流道散热器为例说明本方法，参数不设量纲。该散热器为圆柱形，热源位于圆柱形中间的空腔，散热流道位于管壁内部，属于薄壁流道设计，因此采用二维代替三维。初始设计为平行流道，如图2所示，图中采用外表面黑色贴图代表流道，后面的图片亦采用该说明方式。假设管壁热流密度均匀，产热系数取β＝100，流体密度取为ρ＝1，动力粘度为μ＝1.3×10-3，热容为c＝4.2×103，流体导热率kf＝0.6，固体导热率ks＝400。流体由圆柱一端流至另一端，并在此过程中带走热量，流体流动由压力差驱动，该压力差取为Δp＝20。
实施例中，提取细观结构的方式如图2所示，该细观结构可以通过阵列方式还原流道结构。通过该细观结构提取，可以确定设计域以及边界条件信息，如图3所示，其中为了保证管道连通性，将非设计域设置为流道范围前后各5％长度的区域，该区域恒为流道。在流入边界处，温度设置为T＝0，因为轴向排列有20个细观单元，且假设压力梯度均匀，将入口压力取为pin＝1，出口压力取为pout＝0。上下两侧边界由于位于初始设计流道的中线处，设置边界条件为沿轴向滑移边界条件，即uy＝0，其他边界设置为无滑移壁面边界条件，即u＝0。最大体积分数设置为0.4，设置最大优化迭代步数为200步，迭代收敛条件为：目标函数值上一迭代步之差除以本迭代步目标函数值之商小于10-5，即终止迭代。
如图4所示，初始构型按照步骤(4.1.1)-(4.1.3)所述的方法产生，经过200步迭代以后收敛，由图5的曲线可以看出，迭代过程平稳，这种目标函数趋向于最大化，而体积满足体积分数约束。最后，将优化得出的细观结构进行重新阵列化，得到微流道散热器的最终设计。
本发明还提供了一种微流道散热器，所述微流道散热器是采用如上所述的微流道散热器的细观尺度拓扑优化设计方法设计而成的。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种微流道散热器的细观尺度拓扑优化设计方法，其特征在于，该方法包括以下步骤：
首先，将微流道散热器进行局部特征提取以获得细观单元，并将所述细观单元作为设计对象，进而对所述细观单元进行有限元网格划分；接着，基于所述细观单元的设计域中的有限元网格，通过参数化水平集拓扑优化方法得到细观尺度下的微流道拓扑结构，并将细观尺度下的微流道拓扑结构进行阵列以得到微流道散热器。
2.如权利要求1所述的微流道散热器的细观尺度拓扑优化设计方法，其特征在于：该方法包括以下子步骤：
(1)将微流道散热器进行局部特征提取以获得细观单元，并将所述细观单元作为设计对象，进而对所述细观单元进行有限元网格划分；
(2)将紧支撑径向基函数分布于所述细观单元的设计域中的有限网格节点上，进而得到以水平集函数表示的散热器细观流道拓扑构型；
(3)基于该散热器细观流道拓扑构型得到有限元网格单元密度，继而得到细观流道拓扑构型对应散热器细观单元的物理模型；
(4)基于新得到的物理模型求解得到初始水平集函数后，逐步更新水平集函数的形状以确定最终的固体-流体分布，继而得到细观尺度下的微流道拓扑结构，并进行阵列以得到微流道散热器。
3.如权利要求2所述的微流道散热器的细观尺度拓扑优化设计方法，其特征在于：步骤(1)中，将微流道散热器进行局部特征提取以获得细观单元，并将所述细观单元作为设计对象，根据所述细观单元中原有的微流道形式确定设计域及边界条件，继而根据设计域对微流道散热器的细观单元进行有限网格划分。
4.如权利要求2所述的微流道散热器的细观尺度拓扑优化设计方法，其特征在于：步骤(2)中，将带有扩展系数的紧支撑径向基函数分布于所述细观单元的设计域中的有限网格节点上，并根据网格自适应地调整所述紧支撑径向基函数的支撑半径。
5.如权利要求3所述的微流道散热器的细观尺度拓扑优化设计方法，其特征在于：基于径向基函数差值原理将所有基函数数值相加以得到参数化水平集函数；之后，基于所述参数化水平集函数及水平集理论描述所述设计域的固体及流体材料相，通过改变扩展系数使得所述参数化水平集函数具有不同形状，继而改变所述设计域中固体及流体材料分布，以得到以水平集函数表示的散热器细观流道拓扑构型。
6.如权利要求3所述的微流道散热器的细观尺度拓扑优化设计方法，其特征在于：单个径向基函数的函数值表示为：
其中，
式中，x代表空间坐标；xI代表径向基函数中心点即控制点坐标；dmI用于控制支撑半径；dl为空间坐标为x的点与径向基函数中心点之间的欧氏距离；
支撑半径根据网格单元大小进行自适应调整，具体按照以下公式进行：
式中，xI代表径向基函数中心点坐标；xn代表相邻的基函数中心点坐标；N代表相邻基函数的数量；dmax为常数，通常在2-4范围内取值；n为相邻基函数的标号。
7.如权利要求2-6任一项所述的微流道散热器的细观尺度拓扑优化设计方法，其特征在于：步骤(3)中，基于所述散热器细观流道拓扑构型，采用Ersatz Material模型将水平集函数映射为连续变化的有限元网格单元密度，并采用带有温度方程的Navier-Stokes方程组构造设计域的不可压缩流体非等温强迫对流的物理模型，以表达散热流道中的流体传热行为；接着，根据得到的网格单元密度在各个网格单元中对Navier-Stokes方程组中的流体阻力、导热率和传热率进行固体相与流体相之间的线性插值，使得设计域中不同位置的网格单元分别具有固体或者流体的材料属性，继而得到流道拓扑构型对应散热器细观单元的物理模型。
8.如权利要求7所述的微流道散热器的细观尺度拓扑优化设计方法，其特征在于：带有温度方程的Navier-Stokes方程组形式为：
式中，u为流体速度场；p为流体压力场；μ为流体动力粘度；ρ为流体密度；F为阻力项；c为流体热容；T为温度场；kf为流体导热率；
采用的插值公式为：
F＝-κu＝-[κmax+(κmin-κmax)·H(Φ)]u
式中，u为流体速度场；κmax为固体区域的阻力值，通常为105；κmin为流体区域阻力值，通常为0.01；H(Φ)为单元密度值，与水平集函数Ф相关；
得到的温度方程的插值公式为：
其中，β为产热系数。
9.如权利要求2-6任一项所述的微流道散热器的细观尺度拓扑优化设计方法，其特征在于：基于新得到的物理模型，采用非线性有限元方法对带有温度方程的Navier-Stokes方程组进行求解，以最大化散热效率作为目标函数，结合体积分数约束条件及有限元求解结果采用伴随方法进行敏度分析，并根据敏度分析结果和变密度拓扑优化理论获得初始水平集函数；接着，以该初始水平集函数为起点，根据有限元求解结果及敏度分析结果采用优化准则法逐步更新水平集函数的形状，以确定最终的固体-流体材料分布，继而得到细观尺度下的微流道拓扑结构，进而将细观尺度下的微流道拓扑结构进行阵列以得到完整的微流道散热器。
10.一种微流道散热器，其特征在于：所述微流道散热器是采用权利要求1-9任一项所述的微流道散热器的细观尺度拓扑优化设计方法设计得到的。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110301573.8A CN113094944B (zh) | 2021-03-22 | 2021-03-22 | 一种微流道散热器及其细观尺度拓扑优化设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110301573.8A CN113094944B (zh) | 2021-03-22 | 2021-03-22 | 一种微流道散热器及其细观尺度拓扑优化设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113094944A true CN113094944A (zh) | 2021-07-09 |
| CN113094944B CN113094944B (zh) | 2024-04-19 |
Family
ID=76668881
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110301573.8A Active CN113094944B (zh) | 2021-03-22 | 2021-03-22 | 一种微流道散热器及其细观尺度拓扑优化设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113094944B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114021497A (zh) * | 2021-11-03 | 2022-02-08 | 厦门大学 | 一种基于自动微分的可压缩湍流流体拓扑优化方法 |
| CN114386218A (zh) * | 2021-11-25 | 2022-04-22 | 西安电子科技大学 | 一种基于复杂功率分配的集成微流道的优化方法 |
| CN114564896A (zh) * | 2021-12-21 | 2022-05-31 | 西安交通大学 | 基于多尺度拓扑优化的强制对流微通道换热器设计方法 |
| CN117648778A (zh) * | 2024-01-29 | 2024-03-05 | 地平线(天津)科学技术应用研究有限公司 | 一种单管多翅片排管蒸发器的优化设计方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106547972A (zh) * | 2016-11-04 | 2017-03-29 | 中国科学院长春光学精密机械与物理研究所 | 基于参数化水平集方法的流体管道拓扑优化设计方法 |
| CN107391855A (zh) * | 2017-07-26 | 2017-11-24 | 华中科技大学 | 一种面向多种微观结构的材料结构一体化构建方法 |
| CN111709171A (zh) * | 2020-06-08 | 2020-09-25 | 西安交通大学 | 一种热流强耦合问题的等几何求解及散热拓扑生成方法 |
| CN112417692A (zh) * | 2020-11-24 | 2021-02-26 | 华东交通大学 | 基于载荷不确定性的材料结构多尺度拓扑优化设计方法 |
| US20210073428A1 (en) * | 2019-04-26 | 2021-03-11 | Dalian University Of Technology | Structure topology optimization method based on material-field reduced series expansion |
-
2021
- 2021-03-22 CN CN202110301573.8A patent/CN113094944B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106547972A (zh) * | 2016-11-04 | 2017-03-29 | 中国科学院长春光学精密机械与物理研究所 | 基于参数化水平集方法的流体管道拓扑优化设计方法 |
| CN107391855A (zh) * | 2017-07-26 | 2017-11-24 | 华中科技大学 | 一种面向多种微观结构的材料结构一体化构建方法 |
| US20210073428A1 (en) * | 2019-04-26 | 2021-03-11 | Dalian University Of Technology | Structure topology optimization method based on material-field reduced series expansion |
| CN111709171A (zh) * | 2020-06-08 | 2020-09-25 | 西安交通大学 | 一种热流强耦合问题的等几何求解及散热拓扑生成方法 |
| CN112417692A (zh) * | 2020-11-24 | 2021-02-26 | 华东交通大学 | 基于载荷不确定性的材料结构多尺度拓扑优化设计方法 |
Non-Patent Citations (1)
| Title |
|---|
| 马欣荣 等: "约束条件下微通道拓扑优化研究", 咸阳师范学院学报, no. 04, 25 July 2020 (2020-07-25) * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114021497A (zh) * | 2021-11-03 | 2022-02-08 | 厦门大学 | 一种基于自动微分的可压缩湍流流体拓扑优化方法 |
| CN114386218A (zh) * | 2021-11-25 | 2022-04-22 | 西安电子科技大学 | 一种基于复杂功率分配的集成微流道的优化方法 |
| CN114386218B (zh) * | 2021-11-25 | 2022-09-13 | 西安电子科技大学 | 一种基于复杂功率分配的集成微流道的优化方法 |
| CN114564896A (zh) * | 2021-12-21 | 2022-05-31 | 西安交通大学 | 基于多尺度拓扑优化的强制对流微通道换热器设计方法 |
| CN114564896B (zh) * | 2021-12-21 | 2024-02-09 | 西安交通大学 | 基于多尺度拓扑优化的强制对流微通道换热器设计方法 |
| CN117648778A (zh) * | 2024-01-29 | 2024-03-05 | 地平线(天津)科学技术应用研究有限公司 | 一种单管多翅片排管蒸发器的优化设计方法 |
| CN117648778B (zh) * | 2024-01-29 | 2024-04-23 | 地平线(天津)科学技术应用研究有限公司 | 一种单管多翅片排管蒸发器的优化设计方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113094944B (zh) | 2024-04-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113094944A (zh) | 一种微流道散热器及其细观尺度拓扑优化设计方法 | |
| CN112084591B (zh) | 一种基于三维拓扑优化的散热器冷却通道设计方法 | |
| CN111737835B (zh) | 基于三周期极小曲面的三维多孔散热结构的设计与优化方法 | |
| CN111709171B (zh) | 一种热流强耦合问题的等几何求解及散热拓扑生成方法 | |
| CN111709096B (zh) | 一种强化自然对流换热的异型翅片结构设计方法 | |
| WO2022000132A1 (zh) | 基于三周期极小曲面的三维多孔散热结构的设计与优化方法 | |
| CN111832204B (zh) | 一种热流耦合结构的非迭代式拓扑优化方法 | |
| Jahangiry et al. | Combination of Isogeometric analysis and level-set method in topology optimization of heat-conduction systems | |
| CN112784468B (zh) | 一种用于轻质防隔热承载结构的多尺度拓扑优化方法 | |
| Luo et al. | Multiscale simulation of a novel leaf-vein-inspired gradient porous wick structure | |
| CN115659908B (zh) | 一种印刷电路板换热器的三单元非平衡多孔介质方法 | |
| CN114925572A (zh) | 一种动力电池水冷流板的流道拓扑优化设计方法 | |
| Jiang et al. | Meshless optimization of triply periodic minimal surface based two-fluid heat exchanger | |
| Zhang et al. | A pseudo 3D cooling heat sink model designed by multi-objective topology optimization method | |
| Bianco et al. | Multi-material topology optimization of innovative microchannel heat sinks equipped with metal foams | |
| Iqbal et al. | An enriched finite element method for efficient solutions of transient heat diffusion problems with multiple heat sources | |
| Holka et al. | Density-based topology optimization of a surface cooler in turbulent flow using a continuous adjoint turbulence model | |
| CN113268910A (zh) | 一种重力驱动的自然对流异型热沉结构拓扑优化方法 | |
| Zhang et al. | Radiator optimization design for planar motors based on parametric components | |
| CN117317463B (zh) | 一种基于相变浆体的电池簇冷板流道拓扑优化设计方法 | |
| Yeranee et al. | Design of a serpentine cooling channel under turbulent flow using density-based topology optimization | |
| CN117725709A (zh) | 集成化发射天线散热流道拓扑优化设计方法 | |
| Zhao et al. | Feature-based topology optimization for channel cooling structure | |
| CN116361966A (zh) | 一种基于拓扑优化的脉动热管设计方法 | |
| CN114266184B (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 |