StiffenedGasTwoPhaseFluidProperties.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 "TwoPhaseFluidProperties.h"
#include "LinearInterpolation.h"
#include "NaNInterface.h"

class StiffenedGasFluidProperties;

#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Woverloaded-virtual"

class StiffenedGasTwoPhaseFluidProperties : public TwoPhaseFluidProperties, public NaNInterface
{
public:
  StiffenedGasTwoPhaseFluidProperties(const InputParameters & parameters);

  virtual Real p_critical() const override;
  virtual Real T_triple() const override;
  virtual Real L_fusion() const override;
  virtual Real T_sat(Real pressure) const override;
  virtual Real p_sat(Real temperature) const override;
  virtual Real dT_sat_dp(Real pressure) const override;
  virtual Real sigma_from_T(Real T) const override;
  virtual Real dsigma_dT_from_T(Real T) const override;

  virtual bool supportsPhaseChange() const override { return true; }

protected:
  Real compute_p_sat(const Real & T) const;

  // liquid SGEOS parameters
  const Real _gamma_liquid;
  const Real _cv_liquid;
  const Real _cp_liquid;
  const Real _q_liquid;
  const Real _p_inf_liquid;
  const Real _q_prime_liquid;

  // vapor SGEOS parameters
  const Real _gamma_vapor;
  const Real _cv_vapor;
  const Real _cp_vapor;
  const Real _q_vapor;
  const Real _p_inf_vapor;
  const Real _q_prime_vapor;

  const Real & _T_c;
  const Real & _p_c;
  const Real & _T_triple;
  const Real & _L_fusion;

  const Real & _sigma_A;
  const Real & _sigma_B;
  const Real & _sigma_C;

  const Real & _T_sat_min;
  const Real & _T_sat_max;
  const Real & _p_sat_guess;
  const unsigned int & _n_sat_samples;
  const Real _dT_sat;

  // coefficients for saturation pressure Newton solve
  const Real _A;
  const Real _B;
  const Real _C;
  const Real _D;

  const Real _newton_tol;
  const unsigned int _newton_max_iter;

  // These two vectors store saturation line p(T) information by only calculating
  // once in constructor and then use interpolation to quickly calculate the value
  std::vector<Real> _T_vec;
  std::vector<Real> _p_sat_vec;

  LinearInterpolation _ipol_temp;
  LinearInterpolation _ipol_pressure;

public:
  static InputParameters validParams();
};

#pragma GCC diagnostic pop