Liquid cooling plate flow channel design method and apparatus, storage medium, and electronic device

Abstract

The present application relates to the technical field of heat dissipation, and relates to a liquid cooling plate flow channel design method and apparatus, a storage medium, and an electronic device. The method is applied to an electronic device, and comprises: acquiring a first flow channel distribution of a liquid cooling plate, wherein a first quantity of target values of material densities corresponding to the first flow channel distribution is greater than a threshold, and the target values are between 0 and 1; acquiring an adjustment value of a target parameter, wherein the target parameter comprises at least one of a first parameter for adjusting fluid permeability and a second parameter for filtering the first flow channel distribution; and on the basis of the adjustment value of the target parameter, performing topology optimization on the first flow channel distribution to determine a target flow channel distribution of the liquid cooling plate. The heat dissipation effect of the liquid cooling plate designed according to the method of the present embodiments is closer to the expected target.

Description

Liquid cooling plate flow channel design method and device, storage medium and electronic equipment

[0001] The present application claims priority to the Chinese patent application No. 202411391578.4, filed on September 30, 2024, and entitled "Liquid cooling plate flow channel design method and device, storage medium and electronic equipment", the whole content of which is incorporated herein by reference. TECHNICAL FIELD

[0002] The present application relates to the field of heat dissipation technology, and more particularly, to a liquid cooling plate flow channel design method, a liquid cooling plate flow channel design device, a storage medium and an electronic equipment. BACKGROUND

[0003] With the increasing demand for the performance of electronic equipment, the heat dissipation problem of electronic equipment has gradually become one of the key concerns. The traditional heat dissipation technology is to design a flow channel on a liquid cooling plate, and the heat of the heat generating area of the electronic equipment is taken away by the flow of the fluid in the flow channel.

[0004] However, the heat dissipation effect of the liquid cooling plate obtained by machining the flow channel distribution obtained according to the existing liquid cooling plate flow channel design method cannot meet the expectation. Therefore, it is urgent to provide a new liquid cooling plate flow channel design method to improve the actual heat dissipation effect of the liquid cooling plate. SUMMARY

[0005] The present application provides a liquid cooling plate flow channel design method, a liquid cooling plate flow channel design device, a storage medium and an electronic equipment. The heat dissipation effect of the liquid cooling plate obtained by the method of the present application is closer to the expectation.

[0006] In a first aspect, the present application provides a liquid cooling plate flow channel design method applied to an electronic equipment, comprising:

[0007] obtaining a first flow channel distribution of a liquid cooling plate, wherein a first number of target values of material density corresponding to the first flow channel distribution is greater than a threshold value, and the target value is between 0 and 1;

[0008] obtaining an adjustment value of a target parameter, wherein the target parameter comprises at least one of a first parameter for adjusting fluid permeability and a second parameter for filtering the first flow channel distribution;

[0009] topologically optimizing the first flow channel distribution according to the adjustment value of the target parameter to determine a target flow channel distribution of the liquid cooling plate.

[0010] Optionally, the topologically optimizing the first flow channel distribution according to the adjustment value of the target parameter to determine the target flow channel distribution of the liquid cooling plate comprises:

[0011] updating the first flow channel distribution to a second flow channel distribution according to an adjustment value of the target parameter;

[0012] obtaining a second quantity of target values of the material density corresponding to the second flow channel distribution;

[0013] in a case where the second quantity is less than the threshold value, determining the target flow channel distribution according to the second flow channel distribution.

[0014] Optionally, the first parameter comprises a Darcy penalty function, and the second parameter comprises a hyperbolic tangent projection slope and / or a filtering radius.

[0015] Optionally, the second parameter comprises a hyperbolic tangent projection slope and a filtering radius, and the determining the target flow channel distribution of the liquid cooling plate according to the adjustment value of the target parameter and the topological optimization of the first flow channel distribution comprises:

[0016] topologically optimizing the first flow channel distribution according to an adjustment value of the Darcy penalty function to obtain a third flow channel distribution;

[0017] in a case where a quantity of target values of the material density corresponding to the third flow channel distribution is greater than the threshold value, obtaining a fourth flow channel distribution according to an adjustment value of the hyperbolic tangent projection slope;

[0018] in a case where a quantity of target values of the material density corresponding to the fourth flow channel distribution is greater than the threshold value, obtaining a fifth flow channel distribution according to an adjustment value of the filtering radius;

[0019] determining the target flow channel distribution according to the fifth flow channel distribution.

[0020] Optionally, the determining the target flow channel distribution of the liquid cooling plate according to the adjustment value of the target parameter and the topological optimization of the first flow channel distribution comprises:

[0021] obtaining an objective function, the objective function comprising any one of a heat transfer objective function, a flow objective function and a coupling objective function, the coupling objective function being an objective function obtained by coupling the heat transfer objective function and the flow objective function, the heat transfer objective function being determined based on at least one of an ambient temperature and a heat generation coefficient of a heat source corresponding to the liquid cooling plate, the flow objective function being determined based on at least one of a first fluid flow velocity corresponding to each grid in the flow channel distribution and a fluid domain proportion corresponding to each grid, the first fluid flow velocity comprising a horizontal velocity and a vertical velocity;

