www.mooseframework.org
RichardsPolyLineSink.h
Go to the documentation of this file.
1 //* This file is part of the MOOSE framework
2 //* https://www.mooseframework.org
3 //*
4 //* All rights reserved, see COPYRIGHT for full restrictions
5 //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
6 //*
7 //* Licensed under LGPL 2.1, please see LICENSE for details
8 //* https://www.gnu.org/licenses/lgpl-2.1.html
9 
10 #pragma once
11 
12 #include "DiracKernel.h"
13 #include "LinearInterpolation.h"
14 #include "RichardsSumQuantity.h"
15 #include "RichardsVarNames.h"
16 
17 // Forward Declarations
18 
25 {
26 public:
28 
30 
31  virtual void addPoints();
32  virtual Real computeQpResidual();
33  virtual Real computeQpJacobian();
34 
41  virtual Real computeQpOffDiagJacobian(unsigned int jvar);
42 
43 protected:
50 
53 
55  std::string _point_file;
56 
59 
61  unsigned int _pvar;
62 
65 
68 
70  std::vector<Real> _xs;
71 
73  std::vector<Real> _ys;
74 
76  std::vector<Real> _zs;
77 
83  bool parseNextLineReals(std::ifstream & ifs, std::vector<Real> & myvec);
84 };
Sums into _total This is used, for instance, to record the total mass flowing into a borehole...
Approximates a polyline by a sequence of Dirac Points the mass flux from each Dirac Point is _sink_fu...
virtual Real computeQpResidual()
std::string _point_file
contains rows of the form x y z (space separated)
static InputParameters validParams()
std::vector< Real > _ys
vector of Dirac Points&#39; y positions
This holds maps between pressure_var or pressure_var, sat_var used in RichardsMaterial and kernels...
virtual Real computeQpOffDiagJacobian(unsigned int jvar)
Computes the off-diagonal part of the jacobian Note: at March2014 this is never called since moose do...
RichardsPolyLineSink(const InputParameters &parameters)
LinearInterpolation _sink_func
mass flux = _sink_func as a function of porepressure
virtual Real computeQpJacobian()
unsigned int _pvar
The moose internal variable number of the richards variable of this Dirac Kernel. ...
bool parseNextLineReals(std::ifstream &ifs, std::vector< Real > &myvec)
reads a space-separated line of floats from ifs and puts in myvec
RichardsSumQuantity & _total_outflow_mass
This is used to hold the total fluid flowing into the sink Hence, it is positive for sinks where flui...
const MaterialProperty< std::vector< Real > > & _pp
fluid porepressure (or porepressures in case of multiphase)
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const MaterialProperty< std::vector< std::vector< Real > > > & _dpp_dv
d(porepressure_i)/d(variable_j)
std::vector< Real > _xs
vector of Dirac Points&#39; x positions
std::vector< Real > _zs
vector of Dirac Points&#39; z positions
const InputParameters & parameters() const
const RichardsVarNames & _richards_name_UO
Defines the richards variables in the simulation.