wing2d.py

#!/usr/bin/env python
# -*- coding: utf-8 -*-

def main():
  from bempy2d import geomwing, geomwing, srfmatbc, bcondvel, solvephi, circlefld, \
                      calcphifld, fldmatbc, fldmatbcv, calcvelfld, fieldgrid
  from numpy import array, zeros, savetxt, eye, linspace
  from pylab import plot, show, grid
  nelem = 50
  chord = 1.5
  thickness = 0.3
  u = array([-1.,0.])
  xysrf = zeros((nelem+1,2))
  xnode = geomwing(nelem,chord,thickness)
  xysrf[:nelem,:] = xnode[:,:2]
  xysrf[nelem,:] = xnode[0,:2]
  #print xysrf
  plot(xysrf[:,0],xysrf[:,1])
  
  xnode = geomwing(nelem,chord,thickness)
  xysrf[:nelem,:] = xnode[:,:2]
  xysrf[nelem,:] = xnode[0,:2]
  #print xysrf
  plot(xysrf[:,0],xysrf[:,1])
  #B, C = srfmatbc(xnode)
  #print C
  #print B
  #chisrf = bcondvel(xnode, u)
  #phisrf = solvephi(B,C,chisrf,nelem)
  #print phisrf
  
  grid()
  show()

if( __name__ == '__main__'):
  main()
  