CN-122021179-A - Non-uniform TPMS lattice phase change heat sink design method based on topology optimization density mapping

CN122021179ACN 122021179 ACN122021179 ACN 122021179ACN-122021179-A

Abstract

The invention discloses a topological optimization density mapping-based heterogeneous TPMS lattice phase-change heat sink design method, which is formed by compounding a heterogeneous pore skeleton based on a three-period minimum curved surface and phase-change materials filled in the skeleton, wherein the three-period minimum curved surface unit porosities of the heterogeneous pore skeleton in the horizontal and vertical directions are both obtained by the heat sink topological optimization density mapping method, and the heat sink topological optimization density mapping method is also disclosed.

Inventors

LIU XIANGLEI
LI ZENGYUN
XUAN YIMIN

Assignees

南京航空航天大学

Dates

Publication Date: 20260512
Application Date: 20260225

Claims (10)

1. A design method of a non-uniform TPMS lattice phase-change heat sink based on topology optimization density mapping is characterized in that the non-uniform TPMS lattice phase-change heat sink is formed by compounding a non-uniform pore framework based on three-period minimum curved surfaces and phase-change materials filled in the framework; the three-period minimum curved surface unit porosities of the non-uniform pore skeleton in the horizontal direction and the vertical direction are both obtained by a heat sink topology optimization density mapping method; The heat sink topology optimization density mapping method comprises the following steps: S1, inputting a heat sink topological optimization model and an expression of a three-period minimum curved surface TPMS, and correspondingly designing a domain boundary and initial parameters; S2, generating a density mapping grid for a design domain of the heat sink topological optimization model by using a variable density method to obtain a heat sink topological optimization density field structure; s3, extracting the relative density distribution condition of the material distribution in the unit area according to the topological optimization density field; S4, establishing a mapping relation between the curved surface offset of the three-period minimum curved surface and the relative density in the unit area obtained in the S3, and constructing the three-period minimum curved surface of the design domain to obtain a TPMS mapping lattice phase change heat sink structure; S5, carrying out surface mesh division on the TPMS mapping lattice phase-change heat sink structure, and executing pre-processing operation of a model to obtain a model file of the TPMS mapping lattice phase-change heat sink structure with the surface mesh re-divided; s6, importing the model file into simulation software to carry out finite element body grid division; And S7, setting boundary conditions in simulation software, and performing simulation calculation and analysis to complete the optimal design of the TPMS mapping lattice phase change heat sink structure.
2. The method for designing the heterogeneous TPMS lattice phase-change heat sink based on the topology optimization density mapping, as claimed in claim 1, is characterized in that the framework is made of stainless steel, copper alloy, aluminum alloy, ceramic or resin, and the phase-change material is hexadecylamine or inorganic salt phase-change material.
3. The method for designing a heterogeneous TPMS lattice phase-change heat sink based on topology optimization density mapping as set forth in claim 1, wherein the level set control equation of the three-period minimum curved surface is F (X, Y, Z) =C, wherein X, Y, Z are lengths in three directions in a Cartesian coordinate system, and C is a curved surface offset.
4. The method for designing the heterogeneous TPMS lattice phase-change heat sink based on the topology optimization density map of claim 1, wherein the three-period minimum curved surface comprises PRIMITIVE configuration, gyroid configuration and Diamond configuration.
5. The method for designing the heterogeneous TPMS lattice phase-change heat sink based on the topology optimization density mapping of claim 4, wherein the control equation corresponding to the three-period minimum curved surface is: PRIMITIVE configuration: ; Gyroid configuration: ; Diamond configuration: ; where k=2pi/L.
6. The non-uniform TPMS lattice phase-change heat sink design method based on topology optimization density mapping of claim 1, wherein the heat sink topology optimization model in S1 is as follows: ; ; Wherein C t is thermal compliance, Γ q is the boundary of constant heat flux action in the second class of boundary conditions, In order to satisfy the temperature field of the homogeneous equation, gamma is an artificial density field, k is a heat conductivity coefficient, T is a temperature field, h is a heat convection coefficient, T ∞ is an ambient temperature, gamma is a gradient modulus of a pseudo density field, Q is a bulk heat source, omega is an optimal design domain, ρ e is a pseudo density of materials in the optimal design domain, and phi is an upper limit of a volume fraction of a heat sink structure, which is set to 0.2.
7. The method for designing a heterogeneous TPMS lattice phase-change heat sink based on topology optimization density mapping of claim 6, wherein the boundary range of the design domain in S1 is 。
8. The method for designing a heterogeneous TPMS lattice phase-change heat sink based on topology optimization density mapping of claim 1, wherein in S3, a design domain is discretized into a plurality of unit subareas, and the relative density of each unit subarea is counted respectively.
9. The method for designing the heterogeneous TPMS lattice phase-change heat sink based on the topology optimization density mapping of claim 1, wherein S4 establishes a mapping relation between the curved surface offset and the relative density in the unit area in a three-period minimum curved surface control equation, maps the topology optimization density field into the curved surface offset of the three-period minimum curved surface unit, and constructs the TPMS mapping lattice phase-change heat sink structure.
10. The non-uniform TPMS lattice phase-change heat sink design method based on topology optimization density mapping of claim 1, wherein in S7, a model is imported into simulation software to perform a phase-change heat sink heat transfer characteristic simulation experiment, and the structural phase-change material complete melting time is obtained through given boundary conditions, initial values and constraints.

