Corrected a sign error for forcing function of neumann condition

labmec · May 31, 2024 · 88a50ea · 88a50ea
1 parent eeb9875
commit 88a50ea
Showing 1 changed file with 5 additions and 2 deletions.
diff --git a/Material/DarcyFlow/TPZDarcyFlow.cpp b/Material/DarcyFlow/TPZDarcyFlow.cpp
@@ -86,12 +86,12 @@ void TPZDarcyFlow::ContributeBC(const TPZMaterialDataT<STATE> &data, STATE weigh
         } else if(bc.Type() == 1) {
             v2 = 0.;
             for (int i = 0; i < fDim; i++) {
-                v2 += Flux[i] * normal[i];
+                v2 += -Flux[i] * normal[i];
             }
         } else if(bc.Type() == 2) {
             v2 = 0.;
             for (int i = 0; i < fDim; i++) {
-                v2 += Flux[i] * normal[i];
+                v2 += -Flux[i] * normal[i];
             }
             v2 += bc.Val1()(0,0) * rhs_val[0];
         }
@@ -396,4 +396,7 @@ void TPZDarcyFlow::FillBoundaryConditionDataRequirements(int type, TPZMaterialDa
     if (type == 3 || type == 1) {
         data.fNeedsNormal = true;
     }
+    if (HasForcingFunction()) {
+        data.fNeedsNormal = true;
+    }
 }