21 params.
addParam<MaterialPropertyName>(
"mob_name",
"L",
"The mobility used with the kernel");
22 params.
addParam<std::string>(
"base_name",
"Material property base name");
27 "Range over which order parameters change across an " 28 "interface. By default order parameters are assumed to " 35 _L(getMaterialProperty<
Real>(
"mob_name")),
36 _base_name(isParamValid(
"base_name") ? getParam<
std::string>(
"base_name") +
"_" :
""),
37 _strain(getMaterialPropertyByName<
RankTwoTensor>(_base_name +
"elastic_strain")),
38 _stress(getParam<
Real>(
"stress") / getParam<
Real>(
"op_range"))
57 const Real ds = -
_stress / (grad_norm * grad_norm_sq);
58 const Real dsx = ds * nx;
59 const Real dsy = ds * ny;
60 const Real dsz = ds * nz;
63 _dS(0, 0, 0) = (ny * ny + nz * nz) * dsx;
64 _dS(0, 1, 0) =
_dS(0, 0, 1) = -ny * s - nx * ny * dsx;
65 _dS(0, 1, 1) = 2.0 * nx * s + (nx * nx + nz * nz) * dsx;
66 _dS(0, 2, 0) =
_dS(0, 0, 2) = -nz * s - nx * nz * dsx;
67 _dS(0, 2, 1) =
_dS(0, 1, 2) = -ny * nz * dsx;
68 _dS(0, 2, 2) = 2.0 * nx * s + (nx * nx + ny * ny) * dsx;
71 _dS(1, 0, 0) = 2.0 * ny * s + (ny * ny + nz * nz) * dsy;
72 _dS(1, 1, 0) =
_dS(1, 0, 1) = -nx * s - nx * ny * dsy;
73 _dS(1, 1, 1) = (nx * nx + nz * nz) * dsy;
74 _dS(1, 2, 0) =
_dS(1, 0, 2) = -nx * nz * dsy;
75 _dS(1, 2, 1) =
_dS(1, 1, 2) = -nz * s - ny * nz * dsy;
76 _dS(1, 2, 2) = 2.0 * ny * s + (nx * nx + ny * ny) * dsy;
79 _dS(2, 0, 0) = 2.0 * nz * s + (ny * ny + nz * nz) * dsz;
80 _dS(2, 1, 0) =
_dS(2, 0, 1) = -nx * ny * dsz;
81 _dS(2, 1, 1) = 2.0 * nz * s + (nx * nx + nz * nz) * dsz;
82 _dS(2, 2, 0) =
_dS(2, 0, 2) = -nx * s - nx * nz * dsz;
83 _dS(2, 2, 1) =
_dS(2, 1, 2) = -ny * s - ny * nz * dsz;
84 _dS(2, 2, 2) = (nx * nx + ny * ny) * dsz;
108 const Real ds = -
_stress / (grad_norm * grad_norm_sq);
109 const Real dsx = ds * nx;
110 const Real dsy = ds * ny;
111 const Real dsz = ds * nz;
113 const Real dus = ds * (nx * px + ny * py + pz * nz);
115 const Real b = -3.0 * nx * px - 3.0 * ny * py - 3.0 * nz * pz;
116 const Real dudsx = ds * nx / grad_norm_sq *
b + ds * px;
117 const Real dudsy = ds * ny / grad_norm_sq *
b + ds * py;
118 const Real dudsz = ds * nz / grad_norm_sq *
b + ds * pz;
121 _ddS(0, 0, 0) = (2.0 * ny * py + 2.0 * nz * pz) * dsx + (ny * ny + nz * nz) * dudsx;
123 -py * s - ny * dus - px * ny * dsx - nx * py * dsx - nx * ny * dudsx;
124 _ddS(0, 1, 1) = 2.0 * px * s + 2.0 * nx * dus + (2.0 * nx * px + 2.0 * nz * pz) * dsx +
125 (nx * nx + nz * nz) * dudsx;
127 -pz * s - nz * dus - px * nz * dsx - nx * pz * dsx - nx * nz * dudsx;
128 _ddS(0, 2, 1) =
_ddS(0, 1, 2) = -py * nz * dsx - ny * pz * dsx - ny * nz * dudsx;
129 _ddS(0, 2, 2) = 2.0 * px * s + 2.0 * nx * dus + (2.0 * nx * px + 2.0 * ny * py) * dsx +
130 (nx * nx + ny * ny) * dudsx;
133 _ddS(1, 0, 0) = 2.0 * py * s + 2.0 * ny * dus + (2.0 * ny * py + 2.0 * nz * pz) * dsy +
134 (ny * ny + nz * nz) * dudsy;
136 -px * s - nx * dus - px * ny * dsy - nx * py * dsy - nx * ny * dudsy;
137 _ddS(1, 1, 1) = (2.0 * nx * px + 2.0 * nz * pz) * dsy + (nx * nx + nz * nz) * dudsy;
138 _ddS(1, 2, 0) =
_ddS(1, 0, 2) = -px * nz * dsy - nx * pz * dsy - nx * nz * dudsy;
140 -pz * s - nz * dus - py * nz * dsy - ny * pz * dsy - ny * nz * dudsy;
141 _ddS(1, 2, 2) = 2.0 * py * s + 2.0 * ny * dus + (2.0 * nx * px + 2.0 * ny * py) * dsy +
142 (nx * nx + ny * ny) * dudsy;
145 _ddS(2, 0, 0) = 2.0 * pz * s + 2.0 * nz * dus + (2.0 * ny * py + 2.0 * nz * pz) * dsz +
146 (ny * ny + nz * nz) * dudsz;
147 _ddS(2, 1, 0) =
_ddS(2, 0, 1) = -px * ny * dsz - nx * py * dsz - nx * ny * dudsz;
148 _ddS(2, 1, 1) = 2.0 * pz * s + 2.0 * nz * dus + (2.0 * nx * px + 2.0 * nz * pz) * dsz +
149 (nx * nx + nz * nz) * dudsz;
151 -px * s - nx * dus - px * nz * dsz - nx * pz * dsz - nx * nz * dudsz;
153 -py * s - ny * dus - py * nz * dsz - ny * pz * dsz - ny * nz * dudsz;
154 _ddS(2, 2, 2) = (2.0 * nx * px + 2.0 * ny * py) * dsz + (nx * nx + ny * ny) * dudsz;
ACInterfaceStress(const InputParameters ¶meters)
const MaterialProperty< Real > & _L
Mobility.
const VariableGradient & _grad_u
static InputParameters validParams()
registerMooseObject("PhaseFieldApp", ACInterfaceStress)
static constexpr Real TOLERANCE
const VariablePhiGradient & _grad_phi
virtual Real computeQpJacobian() override
ADRealEigenVector< T, D, asd > sqrt(const ADRealEigenVector< T, D, asd > &)
static InputParameters validParams()
Compute the Allen-Cahn interface stress driving force contribution .
RankThreeTensor _dS
d sigma/d(grad eta), derivative of interface stress tensor with order parameter gradient ...
virtual Real computeQpResidual() override
const Real _stress
interface stress
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const VariableTestGradient & _grad_test
VectorValue< T > doubleContraction(const RankTwoTensorTempl< T > &b) const
RankThreeTensor _ddS
derivative of _dS w.r.t. the finite element coefficients for the Jacobian calculation ...
const MaterialProperty< RankTwoTensor > & _strain