CN114239435B - Numerical calculation method for three-dimensional flat plate water heat pipe - Google Patents

Numerical calculation method for three-dimensional flat plate water heat pipe Download PDF

Info

Publication number
CN114239435B
CN114239435B CN202111564511.2A CN202111564511A CN114239435B CN 114239435 B CN114239435 B CN 114239435B CN 202111564511 A CN202111564511 A CN 202111564511A CN 114239435 B CN114239435 B CN 114239435B
Authority
CN
China
Prior art keywords
interface
liquid
temperature
pressure
vapor
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
CN202111564511.2A
Other languages
Chinese (zh)
Other versions
CN114239435A (en
Inventor
王明军
穆德芳
章静
秋穗正
苏光辉
田文喜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN202111564511.2A priority Critical patent/CN114239435B/en
Publication of CN114239435A publication Critical patent/CN114239435A/en
Application granted granted Critical
Publication of CN114239435B publication Critical patent/CN114239435B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)

Abstract

The invention discloses a numerical calculation method of a three-dimensional flat plate water heat pipe, which comprises the following steps: establishing a three-dimensional geometric model of the flat water heat pipe according to the real design parameters; grid division is carried out on the flat plate water heat pipe geometric model by using grid division software; the grid is led into computational fluid dynamics software, and the boundary condition of the heat pipe is set according to the real situation; setting initial conditions required by transient calculation; programming the UDF to calculate evaporation/condensation mass flow rate, temperature and pressure at the wick-vapor interface, and liquid and vapor densities per time step; and calculating the flowing heat transfer in the three-dimensional flat plate hydrothermal pipe by using computational fluid dynamics software. The method can calculate the operation of the heat pipe with unsteady high heat flux density by using a computational fluid mechanics means, has better stability and convergence in calculation, and provides reference for the design of the newly-appearing heat pipe.

Description

