PorousFlowPorosity.h

Go to the documentation of this file.

//* 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 "PorousFlowPorosityExponentialBase.h"

class PorousFlowPorosity : public PorousFlowPorosityExponentialBase
{
public:
  static InputParameters validParams();

  PorousFlowPorosity(const InputParameters & parameters);

protected:
  virtual Real atNegInfinityQp() const override;
  virtual Real datNegInfinityQp(unsigned pvar) const override;
  virtual Real atZeroQp() const override;
  virtual Real datZeroQp(unsigned pvar) const override;
  virtual Real 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 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 MaterialProperty<Real> * const _vol_strain_qp;

  const MaterialProperty<std::vector<RealGradient>> * const _dvol_strain_qp_dvar;

  const MaterialProperty<Real> * const _pf;

  const MaterialProperty<std::vector<Real>> * const _dpf_dvar;

  const MaterialProperty<Real> * 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;
};