Approximates a line sink by a sequence of Dirac Points. More...

#include <PorousFlowPolyLineSink.h>

Inheritance diagram for PorousFlowPolyLineSink:

Public Member Functions
	PorousFlowPolyLineSink (const InputParameters &parameters)

Protected Types
enum	PorTchoice { PorTchoice::pressure, PorTchoice::temperature }
	whether the flux is a function of pressure or temperature More...

Protected Member Functions
Real	computeQpBaseOutflow (unsigned current_dirac_ptid) const override
	Returns the flux from the line sink (before modification by mobility, etc). Derived classes should override this. More...

void	computeQpBaseOutflowJacobian (unsigned jvar, unsigned current_dirac_ptid, Real &outflow, Real &outflowp) const override
	Calculates the BaseOutflow as well as its derivative wrt jvar. Derived classes should override this. More...

virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

virtual void	addPoints () override
	Add Dirac Points to the borehole. More...

Real	jac (unsigned int jvar)
	Jacobian contribution for the derivative wrt the variable jvar. More...

Real	ptqp () const
	If _p_or_t==0, then returns the quadpoint porepressure, else returns the quadpoint temperature. More...

Real	dptqp (unsigned pvar) const
	If _p_or_t==0, then returns d(quadpoint porepressure)/d(PorousFlow variable), else returns d(quadpoint temperature)/d(PorousFlow variable) More...

bool	parseNextLineReals (std::ifstream &ifs, std::vector< Real > &myvec)
	Reads a space-separated line of floats from ifs and puts in myvec. More...

Protected Attributes
LinearInterpolation	_sink_func
	mass flux = _sink_func as a function of porepressure or temperature More...

const PorousFlowDictator &	_dictator
	PorousFlowDictator UserObject. More...

PorousFlowSumQuantity &	_total_outflow_mass
	This is used to hold the total fluid flowing into the line sink for each time step. More...

const bool	_has_porepressure
	Whether a quadpoint porepressure material exists (for error checking) More...

const bool	_has_temperature
	Whether a quadpoint temperature material exists (for error checking) More...

const bool	_has_mass_fraction
	Whether a mass_fraction material exists (for error checking) More...

const bool	_has_relative_permeability
	Whether a relative permeability material exists (for error checking) More...

const bool	_has_mobility
	Whether enough materials exist to form the mobility (for error checking) More...

const bool	_has_enthalpy
	Whether an enthalpy material exists (for error checking) More...

const bool	_has_internal_energy
	Whether an internal-energy material exists (for error checking) More...

enum PorousFlowLineSink::PorTchoice	_p_or_t

const bool	_use_mass_fraction
	Whether the flux will be multiplied by the mass fraction. More...

const bool	_use_relative_permeability
	Whether the flux will be multiplied by the relative permeability. More...

const bool	_use_mobility
	Whether the flux will be multiplied by the mobility. More...

const bool	_use_enthalpy
	Whether the flux will be multiplied by the enthalpy. More...

const bool	_use_internal_energy
	Whether the flux will be multiplied by the internal-energy. More...

const unsigned int	_ph
	The phase number. More...

const unsigned int	_sp
	The component number (only used if _use_mass_fraction==true) More...

const MaterialProperty< std::vector< Real > > *const	_pp
	Quadpoint pore pressure in each phase. More...

const MaterialProperty< std::vector< std::vector< Real > > > *const	_dpp_dvar
	d(quadpoint pore pressure in each phase)/d(PorousFlow variable) More...

const MaterialProperty< Real > *const	_temperature
	Quadpoint temperature. More...

const MaterialProperty< std::vector< Real > > *const	_dtemperature_dvar
	d(quadpoint temperature)/d(PorousFlow variable) More...

const MaterialProperty< std::vector< Real > > *const	_fluid_density_node
	Fluid density for each phase (at the node) More...

const MaterialProperty< std::vector< std::vector< Real > > > *const	_dfluid_density_node_dvar
	d(Fluid density for each phase (at the node))/d(PorousFlow variable) More...

