Skip to content
Snippets Groups Projects
Commit a40b2af3 authored by wenqing's avatar wenqing Committed by Dmitri Naumov
Browse files

[THM/local assembler] Replaced shape function N_p with N_T in a matrix assembly

parent 7bc09a30
No related branches found
No related tags found
No related merge requests found
...@@ -419,7 +419,7 @@ void ThermoHydroMechanicsLocalAssembler<ShapeFunctionDisplacement, ...@@ -419,7 +419,7 @@ void ThermoHydroMechanicsLocalAssembler<ShapeFunctionDisplacement,
// Add thermo-osmosis effect on KTT // Add thermo-osmosis effect on KTT
KTT.noalias() += KTT.noalias() +=
(dNdx_T.transpose() * effective_thermal_condictivity * dNdx_T + (dNdx_T.transpose() * effective_thermal_condictivity * dNdx_T +
N_p.transpose() * velocity.transpose() * dNdx_p * fluid_density * c_f) * N_T.transpose() * velocity.transpose() * dNdx_T * fluid_density * c_f) *
w - w -
fluid_density * c_f * N_T.transpose() * (dNdx_T * T).transpose() * fluid_density * c_f * N_T.transpose() * (dNdx_T * T).transpose() *
K_pT_thermal_osmosis * dNdx_T * w; K_pT_thermal_osmosis * dNdx_T * w;
......
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