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;
};