const MaterialProperty< std::vector< Real > > *const	_fluid_viscosity
	Viscosity of each component in each phase. More...

const MaterialProperty< std::vector< std::vector< Real > > > *const	_dfluid_viscosity_dvar
	d(Viscosity of each component in each phase)/d(PorousFlow variable) More...

const MaterialProperty< std::vector< Real > > *const	_relative_permeability
	Relative permeability of each phase. More...

const MaterialProperty< std::vector< std::vector< Real > > > *const	_drelative_permeability_dvar
	d(Relative permeability of each phase)/d(PorousFlow variable) More...

const MaterialProperty< std::vector< std::vector< Real > > > *const	_mass_fractions
	Mass fraction of each component in each phase. More...

const MaterialProperty< std::vector< std::vector< std::vector< Real > > > > *const	_dmass_fractions_dvar
	d(Mass fraction of each component in each phase)/d(PorousFlow variable) More...

const MaterialProperty< std::vector< Real > > *const	_enthalpy
	Enthalpy of each phase. More...

const MaterialProperty< std::vector< std::vector< Real > > > *const	_denthalpy_dvar
	d(enthalpy of each phase)/d(PorousFlow variable) More...

const MaterialProperty< std::vector< Real > > *const	_internal_energy
	Internal_Energy of each phase. More...

const MaterialProperty< std::vector< std::vector< Real > > > *const	_dinternal_energy_dvar
	d(internal_energy of each phase)/d(PorousFlow variable) More...

const Real	_line_length
	Line length. This is only used if there is only one borehole point. More...

const RealVectorValue	_line_direction
	Line direction. This is only used if there is only one borehole point. More...

const std::string	_point_file
	File defining the geometry of the borehole. More...

std::vector< Real >	_rs
	Radii of the borehole. More...

std::vector< Real >	_xs
	x points of the borehole More...

std::vector< Real >	_ys
	y points of the borehole More...

std::vector< Real >	_zs
	z points of borehole More...

Point	_bottom_point
	The bottom point of the borehole (where bottom_pressure is defined) More...

std::vector< Real >	_half_seg_len
	0.5*(length of polyline segments between points) More...

Detailed Description

Approximates a line sink by a sequence of Dirac Points.

Definition at line 23 of file PorousFlowPolyLineSink.h.

Member Enumeration Documentation

◆ PorTchoice

enum PorousFlowLineSink::PorTchoice

strongprotectedinherited

whether the flux is a function of pressure or temperature

Enumerator
pressure
temperature

Definition at line 91 of file PorousFlowLineSink.h.

91 { pressure, temperature } _p_or_t;

Constructor & Destructor Documentation

◆ PorousFlowPolyLineSink()

PorousFlowPolyLineSink::PorousFlowPolyLineSink ( const InputParameters & parameters )

Definition at line 41 of file PorousFlowPolyLineSink.C.

   : PorousFlowLineSink(parameters),
     _sink_func(getParam<std::vector<Real>>("p_or_t_vals"), getParam<std::vector<Real>>("fluxes"))
 {
 }

Member Function Documentation

◆ addPoints()

void PorousFlowLineSink::addPoints ( )

overrideprotectedvirtualinherited

Add Dirac Points to the borehole.

Reimplemented from PorousFlowLineGeometry.

Definition at line 217 of file PorousFlowLineSink.C.

 {
   // This function gets called just before the DiracKernel is evaluated
   // so this is a handy place to zero this out.
   _total_outflow_mass.zero();
  
   PorousFlowLineGeometry::addPoints();
 }

◆ computeQpBaseOutflow()

Real PorousFlowPolyLineSink::computeQpBaseOutflow ( unsigned current_dirac_ptid ) const

overrideprotectedvirtual

Returns the flux from the line sink (before modification by mobility, etc). Derived classes should override this.

Implements PorousFlowLineSink.