[0022] topologically optimizing the first flow channel distribution according to a value of the objective function and the adjustment value of the target parameter to determine the target flow channel distribution of the liquid cooling plate.

[0023] Optionally, the obtaining the target function comprises:

[0024] determining, according to at least one of the ambient temperature and the heat generation coefficient, temperature field information corresponding to the flow channel distribution, the temperature field information comprising a temperature at each grid in the flow channel distribution;

[0025] determining, according to a maximum temperature, a minimum temperature and a temperature of a set grid in the temperature field information, the target function;

[0026] Optionally, the obtaining the target function comprises:

[0027] determining, according to at least one of a first fluid flow velocity corresponding to each grid and a fluid domain proportion corresponding to each grid, flow velocity field information corresponding to the flow channel distribution, the flow velocity field information comprising a second fluid flow velocity of each grid in the flow channel distribution;

[0028] determining, according to the second fluid flow velocity, a third fluid flow velocity, the third fluid flow velocity representing a fluid flow velocity corresponding to the flow channel distribution;

[0029] determining, according to a preset maximum fluid flow velocity, a preset minimum fluid flow velocity and the third fluid flow velocity, the target function.

[0030] Optionally, the obtaining the target function comprises:

[0031] obtaining a sum of squares of the horizontal velocity and a sum of squares of the vertical velocity;

[0032] determining, according to a product of the sum of squares of the horizontal velocity and the sum of squares of the vertical velocity and the fluid domain proportion, the target function;

[0033] Optionally, the obtaining the target function comprises:

[0034] the obtaining the target function comprises:

[0035] obtaining an ambient temperature, the heat generation coefficient and a temperature of a set grid in the flow channel distribution;

[0036] determining, according to a product of a difference between the ambient temperature and the temperature of the set grid and the heat generation coefficient, the target function.

[0037] In a second aspect, the embodiment provides a liquid cooling plate flow channel design device, comprising:

[0038] a first obtaining module, configured to obtain a first flow channel distribution of a liquid cooling plate, a first number of target values of material density corresponding to the first flow channel distribution being greater than a threshold value, the target value being between 0 and 1;

[0039] The second acquisition module is configured to acquire an adjustment value of a target parameter, the target parameter including at least one of a first parameter for adjusting fluid permeability, or a second parameter for filtering the first flow channel distribution.

[0040] The determination module is configured to perform topology optimization on the first flow channel distribution according to the adjustment value of the target parameter, and determine a target flow channel distribution of the liquid cooling plate.

[0041] In a third aspect, an embodiment of the present application provides a storage medium having a computer program stored thereon, the computer program being configured to implement the method according to any one of the first aspect of the present application when executed by a processor.

[0042] In a fourth aspect, an embodiment of the present application provides an electronic device including a memory and a processor, the memory being configured to store computer instructions, and the processor being configured to invoke the computer instructions from the memory to execute the method according to any one of the first aspect of the present application.

[0043] The embodiments of the present application adjust the fluid permeability corresponding to the initial first flow channel distribution and / or filter the initial flow channel distribution, perform topology optimization on the region with the material density value between 0 and 1 in the initial flow channel distribution, so that the number of the material density values between 0 and 1 in the target flow channel distribution obtained after optimization meets the threshold value, thereby reducing the deviation between the design result and the processing result. In this way, the heat dissipation effect of the actual processed liquid cooling plate is closer to the expectation, and the actual heat dissipation effect is better than that of the liquid cooling plate in the prior art.

[0044] Other features and advantages of the present application will become apparent from the following detailed description of exemplary embodiments thereof, taken in conjunction with the accompanying drawings. BRIEF DESCRIPTION OF DRAWINGS

[0045] The accompanying drawings incorporated in and forming a part of the specification, illustrate embodiments of the present application and, together with the description, serve to explain the principles of the present application.

[0046] FIG. 1 shows a schematic diagram of a design domain according to an embodiment of the present application.

[0047] FIG. 2 shows a schematic flowchart of a liquid cooling plate flow channel design method according to an embodiment of the present application.

[0048] FIG. 3 shows a schematic diagram of a flow channel distribution according to an embodiment of the present application.

[0049] FIG. 4 shows a schematic flowchart of a liquid cooling plate flow channel design method according to another embodiment of the present application.

[0050] FIG. 5 shows a schematic block diagram of a liquid cooling plate flow channel design device according to an embodiment of the present application.

[0051] FIG. 6 shows a schematic block diagram of an electronic device according to an embodiment of the present application. DETAILED DESCRIPTION

[0052] Various exemplary embodiments of the present application will now be described in detail with reference to the accompanying drawings. Note that the relative arrangement, numerical expressions, and numerical values of components and steps set forth in these embodiments are illustrative only and do not limit the scope of the present application unless specifically stated otherwise.

[0053] The following description of at least one exemplary embodiment is merely illustrative in nature and is in no way intended to limit the scope of the application or its application or uses.

