RichardsPolyLineSink.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 "DiracKernel.h"
#include "LinearInterpolation.h"
#include "RichardsSumQuantity.h"
#include "RichardsVarNames.h"
 
// Forward Declarations
class RichardsPolyLineSink;
 
template <>
InputParameters validParams<RichardsPolyLineSink>();
 
class RichardsPolyLineSink : public DiracKernel
{
public:
  RichardsPolyLineSink(const InputParameters & parameters);
 
  virtual void addPoints();
  virtual Real computeQpResidual();
  virtual Real computeQpJacobian();
 
  virtual Real computeQpOffDiagJacobian(unsigned int jvar);
 
protected:
  RichardsSumQuantity & _total_outflow_mass;
 
  LinearInterpolation _sink_func;
 
  std::string _point_file;
 
  const RichardsVarNames & _richards_name_UO;
 
  unsigned int _pvar;
 
  const MaterialProperty<std::vector<Real>> & _pp;
 
  const MaterialProperty<std::vector<std::vector<Real>>> & _dpp_dv;
 
  std::vector<Real> _xs;
 
  std::vector<Real> _ys;
 
  std::vector<Real> _zs;
 
  bool parseNextLineReals(std::ifstream & ifs, std::vector<Real> & myvec);
};