From 40c576192570b38b917bcd2b9a1ab3aa772a1d65 Mon Sep 17 00:00:00 2001
From: FLE_Knight <ujwalkishore.jinaga@mail.polimi.it>
Date: Fri, 26 Apr 2024 19:08:26 +0200
Subject: [PATCH] corrected
---
NonLinearSolver/materialLaw/mlawNonLinearTVENonLinearTVP.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/NonLinearSolver/materialLaw/mlawNonLinearTVENonLinearTVP.cpp b/NonLinearSolver/materialLaw/mlawNonLinearTVENonLinearTVP.cpp
index 45e52b6c6..a4f98f544 100644
--- a/NonLinearSolver/materialLaw/mlawNonLinearTVENonLinearTVP.cpp
+++ b/NonLinearSolver/materialLaw/mlawNonLinearTVENonLinearTVP.cpp
@@ -583,8 +583,8 @@ void mlawNonLinearTVENonLinearTVP::predictorCorrector_TVP_nonAssociatedFlow_nonL
DdevPhidGamma_RHS(i,j) += ( 0.5*2*_b*Ke*1/q1->_Jacobian*(ptilde+Gamma*dPhiPdGamma)*_I(i,j) - 3*pow(kk,1)*Hb/Cxdev*devPhi(i,j) ); // pow(kk,2) DEBUG
for (int k=0; k<3; k++)
for (int l=0; l<3; l++){
- DdevPhidDgamma_RHS(i,j) += 0.5* 2.*_b/3.*Gamma*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l))*_I(k,l)*Hp*1/q1->_Jacobian;
- DdevPhidGamma_RHS(i,j) += 0.5* ( -3.*Ge_Tensor(i,j,k,l)*1/q1->_Jacobian*devPhi(k,l) + 2.*_b/3.*Gamma*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l))*_I(k,l)*dPhiPdGamma );
+ DdevPhidDgamma_RHS(i,j) += 0.5* 2.*_b/3.*Gamma*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l)*1/q1->_Jacobian)*_I(k,l)*Hp;
+ DdevPhidGamma_RHS(i,j) += 0.5* ( -3.*Ge_Tensor(i,j,k,l)*1/q1->_Jacobian*devPhi(k,l) + 2.*_b/3.*Gamma*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l)*1/q1->_Jacobian)*_I(k,l)*dPhiPdGamma );
DdevPhidGamma_RHS(i,j) += 0.5*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l)*1/q1->_Jacobian)*N(k,l); // DEBUG
}
}
--
GitLab