Approximates a borehole by a sequence of Dirac Points. More...

#include <PorousFlowPeacemanBorehole.h>

Inheritance diagram for PorousFlowPeacemanBorehole:

Public Member Functions
	PorousFlowPeacemanBorehole (const InputParameters &parameters)
	Creates a new PorousFlowPeacemanBorehole This reads the file containing the lines of the form radius x y z that defines the borehole geometry. More...

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

Protected Member Functions
Real	wellConstant (const RealTensorValue &perm, const RealTensorValue &rot, const Real &half_len, const Elem *ele, const Real &rad) const
	Calculates Peaceman's form of the borehole well constant Z Chen, Y Zhang, Well flow models for various numerical methods, Int J Num Analysis and Modeling, 3 (2008) 375-388. More...

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
const Function &	_character
	If positive then the borehole acts as a sink (producion well) for porepressure > borehole pressure, and does nothing otherwise If negative then the borehole acts as a source (injection well) for porepressure < borehole pressure, and does nothing otherwise The flow rate to/from the borehole is multiplied by \|character\|, so usually character = +/- 1. More...

const Real	_p_bot
	Bottomhole pressure of borehole. More...

const RealVectorValue	_unit_weight
	Unit weight of fluid in borehole (for calculating bottomhole pressure at each Dirac Point) More...

const Real	_re_constant
	Borehole constant. More...

const Real	_well_constant
	Well constant. More...

const bool	_has_permeability
	Whether there is a quadpoint permeability material (for error checking) More...

const bool	_has_thermal_conductivity
	Whether there is a quadpoint thermal conductivity material (for error checking) More...

const MaterialProperty< RealTensorValue > &	_perm_or_cond
	Permeability or conductivity of porous material. More...

const MaterialProperty< std::vector< RealTensorValue > > &	_dperm_or_cond_dvar
	d(Permeability)/d(PorousFlow variable) More...

std::vector< RealTensorValue >	_rot_matrix
	Rotation matrix used in well_constant calculation. 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 borehole by a sequence of Dirac Points.

Definition at line 22 of file PorousFlowPeacemanBorehole.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

◆ PorousFlowPeacemanBorehole()

PorousFlowPeacemanBorehole::PorousFlowPeacemanBorehole ( const InputParameters & parameters )

Creates a new PorousFlowPeacemanBorehole This reads the file containing the lines of the form radius x y z that defines the borehole geometry.

It also calculates segment-lengths and rotation matrices needed for computing the borehole well constant

