CN109992846A - The emulation mode of solar cross-season underground pipe accumulation of heat - Google Patents
The emulation mode of solar cross-season underground pipe accumulation of heat Download PDFInfo
- Publication number
- CN109992846A CN109992846A CN201910193763.5A CN201910193763A CN109992846A CN 109992846 A CN109992846 A CN 109992846A CN 201910193763 A CN201910193763 A CN 201910193763A CN 109992846 A CN109992846 A CN 109992846A
- Authority
- CN
- China
- Prior art keywords
- fluid
- temperature
- heat
- heat exchanger
- emulation mode
- 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.)
- Granted
Links
Classifications
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/10—Geothermal energy
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of emulation modes of solar cross-season ground heat exchanger flowing heat transfer numerical simulation, obtain energy equation first;The relevant dimensionless number Nu of fluid flowing heat transfer is calculated, derives to obtain fluid convection coefficient of heat transfer hf;Given fluid domain initial temperature;Crank-Nicolson format is constructed as fluid temperature (F.T.), with fluid convection coefficient of heat transfer hfConstitute the third boundary condition at ground heat exchanger inner wall;The energy equation being derived by is iteratively solved to obtain temperature field using Three-dimensional unstructured grid Finite Volume Method for Air;After each iterative solution, the heat loss of unit length fluid is obtained according to temperature computation near wall, subsequent time fluid temperature (F.T.) is calculated using first body heat balance method;Aforesaid operations are repeated, until the temperature field of zoning reaches stable state;The above results are post-processed again, obtain the profiling temperatures in the region.This method can efficient, Accurate Prediction ground heat exchanger flowing and heat transfer conditions.
Description
Technical field
The present invention relates to solar energy heat utilization field more particularly to a kind of emulation of solar cross-season underground pipe accumulation of heat
Method.
Background technique
In recent years, China's most area is perplexed by haze deeply, seriously threatens national health and living standard.
By investigation, coal-burning stove for heating is one of the main reason for causing haze.But existing " with gas for coal " and " using electricity instead of coal " cost compared with
Height, and solar energy occupies more and more shares as a kind of cleaning, abundant, the cheap energy in energy consumption mix.Too
Positive energy soil cross-season heat-storage system exactly utilizes solar heat to heat circulatory mediator in system loop, is changed by underground pipe
Hot device exchanges circulatory mediator with soil progress caloic, to achieve the purpose that solar energy " summer Chu Dongyong ".
The complicated because being known as many and relationship of ground heat exchanger working performance is influenced, therefore is difficult in the design process
Accurate Prediction is carried out to the heat transfer process between underground pipe and soil.Experimental study is only carried out under existence conditions or establishes demonstration work
Journey research underground pipe soil coupling heat-accumulating process be it is unpractical, not only the test period is very long, invest it is huge, and be directed to single system
System is difficult to carry out optimizing research.In this case, based on by correlation flowing with heat transfer knowledge, suitable mathematics mould is established
Type is great high efficiency and economy to the expansion research of this kind of problems using numerical computation method.
The core of across season soil thermal storage numerical simulation study be fluid interchange process in underground pipe and underground pipe with
Heat transmission process between surrounding soil.Due to needing to model simultaneously to tube fluid region and the outer solid area of pipe,
The calculation amount of entire numerical solution process is very big, and calculating cycle is very long, and domestic and international existing calculation method is difficult to meet practical work
The demand of journey.
Summary of the invention
Based on the problems of prior art, the object of the present invention is to provide a kind of solar cross-season underground pipe accumulation of heats
Emulation mode, can realize under the premise of efficient computational efficiency between flow process in underground pipe and underground pipe soil
The numerical value Accurate Prediction of heat transfer process.
The purpose of the present invention is what is be achieved through the following technical solutions:
Embodiment of the present invention provides a kind of emulation mode of solar cross-season underground pipe accumulation of heat, for changing underground pipe
The accumulation of heat of hot device carries out simulation calculation, comprising the following steps:
Step 1, the relevant parameter of the ground heat exchanger is obtained;
Step 2, applied energy Conservation Relationship and thermally conductive Fourier law obtain the solid area of the ground heat exchanger
Energy equation;
Step 3 is calculated with flowing heat transfer Empirical Equation Gnielinski formula combination tube fluid physical parameter
The relevant dimensionless number Nu of fluid flowing heat transfer is derived from fluid convection coefficient of heat transfer h according to the dimensionless number Nuf;
Step 4, the initial temperature for giving the ground heat exchanger fluid domainTemperature as current time fluidWith subsequent time fluid temperature (F.T.)Initial value;
Step 5, using fluid current time temperatureWith subsequent time temperatureConstruct Crank-Nicolson lattice
Formula constitutes the third class at the inner wall of the pipe face of the ground heat exchanger with fluid convection coefficient of heat transfer hf as fluid temperature (F.T.)
Boundary condition;
Step 6 changes to energy equation obtained in the step 1 using Three-dimensional unstructured grid Finite Volume Method for Air
In generation, solves, and obtains Temperature Distribution;
Step 7, the step 6 every time iteratively solve after, according to the pipeline near wall of the ground heat exchanger
The temperature computation at place obtains the heat loss of unit length fluid, and subsequent time fluid temperature (F.T.) is calculated using first body heat balance method
The operation of step 8, repeating said steps 4~7, until the temperature of the solid area of the ground heat exchanger of calculating
Degree field reaches stable state;
Step 9, the result obtained again to above-mentioned solution post-process, and obtain the profiling temperatures in the region.
As seen from the above technical solution provided by the invention, solar cross-season provided in an embodiment of the present invention is buried
The emulation mode of pipe accumulation of heat, it has the advantage that:
Since the convection transfer rate at fluid near wall is directly calculated by Empirical Equation, and use
Crank-Nicolson format obtains fluid temperature (F.T.), grid dividing is carried out without the fluid domain to ground heat exchanger, without logical
It crosses solution turbulence model and obtains fluid flowing and heat transfer relevant parameter at underground pipe inner wall, effectively reduce the meter of such problem
Calculation amount, largely improves computational efficiency, under the premise of guaranteeing this method high efficiency, realizes to solar cross-season
The Accurate Prediction of soil thermal storage Numerical Temperature.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, required use in being described below to embodiment
Attached drawing be briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for this
For the those of ordinary skill in field, without creative efforts, it can also be obtained according to these attached drawings other
Attached drawing.
Fig. 1 is that the emulation mode process of solar cross-season underground pipe flowing heat transfer provided in an embodiment of the present invention is illustrated
Figure;
The underground pipe outlet temperature that Fig. 2 emulation mode provided in an embodiment of the present invention obtains and turbulence model calculated result
Comparison schematic diagram;
Fig. 3 obtains for what emulation mode provided in an embodiment of the present invention obtained along journey fluid temperature (F.T.) and turbulence model calculation method
To result comparison schematic diagram.
Specific embodiment
Below with reference to particular content of the invention, technical solution in the embodiment of the present invention is clearly and completely retouched
It states, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Based on the present invention
Embodiment, every other embodiment obtained by those of ordinary skill in the art without making creative efforts,
Belong to protection scope of the present invention.The content being not described in detail in the embodiment of the present invention belongs to professional and technical personnel in the field
The well known prior art.
As shown in Figure 1, the embodiment of the present invention provides a kind of emulation mode of solar cross-season underground pipe accumulation of heat, for pair
The accumulation of heat of ground heat exchanger carries out simulation calculation, comprising the following steps:
Step 1, the relevant parameter of the ground heat exchanger is obtained;
Step 2 obtains the solid area applied energy Conservation Relationship of the ground heat exchanger and thermally conductive Fourier law
To energy equation;The solid area refers to: tube wall, drilling well backfilling material and the soil of ground heat exchanger;
Step 3 is calculated based on flowing heat transfer Empirical Equation Gnielinski formula combination tube fluid physical parameter
To the relevant dimensionless number Nu of fluid flowing heat transfer, fluid convection coefficient of heat transfer h is derived from according to the dimensionless number Nuf;
Step 4, the initial temperature for giving the ground heat exchanger fluid domainTemperature as current time fluidWith subsequent time fluid temperature (F.T.)Initial value;
Step 5, using fluid current time temperatureWith subsequent time temperatureConstruct Crank-Nicolson lattice
Formula is as fluid temperature (F.T.), with fluid convection coefficient of heat transfer hfConstitute the third class at the inner wall of the pipe face of the ground heat exchanger
Boundary condition, the third boundary condition refer to: calculating definite condition when conduction process.With this condition, fluid is provided
Convection transfer rate and fluid temperature (F.T.), to carry out thermally conductive solution;
Step 6 changes to energy equation obtained in the step 1 using Three-dimensional unstructured grid Finite Volume Method for Air
In generation, solves, and obtains Temperature Distribution;The step of solution is iterated to energy equation are as follows: discretization 1. is carried out to energy equation;②
Start to carry out energy after given iteration preliminary examination condition, iteration definite condition, iteration convergence condition, greatest iteration step number and time step
Measure the solution of equation;
Step 7, the step 6 every time iteratively solve after, according to the pipeline near wall of the ground heat exchanger
The temperature computation at place obtains the heat loss of unit length fluid, and subsequent time fluid temperature (F.T.) is calculated using first body heat balance method
The operation of step 8, repeating said steps 4~6, until the solid area of the ground heat exchanger calculated
Temperature field reaches stable state;
Step 9, the result obtained again to above-mentioned solution post-process, and obtain the profiling temperatures in the region.
In above method step 1, the relevant parameter for obtaining the ground heat exchanger includes: caliber, branch tube spacing, drilling
Diameter, drilling depth, circulatory mediator flow velocity and inlet temperature.
In the step 2 of the above method,
Energy conservation relation are as follows: net heat flow+body force of the increment rate of thermodynamic energy in micro unit=enter micro unit
The function that micro unit is done with surface force, wherein the micro unit, which refers to, carries out net to the solid area of the ground heat exchanger
Cell cube representated by each grid after lattice division;
The energy equation indicated with fluid specific enthalpy h and temperature T obtained are as follows:
In above formula (1), λ is the thermal coefficient of fluid, ShFor the inner heat source of fluid, Φ is since viscous effect mechanical energy turns
The part of thermal energy, referred to as dissipative function (dissipation function) are turned to, the calculating formula of the dissipative function is as follows:
In above formula (2), pdivU is surface force to fluid micro unit work done, the parameter ignorance;Simultaneously to ideal gas
Body, liquid and solid take h=cpT takes cpFor constant, dissipative function Φ is brought into source item STIn (ST=Sh+ Φ), it obtains:
Solution for soil side solid area simplifies energy equation are as follows:
In above formula (1), λ is the thermal coefficient of fluid;T is temperature;T is the time;U is vector velocity;ρ is density; CPFor
Specific heat capacity is held and (can use constant when calculating);STFor source item, due to being solved for soil side solid, these parameters
For soil.
In the step 3 of the above method, Gnielinski formula are as follows:
In above formula (5), d is pipeline nominal diameter, and l is duct length, and f is the resistance coefficient flowed in pipe, calculation formula
Are as follows:
F=(1.82lgRe-1.64)-2; (6)
For circulatory mediator in pipe are as follows:
According to Nu number definition:
Obtain the convection transfer rate h in the inner wall of the pipe face of ground heat exchangerf。
In the step 5 of the above method,
The third boundary condition in the inner wall of the pipe face of ground heat exchanger are as follows:
In above formula (9),For the temperature gradient at wall surface, TwFor wall surface temperature, TfFor tube fluid temperature, TfAre as follows:
In the step 6 of the above method, the Three-dimensional unstructured grid of the Three-dimensional unstructured grid Finite Volume Method for Air of use
For hexahedral mesh.
In the step 7 of the above method,
The thermally conductive Fourier law of thermal loss at the pipeline near wall of the ground heat exchanger are as follows:
A in above formula (10)pipe_wallFor unit pipe range wall surface area,Ti nFor the temperature of adjacent two nodes at tube wall, δ
For the mesh scale near wall;
First body heat balance method are as follows:
In above formula (11),For first body current time it is interior can, For table on this moment member body
The heat that face imports, For heat derived from this moment lower surface, For wall
Heat derived from face, For subsequent time member body it is interior can, For fluid
Mass flow, titerationFor inner iteration duration;ByDerive the intracorporal fluid temperature (F.T.) of subsequent time member
Ti n+1;Wherein, first body, which refers to, is axially divided into several pieces, list representated by every part for the fluid mass of ground heat exchanger
First body, it is preferred that every part of length takes the 1% of underground pipe buried depth.The number of division is more, and calculated result is more accurate, but therewith
The time of calculating can also increase accordingly.
In the step 9 of the above method, the result obtained to above-mentioned solution is post-processed are as follows: after calculating reaches convergence,
The temperature value of part needed for exporting generates temperature profile figure or line chart, to be analyzed.
In the above method, the step 1 arrives step 9, and the iterative solution for solving soil side conduction process is outer iteration, solves
The iterative solution of fluid temperature (F.T.) is inner iteration.
Method of the invention, due to being not necessarily to carry out grid dividing to fluid domain, without obtaining ground by solving turbulence model
Fluid flowing and heat transfer relevant parameter, realize that the turbulence model without being introduced into numerical solution can solve at pipe laying inner wall
Flowing, largely improves computational efficiency, can exchange heat between flow process in underground pipe and underground pipe soil in pipe
The numerical value Accurate Prediction of journey is a kind of efficiently and accurately test method.
The embodiment of the present invention is specifically described in further detail below.
Fig. 1 is underground pipe flowing heat transfer emulation mode flow diagram provided in an embodiment of the present invention, the emulation mode packet
Include following steps:
Step 1, the relevant parameter of the ground heat exchanger is obtained;
Step 2, first according to heat transfer theory, select by solid area applied energy law of conservation and thermally conductive
Fourier law obtains energy equation;
Choose law of conservation of energy:
According to Conservation Relationship: [increment rate of thermodynamic energy in micro unit]=[into the net heat flow of micro unit]+[volume
The function that power and surface force do micro unit], it can obtain the energy equation indicated with fluid specific enthalpy h and temperature T:
Wherein, λ is the thermal coefficient of fluid, ShFor the inner heat source of fluid, Φ is since viscous effect mechanical energy is converted into
The part of thermal energy, referred to as dissipative function (dissipation function), calculating formula is as follows:
Wherein, pdivU is surface force to fluid micro unit work done, can generally be ignored;Simultaneously to perfect gas, liquid
Body and solid can take h=cpT further takes cpFor constant, and dissipative function Φ is brought into source item STIn (ST=Sh+ Φ),
It obtains:
Solution for soil side solid area can be further simplified energy equation:
Step 3 is calculated based on flowing heat transfer Empirical Equation Gnielinski formula combination tube fluid physical parameter
To the relevant dimensionless number Nu of fluid flowing heat transfer, fluid convection coefficient of heat transfer h is derived to obtainf;
In this step, the flowing heat transfer state of intraductal turbulance is characterized using dimensionless number Nu, Nu number is by intraductal turbulance
Flowing heat transfer correlation Gnielinski formula calculated, Gnielinski formula:
Wherein, d is pipeline nominal diameter, and l is duct length, and f is the resistance coefficient flowed in pipe, calculation formula are as follows:
F=(1.82lgRe-1.64)-2(6)
Particularly, for circulatory mediator in pipe,
According to Nu number definition:
The convection transfer rate h of available underground pipe inner wallf;
Step 4, given fluid domain initial temperatureTemperature as current time fluidWith subsequent time fluid temperature (F.T.)Initial value;
In this step, by grid of the method without fluid domain stated in the present invention, the fluid temperature (F.T.) is imagination
Fluid temperature (F.T.), statement is in this way, the temperature of given fluid domain is so as to subsequent calculating progress below;
Step 5, using fluid current time temperatureWith subsequent time temperatureConstruct Crank-Nicolson format
As fluid temperature (F.T.), with fluid convection coefficient of heat transfer hfConstitute the third boundary condition at ground heat exchanger inner wall;
The third boundary condition in inner wall of the pipe face may be expressed as:
Wherein,For the temperature gradient at wall surface, TwFor wall surface temperature, TfIt, specifically can be with for tube fluid temperature
Expression are as follows:
Step 6 solves energy equation obtained in step 1 using Three-dimensional unstructured grid Finite Volume Method for Air,
Obtain Temperature Distribution;
Three-dimensional unstructured grid employed in the step 5 is hexahedral mesh, to guarantee the convergence of calculating process
Property and stability;
Step 7, step 5 every time iteratively solve after, according to the temperature near wall pass through thermally conductive Fourier law
The heat loss of unit length fluid is calculated, subsequent time fluid temperature (F.T.) is calculated using first body heat balance method
The statement of the thermally conductive Fourier law of thermal loss near wall are as follows:
Wherein Apipe_wallFor unit pipe range wall surface area,Ti nFor the temperature of adjacent two nodes at tube wall, δ is nearly wall
Mesh scale at face;
The statement of first body heat balance method are as follows:
Wherein,For first body current time it is interior can, For the importing of this moment member body upper surface
Heat, For heat derived from this moment lower surface, To be led at wall surface
Heat out, For subsequent time member body it is interior can, For the quality of fluid
Flow, titerationFor inner iteration duration;
ByIt can derive the intracorporal fluid temperature (F.T.) T of subsequent time memberi n+1。
Step 8, the operation for repeating step 4~6, until the temperature field of the zoning reaches stable state;
Step 9, the result obtained again to above-mentioned solution post-process, and obtain the profiling temperatures in the region.
In specific implementation, above-mentioned steps 1 arrive step 9, and the iterative solution for solving soil side conduction process is outer iteration, solve
The iterative solution of fluid temperature (F.T.) is inner iteration.
9 steps are to the heat transfer between underground pipe soil during solar cross-season soil thermal storage described in the present embodiment
Process is simulated, and the specific size and physical parameter for the ground heat exchanger simulated are as shown in Table 1 and Table 2.
The underground pipe outlet temperature and turbulence model calculated result that the emulation mode that Fig. 2 show the embodiment of the present invention obtains
Comparison schematic diagram, Fig. 2 compared the double-U-shaped ground heat exchanger unilateral side U-tube being calculated using the embodiment of the present invention and gone out
The result that fluid temperature (F.T.) at mouthful is calculated with the traditional analog method using turbulence model Standard k-e model.Such as figure
Shown in 2, the result and turbulence model emulated in the embodiment of the present invention is coincide preferably, is demonstrated the present invention and is established emulation side
The accuracy of method.
What Fig. 3 showed that emulation mode in the embodiment of the present invention obtains obtains along journey fluid temperature (F.T.) and conventional numeric calculation method
To result comparison schematic diagram, Fig. 3 compare that different emulation modes obtain along journey fluid temperature (F.T.), as shown in figure 3, the present invention is implemented
It is 0.025% that emulation mode, which obtains result and the maximum deviation of conventional numeric calculation method acquisition result, in example, and average deviation is
0.02%, situation of coincideing between the two is good, further demonstrates the accuracy of the established emulation mode of the present invention.
The content being not described in detail in the embodiment of the present invention belongs to the prior art well known to professional and technical personnel in the field.
In conclusion the method for the embodiment of the present invention using numerical simulation to the flowing heat transfer process of ground heat exchanger into
Row emulation, it may be convenient to the temperature of buried tube fluid and the Temperature Distribution of soil are obtained, without carrying out net to fluid domain
Lattice are divided and are calculated, and can be greatly reduced the calculation amount of solar cross-season soil thermal storage, be shortened the R&D cycle.
The foregoing is only a preferred embodiment of the present invention, but scope of protection of the present invention is not limited thereto,
Within the technical scope of the present disclosure, any changes or substitutions that can be easily thought of by anyone skilled in the art,
It should be covered by the protection scope of the present invention.Therefore, protection scope of the present invention should be with the protection model of claims
Subject to enclosing.
Claims (10)
1. a kind of emulation mode of solar cross-season underground pipe accumulation of heat, which is characterized in that for the storage to ground heat exchanger
Heat carries out simulation calculation, comprising the following steps:
Step 1, the relevant parameter of the ground heat exchanger is obtained;
Step 2, applied energy Conservation Relationship and thermally conductive Fourier law obtain the energy of the solid area of the ground heat exchanger
Measure equation;
Fluid is calculated with flowing heat transfer Empirical Equation Gnielinski formula combination tube fluid physical parameter in step 3
The relevant dimensionless number Nu of flowing heat transfer is derived from fluid convection coefficient of heat transfer h according to the dimensionless number Nuf;
Step 4, the initial temperature for giving the ground heat exchanger fluid domainTemperature as current time fluidWith under
One moment fluid temperature (F.T.)Initial value;
Step 5, using fluid current time temperatureWith subsequent time temperatureConstruct the conduct of Crank-Nicolson format
Fluid temperature (F.T.), with fluid convection coefficient of heat transfer hfConstitute the third class perimeter strip at the inner wall of the pipe face of the ground heat exchanger
Part;
Step 6 is iterated energy equation obtained in the step 1 using Three-dimensional unstructured grid Finite Volume Method for Air and asks
Solution, obtains Temperature Distribution;
Step 7, the step 6 every time iteratively solve after, at the pipeline near wall of the ground heat exchanger
Temperature computation obtains the heat loss of unit length fluid, and subsequent time fluid temperature (F.T.) is calculated using first body heat balance method
The operation of step 8, repeating said steps 4~7, until the temperature field of the solid area of the ground heat exchanger of calculating
Reach stable state;
Step 9, the result obtained again to above-mentioned solution post-process, and obtain the profiling temperatures in the region.
2. the emulation mode of solar cross-season underground pipe accumulation of heat according to claim 1, which is characterized in that the method
In step 1, the relevant parameter for obtaining the ground heat exchanger includes: caliber, branch tube spacing, bore diameter, drilling depth, follows
Ring velocity of medium and inlet temperature.
3. the emulation mode of solar cross-season underground pipe accumulation of heat according to claim 1 or 2, which is characterized in that described
In the step 2 of method,
Energy conservation relation are as follows: the net heat flow+body force and table of the increment rate of thermodynamic energy in micro unit=enter micro unit
The function that face power does micro unit, wherein the micro unit refers to that carrying out grid to the solid area of the ground heat exchanger draws
Cell cube representated by each grid after point;
The energy equation obtained are as follows:
In above formula (4), λ is the thermal coefficient of fluid;T is temperature;T is the time;U is vector velocity;ρ is density;CP is soil
Specific heat capacity;ST is source item.
4. the emulation mode of solar cross-season underground pipe accumulation of heat according to claim 1 or 2, which is characterized in that described
In the step 3 of method, Gnielinski formula are as follows:
In above formula (5), d is pipeline nominal diameter, and l is duct length, and f is the resistance coefficient flowed in pipe, f calculation formula are as follows: f
=(1.82lgRe-1.64)-2; (6)
For circulatory mediator in pipe are as follows:
According to dimensionless number Nu number definition:
Obtain the convection transfer rate h in the inner wall of the pipe face of ground heat exchangerf。
5. the emulation mode of solar cross-season underground pipe accumulation of heat according to claim 1 or 2, which is characterized in that described
In the step 5 of method,
The third boundary condition in the inner wall of the pipe face of ground heat exchanger are as follows:
In above formula (9),For the temperature gradient at wall surface, TwFor wall surface temperature, TfFor tube fluid temperature, TfAre as follows:
6. the emulation mode of solar cross-season underground pipe accumulation of heat according to claim 1 or 2, which is characterized in that described
In the step 6 of method, the Three-dimensional unstructured grid of the Three-dimensional unstructured grid Finite Volume Method for Air used is hexahedral mesh.
7. the emulation mode of solar cross-season underground pipe accumulation of heat according to claim 1 or 2, which is characterized in that described
In the step 7 of method,
The thermally conductive Fourier law of thermal loss at the pipeline near wall of the ground heat exchanger are as follows:
A in above formula (10)pipe_wallFor unit pipe range wall surface area,Ti nFor the temperature of adjacent two nodes at tube wall, δ is close
Mesh scale at wall surface;
First body heat balance method are as follows:
In above formula (11),For first body current time it is interior can, It is led for this moment member body upper surface
The heat entered, For heat derived from this moment lower surface, At wall surface
Derived heat, For subsequent time member body it is interior can, For the quality of fluid
Flow, titerationFor inner iteration duration;ByDerive the intracorporal fluid temperature (F.T.) T of subsequent time memberi n+1;Its
In, the member body, which refers to, is axially divided into several pieces, cell cube representated by every part for the fluid mass of ground heat exchanger.
8. the emulation mode of solar cross-season underground pipe accumulation of heat according to claim 1 or 2, which is characterized in that described
In the step 5 of method, the third boundary condition is definite condition when calculating conduction process.
9. the emulation mode of solar cross-season underground pipe accumulation of heat according to claim 1 or 2, which is characterized in that described
In the step 6 of method, the iterative solution step to energy equation includes:
(1) discretization is carried out to energy equation;
(2) it is opened after giving iteration preliminary examination condition, iteration definite condition, iteration convergence condition, greatest iteration step number and time step
Beginning carries out the solution of energy equation.
10. the emulation mode of solar cross-season underground pipe accumulation of heat according to claim 1 or 2, which is characterized in that described
In method, the step 1 arrives step 9, and the iterative solution for solving soil side conduction process is outer iteration, solves changing for fluid temperature (F.T.)
In generation, is solved to inner iteration.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910193763.5A CN109992846B (en) | 2019-03-14 | 2019-03-14 | Simulation method for solar cross-season buried pipe heat storage |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910193763.5A CN109992846B (en) | 2019-03-14 | 2019-03-14 | Simulation method for solar cross-season buried pipe heat storage |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109992846A true CN109992846A (en) | 2019-07-09 |
CN109992846B CN109992846B (en) | 2022-12-09 |
Family
ID=67130552
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910193763.5A Active CN109992846B (en) | 2019-03-14 | 2019-03-14 | Simulation method for solar cross-season buried pipe heat storage |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109992846B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110489912A (en) * | 2019-08-27 | 2019-11-22 | 北京石油化工学院 | A kind of method of solar cross-season soil thermal storage hierarchy slicing numerical simulation |
CN110610037A (en) * | 2019-08-30 | 2019-12-24 | 中国科学院工程热物理研究所 | Thermodynamic simulation method for counter-flow heat regenerator |
CN110991061A (en) * | 2019-12-10 | 2020-04-10 | 国网浙江省电力有限公司杭州供电公司 | Thermoelectric comprehensive energy simulation planning method |
CN111551288A (en) * | 2020-05-13 | 2020-08-18 | 山东建筑大学 | Prediction method for circulating liquid temperature of middle-deep layer U-shaped heat exchange well |
CN114357838A (en) * | 2022-01-07 | 2022-04-15 | 西安交通大学 | Simulation method of coaxial sleeve type buried pipe heat exchanger with variable flow and variable pipe diameter in different seasons |
CN116562184A (en) * | 2023-04-18 | 2023-08-08 | 中南大学 | Rapid prediction calculation method for mold temperature field in autoclave molding process |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106096293A (en) * | 2016-06-17 | 2016-11-09 | 北京航空航天大学 | A kind of big rotation composite extending arm heat-driven oscillation Forecasting Methodology |
CN106294963A (en) * | 2016-08-04 | 2017-01-04 | 国网上海市电力公司 | Direct-buried cable carrying current calculation method |
JP2017015364A (en) * | 2015-07-06 | 2017-01-19 | 株式会社 トラストプラン | Solar heat underground heat storage snow-melting system and its control method |
CN107192571A (en) * | 2017-06-15 | 2017-09-22 | 西安交通大学 | Horizontal tube bundle and the outer refrigeration working medium falling film evaporation of single tube, pool boiling and condensation phase change heat exchange test device |
-
2019
- 2019-03-14 CN CN201910193763.5A patent/CN109992846B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017015364A (en) * | 2015-07-06 | 2017-01-19 | 株式会社 トラストプラン | Solar heat underground heat storage snow-melting system and its control method |
CN106096293A (en) * | 2016-06-17 | 2016-11-09 | 北京航空航天大学 | A kind of big rotation composite extending arm heat-driven oscillation Forecasting Methodology |
CN106294963A (en) * | 2016-08-04 | 2017-01-04 | 国网上海市电力公司 | Direct-buried cable carrying current calculation method |
CN107192571A (en) * | 2017-06-15 | 2017-09-22 | 西安交通大学 | Horizontal tube bundle and the outer refrigeration working medium falling film evaporation of single tube, pool boiling and condensation phase change heat exchange test device |
Non-Patent Citations (2)
Title |
---|
MÁRCIO SANTOS,ET AL: "Off-design modelling of an organic Rankine cycle micro-CHP: Modular", 《APPLIED THERMAL ENGINEERING》 * |
孙中慧: "太阳能跨季节储热地埋管换热器数值模型适用性研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110489912A (en) * | 2019-08-27 | 2019-11-22 | 北京石油化工学院 | A kind of method of solar cross-season soil thermal storage hierarchy slicing numerical simulation |
CN110489912B (en) * | 2019-08-27 | 2022-12-13 | 北京石油化工学院 | Method for simulating numerical value of solar cross-season soil heat storage layered slice |
CN110610037A (en) * | 2019-08-30 | 2019-12-24 | 中国科学院工程热物理研究所 | Thermodynamic simulation method for counter-flow heat regenerator |
CN110610037B (en) * | 2019-08-30 | 2023-04-18 | 中国科学院工程热物理研究所 | Thermodynamic simulation method for counter-flow heat regenerator |
CN110991061A (en) * | 2019-12-10 | 2020-04-10 | 国网浙江省电力有限公司杭州供电公司 | Thermoelectric comprehensive energy simulation planning method |
CN110991061B (en) * | 2019-12-10 | 2023-07-07 | 国网浙江省电力有限公司杭州供电公司 | Thermoelectric comprehensive energy simulation planning method |
CN111551288A (en) * | 2020-05-13 | 2020-08-18 | 山东建筑大学 | Prediction method for circulating liquid temperature of middle-deep layer U-shaped heat exchange well |
CN114357838A (en) * | 2022-01-07 | 2022-04-15 | 西安交通大学 | Simulation method of coaxial sleeve type buried pipe heat exchanger with variable flow and variable pipe diameter in different seasons |
CN114357838B (en) * | 2022-01-07 | 2024-01-16 | 西安交通大学 | Simulation method of coaxial sleeve type buried pipe heat exchanger with variable-flow pipe diameter |
CN116562184A (en) * | 2023-04-18 | 2023-08-08 | 中南大学 | Rapid prediction calculation method for mold temperature field in autoclave molding process |
CN116562184B (en) * | 2023-04-18 | 2024-02-13 | 中南大学 | Rapid prediction calculation method for mold temperature field in autoclave molding process |
Also Published As
Publication number | Publication date |
---|---|
CN109992846B (en) | 2022-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109992846A (en) | The emulation mode of solar cross-season underground pipe accumulation of heat | |
Eslami-Nejad et al. | Coupling of geothermal heat pumps with thermal solar collectors using double U-tube boreholes with two independent circuits | |
CN105608273B (en) | A kind of system optimizing power battery pack Temperature Distribution based on CFD software | |
Gao et al. | Numerical and experimental investigation on thermal stratification characteristics affected by the baffle plate in thermal storage tank | |
CN107947182B (en) | Dynamic power flow analysis method and dynamic power flow analysis system | |
CN109492325A (en) | A kind of multipotency coupled system tidal current analysis method based on extension energy hub | |
CN106647329A (en) | Heat supply network equivalent modeling method | |
Bouhal et al. | Towards an energy efficiency optimization of solar horizontal storage tanks and circulation pipes integrating evacuated tube collectors through CFD parametric studies | |
CN107025334A (en) | Central heating system heat user integrated dynamic model method for building up and device | |
CN110220405A (en) | Solid heat storage heat transfer control method based on Fourier number consistency | |
CN104732111B (en) | A kind of earth source heat pump it is efficient by when method for numerical simulation | |
CN109028609A (en) | A kind of flat-plate solar heat collector array strings arranged in parallel design method | |
Huang et al. | Optimization of a collector-storage solar air heating system for building heat recovery ventilation preheating in the cold area | |
Afshari | Experimental and numerical study on heat transfer enhancement of home radiators by employing solar cells and fans | |
CN110489912B (en) | Method for simulating numerical value of solar cross-season soil heat storage layered slice | |
CN110781598B (en) | Method for calculating running state of regional comprehensive energy system | |
Abbas et al. | Performance analysis of seasonal soil heat storage system based on numerical simulation and experimental investigation | |
CN112270109A (en) | Method for simulating heating performance of graphite rod in high-temperature carbonization furnace | |
CN111783309A (en) | Dynamic simulation method of steam heating network based on internal conservation | |
CN206609036U (en) | A kind of heating system of solar energy and burnt gas wall hanging furnace coupling application | |
CN206625899U (en) | A kind of distributed natural gas CCHP efficiency accounting system | |
CN203299154U (en) | Variable-frequency heat pump type rock-soil thermophysical property tester | |
CN208859705U (en) | A kind of integral type heating energy-storage system | |
CN104748996B (en) | Power plant's smoke discharging residual heat utilizes system Dynamic Response Simulation experimental apparatus for testing and method | |
CN203981418U (en) | Fibre Optical Sensor becomes power thermal response test instrument |
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 |