[0054] Techniques, methods, and apparatus known to those of ordinary skill in the relevant art can not be discussed in detail herein, but should be considered as part of the specification, where appropriate.

[0055] In all examples shown and discussed herein, any specific values should be interpreted as merely illustrative and not as a limitation. Thus, other examples of the exemplary embodiments can have different values.

[0056] Note that like numbers and letters refer to like elements throughout the several views of the drawings and that, as used herein, the term "and / or" includes any and all combinations of one or more of the associated listed items.

[0057] As the performance requirements of electronic devices continue to increase, the heat dissipation of electronic devices has gradually become one of the key concerns. The traditional heat dissipation technology is to design a flow channel on a liquid cooling plate, and the heat of the heat generating area of the electronic device is taken away by the flow of the fluid in the flow channel.

[0058] However, the heat dissipation effect of the liquid cooling plate obtained by machining the flow channel distribution obtained according to the existing liquid cooling plate flow channel design method cannot achieve the expected effect. Therefore, it is urgent to provide a new liquid cooling plate flow channel design method to improve the actual heat dissipation effect of the liquid cooling plate.

[0059] In order to facilitate understanding of the liquid cooling plate flow channel design method of the present application, the process of liquid cooling plate flow channel design will be briefly introduced first.

[0060] First, the design domain and boundary conditions of the liquid cooling plate are defined, and material properties are assigned to the fluid and solid. The design domain is the area where the flow channel of the liquid cooling plate is designed. The design domain can be a full solid area of a certain size, or a partial solid area and a partial fluid area, or a full fluid area. The fluid domain is the area where the fluid (such as coolant) can flow, and the solid domain is the area where the flow channel support column can be set. The design domain has a specific fluid domain and solid domain, which can be understood as needing to be further optimized and designed based on a specific flow channel. Material properties include density, heat source heat generation coefficient, viscosity, specific heat capacity, etc. The density of the fluid can be set to 1, and the density of the solid can be set to 0. The heat source heat generation coefficient is a physical quantity that describes the heat generation rate per unit area, which is usually used to evaluate the heat generation capacity of electronic devices in electronic equipment under certain conditions. The heat source heat generation coefficient of electronic devices of different power is different. The boundary conditions can include the velocity, pressure, and temperature of the fluid flow at the inlet and outlet of the flow channel.

[0061] Second, an appropriate number of grids are divided for the design domain to ensure that the grids of the flow channel and key areas are fine enough to capture details while avoiding excessive subdivision that wastes computing resources. FIG. 1 shows a simple schematic diagram of a design domain according to an embodiment of the present application. In FIG. 1, the design domain A is divided into 25 grids, and each grid corresponds to an area configured as a fluid domain (the area corresponding to "1" in the figure) or a solid domain (the area corresponding to "0" in the figure).

[0062] Then, the objective function and the constraint conditions are set. The objective function can be determined according to the performance indicators indicated in the design requirements. For example, the objective function can be determined according to at least one of the flow channel heat exchange efficiency, flow velocity, and pressure drop size indicated in the design requirements. The constraint conditions can include the volume rate. The volume rate is used to limit the flow channel volume to a certain proportion of the design domain.

[0063] Finally, topology optimization is automatically performed in the design domain according to the objective function, boundary conditions, and constraint conditions, and the density distribution of the solid domain and the fluid domain is adjusted to obtain an optimized flow channel distribution map. This process usually goes through multiple iterations, so that the performance results of the optimized flow channel distribution meet the design requirements.

[0064] Although the performance result corresponding to the flow channel distribution obtained by the algorithm simulation meets the design requirement, it does not mean that the performance result corresponding to the flow channel actually processed according to the flow channel distribution also meets the design requirement. The applicant found through analysis that compared with a conventional large-size liquid cooling plate, the thickness of the liquid cooling plate of an electronic device is very thin, for example, only about 0.4 mm. In such a thin liquid cooling plate, the number of material density values between 0 and 1 obtained by the flow channel distribution simulated by the topology optimization algorithm is usually relatively large. In the actual processing process, the regions with material density values between 0 and 1 will be processed as solid regions, resulting in a large difference between the actual processed flow channel and the simulated flow channel.

[0065] Based on this, the embodiment of the present application provides a liquid cooling plate flow channel design method applied to an electronic device. As shown in FIG. 2, the method can include steps S110 to S130.

[0066] In step S110, a first flow channel distribution of a liquid cooling plate is obtained, a first number of target values of material density corresponding to the first flow channel distribution is greater than a threshold value, and the target value is between 0 and 1.

[0067] In the embodiment, the first flow channel distribution can be obtained based on a full-solid initial design domain, or can be obtained based on a full-fluid initial design domain, or can be obtained based on a specific initial design domain, or can be directly obtained from an existing flow channel distribution library. The way of obtaining the first flow channel distribution is not limited here.

[0068] The first flow channel distribution represents the path of fluid flow. As an example, FIG. 3 shows a partial schematic diagram of a flow channel distribution, in which the region B with a material density value of 1 is a fluid region, forming a flow channel. The region with a material density value of 0 is a solid region. The region with a material density value between 0 and 1 is a porous medium region. The porous medium region will be processed as a solid region in the actual processing process, so it is necessary to identify the number of material density values between 0 and 1. In the embodiment, the material density value between 0 and 1 is taken as the target value, and the number of target values is first identified. In the case where the number of target values is greater than the threshold value, it is necessary to further optimize the current flow channel distribution to adjust the number of target values to be less than the threshold value.

[0069] The size of the threshold value can be determined according to the proportion of the region corresponding to the target value in the design domain. The specific value of the proportion is not limited in the embodiment, and can be determined according to the difference between the performance of the actual processed flow channel and the performance of the simulated flow channel. When the performance difference is required to be small, a lower proportion can be set.

[0070] In step S120, an adjustment value of a target parameter is obtained, the target parameter including at least one of a first parameter for adjusting fluid permeability and a second parameter for filtering the first flow channel distribution.

[0071] In step S130, the first flow channel distribution is topologically optimized according to the adjustment value of the target parameter to determine a target flow channel distribution of the liquid cooling plate.

[0072] The fluid permeability is an important parameter for describing the difficulty of fluid flow in a porous medium. In the case where the number of target values is greater than the threshold value, the fluid permeability can be adjusted by adjusting the value of the first parameter, so as to update the first flow channel distribution to the second flow channel distribution.

[0073] After obtaining the updated second flow channel distribution, the second flow channel distribution can be further updated by adjusting the second parameter.

[0074] Filtering the flow channel distribution generally means applying certain mathematical or physical rules to the flow channel distribution to smooth and simplify the shape of the flow channel, so that the material density value corresponding to the flow channel distribution approaches 0 or 1, thereby obtaining clear boundaries of fluid and solid domains, and facilitating processing. As an example, the second parameter can be a flow channel length filtering parameter. By limiting the length or width of the flow channel, flow channels that are difficult to process in the processing process are avoided. The second parameter can also be a shape filtering parameter. The geometry of the flow channel is optimized by shape rules, for example, avoiding sharp and pointed corners, so that the flow channel is smoother and more continuous. The second parameter can also be a topological rule filtering parameter, such as hyperbolic tangent projection slope, or filtering radius.

[0075] In the embodiment, the first parameter can include a Darcy penalty function, and the second parameter can include a hyperbolic tangent projection slope and / or a filtering radius. The adjustment range of the Darcy penalty function is 0-10. When the material density value is between 0 and 1, the intermediate density value is penalized, so that the intermediate density value gradually converges to 0 or 1. This can make the topological optimization model of the continuous variable better approximate the optimization model of the 0 or 1 discrete variable. The smaller the value of the parameter, the closer the intermediate density is to 0 or 1. The hyperbolic tangent projection slope is used to filter the overall design domain, and the adjustment range is 1-32. The larger the value of the parameter, the better the filtering effect. The adjustment range of the filtering radius is 0-1 m. By adjusting the value of the parameter, the flow channel distribution can be prevented from presenting a chessboard phenomenon, i.e., presenting more complete solid walls, resulting in disconnected flow channels.

[0076] In the case where the number of target values is greater than the threshold value, the value of the target parameter is adjusted, and the first flow channel distribution is topologically optimized according to the value of the target parameter to obtain the target flow channel distribution in the process of continuous iteration.

[0077] In this embodiment, step S130 can include steps S131-S133.

[0078] In step S131, the first flow channel distribution is updated to a second flow channel distribution according to the adjustment value of the target parameter.

[0079] In step S132, a second number of target values of material density corresponding to the second flow channel distribution is obtained.

[0080] In step S133, the target flow channel distribution is determined according to the second flow channel distribution in a case where the second number of target values of material density is less than a threshold value.

[0081] In this embodiment, the step of determining the target flow channel distribution according to the second flow channel distribution in step S133 can continue to optimize using any flow channel distribution optimization method to determine the target flow channel distribution, for example, using the topology optimization method to obtain the target flow channel distribution.

[0082] The embodiments of the present application adjust the fluid permeability corresponding to the initial first flow channel distribution and / or filter the initial flow channel distribution, so as to perform topology optimization on the region in the initial flow channel distribution where the material density value is between 0 and 1, so that the number of values of material density in the target flow channel distribution obtained after optimization is between 0 and 1 meets the threshold value, thereby reducing the deviation between the design result and the processing result. In this way, the heat dissipation effect of the actual processed liquid cooling plate is closer to the expectation, and the actual heat dissipation effect is better than that of the liquid cooling plate in the prior art.

[0083] In some embodiments, the order of adjustment of the target parameters can be set according to the influence degree of the target parameters on the flow channel topology optimization result. In this embodiment, the first parameter can be adjusted preferentially, and the first flow channel distribution is topology-optimized according to the adjustment value of the first parameter to obtain an updated flow channel distribution. In a case where the number of target values of material density corresponding to the updated flow channel distribution is greater than a threshold value, the second parameter is continued to be adjusted.

