def saveData(flowFields,reader,fileBase,i):
writer = exporters.FluentDataExporterA(reader,fileBase+str(i+1)+"-noslipnew.dat",False,0)
writer.init()
writer.writeScalarField(flowFields.pressure,1)
writer.writeScalarField(flowFields.density,101)
writer.writeVectorField(flowFields.velocity,111)
writer.writeScalarField(flowFields.massFlux,18)
writer.finish()
def saveData(flowFields, reader, fileBase):
writer = exporters.FluentDataExporterA(reader, fileBase + "-prism.dat",
False, 0)
writer.init()
writer.writeScalarField(flowFields.pressure, 1)
writer.writeVectorField(flowFields.velocity, 111)
writer.writeScalarField(flowFields.massFlux, 18)
writer.finish()
#advance(fmodel,numIterations)
t1 = time.time()
if outfile != '/dev/stdout':
print '\nsolution time = %f' % (t1-t0)
print 'solution time = %f' % (t1-t0)
print '\n\npressure integrals\n'
fmodel.printPressureIntegrals()
print '\n\nmomentum flux integrals\n'
fmodel.printMomentumFluxIntegrals()
writer = exporters.FluentDataExporterA(reader,fileBase+"-ibm.dat",False,0)
writer.init()
writer.writeScalarField(flowFields.pressure,1)
writer.writeVectorField(flowFields.velocity,111)
writer.writeScalarField(flowFields.massFlux,18)
writer.finish()
if (atype=='tangent'):
writer = exporters.FluentDataExporterA(reader,fileBase+"-ibm-tangent.dat",False,1)
writer.init()
writer.writeScalarField(flowFields.pressure,1)
writer.writeVectorField(flowFields.velocity,111)
writer.writeScalarField(flowFields.massFlux,18)
writer.finish()
particles = solid.getParticles()
px = solid.getCoordinates()
pv = solid.getVelocities()
geomFields.coordinate[particles] = px
flowFields.velocity[particles] = pv
metricsCalculator.computeIBInterpolationMatrices(particles)
metricsCalculator.computeIBandSolidInterpolationMatrices(particles)
advanceUnsteady(fmodel, meshes, globalTime, numTimeSteps, particles)
t1 = time.time()
print 'solution time = %f' % (t1 - t0)
writer = exporters.FluentDataExporterA(reader,
fileBase + "P=83593-fullbeam.dat",
False, 0)
writer.init()
writer.writeScalarField(flowFields.pressure, 1)
writer.writeVectorField(flowFields.velocity, 111)
writer.writeScalarField(flowFields.massFlux, 18)
writer.finish()
if (atype == 'tangent'):
writer = exporters.FluentDataExporterA(reader,
fileBase + "-prism-tangent.dat",
False, 1)
writer.init()
writer.writeScalarField(flowFields.pressure, 1)
writer.writeVectorField(flowFields.velocity, 111)
writer.writeScalarField(flowFields.massFlux, 18)
bc['specifiedXVelocity']=fx
bc['specifiedYVelocity']=fy
bc['specifiedZVelocity']=fz
bc.bcType = 'VelocityBoundary'
fmodel.init()
initVelocity(geomFields,flowFields.velocity,meshes)
fmodel.computeIBFaceVelocity(particles)
advance(fmodel,particles,numIterations)
t1 = time.time()
if outfile != '/dev/stdout':
print '\nsolution time = %f' % (t1-t0)
writer = exporters.FluentDataExporterA(reader,fileBase+"cube-125k-adapted-withconvection.dat",False,0)
writer.init()
writer.writeScalarField(flowFields.pressure,1)
writer.writeVectorField(flowFields.velocity,111)
writer.writeScalarField(flowFields.massFlux,18)
writer.finish()
if (atype=='tangent'):
writer = exporters.FluentDataExporterA(reader,fileBase+"-prism-tangent.dat",False,1)
writer.init()
writer.writeScalarField(flowFields.pressure,1)
writer.writeVectorField(flowFields.velocity,111)
writer.writeScalarField(flowFields.massFlux,18)
writer.finish()
vc.setVar('viscosity',(1.7894e-5,1))
"""
#import debug
fmodel.init()
dmodel.init()
initIterations=0
#advancenew([fmodel],dmodel,0,initIterations)
advance([fmodel],dmodel,initIterations,numIterations)
#advance([fmodel],dmodel,numIterations)-old
t1 = time.time()
if outfile != '/dev/stdout':
print '\nsolution time = %f' % (t1-t0)
writer = exporters.FluentDataExporterA(reader,fileBase+"-remove.dat",False,0)
writer.init()
writer.writeScalarField(flowFields.pressure,1)
writer.writeScalarField(flowFields.density,101)
writer.writeVectorField(flowFields.velocity,111)
writer.writeScalarField(flowFields.massFlux,18)
writer.finish()
if (atype=='tangent'):
writer = exporters.FluentDataExporterA(reader,fileBase+"-prism-tangent.dat",False,1)
writer.init()
writer.writeScalarField(flowFields.pressure,1)
writer.writeVectorField(flowFields.velocity,111)
writer.writeScalarField(flowFields.massFlux,18)
writer.finish()
#bc.setVar('specifiedPotentialFlux',0)
#bc.bcType = 'SpecifiedCharge'
#bc.setVar('specifiedCharge',350)
elec_options = elec_model.getOptions()
elec_options.electrostatics = 1
elec_options.chargetransport = 0
elec_options.tunneling = 1
elec_options.ibm = 0
elec_model.printBCs()
import ddd
elec_model.init()
#set up permittivity
cells = meshes[0].getCells()
perm = electronicFields.permittivity[cells].asNumPyArray()
perm[:] = 1.0
advanceUnsteady(elec_model, globalTime, numTimeSteps)
t1 = time.time()
print 'solution time = %f' % (t1 - t0)
writer = exporters.FluentDataExporterA(reader, fileBase + "cav32_test.dat",
False, 0)
writer.init()
writer.writeScalarField(electronicFields.charge, 3)
writer.finish()
cells = mesh0.getCells() nCells = cells.getCount() # add particles position and velocities to the main fields particles = solid.getParticles() px = solid.getCoordinates() pt = solid.getTemperatures() geomFields.coordinate[particles]=px thermalFields.temperature[particles]=pt metricsCalculator.computeIBInterpolationMatrices(particles) tmodel.printBCs() tmodel.init() advance(tmodel,particles,numIterations) t1 = time.time() writer = exporters.FluentDataExporterA(reader,fileBase+"cav8_TScalar.dat",False,0) writer.init() writer.writeScalarField(thermalFields.temperature,3) writer.finish()
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()
tmodel.init()
tmodel.advance(1)
writer = exporters.FluentDataExporterA(reader, "test.dat", False, 0)
writer.init()
writer.writeScalarField(thermalFields.temperature, 3)
writer.finish()
