def createModels(self): if self.enablePlateModel == True: self.plateFields = models.PlateFields('plate') self.pmodel = models.PlateModelA(self.geomFields, self.plateFields, self.solidMeshes) self.dmodel = models.PlateDeformationModelA( self.geomFields, self.plateFields, self.solidMeshes) if self.enableElecModel == True: self.elecFields = models.ElectricFields('elec') self.emodel = models.ElectricModelA(self.geomFields, self.elecFields, self.fluidMeshesNew) if self.enableFlowModel == True: self.flowFields = models.FlowFields('flow') self.fmodel = models.FlowModelA(self.geomFields, self.flowFields, self.fluidMeshes) #currently contact model is implemented in ComputeForce.py #if self.enableContactModel == True: # self.contactFields = models.ContactFields('contact') # self.cmodel = models.ContactModelA(self.geomFields,self.contactFields,self.fluidMeshesNew) print 'models are created'
solver_time = zeros(1, dtype='d') solver_maxtime = zeros(1, dtype='d') solver_mintime = zeros(1, dtype='d') solver_start[0] = MPI.Wtime() geomFields = models.GeomFields('geom') fluidMetricsCalculator = models.MeshMetricsCalculatorA(geomFields, fluidMeshes) fluidMetricsCalculator.init() solidBoundaryMeshes = [m.extractBoundaryMesh() for m in solidMeshes] solidBoundaryMetricsCalculator = models.MeshMetricsCalculatorA( geomFields, solidBoundaryMeshes) solidBoundaryMetricsCalculator.init() flowFields = models.FlowFields('flow') fmodel = models.FlowModelA(geomFields, flowFields, fluidMeshes) fluidWalls = [3, 4] fluidInlet = [5] fluidOutlet = [6] circleWalls = [4] bcMap = fmodel.getBCMap() for id in fluidWalls: if id in bcMap: bc = bcMap[id] bc.bcType = 'NoSlipWall' for id in fluidInlet: