RichardsMaterial.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 "Material.h"
 
#include "RichardsVarNames.h"
#include "RichardsDensity.h"
#include "RichardsRelPerm.h"
#include "RichardsSeff.h"
#include "RichardsSat.h"
#include "RichardsSUPG.h"
 
// Forward Declarations
class RichardsMaterial;
 
template <>
InputParameters validParams<RichardsMaterial>();
 
class RichardsMaterial : public Material
{
public:
  RichardsMaterial(const InputParameters & parameters);
 
protected:
  Real _material_por;
 
  const VariableValue & _por_change;
  const VariableValue & _por_change_old;
 
  RealTensorValue _material_perm;
 
  RealVectorValue _material_gravity;
 
  MaterialProperty<Real> & _porosity_old;
  MaterialProperty<Real> & _porosity;
  MaterialProperty<RealTensorValue> & _permeability;
  MaterialProperty<RealVectorValue> & _gravity;
 
  const RichardsVarNames & _richards_name_UO;
  unsigned int _num_p;
 
  std::vector<const RichardsRelPerm *> _material_relperm_UO;
  std::vector<const RichardsSeff *> _material_seff_UO;
  std::vector<const RichardsSat *> _material_sat_UO;
  std::vector<const RichardsDensity *> _material_density_UO;
  std::vector<const RichardsSUPG *> _material_SUPG_UO;
 
  std::vector<const VariableValue *> _perm_change;
 
  virtual void computeProperties();
 
  void computePandSeff();
 
  void computeDerivedQuantities(unsigned int qp);
 
  void compute2ndDerivedQuantities(unsigned int qp);
 
  void zeroSUPG(unsigned int qp);
 
  void computeSUPG();
 
private:
  Real _trace_perm;
 
  std::vector<Real> _material_viscosity;
 
  MaterialProperty<std::vector<Real>> & _pp_old;
 
  MaterialProperty<std::vector<Real>> & _pp;
 
  MaterialProperty<std::vector<std::vector<Real>>> & _dpp_dv;
 
  MaterialProperty<std::vector<std::vector<std::vector<Real>>>> & _d2pp_dv;
 
  MaterialProperty<std::vector<Real>> & _viscosity;
 
  MaterialProperty<std::vector<Real>> & _density_old;
 
  MaterialProperty<std::vector<Real>> & _density;
 
  MaterialProperty<std::vector<std::vector<Real>>> & _ddensity_dv;
 
  MaterialProperty<std::vector<Real>> & _seff_old;
 
  MaterialProperty<std::vector<Real>> & _seff; // effective saturation
 
  MaterialProperty<std::vector<std::vector<Real>>> & _dseff_dv; // d(seff)/dp
 
  MaterialProperty<std::vector<std::vector<std::vector<Real>>>> & _d2seff_dv;
 
  MaterialProperty<std::vector<Real>> & _sat_old;
 
  MaterialProperty<std::vector<Real>> & _sat;
 
  MaterialProperty<std::vector<std::vector<Real>>> & _dsat_dv;
 
  MaterialProperty<std::vector<Real>> & _rel_perm;
 
  MaterialProperty<std::vector<std::vector<Real>>> & _drel_perm_dv;
 
  MaterialProperty<std::vector<Real>> & _mass_old;
 
  MaterialProperty<std::vector<Real>> & _mass;
 
  MaterialProperty<std::vector<std::vector<Real>>> & _dmass;
 
  MaterialProperty<std::vector<RealVectorValue>> & _flux_no_mob;
 
  MaterialProperty<std::vector<std::vector<RealVectorValue>>> & _dflux_no_mob_dv;
 
  MaterialProperty<std::vector<std::vector<RealTensorValue>>> & _dflux_no_mob_dgradv;
 
  MaterialProperty<std::vector<RealVectorValue>> & _flux;
 
  MaterialProperty<std::vector<std::vector<RealVectorValue>>> & _dflux_dv;
 
  MaterialProperty<std::vector<std::vector<RealTensorValue>>> & _dflux_dgradv;
 
  MaterialProperty<std::vector<std::vector<std::vector<RealVectorValue>>>> & _d2flux_dvdv;
 
  MaterialProperty<std::vector<std::vector<std::vector<RealTensorValue>>>> & _d2flux_dgradvdv;
 
  MaterialProperty<std::vector<std::vector<std::vector<RealTensorValue>>>> & _d2flux_dvdgradv;
 
  MaterialProperty<std::vector<RealVectorValue>> & _tauvel_SUPG; // tauSUPG * velSUPG
  MaterialProperty<std::vector<std::vector<RealTensorValue>>> &
      _dtauvel_SUPG_dgradp; // d (_tauvel_SUPG_i)/d(_grad_variable_j)
  MaterialProperty<std::vector<std::vector<RealVectorValue>>> &
      _dtauvel_SUPG_dp; // d (_tauvel_SUPG_i)/d(variable_j)
 
  std::vector<std::vector<std::vector<Real>>> _d2density;
 
  std::vector<std::vector<std::vector<Real>>> _d2rel_perm_dv;
 
  std::vector<const VariableValue *> _pressure_vals;
  std::vector<const VariableValue *> _pressure_old_vals;
  std::vector<const VariableGradient *> _grad_p;
 
  void zero2ndDerivedQuantities(unsigned int qp);
};