Environment - simulationEnvironment: publicData = /home/hudson/saved_omc/libraries/.openmodelica/libraries/TransiEnt 2.0.3-master//Tables/ privateData = /home/hudson/saved_omc/libraries/.openmodelica/libraries/TransiEnt 2.0.3-master//Tables/ superstructureTables = /home/hudson/saved_omc/libraries/.openmodelica/libraries/TransiEnt 2.0.3-master//Tables/superstructure/Tables/ superstructureInput = /home/hudson/saved_omc/libraries/.openmodelica/libraries/TransiEnt 2.0.3-master//Tables/superstructure/RegionInformation/ startTime=0 stopTime=1 tolerance=1e-06 numberOfIntervals=2500 stepSize=0.0004 Regular simulation (gbode: the generic bi-rate ODE solver, with the radauIIA3 method): (cd run_gbode && ../TransiEnt_TransiEnt.Storage.Heat.ATES_L4.Check.Check_ValidityAnalysis -abortSlowSimulation -alarm=240 -lv LOG_STATS -inputPath=.. -s gbode -gbm=radauIIA3) LOG_STDOUT | info | ==================================================================== LOG_STDOUT | info | ATES_L4 validity analysis LOG_STDOUT | info | R_th = 24.4 m, H_a = 60 m, H_a/R_th = 2.5, m_flow_max = 5 kg/s, T_0 = 30 degC, T_inj = 60 degC LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Darcy's law (creeping flow): fulfilled LOG_STDOUT | info | Re_d(r_0) = 0.0622, non-Darcy share E(r_0) = 0.0967 % (fulfilled <= 1 %) LOG_STDOUT | info | E > 10 % for r < 0.00131 m, first control volume ends at r = 1.15 m LOG_STDOUT | info | d = 0.329 mm (estimated from k and n); Ergun (1952), Zeng and Grigg (2006) LOG_STDOUT | info | Local thermal equilibrium: fulfilled LOG_STDOUT | info | d = 0.329 mm (LTE for d <= 7 mm, Gossler et al. 2020) LOG_STDOUT | info | seepage velocity > 1.6 m/d for r < 2.91 m (1.41 % of the plume volume) LOG_STDOUT | info | Darcy velocity > 12 m/d for r < 0.097 m (LTNE > 5 % for d >= 20 mm, Lee et al. 2025) LOG_STDOUT | info | Natural regional groundwater flow: fulfilled LOG_STDOUT | info | N_0 = u*t_cycle/R_th = 0, mobility ratio M = 1.71 -> N = 0 ... 0 (fulfilled <= 0.1, violated > 1) LOG_STDOUT | info | thresholds derived from LT-ATES (Bloemendal and Hartog 2018, Tas et al. 2025), high-temperature bounds after Wheatcraft and Winterberg (1985) LOG_STDOUT | info | Buoyancy switch: WARNING (buoyancy = false) LOG_STDOUT | info | Ra = 26.9 (critical 39.5, Nield and Bejan 2017), q_0 = 0.041 m/d (significant >= 0.05 m/d, Beernink et al. 2024) LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Overall: WARNING LOG_STDOUT | info | ==================================================================== LOG_STDOUT | info | ==================================================================== LOG_STDOUT | info | ATES_L4 validity analysis LOG_STDOUT | info | R_th = 75.5 m, H_a = 21 m, H_a/R_th = 0.28, m_flow_max = 31.4 kg/s, T_0 = 30 degC, T_inj = 60 degC LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Darcy's law (creeping flow): VIOLATED LOG_STDOUT | info | Re_d(r_0) = 85, non-Darcy share E(r_0) = 58.6 % (fulfilled <= 1 %) LOG_STDOUT | info | E > 10 % for r < 1.91 m, first control volume ends at r = 1.15 m LOG_STDOUT | info | d = 25 mm (input d_50); Ergun (1952), Zeng and Grigg (2006) LOG_STDOUT | info | Local thermal equilibrium: VIOLATED LOG_STDOUT | info | d = 25 mm (LTE for d <= 7 mm, Gossler et al. 2020) LOG_STDOUT | info | seepage velocity > 1.6 m/d for r < 43.5 m (33.3 % of the plume volume) LOG_STDOUT | info | Darcy velocity > 12 m/d for r < 1.74 m (LTNE > 5 % for d >= 20 mm, Lee et al. 2025) LOG_STDOUT | info | Natural regional groundwater flow: VIOLATED LOG_STDOUT | info | N_0 = u*t_cycle/R_th = 1.32, mobility ratio M = 1.71 -> N = 1.67 ... 2.26 (fulfilled <= 0.1, violated > 1) LOG_STDOUT | info | thresholds derived from LT-ATES (Bloemendal and Hartog 2018, Tas et al. 2025), high-temperature bounds after Wheatcraft and Winterberg (1985) LOG_STDOUT | info | Buoyancy switch: VIOLATED (buoyancy = false) LOG_STDOUT | info | Ra = 1.65e+03 (critical 39.5, Nield and Bejan 2017), q_0 = 3.74 m/d (significant >= 0.05 m/d, Beernink et al. 2024) LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Overall: VIOLATED LOG_STDOUT | info | ==================================================================== LOG_STDOUT | info | ==================================================================== LOG_STDOUT | info | ATES_L4 validity analysis LOG_STDOUT | info | R_th = 104 m, H_a = 21 m, H_a/R_th = 0.2, m_flow_max = 31.4 kg/s, T_0 = 30 degC, T_inj = 60 degC LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Darcy's law (creeping flow): WARNING LOG_STDOUT | info | Re_d(r_0) = 34, non-Darcy share E(r_0) = 34.6 % (fulfilled <= 1 %) LOG_STDOUT | info | E > 10 % for r < 0.714 m, first control volume ends at r = 1.15 m -> affects only the near-well pressure drop LOG_STDOUT | info | d = 10 mm (input d_50); Ergun (1952), Zeng and Grigg (2006) LOG_STDOUT | info | Local thermal equilibrium: WARNING LOG_STDOUT | info | d = 10 mm (LTE for d <= 7 mm, Gossler et al. 2020) LOG_STDOUT | info | seepage velocity > 1.6 m/d for r < 52.3 m (25.2 % of the plume volume) LOG_STDOUT | info | Darcy velocity > 12 m/d for r < 1.74 m (LTNE > 5 % for d >= 20 mm, Lee et al. 2025) LOG_STDOUT | info | Natural regional groundwater flow: WARNING LOG_STDOUT | info | N_0 = u*t_cycle/R_th = 0.0922, mobility ratio M = 1.71 -> N = 0.116 ... 0.157 (fulfilled <= 0.1, violated > 1) LOG_STDOUT | info | thresholds derived from LT-ATES (Bloemendal and Hartog 2018, Tas et al. 2025), high-temperature bounds after Wheatcraft and Winterberg (1985) LOG_STDOUT | info | Buoyancy switch: fulfilled (buoyancy = true) LOG_STDOUT | info | Ra = 78.5 (critical 39.5, Nield and Bejan 2017), q_0 = 0.265 m/d (significant >= 0.05 m/d, Beernink et al. 2024) LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Overall: WARNING LOG_STDOUT | info | ==================================================================== LOG_STDOUT | info | ==================================================================== LOG_STDOUT | info | ATES_L4 validity analysis LOG_STDOUT | info | R_th = 27.8 m, H_a = 21 m, H_a/R_th = 0.75, m_flow_max = 2.24 kg/s, T_0 = 30 degC, T_inj = 60 degC LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Darcy's law (creeping flow): fulfilled LOG_STDOUT | info | Re_d(r_0) = 0.141, non-Darcy share E(r_0) = 0.219 % (fulfilled <= 1 %) LOG_STDOUT | info | E > 10 % for r < 0.00297 m, first control volume ends at r = 1.15 m LOG_STDOUT | info | d = 0.583 mm (estimated from k and n); Ergun (1952), Zeng and Grigg (2006) LOG_STDOUT | info | Local thermal equilibrium: fulfilled LOG_STDOUT | info | d = 0.583 mm (LTE for d <= 7 mm, Gossler et al. 2020) LOG_STDOUT | info | seepage velocity > 1.6 m/d for r < 3.72 m (1.79 % of the plume volume) LOG_STDOUT | info | Darcy velocity > 12 m/d for r < 0.124 m (LTNE > 5 % for d >= 20 mm, Lee et al. 2025) LOG_STDOUT | info | Natural regional groundwater flow: fulfilled LOG_STDOUT | info | N_0 = u*t_cycle/R_th = 0, mobility ratio M = 1.71 -> N = 0 ... 0 (fulfilled <= 0.1, violated > 1) LOG_STDOUT | info | thresholds derived from LT-ATES (Bloemendal and Hartog 2018, Tas et al. 2025), high-temperature bounds after Wheatcraft and Winterberg (1985) LOG_STDOUT | info | Buoyancy switch: fulfilled (buoyancy = true) LOG_STDOUT | info | Ra = 70.7 (critical 39.5, Nield and Bejan 2017), q_0 = 0.199 m/d (significant >= 0.05 m/d, Beernink et al. 2024) LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Overall: fulfilled LOG_STDOUT | info | ==================================================================== LOG_STDOUT | info | ==================================================================== LOG_STDOUT | info | ATES_L4 validity analysis LOG_STDOUT | info | R_th = 104 m, H_a = 21 m, H_a/R_th = 0.2, m_flow_max = 31.4 kg/s, T_0 = 30 degC, T_inj = 60 degC LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Darcy's law (creeping flow): WARNING LOG_STDOUT | info | Re_d(r_0) = 1.98, non-Darcy share E(r_0) = 2.99 % (fulfilled <= 1 %) LOG_STDOUT | info | E > 10 % for r < 0.0417 m, first control volume ends at r = 1.15 m -> affects only the near-well pressure drop LOG_STDOUT | info | d = 0.583 mm (estimated from k and n); Ergun (1952), Zeng and Grigg (2006) LOG_STDOUT | info | Local thermal equilibrium: fulfilled LOG_STDOUT | info | d = 0.583 mm (LTE for d <= 7 mm, Gossler et al. 2020) LOG_STDOUT | info | seepage velocity > 1.6 m/d for r < 52.3 m (25.2 % of the plume volume) LOG_STDOUT | info | Darcy velocity > 12 m/d for r < 1.74 m (LTNE > 5 % for d >= 20 mm, Lee et al. 2025) LOG_STDOUT | info | Natural regional groundwater flow: fulfilled LOG_STDOUT | info | N_0 = u*t_cycle/R_th = 0, mobility ratio M = 1.71 -> N = 0 ... 0 (fulfilled <= 0.1, violated > 1) LOG_STDOUT | info | thresholds derived from LT-ATES (Bloemendal and Hartog 2018, Tas et al. 2025), high-temperature bounds after Wheatcraft and Winterberg (1985) LOG_STDOUT | info | Buoyancy switch: fulfilled (buoyancy = true) LOG_STDOUT | info | Ra = 70.7 (critical 39.5, Nield and Bejan 2017), q_0 = 0.199 m/d (significant >= 0.05 m/d, Beernink et al. 2024) LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Overall: WARNING LOG_STDOUT | info | ==================================================================== LOG_STDOUT | info | ==================================================================== LOG_STDOUT | info | ATES_L4 validity analysis LOG_STDOUT | info | R_th = 25.6 m, H_a = 21 m, H_a/R_th = 0.82, m_flow_max = 19.4 kg/s, T_0 = 20 degC, T_inj = 60 degC LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Darcy's law (creeping flow): WARNING LOG_STDOUT | info | Re_d(r_0) = 1.84, non-Darcy share E(r_0) = 2.78 % (fulfilled <= 1 %) LOG_STDOUT | info | E > 10 % for r < 0.0257 m, first control volume ends at r = 1.1 m -> affects only the near-well pressure drop LOG_STDOUT | info | d = 0.583 mm (estimated from k and n); Ergun (1952), Zeng and Grigg (2006) LOG_STDOUT | info | Local thermal equilibrium: fulfilled LOG_STDOUT | info | d = 0.583 mm (LTE for d <= 7 mm, Gossler et al. 2020) LOG_STDOUT | info | seepage velocity > 1.6 m/d for r < 32.2 m (100 % of the plume volume) LOG_STDOUT | info | Darcy velocity > 12 m/d for r < 1.07 m (LTNE > 5 % for d >= 20 mm, Lee et al. 2025) LOG_STDOUT | info | Natural regional groundwater flow: WARNING LOG_STDOUT | info | N_0 = u*t_cycle/R_th = 0.0895, mobility ratio M = 2.15 -> N = 0.122 ... 0.192 (fulfilled <= 0.1, violated > 1) LOG_STDOUT | info | thresholds derived from LT-ATES (Bloemendal and Hartog 2018, Tas et al. 2025), high-temperature bounds after Wheatcraft and Winterberg (1985) LOG_STDOUT | info | Buoyancy switch: VIOLATED (buoyancy = false) LOG_STDOUT | info | Ra = 78.7 (critical 39.5, Nield and Bejan 2017), q_0 = 0.206 m/d (significant >= 0.05 m/d, Beernink et al. 2024) LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Overall: VIOLATED LOG_STDOUT | info | ==================================================================== LOG_STDOUT | info | ==================================================================== LOG_STDOUT | info | ATES_L4 validity analysis LOG_STDOUT | info | R_th = 25.6 m, H_a = 21 m, H_a/R_th = 0.82, m_flow_max = 19.4 kg/s, T_0 = 20 degC, T_inj = 60 degC LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Darcy's law (creeping flow): WARNING LOG_STDOUT | info | Re_d(r_0) = 1.84, non-Darcy share E(r_0) = 2.78 % (fulfilled <= 1 %) LOG_STDOUT | info | E > 10 % for r < 0.0257 m, first control volume ends at r = 1.1 m -> affects only the near-well pressure drop LOG_STDOUT | info | d = 0.583 mm (estimated from k and n); Ergun (1952), Zeng and Grigg (2006) LOG_STDOUT | info | Local thermal equilibrium: fulfilled LOG_STDOUT | info | d = 0.583 mm (LTE for d <= 7 mm, Gossler et al. 2020) LOG_STDOUT | info | seepage velocity > 1.6 m/d for r < 32.2 m (100 % of the plume volume) LOG_STDOUT | info | Darcy velocity > 12 m/d for r < 1.07 m (LTNE > 5 % for d >= 20 mm, Lee et al. 2025) LOG_STDOUT | info | Natural regional groundwater flow: WARNING LOG_STDOUT | info | N_0 = u*t_cycle/R_th = 0.0895, mobility ratio M = 2.15 -> N = 0.122 ... 0.192 (fulfilled <= 0.1, violated > 1) LOG_STDOUT | info | thresholds derived from LT-ATES (Bloemendal and Hartog 2018, Tas et al. 2025), high-temperature bounds after Wheatcraft and Winterberg (1985) LOG_STDOUT | info | Buoyancy switch: fulfilled (buoyancy = true) LOG_STDOUT | info | Ra = 78.7 (critical 39.5, Nield and Bejan 2017), q_0 = 0.206 m/d (significant >= 0.05 m/d, Beernink et al. 2024) LOG_STDOUT | info | -------------------------------------------------------------------- LOG_STDOUT | info | Overall: WARNING LOG_STDOUT | info | ==================================================================== LOG_SUCCESS | info | The initialization finished successfully without homotopy method. LOG_STATS | info | ### STATISTICS ### | | | | | timer | | | | | | 0.00260173s reading init.xml | | | | | | 0.000189356s reading info.xml | | | | | | 3.5296e-05s [ 0.1%] pre-initialization | | | | | | 0.024897s [ 91.0%] initialization | | | | | | 4.21e-07s [ 0.0%] steps | | | | | | 2.6053e-05s [ 0.1%] solver (excl. callbacks) | | | | | | 0s [ 0.0%] creating output-file | | | | | | 1.0722e-05s [ 0.0%] event-handling | | | | | | 0.000893804s [ 3.3%] overhead | | | | | | 0.00150191s [ 5.5%] simulation | | | | | | 0.0273652s [100.0%] total | | | | | events | | | | | | 0 state events | | | | | | 0 time events | | | | | solver: euler | | | | | | 2500 steps taken | | | | | | 2500 calls of functionODE | | | | | | 0 evaluations of jacobian | | | | | | 0 error test failures | | | | | | 0 convergence test failures | | | | | | 0s time of jacobian evaluation LOG_SUCCESS | info | The simulation finished successfully.