doxygen/modules/PorousFlowBrineCO2_8h_source.html

//* 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 BrineFluidProperties;

class SinglePhaseFluidProperties;

class Water97FluidProperties;

class PorousFlowBrineCO2;


template <>

InputParameters validParams<PorousFlowBrineCO2>();


class PorousFlowBrineCO2 : public PorousFlowFluidStateMultiComponentBase

{

public:

  PorousFlowBrineCO2(const InputParameters & parameters);


  virtual std::string fluidStateName() const override;


  virtual void thermophysicalProperties(Real pressure,

                                        Real temperature,

                                        Real Xnacl,

                                        Real Z,

                                        unsigned int qp,

                                        std::vector<FluidStateProperties> & fsp) const override;


  void equilibriumMoleFractions(const DualReal & pressure,

                                const DualReal & temperature,

                                const DualReal & Xnacl,

                                DualReal & xco2,

                                DualReal & yh2o) const;


  void equilibriumMassFractions(const DualReal & pressure,

                                const DualReal & temperature,

                                const DualReal & Xnacl,

                                DualReal & Xco2,

                                DualReal & Yh2o) const;


  void massFractions(const DualReal & pressure,

                     const DualReal & temperature,

                     const DualReal & Xnacl,

                     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,

                        const DualReal & Xnacl,

                        std::vector<FluidStateProperties> & fsp) const;


  DualReal saturation(const DualReal & pressure,

                      const DualReal & temperature,

                      const DualReal & Xnacl,

                      const DualReal & Z,

                      std::vector<FluidStateProperties> & fsp) const;


  void twoPhaseProperties(const DualReal & pressure,

                          const DualReal & temperature,

                          const DualReal & Xnacl,

                          const DualReal & Z,

                          unsigned int qp,

                          std::vector<FluidStateProperties> & fsp) const;


  void fugacityCoefficientsLowTemp(const DualReal & pressure,

                                   const DualReal & temperature,

                                   const DualReal & co2_density,

                                   DualReal & fco2,

                                   DualReal & fh2o) const;


  void fugacityCoefficientsHighTemp(const DualReal & pressure,

                                    const DualReal & temperature,

                                    const DualReal & co2_density,

                                    const DualReal & xco2,

                                    const DualReal & yh2o,

                                    DualReal & fco2,

                                    DualReal & fh2o) const;


  DualReal fugacityCoefficientH2OHighTemp(const DualReal & pressure,

                                          const DualReal & temperature,

                                          const DualReal & co2_density,

                                          const DualReal & xco2,

                                          const DualReal & yh2o) const;


  DualReal fugacityCoefficientCO2HighTemp(const DualReal & pressure,

                                          const DualReal & temperature,

                                          const DualReal & co2_density,

                                          const DualReal & xco2,

                                          const DualReal & yh2o) const;


  DualReal activityCoefficientH2O(const DualReal & temperature, const DualReal & xco2) const;


  DualReal activityCoefficientCO2(const DualReal & temperature, const DualReal & xco2) const;


  DualReal activityCoefficient(const DualReal & pressure,

                               const DualReal & temperature,

                               const DualReal & Xnacl) const;


  DualReal activityCoefficientHighTemp(const DualReal & temperature, const DualReal & Xnacl) const;


  DualReal equilibriumConstantH2O(const DualReal & temperature) const;


  DualReal equilibriumConstantCO2(const DualReal & temperature) const;


  DualReal partialDensityCO2(const DualReal & temperature) const;


  virtual Real totalMassFraction(

      Real pressure, Real temperature, Real Xnacl, Real saturation, unsigned int qp) const override;


  unsigned int saltComponentIndex() const { return _salt_component; };


  DualReal henryConstant(const DualReal & temperature, const DualReal & Xnacl) const;


  DualReal enthalpyOfDissolutionGas(const DualReal & temperature, const DualReal & Xnacl) const;


  DualReal enthalpyOfDissolution(const DualReal & temperature) const;


  void equilibriumMoleFractionsLowTemp(const DualReal & pressure,

                                       const DualReal & temperature,

                                       const DualReal & Xnacl,

                                       DualReal & xco2,

                                       DualReal & yh2o) const;


  void solveEquilibriumMoleFractionHighTemp(Real pressure,

                                            Real temperature,

                                            Real Xnacl,

                                            Real co2_density,

                                            Real & xco2,

                                            Real & yh2o) const;


protected:

  void checkVariables(Real pressure, Real temperature) const;


  void smoothCubicInterpolation(

      Real temperature, Real f0, Real df0, Real f1, Real df1, Real & value, Real & deriv) const;


  void funcABHighTemp(Real pressure,

                      Real temperature,

                      Real Xnacl,

                      Real co2_density,

                      Real xco2,

                      Real yh2o,

                      Real & A,

                      Real & B) const;


  void funcABHighTemp(Real pressure,

                      Real temperature,

                      Real Xnacl,

                      Real co2_density,

                      Real xco2,

                      Real yh2o,

                      Real & A,

                      Real & dA_dp,

                      Real & dA_dT,

                      Real & B,

                      Real & dB_dp,

                      Real & dB_dT,

                      Real & dB_dX) const;


  void funcABLowTemp(const DualReal & pressure,

                     const DualReal & temperature,

                     const DualReal & co2_density,

                     DualReal & A,

                     DualReal & B) const;


  const unsigned int _salt_component;

  const BrineFluidProperties & _brine_fp;

  const SinglePhaseFluidProperties & _co2_fp;

  const Real _Mh2o;

  const Real _invMh2o;

  const Real _Mco2;

  const Real _Mnacl;

  const Real _Rbar;

  const Real _Tlower;

  const Real _Tupper;

  const Real _Zmin;

  const std::vector<Real> _co2_henry;

};