fix bugs

YiminJin · YiminJin · commit 60124a101bab · 2025-06-17T13:59:53.000-07:00
diff --git a/source/material_model/grain_size.cc b/source/material_model/grain_size.cc
@@ -198,7 +198,7 @@ namespace aspect
         }
 
       const double strain_rate_dependence = (1.0 - dislocation_creep_exponent[phase_index]) / dislocation_creep_exponent[phase_index];
-      const SymmetricTensor<2,dim> shear_strain_rate = Utilities::Tensors::deviator(shear_strain_rate);
+      const SymmetricTensor<2,dim> shear_strain_rate = Utilities::Tensors::deviator(strain_rate);
       const double second_strain_rate_invariant = std::sqrt(std::max(-Utilities::Tensors::deviatoric_tensor_inv2(shear_strain_rate), 0.));
 
       // If the strain rate is zero, the dislocation viscosity is infinity.

Original file line number	Diff line number	Diff line change
`@@ -198,7 +198,7 @@ namespace aspect`
`198`	`198`	`}`
`199`	`199`
`200`	`200`	`const double strain_rate_dependence = (1.0 - dislocation_creep_exponent[phase_index]) / dislocation_creep_exponent[phase_index];`
`201`		`- const SymmetricTensor<2,dim> shear_strain_rate = Utilities::Tensors::deviator(shear_strain_rate);`
	`201`	`+ const SymmetricTensor<2,dim> shear_strain_rate = Utilities::Tensors::deviator(strain_rate);`
`202`	`202`	`const double second_strain_rate_invariant = std::sqrt(std::max(-Utilities::Tensors::deviatoric_tensor_inv2(shear_strain_rate), 0.));`
`203`	`203`
`204`	`204`	`// If the strain rate is zero, the dislocation viscosity is infinity.`