diff --git a/ProcessLib/HT/MonolithicHTFEM.h b/ProcessLib/HT/MonolithicHTFEM.h index 9bbdb9ce77806298397a8b52ff5d9bad232e7fa1..e86f9dc7da88bba5266a4ba4bd6034d2fbeafbc9 100644 --- a/ProcessLib/HT/MonolithicHTFEM.h +++ b/ProcessLib/HT/MonolithicHTFEM.h @@ -81,9 +81,9 @@ public: auto local_b = MathLib::createZeroedVector<LocalVectorType>( local_b_data, local_matrix_size); - auto Ktt = local_K.template block<temperature_size, temperature_size>( + auto KTT = local_K.template block<temperature_size, temperature_size>( temperature_index, temperature_index); - auto Mtt = local_M.template block<temperature_size, temperature_size>( + auto MTT = local_M.template block<temperature_size, temperature_size>( temperature_index, temperature_index); auto Kpp = local_K.template block<pressure_size, pressure_size>( pressure_index, pressure_index); @@ -178,13 +178,13 @@ public: this->getThermalConductivityDispersivity( vars, porosity, fluid_density, specific_heat_capacity_fluid, velocity, I); - Ktt.noalias() += + KTT.noalias() += (dNdx.transpose() * thermal_conductivity_dispersivity * dNdx + N.transpose() * velocity.transpose() * dNdx * fluid_density * specific_heat_capacity_fluid) * w; Kpp.noalias() += w * dNdx.transpose() * K_over_mu * dNdx; - Mtt.noalias() += + MTT.noalias() += w * this->getHeatEnergyCoefficient(vars, porosity, fluid_density, specific_heat_capacity_fluid) *