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()
