import paraBEM from paraBEM.pan3d import doublet_3_0_n0, doublet_3_0_sphere, doublet_3_0_vsaero from paraBEM.utils import check_path pnt1 = paraBEM.PanelVector3(-0.5, -0.5, 0) pnt2 = paraBEM.PanelVector3(0.5, -0.5, 0) pnt3 = paraBEM.PanelVector3(0.5, 0.5, 0) pnt4 = paraBEM.PanelVector3(-0.5, 0.5, 0) source = paraBEM.Panel3([pnt1, pnt2, pnt3, pnt4]) x = np.linspace(-0, 5, 500) y = [] for xi in x: target = paraBEM.PanelVector3(xi, 0., 0.1) panel_infl = doublet_3_0_vsaero(target, source) point_infl = doublet_3_0_n0(target, source) y.append([panel_infl, point_infl, abs(panel_infl - point_infl)]) y = list(zip(*y)) plt.figure(figsize=(8, 3)) plt.gcf().subplots_adjust(bottom=0.15) plt.ylabel("Einfluss") plt.xlabel("x") plt.plot(x, y[0], label=u"Exakte Lösung") plt.plot(x, y[1], label=u"Näherungslösung") plt.ylim(-0.1, 0.6) plt.xlim(0, None) plt.legend() plt.grid(True) plt.savefig(check_path("results/3d/far_vs_near_doublet.png"))
import paraBEM from paraBEM import pan3d from paraBEM.mesh import mesh_object mesh = mesh_object.from_OBJ("../mesh/box_minimal.obj") case = pan3d.DirichletDoublet0Case3(mesh.panels) case.v_inf = paraBEM.Vector3(1, 0, 0) a = case.panels[0] b = case.panels[1] print(a.center, " ", a.n) print(b.center, " ", b.n) print(pan3d.doublet_3_0_vsaero(a.center, b))
import paraBEM from paraBEM.pan3d import doublet_3_0_n0, doublet_3_0_sphere, doublet_3_0_vsaero from paraBEM.utils import check_path pnt1 = paraBEM.PanelVector3(-0.5, -0.5, 0) pnt2 = paraBEM.PanelVector3(0.5, -0.5, 0) pnt3 = paraBEM.PanelVector3(0.5, 0.5, 0) pnt4 = paraBEM.PanelVector3(-0.5, 0.5, 0) source = paraBEM.Panel3([pnt1, pnt2, pnt3, pnt4]) x = np.linspace(-0, 5, 500) y = [] for xi in x: target = paraBEM.PanelVector3(xi, 0., 0.1) panel_infl = doublet_3_0_vsaero(target, source) point_infl = doublet_3_0_n0(target, source) y.append([panel_infl, point_infl, abs(panel_infl - point_infl)]) y = list(zip(*y)) plt.figure(figsize=(8,3)) plt.gcf().subplots_adjust(bottom=0.15) plt.ylabel("Einfluss") plt.xlabel("x") plt.plot(x, y[0], label=u"Exakte Lösung") plt.plot(x, y[1], label=u"Näherungslösung") plt.ylim(-0.1, 0.6) plt.xlim(0, None) plt.legend() plt.grid(True)