blackoilpolymerparams.hpp

Go to the documentation of this file.

// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
/*
  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 2 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/>.

  Consult the COPYING file in the top-level source directory of this
  module for the precise wording of the license and the list of
  copyright holders.
*/
#ifndef OPM_BLACK_OIL_POLYMER_PARAMS_HPP
#define OPM_BLACK_OIL_POLYMER_PARAMS_HPP

#include <opm/material/common/Tabulated1DFunction.hpp>
#include <opm/material/common/IntervalTabulated2DFunction.hpp>

#include <map>
#include <vector>

namespace Opm {

class EclipseState;

template<class Scalar>
struct BlackOilPolymerParams {
    using TabulatedFunction = Tabulated1DFunction<Scalar>;
    using TabulatedTwoDFunction = IntervalTabulated2DFunction<Scalar>;

    enum AdsorptionBehaviour { Desorption = 1, NoDesorption = 2 };

    template<bool enablePolymer, bool enablePolymerMolarWeight>
    void initFromState(const EclipseState& eclState);

    // a struct containing the constants to calculate polymer viscosity
    // based on Mark-Houwink equation and Huggins equation, the constants are provided
    // by the keyword PLYVMH
    struct PlyvmhCoefficients {
        Scalar k_mh;
        Scalar a_mh;
        Scalar gamma;
        Scalar kappa;
    };

    struct SkprpolyTable {
        double refConcentration;
        TabulatedTwoDFunction table_func;
    };

    std::vector<Scalar> plyrockDeadPoreVolume_{};
    std::vector<Scalar> plyrockResidualResistanceFactor_{};
    std::vector<Scalar> plyrockRockDensityFactor_{};
    std::vector<Scalar> plyrockAdsorbtionIndex_{};
    std::vector<Scalar> plyrockMaxAdsorbtion_{};
    std::vector<TabulatedFunction> plyadsAdsorbedPolymer_{};
    std::vector<TabulatedFunction> plyviscViscosityMultiplierTable_{};
    std::vector<Scalar> plymaxMaxConcentration_{};
    std::vector<Scalar> plymixparToddLongstaff_{};
    std::vector<std::vector<Scalar>> plyshlogShearEffectRefMultiplier_{};
    std::vector<std::vector<Scalar>> plyshlogShearEffectRefLogVelocity_{};
    std::vector<Scalar> shrate_{};
    bool hasShrate_ = false;
    bool hasPlyshlog_ = false;

    std::vector<PlyvmhCoefficients> plyvmhCoefficients_{};
    std::map<int, TabulatedTwoDFunction> plymwinjTables_{};
    std::map<int, TabulatedTwoDFunction> skprwatTables_{};

    std::map<int, SkprpolyTable> skprpolyTables_{};

private:
    void setNumSatRegions(unsigned numRegions);

    void setNumMixRegions(unsigned numRegions, bool enablePolymerMolarWeight);

    void setPlyrock(unsigned satRegionIdx,
                    const Scalar& plyrockDeadPoreVolume,
                    const Scalar& plyrockResidualResistanceFactor,
                    const Scalar& plyrockRockDensityFactor,
                    const Scalar& plyrockAdsorbtionIndex,
                    const Scalar& plyrockMaxAdsorbtion);

    void processPlyrock(const EclipseState& eclState);

    void processPlyads(const EclipseState& eclState);

    template<bool enablePolymerMolarWeight>
    void processPlyvisc(const EclipseState& eclState);

    template<bool enablePolymerMolarWeight>
    void processPlymax(const EclipseState& eclState);

    template<bool enablePolymerMolarWeight>
    void processPlmixpar(const EclipseState& eclState);

    void processPlyshlog(const EclipseState& eclState);

    void processShrate(const EclipseState& eclState);

    void processPlyvmh(const EclipseState& eclState);

    void processPlymwinj(const EclipseState& eclState);

    void processSkprwat(const EclipseState& eclState);

    void processSkprpoly(const EclipseState& eclState);
};

} // namespace Opm

#endif // OPM_BLACK_OIL_POLYMER_PARAMS_HPP