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

  Q2PSaturationFlux(const InputParameters & parameters);

protected:
  virtual Real computeQpResidual() override;

  virtual void computeResidual() override;

  virtual void computeOffDiagJacobian(unsigned int jvar) override;

  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 & _pp;

  const VariableGradient & _grad_pp;

  const VariableValue & _pp_nodal;

  unsigned int _pp_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;
};