RichardsBorehole.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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

// Moose Includes
#include "PeacemanBorehole.h"
#include "RichardsVarNames.h"
#include "RichardsDensity.h"
#include "RichardsRelPerm.h"
#include "RichardsSeff.h"

class RichardsBorehole : public PeacemanBorehole
{
public:
  static InputParameters validParams();

  RichardsBorehole(const InputParameters & parameters);

  virtual void computeResidual();

  virtual Real computeQpResidual();

  virtual void computeJacobian();

  virtual Real computeQpJacobian();

  virtual Real computeQpOffDiagJacobian(unsigned int jvar);

protected:
  const bool _fully_upwind;

  const RichardsVarNames & _richards_name_UO;

  const unsigned int _num_p;

  const unsigned int _pvar;

  const RichardsDensity * const _density_UO;

  const RichardsSeff * const _seff_UO;

  const RichardsRelPerm * const _relperm_UO;

  unsigned int _num_nodes;

  std::vector<Real> _mobility;

  std::vector<std::vector<Real>> _dmobility_dv;

  const MaterialProperty<std::vector<Real>> & _pp;

  const MaterialProperty<std::vector<std::vector<Real>>> & _dpp_dv;

  const MaterialProperty<std::vector<Real>> & _viscosity;

  const MaterialProperty<RealTensorValue> & _permeability;

  const MaterialProperty<std::vector<std::vector<Real>>> & _dseff_dv;

  const MaterialProperty<std::vector<Real>> & _rel_perm;

  const MaterialProperty<std::vector<std::vector<Real>>> & _drel_perm_dv;

  const MaterialProperty<std::vector<Real>> & _density;

  const MaterialProperty<std::vector<std::vector<Real>>> & _ddensity_dv;

  std::vector<const VariableValue *> _ps_at_nodes;

  void prepareNodalValues();

  Real jac(unsigned int wrt_num);
};