HEVPFlowRatePowerLawJ2 Class Reference

This user object classs Computes flow rate based on power law and Direction based on J2. More...

#include <HEVPFlowRatePowerLawJ2.h>

Inheritance diagram for HEVPFlowRatePowerLawJ2:

Public Member Functions
	HEVPFlowRatePowerLawJ2 (const InputParameters &parameters)
virtual bool	computeValue (unsigned int, Real &) const
virtual bool	computeDirection (unsigned int, RankTwoTensor &) const
virtual bool	computeDerivative (unsigned int, const std::string &, Real &) const
virtual bool	computeTensorDerivative (unsigned int, const std::string &, RankTwoTensor &) const

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
RankTwoTensor	computePK2Deviatoric (const RankTwoTensor &, const RankTwoTensor &) const
Real	computeEqvStress (const RankTwoTensor &, const RankTwoTensor &) const

Protected Attributes
Real	_ref_flow_rate
Real	_flow_rate_exponent
Real	_flow_rate_tol
std::string	_strength_prop_name
const std::string	_base_name
const MaterialProperty< Real > &	_strength
std::string	_pk2_prop_name
const MaterialProperty< RankTwoTensor > &	_pk2
const MaterialProperty< RankTwoTensor > &	_ce

Detailed Description

This user object classs Computes flow rate based on power law and Direction based on J2.

Definition at line 24 of file HEVPFlowRatePowerLawJ2.h.

Constructor & Destructor Documentation

HEVPFlowRatePowerLawJ2()

HEVPFlowRatePowerLawJ2::HEVPFlowRatePowerLawJ2

(

const InputParameters &

parameters

)

Definition at line 30 of file HEVPFlowRatePowerLawJ2.C.

  : HEVPFlowRateUOBase(parameters),
    _ref_flow_rate(getParam<Real>("reference_flow_rate")),
    _flow_rate_exponent(getParam<Real>("flow_rate_exponent")),
    _flow_rate_tol(getParam<Real>("flow_rate_tol"))
{
}

Member Function Documentation

computeDerivative()

bool HEVPFlowRatePowerLawJ2::computeDerivative ( unsigned int  qp, const std::string &  coupled_var_name, Real &  val  ) const

virtual

Implements HEVPFlowRateUOBase.

Definition at line 70 of file HEVPFlowRatePowerLawJ2.C.

{
  val = 0.0;
 
  if (_strength_prop_name == coupled_var_name)
  {
    RankTwoTensor pk2_dev = computePK2Deviatoric(_pk2[qp], _ce[qp]);
    Real eqv_stress = computeEqvStress(pk2_dev, _ce[qp]);
    val = -_ref_flow_rate * _flow_rate_exponent *
          std::pow(eqv_stress / _strength[qp], _flow_rate_exponent) / _strength[qp];
  }
 
  return true;
}

computeDirection()

bool HEVPFlowRatePowerLawJ2::computeDirection ( unsigned int  qp, RankTwoTensor &  val  ) const

virtual

Implements HEVPFlowRateUOBase.

Definition at line 57 of file HEVPFlowRatePowerLawJ2.C.

{
  RankTwoTensor pk2_dev = computePK2Deviatoric(_pk2[qp], _ce[qp]);
  Real eqv_stress = computeEqvStress(pk2_dev, _ce[qp]);
 
  val.zero();
  if (eqv_stress > 0.0)
    val = 1.5 / eqv_stress * _ce[qp] * pk2_dev * _ce[qp];
 
  return true;
}

computeEqvStress()

Real HEVPFlowRatePowerLawJ2::computeEqvStress ( const RankTwoTensor &  pk2_dev, const RankTwoTensor &  ce  ) const

protected

Definition at line 135 of file HEVPFlowRatePowerLawJ2.C.

{
  RankTwoTensor sdev = pk2_dev * ce;
  Real val = sdev.doubleContraction(sdev.transpose());
  return std::sqrt(1.5 * val);
}

Referenced by computeDerivative(), computeDirection(), computeTensorDerivative(), and computeValue().

computePK2Deviatoric()

RankTwoTensor HEVPFlowRatePowerLawJ2::computePK2Deviatoric ( const RankTwoTensor &  pk2, const RankTwoTensor &  ce  ) const

protected

Definition at line 128 of file HEVPFlowRatePowerLawJ2.C.

{
  return pk2 - (pk2.doubleContraction(ce) * ce.inverse()) / 3.0;
}

Referenced by computeDerivative(), computeDirection(), computeTensorDerivative(), and computeValue().

computeTensorDerivative()

bool HEVPFlowRatePowerLawJ2::computeTensorDerivative ( unsigned int  qp, const std::string &  coupled_var_name, RankTwoTensor &  val  ) const

virtual

Implements HEVPFlowRateUOBase.

Definition at line 88 of file HEVPFlowRatePowerLawJ2.C.