[0084] Specifically, step S130 can include steps S135-S138.

[0085] In step S135, the first flow channel distribution is topology-optimized according to the adjustment value of the Darcy penalty function to obtain a third flow channel distribution.

[0086] In step S136, in a case where the number of target values of material density corresponding to the third flow channel distribution is greater than a threshold value, a fourth flow channel distribution is obtained according to the adjustment value of the hyperbolic tangent projection slope.

[0087] In step S137, in a case where the number of target values of material density corresponding to the fourth flow channel distribution is greater than a threshold value, a fifth flow channel distribution is obtained according to the adjustment value of the filtering radius.

[0088] In step S138, the target flow channel distribution is determined according to the fifth flow channel distribution.

[0089] In this embodiment, step S138 can include: in a case where the number of target values of the material density corresponding to the fifth flow channel distribution is greater than the threshold value, adjusting the size of the grid corresponding to the fifth flow channel distribution to obtain a sixth flow channel distribution; and in a case where the number of target values of the material density corresponding to the sixth flow channel distribution is less than the threshold value, determining the target flow channel distribution according to the sixth flow channel distribution.

[0090] The embodiments of the present application set an adjustment order priority for the target parameters, and in a case where the flow channel distribution obtained by updating according to the adjustment value of the first parameter meets the requirements, the target parameters do not need to be continuously adjusted, so that the waste of resources can be reduced.

[0091] In some embodiments, the flow channel distribution meeting the design requirements cannot be obtained only by adjusting the value of the target parameter. This is because the accuracy of the value of the objective function used in the related flow channel design scheme is not enough to reflect the real performance. Therefore, the embodiments of the present application establish a new objective function according to the design requirements. Generally, the design requirements include the heat transfer efficiency of the flow channel, the fluid flow rate, etc. In this embodiment, a new heat transfer objective function and a flow objective function are established based on the design requirements of the heat transfer efficiency and the fluid flow rate of the flow channel, so as to improve the accuracy of the value of the objective function.

[0092] For the heat transfer objective function, the related art first simulates the temperature field information corresponding to the flow channel distribution by the finite element calculation method. The temperature field information includes the temperature at each grid in the flow channel distribution. The temperature corresponding to the simulated physical field information is directly used for the calculation of the heat transfer objective function. However, in the actual application of the liquid cooling plate flow channel, the temperature of the liquid cooling plate is also affected by the environmental temperature and the heat production coefficient of the heat source, and the value of the heat transfer objective function calculated directly by using the simulated temperature is obviously not accurate enough. Therefore, in this embodiment, the heat transfer objective function is determined based on at least one of the environmental temperature and the heat production coefficient of the heat source corresponding to the liquid cooling plate, so as to improve the accuracy of the value of the heat transfer objective function.

[0093] For the flow objective function, the related art directly uses the finite element calculation to obtain the velocity field information corresponding to the flow channel distribution to calculate the flow objective function. The velocity field information includes the horizontal velocity at each grid in the flow channel distribution. However, in the simulation calculation of the liquid cooling plate flow channel, the fluid velocity is a two-dimensional variable, including not only the horizontal velocity but also the vertical velocity. The flow objective function value calculated directly using the horizontal velocity is obviously not accurate enough. In addition, in the simulation calculation, each grid in the flow channel distribution is not necessarily all solid regions or fluid regions, but includes partial solid regions and partial fluid regions. Therefore, the flow velocity calculated using all fluid regions as the flow velocity corresponding to the grid including only partial fluid regions is also inaccurate. Therefore, in the embodiment, the flow objective function is determined based on at least one of the first fluid flow velocity corresponding to each grid in the flow channel distribution, including the horizontal velocity and the vertical velocity, and the fluid domain proportion corresponding to each grid, to improve the accuracy of the flow objective function value.

[0094] In the embodiment, step S130 can further include steps S210 to S220.

[0095] In step S210, a target function is obtained, the target function including any one of a heat transfer objective function, a flow objective function, and a coupling objective function, the coupling objective function being a target function obtained by coupling the heat transfer objective function and the flow objective function.

[0096] The coupling objective function can be a target function obtained by weighted summation of the heat transfer objective function and the flow objective function. The formula of the coupling objective function F can be: F=a*A+b*B

[0097] Wherein, A is the heat transfer objective function, a is the weight of the heat transfer objective function, B is the flow objective function, and b is the weight of the flow objective function.

[0098] Optionally, a specific weight can be pre-set for the heat transfer objective function and the flow objective function respectively according to experience to obtain the coupling objective function.

[0099] Optionally, in the topology optimization iteration process, the weight can be updated according to the quality of the current target function value after each iteration. For example, if the heat transfer objective function value has reached a good level in the current iteration, the weight of the heat transfer objective function can be reduced, and the weight of the flow objective function can be increased. Through this method, the algorithm can be more flexible to explore the solution space, avoiding concentrating on a certain specific target too early, thereby helping to find a flow channel distribution closer to the design requirements.

