PorousFlowPorosity.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 "PorousFlowPorosityExponentialBase.h"
 
template <bool is_ad>
class PorousFlowPorosityTempl : public PorousFlowPorosityExponentialBaseTempl<is_ad>
{
public:
  static InputParameters validParams();
 
  PorousFlowPorosityTempl(const InputParameters & parameters);
 
protected:
  virtual GenericReal<is_ad> atNegInfinityQp() const override;
  virtual Real datNegInfinityQp(unsigned pvar) const override;
  virtual GenericReal<is_ad> atZeroQp() const override;
  virtual Real datZeroQp(unsigned pvar) const override;
  virtual GenericReal<is_ad> decayQp() const override;
  virtual Real ddecayQp_dvar(unsigned pvar) const override;
  virtual RealGradient ddecayQp_dgradvar(unsigned pvar) const override;
 
  const bool _mechanical;
 
  const bool _fluid;
 
  const bool _thermal;
 
  const bool _chemical;
 
  const VariableValue & _phi0;
 
  const Real _biot;
 
  const Real _exp_coeff;
 
  const Moose::Functor<Real> * _solid_bulk;
 
  const Real _coeff;
 
  const VariableValue & _t_reference;
 
  const VariableValue & _p_reference;
 
  const unsigned _num_c_ref;
 
  std::vector<const VariableValue *> _c_reference;
 
  const unsigned _num_initial_c;
 
  std::vector<const VariableValue *> _initial_c;
 
  std::vector<Real> _c_weights;
 
  const MaterialProperty<Real> * const _porosity_old;
 
  const GenericMaterialProperty<Real, is_ad> * const _vol_strain_qp;
 
  const MaterialProperty<std::vector<RealGradient>> * const _dvol_strain_qp_dvar;
 
  const GenericMaterialProperty<Real, is_ad> * const _pf;
 
  const MaterialProperty<std::vector<Real>> * const _dpf_dvar;
 
  const GenericMaterialProperty<Real, is_ad> * const _temperature;
 
  const MaterialProperty<std::vector<Real>> * const _dtemperature_dvar;
 
  const MaterialProperty<std::vector<Real>> * const _mineral_conc_old;
 
  const MaterialProperty<std::vector<Real>> * const _reaction_rate;
 
  const MaterialProperty<std::vector<std::vector<Real>>> * const _dreaction_rate_dvar;
 
  const unsigned int _aq_ph;
 
  const MaterialProperty<std::vector<Real>> * const _saturation;
 
  const MaterialProperty<std::vector<std::vector<Real>>> * const _dsaturation_dvar;
 
  usingPorousFlowPorosityExponentialBaseMembers;
};
 
typedef PorousFlowPorosityTempl<false> PorousFlowPorosity;
typedef PorousFlowPorosityTempl<true> ADPorousFlowPorosity;