[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN102663243A - Numerical simulation method of buried tube temperature field of ground source heat pump under thermoosmosis coupling - Google Patents

Numerical simulation method of buried tube temperature field of ground source heat pump under thermoosmosis coupling Download PDF

Info

Publication number
CN102663243A
CN102663243A CN2012100892283A CN201210089228A CN102663243A CN 102663243 A CN102663243 A CN 102663243A CN 2012100892283 A CN2012100892283 A CN 2012100892283A CN 201210089228 A CN201210089228 A CN 201210089228A CN 102663243 A CN102663243 A CN 102663243A
Authority
CN
China
Prior art keywords
over
change step
change
temperature field
temperature
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
Application number
CN2012100892283A
Other languages
Chinese (zh)
Other versions
CN102663243B (en
Inventor
钱华
张磊
郑晓红
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Hongtai Energy Technology Co., Ltd.
Original Assignee
Changshu Nanjing Normal University Development Research Institute Co Ltd
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 Changshu Nanjing Normal University Development Research Institute Co Ltd filed Critical Changshu Nanjing Normal University Development Research Institute Co Ltd
Priority to CN201210089228.3A priority Critical patent/CN102663243B/en
Publication of CN102663243A publication Critical patent/CN102663243A/en
Application granted granted Critical
Publication of CN102663243B publication Critical patent/CN102663243B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to a numerical simulation method of buried tube temperature field of a ground source heat pump under thermoosmosis coupling. The method is used for numerical simulation of the buried tube temperature field of the ground source heat pump under thermoosmosis coupling. The heat exchange between the soil and the ground water is regarded as an internal heat source in an innovation manner, and an internal heat source item in an equation is simplified. The tube groups model condition, rock thermal performance data and iterative step r are initially input; implicit iterative computations are employed; and further, the temperature values of the soil around are obtained under the circumstance of a set seepage velocity and soil parameters. According to the method, changes of temperature distribution of the soil around the large-scale tube groups along with the seepage can be easily obtained; and the method has a low requirement for computer hardware and strong generality, effectively shortens the operation time and enlarges the scale of the calculated tube groups.

Description

Heat is oozed ground buried pipe of ground source heat pump temperature field method for numerical simulation under the coupling
Technical field
The invention belongs to ground source heat pump technology and use and energy-saving field, relate in particular to definite method of earth source heat pump pipe laying crowd's soil moisture field under the situation of considering the influence of soil seepage flow.
Background technology
From some actual engineerings; The seepage action of ground water heat transfer characteristic of pipe laying over the ground has considerable influence, can increase the coefficient of heat transfer of soil and circulating fluid, makes long by the pipe laying length that pure conduction model designed of not considering the seepage flow situation; The waste resource increases cost.Therefore must consider the seepage action of ground water situation during earth-source hot-pump system conceptual design.And for extensive nest of tubes situation, because computer capacity is bigger, general business software simulated time is oversize, is unfavorable for actual engineering reference.The present invention can be after relatively short operation time of this significant improvement; Realize pipe laying and ooze under the coupling situation in heat on every side; Through the implicit expression iterative numerical approach; Calculate the situation that extensive nest of tubes surrounding soil Temperature Distribution changes with seepage flow, for practical engineering design and research provide corresponding reference frame.
Therefore, the present invention proposes a kind of heat and oozes numerical simulation new method in ground buried pipe of ground source heat pump temperature field under the coupling.
Summary of the invention
The purpose of this invention is to provide a kind of heat and ooze ground buried pipe of ground source heat pump temperature field method for numerical simulation under the coupling; The present invention can obtain extensive nest of tubes surrounding soil Temperature Distribution easily with the seepage flow situation of change; Effectively shorten operation time, and enlarge the scale of calculating nest of tubes.
For solving the problems of the technologies described above, technical scheme of the present invention is following:
A kind of heat is oozed ground buried pipe of ground source heat pump temperature field method for numerical simulation under the coupling, it is characterized in that may further comprise the steps:
Step 1: input model condition and ground thermal property data are as the initial parameter of algorithm;
Step 2: the correlation parameter of initialization nest of tubes model and soil comprises: the parameters of mentioning in the step 1;
Step 3: loop initialization pointer τ=1, △ τ=r, r are the change step of time τ, and the variation range of system operation time τ is 1~t, and t moves closing time for the system that sets;
Step 4: judge τ≤t, if then change step 5 over to; If, then do not change step 8 over to;
Step 5: call the heat flux subroutine, calculate the distribution situation of heat flow field;
Step 6: call the temperature interative routine, calculate the distribution situation in temperature field;
Step 7: change τ working time, carry out τ=τ+r, get back to step 4, get into next subcycle;
Step 8: output temperature field evaluation result, and then obtain and set under percolation flow velocity and the soil parameters situation soil temperature value everywhere.
Wherein condition and ground thermal property data comprise that nest of tubes arrangement mode, nest of tubes quantity, the setting of nest of tubes boundary condition, soil initial temperature are t 0, soil thermal conductivity is that λ, specific heat capacity hold to hold for c, specific heat capacity and be c w, soil density is that ρ, constant heat flux value are q, soil moisture content ω and seepage action of ground water speed u.
The step of subroutine heat flux subroutine is following:
Steps A: heat flow field initialization, each grid node place initial assignment are zero;
Step B: loop initialization pointer i=1, the change step of i is 1;
Step C: loop initialization pointer j=1, the change step of j is 1;
Step D: judge j≤N,, if then change step e over to; If not, then change step G over to, N value representation longitudinal grid number obtains during by the heat flow field initialization;
Step e: calculate to judge whether this point is the pipe laying node through computation model, if then change step F over to; If not, then carry out j=j+1, change step D over to;
Step F: record node assignment, carry out j=j+1, change step D over to;
Step G: carry out i=i+1;
Step H: judge i≤M,, if then change step I over to; If not, then change step C over to, M value representation transverse grid number obtains during by the heat flow field initialization;
Step I: the heat flow field data transfer is returned to the group program.
The step of subroutine temperature interative routine is following:
Step a: the temperature field initialization, each grid node place initial assignment is initial ground temperature t 0
Step b: loop initialization pointer i=1, the change step of i is 1;
Step c: judge i≤M,, if then change step 4 over to; If not, then change step h over to, M value representation transverse grid number obtains during by the temperature field initialization;
Steps d: loop initialization pointer j=1, the change step of j is 1;
Step e: judge j≤N,, if then change step f over to; If not, then carry out i=i+1, change step c over to, N value representation longitudinal grid number obtains during by the temperature field initialization;
Step f: calculate this some place temperature value through computation model, Model Calculation can obtain through following formula:
Figure 409754DEST_PATH_IMAGE001
Step g: (τ+r/2) is the temperature field data constantly, carry out j=j+1 then, change step e over to for record;
Step h: loop initialization pointer j=1, the change step of j is 1;
Step I: judge j≤N,, if then change step j over to; If not, then change step n over to, N value representation longitudinal grid number obtains during by the temperature field initialization;
Step j: loop initialization pointer i=1, the change step of i is 1;
Step k: judge i≤M,, if then change step l over to; If not, then carry out j=j+1, change step I over to, M value representation transverse grid number obtains during by the temperature field initialization;
Step l: calculate this some place temperature value through computation model, carry out i=i+1 then, change step j over to; Model Calculation can obtain through following formula:
Figure 670971DEST_PATH_IMAGE002
Step m: (τ+r) is the temperature field data constantly, carry out i=i+1 then, change step k over to for record;
Step n: the temperature field data transfer is returned to master routine.
The present invention proposes a kind of heat and oozes ground buried pipe of ground source heat pump temperature field method for numerical simulation under the coupling, and its program mainly comprises three parts: main working procedure, heat flow field counting subroutine and temperature iterative computation subroutine.
The operation master routine, input model condition and ground thermal property data comprise: the nest of tubes arrangement mode, nest of tubes quantity, the nest of tubes boundary condition is provided with, and the soil initial temperature is t 0, soil thermal conductivity is λ, and it is c that specific heat capacity is held, and it is c that specific heat capacity is held w, soil density is ρ, the constant heat flux value is q, and soil moisture content w, seepage action of ground water speed u, system operation times etc. are as the initial parameter of algorithm; The nest of tubes model, heat flow field and the temperature field that generate in the initialize routine; The time step of setting program iteration, and loop initialization pointer; Call the heat flow field subroutine earlier, calculate to judge through computation model whether this point is the pipe laying node, if, then with hot-fluid parameter assignment in this node, if not, then with zero assignment in this point, obtain the heat flow field distribution situation, and return master routine; Heat flow field and last iteration constantly must be arrived the temperature field to superpose; Then, call the temperature interative routine, utilization formula (1) is carried out the transverse grid implicit iterative, calculates that (τ+r/2) is the temperature field data constantly
Figure 339850DEST_PATH_IMAGE001
(1)
Utilization formula (2) is carried out the longitudinal grid implicit iterative, calculates that (τ+r) is the temperature field data constantly
Figure 333214DEST_PATH_IMAGE002
(2)
Obtain the distribution situation in temperature field, and return master routine; Move to predetermined finish time of output temperature field evaluation result, and then obtain and set under percolation flow velocity and the soil parameters situation soil temperature value everywhere.
Advantage of the present invention is following.
(1) having proposed a kind of new algorithm influences the computing velocity of soil moisture field fast to seepage flow.It is following that two-way implicit algorithm is calculated in the following temperature field of seepage flow influence:
Horizontal iterative computation:
Vertical iterative computation:
Figure 637604DEST_PATH_IMAGE002
It is thus clear that definite needs of soil moisture field are confirmed iteration time step-length r; The present invention adopts the implicit expression computing method, and the gained result is insensitive to the value variable effect of time step r.Earlier given longitudinal grid variable j changes the transverse grid variable i and in scope separately, travels through, and each class value calculates the (temperature value that τ+r/2) is corresponding constantly; Then given transverse grid variable i changes longitudinal grid variable j and in scope separately, travels through, and each class value calculates the (temperature value that τ+r) is corresponding constantly.
(2) this confirms that method oozes to heat that ground buried pipe of ground source heat pump temperature field method for numerical simulation provides a kind of new algorithm under the coupling; Soil and underground water heat are regarded as inner heat exchange; Simplify the thermal source item; Can calculate the soil moisture field under the seepage flow influence at short notice fast, for the design and the research of ground buried pipe of ground source heat pump provides certain reference frame.
Description of drawings
Fig. 1 is that heat of the present invention is oozed ground buried pipe of ground source heat pump temperature field numerical computation method main program block diagram under the coupling.
Fig. 2 is a hot-fluid assignment computing block diagram in pipe laying place in the nest of tubes zone.
Fig. 3 is a temperature field numerical value iterative computation block diagram in the nest of tubes zone.
Embodiment
A kind of heat is oozed ground buried pipe of ground source heat pump temperature field method for numerical simulation under the coupling.
The master routine operation:
Step 1: input model condition and ground thermal property data comprise: the nest of tubes arrangement mode, and nest of tubes quantity, the nest of tubes boundary condition is provided with, and the soil initial temperature is t 0, soil thermal conductivity is λ, and it is c that specific heat capacity is held, and it is c that specific heat capacity is held w, soil density is ρ, the constant heat flux value is q, and soil moisture content ω, seepage action of ground water speed u is as the initial parameter of algorithm;
Step 2: the correlation parameter of initialization nest of tubes model and soil comprises: the parameters of mentioning in the step 1;
Step 3: loop initialization pointer τ=1, △ τ=r, r are the change step of time τ, and the variation range of system operation time τ is 1~t, and t moves closing time for the system that sets;
Step 4: judge τ≤t, if then change step 5 over to; If, then do not change step 8 over to;
Step 5: call the heat flux subroutine, calculate the distribution situation of heat flow field;
Step 6: call the temperature interative routine, calculate the distribution situation in temperature field;
Step 7: change τ working time, carry out τ=τ+r, get back to step 4, get into next subcycle;
Step 8: output temperature field evaluation result, and then obtain and set under percolation flow velocity and the soil parameters situation soil temperature value everywhere.
Hot-fluid subroutine call operation:
Step 1: heat flow field initialization, each grid node place initial assignment are zero;
Step 2: loop initialization pointer i=1, the change step of i is 1;
Step 3: loop initialization pointer j=1, the change step of j is 1;
Step 4: judge j≤N,, if then change step 5 over to; If not, then change step 7 over to, N value representation longitudinal grid number obtains during by the heat flow field initialization;
Step 5: calculate to judge whether this point is the pipe laying node through computation model, if then change step 6 over to; If not, then carry out j=j+1, change step 4 over to;
Step 6: record node assignment, carry out j=j+1, change step 4 over to;
Step 7: carry out i=i+1;
Step 8: judge i≤M,, if then change step 9 over to; If not, then change step 3 over to, M value representation transverse grid number obtains during by the heat flow field initialization;
Step 9: the heat flow field data transfer is returned to the group program.
The temperature interative routine calls operation:
Step 1: the temperature field initialization, each grid node place initial assignment is initial ground temperature t 0
Step 2: loop initialization pointer i=1, the change step of i is 1;
Step 3: judge i≤M,, if then change step 4 over to; If not, then change step 8 over to, M value representation transverse grid number obtains during by the temperature field initialization;
Step 4: loop initialization pointer j=1, the change step of j is 1;
Step 5: judge j≤N,, if then change step 6 over to; If not, then carry out i=i+1, change step 3 over to, N value representation longitudinal grid number obtains during by the temperature field initialization;
Step 6: calculate this some place temperature value through computation model, Model Calculation can obtain through following formula:
Figure 160990DEST_PATH_IMAGE001
Step 7: (τ+r/2) is the temperature field data constantly, carry out j=j+1 then, change step 5 over to for record;
Step 8: loop initialization pointer j=1, the change step of j is 1;
Step 9: judge j≤N,, if then change step 10 over to; If not, then change step 14 over to, N value representation longitudinal grid number obtains during by the temperature field initialization;
Step 10: loop initialization pointer i=1, the change step of i is 1;
Step 11: judge i≤M,, if then change step 12 over to; If not, then carry out j=j+1, change step 9 over to, M value representation transverse grid number obtains during by the temperature field initialization;
Step 12: calculate this some place temperature value through computation model, carry out i=i+1 then, change step 10 over to; Model Calculation can obtain through following formula:
Figure 590834DEST_PATH_IMAGE002
Step 13: (τ+r) is the temperature field data constantly, carry out i=i+1 then, change step 11 over to for record;
Step 14: the temperature field data transfer is returned to master routine.

Claims (4)

1. a heat is oozed ground buried pipe of ground source heat pump temperature field method for numerical simulation under the coupling, it is characterized in that may further comprise the steps:
Step 1: input model condition and ground thermal property data are as the initial parameter of algorithm;
Step 2: the correlation parameter of initialization nest of tubes model and soil comprises: the parameters of mentioning in the step 1;
Step 3: loop initialization pointer τ=1, △ τ=r, r are the change step of time τ, and the variation range of system operation time τ is 1~t, and t moves closing time for the system that sets;
Step 4: judge τ≤t, if then change step 5 over to; If, then do not change step 8 over to;
Step 5: call the heat flux subroutine, calculate the distribution situation of heat flow field;
Step 6: call the temperature interative routine, calculate the distribution situation in temperature field;
Step 7: change τ working time, carry out τ=τ+r, get back to step 4, get into next subcycle;
Step 8: output temperature field evaluation result, and then obtain and set under percolation flow velocity and the soil parameters situation soil temperature value everywhere.
2. heat according to claim 1 is oozed ground buried pipe of ground source heat pump temperature field method for numerical simulation under the coupling, it is characterized in that: said step 1 model condition and ground thermal property data comprise that nest of tubes arrangement mode, nest of tubes quantity, the setting of nest of tubes boundary condition, soil initial temperature are t 0, soil thermal conductivity is that λ, specific heat capacity hold to hold for c, specific heat capacity and be c w, soil density is that ρ, constant heat flux value are q, soil moisture content ω and seepage action of ground water speed u.
3. heat according to claim 1 is oozed ground buried pipe of ground source heat pump temperature field method for numerical simulation under the coupling, it is characterized in that the step of said heat flux subroutine is following:
Steps A: heat flow field initialization, each grid node place initial assignment are zero;
Step B: loop initialization pointer i=1, the change step of i is 1;
Step C: loop initialization pointer j=1, the change step of j is 1;
Step D: judge j≤N,, if then change step e over to; If not, then change step G over to, N value representation longitudinal grid number obtains during by the heat flow field initialization;
Step e: calculate to judge whether this point is the pipe laying node through computation model, if then change step F over to; If not, then carry out j=j+1, change step D over to;
Step F: record node assignment, carry out j=j+1, change step D over to;
Step G: carry out i=i+1;
Step H: judge i≤M,, if then change step I over to; If not, then change step C over to, M value representation transverse grid number obtains during by the heat flow field initialization;
Step I: the heat flow field data transfer is returned to the group program.
4. heat according to claim 1 is oozed ground buried pipe of ground source heat pump temperature field method for numerical simulation under the coupling, it is characterized in that the step of said temperature interative routine is following:
Step a: the temperature field initialization, each grid node place initial assignment is initial ground temperature t 0
Step b: loop initialization pointer i=1, the change step of i is 1;
Step c: judge i≤M,, if then change step 4 over to; If not, then change step h over to, M value representation transverse grid number obtains during by the temperature field initialization;
Steps d: loop initialization pointer j=1, the change step of j is 1;
Step e: judge j≤N,, if then change step f over to; If not, then carry out i=i+1, change step c over to, N value representation longitudinal grid number obtains during by the temperature field initialization;
Step f: calculate this some place temperature value through computation model, Model Calculation can obtain through following formula:
Figure 477174DEST_PATH_IMAGE001
Step g: (τ+r/2) is the temperature field data constantly, carry out j=j+1 then, change step e over to for record;
Step h: loop initialization pointer j=1, the change step of j is 1;
Step I: judge j≤N,, if then change step j over to; If not, then change step n over to, N value representation longitudinal grid number obtains during by the temperature field initialization;
Step j: loop initialization pointer i=1, the change step of i is 1;
Step k: judge i≤M,, if then change step l over to; If not, then carry out j=j+1, change step I over to, M value representation transverse grid number obtains during by the temperature field initialization;
Step l: calculate this some place temperature value through computation model, carry out i=i+1 then, change step j over to; Model Calculation can obtain through following formula:
Figure 2012100892283100001DEST_PATH_IMAGE002
Step m: (τ+r) is the temperature field data constantly, carry out i=i+1 then, change step k over to for record;
Step n: the temperature field data transfer is returned to master routine.
CN201210089228.3A 2012-03-30 2012-03-30 Numerical simulation method of buried tube temperature field of ground source heat pump under thermoosmosis coupling Expired - Fee Related CN102663243B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210089228.3A CN102663243B (en) 2012-03-30 2012-03-30 Numerical simulation method of buried tube temperature field of ground source heat pump under thermoosmosis coupling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210089228.3A CN102663243B (en) 2012-03-30 2012-03-30 Numerical simulation method of buried tube temperature field of ground source heat pump under thermoosmosis coupling

Publications (2)

Publication Number Publication Date
CN102663243A true CN102663243A (en) 2012-09-12
CN102663243B CN102663243B (en) 2015-01-07

Family

ID=46772734

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210089228.3A Expired - Fee Related CN102663243B (en) 2012-03-30 2012-03-30 Numerical simulation method of buried tube temperature field of ground source heat pump under thermoosmosis coupling

Country Status (1)

Country Link
CN (1) CN102663243B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103870706A (en) * 2014-03-28 2014-06-18 中国科学院沈阳应用生态研究所 Method and device for analyzing energy consumption and emission of ground source pump system
WO2016106949A1 (en) * 2014-12-30 2016-07-07 华中科技大学 Method for simulating temperature fields of distributed underground facilities in mountain
CN106815427A (en) * 2017-01-11 2017-06-09 北京航空航天大学 A kind of hot Calculation of Physical Properties method and apparatus of propellant
CN107762495A (en) * 2016-08-18 2018-03-06 中国石油化工股份有限公司 The optimization method of ultra-high water cut stage reservoir model longitudinal grid size
CN113468743A (en) * 2021-06-30 2021-10-01 安徽工业大学 Medium-deep buried pipe fluid temperature field analysis method considering groundwater seepage
CN117744953A (en) * 2024-02-19 2024-03-22 天津大学 Simplified analysis method for shallow soil source buried pipe group of large-scale ground source heat pump system
CN117973008A (en) * 2024-01-12 2024-05-03 天津大学 Ground source heat pump system simulation method and device based on shallow soil source g-DTM model

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008090514A (en) * 2006-09-29 2008-04-17 Nippon Steel Engineering Co Ltd Method and system for predicting performance of soil heat source heat pump system and method for calculating underground temperature
CN102243192A (en) * 2011-04-15 2011-11-16 东南大学 Multifunctional rock-soil body thermal-physical property testing device for ground source heat pump

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008090514A (en) * 2006-09-29 2008-04-17 Nippon Steel Engineering Co Ltd Method and system for predicting performance of soil heat source heat pump system and method for calculating underground temperature
CN102243192A (en) * 2011-04-15 2011-11-16 东南大学 Multifunctional rock-soil body thermal-physical property testing device for ground source heat pump

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
张行洋等: "地源热泵地埋管周围土壤温度特性研究", 《制冷与空调》 *
范蕊等: "热渗耦合作用下地下埋管换热器的传热分析", 《暖通空调》 *
袁艳平等: "地源热泵地埋管换热器传热研究(3):变热流边界条件下单U形地埋管换热器的非稳态传热特性", 《暖通空调》 *
陆观立: "地下水渗流对土壤源热泵地埋管换热特性影响研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103870706A (en) * 2014-03-28 2014-06-18 中国科学院沈阳应用生态研究所 Method and device for analyzing energy consumption and emission of ground source pump system
WO2016106949A1 (en) * 2014-12-30 2016-07-07 华中科技大学 Method for simulating temperature fields of distributed underground facilities in mountain
CN107762495A (en) * 2016-08-18 2018-03-06 中国石油化工股份有限公司 The optimization method of ultra-high water cut stage reservoir model longitudinal grid size
CN106815427A (en) * 2017-01-11 2017-06-09 北京航空航天大学 A kind of hot Calculation of Physical Properties method and apparatus of propellant
CN113468743A (en) * 2021-06-30 2021-10-01 安徽工业大学 Medium-deep buried pipe fluid temperature field analysis method considering groundwater seepage
CN113468743B (en) * 2021-06-30 2023-11-10 安徽工业大学 Medium-deep buried pipe fluid temperature field analysis method considering groundwater seepage
CN117973008A (en) * 2024-01-12 2024-05-03 天津大学 Ground source heat pump system simulation method and device based on shallow soil source g-DTM model
CN117973008B (en) * 2024-01-12 2024-08-16 天津大学 Ground source heat pump system simulation method and device based on shallow soil source g-DTM model
CN117744953A (en) * 2024-02-19 2024-03-22 天津大学 Simplified analysis method for shallow soil source buried pipe group of large-scale ground source heat pump system
CN117744953B (en) * 2024-02-19 2024-04-19 天津大学 Simplified analysis method for shallow soil source buried pipe group of large-scale ground source heat pump system

Also Published As

Publication number Publication date
CN102663243B (en) 2015-01-07

Similar Documents

Publication Publication Date Title
CN102663243B (en) Numerical simulation method of buried tube temperature field of ground source heat pump under thermoosmosis coupling
Rees et al. A three-dimensional numerical model of borehole heat exchanger heat transfer and fluid flow
Fossa The temperature penalty approach to the design of borehole heat exchangers for heat pump applications
Gashti et al. Numerical modelling of thermal regimes in steel energy pile foundations: A case study
Koohi-Fayegh et al. An analytical approach to evaluating the effect of thermal interaction of geothermal heat exchangers on ground heat pump efficiency
Zhang et al. Investigation on heat transfer around buried coils of pile foundation heat exchangers for ground-coupled heat pump applications
Zhou et al. Numerically simulating the thermal behaviors in groundwater wells of groundwater heat pump
Gao et al. Numerical simulation of the thermal interaction between pumping and injecting well groups
CN110968967B (en) Heat transfer coupling simulation order reduction method for underground pipe heat exchanger
CN112597639B (en) Heat transfer dimension reduction analysis method and system for porous buried pipe heat exchanger
Zhang et al. An improved evaluation method for thermal performance of borehole heat exchanger
Li et al. A new model and solutions for a spiral heat exchanger and its experimental validation
Wołoszyn et al. Sensitivity analysis of efficiency thermal energy storage on selected rock mass and grout parameters using design of experiment method
CN102134911B (en) Method for calculating mass concrete water pipe cooling temperature field
CN114117819B (en) Steady-state simulation method for hot steam network
CN104376188B (en) A kind of computational methods of geothermal well heat loss
CN113792462A (en) Method for rapidly acquiring heat extraction energy of intermediate-deep geothermal buried pipe
CN109063230A (en) Different gushing waters or permeable condition go down into a mine the evaluation method of the process of flooding
Li A new constant heat flux model for vertical U-tube ground heat exchangers
Belzile et al. Semi-analytical model for geothermal borefields with independent inlet conditions
CN106777731A (en) A kind of gas ductwork hydraulic calculation method
Ji et al. Development and validation of a three-dimensional numerical model for predicting the ground temperature distribution
CN104732111A (en) Efficient hourly numerical simulation method of ground source heat pump
Gao et al. A numerical study on performance efficiency of a low-temperature horizontal ground-source heat pump system
Tan et al. A Biological Immune Mechanism‐Based Quantum PSO Algorithm and Its Application in Back Analysis for Seepage Parameters

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20170512

Address after: Hengxi Street Shishan road Jiangning District Nanjing city Jiangsu province 210000 No. 49

Patentee after: Nanjing Hongtai Energy Technology Co., Ltd.

Address before: 215500 Changshou City City, Changshu Province, Southeast Economic Development Zone, Suzhou, Jin Road, No. 1, building 8, No.

Patentee before: Changshu Nanjing Normal University Development Research Academy Institute Co., Ltd.

TR01 Transfer of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150107

Termination date: 20180330

CF01 Termination of patent right due to non-payment of annual fee