Testing without Rbf

dart/pyVII/vCoupler.py

@@ -89,7 +89,7 @@ class Coupler:
       error = abs((self.vSolver.Cdt - cdPrev)/self.vSolver.Cdt) if self.vSolver.Cdt != 0 else 1
       
       if error  <= self.couplTol and exitCode==0:
-
+        print('{:>5s}| {:>10s} {:>10s} | {:>8s} {:>6s} {:>8s}'.format('It', 'Cl', 'Cd', 'xtrT', 'xtrB', 'Error'))
         print(ccolors.ANSI_GREEN,'{:>4.0f}| {:>10.5f} {:>10.5f} | {:>8.2f} {:>6.2f} {:>8.2f}\n'.format(couplIter,
         self.iSolver.Cl, self.vSolver.Cdt, xtr[0], xtr[1], np.log10(error)),ccolors.ANSI_RESET)
         return 0

dart/src/wViscSolver.cpp

@@ -220,6 +220,11 @@ int ViscSolver::Run(unsigned int couplIter)
 
                 pointExitCode = tSolver->Integration(iPoint, Sections[iSec][iRegion]);
 
+                if (Sections[iSec][iRegion]->Cf[iPoint] < 0)
+                {
+                    std::cout << "Negative Cf at point " << iPoint << " : " << Sections[iSec][iRegion]->Cf[iPoint] << std::endl;
+                }
+
                 /* if (iRegion == 1 && iPoint == 265)
                 {
                 std::cout << " iSec " << iSec

dart/tests/bli_lowLift.py

@@ -57,9 +57,9 @@ def main():
 
     # define flow variables
     Re = 1e7
-    alpha = 12*math.pi/180
+    alpha = 2*math.pi/180
     #alphaSeq = np.arange(-5, 5.5, 0.5)
-    M_inf = 0.
+    M_inf = 0.715
 
     meshFactor = 1
     CFL0 = 1
@@ -76,8 +76,8 @@ def main():
     # mesh the geometry
     print(ccolors.ANSI_BLUE + 'PyMeshing...' + ccolors.ANSI_RESET)
     tms['msh'].start()
-    pars = {'xLgt' : 5, 'yLgt' : 5, 'msF' : 1, 'msTe' : 0.01, 'msLe' : 0.001}
-    msh, gmshWriter = floD.mesh(dim, 'models/n0012.geo', pars, ['field', 'airfoil', 'downstream'])
+    pars = {'xLgt' : 5, 'yLgt' : 5, 'msF' : 1, 'msTe' : 0.01, 'msLe' : 0.005}
+    msh, gmshWriter = floD.mesh(dim, 'models/rae2822.geo', pars, ['field', 'airfoil', 'downstream'])
     tms['msh'].stop()
 
     # set the problem and add medium, initial/boundary conditions