ACInterfaceKobayashi1 Class Reference

Kernel 1 of 2 for interfacial energy anisotropy in the Allen-Cahn equation as implemented in R. More...

#include <ACInterfaceKobayashi1.h>

Inheritance diagram for ACInterfaceKobayashi1:

Public Member Functions
	ACInterfaceKobayashi1 (const InputParameters &parameters)

Protected Member Functions
virtual RealGradient	precomputeQpResidual ()
	Enum of computeDFDOP inputs. More...
virtual RealGradient	precomputeQpJacobian ()
virtual Real	computeQpOffDiagJacobian (unsigned int jvar)

Protected Attributes
const MaterialProperty< Real > &	_L
	Mobility. More...
const MaterialProperty< Real > &	_dLdop
const MaterialProperty< Real > &	_eps
const MaterialProperty< Real > &	_deps
const MaterialProperty< RealGradient > &	_depsdgrad_op
const MaterialProperty< RealGradient > &	_ddepsdgrad_op
std::vector< const MaterialProperty< Real > * >	_dLdarg
	Mobility derivative w.r.t. other coupled variables. More...
std::vector< const MaterialProperty< Real > * >	_depsdarg
std::vector< const MaterialProperty< Real > * >	_ddepsdarg

Detailed Description

Kernel 1 of 2 for interfacial energy anisotropy in the Allen-Cahn equation as implemented in R.

Kobayashi, Physica D, 63, 410-423 (1993). doi:10.1016/0167-2789(93)90120-P This kernel implements the first two terms on the right side of eq. (3) of the paper.

Definition at line 27 of file ACInterfaceKobayashi1.h.

Constructor & Destructor Documentation

ACInterfaceKobayashi1()

ACInterfaceKobayashi1::ACInterfaceKobayashi1

(

const InputParameters &

parameters

)

Definition at line 36 of file ACInterfaceKobayashi1.C.

  : DerivativeMaterialInterface<JvarMapKernelInterface<KernelGrad>>(parameters),
    _L(getMaterialProperty<Real>("mob_name")),
    _dLdop(getMaterialPropertyDerivative<Real>("mob_name", _var.name())),
    _eps(getMaterialProperty<Real>("eps_name")),
    _deps(getMaterialProperty<Real>("deps_name")),
    _depsdgrad_op(getMaterialProperty<RealGradient>("depsdgrad_op_name")),
    _ddepsdgrad_op(getMaterialProperty<RealGradient>("ddepsdgrad_op_name"))
{
  // Get number of coupled variables
  unsigned int nvar = _coupled_moose_vars.size();
 
  // reserve space for derivatives
  _dLdarg.resize(nvar);
  _depsdarg.resize(nvar);
  _ddepsdarg.resize(nvar);
 
  // Iterate over all coupled variables
  for (unsigned int i = 0; i < nvar; ++i)
  {
    const VariableName iname = _coupled_moose_vars[i]->name();
    _dLdarg[i] = &getMaterialPropertyDerivative<Real>("mob_name", iname);
    _depsdarg[i] = &getMaterialPropertyDerivative<Real>("eps_name", iname);
    _ddepsdarg[i] = &getMaterialPropertyDerivative<Real>("deps_name", iname);
  }
}

Member Function Documentation

computeQpOffDiagJacobian()

Real ACInterfaceKobayashi1::computeQpOffDiagJacobian ( unsigned int  jvar )

protectedvirtual

Definition at line 96 of file ACInterfaceKobayashi1.C.

{
  // get the coupled variable jvar is referring to
  const unsigned int cvar = mapJvarToCvar(jvar);
 
  // Set modified gradient vector
  const RealGradient v(-_grad_u[_qp](1), _grad_u[_qp](0), 0);
 
  // Set off-diagonal jaocbian terms from mobility dependence
  Real dsum =
      _L[_qp] * (_deps[_qp] * (*_depsdarg[cvar])[_qp] * _phi[_j][_qp] * v * _grad_test[_i][_qp]);
  dsum +=
      _L[_qp] * (_eps[_qp] * (*_ddepsdarg[cvar])[_qp] * _phi[_j][_qp] * v * _grad_test[_i][_qp]);
  dsum += (*_dLdarg[cvar])[_qp] * _phi[_j][_qp] * _eps[_qp] * _deps[_qp] * v * _grad_test[_i][_qp];
 
  return dsum;
}

