PorousFlowWaterNCG.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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 : public PorousFlowFluidStateMultiComponentBase
{
public:
  static InputParameters validParams();
 
  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 thermophysicalProperties(const ADReal & pressure,
                                const ADReal & temperature,
                                const ADReal & Xnacl,
                                const ADReal & Z,
                                unsigned int qp,
                                std::vector<FluidStateProperties> & fsp) const override;
  void equilibriumMassFractions(const ADReal & pressure,
                                const ADReal & temperature,
                                ADReal & Xncg,
                                ADReal & Yh2o) const;
 
  void massFractions(const ADReal & pressure,
                     const ADReal & temperature,
                     const ADReal & Z,
                     FluidStatePhaseEnum & phase_state,
                     std::vector<FluidStateProperties> & fsp) const;
 
  void gasProperties(const ADReal & pressure,
                     const ADReal & temperature,
                     std::vector<FluidStateProperties> & fsp) const;
  void liquidProperties(const ADReal & pressure,
                        const ADReal & temperature,
                        std::vector<FluidStateProperties> & fsp) const;
 
  ADReal liquidDensity(const ADReal & pressure, const ADReal & temperature) const;
 
  ADReal gasDensity(const ADReal & pressure,
                    const ADReal & temperature,
                    std::vector<FluidStateProperties> & fsp) const;
 
  ADReal saturation(const ADReal & pressure,
                    const ADReal & temperature,
                    const ADReal & Z,
                    std::vector<FluidStateProperties> & fsp) const;
 
  void twoPhaseProperties(const ADReal & pressure,
                          const ADReal & temperature,
                          const ADReal & Z,
                          unsigned int qp,
                          std::vector<FluidStateProperties> & fsp) const;
 
  ADReal enthalpyOfDissolution(const ADReal & temperature) const;
 
  virtual Real totalMassFraction(
      Real pressure, Real temperature, Real Xnacl, Real saturation, unsigned int qp) const override;
 
protected:
  ADReal moleFractionToMassFraction(const ADReal & 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;
};