print("* wake parallel to the x-axis")
case.create_wake(50, 100, paraBEM.Vector3(1., 0., 0.))
case.run()
writer = CaseToVTK(case, "results/wake/x-parallel")
writer.write_panels()
writer.write_wake_panels()
print("F_l = ", case.cL * case.A_ref * 1.2 * 11**2 / 2.,
"F_w = ", case.cD * case.A_ref * 1.2 * 11**2 / 2.)
#wake parallel to v_inf
print("* wake parallel to v_inf")
case.create_wake(50, 100, case.v_inf)
case.run()
writer = CaseToVTK(case, "results/wake/v_inf-parallel")
writer.write_panels()
writer.write_wake_panels()
print("F_l = ", case.cL * case.A_ref * 1.2 * 11**2 / 2.,
"F_w = ", case.cD * case.A_ref * 1.2 * 11**2 / 2.)
#wake rollup to the x axis
print("* rolled up wake")
case.relax_wake(3, 1)
case.run()
writer = CaseToVTK(case, "results/wake/roll_up")
writer.write_panels()
writer.write_wake_panels()
print("F_l = ", case.cL * case.A_ref * 1.2 * 11**2 / 2.,
"F_w = ", case.cD * case.A_ref * 1.2 * 11**2 / 2.)