Definition at line 48 of file PorousFlowPolyLineSink.C.

 {
   Real outflow = 0.0;
  
   if (current_dirac_ptid > 0)
     // contribution from half-segment "behind" this point (must have >1 point for
     // current_dirac_ptid>0)
     outflow += _half_seg_len[current_dirac_ptid - 1];
  
   if (current_dirac_ptid + 1 < _zs.size() || _zs.size() == 1)
     // contribution from half-segment "ahead of" this point, or we only have one point
     outflow += _half_seg_len[current_dirac_ptid];
  
   return outflow * _test[_i][_qp] * _sink_func.sample(ptqp()) * _rs[current_dirac_ptid];
 }

◆ computeQpBaseOutflowJacobian()

void PorousFlowPolyLineSink::computeQpBaseOutflowJacobian	(	unsigned	jvar,
		unsigned	current_dirac_ptid,
		Real &	outflow,
		Real &	outflowp
	)		const

overrideprotectedvirtual

Calculates the BaseOutflow as well as its derivative wrt jvar. Derived classes should override this.

Implements PorousFlowLineSink.

Definition at line 65 of file PorousFlowPolyLineSink.C.

 {
   outflow = 0.0;
   outflowp = 0.0;
  
   if (_dictator.notPorousFlowVariable(jvar))
     return;
   const unsigned pvar = _dictator.porousFlowVariableNum(jvar);
  
   if (current_dirac_ptid > 0)
     // contribution from half-segment "behind" this point (must have >1 point for
     // current_dirac_ptid>0)
     outflow += _half_seg_len[current_dirac_ptid - 1];
  
   if (current_dirac_ptid + 1 < _zs.size() || _zs.size() == 1)
     // contribution from half-segment "ahead of" this point, or we only have one point
     outflow += _half_seg_len[current_dirac_ptid];
  
   outflowp = outflow * _test[_i][_qp] * _sink_func.sampleDerivative(ptqp()) * dptqp(pvar) *
              _phi[_j][_qp] * _rs[current_dirac_ptid];
   outflow *= _test[_i][_qp] * _sink_func.sample(ptqp()) * _rs[current_dirac_ptid];
 }

◆ computeQpJacobian()

Real PorousFlowLineSink::computeQpJacobian ( )

overrideprotectedvirtualinherited

Definition at line 258 of file PorousFlowLineSink.C.

 {
   return jac(_var.number());
 }

◆ computeQpOffDiagJacobian()

Real PorousFlowLineSink::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtualinherited

Definition at line 264 of file PorousFlowLineSink.C.

 {
   return jac(jvar);
 }

◆ computeQpResidual()

Real PorousFlowLineSink::computeQpResidual ( )

overrideprotectedvirtualinherited

Definition at line 227 of file PorousFlowLineSink.C.

 {
   // Get the ID we initially assigned to this point
   const unsigned current_dirac_ptid = currentPointCachedID();
   Real outflow = computeQpBaseOutflow(current_dirac_ptid);
   if (outflow == 0.0)
     return 0.0;
  
   if (_use_relative_permeability)
     outflow *= (*_relative_permeability)[_i][_ph];
  
   if (_use_mobility)
     outflow *= (*_relative_permeability)[_i][_ph] * (*_fluid_density_node)[_i][_ph] /
                (*_fluid_viscosity)[_i][_ph];
  
   if (_use_mass_fraction)
     outflow *= (*_mass_fractions)[_i][_ph][_sp];
  
   if (_use_enthalpy)
     outflow *= (*_enthalpy)[_i][_ph];
  
   if (_use_internal_energy)
     outflow *= (*_internal_energy)[_i][_ph];
  
   _total_outflow_mass.add(
       outflow * _dt); // this is not thread safe, but DiracKernel's aren't currently threaded
  
   return outflow;
 }

◆ dptqp()

Real PorousFlowLineSink::dptqp ( unsigned pvar ) const

protectedinherited

If _p_or_t==0, then returns d(quadpoint porepressure)/d(PorousFlow variable), else returns d(quadpoint temperature)/d(PorousFlow variable)

Parameters

pvar	The PorousFlow variable number