Numerical calculation method for three-dimensional flat plate water heat pipe
Technical Field
The invention belongs to the technical field of thermal hydraulic numerical calculation, and particularly relates to a numerical calculation method of a three-dimensional flat plate water heat pipe.
Background
A heat pipe is a passive heat transfer device that has found wide use in thermal systems because it provides efficient heat transfer with minimal loss. With the chip-level heat flux reaching 200-300W/cm 2 One of the driving forces in heat pipe development is the need to reliably manage the heat dissipation of electronic systems while maintaining the equipment temperature to specification requirements. The passive nature of its operation and the ability to transfer heat over reasonable distances with only small changes in surface temperature make it particularly attractive in electronic cooling and packaging designs. Heat pipes are also widely used as heat sinks. It uses latent heat exchange to realizeThe heat transfer rate is now high. The increasing power density of electronic chips requires that the performance of these devices be optimized in order to be able to effectively remove heat from the electronic chip while limiting the temperature difference between the chip and the environment. The capillary porous media (wick) employed in the device is a major factor affecting the heat transfer efficiency of the heat pipe and vapor chamber. The wick also determines the maximum heat transfer capacity. The research on the phase-change heat and mass transfer characteristics in the liquid suction core is helpful for optimizing the design of the liquid suction core and improving the performance of the phase-change cooling device.
During transient operation, the system pressure in the heat pipe changes as the liquid/vapor interface in the wick evaporates and condenses. Since the liquid-vapor density is relatively large, a small change in the phase change rate causes a large change in the system pressure. The system pressure in turn changes the interface pressure and thus the saturation temperature by the clausius-claperton equation. The evaporation and condensation rates depend on the interfacial evaporation resistance, which depends on both the interfacial pressure and the volumetric pressure. In addition, the system pressure also directly changes the vapor density through the ideal gas law. These non-linear relationships can cause difficulties in convergence of the numerical formats, especially at high heat additions.
Disclosure of Invention
The invention aims to provide a numerical calculation method of a flat plate water heat pipe, which can calculate the operation of an unsteady high heat flux heat pipe by using a computational fluid mechanical means, has good stability and convergence in calculation, and provides reference for newly-appearing heat pipe design.
In order to achieve the above purpose, the invention adopts the following technical scheme:
a numerical calculation method of a three-dimensional flat plate water heat pipe is characterized by comprising the following steps of: the method comprises the following steps:
step 1: establishing a software solidworks through a geometric model to establish a three-dimensional geometric model of the flat water heat pipe according to a real design, wherein the three-dimensional geometric model comprises a pipe wall, a liquid suction core, a vapor chamber, an evaporation section and a condensation section;
step 2: performing grid division on the flat plate water heat pipe geometric model obtained in the step 1 by using grid division software;
step 3: the grids divided in the step 2 are led into computational fluid dynamics software, an evaporation section and a condensation section of the heat pipe are respectively set as a heat flow boundary condition and a convection boundary condition, other wall surfaces are set as heat insulation wall surfaces, the heat flow density heated by the evaporation section, the heat exchange coefficient of the condensation section and the temperature of a coolant are set according to the actual situation, and the porosity of the liquid absorption core is set;
step 4: setting initial conditions required for transient calculation: setting the initial steam mass as 0, initializing a temperature field, and enabling the working pressure of an initial steam area to be the saturation pressure at the current temperature;
step 5: the user-defined program UDF was used to calculate the gas-liquid interface mass flow rate, the temperature and pressure at the wick-vapor interface, and the liquid and vapor densities per time step:
step 5-1: assuming that the phase change of the liquid working medium and the vapor working medium occurs at the interface between the liquid suction core and the vapor cavity, calculating the temperature T at the interface according to the energy conservation equation at the interface i
The conservation of energy equation at the interface is expressed as:
Figure BDA0003421462880000021
wherein:
λ wick and lambda is v Representing the thermal conductivity, W.m, of the wick and vapor chamber, respectively -1 ·K -1
A i Interface area, m 2
T-temperature, where T i Interface temperature, K;
y-distance in the perpendicular direction to the interface, m;
m i interface Mass transfer, m i < 0 is evaporation, m i Condensation > 0 kg;
h fg latent heat of vaporization, J/kg;
C l -liquid phase constant pressure specific heat, J/(kg·k);
C v -gas phase constant pressure specific heat, J/(kg·k);
assuming two control bodies at two sides of a gas-liquid interface are PW and PV respectively, wherein PW represents one side of a liquid suction core, PV represents one side of a vapor cavity, and dispersing an energy conservation equation to obtain the temperature T at the interface i
Figure BDA0003421462880000031
Wherein:
T PW -the first layer mesh temperature, K, of the wick near the interface;
T PV -the first layer mesh temperature, K, of the vapor chamber near the interface;
Δx-interfacial area, m 2
Calculating the temperature T of the interface according to the formula (2) i
Figure BDA0003421462880000032
Step 5-2: calculating the interface pressure P by using Clausius-Clapeyron equation according to the temperature of the interface obtained in the step 5-1 i
Figure BDA0003421462880000041
Wherein:
P o -reference pressure, pa;
P i interface pressure, pa;
T o -reference temperature, K;
r-gas constant, J/kg;
step 5-3: for gas-liquid interface mass flow rate
Figure BDA0003421462880000042
And (3) performing calculation:
Figure BDA0003421462880000043
wherein:
Figure BDA0003421462880000044
-interface mass flow rate, kg/s;
sigma-adjustment factor; surface tension, N/m;
P v -vapour chamber pressure, pa;
T v -vapor chamber temperature, K;
step 5-4: calculating the working pressure of the vapor region according to an ideal gas formula and the overall mass balance of the vapor region, and calculating the vapor density:
Figure BDA0003421462880000045
wherein:
P op -vapour zone working pressure, pa;
Figure BDA0003421462880000046
-the mass of vapor in the vapor chamber calculated for the last time step, kg;
Figure BDA0003421462880000047
-interface mass flow rate, kg/s;
Δt-time step, s;
r-gas constant;
V cell -vapour chamber pressure, pa;
based on the incompressible assumption of the fluid, the vapor density in the control volume is calculated from the system pressure:
Figure BDA0003421462880000051
step 5-5: calculating the average density ρ of the liquid l
The wick is unsaturated during transients, and the average liquid density is calculated to maintain liquid mass balance, regardless of the problem of unsaturation;
calculating the mass of liquid in the liquid suction core at the current moment:
Figure BDA0003421462880000052
wherein:
M l -mass of liquid in the wick, kg;
Figure BDA0003421462880000053
-the mass of liquid in the wick, kg, calculated over a previous time step;
Figure BDA0003421462880000054
-interface mass flow rate, kg/s;
Δt-time step, s;
from the mass of liquid in the wick, the average density of the liquid is calculated:
Figure BDA0003421462880000055
wherein:
epsilon-wick porosity;
V l wick volume/m 3
Step 6: calculating the flow heat transfer in the three-dimensional flat plate hydrothermal pipe by using computational fluid dynamics software:
step 6-1: respectively solving momentum equations of the liquid suction core and the steam cavity according to the interface mass flow rate, the pressure and the temperature parameters obtained in the step 5 to obtain speed distribution, and correcting the pressure by using a pressure correction equation;
step 6-2: solving an energy conservation equation of the calculation region to obtain a temperature field of the calculation region;
step 6-3: checking whether the calculation converges: and when the residual value is smaller than 0.000001, the calculation is considered to reach convergence, if so, the calculation of the next time step is carried out, otherwise, the step 5 is returned to continue the calculation.
Preferably, in step 2, the calculation domain is gridded by using a pure hexahedral grid in consideration of high requirements on grid quality for fine simulation of the heat pipe.
Compared with the prior art, the invention has the following beneficial effects:
1) According to the invention, the wick and the steam cavity are calculated independently, and the heat pipe is calculated by adding the interaction method of the wick and the steam cavity, so that the internal flow of the wick can be simulated more finely;
2) The numerical format of incompressible limit pressure of the system is calculated by using the whole mass balance, and the stability of a standard sequence program is improved by considering the coupling between the mass flow of an evaporation section/a condensation section, the concentration of an interface substance, the pressure and the system pressure;
3) The model is independent, the method is high in universality, and can be used for calculating heat pipes in other shapes.
Drawings
FIG. 1 is a schematic diagram of a three-dimensional flat heat pipe;
FIG. 2 is a schematic diagram of an energy transfer process within a heat pipe;
fig. 3 is a schematic view of a control body at the interface of a wick and a vapor chamber;
fig. 4 is a flow chart of the method of the present invention.
Detailed Description
The present invention will be described in further detail below with reference to the flow chart of fig. 4, taking a typical three-dimensional flat plate water heat pipe as an example:
the invention discloses a numerical calculation method of a three-dimensional flat plate water heat pipe, which comprises the following steps:
step 1: FIG. 1 is a schematic structural diagram of a three-dimensional flat heat pipe, wherein a software solidworks is established through a geometric model to establish a three-dimensional geometric model of the flat water heat pipe according to a real design, and the three-dimensional geometric model comprises a pipe wall, a liquid suction core, a vapor chamber, an evaporation section and a condensation section;
step 2: grid division is carried out on the flat-plate water heat pipe geometric model obtained in the step 1 by utilizing grid division software, and the embodiment uses a pure hexahedral grid to carry out grid division on a calculation domain in consideration of high requirements on grid quality of fine simulation on the heat pipe;
step 3: the grids divided in the step 2 are led into computational fluid dynamics software, an evaporation section and a condensation section of the heat pipe are respectively set as a heat flow boundary condition and a convection boundary condition, other wall surfaces are set as heat insulation wall surfaces, the heat flow density heated by the evaporation section, the heat exchange coefficient of the condensation section and the temperature of a coolant are set according to the actual situation, and the porosity of the liquid absorption core is set;
step 4: setting initial conditions required for transient calculation: setting the initial steam mass as 0, initializing a temperature field, and enabling the working pressure of an initial steam area to be the saturation pressure at the current temperature;
step 5: the user-defined program UDF was used to calculate the gas-liquid interface mass flow rate, the temperature and pressure at the wick-vapor interface, and the liquid and vapor densities per time step:
step 5-1: assuming that the phase change of the liquid working medium and the vapor working medium occurs at the interface between the liquid suction core and the vapor cavity, calculating the temperature T at the interface according to the energy conservation equation at the interface i
The energy transfer process in the heat pipe is shown in fig. 2, in the evaporation section, the liquid is changed into steam under the heating effect, and the heat is transferred from the evaporation section to the steam cavity; in the condensing section, the steam condenses under the cooling action, and heat is transferred from the steam cavity to the condensing section. The conservation of energy equation at the interface is expressed as:
Figure BDA0003421462880000071
wherein:
λ wick and lambda is v Representing the thermal conductivity, W.m, of the wick and vapor chamber, respectively -1 ·K -1
A i Interface area, m 2
T-temperature, where T i Interface temperature, K;
y-distance in the perpendicular direction to the interface, m;
m i interface Mass transfer, m i < 0 is evaporation, m i Condensation > 0 kg;
h fg latent heat of vaporization, J/kg;
C l -liquid phase constant pressure specific heat, J/(kg·k);
C v -gas phase constant pressure specific heat, J/(kg·k);
FIG. 3 is a schematic diagram of the control body at the interface of the wick and the vapor chamber, assuming two control bodies on both sides of the gas-liquid interface are PW and PV, respectively, where PW represents one side of the wick and PV represents one side of the vapor chamber, and the energy conservation equation is discretized to obtain the temperature T at the interface i
Figure BDA0003421462880000081
Wherein:
T PW -the first layer mesh temperature, K, of the wick near the interface;
T PV -the first layer mesh temperature, K, of the vapor chamber near the interface;
Δx-interfacial area, m 2
Calculating the temperature T of the interface according to the formula (2) i
Figure BDA0003421462880000082
Step 5-2: calculating the interface pressure P by using Clausius-Clapeyron equation according to the temperature of the interface obtained in the step 5-1 i
Figure BDA0003421462880000083
Wherein:
P o -reference pressure, pa;
P i interface pressure, pa;
T o -reference temperature, K;
T i interface temperature, K;
h fg latent heat of vaporization, J/kg;
r-gas constant, J/kg;
step 5-3: for gas-liquid interface mass flow rate
Figure BDA0003421462880000091
And (3) performing calculation:
Figure BDA0003421462880000092
wherein:
Figure BDA0003421462880000093
-interface mass flow rate, kg/s;
sigma-adjustment factor; surface tension, N/m;
P v -vapour chamber pressure, pa;
T v -vapor chamber temperature, K;
step 5-4: calculating the working pressure of the vapor region according to an ideal gas formula and the overall mass balance of the vapor region, and calculating the vapor density:
Figure BDA0003421462880000094
wherein:
P op -vapour zone working pressure, pa;
Figure BDA0003421462880000095
-the mass of vapor in the vapor chamber calculated for the last time step, kg;
Figure BDA0003421462880000096
-interface mass flow rate, kg/s;
Δt-time step, s;
r-gas constant;
V cell -vapour chamber pressure, pa;
based on the incompressible assumption of the fluid, the vapor density in the control volume is calculated from the system pressure:
Figure BDA0003421462880000097
step 5-5: calculating the average density ρ of the liquid l
The wick is typically unsaturated during transients, and the average liquid density is calculated to maintain liquid mass balance, regardless of the problem of unsaturation.
Calculating the mass of liquid in the liquid suction core at the current moment:
Figure BDA0003421462880000101
wherein:
M l -mass of liquid in the wick, kg;
Figure BDA0003421462880000102
-in the wick calculated in the last time stepLiquid mass, kg;
Figure BDA0003421462880000103
-interface mass flow rate, kg/s;
Δt-time step, s;
from the mass of liquid in the wick, the average density of the liquid is calculated:
Figure BDA0003421462880000104
wherein:
epsilon-wick porosity;
V l wick volume/m 3
Step 6: calculating the flow heat transfer in the three-dimensional flat plate hydrothermal pipe by using computational fluid dynamics software:
step 6-1: respectively solving momentum equations of the liquid suction core and the steam cavity according to the parameters of the interface mass flow rate, the pressure, the temperature and the like obtained in the step 5 to obtain speed distribution, and correcting the pressure by using a pressure correction equation;
step 6-2: solving an energy conservation equation of the calculation region to obtain a temperature field of the calculation region;
step 6-3: checking whether the calculation converges: and when the residual value is smaller than 0.000001, the calculation is considered to reach convergence, if so, the calculation of the next time step is carried out, otherwise, the step 5 is returned to continue the calculation.