[0100] In step S220, a topology optimization is performed on the first flow channel distribution according to the value of the target function and the adjustment value of the target parameter to determine the target flow channel distribution of the liquid cooling plate.

[0101] In this embodiment, based on the adjustment value of the target parameter, the value of the first target function corresponding to the first flow channel distribution is calculated, and when the value of the first target function does not meet the preset target, the next iteration is performed to update the flow channel distribution. Then, based on the adjustment value of the target parameter, the second target function value corresponding to the updated flow channel distribution is calculated, and when the second target function value does not meet the preset target or the change of the second target function value compared with the value of the target function corresponding to the previous iteration is greater than a threshold value, the next iteration is continued until the value of the target function meets the preset target or the change of the value of the target function compared with the value of the target function corresponding to the previous iteration is less than the threshold value, and the flow channel distribution at this time is taken as the target flow channel distribution.

[0102] The target function established by the method of the embodiment of the application is more accurate in value when topology optimization is performed, and a flow channel distribution closer to the design requirement can be obtained.

[0103] In some embodiments, step S210 can include steps S211 to S212.

[0104] In step S211, temperature field information corresponding to the flow channel distribution is determined according to at least one of the ambient temperature and the heat source heat generation coefficient, and the temperature field information includes the temperature at each grid in the flow channel distribution.

[0105] In this embodiment, optionally, new temperature field information corresponding to the flow channel distribution can be determined according to at least one of the ambient temperature and the heat source heat generation coefficient. The new temperature field information is used to replace the temperature field information directly simulated by finite element calculation in the existing target function, so as to determine the target function.

[0106] Optionally, in step S212, the target function is determined according to the maximum temperature, the minimum temperature and the temperature of the set grid in the temperature field information obtained in step S211.

[0107] In some embodiments, step S210 can include steps S213 to S214.

[0108] In step S213, the ambient temperature, the heat source heat generation coefficient and the temperature of the set grid in the flow channel distribution are obtained.

[0109] In step S214, the target function is determined according to the product of the difference between the ambient temperature and the temperature of the set grid and the heat source heat generation coefficient.

[0110] In this embodiment, step S214 can take the product of the difference between the ambient temperature and the temperature of the set grid and the heat source heat generation coefficient as the temperature corresponding to each grid, and determine the target function according to the temperature.

[0111] As an example, the formula of the heat transfer objective function A can be: A = (Ta-Tamin) / (Tamax-Tamin)

[0112] wherein Ta = theta * (TQ-T), theta is a heat source heat generation coefficient, TQ is an ambient temperature, and T is a temperature corresponding to each grid obtained by finite element calculation.

[0113] In some embodiments, step S210 can include steps S215-S217.

[0114] Step S215: determining flow rate field information corresponding to the flow channel distribution according to at least one of the first fluid flow velocity corresponding to each grid and the fluid domain proportion corresponding to each grid, the flow rate field information including a second fluid flow velocity of each grid in the flow channel distribution.

[0115] Step S216: determining a third fluid flow velocity according to the second fluid flow velocity, the third fluid flow velocity representing a fluid flow velocity corresponding to the flow channel distribution.

[0116] In this embodiment, the second fluid flow velocity is a fluid flow velocity corresponding to each grid in the flow channel distribution, and the third fluid flow velocity represents a global fluid flow velocity corresponding to the flow channel distribution. As an example, the second fluid flow velocity can be summed to obtain the third fluid flow velocity. As another example, the second fluid flow velocity can be integrated to obtain the third fluid flow velocity.

[0117] Step S217: determining an objective function according to the preset maximum fluid flow velocity, the preset minimum fluid flow velocity, and the third fluid flow velocity.

[0118] In this embodiment, the maximum fluid flow velocity and the minimum fluid flow velocity can be pre-set according to experience for determining the objective function in step S216.

[0119] In some embodiments, step S210 can include steps S218-S219.

[0120] Step S218: obtaining a sum of squares of horizontal velocities and a sum of squares of vertical velocities corresponding to each grid.

[0121] Step S219: determining an objective function according to a product of the sum of squares of horizontal velocities and the sum of squares of vertical velocities and the fluid domain proportion.

[0122] In this embodiment, in step S219, the product of the sum of squares of horizontal velocities and the sum of squares of vertical velocities and the fluid domain proportion can be taken as a second fluid flow velocity corresponding to each grid. According to the second fluid flow velocity, the objective function is determined.

[0123] As an example, the formula of the flow target function B can be: B = (Fa-Famin) / (Famax-Famin)

[0124] wherein, r is the proportion of solid volume in the grid volume, (1-r) is the proportion of fluid volume in the grid volume, i.e. the proportion of fluid domain, u is the horizontal velocity, and v is the vertical velocity.

[0125] Embodiments of the present application also provide a liquid cooling plate flow channel design device. As shown in FIG. 5, the device 100 includes a first acquisition module, a second acquisition module, and a determination module.

[0126] The first acquisition module is configured to acquire a first flow channel distribution of a liquid cooling plate, wherein a first number of target values of material density corresponding to the first flow channel distribution is greater than a threshold value, and the target values are between 0 and 1.

