8ff84920ed
The current Network::Balance class applies default values a little too eagerly and this leads to problems for restarted simulation. We must take care to preserve the "did this item have a value" state also in the restart file lest interactions with TUNING be lost there. To this end, make the minimum timestep size for network balancing be an std::optional<> and use std::nullopt to represent "no value set" in the input file. While here, also preserve negative balancing intervals exactly instead of replacing these values by zero.
108 lines
2.9 KiB
C++
108 lines
2.9 KiB
C++
/*
|
|
Copyright 2018 Statoil ASA.
|
|
|
|
This file is part of the Open Porous Media Project (OPM).
|
|
|
|
OPM is free software: you can redistribute it and/or modify
|
|
it under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation, either version 3 of the License, or
|
|
(at your option) any later version.
|
|
|
|
OPM is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
GNU General Public License for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
|
*/
|
|
|
|
#ifndef OPM_DOUBHEAD_HEADER_INCLUDED
|
|
#define OPM_DOUBHEAD_HEADER_INCLUDED
|
|
|
|
#include <chrono>
|
|
#include <cstddef>
|
|
#include <vector>
|
|
|
|
namespace Opm {
|
|
struct Tuning;
|
|
class Schedule;
|
|
class UDQParams;
|
|
class UnitSystem;
|
|
}
|
|
|
|
namespace Opm { namespace RestartIO {
|
|
|
|
class DoubHEAD
|
|
{
|
|
public:
|
|
struct TimeStamp {
|
|
std::chrono::time_point<std::chrono::system_clock> start;
|
|
std::chrono::duration<double, std::chrono::seconds::period> elapsed;
|
|
};
|
|
|
|
struct guideRate {
|
|
double A;
|
|
double B;
|
|
double C;
|
|
double D;
|
|
double E;
|
|
double F;
|
|
double delay;
|
|
double damping_fact;
|
|
};
|
|
|
|
struct liftOptPar {
|
|
double min_int;
|
|
double incr;
|
|
double min_ec_grad;
|
|
};
|
|
|
|
struct NetBalanceParams {
|
|
explicit NetBalanceParams(const UnitSystem& usys);
|
|
|
|
double balancingInterval;
|
|
double convTolNodPres;
|
|
double convTolTHPCalc;
|
|
double targBranchBalError;
|
|
double maxBranchBalError;
|
|
double minTimeStepSize;
|
|
};
|
|
|
|
DoubHEAD();
|
|
|
|
~DoubHEAD() = default;
|
|
DoubHEAD(const DoubHEAD& rhs) = default;
|
|
DoubHEAD(DoubHEAD&& rhs) = default;
|
|
|
|
DoubHEAD& operator=(const DoubHEAD& rhs) = default;
|
|
DoubHEAD& operator=(DoubHEAD&& rhs) = default;
|
|
|
|
DoubHEAD& tuningParameters(const Tuning& tuning,
|
|
const double cnvT);
|
|
|
|
DoubHEAD& timeStamp(const TimeStamp& ts);
|
|
DoubHEAD& nextStep(const double nextTimeStep);
|
|
|
|
DoubHEAD& drsdt(const Schedule& sched,
|
|
const std::size_t lookup_step,
|
|
const double cnvT);
|
|
|
|
DoubHEAD& udq_param(const UDQParams& udqPar);
|
|
DoubHEAD& guide_rate_param(const guideRate& guide_rp);
|
|
DoubHEAD& lift_opt_param(const liftOptPar& lo_par);
|
|
DoubHEAD& netBalParams(const NetBalanceParams& net_bal_par);
|
|
|
|
const std::vector<double>& data() const
|
|
{
|
|
return this->data_;
|
|
}
|
|
|
|
private:
|
|
std::vector<double> data_;
|
|
};
|
|
|
|
}} // Opm::RestartIO
|
|
|
|
#endif // OPM_DOUBHEAD_HEADER_INCLUDED
|