Definition at line 341 of file PorousFlowLineSink.C.

 {
   return (_p_or_t == PorTchoice::pressure ? (*_dpp_dvar)[_qp][_ph][pvar]
                                           : (*_dtemperature_dvar)[_qp][pvar]);
 }

Referenced by computeQpBaseOutflowJacobian(), and PorousFlowPeacemanBorehole::computeQpBaseOutflowJacobian().

◆ jac()

Real PorousFlowLineSink::jac ( unsigned int jvar )

protectedinherited

Jacobian contribution for the derivative wrt the variable jvar.

Definition at line 270 of file PorousFlowLineSink.C.

 {
   if (_dictator.notPorousFlowVariable(jvar))
     return 0.0;
   const unsigned pvar = _dictator.porousFlowVariableNum(jvar);
  
   Real outflow;
   Real outflowp;
   const unsigned current_dirac_ptid = currentPointCachedID();
   computeQpBaseOutflowJacobian(jvar, current_dirac_ptid, outflow, outflowp);
   if (outflow == 0.0 && outflowp == 0.0)
     return 0.0;
  
   if (_use_relative_permeability)
   {
     const Real relperm_prime = (_i != _j ? 0.0 : (*_drelative_permeability_dvar)[_i][_ph][pvar]);
     outflowp = (*_relative_permeability)[_i][_ph] * outflowp + relperm_prime * outflow;
     outflow *= (*_relative_permeability)[_i][_ph];
   }
  
   if (_use_mobility)
   {
     const Real mob = (*_relative_permeability)[_i][_ph] * (*_fluid_density_node)[_i][_ph] /
                      (*_fluid_viscosity)[_i][_ph];
     const Real mob_prime =
         (_i != _j
              ? 0.0
              : (*_drelative_permeability_dvar)[_i][_ph][pvar] * (*_fluid_density_node)[_i][_ph] /
                        (*_fluid_viscosity)[_i][_ph] +
                    (*_relative_permeability)[_i][_ph] *
                        (*_dfluid_density_node_dvar)[_i][_ph][pvar] / (*_fluid_viscosity)[_i][_ph] -
                    (*_relative_permeability)[_i][_ph] * (*_fluid_density_node)[_i][_ph] *
                        (*_dfluid_viscosity_dvar)[_i][_ph][pvar] /
                        Utility::pow<2>((*_fluid_viscosity)[_i][_ph]));
     outflowp = mob * outflowp + mob_prime * outflow;
     outflow *= mob;
   }
  
   if (_use_mass_fraction)
   {
     const Real mass_fractions_prime =
         (_i != _j ? 0.0 : (*_dmass_fractions_dvar)[_i][_ph][_sp][pvar]);
     outflowp = (*_mass_fractions)[_i][_ph][_sp] * outflowp + mass_fractions_prime * outflow;
     outflow *= (*_mass_fractions)[_i][_ph][_sp];
   }
  
   if (_use_enthalpy)
   {
     const Real enthalpy_prime = (_i != _j ? 0.0 : (*_denthalpy_dvar)[_i][_ph][pvar]);
     outflowp = (*_enthalpy)[_i][_ph] * outflowp + enthalpy_prime * outflow;
     outflow *= (*_enthalpy)[_i][_ph];
   }
  
   if (_use_internal_energy)
   {
     const Real internal_energy_prime = (_i != _j ? 0.0 : (*_dinternal_energy_dvar)[_i][_ph][pvar]);
     outflowp = (*_internal_energy)[_i][_ph] * outflowp + internal_energy_prime * outflow;
     // this multiplication was performed, but the code does not need to know: outflow *=
     // (*_internal_energy)[_i][_ph];
   }
  
   return outflowp;
 }

Referenced by PorousFlowLineSink::computeQpJacobian(), and PorousFlowLineSink::computeQpOffDiagJacobian().

◆ parseNextLineReals()

bool PorousFlowLineGeometry::parseNextLineReals	(	std::ifstream &	ifs,
		std::vector< Real > &	myvec
	)