{
  val.zero();
 
  if (_pk2_prop_name == coupled_var_name)
  {
    RankTwoTensor pk2_dev = computePK2Deviatoric(_pk2[qp], _ce[qp]);
    Real eqv_stress = computeEqvStress(pk2_dev, _ce[qp]);
    Real dflowrate_dseqv = _ref_flow_rate * _flow_rate_exponent *
                           std::pow(eqv_stress / _strength[qp], _flow_rate_exponent - 1.0) /
                           _strength[qp];
 
    RankTwoTensor tau = pk2_dev * _ce[qp];
    RankTwoTensor dseqv_dpk2dev;
 
    dseqv_dpk2dev.zero();
    if (eqv_stress > 0.0)
      dseqv_dpk2dev = 1.5 / eqv_stress * tau * _ce[qp];
 
    RankTwoTensor ce_inv = _ce[qp].inverse();
 
    RankFourTensor dpk2dev_dpk2;
    for (unsigned int i = 0; i < LIBMESH_DIM; ++i)
      for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
        for (unsigned int k = 0; k < LIBMESH_DIM; ++k)
          for (unsigned int l = 0; l < LIBMESH_DIM; ++l)
          {
            dpk2dev_dpk2(i, j, k, l) = 0.0;
            if (i == k && j == l)
              dpk2dev_dpk2(i, j, k, l) = 1.0;
            dpk2dev_dpk2(i, j, k, l) -= ce_inv(i, j) * _ce[qp](k, l) / 3.0;
          }
    val = dflowrate_dseqv * dpk2dev_dpk2.transposeMajor() * dseqv_dpk2dev;
  }
  return true;
}

computeValue()

bool HEVPFlowRatePowerLawJ2::computeValue ( unsigned int  qp, Real &  val  ) const

virtual

Implements HEVPFlowRateUOBase.

Definition at line 39 of file HEVPFlowRatePowerLawJ2.C.

{
  RankTwoTensor pk2_dev = computePK2Deviatoric(_pk2[qp], _ce[qp]);
  Real eqv_stress = computeEqvStress(pk2_dev, _ce[qp]);
  val = std::pow(eqv_stress / _strength[qp], _flow_rate_exponent) * _ref_flow_rate;
 
  if (val > _flow_rate_tol)
  {
#ifdef DEBUG
    mooseWarning(
        "Flow rate greater than ", _flow_rate_tol, " ", val, " ", eqv_stress, " ", _strength[qp]);
#endif
    return false;
  }
  return true;
}

validParams()

InputParameters HEVPFlowRatePowerLawJ2::validParams ( )

static

Definition at line 17 of file HEVPFlowRatePowerLawJ2.C.

{
  InputParameters params = HEVPFlowRateUOBase::validParams();
  params.addParam<Real>(
      "reference_flow_rate", 0.001, "Reference flow rate for rate dependent flow");
  params.addParam<Real>("flow_rate_exponent", 10.0, "Power law exponent in flow rate equation");
  params.addParam<Real>("flow_rate_tol", 1e3, "Tolerance for flow rate");
  params.addClassDescription(
      "User object to evaluate power law flow rate and flow direction based on J2");
 
  return params;
}

Member Data Documentation

_base_name

const std::string HEVPFlowRateUOBase::_base_name

protectedinherited

Definition at line 40 of file HEVPFlowRateUOBase.h.

_ce

const MaterialProperty<RankTwoTensor>& HEVPFlowRateUOBase::_ce

protectedinherited

Definition at line 44 of file HEVPFlowRateUOBase.h.

Referenced by computeDerivative(), computeDirection(), computeTensorDerivative(), and computeValue().

_flow_rate_exponent

Real HEVPFlowRatePowerLawJ2::_flow_rate_exponent

protected

Definition at line 38 of file HEVPFlowRatePowerLawJ2.h.

Referenced by computeDerivative(), computeTensorDerivative(), and computeValue().

_flow_rate_tol

Real HEVPFlowRatePowerLawJ2::_flow_rate_tol

protected

Definition at line 39 of file HEVPFlowRatePowerLawJ2.h.

Referenced by computeValue().

_pk2

const MaterialProperty<RankTwoTensor>& HEVPFlowRateUOBase::_pk2

protectedinherited

Definition at line 43 of file HEVPFlowRateUOBase.h.

Referenced by computeDerivative(), computeDirection(), computeTensorDerivative(), and computeValue().

_pk2_prop_name

std::string HEVPFlowRateUOBase::_pk2_prop_name

protectedinherited

Definition at line 42 of file HEVPFlowRateUOBase.h.

Referenced by computeTensorDerivative().

_ref_flow_rate

Real HEVPFlowRatePowerLawJ2::_ref_flow_rate

protected

Definition at line 37 of file HEVPFlowRatePowerLawJ2.h.

Referenced by computeDerivative(), computeTensorDerivative(), and computeValue().

_strength

const MaterialProperty<Real>& HEVPFlowRateUOBase::_strength

protectedinherited

Definition at line 41 of file HEVPFlowRateUOBase.h.

Referenced by computeDerivative(), computeTensorDerivative(), and computeValue().

_strength_prop_name

std::string HEVPFlowRateUOBase::_strength_prop_name

protectedinherited

Definition at line 39 of file HEVPFlowRateUOBase.h.

Referenced by computeDerivative().

---

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

• tensor_mechanics/include/userobjects/HEVPFlowRatePowerLawJ2.h
• tensor_mechanics/src/userobjects/HEVPFlowRatePowerLawJ2.C