From 4e01911cc40d6ab450d7cf7a264c4bd2f79c8238 Mon Sep 17 00:00:00 2001 From: Ranpeng Li Date: Wed, 5 Feb 2025 14:37:21 +0100 Subject: [PATCH] Add an option to make it crash when look-up table does not cover the range, and manually set it for entropy method --- doc/modules/changes/20250110_ranpengli | 4 + include/aspect/structured_data.h | 6 +- source/material_model/utilities.cc | 2 +- source/structured_data.cc | 31 +++- tests/entropy_initial_lookup.prm | 2 +- tests/entropy_initial_lookup/screen-output | 36 ++-- tests/entropy_initial_lookup/statistics | 32 ++++ tests/entropy_initial_lookup_wb.prm | 2 +- tests/entropy_initial_lookup_wb/screen-output | 18 +- tests/entropy_initial_lookup_wb/statistics | 35 ++++ tests/entropy_plasticity.prm | 2 +- tests/entropy_plasticity/screen-output | 166 ++++-------------- tests/entropy_plasticity/statistics | 6 +- 13 files changed, 174 insertions(+), 168 deletions(-) create mode 100644 doc/modules/changes/20250110_ranpengli create mode 100644 tests/entropy_initial_lookup/statistics create mode 100644 tests/entropy_initial_lookup_wb/statistics diff --git a/doc/modules/changes/20250110_ranpengli b/doc/modules/changes/20250110_ranpengli new file mode 100644 index 00000000000..be9d1ce96c0 --- /dev/null +++ b/doc/modules/changes/20250110_ranpengli @@ -0,0 +1,4 @@ +Changed: The entropy reader throws when the range of the provided look-up table +does not fully cover the entropy-pressure range in the model. +
+(Ranpeng Li, 2025/01/10) diff --git a/include/aspect/structured_data.h b/include/aspect/structured_data.h index af5b56be207..4593c1c7943 100644 --- a/include/aspect/structured_data.h +++ b/include/aspect/structured_data.h @@ -191,10 +191,14 @@ namespace aspect * @param component The index (starting at 0) of the data column to be * returned. The index is therefore less than the number of data * columns in the data file (or specified in the constructor). + * @param crash_if_not_in_range If set to true, the function will throw + * when the requested position is outside the range of the coordinates + * provided by the data file. */ double get_data(const Point &position, - const unsigned int component) const; + const unsigned int component, + const bool crash_if_not_in_range = false) const; /** * Returns the gradient of the function based on the bilinear diff --git a/source/material_model/utilities.cc b/source/material_model/utilities.cc index d106de25a77..d74b59f2134 100644 --- a/source/material_model/utilities.cc +++ b/source/material_model/utilities.cc @@ -837,7 +837,7 @@ namespace aspect EntropyReader::temperature(const double entropy, const double pressure) const { - const double temperature = material_lookup->get_data({entropy,pressure}, 0); + const double temperature = material_lookup->get_data({entropy,pressure}, 0, true); return temperature; } diff --git a/source/structured_data.cc b/source/structured_data.cc index 0caefa1d25a..a436ebbf024 100644 --- a/source/structured_data.cc +++ b/source/structured_data.cc @@ -862,9 +862,38 @@ namespace aspect template double StructuredDataLookup::get_data(const Point &position, - const unsigned int component) const + const unsigned int component, + const bool crash_if_not_in_range) const { Assert(component &x_coordinates = get_interpolation_point_coordinates(0); + + AssertThrow (position[0] >= x_coordinates[0] && position[0] <= x_coordinates[x_coordinates.size()-1], + ExcMessage("The requested position " + + std::to_string(position[0]) + + " is outside the range of the data (minimum value = " + + std::to_string(x_coordinates[0]) + + " , maximum value = " + + std::to_string(x_coordinates[x_coordinates.size()-1]) + + ")." + )); + + const std::vector &y_coordinates = get_interpolation_point_coordinates(1); + + AssertThrow (position[1] >= y_coordinates[0] && position[1] <= y_coordinates[y_coordinates.size()-1], + ExcMessage("The requested position " + + std::to_string(position[1]) + + " is outside the range of the data (minimum value = " + + std::to_string(y_coordinates[0]) + + " , maximum value = " + + std::to_string(y_coordinates[y_coordinates.size()-1]) + + ")." + )); + } + return data[component]->value(position); } diff --git a/tests/entropy_initial_lookup.prm b/tests/entropy_initial_lookup.prm index 34f1f0b032d..2df40382e94 100644 --- a/tests/entropy_initial_lookup.prm +++ b/tests/entropy_initial_lookup.prm @@ -5,7 +5,7 @@ set Dimension = 2 set Use years in output instead of seconds = true set End time = 0 set Nonlinear solver scheme = iterated Advection and Stokes -set Surface pressure = 0 +set Surface pressure = 25000 set Adiabatic surface temperature = 1600.0 subsection Formulation diff --git a/tests/entropy_initial_lookup/screen-output b/tests/entropy_initial_lookup/screen-output index f411dec44ee..903e218ae79 100644 --- a/tests/entropy_initial_lookup/screen-output +++ b/tests/entropy_initial_lookup/screen-output @@ -10,56 +10,56 @@ Number of degrees of freedom: 3,217 (1,206+202+603+603+603) Solving entropy system ... 3 iterations. Copying properties into prescribed compositional field density_field... done. Solving Stokes system (GMG)... 1000+3 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000495634, 4.77467e-10, 0, 5.9652 - Relative nonlinear residual (total system) after nonlinear iteration 1: 5.9652 + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000495626, 4.77402e-10, 0, 5.96521 + Relative nonlinear residual (total system) after nonlinear iteration 1: 5.96521 Copying properties into prescribed temperature field... done. Solving temperature system... 3 iterations. Solving entropy system ... 3 iterations. Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 600+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 8.51534e-10, 5.64728e-10, 0, 0.0685014 - Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0685014 + Solving Stokes system (GMG)... 676+0 iterations. + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 8.51395e-10, 5.64621e-10, 0, 0.0685012 + Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0685012 Copying properties into prescribed temperature field... done. Solving temperature system... 3 iterations. Solving entropy system ... 3 iterations. Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 500+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.01722e-09, 4.2244e-10, 0, 0.0168383 - Relative nonlinear residual (total system) after nonlinear iteration 3: 0.0168383 + Solving Stokes system (GMG)... 550+0 iterations. + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.01699e-09, 4.22363e-10, 0, 0.0168382 + Relative nonlinear residual (total system) after nonlinear iteration 3: 0.0168382 Copying properties into prescribed temperature field... done. Solving temperature system... 3 iterations. Solving entropy system ... 3 iterations. Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 286+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 7.51013e-10, 3.95146e-10, 0, 0.00246445 + Solving Stokes system (GMG)... 245+0 iterations. + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 7.509e-10, 3.95043e-10, 0, 0.00246445 Relative nonlinear residual (total system) after nonlinear iteration 4: 0.00246445 Copying properties into prescribed temperature field... done. Solving temperature system... 3 iterations. Solving entropy system ... 2 iterations. Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 194+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 7.21895e-10, 4.19474e-10, 0, 0.000253486 - Relative nonlinear residual (total system) after nonlinear iteration 5: 0.000253486 + Solving Stokes system (GMG)... 190+0 iterations. + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 7.21723e-10, 4.19413e-10, 0, 0.00025346 + Relative nonlinear residual (total system) after nonlinear iteration 5: 0.00025346 Copying properties into prescribed temperature field... done. Solving temperature system... 3 iterations. Solving entropy system ... 3 iterations. Copying properties into prescribed compositional field density_field... done. Solving Stokes system (GMG)... 76+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 7.79118e-10, 3.18101e-10, 0, 2.03324e-05 - Relative nonlinear residual (total system) after nonlinear iteration 6: 2.03324e-05 + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 7.79017e-10, 3.17962e-10, 0, 2.03204e-05 + Relative nonlinear residual (total system) after nonlinear iteration 6: 2.03204e-05 Copying properties into prescribed temperature field... done. Solving temperature system... 3 iterations. Solving entropy system ... 3 iterations. Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 55+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 5.78613e-10, 3.56974e-10, 0, 1.35817e-06 - Relative nonlinear residual (total system) after nonlinear iteration 7: 1.35817e-06 + Solving Stokes system (GMG)... 56+0 iterations. + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 5.78371e-10, 3.56725e-10, 0, 1.35606e-06 + Relative nonlinear residual (total system) after nonlinear iteration 7: 1.35606e-06 Postprocessing: diff --git a/tests/entropy_initial_lookup/statistics b/tests/entropy_initial_lookup/statistics new file mode 100644 index 00000000000..4714d1c15bc --- /dev/null +++ b/tests/entropy_initial_lookup/statistics @@ -0,0 +1,32 @@ +# 1: Time step number +# 2: Time (years) +# 3: Time step size (years) +# 4: Number of mesh cells +# 5: Number of Stokes degrees of freedom +# 6: Number of temperature degrees of freedom +# 7: Number of degrees of freedom for all compositions +# 8: Number of nonlinear iterations +# 9: Iterations for temperature solver +# 10: Iterations for composition solver 1 +# 11: Iterations for composition solver 2 +# 12: Iterations for Stokes solver +# 13: Velocity iterations in Stokes preconditioner +# 14: Schur complement iterations in Stokes preconditioner +# 15: Visualization file name +# 16: RMS velocity (m/year) +# 17: Max. velocity (m/year) +# 18: Minimal temperature (K) +# 19: Average temperature (K) +# 20: Maximal temperature (K) +# 21: Average nondimensional temperature (K) +# 22: Outward mass flux through boundary with indicator 0 ("left") (kg/yr) +# 23: Outward mass flux through boundary with indicator 1 ("right") (kg/yr) +# 24: Outward mass flux through boundary with indicator 2 ("bottom") (kg/yr) +# 25: Outward mass flux through boundary with indicator 3 ("top") (kg/yr) +# 26: Minimal value for composition entropy +# 27: Maximal value for composition entropy +# 28: Global mass for composition entropy +# 29: Minimal value for composition density_field +# 30: Maximal value for composition density_field +# 31: Global mass for composition density_field +0 0.000000000000e+00 0.000000000000e+00 100 1408 603 1206 7 25 20 4294967289 2790 3108 2827 output-entropy_initial_lookup/solution/solution-00000 7.88044505e-03 9.98412057e-03 1.59574424e+03 1.60014975e+03 1.60309788e+03 1.16764184e-03 0.00000000e+00 0.00000000e+00 4.08138251e+05 -4.08137893e+05 2.19666723e+03 2.53347698e+03 3.95671887e+13 3.13952225e+03 4.77911946e+03 6.90810628e+13 diff --git a/tests/entropy_initial_lookup_wb.prm b/tests/entropy_initial_lookup_wb.prm index 65333180d21..6290aca29b1 100644 --- a/tests/entropy_initial_lookup_wb.prm +++ b/tests/entropy_initial_lookup_wb.prm @@ -8,7 +8,7 @@ set End time = 0 set Output directory = output_3 set Timing output frequency = 10 set Pressure normalization = surface -set Surface pressure = 0 +set Surface pressure = 25000 set Adiabatic surface temperature = 1573.0 set Resume computation = false set Nonlinear solver scheme = iterated Advection and Stokes diff --git a/tests/entropy_initial_lookup_wb/screen-output b/tests/entropy_initial_lookup_wb/screen-output index 762ab6119a2..25384dfc0ad 100644 --- a/tests/entropy_initial_lookup_wb/screen-output +++ b/tests/entropy_initial_lookup_wb/screen-output @@ -39,36 +39,36 @@ Number of degrees of freedom: 47,413 (18,018+2,368+9,009+9,009+9,009) Solving entropy system ... 0 iterations. Copying properties into prescribed compositional field density_field... done. Solving Stokes system (GMG)... 33+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.60794e-16, 1.69272e-16, 0, 0.00433026 - Relative nonlinear residual (total system) after nonlinear iteration 1: 0.00433026 + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.5721e-16, 1.61845e-16, 0, 0.00433027 + Relative nonlinear residual (total system) after nonlinear iteration 1: 0.00433027 Solving temperature system... 0 iterations. Solving entropy system ... 0 iterations. Copying properties into prescribed compositional field density_field... done. Solving Stokes system (GMG)... 28+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.60794e-16, 1.69272e-16, 0, 0.00148712 - Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00148712 + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.5721e-16, 1.61845e-16, 0, 0.00148713 + Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00148713 Solving temperature system... 0 iterations. Solving entropy system ... 0 iterations. Copying properties into prescribed compositional field density_field... done. Solving Stokes system (GMG)... 23+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.60794e-16, 1.69272e-16, 0, 6.99073e-05 - Relative nonlinear residual (total system) after nonlinear iteration 3: 6.99073e-05 + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.5721e-16, 1.61845e-16, 0, 6.99077e-05 + Relative nonlinear residual (total system) after nonlinear iteration 3: 6.99077e-05 Solving temperature system... 0 iterations. Solving entropy system ... 0 iterations. Copying properties into prescribed compositional field density_field... done. Solving Stokes system (GMG)... 18+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.60794e-16, 1.69272e-16, 0, 5.0992e-06 - Relative nonlinear residual (total system) after nonlinear iteration 4: 5.0992e-06 + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.5721e-16, 1.61845e-16, 0, 5.09921e-06 + Relative nonlinear residual (total system) after nonlinear iteration 4: 5.09921e-06 Postprocessing: Writing graphical output: output-entropy_initial_lookup_wb/solution/solution-00003 RMS, max velocity: 0.0138 m/year, 0.124 m/year Temperature min/avg/max: 273 K, 2184 K, 3500 K - Mass fluxes through boundary parts: 0.004614 kg/yr, 0.06578 kg/yr, 0 kg/yr, 2.212e-08 kg/yr + Mass fluxes through boundary parts: 0.004614 kg/yr, 0.06577 kg/yr, 0 kg/yr, 2.203e-08 kg/yr Compositions min/max/mass: 577.9/2758/1.089e+17 // 3228/5592/1.939e+17 Termination requested by criterion: end time diff --git a/tests/entropy_initial_lookup_wb/statistics b/tests/entropy_initial_lookup_wb/statistics new file mode 100644 index 00000000000..1af85027c6b --- /dev/null +++ b/tests/entropy_initial_lookup_wb/statistics @@ -0,0 +1,35 @@ +# 1: Time step number +# 2: Time (years) +# 3: Time step size (years) +# 4: Number of mesh cells +# 5: Number of Stokes degrees of freedom +# 6: Number of temperature degrees of freedom +# 7: Number of degrees of freedom for all compositions +# 8: Number of nonlinear iterations +# 9: Visualization file name +# 10: RMS velocity (m/year) +# 11: Max. velocity (m/year) +# 12: Minimal temperature (K) +# 13: Average temperature (K) +# 14: Maximal temperature (K) +# 15: Average nondimensional temperature (K) +# 16: Outward mass flux through boundary with indicator 0 ("bottom") (kg/yr) +# 17: Outward mass flux through boundary with indicator 1 ("top") (kg/yr) +# 18: Outward mass flux through boundary with indicator 2 ("west") (kg/yr) +# 19: Outward mass flux through boundary with indicator 3 ("east") (kg/yr) +# 20: Minimal value for composition entropy +# 21: Maximal value for composition entropy +# 22: Global mass for composition entropy +# 23: Minimal value for composition density_field +# 24: Maximal value for composition density_field +# 25: Global mass for composition density_field +# 26: Iterations for temperature solver +# 27: Iterations for composition solver 1 +# 28: Iterations for composition solver 2 +# 29: Iterations for Stokes solver +# 30: Velocity iterations in Stokes preconditioner +# 31: Schur complement iterations in Stokes preconditioner +0 0.000000000000e+00 0.000000000000e+00 384 3739 1649 3298 0 output-entropy_initial_lookup_wb/solution/solution-00000 0.00000000e+00 0.00000000e+00 2.73000000e+02 2.17080114e+03 3.50000000e+03 5.88100755e-01 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 5.77925312e+02 2.75762818e+03 1.07732840e+17 0.00000000e+00 0.00000000e+00 0.00000000e+00 0 0 0 0 0 0 +0 0.000000000000e+00 0.000000000000e+00 684 7217 3189 6378 0 output-entropy_initial_lookup_wb/solution/solution-00001 0.00000000e+00 0.00000000e+00 2.73000000e+02 2.18276710e+03 3.50000000e+03 5.91808832e-01 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 5.77925312e+02 2.75762818e+03 1.08621269e+17 0.00000000e+00 0.00000000e+00 0.00000000e+00 0 0 0 0 0 0 +0 0.000000000000e+00 0.000000000000e+00 1023 10982 4851 9702 0 output-entropy_initial_lookup_wb/solution/solution-00002 0.00000000e+00 0.00000000e+00 2.73000000e+02 2.18453715e+03 3.50000000e+03 5.92357344e-01 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 5.77925312e+02 2.75762818e+03 1.08867722e+17 0.00000000e+00 0.00000000e+00 0.00000000e+00 0 0 0 0 0 0 +0 0.000000000000e+00 0.000000000000e+00 1950 20386 9009 18018 4 output-entropy_initial_lookup_wb/solution/solution-00003 1.37687774e-02 1.23769331e-01 2.73000000e+02 2.18434812e+03 3.50000000e+03 5.92298767e-01 4.61392016e-03 6.57718010e-02 0.00000000e+00 2.21121666e-08 5.77925312e+02 2.75762818e+03 1.08896595e+17 3.22822228e+03 5.59177718e+03 1.93948135e+17 0 0 4294967292 98 106 106 diff --git a/tests/entropy_plasticity.prm b/tests/entropy_plasticity.prm index f699984b4bb..4562225e0df 100644 --- a/tests/entropy_plasticity.prm +++ b/tests/entropy_plasticity.prm @@ -108,7 +108,7 @@ subsection Initial composition model set Coordinate system = spherical set Variable names = r,phi set Function expression = 2535.08 + 10 * sin(2*phi)*sin(pi*(r-6371000)/2890000) - 1878.618 * erfc((6371000-r)/(2*sqrt(1.152e-6*5e7*31557600))) + \ - 486.368 * max(erfc((r-3481000)/(2*sqrt(1.152e-6*5e8*31557600))), exp(-((phi - 3/4*pi)^2)/(2*(1/35*pi)^2) - ((r-3481000)^2)/(2*400000^2)) );0 + 100 * max(erfc((r-3481000)/(2*sqrt(1.152e-6*5e8*31557600))), exp(-((phi - 3/4*pi)^2)/(2*(1/35*pi)^2) - ((r-3481000)^2)/(2*400000^2)) );0 #set Function expression = 2687.748 + 10 * sin(2*phi)*sin(pi*(r-6371000)/2890000) - 2031.286 * erfc((6371000-r)/(2*sqrt(1.152e-6*5e7*31557600))) + 333.7 * erfc((r-3481000)/(2*sqrt(1.152e-6*5e8*31557600)));0 end diff --git a/tests/entropy_plasticity/screen-output b/tests/entropy_plasticity/screen-output index f2efd0b9f34..746cd724a5f 100644 --- a/tests/entropy_plasticity/screen-output +++ b/tests/entropy_plasticity/screen-output @@ -9,177 +9,81 @@ Number of degrees of freedom: 17,184 (6,528+864+3,264+3,264+3,264) Solving temperature system... 0 iterations. Solving entropy system ... 0 iterations. Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 97+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.63904e-16, 1.61442e-16, 0, 0.00367688 - Relative nonlinear residual (total system) after nonlinear iteration 1: 0.00367688 + Solving Stokes system (GMG)... 75+0 iterations. + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.61151e-16, 1.59986e-16, 0, 0.00372447 + Relative nonlinear residual (total system) after nonlinear iteration 1: 0.00372447 Copying properties into prescribed temperature field... done. Solving temperature system... 0 iterations. Solving entropy system ... 0 iterations. Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 72+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.63904e-16, 1.61442e-16, 0, 0.00195243 - Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00195243 + Solving Stokes system (GMG)... 52+0 iterations. + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.61151e-16, 1.59986e-16, 0, 0.000935405 + Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000935405 Copying properties into prescribed temperature field... done. Solving temperature system... 0 iterations. Solving entropy system ... 0 iterations. Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 64+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.63904e-16, 1.61442e-16, 0, 8.95835e-05 - Relative nonlinear residual (total system) after nonlinear iteration 3: 8.95835e-05 + Solving Stokes system (GMG)... 44+0 iterations. + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.61151e-16, 1.59986e-16, 0, 5.67375e-05 + Relative nonlinear residual (total system) after nonlinear iteration 3: 5.67375e-05 Copying properties into prescribed temperature field... done. Solving temperature system... 0 iterations. Solving entropy system ... 0 iterations. Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 49+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.63904e-16, 1.61442e-16, 0, 5.0626e-06 - Relative nonlinear residual (total system) after nonlinear iteration 4: 5.0626e-06 + Solving Stokes system (GMG)... 32+0 iterations. + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.61151e-16, 1.59986e-16, 0, 2.72951e-06 + Relative nonlinear residual (total system) after nonlinear iteration 4: 2.72951e-06 Postprocessing: Writing graphical output: output-entropy_plasticity/solution/solution-00000 - RMS, max velocity: 0.039 m/year, 0.467 m/year - Temperature min/avg/max: 321 K, 2008 K, 3448 K - Mass fluxes through boundary parts: -1.545e-06 kg/yr, 9.844e-07 kg/yr - Compositions min/max/mass: 600/2976/2.237e+17 // 3514/5130/3.852e+17 - -*** Timestep 1: t=199658 years, dt=199658 years - Copying properties into prescribed temperature field... done. - Solving temperature system... 7 iterations. - Solving entropy system ... 11 iterations. - Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 78+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.61476e-05, 0.00115486, 0, 0.000363459 - Relative nonlinear residual (total system) after nonlinear iteration 1: 0.00115486 - - Copying properties into prescribed temperature field... done. - Solving temperature system... 9 iterations. - Solving entropy system ... 11 iterations. - Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 76+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00224359, 0.000328539, 0, 0.000282533 - Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00224359 - - Copying properties into prescribed temperature field... done. - Solving temperature system... 8 iterations. - Solving entropy system ... 11 iterations. - Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 70+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000631761, 9.10519e-05, 0, 7.59395e-05 - Relative nonlinear residual (total system) after nonlinear iteration 3: 0.000631761 + RMS, max velocity: 0.0117 m/year, 0.0483 m/year + Temperature min/avg/max: 321 K, 1974 K, 2582 K + Mass fluxes through boundary parts: -6.708e-07 kg/yr, 2.047e-07 kg/yr + Compositions min/max/mass: 600/2950/2.228e+17 // 3514/5130/3.853e+17 +*** Timestep 1: t=400000 years, dt=400000 years Copying properties into prescribed temperature field... done. Solving temperature system... 8 iterations. - Solving entropy system ... 10 iterations. - Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 62+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000176959, 2.28632e-05, 0, 1.79573e-05 - Relative nonlinear residual (total system) after nonlinear iteration 4: 0.000176959 - - Copying properties into prescribed temperature field... done. - Solving temperature system... 7 iterations. - Solving entropy system ... 9 iterations. - Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 57+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 4.39517e-05, 5.25597e-06, 0, 4.016e-06 - Relative nonlinear residual (total system) after nonlinear iteration 5: 4.39517e-05 - - Copying properties into prescribed temperature field... done. - Solving temperature system... 7 iterations. Solving entropy system ... 9 iterations. Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 50+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 9.96903e-06, 1.1388e-06, 0, 8.54258e-07 - Relative nonlinear residual (total system) after nonlinear iteration 6: 9.96903e-06 - - - Postprocessing: - RMS, max velocity: 0.0442 m/year, 0.57 m/year - Temperature min/avg/max: 321 K, 2007 K, 3448 K - Mass fluxes through boundary parts: -1.65e-06 kg/yr, 1.007e-06 kg/yr - Compositions min/max/mass: 600/2986/2.237e+17 // 3514/5130/3.852e+17 - -*** Timestep 2: t=384366 years, dt=184708 years - Copying properties into prescribed temperature field... done. - Solving temperature system... 7 iterations. - Solving entropy system ... 11 iterations. - Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 76+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.64282e-05, 0.00048469, 0, 0.000319382 - Relative nonlinear residual (total system) after nonlinear iteration 1: 0.00048469 + Solving Stokes system (GMG)... 56+0 iterations. + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 4.27436e-05, 0.000983831, 0, 0.00038342 + Relative nonlinear residual (total system) after nonlinear iteration 1: 0.000983831 Copying properties into prescribed temperature field... done. Solving temperature system... 9 iterations. - Solving entropy system ... 10 iterations. - Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 67+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00317547, 9.18909e-05, 0, 8.97411e-05 - Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00317547 - - Copying properties into prescribed temperature field... done. - Solving temperature system... 8 iterations. - Solving entropy system ... 9 iterations. - Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 60+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00017361, 1.19549e-05, 0, 1.73366e-05 - Relative nonlinear residual (total system) after nonlinear iteration 3: 0.00017361 - - Copying properties into prescribed temperature field... done. - Solving temperature system... 7 iterations. Solving entropy system ... 8 iterations. Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 49+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.38815e-05, 9.61741e-07, 0, 1.16635e-06 - Relative nonlinear residual (total system) after nonlinear iteration 4: 2.38815e-05 - - Copying properties into prescribed temperature field... done. - Solving temperature system... 6 iterations. - Solving entropy system ... 6 iterations. - Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 34+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.89356e-06, 5.36067e-08, 0, 5.50108e-08 - Relative nonlinear residual (total system) after nonlinear iteration 5: 1.89356e-06 - - - Postprocessing: - RMS, max velocity: 0.0502 m/year, 0.818 m/year - Temperature min/avg/max: 321 K, 2006 K, 3448 K - Mass fluxes through boundary parts: -1.668e-06 kg/yr, 1.067e-06 kg/yr - Compositions min/max/mass: 600/2986/2.236e+17 // 3513/5130/3.852e+17 - -*** Timestep 3: t=400000 years, dt=15634.2 years - Copying properties into prescribed temperature field... done. - Solving temperature system... 7 iterations. - Solving entropy system ... 7 iterations. - Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 63+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.46452e-05, 2.5067e-05, 0, 7.81458e-06 - Relative nonlinear residual (total system) after nonlinear iteration 1: 2.5067e-05 + Solving Stokes system (GMG)... 47+0 iterations. + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00196115, 3.87599e-05, 0, 3.99673e-05 + Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00196115 Copying properties into prescribed temperature field... done. Solving temperature system... 8 iterations. - Solving entropy system ... 7 iterations. + Solving entropy system ... 6 iterations. Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 54+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00030441, 3.82669e-06, 0, 7.23454e-06 - Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00030441 + Solving Stokes system (GMG)... 39+0 iterations. + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 7.9364e-05, 2.26612e-06, 0, 4.94645e-06 + Relative nonlinear residual (total system) after nonlinear iteration 3: 7.9364e-05 Copying properties into prescribed temperature field... done. Solving temperature system... 7 iterations. - Solving entropy system ... 6 iterations. + Solving entropy system ... 5 iterations. Copying properties into prescribed compositional field density_field... done. - Solving Stokes system (GMG)... 42+0 iterations. - Relative nonlinear residuals (temperature, compositional fields, Stokes system): 7.71754e-06, 3.59224e-07, 0, 9.10599e-07 - Relative nonlinear residual (total system) after nonlinear iteration 3: 7.71754e-06 + Solving Stokes system (GMG)... 30+0 iterations. + Relative nonlinear residuals (temperature, compositional fields, Stokes system): 4.49419e-06, 1.64894e-07, 0, 3.88155e-07 + Relative nonlinear residual (total system) after nonlinear iteration 4: 4.49419e-06 Postprocessing: - RMS, max velocity: 0.0507 m/year, 0.841 m/year - Temperature min/avg/max: 321 K, 2006 K, 3448 K - Mass fluxes through boundary parts: -1.625e-06 kg/yr, 1.076e-06 kg/yr - Compositions min/max/mass: 600/2989/2.236e+17 // 3513/5130/3.852e+17 + RMS, max velocity: 0.0115 m/year, 0.0485 m/year + Temperature min/avg/max: 321 K, 1973 K, 2582 K + Mass fluxes through boundary parts: -6.514e-07 kg/yr, 1.937e-07 kg/yr + Compositions min/max/mass: 600/2950/2.227e+17 // 3516/5130/3.853e+17 Termination requested by criterion: end time diff --git a/tests/entropy_plasticity/statistics b/tests/entropy_plasticity/statistics index 34aa15253e3..7879856e5f3 100644 --- a/tests/entropy_plasticity/statistics +++ b/tests/entropy_plasticity/statistics @@ -27,7 +27,5 @@ # 27: Minimal value for composition density_field # 28: Maximal value for composition density_field # 29: Global mass for composition density_field -0 0.000000000000e+00 0.000000000000e+00 768 7392 3264 6528 4 0 0 4294967292 278 286 286 output-entropy_plasticity/solution/solution-00000 3.89724518e-02 4.67399599e-01 3.21018877e+02 2.00757498e+03 3.44750753e+03 4.47006654e-01 -1.54485848e-06 9.84396700e-07 6.00000000e+02 2.97625563e+03 2.23741621e+17 3.51431122e+03 5.13045059e+03 3.85216157e+17 -1 1.996576577258e+05 1.996576577258e+05 768 7392 3264 6528 6 46 61 4294967290 387 399 399 "" 4.41762729e-02 5.69736081e-01 3.21018877e+02 2.00684655e+03 3.44750753e+03 4.46813590e-01 -1.64951380e-06 1.00659084e-06 6.00000000e+02 2.98563135e+03 2.23690620e+17 3.51367980e+03 5.13045059e+03 3.85225361e+17 -2 3.843657524377e+05 1.847080947119e+05 768 7392 3264 6528 5 37 44 4294967291 281 291 291 "" 5.01562464e-02 8.18084445e-01 3.21018877e+02 2.00617889e+03 3.44750753e+03 4.46636632e-01 -1.66839436e-06 1.06658731e-06 6.00000000e+02 2.98552383e+03 2.23643652e+17 3.51305713e+03 5.13045059e+03 3.85233771e+17 -3 4.000000000000e+05 1.563424756231e+04 768 7392 3264 6528 3 22 20 4294967293 156 162 162 "" 5.06901714e-02 8.40850160e-01 3.21018877e+02 2.00612466e+03 3.44750753e+03 4.46622259e-01 -1.62480117e-06 1.07582833e-06 6.00000000e+02 2.98888782e+03 2.23639684e+17 3.51300258e+03 5.13045059e+03 3.85234480e+17 +0 0.000000000000e+00 0.000000000000e+00 768 7392 3264 6528 4 0 0 4294967292 199 207 207 output-entropy_plasticity/solution/solution-00000 1.16744285e-02 4.83332658e-02 3.21018877e+02 1.97425661e+03 2.58226935e+03 4.38175918e-01 -6.70759265e-07 2.04705360e-07 6.00000000e+02 2.95000000e+03 2.22829153e+17 3.51431122e+03 5.13045059e+03 3.85325153e+17 +1 4.000000000000e+05 4.000000000000e+05 768 7392 3264 6528 4 32 28 4294967292 168 176 176 "" 1.14747703e-02 4.84865775e-02 3.21018877e+02 1.97279434e+03 2.58226935e+03 4.37788354e-01 -6.51375447e-07 1.93672617e-07 6.00000000e+02 2.95000000e+03 2.22727088e+17 3.51550463e+03 5.13045059e+03 3.85343592e+17