precomputeQpJacobian()

RealGradient ACInterfaceKobayashi1::precomputeQpJacobian ( )

protectedvirtual

Definition at line 74 of file ACInterfaceKobayashi1.C.

{
  // Set modified gradient vector
  const RealGradient v(-_grad_u[_qp](1), _grad_u[_qp](0), 0);
 
  // dvdgrad_op*_grad_phi
  const RealGradient dv(-_grad_phi[_j][_qp](1), _grad_phi[_j][_qp](0), 0);
 
  // Derivative of epsilon wrt nodal op values
  Real depsdop_i = _depsdgrad_op[_qp] * _grad_phi[_j][_qp];
  Real ddepsdop_i = _ddepsdgrad_op[_qp] * _grad_phi[_j][_qp];
  ;
 
  // Set the Jacobian
  RealGradient jac1 = _eps[_qp] * _deps[_qp] * dv;
  RealGradient jac2 = _deps[_qp] * depsdop_i * v;
  RealGradient jac3 = _eps[_qp] * ddepsdop_i * v;
 
  return _L[_qp] * (jac1 + jac2 + jac3);
}

precomputeQpResidual()

RealGradient ACInterfaceKobayashi1::precomputeQpResidual ( )

protectedvirtual

Enum of computeDFDOP inputs.

Definition at line 64 of file ACInterfaceKobayashi1.C.

{
  // Set modified gradient vector
  const RealGradient v(-_grad_u[_qp](1), _grad_u[_qp](0), 0);
 
  // Define anisotropic interface residual
  return _eps[_qp] * _deps[_qp] * _L[_qp] * v;
}

Member Data Documentation

_ddepsdarg

std::vector<const MaterialProperty<Real> *> ACInterfaceKobayashi1::_ddepsdarg

protected

Definition at line 49 of file ACInterfaceKobayashi1.h.

Referenced by ACInterfaceKobayashi1(), and computeQpOffDiagJacobian().

_ddepsdgrad_op

const MaterialProperty<RealGradient>& ACInterfaceKobayashi1::_ddepsdgrad_op

protected

Definition at line 44 of file ACInterfaceKobayashi1.h.

Referenced by precomputeQpJacobian().

_deps

const MaterialProperty<Real>& ACInterfaceKobayashi1::_deps

protected

Definition at line 42 of file ACInterfaceKobayashi1.h.

Referenced by computeQpOffDiagJacobian(), precomputeQpJacobian(), and precomputeQpResidual().

_depsdarg

std::vector<const MaterialProperty<Real> *> ACInterfaceKobayashi1::_depsdarg

protected

Definition at line 48 of file ACInterfaceKobayashi1.h.

Referenced by ACInterfaceKobayashi1(), and computeQpOffDiagJacobian().

_depsdgrad_op

const MaterialProperty<RealGradient>& ACInterfaceKobayashi1::_depsdgrad_op

protected

Definition at line 43 of file ACInterfaceKobayashi1.h.

Referenced by precomputeQpJacobian().

_dLdarg

std::vector<const MaterialProperty<Real> *> ACInterfaceKobayashi1::_dLdarg

protected

Mobility derivative w.r.t. other coupled variables.

Definition at line 47 of file ACInterfaceKobayashi1.h.

Referenced by ACInterfaceKobayashi1(), and computeQpOffDiagJacobian().

_dLdop

const MaterialProperty<Real>& ACInterfaceKobayashi1::_dLdop

protected

Definition at line 40 of file ACInterfaceKobayashi1.h.

_eps

const MaterialProperty<Real>& ACInterfaceKobayashi1::_eps

protected

Definition at line 41 of file ACInterfaceKobayashi1.h.

Referenced by computeQpOffDiagJacobian(), precomputeQpJacobian(), and precomputeQpResidual().

_L

const MaterialProperty<Real>& ACInterfaceKobayashi1::_L

protected

Mobility.

Definition at line 39 of file ACInterfaceKobayashi1.h.

Referenced by computeQpOffDiagJacobian(), precomputeQpJacobian(), and precomputeQpResidual().

---

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

• phase_field/include/kernels/ACInterfaceKobayashi1.h
• phase_field/src/kernels/ACInterfaceKobayashi1.C