Description

Non-uniform TPMS lattice phase change heat sink design method based on topology optimization density mapping Technical Field The invention relates to the technical field of phase-change heat storage, in particular to a non-uniform TPMS lattice phase-change heat sink design method based on topology optimization density mapping. Background With the rapid development of modern electronic components to high integration and high power density, power semiconductor devices are increasingly widely applied in modern electronic systems, and the problem of local extremely high heat flux density caused by the increasingly wide application of power semiconductor devices becomes a key bottleneck for restricting the performance and reliability of equipment. The traditional active heat dissipation modes such as air cooling and liquid cooling are difficult to be applied under the working condition of limited space or passive working condition, and the phase change heat storage technology provides an ideal passive heat management scheme for solving intermittent high thermal shock by virtue of the characteristic of high latent heat density and constant temperature heat storage of the phase change material. However, the inherent low thermal conductivity of common organic phase change materials severely restricts the rapid transport of heat, so that the thermal response inside the phase change heat sink is delayed, and the heat of a heat source cannot be timely led out. To solve this problem, the introduction of a high thermal conductivity metal framework into PCMs matrices to form composite materials has become the dominant technological route. The three-period minimum curved surface (TPMS) structure has zero average curvature, high specific surface area and full-communication pore characteristics, and is regarded as a potential reinforced heat transfer framework. However, the prior art is deficient in addressing the localized high heat flow conditions typical of electronic devices. Traditional TPMS skeleton design is mostly limited to uniform porosity or simple linear gradient distribution, and the homogenization design cannot be matched with the heat load characteristics of a bottom local heat source, so that the heat source area has overlarge diffusion heat resistance, and heat accumulation phenomenon occurs. Disclosure of Invention In order to overcome the defects of the background technology, the invention provides a non-uniform TPMS lattice phase-change heat sink design method based on topology optimization density mapping; the method takes a given design domain, constraint conditions and initial conditions as inputs, takes the minimum thermal flexibility as the topological optimization of an optimization target, carries out calculation iteration on the structure in the given design domain, finally obtains a reasonable material density field model, further uses a density mapping method to map by utilizing the grid relative density of a heat sink topological structure and the curved surface offset of a three-period minimum curved surface, and designs the TPMS mapping lattice phase change heat sink structure. According to the design method of the non-uniform TPMS lattice phase-change heat sink based on the topological optimization density mapping, the non-uniform TPMS lattice phase-change heat sink is formed by compounding a non-uniform pore framework based on a three-period minimum curved surface and phase-change materials filled in the framework; the three-period minimum curved surface unit porosities of the non-uniform pore skeleton in the horizontal direction and the vertical direction are both obtained by a heat sink topology optimization density mapping method; The heat sink topology optimization density mapping method comprises the following steps: S1, inputting a heat sink topological optimization model and an expression of a three-period minimum curved surface TPMS, and correspondingly designing a domain boundary and initial parameters; S2, generating a density mapping grid for a design domain of the heat sink topological optimization model by using a variable density method to obtain a heat sink topological optimization density field structure; s3, extracting the relative density distribution condition of the material distribution in the unit area according to the topological optimization density field; S4, establishing a mapping relation between the curved surface offset of the three-period minimum curved surface and the relative density in the unit area obtained in the S3, and constructing the three-period minimum curved surface of the design domain to obtain a TPMS mapping lattice phase change heat sink structure; S5, carrying out surface mesh division on the TPMS mapping lattice phase-change heat sink structure, and executing pre-processing operation of a model to obtain a model file of the TPMS mapping lattice phase-change heat sink structure with the surface mesh re-divided; s6, importing the model file int