Skip to content
Snippets Groups Projects
Commit 74524da7 authored by Thomas Nagel's avatar Thomas Nagel Committed by Dmitri Naumov
Browse files

Removed pressure coupling into displacement when flow in matrix is deactivated.

parent d81d6c23
No related branches found
No related tags found
No related merge requests found
...@@ -249,16 +249,14 @@ void HydroMechanicsLocalAssemblerMatrix<ShapeFunctionDisplacement, ...@@ -249,16 +249,14 @@ void HydroMechanicsLocalAssemblerMatrix<ShapeFunctionDisplacement,
rhs_u.noalias() -= -H_u.transpose() * rho * gravity_vec * ip_w; rhs_u.noalias() -= -H_u.transpose() * rho * gravity_vec * ip_w;
// //
// displacement equation, pressure part // pressure equation, pressure part and displacement equation, pressure
// // part
Kup.noalias() += B.transpose() * alpha * identity2 * N_p * ip_w;
//
// pressure equation, pressure part.
// //
if (!_process_data.deactivate_matrix_in_flow) // Only for hydraulically if (!_process_data.deactivate_matrix_in_flow) // Only for hydraulically
// active matrix // active matrix
{ {
Kup.noalias() += B.transpose() * alpha * identity2 * N_p * ip_w;
double const k_over_mu = double const k_over_mu =
_process_data.intrinsic_permeability(t, x_position)[0] / _process_data.intrinsic_permeability(t, x_position)[0] /
_process_data.fluid_viscosity(t, x_position)[0]; _process_data.fluid_viscosity(t, x_position)[0];
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment