n line

BattModels · Jan 29, 2019 · 75302c3 · 75302c3
1 parent bdd3260
commit 75302c3
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/kernels/KineticsLiqCrysVarNeg.C b/src/kernels/KineticsLiqCrysVarNeg.C
@@ -60,7 +60,7 @@ Real
 KineticsLiqCrysVarNeg::computeQpJacobian()
 {
   Real divn = _grad_n[_qp].tr();
-  Real gradderiv =  -2.0 * (_n[_qp] * _gradeta[_qp]) * divn -2.0 * _n[_qp] * (_grad_n[_qp] * _gradeta[_qp] + _grad_gradeta[_qp] * _n[_qp]); //There is a divergence in 1st term, not implemented yet                                                               
+  Real gradderiv =  -2.0 * (_n[_qp] * _gradeta[_qp]) * divn -2.0 * _n[_qp] * (_grad_n[_qp] * _gradeta[_qp] + _grad_gradeta[_qp] * _n[_qp]); //There is a divergence in 1st term, not implemented yet                                                
   Real epen = (- _dhdu[_qp]) * (_constfactor * _grad_n[_qp].contract(_grad_n[_qp]) + _penaltyconstfactor * (_n[_qp] * _gradeta[_qp]) * (_n[_qp] * _gradeta[_qp]) ) + _penaltyconstfactor * (1.0 - _h[_qp]) * gradderiv;
   Real pepenpuj = -_d2hdu2[_qp] * _phi[_j][_qp] * (_constfactor * _grad_n[_qp].contract(_grad_n[_qp]) + _penaltyconstfactor * (_n[_qp] * _gradeta[_qp]) * (_n[_qp] * _gradeta[_qp])) - _dhdu[_qp] * _phi[_j][_qp] * _penaltyconstfactor *  gradderiv;
   return -pepenpuj * exp(-epen) * _Fbv[_qp] * _test[_i][_qp] + exp( epen ) * _dFbvdu[_qp] * _phi[_j][_qp] * _test[_i][_qp];