From ca62dcef4305a8eb3fbfad12e2f039e2f9ac6d17 Mon Sep 17 00:00:00 2001
From: AmauryBilocq <amaury.bilocq@gmail.com>
Date: Wed, 30 Oct 2019 18:41:17 +0100
Subject: [PATCH] Trial to find the stagnation point
---
flow/viscous/solver.py | 10 +++++++++-
1 file changed, 9 insertions(+), 1 deletion(-)
diff --git a/flow/viscous/solver.py b/flow/viscous/solver.py
index c6b5a9cc..002f585a 100644
--- a/flow/viscous/solver.py
+++ b/flow/viscous/solver.py
@@ -20,6 +20,14 @@ class Solver():
def run(self):
''' Solve BLE
'''
+ # define tangential velocity in 2D
+ vt = self.v[:,0] + self.v[:,1]
+ # define stagnation point
+ stag = np.argmin(abs(self.v[:,1])) #define the stagnation point
+ print 'Stag is'
+ print stag
+ print len(self.boundary.nodes)
+
print '--- I am a fake BL solver and I am setting dummy blowing velocities! ---'
# 6th order fitting of delta_star for naca 0012 at alpha = 0, from XFoil
a = 0.000054925976379640116
@@ -32,4 +40,4 @@ class Solver():
for i in range(0, self.nE):
x = 0.5*(self.boundary.groups[0].tag.elems[i].nodes[0].pos.x[0]+self.boundary.groups[0].tag.elems[i].nodes[1].pos.x[0])
# CAUTION: a positive velocity will result in a suction!
- self.u[i] = -(b + 2*c*x + 3*d*x*x* + 4*e*x*x*x + 5*f*x*x*x*x + 6*g*x*x*x*x*x) # ue = d(delta_star) / dx
+ self.u[i] = -(b + 2*c*x + 3*d*x*x + 4*e*x*x*x + 5*f*x*x*x*x + 6*g*x*x*x*x*x) # ue = d(delta_star) / dx
--
GitLab