import sys import numpy import dakota_utils from build_mesh import build_mesh from apame_utils import run_case x,varnames,=dakota_utils.read_input(sys.argv[1]) filename = 'wing.vtp' build_mesh(filename, semi_span=x[0], root_chord=x[1], tip_chord=x[2], root_m=x[3], root_p=x[4], tip_m=x[5], tip_p=x[6], t=0.1, nu=50, nv=21) Patm = 101325.#Pa rho = 1.225#kg.m-3 c = numpy.sqrt(1.4*Patm/rho) Mach = 0.17 airspeed = c*Mach run_case(filename, wake_length=20., alpha=[x[7]],
tip_chord=1./1.1 root_m=0.01 root_p=0.4 tip_m=0.01 tip_p=0.4 t = 0.11 filename = 'wing.vtp' build_mesh(filename, semi_span, root_chord, tip_chord, root_m, root_p, tip_m, tip_p, t=t, nu=50, nv=41) Patm = 101325.#Pa rho = 1.225#kg.m-3 c = numpy.sqrt(1.4*Patm/rho) Mach = 0.15 airspeed = c*Mach mu=1.5e-5 Re = rho*airspeed*root_chord/mu print "Re=",Re print "airspeed=",airspeed*3.6