protectedinherited

Reads a space-separated line of floats from ifs and puts in myvec.

Definition at line 103 of file PorousFlowLineGeometry.C.

 {
   std::string line;
   myvec.clear();
   bool gotline(false);
   if (getline(ifs, line))
   {
     gotline = true;
  
     // Harvest floats separated by whitespace
     std::istringstream iss(line);
     Real f;
     while (iss >> f)
     {
       myvec.push_back(f);
     }
   }
   return gotline;
 }

Referenced by PorousFlowLineGeometry::PorousFlowLineGeometry().

◆ ptqp()

Real PorousFlowLineSink::ptqp ( ) const

protectedinherited

If _p_or_t==0, then returns the quadpoint porepressure, else returns the quadpoint temperature.

Definition at line 335 of file PorousFlowLineSink.C.

 {
   return (_p_or_t == PorTchoice::pressure ? (*_pp)[_qp][_ph] : (*_temperature)[_qp]);
 }

Referenced by computeQpBaseOutflow(), PorousFlowPeacemanBorehole::computeQpBaseOutflow(), computeQpBaseOutflowJacobian(), and PorousFlowPeacemanBorehole::computeQpBaseOutflowJacobian().

Member Data Documentation

◆ _bottom_point

Point PorousFlowLineGeometry::_bottom_point

protectedinherited

The bottom point of the borehole (where bottom_pressure is defined)

Definition at line 61 of file PorousFlowLineGeometry.h.

Referenced by PorousFlowPeacemanBorehole::computeQpBaseOutflow(), PorousFlowPeacemanBorehole::computeQpBaseOutflowJacobian(), and PorousFlowLineGeometry::PorousFlowLineGeometry().

◆ _denthalpy_dvar

const MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowLineSink::_denthalpy_dvar

protectedinherited

d(enthalpy of each phase)/d(PorousFlow variable)

Definition at line 154 of file PorousFlowLineSink.h.

◆ _dfluid_density_node_dvar

const MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowLineSink::_dfluid_density_node_dvar

protectedinherited

d(Fluid density for each phase (at the node))/d(PorousFlow variable)

Definition at line 130 of file PorousFlowLineSink.h.

◆ _dfluid_viscosity_dvar

const MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowLineSink::_dfluid_viscosity_dvar

protectedinherited

d(Viscosity of each component in each phase)/d(PorousFlow variable)

Definition at line 136 of file PorousFlowLineSink.h.

◆ _dictator

const PorousFlowDictator& PorousFlowLineSink::_dictator

protectedinherited

PorousFlowDictator UserObject.

Definition at line 60 of file PorousFlowLineSink.h.

Referenced by computeQpBaseOutflowJacobian(), PorousFlowPeacemanBorehole::computeQpBaseOutflowJacobian(), PorousFlowLineSink::jac(), and PorousFlowLineSink::PorousFlowLineSink().

◆ _dinternal_energy_dvar

const MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowLineSink::_dinternal_energy_dvar

protectedinherited

d(internal_energy of each phase)/d(PorousFlow variable)

Definition at line 160 of file PorousFlowLineSink.h.

◆ _dmass_fractions_dvar

const MaterialProperty<std::vector<std::vector<std::vector<Real> > > >* const PorousFlowLineSink::_dmass_fractions_dvar

protectedinherited

d(Mass fraction of each component in each phase)/d(PorousFlow variable)

Definition at line 148 of file PorousFlowLineSink.h.

◆ _dpp_dvar

const MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowLineSink::_dpp_dvar

protectedinherited

d(quadpoint pore pressure in each phase)/d(PorousFlow variable)

Definition at line 118 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::dptqp().

◆ _drelative_permeability_dvar

const MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowLineSink::_drelative_permeability_dvar

protectedinherited

d(Relative permeability of each phase)/d(PorousFlow variable)

Definition at line 142 of file PorousFlowLineSink.h.

◆ _dtemperature_dvar