[0127] The second acquisition module is configured to acquire an adjustment value of a target parameter, wherein the target parameter includes at least one of a first parameter for adjusting fluid permeability and a second parameter for filtering the first flow channel distribution.

[0128] The determination module is configured to perform topological optimization on the first flow channel distribution according to the adjustment value of the target parameter, and determine a target flow channel distribution of the liquid cooling plate.

[0129] Embodiments of the present application also provide a storage medium having a computer program stored thereon, wherein the computer program is executed by a processor to implement the method described in any one of the above method embodiments.

[0130] Embodiments of the present application also provide an electronic device, as shown in FIG. 6, the electronic device 200 includes a memory 210 and a processor 220,

[0131] The memory 210 is configured to store computer instructions, and the processor 220 is configured to call the computer instructions from the memory 210 to execute the method described in any one of the above method embodiments.

[0132] Each of the embodiments of the present application adopts a progressive manner for description, and the same or similar parts of each embodiment can be referred to each other. Each embodiment focuses on the difference from other embodiments. In particular, for the device and equipment embodiments, since they are basically similar to the method embodiments, the description is relatively simple, and the relevant parts can be referred to the part of the method embodiment.

[0133] The above-described embodiments of the application have several aspects, no single one of which is solely responsible for the application's desirable attributes. Without limiting the scope of the application as expressed by the claims that follow, some further embodiments make these aspects even more useful. Other embodiments fall within the scope of the appended claims.

[0134] The present application can be a system, a method, and / or a computer program product. The computer program product can include a computer readable storage medium (or media) having computer readable program instructions thereon for causing a processor to carry out aspects of the present application.

[0135] The computer readable storage medium can be a tangible device that can retain and store instructions for use by an instruction execution device. The computer readable storage medium can be, for example, but is not limited to, an electronic storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. A non-exhaustive list of more specific examples of the computer readable storage medium includes the following: a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanically encoded device such as punch-cards or punched tape, a

[0136] The computer readable program instructions described herein can be downloaded to respective computing / processing devices from a computer readable storage medium or to an external computer or external storage device via a network, for example, the Internet, a local area network, a wide area network and / or a wireless network. The network can comprise copper transmission cables, optical transmission fibers, wireless transmission, routers, firewalls, switches, gateway computers and / or edge servers. A network adapter card or network interface in each computing / processing device receives computer readable program instructions from the network and forwards the computer readable program instructions for storage in a computer readable storage medium within the respective computing / processing device.

[0137] Computer readable program instructions for carrying out operations of the present application can be assembler instructions, instruction-set-architecture (ISA) instructions, machine instructions, machine dependent instructions, microcode, firmware instructions, state-setting data, or either source code or object code written in any combination of one or more programming languages, including an object oriented programming language such as Smalltalk, C++ or the like and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The computer readable program instructions can execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer can be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection can be made to an external computer (for example, through the Internet using an Internet Service Provider). In some embodiments, electronic circuitry including, for example, programmable logic circuitry, field-programmable gate array (FPGA), or programmable logic array (PLA) can execute the computer readable program instructions by utilizing state information of the computer readable program instructions to personalize the electronic circuitry, in order to perform aspects of the present application.

[0138] The computer readable program instructions can also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide processes for implementing the functions / acts specified in the flowchart and / or block diagram block or blocks.

[0139] These computer readable program instructions can be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions / acts specified in the flowchart and / or block diagram block or blocks. These computer readable program instructions can also be stored in a computer readable storage medium that can include random access memory (RAM), read only memory (ROM), electrically erasable programmable read only memory (EEPROM), flash memory or other data storage device. When the computer readable program instructions are loaded into the computer and other programmable data processing apparatus, a series of operational steps are implemented that provide processes such that the computer or other programmable apparatus provide processes for implementing the functions / acts specified in the flowchart and / or block diagram block or blocks.

[0140] The computer readable program instructions can also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable data processing apparatus or other device to produce a computer implemented process such that the instructions which execute on the computer, other programmable data processing apparatus, or other device implement the functions / acts specified in the flowchart and / or block diagram block or blocks.

[0141] The computer readable program instructions can also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable data processing apparatus or other device to produce a computer implemented process such that the instructions which execute on the computer, other programmable data processing apparatus, or other device implement the functions / acts specified in the flowchart and / or block diagram block or blocks.

[0142] Embodiments of the present application have been described above, and the description is intended to be illustrative, and not restrictive, of the various embodiments of the present application. Many modifications and variations of the described embodiments of the present application are possible, given the benefit of the present disclosure, without departing from the scope and spirit of the described embodiments of the present application. The scope of the present application is defined by the appended claims.

Claims

1. A liquid-cooled plate flow channel design method, characterized in that, The application is applied to an electronic device, comprising: obtaining a first flow channel distribution of a liquid cooling plate, a first quantity of target values of material density corresponding to the first flow channel distribution being greater than a threshold value, the target values being between 0 and 1; obtaining an adjustment value of a target parameter, the target parameter comprising at least one of a first parameter for adjusting fluid permeability and a second parameter for filtering the first flow channel distribution; topologically optimizing the first flow channel distribution according to the adjustment value of the target parameter to determine a target flow channel distribution of the liquid cooling plate.