Claims (2)

1. A numerical calculation method of a three-dimensional flat plate water heat pipe is characterized by comprising the following steps of: the method comprises the following steps:
step 1: establishing a software solidworks through a geometric model to establish a three-dimensional geometric model of the flat water heat pipe according to a real design, wherein the three-dimensional geometric model comprises a pipe wall, a liquid suction core, a vapor chamber, an evaporation section and a condensation section;
step 2: performing grid division on the flat plate water heat pipe geometric model obtained in the step 1 by using grid division software;
step 3: the grids divided in the step 2 are led into computational fluid dynamics software, an evaporation section and a condensation section of the heat pipe are respectively set as a heat flow boundary condition and a convection boundary condition, other wall surfaces are set as heat insulation wall surfaces, the heat flow density heated by the evaporation section, the heat exchange coefficient of the condensation section and the temperature of a coolant are set according to the actual situation, and the porosity of the liquid absorption core is set;
step 4: setting initial conditions required for transient calculation: setting the initial steam mass as 0, initializing a temperature field, and enabling the working pressure of an initial steam area to be the saturation pressure at the current temperature;
step 5: the user-defined program UDF was used to calculate the gas-liquid interface mass flow rate, the temperature and pressure at the wick-vapor interface, and the liquid and vapor densities for each time step:
step 5-1: assuming that the phase change of the liquid working medium and the vapor working medium occurs at the interface between the liquid suction core and the vapor cavity, calculating the temperature T at the interface according to the energy conservation equation at the interface i
The conservation of energy equation at the interface is expressed as:
Figure FDA0004097585420000011
wherein:
λ wick and lambda is v Representing the thermal conductivity, W.m, of the wick and vapor chamber, respectively -1 ·K -1
A i Interface area, m 2
T-temperature, where T i Interface temperature, K;
y-distance in the perpendicular direction to the interface, m;
m i interface Mass transfer, m i < 0 is evaporation, m i Condensation > 0 kg;
h fg latent heat of vaporization, J/kg;
C l liquid phase constant pressure specific heat, J/(kg)K);
C v -gas phase constant pressure specific heat, J/(kg·k);
assuming two control bodies at two sides of a gas-liquid interface are PW and PV respectively, wherein PW represents one side of a liquid suction core, PV represents one side of a vapor cavity, and dispersing an energy conservation equation to obtain the temperature T at the interface i
Figure FDA0004097585420000021
Wherein:
T PW -the first layer mesh temperature, K, of the wick near the interface;
T PV -the first layer mesh temperature, K, of the vapor chamber near the interface;
Δx-interfacial area, m 2
Calculating the temperature T of the interface according to the formula (2) i
Figure FDA0004097585420000022
Figure FDA0004097585420000023
Step 5-2: calculating the interface pressure P by using Clausius-Clapeyron equation according to the temperature of the interface obtained in the step 5-1 i
Figure FDA0004097585420000024
Wherein:
P o -reference pressure, pa;
P i interface pressure, pa;
T o -reference temperature, K;
r-gas constant, J/kg;
step 5-3: for gas-liquid interface mass flow rate
Figure FDA0004097585420000031
And (3) performing calculation:
Figure FDA0004097585420000032
wherein:
Figure FDA0004097585420000033
-interface mass flow rate, kg/s;
sigma-adjustment factor; surface tension, N/m;
P v -vapour chamber pressure, pa;
T v -vapor chamber temperature, K;
step 5-4: calculating the working pressure of the vapor region according to an ideal gas formula and the overall mass balance of the vapor region, and calculating the vapor density:
Figure FDA0004097585420000034
wherein:
P op -vapour zone working pressure, pa;
Figure FDA0004097585420000035
-the mass of vapor in the vapor chamber calculated for the last time step, kg;
Figure FDA0004097585420000036
-interface mass flow rate, kg/s;
Δt-time step, s;
r-gas constant;
V cell -vapour chamber pressure, pa;
based on the incompressible assumption of the fluid, the vapor density in the control volume is calculated from the system pressure:
Figure FDA0004097585420000037
step 5-5: calculating the average density ρ of the liquid l
The wick is unsaturated during transients, and the average liquid density is calculated to maintain liquid mass balance, regardless of the problem of unsaturation;
calculating the mass of liquid in the liquid suction core at the current moment:
Figure FDA0004097585420000041
wherein:
M l -mass of liquid in the wick, kg;
Figure FDA0004097585420000042
-the mass of liquid in the wick, kg, calculated over a previous time step;
Figure FDA0004097585420000043
-interface mass flow rate, kg/s;
Δt-time step, s;
from the mass of liquid in the wick, the average density of the liquid is calculated:
Figure FDA0004097585420000044
wherein:
epsilon-wick porosity;
V l wick volume/m 3
Step 6: calculating the flow heat transfer in the three-dimensional flat plate hydrothermal pipe by using computational fluid dynamics software:
step 6-1: respectively solving momentum equations of the liquid suction core and the steam cavity according to the interface mass flow rate, the pressure and the temperature parameters obtained in the step 5 to obtain speed distribution, and correcting the pressure by using a pressure correction equation;
step 6-2: solving an energy conservation equation of the calculation region to obtain a temperature field of the calculation region;
step 6-3: checking whether the calculation converges: and when the residual value is smaller than 0.000001, the calculation is considered to reach convergence, if so, the calculation of the next time step is carried out, otherwise, the step 5 is returned to continue the calculation.
2. The method for calculating the numerical value of the three-dimensional flat plate water heat pipe according to claim 1, which is characterized in that: in step 2, the calculation domain is grid-divided by using a pure hexahedral grid in consideration of high requirements on grid quality in fine simulation of the heat pipe.
CN202111564511.2A 2021-12-20 2021-12-20 Numerical calculation method for three-dimensional flat plate water heat pipe Active CN114239435B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111564511.2A CN114239435B (en) 2021-12-20 2021-12-20 Numerical calculation method for three-dimensional flat plate water heat pipe

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111564511.2A CN114239435B (en) 2021-12-20 2021-12-20 Numerical calculation method for three-dimensional flat plate water heat pipe

