From 0abf06b73b1d75c74d68e2b8a0b069e024ede195 Mon Sep 17 00:00:00 2001
From: Paul Dechamps <paulzer@Pauls-MacBook-Pro.local>
Date: Fri, 15 Apr 2022 14:20:01 +0200
Subject: [PATCH] Testing without Rbf
---
dart/pyVII/vCoupler.py | 2 +-
dart/src/wViscSolver.cpp | 5 +++++
dart/tests/bli_lowLift.py | 8 ++++----
3 files changed, 10 insertions(+), 5 deletions(-)
diff --git a/dart/pyVII/vCoupler.py b/dart/pyVII/vCoupler.py
index 449f2a2..e0d5813 100644
--- a/dart/pyVII/vCoupler.py
+++ b/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
diff --git a/dart/src/wViscSolver.cpp b/dart/src/wViscSolver.cpp
index 586cdac..84f5388 100644
--- a/dart/src/wViscSolver.cpp
+++ b/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
diff --git a/dart/tests/bli_lowLift.py b/dart/tests/bli_lowLift.py
index 78b802f..eddabfa 100644
--- a/dart/tests/bli_lowLift.py
+++ b/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
--
GitLab