2. The method of claim 1, wherein, The topologically optimizing the first flow channel distribution according to the adjustment value of the target parameter to determine a target flow channel distribution of the liquid cooling plate comprises: updating the first flow channel distribution to a second flow channel distribution according to the adjustment value of the target parameter; obtaining a second quantity of target values of material density corresponding to the second flow channel distribution; in a case where the second quantity is less than the threshold value, determining the target flow channel distribution according to the second flow channel distribution.

3. The method of claim 1, wherein, The first parameter comprises a Darcy penalty function, and the second parameter comprises a hyperbolic tangent projection slope and / or a filtering radius.

4. The method of claim 3, wherein, The second parameter comprises a hyperbolic tangent projection slope and a filtering radius, and the topologically optimizing the first flow channel distribution according to the adjustment value of the target parameter to determine a target flow channel distribution of the liquid cooling plate comprises: topologically optimizing the first flow channel distribution according to the adjustment value of the Darcy penalty function to obtain a third flow channel distribution; in a case where a quantity of target values of material density corresponding to the third flow channel distribution is greater than the threshold value, obtaining a fourth flow channel distribution according to the adjustment value of the hyperbolic tangent projection slope; in a case where a quantity of target values of material density corresponding to the fourth flow channel distribution is greater than the threshold value, obtaining a fifth flow channel distribution according to the adjustment value of the filtering radius; determining the target flow channel distribution according to the fifth flow channel distribution.

5. The method according to any one of claims 1 to 4, characterized in that, The topologically optimizing the first flow channel distribution according to the adjustment value of the target parameter to determine a target flow channel distribution of the liquid cooling plate comprises: obtaining an objective function, the objective function comprising any one of a heat transfer objective function, a flow objective function and a coupling objective function, the coupling objective function being an objective function obtained by coupling the heat transfer objective function and the flow objective function, the heat transfer objective function being determined based on at least one of an ambient temperature and a heat generation coefficient of a heat source corresponding to the liquid cooling plate, the flow objective function being determined based on at least one of a first fluid flow velocity corresponding to each grid in the flow channel distribution and a fluid domain proportion corresponding to each grid, the first fluid flow velocity comprising a horizontal velocity and a vertical velocity; topologically optimizing the first flow channel distribution according to a value of the objective function and the adjustment value of the target parameter to determine a target flow channel distribution of the liquid cooling plate.

6. The method of claim 5, wherein, The obtaining the objective function comprises: determining temperature field information corresponding to the flow channel distribution according to at least one of the ambient temperature and the heat generation coefficient, the temperature field information comprising a temperature at each grid in the flow channel distribution; Determine the objective function according to the maximum temperature, the minimum temperature, and the temperature of the set grid in the temperature field information; Or, the obtaining of the objective function comprises: Determine the flow velocity field information corresponding to the flow channel distribution according to at least one of the first fluid flow velocity corresponding to each grid and the fluid domain proportion corresponding to each grid, wherein the flow velocity field information comprises a second fluid flow velocity of each grid in the flow channel distribution; Determine a third fluid flow velocity according to the second fluid flow velocity, wherein the third fluid flow velocity represents the fluid flow velocity corresponding to the flow channel distribution; Determine the objective function according to the preset maximum fluid flow velocity, the preset minimum fluid flow velocity, and the third fluid flow velocity.

7. The method of claim 5, wherein, The obtaining of the objective function comprises: Obtain the sum of the squares of the horizontal velocity and the sum of the squares of the vertical velocity; Determine the objective function according to the product of the sum of the squares of the horizontal velocity and the sum of the squares of the vertical velocity and the fluid domain proportion; Or, The obtaining of the objective function comprises: Obtain the ambient temperature, the heat generation coefficient, and the temperature of the set grid in the flow channel distribution; Determine the objective function according to the product of the difference between the ambient temperature and the temperature of the set grid and the heat generation coefficient.

8. A liquid cold plate flow channel design apparatus, characterized by, Comprise: The first acquisition module is used for acquiring a first flow channel distribution of a liquid cooling plate, wherein a first number of target values of material density corresponding to the first flow channel distribution is greater than a threshold value, and the target value is between 0 and 1; The second acquisition module is used for acquiring an adjustment value of a target parameter, wherein the target parameter comprises at least one of a first parameter for adjusting fluid permeability and a second parameter for filtering the first flow channel distribution; The determination module is used for performing topological optimization on the first flow channel distribution according to the adjustment value of the target parameter to determine a target flow channel distribution of the liquid cooling plate.

9. A storage medium, characterized by A computer program is stored thereon, and the computer program implements the method in any one of claims 1 to 7 when executed by a processor.

10. An electronic device, comprising: Comprise a memory and a processor, The memory is used for storing computer instructions, and the processor is used for calling the computer instructions from the memory to execute the method in any one of claims 1 to 7.

Images