CN105857644A - Optimized design method of heat pipe radiator - Google Patents
Optimized design method of heat pipe radiator Download PDFInfo
- Publication number
- CN105857644A CN105857644A CN201610188808.6A CN201610188808A CN105857644A CN 105857644 A CN105857644 A CN 105857644A CN 201610188808 A CN201610188808 A CN 201610188808A CN 105857644 A CN105857644 A CN 105857644A
- Authority
- CN
- China
- Prior art keywords
- radiator
- heat
- heat pipe
- reduced
- conduction model
- 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.)
- Pending
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/46—Arrangements or adaptations of devices for control of environment or living conditions
- B64G1/50—Arrangements or adaptations of devices for control of environment or living conditions for temperature control
- B64G1/506—Heat pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Biodiversity & Conservation Biology (AREA)
- General Health & Medical Sciences (AREA)
- Toxicology (AREA)
- Environmental & Geological Engineering (AREA)
- Health & Medical Sciences (AREA)
- Aviation & Aerospace Engineering (AREA)
- Environmental Sciences (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)
Abstract
The invention provides an optimized design method of a heat pipe radiator. The method includes the following steps that S1, a two-dimensional steady-state heat conduction model of the radiator is established after simplification; S2, the temperature field distribution of the radiator is obtained through the two-dimensional steady-state heat conduction model; and S3, the heat conduction coefficient of the corresponding position is changed by pre-burying heat pipes in the maximum temperature field temperature gradient position. The method further includes the step of repeating the step S2 and the step S3 until the isothermality design objective requirement of the radiator is met. According to the method, the local heat conduction coefficient is changed by pre-burying the heat pipes in the maximum temperature gradient position of the radiator, so that the heat conduction capacity is improved, and isothermality design of the whole radiator is achieved. According to the method, human intervention is reduced, and the method can be used for optimized design of radiators in any geometrical shape and with any heat source distribution mode.
Description
Technical field
The present invention relates to spacecraft thermal control, in particular it relates to the Optimization Design of a kind of heat pipe radiator.
Background technology
Along with the development of space tasks demand, star occurring the unit of big heat consumption, high-temperature stability, such as certain swashs
Optical radar load, its heat consumption is up to more than 470W, and the stability requirement of temperature is within ± 2 DEG C.It is thus desirable to one
Individual single opaco radiator is as radiating surface.In order to improve the radiating efficiency of radiator, it is necessary to improve radiator
Temperature stability.
Heat pipe radiator is a kind of better method of radiator isothermalization design.It is that heat pipe is embedded in cellular board
In, therefore radiator has cellular board high intensity, high rigidity and the feature of light weight and the high thermal conductivity of heat pipe concurrently.
The optimization design of heat pipe radiator is i.e. in the case of heat pipe radiator area and weight are certain, pre-buried by optimizing
The permutation and combination of heat pipe so that radiator temperature is more uniformly distributed, the heat loss through radiation amount of unit are is bigger.In reality
In the middle of engineer applied, the most by virtue of experience carry out pre-buried heat pipe, but for the radiator of random geometry, often
Difficulty is bigger.
Summary of the invention
For defect of the prior art, it is an object of the invention to provide the optimization design side of a kind of heat pipe radiator
Method, by high temperature or the high heat flux of radiator local, by pre-buried heat pipe, is transferred to whole radiator efficiently,
To ensure that the maximum temperature difference on radiator is not more than a certain permissible value, it is achieved the isothermalization design of radiator.
According to the Optimization Design of the heat pipe radiator that the present invention provides, comprise the steps:
Step S1: set up the two-dimensional steady-state conduction model of radiator;
Step S2: obtained the thermo parameters method of radiator by described two-dimensional steady-state conduction model;
Step S3: the thermograde maximum in described temperature field changes the thermal conductivity factor of relevant position by pre-buried heat pipe.
Preferably, also comprise the steps:
-repeat step S2 to step S3, until meet radiator thermograde less than or equal to isothermal design objective
Requirement.
Preferably, described two-dimensional steady-state conduction model particularly as follows:
Wherein, λ is the thermal conductivity factor of material;T is temperature;X, y are the plane coordinates of radiator;For unit
The endogenous pyrogen of volume, is inputted hot-fluid Φ by loadinWith heat loss through radiation Orbital heat flux ΦoutTwo parts form, load input heat
Stream ΦinRelevant to heat input;Endogenous pyrogenWith heat loss through radiation Orbital heat flux ΦoutExpression formula as follows:
Φout=-ε σ AT4
Wherein, V represents thermal source volume;A is radiator area;ε is emissivity;σ is that this fence-Boltzmann of making a mistake is normal
Number.
Preferably, described two-dimensional steady-state conduction model is set up by hypothesis below:
Radiator sunny slope is reduced to insulating surface;
Four borders of radiator are reduced to adiabatic boundary;
Radiator less than the threshold value set, is then reduced to the two dimension of one side heat radiation by the sunny slope of radiator and opaco thermal resistance
Plane;
The load input heat consumption of radiator surface and the external heat loss through radiation of radiator are all reduced to the endogenous pyrogen of radiator;
By with upper type, the two-dimensional steady-state conduction model being reduced to endogenous pyrogen, adiabatic boundary of radiator.
Compared with prior art, the present invention has a following beneficial effect:
1, the present invention is at radiator thermograde maximum, changes localized thermal conductivity by pre-buried heat pipe, with
Strengthen capacity of heat transmission, and then realize the sammingization design of whole radiator.
2, the present invention can decrease human intervention, can be used for random geometry, the radiation of any thermal source distribution
The optimization design of device.
Accompanying drawing explanation
The detailed description made non-limiting example with reference to the following drawings by reading, other of the present invention is special
Levy, purpose and advantage will become more apparent upon:
Fig. 1 is radiator and the schematic diagram in high heat consumption region in the present invention;
Fig. 2 is the distribution schematic diagram in the temperature field optimizing previous irradiation device in the present invention;
Fig. 3 is the distribution schematic diagram of heat pipe in the present invention;
Fig. 4 be the present invention optimizes after the distribution schematic diagram in temperature field of radiator;
Fig. 5 is the flow chart of steps of the Optimization Design of heat pipe radiator in the present invention.
Detailed description of the invention
Below in conjunction with specific embodiment, the present invention is described in detail.Following example will assist in those skilled in the art
Member is further appreciated by the present invention, but limits the present invention the most in any form.It should be pointed out that, the common skill to this area
For art personnel, without departing from the inventive concept of the premise, it is also possible to make some deformation and improvement.These broadly fall into
Protection scope of the present invention.
In the present embodiment, the Optimization Design of the heat pipe radiator that the present invention provides, comprise the steps:
Step S1: set up the two-dimensional steady-state conduction model of radiator;
Step S2: obtained the thermo parameters method of radiator by described two-dimensional steady-state conduction model;
Step S3: the thermograde maximum in described temperature field changes the thermal conductivity factor of relevant position by pre-buried heat pipe.
The Optimization Design of the heat pipe radiator that the present invention provides, also comprises the steps:
-repetition step S2 is to step S3, until the thermograde meeting radiator is wanted less than or equal to isothermal design objective
Ask.Such as requiring to be 10 centimetres when isothermal design objective, 20 DEG C, then repetition step S2 is to step S3, until 10
In the inside, thermograde is less than or equal to 20 degrees Celsius.
Described two-dimensional steady-state conduction model particularly as follows:
Wherein, λ is the thermal conductivity factor of material;T is temperature;X, y are the plane coordinates of radiator;For unit
The endogenous pyrogen of volume, is inputted hot-fluid Φ by loadinWith heat loss through radiation Orbital heat flux ΦoutTwo parts form, load input heat
Stream ΦinRelevant to heat input;Endogenous pyrogenWith heat loss through radiation Orbital heat flux ΦoutExpression formula as follows:
Φout=-ε σ AT4 (3)
Wherein, V represents thermal source volume;A is radiator area;ε is emissivity;σ is that this fence-Boltzmann of making a mistake is normal
Number.
The thermal control of the high-power load on star is typically with outer patch heat pipe, LHP or pump and drives the method for fluid circuit to carry
The heat consumption of lotus is transferred on radiator.Radiator is generally made up of lead-covering and aluminium honeycomb, its sunny slope cladding multilayer, the back of the body
Sunny side spray-coated white paint is dispelled the heat to the cold space radiation of 4K.
Therefore, described two-dimensional steady-state conduction model is set up by hypothesis below:
Radiator sunny slope cladding multilayer, can be reduced to insulating surface;
Four border heat exchanges of radiator are the least, can be reduced to adiabatic boundary;
Sunny slope and the opaco thermal resistance of radiator are less, and therefore radiator can be reduced to the two dimensional surface of one side heat radiation.
Local load's input heat consumption of radiator surface and the external heat loss through radiation of radiator all can be reduced to the interior of radiator
Thermal source.By with upper type, the two-dimensional steady-state conduction model being reduced to endogenous pyrogen, adiabatic boundary of radiator.
More specifically, it is assumed that the thermal source heat consumption at regional area is up to 470W, according to heat radiation equation,
Φ=A ε σ T4 (4)
It is 2.5m that preresearch estimates obtains the area requirements of radiator2(long 2m, wide 1.25 meters);
Φ is radiant heat, and A is radiator area, and ε is emissivity, and σ is this fence-Boltzmann constant of making a mistake, and T is temperature.
Numerical solution by solving equation (1), it is thus achieved that the initially temperature profile of radiator during the most pre-buried heat pipe, as
Shown in Fig. 2, it can be seen that the thermograde of radiator reaches more than 125 degree.According to the present invention, when the maximum on radiator
When thermograde requires less than 20 degree, the distribution map of pre-buried heat pipe as it is shown on figure 3, in figure red line represent pre-buried heat pipe.
The temperature profile of radiator after optimization, as shown in fig. 4, it can be seen that maximum temperature gradient is less than 20 degree, meets heat
The optimization design requirement of pipe radiator.In the middle of practical implementation, the pre-buried heat pipe of Fig. 3 can be modified, make
Reach equal heat-transfer capability.
Above the specific embodiment of the present invention is described.It is to be appreciated that the invention is not limited in
Stating particular implementation, those skilled in the art can make various deformation or amendment within the scope of the claims,
This has no effect on the flesh and blood of the present invention.
Claims (4)
1. the Optimization Design of a heat pipe radiator, it is characterised in that comprise the steps:
Step S1: simplify and set up the two-dimensional steady-state conduction model of radiator;
Step S2: obtained the thermo parameters method of radiator by described two-dimensional steady-state conduction model;
Step S3: the thermograde maximum in described temperature field changes the thermal conductivity factor of relevant position by pre-buried heat pipe.
The Optimization Design of heat pipe radiator the most according to claim 1, it is characterised in that also include as follows
Step:
-repetition step S2 is to step S3, until the thermograde meeting radiator is wanted less than or equal to isothermal design objective
Ask.
The Optimization Design of heat pipe radiator the most according to claim 1, it is characterised in that described two stability maintenances
State conduction model particularly as follows:
Wherein, λ is the thermal conductivity factor of material;T is temperature;X, y are the plane coordinates of radiator;For unit
The endogenous pyrogen of volume, is inputted hot-fluid Φ by loadinWith heat loss through radiation Orbital heat flux ΦoutTwo parts form, load input heat
Stream ΦinRelevant to heat input;Endogenous pyrogenWith heat loss through radiation Orbital heat flux ΦoutExpression formula as follows:
Φout=-ε σ AT4
Wherein, V represents thermal source volume;A is radiator area;ε is emissivity;σ is that this fence-Boltzmann of making a mistake is normal
Number.
The Optimization Design of heat pipe radiator the most according to claim 1, it is characterised in that described two stability maintenances
State conduction model is set up by hypothesis below:
Radiator sunny slope is reduced to insulating surface;
Four borders of radiator are reduced to adiabatic boundary;
Radiator less than the threshold value set, is then reduced to the two dimension of one side heat radiation by the sunny slope of radiator and opaco thermal resistance
Plane;
The load input heat consumption of radiator surface and the external heat loss through radiation of radiator are all reduced to the endogenous pyrogen of radiator;
By with upper type, the two-dimensional steady-state conduction model being reduced to endogenous pyrogen, adiabatic boundary of radiator.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610188808.6A CN105857644A (en) | 2016-03-29 | 2016-03-29 | Optimized design method of heat pipe radiator |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610188808.6A CN105857644A (en) | 2016-03-29 | 2016-03-29 | Optimized design method of heat pipe radiator |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105857644A true CN105857644A (en) | 2016-08-17 |
Family
ID=56626627
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610188808.6A Pending CN105857644A (en) | 2016-03-29 | 2016-03-29 | Optimized design method of heat pipe radiator |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105857644A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110987503A (en) * | 2019-12-19 | 2020-04-10 | 吉林建筑大学 | Heat pipe cooling system based on porous heat conduction material |
CN112113448A (en) * | 2020-09-23 | 2020-12-22 | 安徽农业大学 | Method for determining equivalent heat transfer coefficient of heat pipe embedded in oscillating plow breast |
US20220073217A1 (en) * | 2020-09-10 | 2022-03-10 | Brigham Young University (Byu) | Dynamic radiative thermal management of spacecraft |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008243063A (en) * | 2007-03-28 | 2008-10-09 | Mazda Motor Corp | Design support system for vehicle |
CN105160076A (en) * | 2015-08-07 | 2015-12-16 | 中国石油天然气股份有限公司 | Calculation method and device for annular air injection heat insulation parameters |
-
2016
- 2016-03-29 CN CN201610188808.6A patent/CN105857644A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008243063A (en) * | 2007-03-28 | 2008-10-09 | Mazda Motor Corp | Design support system for vehicle |
CN105160076A (en) * | 2015-08-07 | 2015-12-16 | 中国石油天然气股份有限公司 | Calculation method and device for annular air injection heat insulation parameters |
Non-Patent Citations (1)
Title |
---|
程雪涛等: "空间辐射器的等温化设计", 《工程热物理学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110987503A (en) * | 2019-12-19 | 2020-04-10 | 吉林建筑大学 | Heat pipe cooling system based on porous heat conduction material |
US20220073217A1 (en) * | 2020-09-10 | 2022-03-10 | Brigham Young University (Byu) | Dynamic radiative thermal management of spacecraft |
CN112113448A (en) * | 2020-09-23 | 2020-12-22 | 安徽农业大学 | Method for determining equivalent heat transfer coefficient of heat pipe embedded in oscillating plow breast |
CN112113448B (en) * | 2020-09-23 | 2021-07-09 | 安徽农业大学 | Method for determining equivalent heat transfer coefficient of heat pipe embedded in oscillating plow breast |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Basak et al. | Entropy generation during natural convection in a porous cavity: effect of thermal boundary conditions | |
Chen | Progress in study on constructal theory and its applications | |
Villar et al. | Numerical 3-D heat flux simulations on flat plate solar collectors | |
Abreu et al. | An experimental study of two-phase closed thermosyphons for compact solar domestic hot-water systems | |
Kabeel et al. | Modified mathematical model for evaluating the performance of water-in-glass evacuated tube solar collector considering tube shading effect | |
Al Siyabi et al. | Experimental and numerical thermal analysis of multi-layered microchannel heat sink for concentrating photovoltaic application | |
CN105857644A (en) | Optimized design method of heat pipe radiator | |
Kaluri et al. | Entropy generation minimization versus thermal mixing due to natural convection in differentially and discretely heated square cavities | |
Nardini et al. | Heat transfer experiment on natural convection in a square cavity with discrete sources | |
Talukdar et al. | Finite volume method in 3-D curvilinear coordinates with multiblocking procedure for radiative transport problems | |
Yasin et al. | Impact of uniform and non-uniform heated rods on free convective flow inside a porous enclosure: finite element analysis | |
Wang et al. | Combined heat transfer by natural convection–conduction and surface radiation in an open cavity under constant heat flux heating | |
CN103434659A (en) | High-temperature-uniformity heat radiator for navigation satellite | |
Aly et al. | Double-diffusive convection of solid particles in a porous X-shaped cavity filled with a nanofluid | |
Dannelley et al. | Natural convection heat transfer from fractal-like fins | |
Kaluri et al. | Bejan's heatlines and numerical visualization of heat flow and thermal mixing in various differentially heated porous square cavities | |
Assari et al. | An energy and exergy analysis of water and air with different passage in a solar collector | |
Singh et al. | Role of entropy generation on thermal management during natural convection in a tilted square cavity with isothermal and non-isothermal hot walls | |
Rchida et al. | Two dimensions modelization of radiative exchange in polar bears’ Fur | |
CN107995825A (en) | The high heat consumption unit lightweight cooling system of satellite based on heat pipe | |
Limia et al. | Thermal modeling and efficiency of a dual-stage sodium heat engine | |
CN105744804A (en) | Temperature control system of high-heat consumption and high-stability single machine | |
Wright | Application of the thermodynamics of radiation to dyson spheres as work extractors and computational engines and their observational consequences | |
Souza et al. | Constructal design of high-conductivity inserts | |
Nesterov et al. | Experimental investigations of flat T-shaped copper and titanium heat pipes |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160817 |