Publications (2)

Publication Number Publication Date
CN114239435A CN114239435A (en) 2022-03-25
CN114239435B true CN114239435B (en) 2023-05-16

Family

ID=80759505

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111564511.2A Active CN114239435B (en) 2021-12-20 2021-12-20 Numerical calculation method for three-dimensional flat plate water heat pipe

Country Status (1)

Country Link
CN (1) CN114239435B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107423459A (en) * 2017-03-21 2017-12-01 哈尔滨工程大学 A kind of heat exchanger porous media model porosity and Permeability Parameters processing method based on CAD software
CN107704674A (en) * 2017-09-26 2018-02-16 吉林省电力科学研究院有限公司 The method for numerical simulation of air cooling tubes condenser water vapor condensation process

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190276906A1 (en) * 2011-03-30 2019-09-12 Macrae Technologies, Inc. High heat flux regime coolers
NO337063B1 (en) * 2014-01-24 2016-01-11 Ledaflow Tech Da Method for Transient Quasi-Three-Dimensional Simulation of Multiphase Fluid Flow in Pipelines
CN108694274A (en) * 2018-04-18 2018-10-23 青岛海聚仿真软件技术有限公司 Condense the Numeral Emulation System and its method of phase transformation secondary development analysis
CN109670216B (en) * 2018-11-30 2023-07-04 中国船舶重工集团公司第七一九研究所 Passive residual heat removal condenser position optimization design method based on CFD technology
CN112597640B (en) * 2020-12-08 2024-08-09 清华大学 Heat pipe simulation method and device and electronic equipment
CN113434988B (en) * 2021-06-11 2022-10-28 西安交通大学 Heat transfer characteristic analysis method for heat pipe in nuclear power system
CN113609795B (en) * 2021-08-02 2022-12-09 西安交通大学 Three-dimensional calculation method for thermal hydraulic characteristics of high-temperature gas cooled reactor spiral tube steam generator

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107423459A (en) * 2017-03-21 2017-12-01 哈尔滨工程大学 A kind of heat exchanger porous media model porosity and Permeability Parameters processing method based on CAD software
CN107704674A (en) * 2017-09-26 2018-02-16 吉林省电力科学研究院有限公司 The method for numerical simulation of air cooling tubes condenser water vapor condensation process

