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 : public Material
{
public:
  static InputParameters validParams();

  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;

  const 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);
};