12#ifndef DUMUX_BRINE_AIR_FLUID_SYSTEM_HH 
   13#define DUMUX_BRINE_AIR_FLUID_SYSTEM_HH 
   41template<
bool fastButSimplifiedRelations = false>
 
   57template <
class Scalar,
 
   59          class Policy = BrineAirDefaultPolicy<>>
 
   61: 
public Base<Scalar, BrineAir<Scalar, H2Otype, Policy>>
 
  104    struct BrineAdapterPolicy
 
  106        using FluidSystem = 
Brine;
 
  108        static constexpr int phaseIdx(
int brinePhaseIdx) { 
return liquidPhaseIdx; }
 
  109        static constexpr int compIdx(
int brineCompIdx)
 
  111            switch (brineCompIdx)
 
  120    template<
class Flu
idState>
 
  135        assert(0 <= phaseIdx && phaseIdx < 
numPhases);
 
  141        DUNE_THROW(Dune::InvalidStateException, 
"Invalid phase index " << phaseIdx);
 
 
  154    static constexpr bool isGas(
int phaseIdx)
 
  156        assert(0 <= phaseIdx && phaseIdx < 
numPhases);
 
 
  176        assert(0 <= phaseIdx && phaseIdx < 
numPhases);
 
 
  194        assert(0 <= phaseIdx && phaseIdx < 
numPhases);
 
 
  209        assert(0 <= phaseIdx && phaseIdx < 
numPhases);
 
 
  222        assert(0 <= phaseIdx && phaseIdx < 
numPhases);
 
 
  241            case H2OIdx: 
return H2O::name();
 
  245        DUNE_THROW(Dune::InvalidStateException, 
"Invalid component index " << compIdx);
 
 
  257            case H2OIdx: 
return H2O::molarMass();
 
  261        DUNE_THROW(Dune::InvalidStateException, 
"Invalid component index " << compIdx);
 
 
  270    template <
class Flu
idState>
 
  277            if (!Policy::useKelvinVaporPressure())
 
  283                const auto pc =  (fluidState.wettingPhase() == (
int) 
liquidPhaseIdx)
 
  292            DUNE_THROW(Dune::NotImplemented, 
"NaCl::vaporPressure(t)");
 
  294            DUNE_THROW(Dune::NotImplemented, 
"Invalid component index " << compIdx);
 
 
  330        std::cout << 
"The brine-air fluid system was configured with the following policy:\n";
 
  331        std::cout << 
" - use brine density as liquid mixture density: " << std::boolalpha << Policy::useBrineDensityAsLiquidMixtureDensity() << 
"\n";
 
  332        std::cout << 
" - use ideal gas density: " << std::boolalpha << Policy::useIdealGasDensity() << std::endl;
 
  334        if (H2O::isTabulated)
 
  335            H2O::init(tempMin, tempMax, nTemp, pressMin, pressMax, nPress);
 
 
  352    template <
class Flu
idState>
 
  355        assert(0 <= phaseIdx && phaseIdx < 
numPhases);
 
  357        const auto T = fluidState.temperature(phaseIdx);
 
  358        const auto p = fluidState.pressure(phaseIdx);
 
  363            if (Policy::useBrineDensityAsLiquidMixtureDensity())
 
  376            if (Policy::useIdealGasDensity())
 
  381            return H2O::gasDensity(T, fluidState.partialPressure(
phase1Idx, 
H2OIdx))
 
  385            DUNE_THROW(Dune::InvalidStateException, 
"Invalid phase index " << phaseIdx);
 
 
  390    template <
class Flu
idState>
 
  397            const Scalar T = fluidState.temperature(phaseIdx);
 
  400            if (Policy::useIdealGasDensity())
 
  405            return H2O::gasMolarDensity(T, fluidState.partialPressure(
phase1Idx, 
H2OIdx))
 
  409            DUNE_THROW(Dune::InvalidStateException, 
"Invalid phase index " << phaseIdx);
 
 
  423    template <
class Flu
idState>
 
  426        assert(0 <= phaseIdx && phaseIdx < 
numPhases);
 
  431            return Air::gasViscosity(fluidState.temperature(phaseIdx), fluidState.pressure(phaseIdx));
 
 
  456    template <
class Flu
idState>
 
  459        assert(0 <= phaseIdx && phaseIdx < 
numPhases);
 
  462        Scalar T = fluidState.temperature(phaseIdx);
 
  463        Scalar p = fluidState.pressure(phaseIdx);
 
  477                return H2O::vaporPressure(T)/p;
 
  479            else if (compIdx == 
AirIdx)
 
  487        DUNE_THROW(Dune::InvalidStateException, 
"Invalid phase index " << phaseIdx);
 
 
  492    template <
class Flu
idState>
 
  497        DUNE_THROW(Dune::NotImplemented, 
"Diffusion coefficients");
 
 
  502    template <
class Flu
idState>
 
  508        assert(0 <= phaseIdx && phaseIdx < 
numPhases);
 
  512        const auto T = fluidState.temperature(phaseIdx);
 
  513        const auto p = fluidState.pressure(phaseIdx);
 
  515        if (compIIdx > compJIdx)
 
  516            std::swap(compIIdx, compJIdx);
 
  525                DUNE_THROW(Dune::NotImplemented, 
"Binary diffusion coefficient of components " 
  526                                                 << compIIdx << 
" and " << compJIdx
 
  527                                                 << 
" in phase " << phaseIdx);
 
  541                DUNE_THROW(Dune::NotImplemented, 
"Binary diffusion coefficient of components " 
  542                                                 << compIIdx << 
" and " << compJIdx
 
  543                                                 << 
" in phase " << phaseIdx);
 
  546        DUNE_THROW(Dune::InvalidStateException, 
"Invalid phase index " << phaseIdx);
 
 
  569    template <
class Flu
idState>
 
  572        assert(0 <= phaseIdx && phaseIdx < 
numPhases);
 
  574        const Scalar T = fluidState.temperature(phaseIdx);
 
  575        const Scalar p = fluidState.pressure(phaseIdx);
 
  587        DUNE_THROW(Dune::InvalidStateException, 
"Invalid phase index " << phaseIdx);
 
 
  596    template <
class Flu
idState>
 
  603            DUNE_THROW(Dune::NotImplemented, 
"The component enthalpies in the liquid phase are not implemented.");
 
  607            if (componentIdx == 
H2OIdx)
 
  608                return H2O::gasEnthalpy(T, p);
 
  609            else if (componentIdx == 
AirIdx)
 
  611            else if (componentIdx == 
NaClIdx)
 
  612                DUNE_THROW(Dune::InvalidStateException, 
"Implementation assumes NaCl not to be present in gas phase");
 
  613            DUNE_THROW(Dune::InvalidStateException, 
"Invalid component index " << componentIdx);
 
  615        DUNE_THROW(Dune::InvalidStateException, 
"Invalid phase index " << phaseIdx);
 
 
  630    template <
class Flu
idState>
 
  638            + H2O::gasThermalConductivity(fluidState.temperature(phaseIdx), fluidState.pressure(phaseIdx)) * fluidState.massFraction(
gasPhaseIdx, 
H2OIdx);
 
  640        DUNE_THROW(Dune::InvalidStateException, 
"Invalid phase index " << phaseIdx);
 
 
  654    template <
class Flu
idState>
 
  657        const Scalar T = fluidState.temperature(phaseIdx);
 
  658        const Scalar p = fluidState.pressure(phaseIdx);
 
  668        DUNE_THROW(Dune::InvalidStateException, 
"Invalid phase index " << phaseIdx);
 
 
 
