def createPoints(self):
pointList = cylinder(150)
self.matCr = MatCreator(pointList,self.flows)
self.vortexPList = self.matCr.getMatVP()
self.vortexBlobPoints = []
self.vortexPanelPoints = []
for element in self.vortexPList:
vPanel = element.position*(1+element.length/math.pi/abs(element.position))
self.vortexBlobPoints.append(vPanel)
self.vortexPanelPoints.append(1.01*element.position)
self.currentFlowArray =[]
self.currentFlowArray.extend(self.vortexPList)
self.currentFlowArray.extend(self.flows)
self.sim.updateElements([],self.currentFlowArray)
def createPoints(self):
pointList = cylinder(150)
self.matCr = MatCreator(pointList, self.flows)
self.vortexPList = self.matCr.getMatVP()
self.vortexBlobPoints = []
self.vortexPanelPoints = []
for element in self.vortexPList:
vPanel = element.position * (
1 + element.length / math.pi / abs(element.position))
self.vortexBlobPoints.append(vPanel)
self.vortexPanelPoints.append(1.01 * element.position)
self.currentFlowArray = []
self.currentFlowArray.extend(self.vortexPList)
self.currentFlowArray.extend(self.flows)
self.sim.updateElements([], self.currentFlowArray)
import flows as f
from simulator import Simulator
import plotter as p
def addTracerPoint(pFlowArray):
for i in range(-20, 20, 5):
pos = complex(-2, i / 10.0)
pFlowArray.append(f.Tracer(pos, 1, 1, False))
if __name__ == '__main__':
print('Flow pass a circular cylinder')
pFlowArray = []
pFlowArray.append(f.Uniform(complex(-100, 0), 1, 1, True))
# pFlowArray.append(f.Vortex(complex(-1.010,0),1,1,False))
pointlist = cylinder(40)
matCr = MatCreator(pointlist, pFlowArray)
#addTracerPoint(pFlowArray)
length = len(pFlowArray)
pFlowArray.extend(matCr.getMatVP())
sim = Simulator(timestep=0.1)
sim.updateElements(pFlowArray)
data, colors = sim.run(0.2, matCr.updateFlows)
plotter = p.ParticlePlotter((-2, 2), (-2, 2))
plotter.update(data, colors)
filename = "temp"
#plotter.run(filename,True,False,False)
plotter.run(filename, False, True, False)
import flows as f
from simulator import Simulator
import plotter as p
def addTracerPoint(pFlowArray):
for i in range(-20, 20, 5):
pos = complex(-2, i / 10.0)
pFlowArray.append(f.Tracer(pos, 1, 1, False))
if __name__ == "__main__":
print("Flow pass a circular cylinder")
pFlowArray = []
pFlowArray.append(f.Uniform(complex(-100, 0), 1, 1, True))
# pFlowArray.append(f.Vortex(complex(-1.010,0),1,1,False))
pointlist = cylinder(40)
matCr = MatCreator(pointlist, pFlowArray)
# addTracerPoint(pFlowArray)
length = len(pFlowArray)
pFlowArray.extend(matCr.getMatVP())
sim = Simulator(timestep=0.1)
sim.updateElements(pFlowArray)
data, colors = sim.run(0.2, matCr.updateFlows)
plotter = p.ParticlePlotter((-2, 2), (-2, 2))
plotter.update(data, colors)
filename = "temp"
# plotter.run(filename,True,False,False)
plotter.run(filename, False, True, False)
