This material uses a strain-dependent diffusivity tensor that accommodates strain-induced anisotropy: $$\mathbf{d}=d_{0}\mathbf{I}+\sum\limits _{a=1}^{3}d_{1}^{a}\mathbf{m}_{a}+d_{2}^{a}\left(\mathbf{m}_{a}\cdot\mathbf{b}+\mathbf{b}\cdot\mathbf{m}_{a}\right),$$ where, $$\begin{aligned}d_{0} & =d_{0r}\left(\frac{J-\varphi_{r}^{s}}{1-\varphi_{r}^{s}}\right)^{\alpha_{^{0}}}e^{M_{^{0}}\left(J^{2}-1\right)/2}\\ d_{1}^{a} & =\frac{d_{1r}^{a}}{J^{2}}\left(\frac{J-\varphi_{r}^{s}}{1-\varphi_{r}^{s}}\right)^{\alpha_{^{a}}}e^{M_{^{a}}\left(J^{2}-1\right)/2}\\ d_{2}^{a} & =\frac{d_{2r}^{a}}{2J^{4}}\left(\frac{J-\varphi_{r}^{s}}{1-\varphi_{r}^{s}}\right)^{\alpha_{^{a}}}e^{M_{a}\left(J^{2}-1\right)/2} \end{aligned} ,\quad a=1,2,3$$ $J$ is the Jacobian of the deformation, i.e. $J=\det\mathbf{F}$ where $\mathbf{F}$ is the deformation gradient. $\mathbf{m}_{a}$ are second order tensors representing the spatial structural tensors describing the orthogonal planes of symmetry, given by $$\mathbf{m}_{a}=\mathbf{F}\cdot\left(\mathbf{V}_{a}\otimes\mathbf{V}_{a}\right)\cdot\mathbf{F}^{T},\quad a=1-3,$$ where $\mathbf{V}_{a}$ are orthonormal vectors normal to the planes of symmetry (defined as described in Section 4.1.1↑). Note that the diffusivity in the reference state ($\mathbf{F}=\mathbf{I})$ is given by $\mathbf{d}=d_{0r}\mathbf{I}+\sum\limits _{a=1}^{3}\left(d_{1r}^{a}+d_{2r}^{a}\right)\mathbf{V}_{a}\otimes\mathbf{V}_{a}$ .

Example:

<diffusivity name="Diffusivity" type="diff-ref-ortho">
  <phi0>0.2</phi0>
  <free_diff>0.005</free_diff>
  <diff00>0.001</diff00>
  <diff1>0.01, 0.02, 0.03</diff1>
  <diff2>0.001, 0.002, 0.003</diff2>
  <M0>0.5</M0>
  <M>1.5, 2.0, 2.5</M>
  <alpha0>1.5</alpha0>
  <alpha>2, 2.5, 3</alpha>
</diffusivity>