My Project
EclSolidEnergyLawMultiplexer.hpp
Go to the documentation of this file.
1 // -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
2 // vi: set et ts=4 sw=4 sts=4:
3 /*
4  This file is part of the Open Porous Media project (OPM).
5 
6  OPM is free software: you can redistribute it and/or modify
7  it under the terms of the GNU General Public License as published by
8  the Free Software Foundation, either version 2 of the License, or
9  (at your option) any later version.
10 
11  OPM is distributed in the hope that it will be useful,
12  but WITHOUT ANY WARRANTY; without even the implied warranty of
13  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  GNU General Public License for more details.
15 
16  You should have received a copy of the GNU General Public License
17  along with OPM. If not, see <http://www.gnu.org/licenses/>.
18 
19  Consult the COPYING file in the top-level source directory of this
20  module for the precise wording of the license and the list of
21  copyright holders.
22 */
27 #ifndef OPM_ECL_SOLID_ENERGY_LAW_MULTIPLEXER_HPP
28 #define OPM_ECL_SOLID_ENERGY_LAW_MULTIPLEXER_HPP
29 
31 
32 #include "EclHeatcrLaw.hpp"
33 #include "EclSpecrockLaw.hpp"
34 #include "NullSolidEnergyLaw.hpp"
35 
37 
38 namespace Opm
39 {
45 template <class ScalarT,
46  class FluidSystem,
47  class ParamsT = EclSolidEnergyLawMultiplexerParams<ScalarT>>
49 {
50  enum { numPhases = FluidSystem::numPhases };
51 
55 
56 public:
57  typedef ParamsT Params;
58  typedef typename Params::Scalar Scalar;
59 
63  template <class FluidState, class Evaluation = typename FluidState::Scalar>
64  static Evaluation solidInternalEnergy(const Params& params, const FluidState& fluidState)
65  {
66  switch (params.solidEnergyApproach()) {
67  case Params::heatcrApproach:
68  // relevant ECL keywords: HEATCR, HEATCRT and STCOND
69  return HeatcrLaw::solidInternalEnergy(params.template getRealParams<Params::heatcrApproach>(), fluidState);
70 
71  case Params::specrockApproach:
72  // relevant ECL keyword: SPECROCK
73  return SpecrockLaw::solidInternalEnergy(params.template getRealParams<Params::specrockApproach>(), fluidState);
74 
75  case Params::nullApproach:
76  // (no relevant ECL keyword)
77  return NullLaw::solidInternalEnergy(0, fluidState);
78 
79  default:
80  throw std::logic_error("Invalid solid energy approach: "+std::to_string(int(params.solidEnergyApproach())));
81  }
82  }
83 };
84 } // namespace Opm
85 
86 #endif
Implements the volumetric interior energy relations of rock used by ECL.
The default implementation of a parameter object for the ECL thermal law.
Implements the volumetric interior energy relations of rock used by ECL.
A number of commonly used algebraic functions for the localized OPM automatic differentiation (AD) fr...
Implements a solid energy storage law which just returns 0.
Implements the volumetric interior energy relations of rock used by ECL.
Definition: EclHeatcrLaw.hpp:47
static Evaluation solidInternalEnergy(const Params &params, const FluidState &fluidState)
Given a fluid state, compute the volumetric internal energy of the rock [W/m^3].
Definition: EclHeatcrLaw.hpp:56
Provides the energy storage relation of rock.
Definition: EclSolidEnergyLawMultiplexer.hpp:49
static Evaluation solidInternalEnergy(const Params &params, const FluidState &fluidState)
Given a fluid state, compute the volumetric internal energy of the rock [W/m^3].
Definition: EclSolidEnergyLawMultiplexer.hpp:64
Implements the volumetric interior energy relations of rock used by ECL.
Definition: EclSpecrockLaw.hpp:46
static Evaluation solidInternalEnergy(const Params &params, const FluidState &fluidState)
Given a fluid state, compute the volumetric internal energy of the rock [W/m^3].
Definition: EclSpecrockLaw.hpp:55
Implements a solid energy storage law which just returns 0.
Definition: NullSolidEnergyLaw.hpp:41
static Evaluation solidInternalEnergy(const Params &, const FluidState &)
Given a fluid state, compute the volumetric internal energy of the solid matrix [W/m^3].
Definition: NullSolidEnergyLaw.hpp:53