Q2PPorepressureFlux.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 "Kernel.h"
#include "RichardsDensity.h"
#include "RichardsRelPerm.h"
#include "Material.h"
 
// Forward Declarations
class Q2PPorepressureFlux;
 
template <>
InputParameters validParams<Q2PPorepressureFlux>();
 
class Q2PPorepressureFlux : public Kernel
{
public:
  Q2PPorepressureFlux(const InputParameters & parameters);
 
protected:
  virtual Real computeQpResidual() override;
 
  virtual void computeResidual() override;
 
  virtual void computeOffDiagJacobian(MooseVariableFEBase & jvar) override;
  using Kernel::computeOffDiagJacobian;
 
  virtual void computeJacobian() override;
 
  Real computeQpJac(unsigned int dvar);
 
  void upwind(bool compute_res, bool compute_jac, unsigned int jvar);
 
  void prepareNodalValues();
 
  const RichardsDensity & _density;
 
  const VariableValue & _sat;
 
  unsigned int _sat_var;
 
  const RichardsRelPerm & _relperm;
 
  Real _viscosity;
 
  const MaterialProperty<RealVectorValue> & _gravity;
 
  const MaterialProperty<RealTensorValue> & _permeability;
 
  unsigned int _num_nodes;
 
  std::vector<Real> _mobility;
 
  std::vector<Real> _dmobility_dp;
 
  std::vector<Real> _dmobility_ds;
};