Definition at line 69 of file PorousFlowPeacemanBorehole.C.

   : PorousFlowLineSink(parameters),
     _character(getFunction("character")),
     _p_bot(getParam<Real>("bottom_p_or_t")),
     _unit_weight(getParam<RealVectorValue>("unit_weight")),
     _re_constant(getParam<Real>("re_constant")),
     _well_constant(getParam<Real>("well_constant")),
     _has_permeability(
         hasMaterialProperty<RealTensorValue>("PorousFlow_permeability_qp") &&
         hasMaterialProperty<std::vector<RealTensorValue>>("dPorousFlow_permeability_qp_dvar")),
     _has_thermal_conductivity(
         hasMaterialProperty<RealTensorValue>("PorousFlow_thermal_conductivity_qp") &&
         hasMaterialProperty<std::vector<RealTensorValue>>(
             "dPorousFlow_thermal_conductivity_qp_dvar")),
     _perm_or_cond(_p_or_t == PorTchoice::pressure
                       ? getMaterialProperty<RealTensorValue>("PorousFlow_permeability_qp")
                       : getMaterialProperty<RealTensorValue>("PorousFlow_thermal_conductivity_qp")),
     _dperm_or_cond_dvar(
         _p_or_t == PorTchoice::pressure
             ? getMaterialProperty<std::vector<RealTensorValue>>("dPorousFlow_permeability_qp_dvar")
             : getMaterialProperty<std::vector<RealTensorValue>>(
                   "dPorousFlow_thermal_conductivity_qp_dvar"))
 {
   if (_p_or_t == PorTchoice::pressure && !_has_permeability)
     mooseError("PorousFlowPeacemanBorehole: You have specified function_of=porepressure, but you "
                "do not have a quadpoint permeability material");
   if (_p_or_t == PorTchoice::temperature && !_has_thermal_conductivity)
     mooseError("PorousFlowPeacemanBorehole: You have specified function_of=temperature, but you do "
                "not have a quadpoint thermal_conductivity material");
  
   // construct the rotation matrix needed to rotate the permeability
   const unsigned int num_pts = _zs.size();
   _rot_matrix.resize(std::max(num_pts - 1, (unsigned)1));
   for (unsigned int i = 0; i + 1 < num_pts; ++i)
   {
     const RealVectorValue v2(_xs[i + 1] - _xs[i], _ys[i + 1] - _ys[i], _zs[i + 1] - _zs[i]);
     _rot_matrix[i] = RotationMatrix::rotVecToZ(v2);
   }
   if (num_pts == (unsigned)1)
     _rot_matrix[0] = RotationMatrix::rotVecToZ(_line_direction);
 }

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 PorousFlowPeacemanBorehole::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 216 of file PorousFlowPeacemanBorehole.C.

 {
   const Real character = _character.value(_t, _q_point[_qp]);
   if (character == 0.0)
     return 0.0;
  
   const Real bh_pressure = _p_bot + _unit_weight * (_q_point[_qp] - _bottom_point);
   const Real pp = ptqp();
  
   Real outflow = 0.0; // this is the flow rate from porespace out of the system
  
   if (current_dirac_ptid > 0)
   // contribution from half-segment "behind" this point (must have >1 point for
   // current_dirac_ptid>0)
   {
     if ((character < 0.0 && pp < bh_pressure) || (character > 0.0 && pp > bh_pressure))
     {
       // injection, so outflow<0 || production, so outflow>0
       const Real wc = wellConstant(_perm_or_cond[_qp],
                                    _rot_matrix[current_dirac_ptid - 1],
                                    _half_seg_len[current_dirac_ptid - 1],
                                    _current_elem,
                                    _rs[current_dirac_ptid]);
       outflow += wc * (pp - bh_pressure);
     }
   }
  
   if (current_dirac_ptid + 1 < _zs.size() || _zs.size() == 1)
   // contribution from half-segment "ahead of" this point, or we only have one point
   {
     if ((character < 0.0 && pp < bh_pressure) || (character > 0.0 && pp > bh_pressure))
     {
       // injection, so outflow<0 || // production, so outflow>0
       const Real wc = wellConstant(_perm_or_cond[_qp],
                                    _rot_matrix[current_dirac_ptid],
                                    _half_seg_len[current_dirac_ptid],
                                    _current_elem,
                                    _rs[current_dirac_ptid]);
       outflow += wc * (pp - bh_pressure);
     }
   }
  
   return outflow * _test[_i][_qp] * std::abs(character);
 }

◆ computeQpBaseOutflowJacobian()

