def naca_5_digit(des, npointside):
npointside_p = xi.copy_intp(npointside)
xout_a = xi.new_doublea(2 * npointside)
zout_a = xi.new_doublea(2 * npointside)
nout_p = xi.new_intp()
stat_p = xi.new_intp()
xi.naca_5_digit(des, npointside_p, xout_a, zout_a, nout_p, stat_p)
nout = xi.intp_value(nout_p)
stat = xi.intp_value(stat_p)
if stat != 0:
nout = 1
xout = nout * [0]
zout = nout * [0]
for i in range(nout):
xout[i] = xi.doublea_getitem(xout_a, i)
zout[i] = xi.doublea_getitem(zout_a, i)
xi.delete_intp(npointside_p)
xi.delete_doublea(xout_a)
xi.delete_doublea(zout_a)
xi.delete_intp(nout_p)
xi.delete_intp(stat_p)
return xout, zout, nout, stat
def naca_4_digit(camber, xcamber, thick, npointside):
camber_p = xi.copy_doublep(camber)
xcamber_p = xi.copy_doublep(xcamber)
thick_p = xi.copy_doublep(thick)
npointside_p = xi.copy_intp(npointside)
xout_a = xi.new_doublea(2 * npointside)
zout_a = xi.new_doublea(2 * npointside)
nout_p = xi.new_intp()
xi.naca_4_digit(camber_p, xcamber_p, thick_p, npointside_p, xout_a, zout_a,
nout_p)
nout = xi.intp_value(nout_p)
xout = nout * [0]
zout = nout * [0]
for i in range(nout):
xout[i] = xi.doublea_getitem(xout_a, i)
zout[i] = xi.doublea_getitem(zout_a, i)
xi.delete_doublep(camber_p)
xi.delete_doublep(xcamber_p)
xi.delete_doublep(thick_p)
xi.delete_intp(npointside_p)
xi.delete_doublea(xout_a)
xi.delete_doublea(zout_a)
xi.delete_intp(nout_p)
return xout, zout, nout
def xfoil_speccl(xdg, cl_spec):
cl_spec_p = xi.copy_doublep(cl_spec)
alpha_p = xi.new_doublep()
lift_p = xi.new_doublep()
drag_p = xi.new_doublep()
moment_p = xi.new_doublep()
converged_p = xi.new_boolp()
stat_p = xi.new_intp()
xi.xfoil_speccl(xdg, cl_spec_p, alpha_p, lift_p, drag_p, moment_p,
converged_p, stat_p)
alpha = xi.doublep_value(alpha_p)
lift = xi.doublep_value(lift_p)
drag = xi.doublep_value(drag_p)
moment = xi.doublep_value(moment_p)
converged = xi.boolp_value(converged_p)
stat = xi.intp_value(stat_p)
xi.delete_doublep(cl_spec_p)
xi.delete_doublep(alpha_p)
xi.delete_doublep(lift_p)
xi.delete_doublep(drag_p)
xi.delete_doublep(moment_p)
xi.delete_boolp(converged_p)
xi.delete_intp(stat_p)
return alpha, lift, drag, moment, converged, stat
def xfoil_modify_tegap(xdg, gap, blendloc):
stat_p = xi.new_intp()
gap_p = xi.copy_doublep(gap)
blendloc_p = xi.copy_doublep(blendloc)
npointout_p = xi.new_intp()
xi.xfoil_modify_tegap(xdg, gap_p, blendloc_p, npointout_p, stat_p)
npointout = xi.intp_value(npointout_p)
stat = xi.intp_value(stat_p)
xi.delete_doublep(gap_p)
xi.delete_doublep(blendloc_p)
xi.delete_intp(npointout_p)
xi.delete_intp(stat_p)
return npointout, stat
def xfoil_get_wakepoints(xdg):
nwake_p = xi.new_intp()
xi.xfoil_get_wakepoints(xdg, nwake_p)
nwake = xi.intp_value(nwake_p)
xi.delete_intp(nwake_p)
return nwake
def xfoil_smooth_paneling(xdg):
stat_p = xi.new_intp()
xi.xfoil_smooth_paneling(xdg, stat_p)
stat = xi.intp_value(stat_p)
xi.delete_intp(stat_p)
return stat
def xfoil_apply_flap_deflection(xdg, xflap, zflap, z_flap_spec, degrees):
xflap_p = xi.copy_doublep(xflap)
zflap_p = xi.copy_doublep(zflap)
z_flap_spec_p = xi.copy_intp(z_flap_spec)
degrees_p = xi.copy_doublep(degrees)
npointout_p = xi.new_intp()
stat_p = xi.new_intp()
xi.xfoil_apply_flap_deflection(xdg, xflap_p, zflap_p, z_flap_spec_p,
degrees_p, npointout_p, stat_p)
npointout = xi.intp_value(npointout_p)
stat = xi.intp_value(stat_p)
xi.delete_doublep(xflap_p)
xi.delete_doublep(zflap_p)
xi.delete_intp(z_flap_spec_p)
xi.delete_intp(npointout_p)
xi.delete_intp(stat_p)
return npointout, stat
def xfoil_get_current_airfoil(xdg, npoint):
xout_a = xi.new_doublea(npoint)
zout_a = xi.new_doublea(npoint)
npoint_p = xi.copy_intp(npoint)
stat_p = xi.new_intp()
xi.xfoil_get_current_airfoil(xdg, xout_a, zout_a, npoint_p, stat_p)
xout = npoint * [0]
zout = npoint * [0]
for i in range(npoint):
xout[i] = xi.doublea_getitem(xout_a, i)
zout[i] = xi.doublea_getitem(zout_a, i)
stat = xi.intp_value(stat_p)
xi.delete_doublea(xout_a)
xi.delete_doublea(zout_a)
xi.delete_intp(npoint_p)
xi.delete_intp(stat_p)
return xout, zout, stat
def xfoil_geometry_info(xdg):
maxt_p = xi.new_doublep()
xmaxt_p = xi.new_doublep()
maxc_p = xi.new_doublep()
xmaxc_p = xi.new_doublep()
stat_p = xi.new_intp()
xi.xfoil_geometry_info(xdg, maxt_p, xmaxt_p, maxc_p, xmaxc_p, stat_p)
maxt = xi.doublep_value(maxt_p)
xmaxt = xi.doublep_value(xmaxt_p)
maxc = xi.doublep_value(maxc_p)
xmaxc = xi.doublep_value(xmaxc_p)
stat = xi.intp_value(stat_p)
xi.delete_doublep(maxt_p)
xi.delete_doublep(xmaxt_p)
xi.delete_doublep(maxc_p)
xi.delete_doublep(xmaxc_p)
xi.delete_intp(stat_p)
return maxt, xmaxt, maxc, xmaxc, stat
def run_xfoil(npointin, xin, zin, geom_opts, noppoint, operating_points,
op_modes, reynolds_numbers, mach_numbers, use_flap, x_flap,
z_flap, z_flap_spec, flap_degrees, xfoil_opts, reinitialize,
fix_unconverged):
npointin_p = xi.copy_intp(npointin)
xin_a = xi.new_doublea(npointin)
zin_a = xi.new_doublea(npointin)
noppoint_p = xi.copy_intp(noppoint)
operating_points_a = xi.new_doublea(noppoint)
op_modes_a = xi.new_inta(noppoint)
reynolds_numbers_a = xi.new_doublea(noppoint)
mach_numbers_a = xi.new_doublea(noppoint)
use_flap_p = xi.copy_boolp(use_flap)
if use_flap:
x_flap_p = xi.copy_doublep(x_flap)
z_flap_p = xi.copy_doublep(z_flap)
z_flap_spec_p = xi.copy_intp(z_flap_spec)
else:
x_flap_p = xi.new_doublep()
z_flap_p = xi.new_doublep()
z_flap_spec_p = xi.new_intp()
flap_degrees_a = xi.new_doublea(noppoint)
reinitialize_p = xi.copy_boolp(reinitialize)
fix_unconverged_p = xi.copy_boolp(fix_unconverged)
lift_a = xi.new_doublea(noppoint)
drag_a = xi.new_doublea(noppoint)
moment_a = xi.new_doublea(noppoint)
viscrms_a = xi.new_doublea(noppoint)
alpha_a = xi.new_doublea(noppoint)
xtrt_a = xi.new_doublea(noppoint)
xtrb_a = xi.new_doublea(noppoint)
stat_p = xi.new_intp()
for i in range(npointin):
xi.doublea_setitem(xin_a, i, xin[i])
xi.doublea_setitem(zin_a, i, zin[i])
for i in range(noppoint):
xi.doublea_setitem(operating_points_a, i, operating_points[i])
xi.inta_setitem(op_modes_a, i, op_modes[i])
xi.doublea_setitem(reynolds_numbers_a, i, reynolds_numbers[i])
xi.doublea_setitem(mach_numbers_a, i, mach_numbers[i])
if use_flap:
xi.doublea_setitem(flap_degrees_a, i, flap_degrees[i])
xi.run_xfoil(npointin_p, xin_a, zin_a, geom_opts, noppoint_p,
operating_points_a, op_modes_a, reynolds_numbers_a,
mach_numbers_a, use_flap_p, x_flap_p, z_flap_p, z_flap_spec_p,
flap_degrees_a, xfoil_opts, reinitialize_p, fix_unconverged_p,
lift_a, drag_a, moment_a, viscrms_a, alpha_a, xtrt_a, xtrb_a,
stat_p)
stat = xi.intp_value(stat_p)
lift = noppoint * [0]
drag = noppoint * [0]
moment = noppoint * [0]
viscrms = noppoint * [0]
alpha = noppoint * [0]
xtrt = noppoint * [0]
xtrb = noppoint * [0]
for i in range(noppoint):
lift[i] = xi.doublea_getitem(lift_a, i)
drag[i] = xi.doublea_getitem(drag_a, i)
moment[i] = xi.doublea_getitem(moment_a, i)
viscrms[i] = xi.doublea_getitem(viscrms_a, i)
alpha[i] = xi.doublea_getitem(alpha_a, i)
xtrt[i] = xi.doublea_getitem(xtrt_a, i)
xtrb[i] = xi.doublea_getitem(xtrb_a, i)
xi.delete_intp(npointin_p)
xi.delete_doublea(xin_a)
xi.delete_doublea(zin_a)
xi.delete_intp(noppoint_p)
xi.delete_doublea(operating_points_a)
xi.delete_inta(op_modes_a)
xi.delete_doublea(reynolds_numbers_a)
xi.delete_doublea(mach_numbers_a)
xi.delete_boolp(use_flap_p)
xi.delete_doublep(x_flap_p)
xi.delete_doublep(z_flap_p)
xi.delete_intp(z_flap_spec_p)
xi.delete_doublea(flap_degrees_a)
xi.delete_doublea(lift_a)
xi.delete_doublea(drag_a)
xi.delete_doublea(moment_a)
xi.delete_doublea(viscrms_a)
xi.delete_doublea(alpha_a)
xi.delete_doublea(xtrt_a)
xi.delete_doublea(xtrb_a)
xi.delete_intp(stat_p)
return lift, drag, moment, viscrms, alpha, xtrt, xtrb, stat
