more well state output

Author: GitPaean

opm/simulators/wells/MSWellHelpers.cpp

@@ -167,9 +167,6 @@ applyUMFPack(Dune::UMFPack<MatrixType>& linsolver,
         // Solve
         linsolver.apply(y, x, res);
 
-        const std::string msg {"tesing numerical problem in UMFPack"};
-        OPM_THROW_NOLOG(NumericalProblem, msg);
-
         // Checking if there is any inf or nan in y
         // it will be the solution before we find a way to catch the singularity of the matrix
         for (std::size_t i_block = 0; i_block < y.size(); ++i_block) {

opm/simulators/wells/MultisegmentWell.hpp

@@ -349,7 +349,8 @@ namespace Opm {
 
         FSInfo getFirstPerforationFluidStateInfo(const Simulator& simulator) const;
 
-        void outputDebugInfoNumericalProblem(DeferredLogger& deferred_logger) const;
+        void outputDebugInfoNumericalProblem(const SingleWellState<Scalar, IndexTraits>& ws,
+                                             DeferredLogger& deferred_logger) const;
     };
 
 }

opm/simulators/wells/MultisegmentWell_impl.hpp

@@ -283,7 +283,8 @@ namespace Opm
             // this is the process with rank zero)
             deferred_logger.problem("In MultisegmentWell::recoverWellSolutionAndUpdateWellState for well "
                                     + this->name() +": "+exp.what());
-            this->outputDebugInfoNumericalProblem(deferred_logger);
+            const auto& ws = well_state.well(this->index_of_well_);
+            this->outputDebugInfoNumericalProblem(ws, deferred_logger);
             throw;
         }
     }
@@ -633,7 +634,8 @@ namespace Opm
             // this is the process with rank zero)
             deferred_logger.problem("In MultisegmentWell::solveEqAndUpdateWellState for well "
                                     + this->name() +": "+exp.what());
-            this->outputDebugInfoNumericalProblem(deferred_logger);
+            const auto& ws = well_state.well(this->index_of_well_);
+            this->outputDebugInfoNumericalProblem(ws, deferred_logger);
             throw;
         }
     }
@@ -1609,7 +1611,8 @@ namespace Opm
                 // this is the process with rank zero)
                 deferred_logger.problem("In MultisegmentWell::iterateWellEqWithControl for well "
                                         + this->name() +": "+exp.what());
-                this->outputDebugInfoNumericalProblem(deferred_logger);
+                const auto& ws = well_state.well(this->index_of_well_);
+                this->outputDebugInfoNumericalProblem(ws, deferred_logger);
                 throw;
             }
         }
@@ -1791,7 +1794,8 @@ namespace Opm
                 // this is the process with rank zero)
                 deferred_logger.problem("In MultisegmentWell::iterateWellEqWithSwitching for well "
                                         + this->name() +": "+exp.what());
-                this->outputDebugInfoNumericalProblem(deferred_logger);
+                const auto& ws = well_state.well(this->index_of_well_);
+                this->outputDebugInfoNumericalProblem(ws, deferred_logger);
                 throw;
             }
         }
@@ -2386,9 +2390,11 @@ namespace Opm
     template <typename TypeTag>
     void
     MultisegmentWell<TypeTag>::
-    outputDebugInfoNumericalProblem(DeferredLogger& deferred_logger) const
+    outputDebugInfoNumericalProblem(const SingleWellState<Scalar, IndexTraits>& ws,
+                                    DeferredLogger& deferred_logger) const
     {
-        std::string msg = this->primary_variables_.debugInfo();
+        std::string msg = ws.debugInfo();
+        msg += this->primary_variables_.debugInfo();
         deferred_logger.debug(msg);
     }
 

opm/simulators/wells/SegmentState.cpp

@@ -31,6 +31,8 @@
 #include <stdexcept>
 #include <vector>
 
+#include <fmt/format.h>
+
 namespace {
 
 std::vector<int> make_segment_number(const Opm::WellSegments& segments)
@@ -147,6 +149,27 @@ bool SegmentState<Scalar>::operator==(const SegmentState& rhs) const
            this->m_segment_number == rhs.m_segment_number;
 }
 
+template<class Scalar>
+std::string
+SegmentState<Scalar>::debugInfo() const
+{
+    if (this->empty()) {
+        return "";
+    }
+    const std::size_t num_phases = this->rates.size() / this->size();
+    std::string info = "SegmentState:\n";
+    for (std::size_t i = 0; i < this->size(); ++i) {
+        info += fmt::format("  Segment {:4}: Segment number {:4}, Pressure = {:8.2e} Pa, Rate: ",
+                            i, this->m_segment_number[i], this->pressure[i]);
+        for (std::size_t p = 0; p < num_phases; ++p) {
+            info += fmt::format(" {: 8.2e}", this->rates[i * num_phases + p]);
+        }
+        info += fmt::format("\n");
+    }
+    return info;
+}
+
+
 template class SegmentState<double>;
 
 #if FLOW_INSTANTIATE_FLOAT

opm/simulators/wells/SegmentState.hpp

@@ -22,6 +22,7 @@
 
 #include <cstddef>
 #include <vector>
+#include <string>
 
 namespace Opm {
 class WellSegments;
@@ -104,6 +105,7 @@ class SegmentState
     std::vector<Scalar> pressure_drop_hydrostatic;
     std::vector<Scalar> pressure_drop_accel;
 
+    std::string debugInfo() const;
 private:
     std::vector<int>    m_segment_number;
 };

opm/simulators/wells/SingleWellState.cpp

@@ -423,10 +423,34 @@ SingleWellState<Scalar, IndexTraits>::
 debugInfo() const
 {
     std::string info = "Well name: " + this->name + " well state:\n";
-    fmt::format_to(std::back_inserter(info), " type: {}, staus: {}, control type: \n",
+    fmt::format_to(std::back_inserter(info), " type: {}, staus: {}, control type: {}, BHP: {:8.2e} Pa, THP: {:8.2e} Pa\n",
                    this->producer? "Producer" : " Injector", WellStatus2String(this->status),
                    this->producer ? WellProducerCMode2String(this->production_cmode)
-                                  : WellInjectorCMode2String(this->injection_cmode));
+                                  : WellInjectorCMode2String(this->injection_cmode),
+                   this->bhp / Opm::unit::barsa,
+                   this->thp / Opm::unit::barsa);
+
+    fmt::format_to(std::back_inserter(info), "  Surface rates (m3/s): ");
+    for (const auto& r : this->surface_rates) {
+        fmt::format_to(std::back_inserter(info), " {:8.2e}", r);
+    }
+    fmt::format_to(std::back_inserter(info), "\n");
+    fmt::format_to(std::back_inserter(info), "  Connection pressures and connection rates:\n");
+    for (std::size_t perf = 0; perf < this->perf_data.size(); perf++) {
+        fmt::format_to(std::back_inserter(info),
+                       "  Connection {:4}: Pressure: {:8.2e} Pa, Rate:",
+                       perf,
+                       this->perf_data.pressure[perf]);
+        const auto& connection_rates = this->perf_data.phase_rates;
+        const std::size_t num_phases = this->pu.numActivePhases();
+
+        for (std::size_t p = 0; p < num_phases; ++p) {
+            fmt::format_to(std::back_inserter(info), " {: 8.2e} m3/s", connection_rates[perf * num_phases + p]);
+        }
+        fmt::format_to(std::back_inserter(info), "\n");
+    }
+
+    info += this->segments.debugInfo();
 
     return info;
 }