const MaterialProperty<std::vector<Real> >* const PorousFlowLineSink::_dtemperature_dvar

protectedinherited

d(quadpoint temperature)/d(PorousFlow variable)

Definition at line 124 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::dptqp().

◆ _enthalpy

const MaterialProperty<std::vector<Real> >* const PorousFlowLineSink::_enthalpy

protectedinherited

Enthalpy of each phase.

Definition at line 151 of file PorousFlowLineSink.h.

◆ _fluid_density_node

const MaterialProperty<std::vector<Real> >* const PorousFlowLineSink::_fluid_density_node

protectedinherited

Fluid density for each phase (at the node)

Definition at line 127 of file PorousFlowLineSink.h.

◆ _fluid_viscosity

const MaterialProperty<std::vector<Real> >* const PorousFlowLineSink::_fluid_viscosity

protectedinherited

Viscosity of each component in each phase.

Definition at line 133 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::jac().

◆ _half_seg_len

std::vector<Real> PorousFlowLineGeometry::_half_seg_len

protectedinherited

0.5*(length of polyline segments between points)

Definition at line 64 of file PorousFlowLineGeometry.h.

Referenced by computeQpBaseOutflow(), PorousFlowPeacemanBorehole::computeQpBaseOutflow(), computeQpBaseOutflowJacobian(), PorousFlowPeacemanBorehole::computeQpBaseOutflowJacobian(), and PorousFlowLineGeometry::PorousFlowLineGeometry().

◆ _has_enthalpy

const bool PorousFlowLineSink::_has_enthalpy

protectedinherited

Whether an enthalpy material exists (for error checking)

Definition at line 85 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::PorousFlowLineSink().

◆ _has_internal_energy

const bool PorousFlowLineSink::_has_internal_energy

protectedinherited

Whether an internal-energy material exists (for error checking)

Definition at line 88 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::PorousFlowLineSink().

◆ _has_mass_fraction

const bool PorousFlowLineSink::_has_mass_fraction

protectedinherited

Whether a mass_fraction material exists (for error checking)

Definition at line 76 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::PorousFlowLineSink().

◆ _has_mobility

const bool PorousFlowLineSink::_has_mobility

protectedinherited

Whether enough materials exist to form the mobility (for error checking)

Definition at line 82 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::PorousFlowLineSink().

◆ _has_porepressure

const bool PorousFlowLineSink::_has_porepressure

protectedinherited

Whether a quadpoint porepressure material exists (for error checking)

Definition at line 70 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::PorousFlowLineSink().

◆ _has_relative_permeability

const bool PorousFlowLineSink::_has_relative_permeability

protectedinherited

Whether a relative permeability material exists (for error checking)

Definition at line 79 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::PorousFlowLineSink().

◆ _has_temperature

const bool PorousFlowLineSink::_has_temperature

protectedinherited

Whether a quadpoint temperature material exists (for error checking)

Definition at line 73 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::PorousFlowLineSink().

◆ _internal_energy

const MaterialProperty<std::vector<Real> >* const PorousFlowLineSink::_internal_energy

protectedinherited

Internal_Energy of each phase.

Definition at line 157 of file PorousFlowLineSink.h.

◆ _line_direction

const RealVectorValue PorousFlowLineGeometry::_line_direction

protectedinherited

Line direction. This is only used if there is only one borehole point.

Definition at line 39 of file PorousFlowLineGeometry.h.

Referenced by PorousFlowPeacemanBorehole::PorousFlowPeacemanBorehole().

◆ _line_length

const Real PorousFlowLineGeometry::_line_length

protectedinherited

Line length. This is only used if there is only one borehole point.

Definition at line 36 of file PorousFlowLineGeometry.h.

Referenced by PorousFlowLineGeometry::PorousFlowLineGeometry().

◆ _mass_fractions

const MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowLineSink::_mass_fractions

protectedinherited

Mass fraction of each component in each phase.

Definition at line 145 of file PorousFlowLineSink.h.

◆ _p_or_t

