Revert Dec commit 4e6928a2: related to heat flux depending on electric potential

NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp

@@ -523,7 +523,10 @@ void mlawLinearElecMagTherMech::constitutive(
             dThermalEMFieldSourcedGradV(m) += (-dedgradV(m,n) * (An(n) - A0(n)) * invdh);
 
             if (_useFluxT)
+            {
                 dThermalEMFieldSourcedGradV(m) += (- dedgradV(m,n) * gradV(n));
+                dThermalEMFieldSourcedA(m) += (-dfluxDdA(m,n) * gradV(n));
+            }
         }
         dThermalEMFieldSourcedA(m) += (-fluxD(m) * invdh);
         dThermalEMFieldSourcedB(m) += /*(H(m) * invdh)*/ 0.0;

NonLinearSolver/materialLaw/mlawLinearElecTherMech.cpp

@@ -206,8 +206,8 @@ if(evaluateCurlField)
     STensorOperation::zero(dkdT);
     STensorOperation::zero(dseebeckdT);
     STensorOperation::zero(djedV);
-    //STensorOperation::zero(dqdV);
-    //STensorOperation::zero(dqdF); // dependency of flux with deformation gradient
+    STensorOperation::zero(dqdV);
+    STensorOperation::zero(dqdF); // dependency of flux with deformation gradient
     STensorOperation::zero(djedF);
     STensorOperation::zero(dPdV);
 
@@ -257,7 +257,7 @@ if(evaluateCurlField)
 
 
 
-    /*for(int i=0;i<3;i++)
+    for(int i=0;i<3;i++)
      {
        fluxT(i)=0.;
        for(int j=0;j<3;j++)
@@ -265,7 +265,7 @@ if(evaluateCurlField)
 	 fluxT(i)+=(_k(i,j)+_l0(i,j)*_seebeck*_seebeck*T)*gradT(j)+_l0(i,j)*_seebeck*T*gradV(j);
          fluxjy(i)+=(_k(i,j)+_l0(i,j)*_seebeck*_seebeck*T+_l0(i,j)*_seebeck*V)*gradT(j)+(_l0(i,j)*V+_l0(i,j)*_seebeck*T)*gradV(j);
        }
-     }*/
+     }
 
      for(int i=0;i<3;i++)
       {
@@ -319,7 +319,7 @@ if(evaluateCurlField)
 
     djedgradV =_l0;
 
-    /*for(int i=0;i<3;i++)
+    for(int i=0;i<3;i++)
     {
       dqdT(i)=0.;
       for(int j=0;j<3;j++)
@@ -347,7 +347,7 @@ if(evaluateCurlField)
 	djydgradTdT(i,k)+=_l0(i,k)*_seebeck*_seebeck;
 	djydgradTdV(i,k)+=_l0(i,k)*_seebeck;
       }
-    }*/
+    }
    double dh;
    double dcpdt=0.;
    dh = getTimeStep();
@@ -373,7 +373,7 @@ if(evaluateCurlField)
 
 // to evaluate q or jy depending on je
 // in TM problem with EM coupling
-    {
+    /*{
         double dseebeckdT;
         static STensor3 dldT,dkdT;//,k1,k2,l10,l20;
 
@@ -428,7 +428,7 @@ if(evaluateCurlField)
                 djydgradTdV(i,k)+=_l0(i,k)*_seebeck;
             }
         }
-    }
+    }*/
 }
   
   double mlawLinearElecTherMech::deformationEnergy(STensor3 defo,const SVector3 &gradV,const SVector3 &gradT,const double T, const double V) const