Skip to content

Commit c18e817

Browse files
ogsbotogsbot
committed
Update
1 parent 4927e47 commit c18e817

20 files changed

+16762
-0
lines changed

‎Parabolic/T/3D_3BHEs_array/3bhes_1U.gml‎

Lines changed: 43 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,43 @@
1+
<?xml version="1.0" encoding="ISO-8859-1"?>
2+
<?xml-stylesheet type="text/xsl" href="OpenGeoSysGLI.xsl"?>
3+
4+
<OpenGeoSysGLI xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://www.opengeosys.org/images/xsd/OpenGeoSysGLI.xsd" xmlns:ogs="http://www.opengeosys.org">
5+
<name>3bhes_1U</name>
6+
<points>
7+
<point id="0" x="-25.0" y="0.0" z="0.0" name="TOP_left_down"/>
8+
<point id="1" x="-25.0" y="50.0" z="0.0" name="TOP_left_upper"/>
9+
<point id="2" x="25.0" y="50.0" z="0.0" name="TOP_right_upper"/>
10+
<point id="3" x="25.0" y="0.0" z="0.0" name="TOP_right_down"/>
11+
<point id="4" x="-25.0" y="0.0" z="-72.0" name="bottom_left_down"/>
12+
<point id="5" x="-25.0" y="50.0" z="-72.0" name="bottom_left_upper"/>
13+
<point id="6" x="25.0" y="50.0" z="-72.0" name="bottom_right_uppert"/>
14+
<point id="7" x="25.0" y="0.0" z="-72.0" name="bottom_right_down"/>
15+
16+
<point id="8" x="0.0" y="25.0" z="-52.0" name="BHE1_BOTTOM"/>
17+
<point id="9" x="0.0" y="25.0" z="-2.0" name="BHE1_TOP"/>
18+
<point id="10" x="-6.0" y="25.0" z="-52.0" name="BHE2_BOTTOM"/>
19+
<point id="11" x="-6.0" y="25.0" z="-2.0" name="BHE2_TOP"/>
20+
<point id="12" x="0.0" y="25.0" z="-52.0" name="BHE3_BOTTOM"/>
21+
<point id="13" x="0.0" y="25.0" z="-2.0" name="BHE3_TOP"/>
22+
</points>
23+
<polylines>
24+
<polyline id="0" name="BHE_1">
25+
<pnt>9</pnt>
26+
<pnt>8</pnt>
27+
</polyline>
28+
<polyline id="1" name="BHE_2">
29+
<pnt>11</pnt>
30+
<pnt>10</pnt>
31+
</polyline>
32+
<polyline id="2" name="BHE_3">
33+
<pnt>13</pnt>
34+
<pnt>12</pnt>
35+
</polyline>
36+
</polylines>
37+
<surfaces>
38+
<surface id="0" name="top">
39+
<element p1="0" p2="1" p3="2"/>
40+
<element p1="0" p2="3" p3="2"/>
41+
</surface>
42+
</surfaces>
43+
</OpenGeoSysGLI>

‎Parabolic/T/3D_3BHEs_array/3bhes_1U.prj‎

