def read(self):
beamReader = FluentCase(self.beamCaseFile)
beamReader.read()
print 'beam mesh read in. Done!'
print '------------------------------------------------------------'
self.solidMeshes = beamReader.getMeshList()
self.geomFields = models.GeomFields('geom')
self.solidMetricsCalculator = models.MeshMetricsCalculatorA(
self.geomFields, self.solidMeshes)
self.solidMetricsCalculator.init()
fluidReader = FluentCase(self.fluidCaseFile)
fluidReader.read()
print 'fluid mesh read in. Done!'
print '------------------------------------------------------------'
self.fluidMeshes = fluidReader.getMeshList()
self.fluidMeshesNew = self.fluidMeshes
self.fluidMetricsCalculator = models.MeshMetricsCalculatorA(
self.geomFields, self.fluidMeshes)
self.fluidMetricsCalculator.init()
self.solidBoundaryMeshes = [
m.extrude(1, self.beam_thickness, True) for m in self.solidMeshes
]
self.solidBoundaryMetricsCalculator = models.MeshMetricsCalculatorA(
self.geomFields, self.solidBoundaryMeshes)
self.solidBoundaryMetricsCalculator.init()
def __init__(self, beamMesh, backgroundMesh, beamThickness, initialGap):
## Read in 2d Solid Mesh
self.beam_thickness = beamThickness
self.Gap = initialGap
beamReader = FluentCase(beamMesh)
beamReader.read()
print "read solid mesh"
self.solidMeshes = beamReader.getMeshList()
self.geomFields = models.GeomFields('geom')
self.solidMetricsCalculator = models.MeshMetricsCalculatorA(
self.geomFields, self.solidMeshes)
self.solidMetricsCalculator.init()
## Define plate and deformation model
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)
bcMap = self.pmodel.getBCMap()
## Apply a default Boundary Condition
#for i, bc in bcMap.iteritems():
#bc.bcType = 'SpecifiedTraction'
## Read in 3d Background Mesh
fluidReader = FluentCase(backgroundMesh)
fluidReader.read()
self.fluidMeshes = fluidReader.getMeshList()
self.fluidMetricsCalculator = models.MeshMetricsCalculatorA(
self.geomFields, self.fluidMeshes)
self.fluidMetricsCalculator.init()
## Define electric model
self.elecFields = models.ElectricFields('elec')
self.emodel = models.ElectricModelA(self.geomFields, self.elecFields,
self.fluidMeshes)
bcMap = self.emodel.getBCMap()
## Apply Default boundary conditions
for i, bc in bcMap.iteritems():
bc.bcType = "Symmetry"
self.solidBoundaryMeshes = [
m.extrude(1, beamThickness, True) for m in self.solidMeshes
]
self.solidBoundaryMetricsCalculator = models.MeshMetricsCalculatorA(
self.geomFields, self.solidBoundaryMeshes)
self.solidBoundaryMetricsCalculator.init()
def __init__(self, beamMesh, backgroundMesh, beamThickness, initialGap):
## Read in 2d Solid Mesh
self.beam_thickness = beamThickness
self.Gap = initialGap
beamReader = FluentCase(beamMesh)
beamReader.read()
print "read solid mesh"
self.solidMeshes = beamReader.getMeshList()
self.geomFields = models.GeomFields('geom')
self.solidMetricsCalculator = models.MeshMetricsCalculatorA(self.geomFields,self.solidMeshes)
self.solidMetricsCalculator.init()
## Define plate and deformation model
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)
bcMap = self.pmodel.getBCMap()
## Apply a default Boundary Condition
#for i, bc in bcMap.iteritems():
#bc.bcType = 'SpecifiedTraction'
## Read in 3d Background Mesh
fluidReader = FluentCase(backgroundMesh)
fluidReader.read();
self.fluidMeshes = fluidReader.getMeshList()
self.fluidMetricsCalculator = models.MeshMetricsCalculatorA(self.geomFields,self.fluidMeshes)
self.fluidMetricsCalculator.init()
## Define electric model
self.elecFields = models.ElectricFields('elec')
self.emodel = models.ElectricModelA(self.geomFields,self.elecFields,self.fluidMeshes)
bcMap = self.emodel.getBCMap()
## Apply Default boundary conditions
for i, bc in bcMap.iteritems():
bc.bcType = "Symmetry"
self.solidBoundaryMeshes = [m.extrude(1, beamThickness, True) for m in self.solidMeshes]
self.solidBoundaryMetricsCalculator = models.MeshMetricsCalculatorA(self.geomFields,self.solidBoundaryMeshes)
self.solidBoundaryMetricsCalculator.init()
def read(self):
beamReader = FluentCase(self.beamCaseFile)
beamReader.read()
print 'beam mesh read in. Done!'
print '------------------------------------------------------------'
self.solidMeshes = beamReader.getMeshList()
self.geomFields = models.GeomFields('geom')
self.solidMetricsCalculator = models.MeshMetricsCalculatorA(self.geomFields,self.solidMeshes)
self.solidMetricsCalculator.init()
fluidReader = FluentCase(self.fluidCaseFile)
fluidReader.read()
print 'fluid mesh read in. Done!'
print '------------------------------------------------------------'
self.fluidMeshes = fluidReader.getMeshList()
self.fluidMeshesNew = self.fluidMeshes
self.fluidMetricsCalculator = models.MeshMetricsCalculatorA(self.geomFields,self.fluidMeshes)
self.fluidMetricsCalculator.init()
self.solidBoundaryMeshes = [m.extrude(1, self.beam_thickness, True) for m in self.solidMeshes]
self.solidBoundaryMetricsCalculator = models.MeshMetricsCalculatorA(self.geomFields,self.solidBoundaryMeshes)
self.solidBoundaryMetricsCalculator.init()
fmodel.updateTime()
# change as needed
if __name__ == '__main__' and fileBase is None:
if len(sys.argv) != 2:
usage()
fileBase = sys.argv[1]
filename = fileBase + "fvm_coupling_hexa10x4x4_80x32x32.cas"
#filename = fileBase+"fvm_coupling_hexa10x4x4_120x48x48.cas"
#filename = fileBase+"fvm_coupling_hexa10x4x4_160x64x64.cas"
reader = FluentCase(filename)
print "reading %s", filename
#import debug
reader.read()
meshes = reader.getMeshList()
mesh0 = meshes[0]
import time
t0 = time.time()
geomFields = fvm.models.GeomFields('geom')
metricsCalculator = fvm.models.MeshMetricsCalculatorA(geomFields, meshes)
f.write(str(cellNodes[n](i, node) + 1) + " ")
f.write("\n")
f.write("\n")
f.close()
parser = OptionParser()
parser.set_defaults(type='quad')
parser.add_option("--type", help="'quad'[default], 'tri', 'hexa', or 'tetra'")
parser.add_option("--xdmf", action='store_true', help="Dump data in xdmf")
parser.add_option("--time",
"-t",
action='store_true',
help="Print timing information.")
(options, args) = parser.parse_args()
reader0 = FluentCase(sys.argv[1])
reader1 = FluentCase(sys.argv[2])
#import debug
reader0.read()
reader1.read()
meshes0 = reader0.getMeshList()
meshes1 = reader1.getMeshList()
#for mesh in meshes:
# mesh.getCells().clearGatherScatterMaps()
mesh0 = meshes0[0]
mesh1 = meshes1[0]
nmesh = 2
#fileBase = "/home/yildirim/memosa/src/fvm/test/cav32"
# change as needed
outfile = None
if __name__ == '__main__' and fileBase is None:
if len(sys.argv) < 2:
usage()
fileBase = sys.argv[1]
if len(sys.argv) == 3:
outfile = sys.argv[2]
if outfile == None:
outfile = fileBase + "-prism.dat"
reader = FluentCase(fileBase + ".cas")
#import debug
reader.read()
import sys
fluent_meshes = reader.getMeshList()
import time
t0 = time.time()
nmesh = MPI.COMM_WORLD.Get_size()
#print "nmesh = ", nmesh
#npart = fvmparallel.IntVector(1,nmesh) #total of distributed meshes
#etype = fvmparallel.IntVector(1,1) #triangle
#!/usr/bin/env python
import fvm
fvm.set_atype('double')
import fvm.fvmbaseExt as fvmbaseExt
import os
from mpi4py import MPI
from FluentCase import FluentCase
fileBase = os.getcwd()
coordFile = os.path.join(fileBase, "Grid_Coord.dat")
velocityFile = os.path.join(fileBase, "Grid_Velocity.dat")
infile = os.path.join(fileBase, "3D-cantilever.cas")
outfile = '/dev/stdout'
reader = FluentCase(infile)
reader.read();
meshes = reader.getMeshList()
mesh0 = meshes[0]
geomFields = fvm.models.GeomFields('geom')
metricsCalculator = fvm.models.MeshMetricsCalculatorA(geomFields,meshes)
metricsCalculator.init()
if fvm.atype == 'tangent':
metricsCalculator.setTangentCoords(0,7,1)
flowFields = fvm.models.FlowFields('flow')
fmodel = fvm.models.FlowModelA(geomFields,flowFields,meshes)
reader.importFlowBCs(fmodel)
import sys
sys.setdlopenflags(0x100|0x2)
import fvm.fvmbaseExt as fvmbaseExt
import fvm.importers as importers
import fvm.models_atyped_double as models
import fvm.exporters_atyped_double as exporters
from FluentCase import FluentCase
#fvmbaseExt.enableDebug("cdtor")
from mpi4py import MPI
reader = FluentCase(sys.argv[1])
reader.read();
meshes_case = reader.getMeshList()
# add double shell to mesh between two materials
# When creating double shell mesh, the mesh it is created from
# is called the 'parent' mesh and the mesh that is passed as an argument
# is called the 'other' mesh
#
# The phi values at the interface between the two meshes are related as follows:
# Phi_other = A * Phi_parent + B
interfaceID = 9
shellmesh = meshes_case[1].createDoubleShell(interfaceID, meshes_case[0], interfaceID)
meshes = [meshes_case[0], meshes_case[1], shellmesh]
geomFields = models.GeomFields('geom')
metricsCalculator = models.MeshMetricsCalculatorA(geomFields,meshes)
metricsCalculator.init()
timeStep = 10
numTimeSteps = 10
saveFrequency = 2
parser = optparse.OptionParser()
parser.add_option("--volt", type=float)
parser.add_option("--type", help="'tri'[default], 'quad', 'hexa', or 'tetra'")
parser.add_option("--time",
"-t",
action='store_true',
help="Print timing information.")
(options, args) = parser.parse_args()
print args[0]
print args[1]
fluidReader = FluentCase(args[0])
solidReader = FluentCase(args[1])
fluidReader.read()
solidReader.read()
fluidMeshes0 = fluidReader.getMeshList()
solidMeshes = solidReader.getMeshList()
nodeCoord = solidMeshes[0].getNodeCoordinates().asNumPyArray()
nodeCoord[:, :] *= 0.5
#paritioning
nmesh = 1
npart = [MPI.COMM_WORLD.Get_size()]
print "options.typeFluid = ", options.type
if ((i + 1) % 1000) == 0:
saveData(flowFields, reader, fileBase, i)
outfile = None
if __name__ == '__main__' and fileBase is None:
if len(sys.argv) < 2:
usage()
fileBase = sys.argv[1]
if len(sys.argv) == 3:
outfile = sys.argv[2]
if outfile == None:
outfile = fileBase + "-memosanew.dat"
reader = FluentCase(sys.argv[1])
#caseBase = fileBase + "Noslipbc_igc"
#reader = FluentCase(caseBase+".cas")
#import debug
reader.read()
meshes = reader.getMeshList()
import time
t0 = time.time()
geomFields = fvm.models.GeomFields('geom')
metricsCalculator = fvm.models.MeshMetricsCalculatorA(geomFields, meshes)
metricsCalculator.init()
import sys
sys.setdlopenflags(0x100 | 0x2)
import fvm.fvmbaseExt as fvmbaseExt
import fvm.importers as importers
import fvm.models_atyped_double as models
import fvm.exporters_atyped_double as exporters
from FluentCase import FluentCase
#fvmbaseExt.enableDebug("cdtor")
import time
reader = FluentCase("../test/FullBatterySimple.cas")
reader.read()
meshes_case = reader.getMeshList()
# add double shell to mesh between two materials
# When creating double shell mesh, the mesh it is created from
# is called the 'parent' mesh and the mesh that is passed as an argument
# is called the 'other' mesh
#
# parent has to be electrolyte, other has to be solid electrode
# for B-V equations to work
#interfaceID1 = 11
shellMeshAnode = meshes_case[2].createDoubleShell(6, meshes_case[1], 17)
#interfaceID2 = 17
shellMeshCathode = meshes_case[2].createDoubleShell(5, meshes_case[0], 16)
meshes = [
meshes_case[0], meshes_case[1], meshes_case[2], shellMeshAnode,
shellMeshCathode
writer.finish()
timeStep = 10
numTimeSteps = 10
saveFrequency = 2
parser = optparse.OptionParser()
parser.add_option("--volt", type=float)
parser.add_option("--type", help="'tri'[default], 'quad', 'hexa', or 'tetra'")
parser.add_option("--time","-t",action='store_true',help="Print timing information.")
(options, args) = parser.parse_args()
print args[0]
print args[1]
fluidReader = FluentCase(args[0])
solidReader = FluentCase(args[1])
fluidReader.read()
solidReader.read()
fluidMeshes0 = fluidReader.getMeshList()
solidMeshes = solidReader.getMeshList()
nodeCoord = solidMeshes[0].getNodeCoordinates().asNumPyArray()
nodeCoord[:,:] *=0.5
#paritioning
nmesh = 1
npart = [MPI.COMM_WORLD.Get_size()]
import sys
sys.setdlopenflags(0x100|0x2)
import fvmbaseExt
import importers
import models_atyped_double as models
import exporters_atyped_double as exporters
from FluentCase import FluentCase
#fvmbaseExt.enableDebug("cdtor")
reader = FluentCase("/home/sm/olda/data/cond.cas")
#import debug
reader.read();
meshes = reader.getMeshList()
geomFields = models.GeomFields('geom')
metricsCalculator = models.MeshMetricsCalculatorA(geomFields,meshes)
metricsCalculator.init()
thermalFields = models.ThermalFields('therm')
tmodel = models.ThermalModelA(geomFields,thermalFields,meshes)
bcmap = tmodel.getBCMap()
bc2 = bcmap[2]
reader.importThermalBCs(tmodel)
tmodel.printBCs()
yvela[i] = Ur * sin(alfa) * sin(beta) + Ualfa * cos(alfa) * sin(beta)
zvela[i] = Ur * cos(alfa) - Ualfa * sin(alfa)
#pdb.set_trace()
fx[fg.site] = xvel
fy[fg.site] = yvel
fz[fg.site] = zvel
return fx,fy,fz
outfile = fileBase+"-prism.dat"
reader = FluentCase(fileBase+"cube-125k-adapted-1.cas")
#import debug
reader.read();
meshes = reader.getMeshList()
import time
t0 = time.time()
geomFields = models.GeomFields('geom')
metricsCalculator = models.MeshMetricsCalculatorA(geomFields,meshes)
metricsCalculator.init()
if atype == 'tangent':
f.write( str(cellNodes[n](i,node)+1) + " ")
f.write("\n")
f.write("\n")
f.close()
parser = OptionParser()
parser.set_defaults(type='quad')
parser.add_option("--type", help="'quad'[default], 'tri', 'hexa', or 'tetra'")
parser.add_option("--xdmf", action='store_true', help="Dump data in xdmf")
parser.add_option("--time","-t",action='store_true',help="Print timing information.")
(options, args) = parser.parse_args()
nmesh = 1
fileBase0 = "plate_steady_"
reader0 = FluentCase(sys.argv[1])
#import debug
reader0.read();
meshes0 = reader0.getMeshList()
mesh0 = meshes0[0]
nmesh = 1
meshes = []
meshes.append(mesh0)
import time
t0 = time.time()
from numpy import *
from mpi4py import MPI
from FluentCase import FluentCase
numIterations = 10
def usage():
print __doc__
sys.exit(1)
if len(sys.argv) < 2:
usage()
casfile = sys.argv[1]
reader = FluentCase(casfile)
reader.read()
fluent_meshes = reader.getMeshList()
nmesh = MPI.COMM_WORLD.Get_size()
#print "nmesh = ", nmesh
#npart = fvmparallel.IntVector(1,nmesh) #total of distributed meshes
#etype = fvmparallel.IntVector(1,1) #triangle
npart = [nmesh]
etype = [2]
part_mesh = fvmparallel.PartMesh(fluent_meshes, npart, etype)
part_mesh.setWeightType(0)
part_mesh.setNumFlag(0)
import sys
sys.setdlopenflags(0x100|0x2)
import fvm.fvmbaseExt as fvmbaseExt
import fvm.importers as importers
import fvm.models_atyped_double as models
import fvm.exporters_atyped_double as exporters
from FluentCase import FluentCase
#fvmbaseExt.enableDebug("cdtor")
reader = FluentCase("../test/RadialTransfer.cas")
reader.read();
meshes_case = reader.getMeshList()
# add double shell to mesh between two materials
# When creating double shell mesh, the mesh it is created from
# is called the 'parent' mesh and the mesh that is passed as an argument
# is called the 'other' mesh
#
# parent has to be electrolyte, other has to be solid electrode
# for B-V equations to work
interfaceID = 4
shellmesh = meshes_case[1].createDoubleShell(interfaceID, meshes_case[0], interfaceID)
meshes = [meshes_case[0], meshes_case[1], shellmesh]
geomFields = models.GeomFields('geom')
metricsCalculator = models.MeshMetricsCalculatorA(geomFields,meshes)
metricsCalculator.init()
anchors = [4, 5, 6, 7]
numTimeSteps = 1
globalTime = 0
globalCount = 0
timeStep = 5e-8
saveFrequency = 50
initialTransient = False
probeIndex = 50
### ===================== mesh read ===============================================###
fileBase = "./"
### 2D plate mesh
beamReader = FluentCase(sys.argv[2])
beamReader.read();
solidMeshes = beamReader.getMeshList()
geomFields = models.GeomFields('geom')
solidMetricsCalculator = models.MeshMetricsCalculatorA(geomFields,solidMeshes)
solidMetricsCalculator.init()
### 3D fluid mesh
fluidReader = FluentCase(sys.argv[1])
fluidReader.read();
fluidMeshes = fluidReader.getMeshList()
fluidMetricsCalculator = models.MeshMetricsCalculatorA(geomFields,fluidMeshes)
fluidMetricsCalculator.init()
nodes = fluidMeshes[0].getNodes()
xn = fluidMeshes[0].getNodeCoordinates().asNumPyArray()
import sys
sys.setdlopenflags(0x100 | 0x2)
import fvm.fvmbaseExt as fvmbaseExt
import fvm.importers as importers
import fvm.models_atyped_double as models
import fvm.exporters_atyped_double as exporters
from FluentCase import FluentCase
#fvmbaseExt.enableDebug("cdtor")
reader = FluentCase("../test/SpeciesTest.cas")
#import debug
reader.read()
meshes = reader.getMeshList()
print "here"
geomFields = models.GeomFields('geom')
metricsCalculator = models.MeshMetricsCalculatorA(geomFields, meshes)
metricsCalculator.init()
nSpecies = 2
smodel = models.SpeciesModelA(geomFields, meshes, nSpecies)
#Species 1 boundary conditions
SN = 0
bcmap1 = smodel.getBCMap(SN)
vcmap1 = smodel.getVCMap(SN)
import sys
sys.setdlopenflags(0x100 | 0x2)
import fvm.fvmbaseExt as fvmbaseExt
import fvm.importers as importers
import fvm.models_atyped_double as models
import fvm.exporters_atyped_double as exporters
from FluentCase import FluentCase
#fvmbaseExt.enableDebug("cdtor")
reader = FluentCase("../test/TwoMaterialTest.cas")
#import debug
reader.read()
meshes = reader.getMeshList()
geomFields = models.GeomFields('geom')
metricsCalculator = models.MeshMetricsCalculatorA(geomFields, meshes)
metricsCalculator.init()
nSpecies = 1
# species conditions
smodel = models.SpeciesModelA(geomFields, meshes, nSpecies)
bcmap = smodel.getBCMap(0)
vcmap = smodel.getVCMap(0)
vcRightZone = vcmap[0]
vcLeftZone = vcmap[1]
# change as needed
outfile = None
if __name__ == '__main__' and fileBase is None:
if len(sys.argv) < 2:
usage()
fileBase = sys.argv[1]
if len(sys.argv) == 3:
outfile = sys.argv[2]
if outfile == None:
outfile = fileBase+"-prism.dat"
reader = FluentCase(fileBase+".cas")
#import debug
reader.read();
import sys
fluent_meshes = reader.getMeshList()
import time
t0 = time.time()
nmesh = MPI.COMM_WORLD.Get_size()
#print "nmesh = ", nmesh
#npart = fvmparallel.IntVector(1,nmesh) #total of distributed meshes
#etype = fvmparallel.IntVector(1,1) #triangle
numTimeSteps = 2 globalTime = 0 globalCount = 0 timeStep = 2e-8 saveFrequency = 1 initialTransient = False ### ===================== mesh read ===============================================### fileBase_input = "/home/yildirim/memosa/src/fvm/scripts/cantilever3D_coupling/" fileBase_output = "./" + str(int(-applied_voltage)) + "/" print fileBase_input + "fluid_3D_new.cas" ### 3D fluid mesh fluidReader = FluentCase(fileBase_input + "fluid_3D_new.cas") fluidReader.read() fluent_meshes_fluid = fluidReader.getMeshList() # paritioning nmesh = 1 npart = [MPI.COMM_WORLD.Get_size()] print "options.typeFluid = ", options.typeFluid etypeFluid = [etype[options.typeFluid]] # partMesh constructor and setTypes part_mesh_fluid = fvmparallel.MeshPartitioner(fluent_meshes_fluid, npart, etypeFluid) part_mesh_fluid.setWeightType(0) part_mesh_fluid.setNumFlag(0) # actions part_mesh_fluid.isDebug(0) part_mesh_fluid.partition() part_mesh_fluid.mesh()
levels = int(sys.argv[5])
T1 = float(sys.argv[6])
T2 = float(sys.argv[7])
dimension = 3
eVtoJoule = 1.60217646e-19
def usage():
print "Usage: %s filebase [outfilename]" % sys.argv[0]
print "Where filebase.cas is a Fluent case file."
print "Output will be in filebase-prism.dat if it is not specified."
sys.exit(1)
reader = FluentCase(fileBase + filename + extension)
reader.read()
meshes = reader.getMeshList()
geomFields = fvm.models.GeomFields('geom')
metricsCalculator = fvm.models.MeshMetricsCalculatorA(geomFields, meshes)
metricsCalculator.init()
ntheta = 10
nphi = 20
tau = 1e-2
vg = 1000
K_space = pa.KspaceA(fileBase + BZfile + BZdisc + ".txt", dimension, 1)
#K_space=pa.KspaceA(1,tau,vg,4e14,ntheta,nphi)
kArray = (K_space.getHollandConductivity((T1 + T2) / 2)).asNumPyArray()
kArray = kArray * eVtoJoule
f.write(str(cellNodes[n](i, node) + 1) + " ")
f.write("\n")
f.write("\n")
f.close()
parser = OptionParser()
parser.set_defaults(type='quad')
parser.add_option("--type", help="'quad'[default], 'tri', 'hexa', or 'tetra'")
parser.add_option("--xdmf", action='store_true', help="Dump data in xdmf")
parser.add_option("--time",
"-t",
action='store_true',
help="Print timing information.")
(options, args) = parser.parse_args()
reader0 = FluentCase(fileBase0 + ".cas")
reader1 = FluentCase(fileBase1 + ".cas")
#import debug
reader0.read()
reader1.read()
meshes0 = reader0.getMeshList()
meshes1 = reader1.getMeshList()
#for mesh in meshes:
# mesh.getCells().clearGatherScatterMaps()
mesh0 = meshes0[0]
mesh1 = meshes1[0]
nmesh = 2
print " "
def usage():
print "Usage: %s filebase [outfilename]" % sys.argv[0]
print "Where filebase.cas is a Fluent case file."
print "Output will be in filebase-prism.dat if it is not specified."
sys.exit(1)
############################################################
#Read the mesh, and extract the information into the MEMOSA#
#data structures #
############################################################
print "1"
reader = FluentCase(
"/home/brad/memosa/src/fvm/src/modules/phononbase/testmesh/Toy2.cas")
reader.read()
print "2"
meshes = reader.getMeshList()
print "3"
geomFields = fvm.models.GeomFields('geom')
metricsCalculator = fvm.models.MeshMetricsCalculatorA(geomFields, meshes)
metricsCalculator.init()
###################################################################
#Specify the group velocity and the scattering rate. Since #
#the only parameter that matters is the Knudsen number (=vg*tau/L)#
#they needn't be physical values, so long as the desired Knudsen #
#is obtained. It is also assumed that the Brillioun zone is #
#isotropic and the group velocity is the same direction as the #
#K-space vector. Because of this, to discretize k-space, you #
'quad' : 'FEQUADRILATERAL',
'tetra' : 'FETETRAHEDRON',
'hexa' : 'FEBRICK'
}
xtype = {
'tri' : 'Triangle',
'quad' : 'Quadrilateral',
'tetra' : 'Tetrahedron',
'hexa' : 'Hexahedron'
}
parser = OptionParser()
parser.set_defaults(type='tri')
parser.add_option("--type", help="'tri'[default], 'quad', 'hexa', or 'tetra'")
parser.add_option("--xdmf", action='store_true', help="Dump data in xdmf")
parser.add_option("--time","-t",action='store_true',help="Print timing information.")
(options, args) = parser.parse_args()
if len(args) != 1:
usage()
numIterations = 10
reader = FluentCase(args[0])
reader.read()
#import ddd
meshes_fluent = reader.getMeshList()
mesh_assembler = fvmbaseExt.MeshAssembler( meshes_fluent )
mesh_assembler.debug_faceCells()
outfile = sys.argv[2]
if outfile == None:
outfile = fileBase+"-tecplt.dat"
#import debug
parser = OptionParser()
parser.set_defaults(type='tri')
parser.add_option("--type", help="'tri'[default], 'quad', 'hexa', or 'tetra'")
parser.add_option("--xdmf", action='store_true', help="Dump data in xdmf")
parser.add_option("--time","-t",action='store_true',help="Print timing information.")
(options, args) = parser.parse_args()
if len(args) != 1:
usage()
reader = FluentCase(args[0])
#reader = FluentCase(fileBase+".cas")
reader.read();
t0 = time.time()
meshes_fluent = reader.getMeshList()
nmesh = 1
npart = [MPI.COMM_WORLD.Get_size()]
etype = [etype[options.type]]
if not MPI.COMM_WORLD.Get_rank():
print "parmesh is processing"
#Partitinoing of spatial mesh
#partMesh constructor and setTypes
part_mesh = fvmparallel.MeshPartitioner( meshes_fluent, npart, etype );
pd.timeStep = 1.e-8
pd.saveFrequency = 1
pd.fxSum = 0
pd.fySum = 0
pd.fzSum = 0
pd.initialTransient = False
pd.probeIndex = 671699
pd.probeFile = open(fileBase + "tipDisplacement-se.dat", "w")
pd.forceFile = open(fileBase + "beamForce-se.dat", "w")
pd.velocityFile = open(fileBase + "tipVelocity-se.dat", "w")
### read in meshes ###
beamFile = fileBase + 'Gen2_IBM_beam_12_1200.cas'
fluidFile = fileBase + 'Gen2_IBM_backgroundmesh_C2D2_wo_sub.cas'
fluidReader = FluentCase(sys.argv[1])
beamReader = FluentCase(sys.argv[2])
beamReader.read()
fluidReader.read()
pd.solidMeshes = beamReader.getMeshList()
pd.fluidMeshes = fluidReader.getMeshList()
pd.solidBoundaryMeshes = [m.extractBoundaryMesh() for m in pd.solidMeshes]
### geometry field ###
pd.geomFields = models.GeomFields('geom')
pd.solidMetricsCalculator = models.MeshMetricsCalculatorA(
pd.geomFields, pd.solidMeshes)
pd.fluidMetricsCalculator = models.MeshMetricsCalculatorA(
parser = OptionParser()
parser.set_defaults(type='tri')
parser.add_option("--type", help="'tri'[default], 'quad', 'hexa', or 'tetra'")
parser.add_option("--xdmf", action='store_true', help="Dump data in xdmf")
parser.add_option("--time",
"-t",
action='store_true',
help="Print timing information.")
(options, args) = parser.parse_args()
if len(args) != 1:
usage()
numIterations = 10
reader = FluentCase(args[0])
reader.read()
fluent_meshes = reader.getMeshList()
nmesh = 1
npart = [MPI.COMM_WORLD.Get_size()]
etype = [etype[options.type]]
if options.time:
# time profile for partmesh
part_mesh_time = zeros(1, dtype='d')
part_mesh_start = zeros(1, dtype='d')
part_mesh_end = zeros(1, dtype='d')
part_mesh_maxtime = zeros(1, dtype='d')
part_mesh_mintime = zeros(1, dtype='d')
part_mesh_start[0] = MPI.Wtime()
for i in range(0,ncell):
nnodes_per_cell = cellNodes[n].getCount(i)
for node in range(0,nnodes_per_cell):
f.write( str(cellNodes[n](i,node)+1) + " ")
f.write("\n")
f.write("\n")
f.close()
parser = OptionParser()
parser.set_defaults(type='quad')
parser.add_option("--type", help="'quad'[default], 'tri', 'hexa', or 'tetra'")
parser.add_option("--xdmf", action='store_true', help="Dump data in xdmf")
parser.add_option("--time","-t",action='store_true',help="Print timing information.")
(options, args) = parser.parse_args()
reader0 = FluentCase(fileBase0+".cas")
reader1 = FluentCase(fileBase1+".cas")
#import debug
reader0.read();
reader1.read();
meshes0 = reader0.getMeshList()
meshes1 = reader1.getMeshList()
#for mesh in meshes:
# mesh.getCells().clearGatherScatterMaps()
mesh0 = meshes0[0]
mesh1 = meshes1[0]
parser = OptionParser()
parser.set_defaults(type='tri')
parser.add_option("--type", help="'tri'[default], 'quad', 'hexa', or 'tetra'")
parser.add_option("--xdmf", action='store_true', help="Dump data in xdmf")
parser.add_option("--time","-t",action='store_true',help="Print timing information.")
(options, args) = parser.parse_args()
if len(args) != 1:
usage()
pd = ProblemDefinition()
pd.fileBase = "."
reader = FluentCase(args[0])
reader.read();
pd.meshes = reader.getMeshList()
nmesh = 1
npart = [MPI.COMM_WORLD.Get_size()]
etype = [etype[options.type]]
if not MPI.COMM_WORLD.Get_rank():
print "parmesh is processing"
if options.time:
# time profile for partmesh
part_mesh_time = zeros(1,dtype='d')
part_mesh_start = zeros(1, dtype='d')
numTimeSteps = 2 globalTime = 0 globalCount = 0 timeStep = 2e-8 saveFrequency = 1 initialTransient = False probeIndex = 0 ### ===================== mesh read ===============================================### fileBase_input = "/home/yildirim/memosa/src/fvm/scripts/cantilever3D_coupling/" fileBase_output = "./" + str(int(-applied_voltage)) + "/" ### 2D plate mesh beamReader = FluentCase(fileBase_input+"beam_2D.cas") beamReader.read() ##paritioning #nmesh = 1 #npart = [MPI.COMM_WORLD.Get_size()] #print "options folud.type = ", options.type #etype = [etype[options.type]] ##partMesh constructor and setTypes #part_mesh = fvmparallel.MeshPartitioner( fluentMeshes, npart, etype ); #part_mesh.setWeightType(0); #part_mesh.setNumFlag(0); ##actions #part_mesh.isDebug(0) #part_mesh.partition() #part_mesh.mesh()
pd.timeStep = 1.e-8
pd.saveFrequency = 1
pd.fxSum = 0
pd.fySum = 0
pd.fzSum = 0
pd.initialTransient=False
pd.probeIndex = 671699
pd.probeFile = open(fileBase + "tipDisplacement-se.dat", "w")
pd.forceFile = open(fileBase + "beamForce-se.dat", "w")
pd.velocityFile = open(fileBase + "tipVelocity-se.dat", "w")
### read in meshes ###
beamFile = fileBase + 'Gen2_IBM_beam_12_1200.cas'
fluidFile = fileBase + 'Gen2_IBM_backgroundmesh_C2D2_wo_sub.cas'
fluidReader = FluentCase(sys.argv[1])
beamReader = FluentCase(sys.argv[2])
beamReader.read()
fluidReader.read()
pd.solidMeshes = beamReader.getMeshList()
pd.fluidMeshes = fluidReader.getMeshList()
pd.solidBoundaryMeshes = [m.extractBoundaryMesh() for m in pd.solidMeshes]
### geometry field ###
pd.geomFields = models.GeomFields('geom')
pd.solidMetricsCalculator = models.MeshMetricsCalculatorA(pd.geomFields,
pd.solidMeshes)
pd.fluidMetricsCalculator = models.MeshMetricsCalculatorA(pd.geomFields,
help="Print timing information.")
(options, args) = parser.parse_args()
outfile = None
fileBase = None
if __name__ == '__main__' and fileBase is None:
if len(sys.argv) < 2:
usage()
fileBase = sys.argv[1]
if len(sys.argv) == 3:
outfile = sys.argv[2]
if outfile == None:
outfile = fileBase + "-prism.dat"
reader = FluentCase(fileBase)
#import debug
reader.read()
meshes = reader.getMeshList()
mesh = meshes[0]
nmesh = 1
print outfile
import time
t0 = time.time()
geomFields = models.GeomFields('geom')
metricsCalculator = models.MeshMetricsCalculatorA(geomFields, meshes)
metricsCalculator.init()
numIterations =50
fileBase = "/home/lin/work/app-memosa/src/fvm/verification/IBM_convergence/"
def advance(tmodel,particles,niter):
for i in range(0,niter):
try:
tmodel.computeIBFaceTemperature(particles)
tmodel.advance(1)
except KeyboardInterrupt:
break
outfile = fileBase+"-prism.dat"
reader = FluentCase(fileBase+"cav64_TScalar.cas")
mpmFileName = fileBase + "MPMs_test1.dat"
#import ddd
reader.read();
meshes = reader.getMeshList()
import time
t0 = time.time()
geomFields = models.GeomFields('geom')
metricsCalculator = models.MeshMetricsCalculatorA(geomFields,meshes)
f.write("</Grid>\n")
f.close()
parser = OptionParser()
parser.set_defaults(type='tri')
parser.add_option("--type", help="'tri'[default], 'quad', 'hexa', or 'tetra'")
parser.add_option("--xdmf", action='store_true', help="Dump data in xdmf")
parser.add_option("--time","-t",action='store_true',help="Print timing information.")
(options, args) = parser.parse_args()
if len(args) != 1:
usage()
numIterations = 10
reader = FluentCase(args[0])
reader.read()
fluent_meshes = reader.getMeshList()
nmesh = 1
npart = [MPI.COMM_WORLD.Get_size()]
etype = [etype[options.type]]
if options.time:
# time profile for partmesh
part_mesh_time = zeros(1,dtype='d')
part_mesh_start = zeros(1, dtype='d')
part_mesh_end = zeros(1, dtype='d')
part_mesh_maxtime = zeros(1,dtype='d')
part_mesh_mintime = zeros(1, dtype='d')
part_mesh_start[0] = MPI.Wtime()
pd.saveFrequency = 2 pd.fxSum = 0 pd.fySum = 0 pd.fzSum = 0 pd.initialTransient=False #pd.probeIndex = 7106 #tip location at proc2 for 3processors run ### read in meshes ### #beamFile = fileBase + 'frogleg_device_252k.cas' #fluidFile = fileBase + 'frogleg_fluid_300k.cas' beamFile = fileBase + 'beam_500x10.cas' fluidFile = fileBase + 'fluid_500x40.cas' beamReader = FluentCase(beamFile) beamReader.read() #all fluid process will be aware of whole meshes pd.solidMeshes = beamReader.getMeshList() #partitioning fluidReader = FluentCase(fluidFile) fluidReader.read() nmesh = 1 npart = [MPI.COMM_WORLD.Get_size()] etype = [etype['quad']] #partMesh constructor and setTypes fluent_meshes = fluidReader.getMeshList() part_mesh = fvmparallel.MeshPartitioner( fluent_meshes, npart, etype ); part_mesh.setWeightType(0); part_mesh.setNumFlag(0);
numIterations =20
fileBase = "/home/lin/work/app-memosa/src/fvm/verification/ElectroStatics/"
def advance(elec_model,niter):
for i in range(0,niter):
try:
stopFlag=elec_model.advance(1)
if stopFlag == 1:
break
except KeyboardInterrupt:
break
reader = FluentCase(fileBase + "cav32_elec.cas")
reader.read()
#pdb.set_trace()
meshes = reader.getMeshList()
mesh0 = meshes[0]
import time
t0 = time.time()
geomFields = models.GeomFields('geom')
metricsCalculator = models.MeshMetricsCalculatorA(geomFields,meshes)
print "Output will be in filebase-cellmark.dat if it is not specified."
sys.exit(1)
outfile = None
if __name__ == '__main__' and fileBase is None:
if len(sys.argv) < 2:
usage()
fileBase = sys.argv[1]
if len(sys.argv) == 3:
outfile = sys.argv[2]
if outfile == None:
outfile = fileBase + "output.dat"
reader = FluentCase(fileBase + "cav32-new.cas")
reader.read()
meshes = reader.getMeshList()
mesh0 = meshes[0]
geomFields = models.GeomFields('geom')
metricsCalculator = models.MeshMetricsCalculatorA(geomFields, meshes)
metricsCalculator.init()
if fvm.atype == 'tangent':
metricsCalculator.setTangentCoords(0, 7, 1)
# change as needed
outfile = None
if __name__ == '__main__' and fileBase is None:
if len(sys.argv) < 2:
usage()
fileBase = sys.argv[1]
if len(sys.argv) == 3:
outfile = sys.argv[2]
if outfile == None:
outfile = fileBase + "-prism.dat"
reader = FluentCase(fileBase + ".cas")
#import debug
reader.read()
meshes = reader.getMeshList()
import time
t0 = time.time()
geomFields = models.GeomFields('geom')
metricsCalculator = models.MeshMetricsCalculatorA(geomFields, meshes)
metricsCalculator.init()
if atype == 'tangent':
nnodes_per_cell = cellNodes[n].getCount(i)
for node in range(0,nnodes_per_cell):
f.write( str(cellNodes[n](i,node)+1) + " ")
f.write("\n")
f.write("\n")
f.close()
# change as needed
#import debug
reader = FluentCase(sys.argv[1])
reader.read()
fluent_meshes = reader.getMeshList()
nmesh = 1
import time
t0 = time.time()
#print "nmesh = ", nmesh
#npart = fvmparallel.IntVector(1,nmesh) #total of distributed meshes
#etype = fvmparallel.IntVector(1,1) #triangle
npart = [MPI.COMM_WORLD.Get_size()]
numTimeSteps = 2
globalTime = 0
globalCount = 0
timeStep = 2e-8
saveFrequency = 1
initialTransient = False
### ===================== mesh read ===============================================###
fileBase_input = "/home/yildirim/memosa/src/fvm/scripts/cantilever3D_coupling/"
fileBase_output = "./" + str(int(-applied_voltage)) + "/"
print fileBase_input + "fluid_3D_new.cas"
### 3D fluid mesh
fluidReader = FluentCase(fileBase_input + "fluid_3D_new.cas")
fluidReader.read()
fluent_meshes_fluid = fluidReader.getMeshList()
#paritioning
nmesh = 1
npart = [MPI.COMM_WORLD.Get_size()]
print "options.typeFluid = ", options.typeFluid
etypeFluid = [etype[options.typeFluid]]
#partMesh constructor and setTypes
part_mesh_fluid = fvmparallel.MeshPartitioner(fluent_meshes_fluid, npart,
etypeFluid)
part_mesh_fluid.setWeightType(0)
part_mesh_fluid.setNumFlag(0)
#actions
part_mesh_fluid.isDebug(0)
part_mesh_fluid.partition()
f.write( str(cellNodes[n](i,node)+1) + " ")
f.write("\n")
f.write("\n")
f.close()
parser = OptionParser()
parser.set_defaults(type='quad')
parser.add_option("--type", help="'quad'[default], 'tri', 'hexa', or 'tetra'")
parser.add_option("--xdmf", action='store_true', help="Dump data in xdmf")
parser.add_option("--time","-t",action='store_true',help="Print timing information.")
(options, args) = parser.parse_args()
nmesh = 1
fileBase0 = 'plate_creep_deformation_'
reader0 = FluentCase(sys.argv[1])
#import debug
reader0.read();
meshes0 = reader0.getMeshList()
mesh0 = meshes0[0]
nmesh = 1
meshes = []
meshes.append(mesh0)
import time
t0 = time.time()
fmodel.updateTime()
# change as needed
if __name__ == '__main__' and fileBase is None:
if len(sys.argv) != 2:
usage()
fileBase = sys.argv[1]
filename = fileBase+"fvm_coupling_hexa10x4x4_80x32x32.cas"
#filename = fileBase+"fvm_coupling_hexa10x4x4_120x48x48.cas"
#filename = fileBase+"fvm_coupling_hexa10x4x4_160x64x64.cas"
reader = FluentCase(filename)
print "reading %s", filename
#import debug
reader.read();
meshes = reader.getMeshList()
mesh0 = meshes[0]
import time
t0 = time.time()
geomFields = fvm.models.GeomFields('geom')
metricsCalculator = fvm.models.MeshMetricsCalculatorA(geomFields,meshes)
#!/usr/bin/env python
import fvm
fvm.set_atype('double')
import fvm.fvmbaseExt as fvmbaseExt
import fvm.importers as importers
from FluentCase import FluentCase
from mpi4py import MPI
from fvm.fvmbaseExt import vectorInt
from fvm.fvmbaseExt import VecD3
reader = FluentCase("cav32.cas")
reader.read()
meshes = reader.getMeshList()
geomFields = fvm.models.GeomFields('geom')
metricsCalculator = fvm.models.MeshMetricsCalculatorA(geomFields, meshes)
metricsCalculator.init()
flowFields = fvm.models.FlowFields('flow')
fmodel = fvm.models.FlowModelA(geomFields, flowFields, meshes)
reader.importFlowBCs(fmodel)
octree = fvmbaseExt.Octree()
mesh0 = meshes[0]
cells = mesh0.getCells()
nCells = cells.getCount()
cellCentroid = geomFields.coordinate[cells].asNumPyArray()
"""
cellCoord = open(fileBase+"cellCentroid.dat", "w")
for c in range(0, nCells):
cellCoord.write("%lf\t%lf\t%lf\n" % (cellCentroid[c][0], cellCentroid[c][1], cellCentroid[c][2]))
cellCoord.close()
"""
