boundaryNodes = mesh.boundaryNodes neighborsNodes = mesh.neighborsNodes neighborsNodesALE = mesh.neighborsNodesALE neighborsElements = mesh.neighborsElements minLengthMesh = mesh.minLengthMesh FreedomDegree = mesh.FreedomDegree numPhysical = mesh.numPhysical Re = 1.0 Sc = 1.0 CFL = 0.5 dt = float(CFL*minLengthMesh) #dt = 0.1 #SL elif polynomial_option == 2: mesh = importMSH.Mini2D(pathMSHFile, mshFileName) numNodes = mesh.numNodes numVerts = mesh.numVerts numElements = mesh.numElements x = mesh.x y = mesh.y IEN = mesh.IEN boundaryEdges = mesh.boundaryEdges boundaryNodes = mesh.boundaryNodes neighborsNodes = mesh.neighborsNodes neighborsNodesALE = mesh.neighborsNodesALE neighborsNodesPressure = mesh.neighborsNodesPressure neighborsElements = mesh.neighborsElements minLengthMesh = mesh.minLengthMesh velocityFreedomDegree = mesh.velocityFreedomDegree
pathMSHFile = searchMSH.Find(mshFileName) if pathMSHFile == 'File not found': sys.exit() mesh = importMSH.Linear2D(pathMSHFile, mshFileName) # Mini Element elif polynomial_option == 2: mshFileName = 'linearHalfPoiseuille.msh' equation_number = 3 pathMSHFile = searchMSH.Find(mshFileName) if pathMSHFile == 'File not found': sys.exit() mesh = importMSH.Mini2D(pathMSHFile, mshFileName, equation_number) mesh.coord() mesh.ien() # Quad Element elif polynomial_option == 3: mshFileName = 'quadHalfPoiseuille_quad.msh' equation_number = 3 pathMSHFile = searchMSH.Find(mshFileName) if pathMSHFile == 'File not found': sys.exit() mesh = importMSH.Quad2D(pathMSHFile, mshFileName, equation_number) mesh.coord() mesh.ien()