nCellsZI = 40 #number of cells in inlet z direction # #----------------------------------------------------------------------- # mesh grading - not implemented #~ grDX, grDY, grDZ = 1.0, 1.0, 1.0 #grading in dense block #~ grSX, grSY, grSZ = 1.0, 1.0, 1.0 #gradins in sparse block # scale (to meters) mScale = 1 #----------------------------------------------------------------------- # FUNCTION CALL #----------------------------------------------------------------------- print 'WRITING BLOCKMESHDICT======================\n\n' fblockMeshGen(caseDir, x0, y0, z0, cL, cW, cH, iH, tL, dbW, nDiv, d0, nCellsXDL, nCellsXDS, nCellsYD, nCellsZD, nCellsXSL, nCellsXSS, nCellsYS, nCellsZS, nCellsZI, mScale) print 'DONE=======================================\n\n' #BC FILES MODIFICATION================================================== #----------------------------------------------------------------------- # FUNCTION CALL #----------------------------------------------------------------------- print 'ADJUSTING BC===============================\n\n' README.write('\n 0.org==============================================\n') #~ finletBCWriter(a0,h0,eps,aaRat,beta,rho,g,mu,sigma,nCellsW,nCellsL) #----------------------------------------------------------------------- # U #----------------------------------------------------------------------- # # Boundary conditions for the velocity field
README.write('nCores \t = \t ' + repr(nCores) + '\n') README.write('startTime \t = \t ' + repr(startTime) + ' s\n') README.write('endTime \t = \t ' + repr(endTime) + ' s\n') #BLOCKMESHDICT FILE GENERATION========================================== # -- prepara variables geomSize= [aG,lG,hG] #variable with geometry dimensions cellSize= [dA,dL,dH] #variable with cell dimensions # -- additional parameters spGrad = 5 #grading intensity in zDir of sparse block mScale = 1 #conversion to metres #----------------------------------------------------------------------- # FUNCTION CALL #----------------------------------------------------------------------- print 'WRITING BLOCKMESHDICT======================\n\n' fblockMeshGen(caseDir,h0,geomSize,cellSize, spGrad,mScale) print 'DONE=======================================\n\n' #BC FILES MODIFICATION================================================== #----------------------------------------------------------------------- # FUNCTION CALL #----------------------------------------------------------------------- print 'ADJUSTING BC===============================\n\n' README.write('\n 0.org==============================================\n') #----------------------------------------------------------------------- # alpha.liquid #----------------------------------------------------------------------- # # Boundary conditions for the velocity field #
# ----------------------------------------------------------------------- # FUNCTION CALL # ----------------------------------------------------------------------- print "WRITING BLOCKMESHDICT======================\n\n" fblockMeshGen( caseDir, x0, y0, z0, cL, cW, cH, iH, tL, dbW, nDiv, d0, nCellsXDL, nCellsXDS, nCellsYD, nCellsZD, nCellsXSL, nCellsXSS, nCellsYS, nCellsZS, nCellsZI, mScale, ) print "DONE=======================================\n\n" # BC FILES MODIFICATION==================================================