RichardsPiecewiseLinearSink.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 "IntegratedBC.h"
#include "LinearInterpolation.h"
#include "Function.h"
#include "RichardsVarNames.h"
#include "RichardsDensity.h"
#include "RichardsRelPerm.h"
#include "RichardsSeff.h"
 
// Forward Declarations
class RichardsPiecewiseLinearSink;
 
template <>
InputParameters validParams<RichardsPiecewiseLinearSink>();
 
class RichardsPiecewiseLinearSink : public IntegratedBC
{
public:
  RichardsPiecewiseLinearSink(const InputParameters & parameters);
 
protected:
  virtual void computeResidual() override;
 
  virtual Real computeQpResidual() override;
 
  virtual void computeJacobian() override;
 
  virtual Real computeQpJacobian() override;
 
  virtual void computeJacobianBlock(MooseVariableFEBase & jvar) override;
  using IntegratedBC::computeJacobianBlock;
 
  virtual Real computeQpOffDiagJacobian(unsigned int jvar) override;
 
  bool _use_mobility;
 
  bool _use_relperm;
 
  bool _fully_upwind;
 
  LinearInterpolation _sink_func;
 
  const Function & _m_func;
 
  const RichardsVarNames & _richards_name_UO;
 
  unsigned int _num_p;
 
  unsigned int _pvar;
 
  const RichardsDensity * _density_UO;
 
  const RichardsSeff * _seff_UO;
 
  const RichardsRelPerm * _relperm_UO;
 
  const PostprocessorValue & _area_pp;
 
  unsigned int _num_nodes;
 
  std::vector<Real> _nodal_density;
 
  std::vector<std::vector<Real>> _dnodal_density_dv;
 
  std::vector<Real> _nodal_relperm;
 
  std::vector<std::vector<Real>> _dnodal_relperm_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);
};