void PorousFlowPeacemanBorehole::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 262 of file PorousFlowPeacemanBorehole.C.

 {
   outflow = 0.0;
   outflowp = 0.0;
  
   const Real character = _character.value(_t, _q_point[_qp]);
   if (character == 0.0)
     return;
  
   if (_dictator.notPorousFlowVariable(jvar))
     return;
   const unsigned pvar = _dictator.porousFlowVariableNum(jvar);
  
   const Real bh_pressure = _p_bot + _unit_weight * (_q_point[_qp] - _bottom_point);
   const Real pp = ptqp();
   const Real pp_prime = dptqp(pvar) * _phi[_j][_qp];
  
   if (current_dirac_ptid > 0)
   // contribution from half-segment "behind" this point
   {
     if ((character < 0.0 && pp < bh_pressure) || (character > 0.0 && pp > bh_pressure))
     {
       // injection, so outflow<0 || // production, so outflow>0
       const Real wc = wellConstant(_perm_or_cond[_qp],
                                    _rot_matrix[current_dirac_ptid - 1],
                                    _half_seg_len[current_dirac_ptid - 1],
                                    _current_elem,
                                    _rs[current_dirac_ptid]);
       outflowp += wc * pp_prime;
       outflow += wc * (pp - bh_pressure);
     }
   }
  
   if (current_dirac_ptid < _zs.size() - 1 || _zs.size() == 1)
   // contribution from half-segment "ahead of" this point
   {
     if ((character < 0.0 && pp < bh_pressure) || (character > 0.0 && pp > bh_pressure))
     {
       // injection, so outflow<0 || // production, so outflow>0
       const Real wc = wellConstant(_perm_or_cond[_qp],
                                    _rot_matrix[current_dirac_ptid],
                                    _half_seg_len[current_dirac_ptid],
                                    _current_elem,
                                    _rs[current_dirac_ptid]);
       outflowp += wc * pp_prime;
       outflow += wc * (pp - bh_pressure);
     }
   }
  
   outflowp *= _test[_i][_qp] * std::abs(character);
   outflow *= _test[_i][_qp] * std::abs(character);
 }

◆ 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 PorousFlowPolyLineSink::computeQpBaseOutflowJacobian(), and 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 PorousFlowPolyLineSink::computeQpBaseOutflow(), computeQpBaseOutflow(), PorousFlowPolyLineSink::computeQpBaseOutflowJacobian(), and computeQpBaseOutflowJacobian().

◆ wellConstant()

Real PorousFlowPeacemanBorehole::wellConstant	(	const RealTensorValue &	perm,
		const RealTensorValue &	rot,
		const Real &	half_len,
		const Elem *	ele,
		const Real &	rad
	)		const

protected

Calculates Peaceman's form of the borehole well constant Z Chen, Y Zhang, Well flow models for various numerical methods, Int J Num Analysis and Modeling, 3 (2008) 375-388.

