From 1c623f1f221d546be9a6e1b627a6892c2e77a29d Mon Sep 17 00:00:00 2001 From: kfbeckers Date: Wed, 24 Jul 2024 18:14:00 -0600 Subject: [PATCH] changes invertlaplace method to Stehfest Algorithm to make Gringarten model more stable --- src/geophires_x/MPFReservoir.py | 2 +- tests/examples/Fervo_Norbeck_Latimer_2023.out | 10 +++++----- 2 files changed, 6 insertions(+), 6 deletions(-) diff --git a/src/geophires_x/MPFReservoir.py b/src/geophires_x/MPFReservoir.py index 72fdde64..de466381 100644 --- a/src/geophires_x/MPFReservoir.py +++ b/src/geophires_x/MPFReservoir.py @@ -88,7 +88,7 @@ def Calculate(self, model: Model): Twnd = [] try: for t in range(1, len(model.reserv.timevector.value)): - Twnd = Twnd + [float(invertlaplace(fp, td[t], method='talbot'))] + Twnd = Twnd + [float(invertlaplace(fp, td[t], method='stehfest'))] except: print( "Error: GEOPHIRES could not execute numerical inverse laplace calculation for reservoir model 1. Simulation will abort.") diff --git a/tests/examples/Fervo_Norbeck_Latimer_2023.out b/tests/examples/Fervo_Norbeck_Latimer_2023.out index 53bb4b03..0341a33a 100644 --- a/tests/examples/Fervo_Norbeck_Latimer_2023.out +++ b/tests/examples/Fervo_Norbeck_Latimer_2023.out @@ -4,10 +4,10 @@ Simulation Metadata ---------------------- - GEOPHIRES Version: 3.4.42 - Simulation Date: 2024-07-17 - Simulation Time: 11:19 - Calculation Time: 0.353 sec + GEOPHIRES Version: 3.5.0 + Simulation Date: 2024-07-24 + Simulation Time: 18:06 + Calculation Time: 1.152 sec ***SUMMARY OF RESULTS*** @@ -144,7 +144,7 @@ Simulation Metadata 6 1.0129 184.12 0.6919 2.1612 8.8239 7 1.0129 184.13 0.6920 2.1612 8.8237 8 1.0108 183.75 0.6925 2.1430 8.7720 - 9 1.0019 182.13 0.6948 2.0640 8.5463 + 9 1.0019 182.13 0.6948 2.0640 8.5464 10 0.9771 177.63 0.7009 1.8532 7.9258