#==============================================================================
# Initialization of the constants
#==============================================================================
# Control Parameters
windspeed = np.array([[1.], [0.]]) # x-dir and y-dir respectively
n_panels = 100
#==============================================================================
# Initialization of bodies
#==============================================================================
# Creating bodies
tower = body(name = 'tower',
shape = ('cylinder', n_panels),
chord = 1,
local_pitch = 0.,
pivot_point = [0.5, 0.])
# Creating multi-body
cylinder = multiBody(dict(body = tower,
location = [0., 0.],
global_pitch = 0.))
#==============================================================================
# Defining Vortex and scan points
#==============================================================================
# Plotting: Mesh grid
x,y = meshgrid(linspace(-1,1,100),linspace(-1,1,100))
mesh = array([concatenate(x),concatenate(y)])
#==============================================================================
# Initialization of the constants
#==============================================================================
# Control Parameters
windspeed = np.array([[1.], [0.]]) # x-dir and y-dir respectively
n_panels = 10
#==============================================================================
# Initialization of bodies
#==============================================================================
# Creating bodies
airfoil = body(name = 'airfoil',
shape = ('geometries/NACA0012.txt', False),
chord = 1.,
local_pitch = 0.,
pivot_point = [0.25, 0.])
tower = body(name = 'tower',
shape = ('cylinder', n_panels),
chord = 1.,#sqrt(1.+1.),
local_pitch = 0.,
pivot_point = [0.5, 0.])
box = body(name = 'block',
shape = ('block', n_panels),
chord = 1.,
local_pitch = 0.,
pivot_point = [0., 0.])