Definition at line 112 of file PorousFlowPeacemanBorehole.C.

 {
   if (_well_constant > 0)
     return _well_constant;
  
   // rot_perm has its "2" component lying along the half segment.
   // We want to determine the eigenvectors of rot(0:1, 0:1), since, when
   // rotated back to the original frame we will determine the element
   // lengths along these directions
   const RealTensorValue rot_perm = (rot * perm) * rot.transpose();
   const Real trace2D = rot_perm(0, 0) + rot_perm(1, 1);
   const Real det2D = rot_perm(0, 0) * rot_perm(1, 1) - rot_perm(0, 1) * rot_perm(1, 0);
   const Real sq = std::sqrt(std::max(0.25 * trace2D * trace2D - det2D,
                                      0.0)); // the std::max accounts for wierdo precision loss
   const Real eig_val1 = 0.5 * trace2D + sq;
   const Real eig_val2 = 0.5 * trace2D - sq;
   RealVectorValue eig_vec1, eig_vec2;
   if (sq > std::abs(trace2D) * 1E-7) // matrix is not a multiple of the identity (1E-7 accounts for
                                      // precision in a crude way)
   {
     if (rot_perm(1, 0) != 0)
     {
       eig_vec1(0) = eig_val1 - rot_perm(1, 1);
       eig_vec1(1) = rot_perm(1, 0);
       eig_vec2(0) = eig_val2 - rot_perm(1, 1);
       eig_vec2(1) = rot_perm(1, 0);
     }
     else if (rot_perm(0, 1) != 0)
     {
       eig_vec1(0) = rot_perm(0, 1);
       eig_vec1(1) = eig_val1 - rot_perm(0, 0);
       eig_vec2(0) = rot_perm(0, 1);
       eig_vec2(1) = eig_val2 - rot_perm(0, 0);
     }
     else // off diagonal terms are both zero
     {
       eig_vec1(0) = 1.0;
       eig_vec2(1) = 1.0;
     }
   }
   else // matrix is basically a multiple of the identity
   {
     eig_vec1(0) = 1.0;
     eig_vec2(1) = 1.0;
   }
  
   // finally, rotate these to original frame and normalise
   eig_vec1 = rot.transpose() * eig_vec1;
   eig_vec1 /= std::sqrt(eig_vec1 * eig_vec1);
   eig_vec2 = rot.transpose() * eig_vec2;
   eig_vec2 /= std::sqrt(eig_vec2 * eig_vec2);
  
   // find the "length" of the element in these directions
   // TODO - maybe better to use variance than max&min
   Real max1 = eig_vec1 * ele->point(0);
   Real max2 = eig_vec2 * ele->point(0);
   Real min1 = max1;
   Real min2 = max2;
   Real proj;
   for (unsigned int i = 1; i < ele->n_nodes(); i++)
   {
     proj = eig_vec1 * ele->point(i);
     max1 = (max1 < proj) ? proj : max1;
     min1 = (min1 < proj) ? min1 : proj;
  
     proj = eig_vec2 * ele->point(i);
     max2 = (max2 < proj) ? proj : max2;
     min2 = (min2 < proj) ? min2 : proj;
   }
   const Real ll1 = max1 - min1;
   const Real ll2 = max2 - min2;
  
   Real r0;
   if (eig_val1 <= 0.0)
     r0 = _re_constant * ll1;
   else if (eig_val2 <= 0.0)
     r0 = _re_constant * ll2;
   else
     r0 = _re_constant * std::sqrt(std::sqrt(eig_val1 / eig_val2) * std::pow(ll2, 2) +
                                   std::sqrt(eig_val2 / eig_val1) * std::pow(ll1, 2)) /
          (std::pow(eig_val1 / eig_val2, 0.25) + std::pow(eig_val2 / eig_val1, 0.25));
  
   const Real effective_perm = (det2D >= 0.0 ? std::sqrt(det2D) : 0.0);
  
   const Real halfPi = acos(0.0);
  
   if (r0 <= rad)
     mooseError("The effective element size (about 0.2-times-true-ele-size) for an element "
                "containing a Peaceman-type borehole must be (much) larger than the borehole radius "
                "for the Peaceman formulation to be correct.  Your element has effective size ",
                r0,
                " and the borehole radius is ",
                rad,
                "\n");
  
   return 4 * halfPi * effective_perm * half_len / std::log(r0 / rad);
 }

Referenced by computeQpBaseOutflow(), and 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 computeQpBaseOutflow(), computeQpBaseOutflowJacobian(), and PorousFlowLineGeometry::PorousFlowLineGeometry().

◆ _character

const Function& PorousFlowPeacemanBorehole::_character

protected

If positive then the borehole acts as a sink (producion well) for porepressure > borehole pressure, and does nothing otherwise If negative then the borehole acts as a source (injection well) for porepressure < borehole pressure, and does nothing otherwise The flow rate to/from the borehole is multiplied by |character|, so usually character = +/- 1.

Definition at line 43 of file PorousFlowPeacemanBorehole.h.

Referenced by computeQpBaseOutflow(), and computeQpBaseOutflowJacobian().

◆ _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 PorousFlowPolyLineSink::computeQpBaseOutflowJacobian(), 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.

◆ _dperm_or_cond_dvar

const MaterialProperty<std::vector<RealTensorValue> >& PorousFlowPeacemanBorehole::_dperm_or_cond_dvar

protected

d(Permeability)/d(PorousFlow variable)

Definition at line 67 of file PorousFlowPeacemanBorehole.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 PorousFlowPolyLineSink::computeQpBaseOutflow(), computeQpBaseOutflow(), PorousFlowPolyLineSink::computeQpBaseOutflowJacobian(), 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_permeability

const bool PorousFlowPeacemanBorehole::_has_permeability

protected

Whether there is a quadpoint permeability material (for error checking)

Definition at line 58 of file PorousFlowPeacemanBorehole.h.

Referenced by PorousFlowPeacemanBorehole().

◆ _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().

◆ _has_thermal_conductivity

const bool PorousFlowPeacemanBorehole::_has_thermal_conductivity

