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 PorousFlowPeacemanBorehole : public PorousFlowLineSink
{
public:
  static InputParameters validParams();
 
  PorousFlowPeacemanBorehole(const InputParameters & parameters);
 
  void virtual initialSetup() override;
 
protected:
  const Function & _character;
 
  const Function & _p_bot;
 
  const RealVectorValue _unit_weight;
 
  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;
};