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)])
 
#==============================================================================
# Calculating the induced Velocity
#==============================================================================
V_external           = array([repeat(windspeed[0],shape(cylinder.collocationPoint)[1]),
                              repeat(windspeed[1],shape(cylinder.collocationPoint)[1])])
Example #2
0
# Creating multi-body
windturbine = multiBody(dict(body         = tower,
                             location     = [0.,0.],
                             global_pitch = 0.),
                             
                        dict(body         = airfoilA,
                             location     = [0.,2.],
                             global_pitch = 0.),  
                        
                        dict(body         = airfoilB,
                             location     = [0.,-2.], 
                             global_pitch = 180.))
'''                      
                                   
cylinder = multiBody(dict(body         = tower,
                          location     = [0., 0.],
                          global_pitch = -360./(2*n_panels))) 

#block   = multiBody(dict(body         = box,
#                         location     = [0., 0.],
#                         global_pitch = 0.))
                                                             
#==============================================================================
# Defining Vortex and scan points
#==============================================================================
# Vortex Field
#vort         = vortex([-1.],[0.],[1.])

# Plotting: Mesh grid
x,y = meshgrid(linspace(-10,10,200),linspace(-10,10,200))
mesh = array([concatenate(x),concatenate(y)])