CN105857644A - Optimized design method of heat pipe radiator - Google Patents

Optimized design method of heat pipe radiator Download PDF

Info

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
Application number
CN201610188808.6A
Other languages
Chinese (zh)
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.)
Shanghai Institute of Satellite Engineering
Original Assignee
Shanghai Institute of Satellite Engineering
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 Shanghai Institute of Satellite Engineering filed Critical Shanghai Institute of Satellite Engineering
Priority to CN201610188808.6A priority Critical patent/CN105857644A/en
Publication of CN105857644A publication Critical patent/CN105857644A/en
Pending legal-status Critical Current

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/46Arrangements or adaptations of devices for control of environment or living conditions
    • B64G1/50Arrangements or adaptations of devices for control of environment or living conditions for temperature control
    • B64G1/506Heat pipes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • 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

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

A kind of Optimization Design of heat pipe radiator
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:
λ ∂ 2 T ∂ x 2 + λ ∂ 2 T ∂ y 2 + Φ · = 0
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:
Φ · = ( Φ i n + Φ o u t ) / V
Φ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:
λ ∂ 2 T ∂ x 2 + λ ∂ 2 T ∂ y 2 + Φ · = 0 - - - ( 1 )
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:
Φ · = ( Φ i n + Φ o u t ) / V - - - ( 2 )
Φ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:
λ ∂ 2 T ∂ x 2 + λ ∂ 2 T ∂ y 2 + Φ · = 0
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:
Φ · = ( Φ i n + Φ o u t ) / V
Φ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.
CN201610188808.6A 2016-03-29 2016-03-29 Optimized design method of heat pipe radiator Pending CN105857644A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
程雪涛等: "空间辐射器的等温化设计", 《工程热物理学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
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