(cd run_fmpy-fmi3-cs-cvode && python3 -m fmpy simulate --output-file ../TransiEnt_TransiEnt.Storage.Heat.ATES_L4.Check.Check_ValidityAnalysis_fmpy-fmi3-cs-cvode_res.csv --start-time 0 --stop-time 1 --timeout 240 --relative-tolerance 1e-06 --interface-type CoSimulation --output-interval 0.0004 ../TransiEnt_TransiEnt_Storage_Heat_ATES_L4_Check_Check_ValidityAnalysis_fmpy-fmi3-cs-cvode.fmu) 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 | ====================================================================