enum PorousFlowLineSink::PorTchoice PorousFlowLineSink::_p_or_t

protectedinherited

Referenced by PorousFlowLineSink::dptqp(), PorousFlowLineSink::PorousFlowLineSink(), PorousFlowPeacemanBorehole::PorousFlowPeacemanBorehole(), and PorousFlowLineSink::ptqp().

◆ _ph

const unsigned int PorousFlowLineSink::_ph

protectedinherited

The phase number.

Definition at line 109 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::computeQpResidual(), PorousFlowLineSink::dptqp(), PorousFlowLineSink::jac(), PorousFlowLineSink::PorousFlowLineSink(), and PorousFlowLineSink::ptqp().

◆ _point_file

const std::string PorousFlowLineGeometry::_point_file

protectedinherited

File defining the geometry of the borehole.

Each row has format weight x y z and the list of such points defines a polyline that is the line sink

Definition at line 46 of file PorousFlowLineGeometry.h.

Referenced by PorousFlowLineGeometry::PorousFlowLineGeometry().

◆ _pp

const MaterialProperty<std::vector<Real> >* const PorousFlowLineSink::_pp

protectedinherited

Quadpoint pore pressure in each phase.

Definition at line 115 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::ptqp().

◆ _relative_permeability

const MaterialProperty<std::vector<Real> >* const PorousFlowLineSink::_relative_permeability

protectedinherited

Relative permeability of each phase.

Definition at line 139 of file PorousFlowLineSink.h.

◆ _rs

std::vector<Real> PorousFlowLineGeometry::_rs

protectedinherited

Radii of the borehole.

Definition at line 49 of file PorousFlowLineGeometry.h.

Referenced by computeQpBaseOutflow(), PorousFlowPeacemanBorehole::computeQpBaseOutflow(), computeQpBaseOutflowJacobian(), PorousFlowPeacemanBorehole::computeQpBaseOutflowJacobian(), and PorousFlowLineGeometry::PorousFlowLineGeometry().

◆ _sink_func

LinearInterpolation PorousFlowPolyLineSink::_sink_func

protected

mass flux = _sink_func as a function of porepressure or temperature

Definition at line 30 of file PorousFlowPolyLineSink.h.

Referenced by computeQpBaseOutflow(), and computeQpBaseOutflowJacobian().

◆ _sp

const unsigned int PorousFlowLineSink::_sp

protectedinherited

The component number (only used if _use_mass_fraction==true)

Definition at line 112 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::computeQpResidual(), PorousFlowLineSink::jac(), and PorousFlowLineSink::PorousFlowLineSink().

◆ _temperature

const MaterialProperty<Real>* const PorousFlowLineSink::_temperature

protectedinherited

Quadpoint temperature.

Definition at line 121 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::ptqp().

◆ _total_outflow_mass

PorousFlowSumQuantity& PorousFlowLineSink::_total_outflow_mass

protectedinherited

This is used to hold the total fluid flowing into the line sink for each time step.

Hence, it is positive for production wells where fluid is flowing from porespace into the line sink (and hence removed from the model)

Definition at line 67 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::addPoints(), PorousFlowLineSink::computeQpResidual(), and PorousFlowLineSink::PorousFlowLineSink().

◆ _use_enthalpy

const bool PorousFlowLineSink::_use_enthalpy

protectedinherited

Whether the flux will be multiplied by the enthalpy.

Definition at line 103 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::computeQpResidual(), PorousFlowLineSink::jac(), and PorousFlowLineSink::PorousFlowLineSink().

◆ _use_internal_energy

const bool PorousFlowLineSink::_use_internal_energy

protectedinherited

Whether the flux will be multiplied by the internal-energy.

Definition at line 106 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::computeQpResidual(), PorousFlowLineSink::jac(), and PorousFlowLineSink::PorousFlowLineSink().

◆ _use_mass_fraction

const bool PorousFlowLineSink::_use_mass_fraction

protectedinherited

Whether the flux will be multiplied by the mass fraction.