protected

Whether there is a quadpoint thermal conductivity material (for error checking)

Definition at line 61 of file PorousFlowPeacemanBorehole.h.

Referenced by PorousFlowPeacemanBorehole().

◆ _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().

◆ _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_bot

const Real PorousFlowPeacemanBorehole::_p_bot

protected

Bottomhole pressure of borehole.

Definition at line 46 of file PorousFlowPeacemanBorehole.h.

Referenced by computeQpBaseOutflow(), and computeQpBaseOutflowJacobian().

◆ _p_or_t

enum PorousFlowLineSink::PorTchoice PorousFlowLineSink::_p_or_t

protectedinherited

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

◆ _perm_or_cond

const MaterialProperty<RealTensorValue>& PorousFlowPeacemanBorehole::_perm_or_cond

protected

Permeability or conductivity of porous material.

Definition at line 64 of file PorousFlowPeacemanBorehole.h.

Referenced by computeQpBaseOutflow(), and computeQpBaseOutflowJacobian().

◆ _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().

◆ _re_constant

const Real PorousFlowPeacemanBorehole::_re_constant

protected

Borehole constant.

Definition at line 52 of file PorousFlowPeacemanBorehole.h.

◆ _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.

◆ _rot_matrix

std::vector<RealTensorValue> PorousFlowPeacemanBorehole::_rot_matrix

protected

Rotation matrix used in well_constant calculation.

Definition at line 70 of file PorousFlowPeacemanBorehole.h.

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

◆ _rs

std::vector<Real> PorousFlowLineGeometry::_rs

protectedinherited

Radii of the borehole.

Definition at line 49 of file PorousFlowLineGeometry.h.

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

◆ _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().

◆ _unit_weight

const RealVectorValue PorousFlowPeacemanBorehole::_unit_weight

protected

Unit weight of fluid in borehole (for calculating bottomhole pressure at each Dirac Point)

Definition at line 49 of file PorousFlowPeacemanBorehole.h.

Referenced by computeQpBaseOutflow(), and computeQpBaseOutflowJacobian().

◆ _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().

◆ _well_constant

const Real PorousFlowPeacemanBorehole::_well_constant

protected

Well constant.

Definition at line 55 of file PorousFlowPeacemanBorehole.h.

◆ _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().

◆ _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().

◆ _zs

std::vector<Real> PorousFlowLineGeometry::_zs

protectedinherited

z points of borehole

Definition at line 58 of file PorousFlowLineGeometry.h.

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

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

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

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Detailed Description

Member Enumeration Documentation

◆ PorTchoice

Constructor & Destructor Documentation

◆ PorousFlowPeacemanBorehole()

Member Function Documentation

◆ addPoints()

◆ computeQpBaseOutflow()

◆ computeQpBaseOutflowJacobian()

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

◆ dptqp()

◆ jac()

◆ parseNextLineReals()

◆ ptqp()

◆ wellConstant()

Member Data Documentation

◆ _bottom_point

◆ _character

◆ _denthalpy_dvar

◆ _dfluid_density_node_dvar

◆ _dfluid_viscosity_dvar

◆ _dictator

◆ _dinternal_energy_dvar

◆ _dmass_fractions_dvar

◆ _dperm_or_cond_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_permeability

◆ _has_porepressure

◆ _has_relative_permeability

◆ _has_temperature

◆ _has_thermal_conductivity

◆ _internal_energy

◆ _line_direction

◆ _line_length

◆ _mass_fractions

◆ _p_bot

◆ _p_or_t

◆ _perm_or_cond

◆ _ph

◆ _point_file

◆ _pp

◆ _re_constant

◆ _relative_permeability

◆ _rot_matrix

◆ _rs

◆ _sp

◆ _temperature

◆ _total_outflow_mass

◆ _unit_weight

◆ _use_enthalpy

◆ _use_internal_energy

◆ _use_mass_fraction

◆ _use_mobility

◆ _use_relative_permeability

◆ _well_constant

◆ _xs

◆ _ys

◆ _zs