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==================================================