CN105138727B - Modeling method of the underground power station based on the vertical enclosed busbar nature hot pressing of length - Google Patents
Modeling method of the underground power station based on the vertical enclosed busbar nature hot pressing of length Download PDFInfo
- Publication number
- CN105138727B CN105138727B CN201510437398.XA CN201510437398A CN105138727B CN 105138727 B CN105138727 B CN 105138727B CN 201510437398 A CN201510437398 A CN 201510437398A CN 105138727 B CN105138727 B CN 105138727B
- Authority
- CN
- China
- Prior art keywords
- vertical
- equation
- power station
- length
- coefficient
- 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
Links
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Air Conditioning Control Device (AREA)
Abstract
It is that physical and mathematical modeling is carried out by the thermal balance state in the shaft space to placing enclosed busbar, to calculate the process of Temperature Distribution in vertical shaft, specific steps include the invention discloses modeling method of the underground power station based on the vertical enclosed busbar nature hot pressing of length：1) according to mass conservation law, the law of conservation of momentum and law of conservation of energy, governing equation is established to the thermal balance in shaft space；2) to the value of pressure p carry out it is assumed that solve the equation of momentum in air vertical direction average speed components ui；3) u obtained by calculationiAnd continuity equation, the value of pair pressure p assumed are modified；4) average speed components u is corrected by revised pressure pi, repeat step 1)~4) and convergence solution is obtained, to obtain detailed temperature field distribution in vertical shaft.
Description
Technical field
The invention belongs to hydraulic and hydroelectric engineering field of electromechanical technology, and long vertical closing is based in particular to a kind of underground power station
The modeling method of busbar nature hot pressing.
Background technology
Large Underground Power Station grows vertical enclosed busbar Temperature Distribution and heat dissipation problem is the guarantee of bus run safety.With
Domestic hydroelectric project is grown rapidly, and application of the long vertical isolated-phase enclosed bus (IPB) in underground power station is also more and more extensive,
But to the IPB of different capabilities and vertical height, without relevant design, manufacturer's standard or regulation.The cooling of each underground power station
Mode is cooled down by gravity-flow ventilation, such as Pengshui Hydropower Station, water cloth a strip of land between hills power station, enclosed busbar rated current are respectively
14kA, 16kA, the vertical IPB of the above-mentioned water power head of a station be one hole lead-out mode of a machine, and gravity-flow ventilation, forced ventilation is respectively adopted
Mode meets the requirement of IPB thermal balances.
The increase of the vertical specified heat flow of enclosed busbar and the complication of arrangement are grown with underground power station, by tradition
The Temperature Distribution of mode calculates and radiating mode is difficult to meet the requirement of Large Underground Power Station enclosed busbar.With Three Gorges underground power station
For 700MW grades long vertical high current isolated-phase enclosed bus, rated current reaches 26kA, vertical height 80m, using two machines, one hole
Mode, be to have been put into operation in the world with the maximum long vertical section high current of rated current, rated capacity in project under construction from mutually sealing
The arrangement of close nut line, two machines, one hole is also that the world is pioneering.The IPB calorific values of Three Gorges underground power station are huge, and IPB is in rated capacity
Heat dissipation capacity reaches about 1kW/m phases under operating condition, and the calorific value in IPB vertical shafts along every meter of short transverse reaches 6kW.It is so big
Calorific value cannot such as exclude in time, and IPB temperature is caused quickly to increase, and substantially exceed permission temperature, not only caused by IPB be lost
It increases, and the temperature capacity beyond IPB, seriously affects the structurally and electrically insulation performance of busbar, even result in IPB forfeitures
Working performance is unable to operate normally.
However, heavy current and special arrangement make the nature in IPB the and IPB vertical shafts of the Large Underground Power Station
Ventilation has the characteristics that heat-flash source continuous heating in restricted clearance, as busbar persistently causes heat to air heat effect in vertical shaft
Air gathers on vertical shaft top, since the free convection that hot pressing generates is not enough to take away heat, by the mode of free convection
Heat dissipation and the requirement of Temperature Distribution, which whether can be fully met, need to carry out detailed calculating analysis.
The heat-flash stream of vertical enclosed busbar and the complexity of arrangement are grown for large-scale power station, needs the heat according to different IPB
Source strength, the factors detailed analysis such as different shaft height calculate, and find systematic IPB thermal balances rule, are similar from now on
The Temperature Distribution and advanced reliable radiating mode for ensureing IPB safe and reliable operations in engineering provide guidance.
Invention content
In order to solve, current Temperature Distribution in a traditional way calculates and radiating mode is difficult to meet Large Underground Power Station envelope
The problem of requirement of close nut line, the present invention provides a kind of can go out IPB temperature point according to heat source strengths different IPB, detailed analysis
Modeling method of the underground power station of cloth based on the vertical enclosed busbar nature hot pressing of length.
To achieve the above object, modeling of the underground power station based on the vertical enclosed busbar nature hot pressing of length designed by the present invention
Method is to carry out physical and mathematical modeling by the thermal balance state in the shaft space to placing enclosed busbar, to calculate vertical shaft
The process of interior Temperature Distribution, is characterized in that, specific steps include：
1) according to mass conservation law, the law of conservation of momentum and law of conservation of energy, to each point in shaft space
Thermal balance state establishes governing equation：
Continuity equation：
The equation of momentum：
Energy equation：
In formula：ui：Average speed components of the air in vertical direction；uj：Air average speed components in the horizontal direction；
xi：Vertical height：t：Time；ρ：Atmospheric density；p：Pressure；v：The laminar flow coefficient of viscosity；vt：Turbulence factor；gi：Vertical direction weight
Power acceleration；β：Air thermal expansion coefficient；T：At hoistway entrance and exit mean temperature；T∞：Environment temperature；T：Practical temperature
Degree；Γ：Generalized diffusion process coefficient；I：Radiation intensity；
2) to the value of pressure p carry out it is assumed that solve the equation of momentum in air vertical direction average speed components ui；
3) u obtained by calculationiAnd continuity equation, the value of pair pressure p assumed are modified so that the air gone out
In the average speed components u of vertical directioniMeet continuity equation；
4) by revised pressure p and air vertical direction average speed components uiValue substitute into energy equation, ask
Solve actual temperature T；
5) step 1)~4 are repeated) until calculating the actual temperature T of each point in vertical shaft.
Preferably, in the energy equation of the step 1), generalized diffusion process coefficient Γ and Prandtl number PrWith σTThere are following passes
It is formula：
Γ=v/Pr+vt/σT
In formula：v：The laminar flow coefficient of viscosity；vt：Turbulence factor；σT：Empirical takes 0.9~1；Pr：Prandtl number.
Preferably, in the energy equation of the step 1), radiation intensity I is calculated using multi-surface radiation patterns.
Most preferably, the computation model of the radiation intensity I is：
In formula：Position vector；s：Vertical length；Direction vector；a：Absorption coefficient；σs：Scattering coefficient；n：Refraction
Coefficient；σ：Stefan-Boltzmann constant；T：Actual temperature；Φ：Phase function；Ω'：Space multistory angle.
The present invention passes through to thermal balance and temperature under the conditions of the vertical enclosed busbar gravity-flow ventilation of Large Underground Power Station heat-flash stream length
Degree distribution carries out numerical computations and field measurement comparative study, temperature under the conditions of obtaining the gravity-flow ventilation in the restricted clearance of heat-flash source
The regularity of distribution and feature are spent, and analyzes vertical shaft size, influence of the factors such as heat source strength, busbar height to IPB Temperature Distributions, from
And the heat dissipation strategy that parameter and reply different situations that crucial controlling influences are played on its Temperature Distribution can be found out.
The present invention is by growing underground power station the thermal balance research of vertical IPB and putting into practice, find and propose influence and is long vertical
The thermally equilibrated key parameter of straight enclosed busbar：Heat source strength.Due to the high hot-fluid of Large Underground Power Station grow vertical enclosed busbar by
Heat-flash source heat-dissipating in limit space shows the complexity of the heat dissipations such as its convection current, radiation；Meanwhile the heat flow rate per unit area and consumption of busbar
Damage, the geometric parameter of busbar vertical shaft, the arrangement of busbar in the shaft in its restricted clearance convection current and radiation dissipate
Thermosetting influences, to affect IPB Temperature Distributions and the temperature difference in vertical height direction.It is computed analysis, busbar heat source
Intensity is the most important influence factor of thermal balance of busbar.
The present invention is distributed different, the different location busbar in restricted clearance according to the busbar in isolated-phase enclosed bus vertical shaft
Local temperature and temperature distributing rule along vertical height direction complete three-dimensional Numerical-Mode in different Variation Features, will be passed through
Quasi-, under conditions of providing different operating parameters, analysis high current grows vertical isolated-phase enclosed bus temperature field, the single mother of different location
Line vertical temperature distribution rule imports and exports difference variation rule, and then the high current explored and recognized under mixed ventilation pattern is long
Hot-fluid in vertical isolated-phase enclosed bus vertical shaft and Temperature Distribution performance, and its different operating parameters flow it shadow with heat exchange
Ring rule.
Description of the drawings
Fig. 1 is IPB vertical shaft sectional views.
Fig. 2 is IPB vertical shaft plan views.
Specific implementation mode
Below in conjunction with the drawings and specific embodiments, the present invention is described in further detail.
1) thermal balance in IPB vertical shafts
The heat dissipation of IPB conductors is carried out by the heat radiation between conductor and shell, shell inner air convection, and shell is then led
If passing through the heat loss through convection of air outside heat radiation, shell.Heat radiation is determined by radiation coefficient (for constant), in IPB vertical direction
Upper no difference, therefore influence the convection current of the IPB mainly inside and outside air of IPB of Temperature Distribution in vertical direction.
PM=QMF+QMD
PM+PK=QKF+QKD
PMFor the power attenuation of conductor, W/ (m phases)；
QMFIt is IPB to the heat loss through radiation of shell；
QMDConvection current heat transfer between conductor and shell；
PKFor the power attenuation of conductor, shell, W/ (m phases)；
QKFFor the heat loss through radiation of shell；
QKDFor the heat loss through convection of air outside shell and shell.
The fever of IPB itself is since caused by the loss of conductor and shell, conductor and shell directly have convection current, radiation
Its mechanism of conducting heat is complex, it is therefore desirable to which model carries out following simplify：
(1) due to consideration that the hot property of air, viscosity, thus consider conductor by 100% heat transfer to shell, and
The length direction of conductor ignores loss of the conductor with gas itself.
(2) in each section of conductor of IPB vertical direction, shell, hole wall only between this section of object generation heat transfer, each section without heat
Conduction, ignores the thermal conduction resistance of IPB conductors, case material.
2) physical model of IPB vertical shafts
The physical model of IPB vertical shafts is as depicted in figs. 1 and 2, and Fig. 1 is that IPB vertical shaft sectional views Fig. 2 is IPB vertical shaft plan views.
IPB and its vertical shaft can be reduced to heat-generating cylindrical body and be placed in the enclosure space of vertical shaft formation, and forming one has endogenous pyrogen
Restricted clearance heat exchange models.
3) mathematical model of IPB vertical shafts
The heat dissipation problem of underground power station IPB is a complicated heat transfer process, if solved using conventional method, no
Only formula is complicated, and needs to make a large amount of simplified to mathematical model and it is assumed that necessarily cause result of calculation deviation practical.With
The development of computer technology can increase substantially calculating essence for being solved using numerical computation method under complex situations
Degree.
Specific steps include：
1) according to mass conservation law, the law of conservation of momentum and law of conservation of energy, to each point in shaft space
Thermal balance state establishes governing equation：
Continuity equation：
The equation of momentum：
Energy equation：
In formula：ui：Average speed components of the air in vertical direction；uj：Air average speed components in the horizontal direction；
xi：Vertical height：t：Time；ρ：Atmospheric density；p：Pressure；v：The laminar flow coefficient of viscosity；vt：Turbulence factor；gi：Vertical direction weight
Power acceleration；β：Air thermal expansion coefficient；T：At hoistway entrance and exit mean temperature；T∞：Environment temperature；T：Practical temperature
Degree；Γ：Generalized diffusion process coefficient；I：Radiation intensity；
2) to the value of pressure p carry out it is assumed that solve the equation of momentum in air vertical direction average speed components ui；
3) u obtained by calculationiAnd continuity equation, the value of pair pressure p assumed are modified so that the air gone out
In the average speed components u of vertical directioniMeet continuity equation；
4) by revised pressure p and air vertical direction average speed components uiValue substitute into energy equation, ask
Solve actual temperature T；
5) step 1)~4 are repeated) until calculating the actual temperature T of each point in vertical shaft.
Wherein, in the energy equation of step 1), generalized diffusion process coefficient Γ and Prandtl number PrWith σTThere are following relational expressions：
Γ=v/Pr+vt/σT
In formula：v：The laminar flow coefficient of viscosity；vt：Turbulence factor；σT：Empirical takes 0.9~1；Pr：Prandtl number.
Prandtl number Pr can be according to formula Pr=cpV/ λ are found out, in formula：cpFor specific heat at constant pressure, λ is thermal coefficient.
The essence of underground power station IPB heat dissipation problems really has with the non-isothermal turbulent-flow heat-exchanging of endogenous pyrogen, external heat radiation
Under comprehensive function, the temperature field of the air in vertical shaft restricted clearance, velocity field reach the process accordingly balanced.This limited sky
Interior air-flow is by adding thermogenetic buoyancy lift to be the thermal current effect from bottom to top of power by IPB shells.Numerical-Mode
Analog model uses Standard law of wall method using the K-ε two-equation models for considering buoyancy lift, wall surface.It is solved using SIMPLE algorithms
Discrete, discrete equation use QUICK formats to prevent pseudo- diffusion, use the progress of multilist surface radiation (S2S) model
Radiation calculates.
The computation model of radiation intensity I is：
In formula：Position vector；s：Vertical length；Direction vector；a：Absorption coefficient；σs：Scattering coefficient；n：Refraction
Coefficient；σ：Stefan-Boltzmann constant；T：Actual temperature；Φ：Phase function；Ω'：Space multistory angle.Space multistory angle
The parameters such as Ω ', absorption coefficient a have model structure and material character to determine.
The content not being described in detail in this specification belongs to the prior art well known to professional and technical personnel in the field.
Claims (4)
1. modeling method of the underground power station based on the vertical enclosed busbar nature hot pressing of length is by the vertical shaft to placing enclosed busbar
Thermal balance state in space carries out physical and mathematical modeling, to calculate the process of Temperature Distribution in vertical shaft, it is characterised in that specific
Step includes：The physical model of IPB vertical shafts is placed in heat-generating cylindrical body in the enclosure space of vertical shaft formation, forms one
Restricted clearance heat exchange models with endogenous pyrogen；
1) flat to the heat that each is put in shaft space according to mass conservation law, the law of conservation of momentum and law of conservation of energy
Weighing apparatus situation establishes governing equation：
Continuity equation：
The equation of momentum：
Energy equation：
In formula：ui：Average speed components of the air in vertical direction；uj：Air average speed components in the horizontal direction；xi：It hangs down
Straight height：xj：Horizontal length：t：Time；ρ：Atmospheric density；p：Pressure；v：The laminar flow coefficient of viscosity；vt：Turbulence factor；gi：Vertically
Direction acceleration of gravity；β：Air thermal expansion coefficient；At hoistway entrance and exit mean temperature；T∞：Environment temperature；T：
Actual temperature；Γ：Generalized diffusion process coefficient；I：Radiation intensity；
2) to the value of pressure p carry out it is assumed that solve the equation of momentum in air vertical direction average speed components ui；
3) u obtained by calculationiAnd continuity equation, the value of pair pressure p assumed are modified so that the air gone out is hanging down
Histogram to average speed components uiMeet continuity equation；
4) by revised pressure p and air vertical direction average speed components uiValue substitute into energy equation, solve practical
Temperature T；
5) step 1)~4 are repeated) until calculating the actual temperature T of each point in vertical shaft.
2. modeling method of the underground power station according to claim 1 based on the vertical enclosed busbar nature hot pressing of length, feature
It is：In the energy equation of the step 1), generalized diffusion process coefficient Γ and Prandtl number PrWith σTThere are following relational expressions：
Γ=v/Pr+vt/σT
In formula：v：The laminar flow coefficient of viscosity；vt：Turbulence factor；σT：Empirical takes 0.9~1；Pr：Prandtl number.
3. modeling method of the underground power station according to claim 1 or 2 based on the vertical enclosed busbar nature hot pressing of length, special
Sign is：In the energy equation of the step 1), radiation intensity I is calculated using multi-surface radiation patterns.
4. modeling method of the underground power station according to claim 3 based on the vertical enclosed busbar nature hot pressing of length, feature
It is：The computation model of the radiation intensity I is：
In formula：Position vector；s：Vertical length；Direction vector；a：Absorption coefficient；σs：Scattering coefficient；n：Refraction coefficient；σ：
Stefan-Boltzmann constant；T：Actual temperature；Φ：Phase function；Ω'：Space multistory angle.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510437398.XA CN105138727B (en) | 2015-07-23 | 2015-07-23 | Modeling method of the underground power station based on the vertical enclosed busbar nature hot pressing of length |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510437398.XA CN105138727B (en) | 2015-07-23 | 2015-07-23 | Modeling method of the underground power station based on the vertical enclosed busbar nature hot pressing of length |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105138727A CN105138727A (en) | 2015-12-09 |
| CN105138727B true CN105138727B (en) | 2018-10-16 |
Family
ID=54724074
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510437398.XA Active CN105138727B (en) | 2015-07-23 | 2015-07-23 | Modeling method of the underground power station based on the vertical enclosed busbar nature hot pressing of length |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105138727B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105912505B (en) * | 2016-04-28 | 2018-05-11 | 河海大学 | A kind of gaseous state based on power-law function describes method |
| CN105893325A (en) * | 2016-06-03 | 2016-08-24 | 江西理工大学 | Method for judging stability of metal mine artificial pillar |
| CN110222954B (en) * | 2019-05-20 | 2022-12-23 | 中国神华能源股份有限公司 | Temperature analysis method and device for long vertical section phase-isolated enclosed bus of thermal power plant |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN204333850U (en) * | 2014-12-25 | 2015-05-13 | 长江勘测规划设计研究有限责任公司 | Long vertical heavy-current enclosed busbar thermal equilibrium control system and device |
-
2015
- 2015-07-23 CN CN201510437398.XA patent/CN105138727B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN204333850U (en) * | 2014-12-25 | 2015-05-13 | 长江勘测规划设计研究有限责任公司 | Long vertical heavy-current enclosed busbar thermal equilibrium control system and device |
Non-Patent Citations (7)
| Title |
|---|
| 《Analyses of Three-Dimensional Eddy Current Field and Thermal Problems in an Isolated Phase Bus》;S. L. Ho等;《IEEE Transactions on Magnetics》;20030513;第39卷(第3期);第1515-1518页 * |
| 《Seitevare水电站长垂直封闭母线设计考察与研究》;梁波等;《人民长江》;20061120;第37卷(第11期);第38-40页 * |
| 《垂直布置的母线温度场计算及分析》;张兴娟等;《中国电机工程学报》;20000830;第20卷(第8期);第11-13页 * |
| 《大电流长垂直离相封闭母线热平衡模拟及实测分析》;郭建辉等;《暖通空调》;20150215;第45卷(第2期);第48-49页 * |
| 《封闭母线导体热平衡计算》;赵丽娟等;《水力发电》;20070612;第33卷(第6期);第58-59、68页 * |
| 《建筑一体化太阳能烟囱的通风性能研究》;杨启容等;《太阳能学报》;20100728;第31卷(第7期);第874-875页第1.2节 * |
| 《竖直通道内逆浮升力方向的混合对流传热研究》;解金海;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20080315(第3期);第66页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105138727A (en) | 2015-12-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yang et al. | Effect of fin number on the melting phase change in a horizontal finned shell-and-tube thermal energy storage unit | |
| Kasaeian et al. | 3D simulation of solar chimney power plant considering turbine blades | |
| Al-Kayiem et al. | Mathematical analysis of the influence of the chimney height and collector area on the performance of a roof top solar chimney | |
| Kamkari et al. | Experimental investigation of phase change material melting in rectangular enclosures with horizontal partial fins | |
| CN105138727B (en) | Modeling method of the underground power station based on the vertical enclosed busbar nature hot pressing of length | |
| Guo et al. | Numerical analysis of the optimal turbine pressure drop ratio in a solar chimney power plant | |
| Kent | Numerical analysis of laminar natural convection in isosceles triangular enclosures for cold base and hot inclined walls | |
| Gao et al. | Performance prediction of wet cooling tower using artificial neural network under cross-wind conditions | |
| CN107871037A (en) | A kind of 35kV dry reactors internal temperature evaluation method | |
| CN110331946A (en) | Polar region drilling platforms derrick attemperator and its design method | |
| Hu et al. | Natural convection heat transfer for eutectic binary nitrate salt based Al2O3 nanocomposites in solar power systems | |
| CN105006791B (en) | Thermal balance temperature difference control method based on the natural hot pressing of long vertical enclosed busbar | |
| CN106529010A (en) | Method for designing anti-condensation ring main unit shell by using finite element model | |
| CN106202873B (en) | A kind of appraisal procedure of the aluminium cell three-dimensional stove side based on heat-flow coupling | |
| Lu et al. | Numerical simulation of temperature field in rotor of large turbo generator with air-coolant. | |
| CN105116727B (en) | Optimal control method based on long vertical enclosed busbar force ventilation amount | |
| Yu et al. | Design analysis and improvement of cooler in positive-pressure explosion-proof low-speed high-capacity induction motors | |
| Li et al. | Calculation of thermal performance in amorphous core dry-type transformers | |
| Xu et al. | Modeling of natural ventilation in built-in photovoltaic-Trombe wall | |
| CN110222954B (en) | Temperature analysis method and device for long vertical section phase-isolated enclosed bus of thermal power plant | |
| Talabi et al. | Flow field and Reynolds stress distribution in low turbulence natural convection in a triangular cavity | |
| Preetham et al. | Numerical Simulation of Natural Convection In a North-Light Roof Under Winter-Day Boundary Conditions | |
| Kant et al. | Numerical simulation model of solar air heater using ansys-14 | |
| Xin et al. | Experimental and numerical analysis of solar enhanced natural draft dry cooling tower | |
| Lv et al. | A novel layout of the heat-sink base with a high cooling efficiency for polysilicon reduction furnace |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |