PorousFlowWaterNCG.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#pragma once
 
#include "PorousFlowFluidStateMultiComponentBase.h"
 
class SinglePhaseFluidProperties;
class Water97FluidProperties;
class PorousFlowWaterNCG;
 
template <>
InputParameters validParams<PorousFlowWaterNCG>();
 
class PorousFlowWaterNCG : public PorousFlowFluidStateMultiComponentBase
{
public:
  PorousFlowWaterNCG(const InputParameters & parameters);
 
  virtual std::string fluidStateName() const override;
 
  void thermophysicalProperties(Real pressure,
                                Real temperature,
                                Real Xnacl,
                                Real Z,
                                unsigned int qp,
                                std::vector<FluidStateProperties> & fsp) const override;
  void equilibriumMassFractions(const DualReal & pressure,
                                const DualReal & temperature,
                                DualReal & Xncg,
                                DualReal & Yh2o) const;
 
  void massFractions(const DualReal & pressure,
                     const DualReal & temperature,
                     const DualReal & Z,
                     FluidStatePhaseEnum & phase_state,
                     std::vector<FluidStateProperties> & fsp) const;
 
  void gasProperties(const DualReal & pressure,
                     const DualReal & temperature,
                     std::vector<FluidStateProperties> & fsp) const;
  void liquidProperties(const DualReal & pressure,
                        const DualReal & temperature,
                        std::vector<FluidStateProperties> & fsp) const;
 
  DualReal liquidDensity(const DualReal & pressure, const DualReal & temperature) const;
 
  DualReal gasDensity(const DualReal & pressure,
                      const DualReal & temperature,
                      std::vector<FluidStateProperties> & fsp) const;
 
  DualReal saturation(const DualReal & pressure,
                      const DualReal & temperature,
                      const DualReal & Z,
                      std::vector<FluidStateProperties> & fsp) const;
 
  void twoPhaseProperties(const DualReal & pressure,
                          const DualReal & temperature,
                          const DualReal & Z,
                          unsigned int qp,
                          std::vector<FluidStateProperties> & fsp) const;
 
  DualReal enthalpyOfDissolution(const DualReal & temperature) const;
 
  virtual Real totalMassFraction(
      Real pressure, Real temperature, Real Xnacl, Real saturation, unsigned int qp) const override;
 
protected:
  DualReal moleFractionToMassFraction(const DualReal & xmol) const;
 
  void checkVariables(Real temperature) const;
 
  const SinglePhaseFluidProperties & _water_fp;
  const Water97FluidProperties & _water97_fp;
  const SinglePhaseFluidProperties & _ncg_fp;
  const Real _Mh2o;
  const Real _Mncg;
  const Real _water_triple_temperature;
  const Real _water_critical_temperature;
  const std::vector<Real> _ncg_henry;
};