# -*- coding: utf-8 -*- import numpy as np import paraBEM from paraBEM.vtk_export import VtkWriter from paraBEM import PanelVector2, Vector2, Panel2 from paraBEM.airfoil import Airfoil from paraBEM.pan2d import NeumannSource0Case2 as Case from paraBEM.utils import check_path airfoil = Airfoil.joukowsky(m=-0.1 + 0.1j) airfoil.numpoints = 50 alpha = np.deg2rad(10) # panelmethode case = Case(airfoil.panels) case.v_inf = Vector2(np.cos(alpha), np.sin(alpha)) case.run() print(np.array(case.matrix.values)) nx = 300 ny = 300 space_x = np.linspace(-1, 2, nx) space_y = np.linspace(-0.2, 0.2, ny) vec = lambda x: paraBEM.Vector2(x[0], x[1]) vec3 = lambda x: [x[0], x[1], 0] grid = [paraBEM.Vector2(x, y) for y in space_y for x in space_x] velocity = list(map(vec3, map(case.off_body_velocity, grid))) vel1 = [(i[0]**2 + i[1]**2)**(0.5) for i in velocity]
# -*- coding: utf-8 -*- import numpy as np import paraBEM from paraBEM.vtk_export import VtkWriter from paraBEM import PanelVector2, Vector2, Panel2 from paraBEM.airfoil import Airfoil from paraBEM.pan2d import NeumannSource0Case2 as Case from paraBEM.utils import check_path airfoil = Airfoil.joukowsky(m=-0.1 +0.1j) airfoil.numpoints = 50 alpha = np.deg2rad(10) # panelmethode case = Case(airfoil.panels) case.v_inf = Vector2(np.cos(alpha), np.sin(alpha)) case.run() print(np.array(case.matrix.values)) nx = 300 ny = 300 space_x = np.linspace(-1, 2, nx) space_y = np.linspace(-0.2, 0.2, ny) vec = lambda x: paraBEM.Vector2(x[0], x[1]) vec3 = lambda x: [x[0], x[1], 0] grid = [paraBEM.Vector2(x, y) for y in space_y for x in space_x] velocity = list(map(vec3, map(case.off_body_velocity, grid))) vel1 = [(i[0]**2 + i[1]**2)**(0.5) for i in velocity]
from paraBEM.vtk_export import CaseToVTK
from paraBEM.pan3d import DirichletDoublet0Source0Case3 as Case
from paraBEM.airfoil import Airfoil
def rib3d(airfoil, y_pos):
out = [paraBEM.PanelVector3(coo[0], y_pos, coo[1]) for coo in airfoil.coordinates[:-1]]
out.append(out[0])
out[0].wake_vertex = True
return out
n_x = 50
n_y = 10
a = Airfoil.joukowsky(-0.01+1j)
a.numpoints = n_x
print(a.coordinates)
ribs = [rib3d(a, y) for y in np.linspace(-5, 5, n_y)]
panels = []
for i in range(n_y)[:-1]:
for j in range(n_x):
panels.append(paraBEM.Panel3([ribs[i][j], ribs[i + 1][j], ribs[i + 1][j + 1], ribs[i][j + 1]]))
te = [rib[0] for rib in ribs]
print(te)
case = Case(panels, te)
case.farfield = 5
case.v_inf = paraBEM.Vector3(1, 0, 0.0)
case.create_wake(length=10000, count=3) # length, count
case.run()
