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:
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:
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
(1)
Utilization formula (2) is carried out the longitudinal grid implicit iterative, calculates that (τ+r) is the temperature field data constantly
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:
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.
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:
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:
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.