diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h index c1368c241a0eee7a44f0bc41d4a88730db8d58d5..15288e84d58318dc9cc48768092678f7d098ec96 100644 --- a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h +++ b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h @@ -587,17 +587,6 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement, double const chi_S_L = chi(S_L); double const chi_S_L_prev = chi(S_L_prev); - double const k_rel = - medium->property(MPL::PropertyType::relative_permeability) - .template value<double>(variables, x_position, t, dt); - auto const mu = liquid_phase.property(MPL::PropertyType::viscosity) - .template value<double>(variables, x_position, t, dt); - auto const K_intrinsic = MPL::formEigenTensor<DisplacementDim>( - solid_phase.property(MPL::PropertyType::permeability) - .value(variables, x_position, t, dt)); - - GlobalDimMatrixType const rho_Ki_over_mu = K_intrinsic * rho_LR / mu; - variables[static_cast<int>( MPL::Variable::effective_pore_pressure_rate)] = (chi_S_L * (-p_cap_ip) - @@ -626,6 +615,17 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement, dt; } + double const k_rel = + medium->property(MPL::PropertyType::relative_permeability) + .template value<double>(variables, x_position, t, dt); + auto const mu = liquid_phase.property(MPL::PropertyType::viscosity) + .template value<double>(variables, x_position, t, dt); + auto const K_intrinsic = MPL::formEigenTensor<DisplacementDim>( + solid_phase.property(MPL::PropertyType::permeability) + .value(variables, x_position, t, dt)); + + GlobalDimMatrixType const rho_Ki_over_mu = K_intrinsic * rho_LR / mu; + // // displacement equation, displacement part //