Adapter class for fluid states with different indices.
A simple class for the air fluid properties.
Binary coefficients for water and air.
Definition h2o_air.hh:24
static Scalar henry(Scalar temperature)
Henry coefficient  for air in liquid water.
Definition h2o_air.hh:35
static Scalar gasDiffCoeff(Scalar temperature, Scalar pressure)
Binary diffusion coefficient  for molecular water and air.
Definition h2o_air.hh:55
static Scalar liquidDiffCoeff(Scalar temperature, Scalar pressure)
Diffusion coefficient  for molecular nitrogen in liquid water.
Definition h2o_air.hh:88
A class for the air fluid properties.
Definition air.hh:33
static Scalar gasDensity(Scalar temperature, Scalar pressure)
The density  of Air at a given pressure and temperature.
Definition air.hh:71
static constexpr Scalar molarMass()
The molar mass in  of Air.
Definition air.hh:48
static const Scalar gasHeatCapacity(Scalar temperature, Scalar pressure)
Specific isobaric heat capacity  of pure air.
Definition air.hh:299
static Scalar gasViscosity(Scalar temperature, Scalar pressure)
The dynamic viscosity  of Air at a given pressure and temperature.
Definition air.hh:179
static Scalar gasThermalConductivity(Scalar temperature, Scalar pressure)
Thermal conductivity  of air.
Definition air.hh:337
static constexpr bool gasIsIdeal()
Returns true, the gas phase is assumed to be ideal.
Definition air.hh:95
static Scalar gasEnthalpy(Scalar temperature, Scalar pressure)
Specific enthalpy of Air  with 273.15  as basis.
Definition air.hh:262
static std::string name()
A human readable name for Air.
Definition air.hh:40
static Scalar gasMolarDensity(Scalar temperature, Scalar pressure)
The molar density of air in , depending on pressure and temperature.
Definition air.hh:83
A class for the NaCl properties.
Definition nacl.hh:34
static std::string name()
A human readable name for the NaCl.
Definition nacl.hh:39
static constexpr Scalar molarMass()
The molar mass of NaCl in .
Definition nacl.hh:47
Tabulates all thermodynamic properties of a given component.
Definition tabulatedcomponent.hh:672
static constexpr Scalar R
The ideal gas constant .
Definition constants.hh:32
Adapter class for fluid states with different indices.
Definition adapter.hh:32
Fluid system base class.
Definition fluidsystems/base.hh:32
ScalarType Scalar
export the scalar type
Definition fluidsystems/base.hh:35
A compositional two-phase fluid system with a liquid and a gaseous phase and ,  and  (dissolved miner...
Definition brineair.hh:62
static Scalar componentEnthalpy(const FluidState &fluidState, int phaseIdx, int componentIdx)
Returns the specific enthalpy  of a component in a specific phase.
Definition brineair.hh:597
static Scalar density(const FluidState &fluidState, int phaseIdx)
Given a phase's composition, temperature, pressure, and the partial pressures of all components,...
Definition brineair.hh:353
static Scalar thermalConductivity(const FluidState &fluidState, int phaseIdx)
Thermal conductivity of a fluid phase .
Definition brineair.hh:631
static constexpr int comp0Idx
Definition brineair.hh:99
static constexpr bool isCompressible(int phaseIdx)
Returns true if and only if a fluid phase is assumed to be compressible.
Definition brineair.hh:192
static constexpr int liquidPhaseIdx
Definition brineair.hh:87
static constexpr int AirIdx
Definition brineair.hh:97
static constexpr int NaClIdx
Definition brineair.hh:98
static void init()
Initialize the fluid system's static parameters generically.
Definition brineair.hh:306
static constexpr int numPhases
Definition brineair.hh:84
static Scalar heatCapacity(const FluidState &fluidState, int phaseIdx)
Specific isobaric heat capacity of a fluid phase. .
Definition brineair.hh:655
BinaryCoeff::H2O_Air H2O_Air
Definition brineair.hh:76
static constexpr int numComponents
Definition brineair.hh:85
NullParameterCache ParameterCache
Definition brineair.hh:79
static Scalar binaryDiffusionCoefficient(const FluidState &fluidState, int phaseIdx, int compIIdx, int compJIdx)
Given a phase's composition, temperature and pressure, return the binary diffusion coefficient  for c...
Definition brineair.hh:503
H2Otype H2O
Definition brineair.hh:68
static Scalar viscosity(const FluidState &fluidState, int phaseIdx)
Calculate the dynamic viscosity of a fluid phase .
Definition brineair.hh:424
static constexpr int comp2Idx
Definition brineair.hh:101
static Scalar enthalpy(const FluidState &fluidState, int phaseIdx)
Given a phase's composition, temperature and pressure, return its specific enthalpy .
Definition brineair.hh:570
static Scalar molarDensity(const FluidState &fluidState, int phaseIdx)
Calculate the molar density  of a fluid phase.
Definition brineair.hh:391
static constexpr int getMainComponent(int phaseIdx)
Get the main component of a given phase if possible.
Definition brineair.hh:220
Components::Air< Scalar > Air
Definition brineair.hh:69
static constexpr bool isGas(int phaseIdx)
Return whether a phase is gaseous.
Definition brineair.hh:154
static Scalar fugacityCoefficient(const FluidState &fluidState, int phaseIdx, int compIdx)
Returns the fugacity coefficient  of a component in a phase.
Definition brineair.hh:457
static void init(Scalar tempMin, Scalar tempMax, unsigned nTemp, Scalar pressMin, Scalar pressMax, unsigned nPress)
Initialize the fluid system's static parameters using problem specific temperature and pressure range...
Definition brineair.hh:327
static Scalar molarMass(int compIdx)
Return the molar mass of a component in .
Definition brineair.hh:252
static Scalar vaporPressure(const FluidState &fluidState, int compIdx)
Vapor pressure of a component .
Definition brineair.hh:271
static std::string componentName(int compIdx)
Return the human readable name of a component.
Definition brineair.hh:236
static std::string phaseName(int phaseIdx)
Return the human readable name of a fluid phase.
Definition brineair.hh:133
static constexpr int comp1Idx
Definition brineair.hh:100
static constexpr int phase1Idx
Definition brineair.hh:91
static constexpr int phase0Idx
Definition brineair.hh:90
static constexpr bool isMiscible()
Returns whether the fluids are miscible.
Definition brineair.hh:147
static Scalar diffusionCoefficient(const FluidState &fluidState, int phaseIdx, int compIdx)
Calculate the binary molecular diffusion coefficient for a component in a fluid phase .
Definition brineair.hh:493
static bool isIdealGas(int phaseIdx)
Returns true if and only if a fluid phase is assumed to be an ideal gas.
Definition brineair.hh:207
static constexpr int gasPhaseIdx
Definition brineair.hh:88
static bool isIdealMixture(int phaseIdx)
Returns true if and only if a fluid phase is assumed to be an ideal mixture.
Definition brineair.hh:174
Components::NaCl< Scalar > NaCl
Definition brineair.hh:70
Dumux::FluidSystems::Brine< Scalar, H2Otype > Brine
Definition brineair.hh:73
static constexpr int H2OIdx
Definition brineair.hh:96
A compositional single phase fluid system consisting of two components, which are H2O and NaCl.
Definition fluidsystems/brine.hh:34
static Scalar thermalConductivity(const FluidState &fluidState, int phaseIdx)
Definition fluidsystems/brine.hh:434
static Scalar viscosity(const FluidState &fluidState, int phaseIdx=liquidPhaseIdx)
Definition fluidsystems/brine.hh:275
static constexpr int H2OIdx
Definition fluidsystems/brine.hh:48
static Scalar heatCapacity(const FluidState &fluidState, int phaseIdx)
Definition fluidsystems/brine.hh:450
static Scalar enthalpy(const FluidState &fluidState, int phaseIdx)
Definition fluidsystems/brine.hh:335
static Scalar vaporPressure(const FluidState &fluidState, int compIdx)
Definition fluidsystems/brine.hh:306
static constexpr int NaClIdx
Definition fluidsystems/brine.hh:49
static Scalar molarDensity(const FluidState &fluidState, int phaseIdx=liquidPhaseIdx)
Definition fluidsystems/brine.hh:372
static constexpr int liquidPhaseIdx
Definition fluidsystems/brine.hh:46
static bool isCompressible(int phaseIdx=liquidPhaseIdx)
Definition fluidsystems/brine.hh:111
static Scalar density(const FluidState &fluidState, int phaseIdx=liquidPhaseIdx)
Definition fluidsystems/brine.hh:217
static Scalar binaryDiffusionCoefficient(const FluidState &fluidState, int phaseIdx, int compIIdx, int compJIdx)
Definition fluidsystems/brine.hh:401
Relations valid for an ideal gas.
Definition idealgas.hh:25
static constexpr Scalar density(Scalar avgMolarMass, Scalar temperature, Scalar pressure)
The density of the gas in , depending on pressure, temperature and average molar mass of the gas.
Definition idealgas.hh:37
static constexpr Scalar molarDensity(Scalar temperature, Scalar pressure)
The molar density of the gas , depending on pressure and temperature.
Definition idealgas.hh:58
The a parameter cache which does nothing.
Definition nullparametercache.hh:22
Some exceptions thrown in DuMux
A fluid system for brine, i.e. H2O with dissolved NaCl.
Material properties of pure water .
Binary coefficients for water and air.
Relations valid for an ideal gas.
Material properties of pure salt .
A collection of input/output field names for common physical quantities.
std::string gaseousPhase() noexcept
I/O name of gaseous phase.
Definition name.hh:111
std::string liquidPhase() noexcept
I/O name of liquid phase.
Definition name.hh:107
Policy for the brine-air fluid system.
Definition brineair.hh:43
static constexpr bool useIdealGasDensity()
Definition brineair.hh:45
static constexpr bool useKelvinVaporPressure()
Definition brineair.hh:46
static constexpr bool useBrineDensityAsLiquidMixtureDensity()
Definition brineair.hh:44
Tabulates all thermodynamic properties of a given untabulated chemical species.