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)