Also Published As

Publication number Publication date
CN114239435A (en) 2022-03-25

Similar Documents

Publication Publication Date Title
CN111725157B (en) Intelligent phase change cooling structure based on shape memory material
CN111947219B (en) Compact type step heat storage and supply system and method based on natural circulation heat exchange
Kandilli et al. Optimisation design and operation parameters of a photovoltaic thermal system integrated with natural zeolite
JP7540089B2 (en) A method for calculating the amount of molten salt storage that supports peak adjustment of thermal power generation units using molten salt thermal storage
CN102997729A (en) Heat pipe radiator of phase change drive loop
CN108808160A (en) High intensity heat transfer structure for cooling down power battery
CN102419123B (en) Anti-gravity bubbling loop heat pipe
CN113959244B (en) Double-evaporator condenser loop heat pipe
Ling et al. Experimental investigation of loop heat pipe with novel interlaced microchannel condenser
CN114239435B (en) Numerical calculation method for three-dimensional flat plate water heat pipe
Jose et al. Recent advances in loop heat pipes: A review
CN113013120A (en) Heat dissipation device and electronic equipment
CN115395132A (en) Active control battery liquid cooling device and working method thereof
CN115854764A (en) Integrated finned tube phase-change heat storage heat exchanger
Al-rifai et al. Numerical Study on the Dominant Condensation Mechanism in the Cross-Flow Transport Membrane Condenser
Alkhafaj et al. Optimization of wickless heat pipe heat exchanger using R1234-yf and ethanol as working fluids
Momeni et al. Thermal Response Analysis of a Cross-Flow PCM Heat Exchanger Based on Air and Liquid Flow
CN105486133A (en) Heat pipe flue gas waste heat recycling device and working medium
CN221409554U (en) Circulation evaporation radiator
CN217358194U (en) Steam solid heat storage system
Wei et al. Numerical study on electrical and thermal performance of a novel concentrated Photovoltaic/Thermal system
CN103763888B (en) A kind of heat-pipe radiator and heat pipe radiating system
CN204678743U (en) Evaporimeter and Split hot heat exchange of heat pipe
Li et al. A Parametric Study on the Thermal Homogeneity of Vapor Chamber in a Battery Cooling System
Patel Experimental Investigation Of Solar Air-Water Heater With Heat Exchanger For Industrial Process Heating

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