Definition at line 94 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::computeQpResidual(), PorousFlowLineSink::jac(), and PorousFlowLineSink::PorousFlowLineSink().

◆ _use_mobility

const bool PorousFlowLineSink::_use_mobility

protectedinherited

Whether the flux will be multiplied by the mobility.

Definition at line 100 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::computeQpResidual(), PorousFlowLineSink::jac(), and PorousFlowLineSink::PorousFlowLineSink().

◆ _use_relative_permeability

const bool PorousFlowLineSink::_use_relative_permeability

protectedinherited

Whether the flux will be multiplied by the relative permeability.

Definition at line 97 of file PorousFlowLineSink.h.

Referenced by PorousFlowLineSink::computeQpResidual(), PorousFlowLineSink::jac(), and PorousFlowLineSink::PorousFlowLineSink().

◆ _xs

std::vector<Real> PorousFlowLineGeometry::_xs

protectedinherited

x points of the borehole

Definition at line 52 of file PorousFlowLineGeometry.h.

Referenced by PorousFlowLineGeometry::addPoints(), PorousFlowLineGeometry::PorousFlowLineGeometry(), and PorousFlowPeacemanBorehole::PorousFlowPeacemanBorehole().

◆ _ys

std::vector<Real> PorousFlowLineGeometry::_ys

protectedinherited

y points of the borehole

Definition at line 55 of file PorousFlowLineGeometry.h.

Referenced by PorousFlowLineGeometry::addPoints(), PorousFlowLineGeometry::PorousFlowLineGeometry(), and PorousFlowPeacemanBorehole::PorousFlowPeacemanBorehole().

◆ _zs

std::vector<Real> PorousFlowLineGeometry::_zs

protectedinherited

z points of borehole

Definition at line 58 of file PorousFlowLineGeometry.h.

Referenced by PorousFlowLineGeometry::addPoints(), computeQpBaseOutflow(), PorousFlowPeacemanBorehole::computeQpBaseOutflow(), computeQpBaseOutflowJacobian(), PorousFlowPeacemanBorehole::computeQpBaseOutflowJacobian(), PorousFlowLineGeometry::PorousFlowLineGeometry(), and PorousFlowPeacemanBorehole::PorousFlowPeacemanBorehole().

The documentation for this class was generated from the following files:

porous_flow/include/dirackernels/PorousFlowPolyLineSink.h
porous_flow/src/dirackernels/PorousFlowPolyLineSink.C

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Detailed Description

Member Enumeration Documentation

◆ PorTchoice

Constructor & Destructor Documentation

◆ PorousFlowPolyLineSink()

Member Function Documentation

◆ addPoints()

◆ computeQpBaseOutflow()

◆ computeQpBaseOutflowJacobian()

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

◆ dptqp()

◆ jac()

◆ parseNextLineReals()

◆ ptqp()

Member Data Documentation

◆ _bottom_point

◆ _denthalpy_dvar

◆ _dfluid_density_node_dvar

◆ _dfluid_viscosity_dvar

◆ _dictator

◆ _dinternal_energy_dvar

◆ _dmass_fractions_dvar

◆ _dpp_dvar

◆ _drelative_permeability_dvar

◆ _dtemperature_dvar

◆ _enthalpy

◆ _fluid_density_node

◆ _fluid_viscosity

◆ _half_seg_len

◆ _has_enthalpy

◆ _has_internal_energy

◆ _has_mass_fraction

◆ _has_mobility

◆ _has_porepressure

◆ _has_relative_permeability

◆ _has_temperature

◆ _internal_energy

◆ _line_direction

◆ _line_length

◆ _mass_fractions

◆ _p_or_t

◆ _ph

◆ _point_file

◆ _pp

◆ _relative_permeability

◆ _rs

◆ _sink_func

◆ _sp

◆ _temperature

◆ _total_outflow_mass

◆ _use_enthalpy

◆ _use_internal_energy

◆ _use_mass_fraction

◆ _use_mobility

◆ _use_relative_permeability

◆ _xs

◆ _ys

◆ _zs