00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00027 #ifndef OPM_EFF_TO_ABS_LAW_PARAMS_HPP
00028 #define OPM_EFF_TO_ABS_LAW_PARAMS_HPP
00029
00030 #include <cassert>
00031
00032 #include <opm/material/common/EnsureFinalized.hpp>
00033
00034 namespace Opm {
00042 template <class EffLawParamsT, int numPhases>
00043 class EffToAbsLawParams : public EffLawParamsT
00044 {
00045 typedef EffLawParamsT EffLawParams;
00046 typedef typename EffLawParams::Traits::Scalar Scalar;
00047
00048 public:
00049 typedef typename EffLawParams::Traits Traits;
00050
00051
00052 EffToAbsLawParams()
00053 : EffLawParams()
00054 {
00055 for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx)
00056 residualSaturation_[phaseIdx] = 0.0;
00057 }
00058
00063 void finalize()
00064 {
00065 sumResidualSaturations_ = 0.0;
00066 for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx)
00067 sumResidualSaturations_ += residualSaturation_[phaseIdx];
00068
00069 EffLawParams::finalize();
00070 }
00071
00075 Scalar residualSaturation(unsigned phaseIdx) const
00076 { EnsureFinalized::check(); return residualSaturation_[phaseIdx]; }
00077
00081 Scalar sumResidualSaturations() const
00082 { EnsureFinalized::check(); return sumResidualSaturations_; }
00083
00087 void setResidualSaturation(unsigned phaseIdx, Scalar value)
00088 { residualSaturation_[phaseIdx] = value; }
00089
00090 private:
00091
00092 Scalar residualSaturation_[numPhases];
00093 Scalar sumResidualSaturations_;
00094 };
00095
00096 }
00097
00098 #endif