Lines changed: 350 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,350 @@
1+
<?xml version="1.0" encoding="ISO-8859-1"?>
2+
<OpenGeoSysProject>
3+
<mesh>3bhes_1U.vtu</mesh>
4+
<geometry>3bhes_1U.gml</geometry>
5+
<python_script>bcs_tespy.py</python_script>
6+
<processes>
7+
<process>
8+
<name>HeatTransportBHE</name>
9+
<type>HEAT_TRANSPORT_BHE</type>
10+
<integration_order>2</integration_order>
11+
<process_variables>
12+
<process_variable>temperature_soil</process_variable>
13+
<process_variable>temperature_BHE1</process_variable>
14+
<process_variable>temperature_BHE2</process_variable>
15+
<process_variable>temperature_BHE3</process_variable>
16+
</process_variables>
17+
<borehole_heat_exchangers>
18+
<borehole_heat_exchanger>
19+
<type>1U</type>
20+
<bhe_if_use_python_bc>true</bhe_if_use_python_bc>
21+
<flow_and_temperature_control>
22+
<type>TemperatureCurveConstantFlow</type>
23+
<flow_rate>2.0e-4</flow_rate>
24+
<temperature_curve>inflow_temperature</temperature_curve>
25+
</flow_and_temperature_control>
26+
<borehole>
27+
<length>50.0</length>
28+
<diameter>0.126</diameter>
29+
</borehole>
30+
<grout>
31+
<density>2190.0</density>
32+
<porosity>0.0</porosity>
33+
<heat_capacity>1735.160</heat_capacity>
34+
<thermal_conductivity>0.806</thermal_conductivity>
35+
</grout>
36+
<pipes>
37+
<inlet>
38+
<diameter> 0.013665</diameter>
39+
<wall_thickness>0.003035</wall_thickness>
40+
<wall_thermal_conductivity>0.39</wall_thermal_conductivity>
41+
</inlet>
42+
<outlet>
43+
<diameter>0.013665</diameter>
44+
<wall_thickness>0.003035</wall_thickness>
45+
<wall_thermal_conductivity>0.39</wall_thermal_conductivity>
46+
</outlet>
47+
<distance_between_pipes>0.053</distance_between_pipes>
48+
<longitudinal_dispersion_length>0.001</longitudinal_dispersion_length>
49+
</pipes>
50+
<refrigerant>
51+
<density>992.92</density>
52+
<viscosity>0.00067418</viscosity>
53+
<specific_heat_capacity>4198</specific_heat_capacity>
54+
<thermal_conductivity>0.62863</thermal_conductivity>
55+
<reference_temperature>298.15</reference_temperature>
56+
</refrigerant>
57+
</borehole_heat_exchanger>
58+
<borehole_heat_exchanger>
59+
<type>1U</type>
60+
<bhe_if_use_python_bc>true</bhe_if_use_python_bc>
61+
<flow_and_temperature_control>
62+
<type>TemperatureCurveConstantFlow</type>
63+
<flow_rate>2.0e-4</flow_rate>
64+
<temperature_curve>inflow_temperature</temperature_curve>
65+
</flow_and_temperature_control>
66+
<borehole>
67+
<length>50.0</length>
68+
<diameter>0.126</diameter>
69+
</borehole>
70+
<grout>
71+
<density>2190.0</density>
72+
<porosity>0.0</porosity>
73+
<heat_capacity>1735.160</heat_capacity>
74+
<thermal_conductivity>0.806</thermal_conductivity>
75+
</grout>
76+
<pipes>
77+
<inlet>
78+
<diameter> 0.013665</diameter>
79+
<wall_thickness>0.003035</wall_thickness>
80+
<wall_thermal_conductivity>0.39</wall_thermal_conductivity>
81+
</inlet>
82+
<outlet>
83+
<diameter>0.013665</diameter>
84+
<wall_thickness>0.003035</wall_thickness>
85+
<wall_thermal_conductivity>0.39</wall_thermal_conductivity>
86+
</outlet>
87+
<distance_between_pipes>0.053</distance_between_pipes>
88+
<longitudinal_dispersion_length>0.001</longitudinal_dispersion_length>
89+
</pipes>
90+
<refrigerant>
91+
<density>992.92</density>
92+
<viscosity>0.00067418</viscosity>
93+
<specific_heat_capacity>4198</specific_heat_capacity>
94+
<thermal_conductivity>0.62863</thermal_conductivity>
95+
<reference_temperature>298.15</reference_temperature>
96+
</refrigerant>
97+
</borehole_heat_exchanger>
98+
<borehole_heat_exchanger>
99+
<type>1U</type>
100+
<bhe_if_use_python_bc>true</bhe_if_use_python_bc>
101+
<flow_and_temperature_control>
102+
<type>TemperatureCurveConstantFlow</type>
103+
<flow_rate>2.0e-4</flow_rate>
104+
<temperature_curve>inflow_temperature</temperature_curve>
105+
</flow_and_temperature_control>
106+
<borehole>
107+
<length>50.0</length>
108+
<diameter>0.126</diameter>
109+
</borehole>
110+
<grout>
111+
<density>2190.0</density>
112+
<porosity>0.0</porosity>
113+
<heat_capacity>1735.160</heat_capacity>
114+
<thermal_conductivity>0.806</thermal_conductivity>
115+
</grout>
116+
<pipes>
117+
<inlet>
118+
<diameter> 0.013665</diameter>
119+
<wall_thickness>0.003035</wall_thickness>
120+
<wall_thermal_conductivity>0.39</wall_thermal_conductivity>
121+
</inlet>
122+
<outlet>
123+
<diameter>0.013665</diameter>
124+
<wall_thickness>0.003035</wall_thickness>
125+
<wall_thermal_conductivity>0.39</wall_thermal_conductivity>
126+
</outlet>
127+
<distance_between_pipes>0.053</distance_between_pipes>
128+
<longitudinal_dispersion_length>0.001</longitudinal_dispersion_length>
129+
</pipes>
130+
<refrigerant>
131+
<density>992.92</density>
132+
<viscosity>0.00067418</viscosity>
133+
<specific_heat_capacity>4198</specific_heat_capacity>
134+
<thermal_conductivity>0.62863</thermal_conductivity>
135+
<reference_temperature>298.15</reference_temperature>
136+
</refrigerant>
137+
</borehole_heat_exchanger>
138+
</borehole_heat_exchangers>
139+
</process>
140+
</processes>
141+
<media>
142+
<medium id="0">
143+
<phases>
144+
<phase>
145+
<type>AqueousLiquid</type>
146+
<properties>
147+
<property>
148+
<name>phase_velocity</name>
149+
<type>Constant</type>
150+
<value>0 0 0</value>
151+
</property>
152+
<property>
153+
<name>specific_heat_capacity</name>
154+
<type>Constant</type>
155+
<value>4068</value>
156+
</property>
157+
<property>
158+
<name>density</name>
159+
<type>Constant</type>
160+
<value>992.92</value>
161+
</property>
162+
</properties>
163+
</phase>
164+
<phase>
165+
<type>Solid</type>
166+
<properties>
167+
<property>
168+
<name>specific_heat_capacity</name>
169+
<type>Constant</type>
170+
<value>1780</value>
171+
</property>
172+
<property>
173+
<name>density</name>
174+
<type>Constant</type>
175+
<value>1120</value>
176+
</property>
177+
</properties>
178+
</phase>
179+
<phase>
180+
<type>Gas</type>
181+
<properties>
182+
<property>
183+
<name>specific_heat_capacity</name>
184+
<type>Constant</type>
185+
<value>1000</value>
186+
</property>
187+
<property>
188+
<name>density</name>
189+
<type>Constant</type>
190+
<value>2500</value>
191+
</property>
192+
</properties>
193+
</phase>
194+
</phases>
195+
<properties>
196+
<property>
197+
<name>porosity</name>
198+
<type>Constant</type>
199+
<value>0</value>
200+
</property>
201+
<property>
202+
<name>thermal_conductivity</name>
203+
<type>Constant</type>
204+
<value>2.4</value>
205+
</property>
206+
<property>
207+
<name>thermal_longitudinal_dispersivity</name>
208+
<type>Constant</type>
209+
<value>0</value>
210+
</property>
211+
<property>
212+
<name>thermal_transversal_dispersivity</name>
213+
<type>Constant</type>
214+
<value>0</value>
215+
</property>
216+
</properties>
217+
</medium>
218+
</media>
219+
<time_loop>
220+
<processes>
221+
<process ref="HeatTransportBHE">
222+
<nonlinear_solver>basic_picard</nonlinear_solver>
223+
<convergence_criterion>
224+
<type>DeltaX</type>
225+
<norm_type>NORM2</norm_type>
226+
<reltol>1.e-5</reltol>
227+
</convergence_criterion>
228+
<time_discretization><type>BackwardEuler</type></time_discretization>
229+
<time_stepping>
230+
<type>FixedTimeStepping</type>
231+
<t_initial> 0.0 </t_initial>
232+
<!-- use the following for full simulation
233+
<t_end> 15552000 </t_end>
234+
-->
235+
<t_end> 7200 </t_end>
236+
<timesteps>
237+
<!-- use the following for full simulation
238+
<pair><repeat>10</repeat><delta_t>720</delta_t></pair>
239+
<pair><repeat>4</repeat><delta_t>3600</delta_t></pair>
240+
<pair><repeat>1</repeat><delta_t>10800</delta_t></pair>
241+
-->
242+
<pair>
243+
<repeat>10</repeat>
244+
<delta_t>720</delta_t>
245+
</pair>
246+
</timesteps>
247+
</time_stepping>
248+
</process>
249+
</processes>
250+
<output>
251+
<type>VTK</type>
252+
<prefix>3bhes_1U</prefix>
253+
<timesteps>
254+
<!-- use the following for full simulation
255+
<pair><repeat> 1</repeat><each_steps> 92 </each_steps></pair>
256+
<pair><repeat> 1</repeat><each_steps> 80 </each_steps></pair>
257+
-->
258+
<pair>
259+
<repeat> 10</repeat>
260+
<each_steps> 1 </each_steps>
261+
</pair>
262+
</timesteps>
263+
<variables>
264+
<variable>temperature_soil</variable>
265+
<variable>temperature_BHE1</variable>
266+
<variable>temperature_BHE2</variable>
267+
<variable>temperature_BHE3</variable>
268+
</variables>
269+
</output>
270+
</time_loop>
271+
<parameters>
272+
<parameter>
273+
<name>T0</name>
274+
<type>Constant</type>
275+
<value>315.15</value>
276+
</parameter>
277+
<parameter>
278+
<name>T0_BHE</name>
279+
<type>Constant</type>
280+
<values>305.36 305.13 305.23 305.115</values>
281+
</parameter>
282+
</parameters>
283+
<process_variables>
284+
<process_variable>
285+
<name>temperature_soil</name>
286+
<components>1</components>
287+
<order>1</order>
288+
<initial_condition>T0</initial_condition>
289+
<boundary_conditions>
290+
<boundary_condition>
291+
<geometrical_set>3bhes_1U</geometrical_set>
292+
<geometry>top</geometry>
293+
<type>Dirichlet</type>
294+
<parameter>T0</parameter>
295+
</boundary_condition>
296+
</boundary_conditions>
297+
</process_variable>
298+
<process_variable>
299+
<name>temperature_BHE1</name>
300+
<components>4</components>
301+
<order>1</order>
302+
<initial_condition>T0_BHE</initial_condition>
303+
</process_variable>
304+
<process_variable>
305+
<name>temperature_BHE2</name>
306+
<components>4</components>
307+
<order>1</order>
308+
<initial_condition>T0_BHE</initial_condition>
309+
</process_variable>
310+
<process_variable>
311+
<name>temperature_BHE3</name>
312+
<components>4</components>
313+
<order>1</order>
314+
<initial_condition>T0_BHE</initial_condition>
315+
</process_variable>
316+
</process_variables>
317+
<nonlinear_solvers>
318+
<nonlinear_solver>
319+
<name>basic_picard</name>
320+
<type>Picard</type>
321+
<max_iter>100</max_iter>
322+
<linear_solver>general_linear_solver</linear_solver>
323+
</nonlinear_solver>
324+
</nonlinear_solvers>
325+
<linear_solvers>
326+
<linear_solver>
327+
<name>general_linear_solver</name>
328+
<lis>-i cg -p jacobi -tol 1e-16 -maxiter 10000</lis>
329+
<eigen>
330+
<solver_type>BiCGSTAB</solver_type>
331+
<precon_type>ILUT</precon_type>
332+
<max_iteration_step>1000</max_iteration_step>
333+
<error_tolerance>1e-16</error_tolerance>
334+
</eigen>
335+
<petsc>
336+
<prefix>gw</prefix>
337+
<parameters>-gw_ksp_type cg -gw_pc_type bjacobi -gw_ksp_rtol 1e-16 -gw_ksp_max_it 10000</parameters>
338+
</petsc>
339+
</linear_solver>
340+
</linear_solvers>
341+
<curves>
342+
<curve>
343+
<name>inflow_temperature</name>
344+
<coords>0 933120000
345+
</coords>
346+
<values>295.3611111 295.3611111
347+
</values>
348+
</curve>
349+
</curves>
350+
</OpenGeoSysProject>

0 commit comments

Comments
 (0)