PorousFlowPeacemanBorehole.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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 "PorousFlowLineSink.h"
 
class SinglePhaseFluidProperties;
class MooseVariableFieldBase;
 
namespace libMesh
{
class System;
}
 
class PorousFlowPeacemanBorehole : public PorousFlowLineSink
{
public:
  static InputParameters validParams();
 
  PorousFlowPeacemanBorehole(const InputParameters & parameters);
 
  void virtual initialSetup() override;
  void virtual residualSetup() override;
  void virtual jacobianSetup() override;
 
protected:
  const Function & _character;
 
  const Function & _p_bot;
 
  const RealVectorValue _unit_weight;
 
  const bool _use_density_from_temperature;
 
  const SinglePhaseFluidProperties * const _fp;
 
  const MooseVariableFieldBase * const _temperature_var;
 
  const libMesh::System * const _temperature_system;
 
  const unsigned int _temperature_var_number;
 
  const RealVectorValue _gravity;
 
  const Real _density_reference_pressure;
 
  const Real _t_c2k;
 
  std::vector<Real> _bh_pressure;
 
  void computeWellborePressures();
 
  Real wellborePressure(unsigned current_dirac_ptid) const;
 
  const Real _re_constant;
 
  const Real _well_constant;
 
  const bool _has_permeability;
 
  const bool _has_thermal_conductivity;
 
  const MaterialProperty<RealTensorValue> & _perm_or_cond;
 
  const MaterialProperty<std::vector<RealTensorValue>> & _dperm_or_cond_dvar;
 
  std::vector<RealTensorValue> _rot_matrix;
 
  Real wellConstant(const RealTensorValue & perm,
                    const RealTensorValue & rot,
                    const Real & half_len,
                    const Elem * ele,
                    const Real & rad) const;
 
  Real computeQpBaseOutflow(unsigned current_dirac_ptid) const override;
  void computeQpBaseOutflowJacobian(unsigned jvar,
                                    unsigned current_dirac_ptid,
                                    Real & outflow,
                                    Real & outflowp) const override;
};