#include <GBEvolutionBase.h>

Inheritance diagram for GBEvolutionBase:

Public Member Functions
	GBEvolutionBase (const InputParameters &parameters)

Protected Member Functions
virtual void	computeQpProperties ()

Protected Attributes
Real	_f0s

Real	_wGB

Real	_length_scale

Real	_time_scale

Real	_GBmob0

Real	_Q

Real	_GBMobility

Real	_molar_vol

const VariableValue &	_T

MaterialProperty< Real > &	_sigma

MaterialProperty< Real > &	_M_GB

MaterialProperty< Real > &	_kappa

MaterialProperty< Real > &	_gamma

MaterialProperty< Real > &	_L

MaterialProperty< Real > *	_dLdT

MaterialProperty< Real > &	_l_GB

MaterialProperty< Real > &	_mu

MaterialProperty< Real > &	_entropy_diff

MaterialProperty< Real > &	_molar_volume

MaterialProperty< Real > &	_act_wGB

const Real	_kb

const Real	_JtoeV

Detailed Description

Definition at line 21 of file GBEvolutionBase.h.

Constructor & Destructor Documentation

◆ GBEvolutionBase()

GBEvolutionBase::GBEvolutionBase ( const InputParameters & parameters )

Definition at line 36 of file GBEvolutionBase.C.

   : DerivativeMaterialInterface<Material>(parameters),
     _f0s(getParam<Real>("f0s")),
     _wGB(getParam<Real>("wGB")),
     _length_scale(getParam<Real>("length_scale")),
     _time_scale(getParam<Real>("time_scale")),
     _GBmob0(getParam<Real>("GBmob0")),
     _Q(getParam<Real>("Q")),
     _GBMobility(getParam<Real>("GBMobility")),
     _molar_vol(getParam<Real>("molar_volume")),
     _T(coupledValue("T")),
     _sigma(declareProperty<Real>("sigma")),
     _M_GB(declareProperty<Real>("M_GB")),
     _kappa(declareProperty<Real>("kappa_op")),
     _gamma(declareProperty<Real>("gamma_asymm")),
     _L(declareProperty<Real>("L")),
     _dLdT(parameters.hasDefaultCoupledValue("T")
               ? nullptr
               : &declarePropertyDerivative<Real>("L", getVar("T", 0)->name())),
     _l_GB(declareProperty<Real>("l_GB")),
     _mu(declareProperty<Real>("mu")),
     _entropy_diff(declareProperty<Real>("entropy_diff")),
     _molar_volume(declareProperty<Real>("molar_volume")),
     _act_wGB(declareProperty<Real>("act_wGB")),
     _kb(8.617343e-5),     // Boltzmann constant in eV/K
     _JtoeV(6.24150974e18) // Joule to eV conversion
 {
   if (_GBMobility == -1 && _GBmob0 == 0)
     mooseError("Either a value for GBMobility or for GBmob0 and Q must be provided");
 }

Member Function Documentation

◆ computeQpProperties()

void GBEvolutionBase::computeQpProperties ( )

protectedvirtual

Reimplemented in GBEvolution.

Definition at line 68 of file GBEvolutionBase.C.

 {
   const Real length_scale4 = _length_scale * _length_scale * _length_scale * _length_scale;
  
   // GB mobility Derivative
   Real dM_GBdT;
  
   if (_GBMobility < 0)
   {
     // Convert to lengthscale^4/(eV*timescale);
     const Real M0 = _GBmob0 * _time_scale / (_JtoeV * length_scale4);
  
     _M_GB[_qp] = M0 * std::exp(-_Q / (_kb * _T[_qp]));
     dM_GBdT = _M_GB[_qp] * _Q / (_kb * _T[_qp] * _T[_qp]);
   }
   else
   {
     // Convert to lengthscale^4/(eV*timescale)
     _M_GB[_qp] = _GBMobility * _time_scale / (_JtoeV * length_scale4);
     dM_GBdT = 0.0;
   }
  
   // in the length scale of the system
   _l_GB[_qp] = _wGB;
  
   _L[_qp] = 4.0 / 3.0 * _M_GB[_qp] / _l_GB[_qp];
   if (_dLdT)
     (*_dLdT)[_qp] = 4.0 / 3.0 * dM_GBdT / _l_GB[_qp];
   _kappa[_qp] = 3.0 / 4.0 * _sigma[_qp] * _l_GB[_qp];
   _gamma[_qp] = 1.5;
   _mu[_qp] = 3.0 / 4.0 * 1.0 / _f0s * _sigma[_qp] / _l_GB[_qp];
  
   // J/(K mol) converted to eV(K mol)
   _entropy_diff[_qp] = 8.0e3 * _JtoeV;
  
   // m^3/mol converted to ls^3/mol
   _molar_volume[_qp] = _molar_vol / (_length_scale * _length_scale * _length_scale);
   _act_wGB[_qp] = 0.5e-9 / _length_scale;
 }