# Absorption spectrum
aabs = Property.Property(value=ex_em_import.getAbs(),
propID=PropertyID.PID_AsorptionSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(aabs)
hmri = Property.Property(1.55,
PropertyID.PID_RefractiveIndex,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_CONE)
mieApp.setProperty(hmri)
logger.info('Properties set!')
# Solve
mieApp.solveStep(0)
tracerApp.solveStep(0, runInBackground=False)
comsolApp.solveStep(0)
# Plot data to file
logger.info("Saving vtk")
v = fTemp.field2VTKData()
v.tofile('testTemperature.vtk')
v = fHeat.field2VTKData()
v.tofile('testHeat.vtk')
# Absorption spectrum
aabs = Property.Property(value=ex_em_import.getAbs(),
propID=PropertyID.PID_AsorptionSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(aabs)
hmri = Property.Property(1.55, PropertyID.PID_RefractiveIndex,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_CONE)
mieApp.setProperty(hmri)
logger.info('Properties set!')
# Solve
mieApp.solveStep(0)
tracerApp.solveStep(0, runInBackground=False)
comsolApp.solveStep(0)
# Plot data to file
logger.info("Saving vtk")
v = fTemp.field2VTKData()
v.tofile('testTemperature.vtk')
v = fHeat.field2VTKData()
v.tofile('testHeat.vtk')
c = {}
C = np.loadtxt('../../../mmp_tracer_api/data/Abs_GREEN.dat')
c['wavelengths'] = C[:, 0]
c['intensities'] = C[:, 1]
aabs = Property.Property(value=c,
propID=PropertyID.PID_AsorptionSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(aabs)
# Solve
times = np.arange(0, 3)
for t in times:
mieApp.solveStep(t)
tracerApp.solveStep(t, runInBackground=False)
comsolApp.solveStep(t)
fHeat = tracerApp.getField(FieldID.FID_Thermal_absorption_volume, 1.5)
pDens = tracerApp.getProperty(propID=PropertyID.PID_ParticleNumberDensity,
objectID=objID.OBJ_CONE,
time=1.5)
# print('%.10f' % pDens.value)
# Plot data to file
logger.info("Saving vtk")
v = fHeat.field2VTKData()
v.tofile('testHeat.vtk')
C = np.loadtxt('../../../mmp_tracer_api/data/Abs_GREEN.dat')
c['wavelengths'] = C[:, 0]
c['intensities'] = C[:, 1]
aabs = Property.Property(value=c,
propID=PropertyID.PID_AsorptionSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(aabs)
# Solve
times = np.arange(0, 3)
for t in times:
mieApp.solveStep(t)
tracerApp.solveStep(t, runInBackground=False)
comsolApp.solveStep(t)
fHeat = tracerApp.getField(FieldID.FID_Thermal_absorption_volume, 1.5)
pDens = tracerApp.getProperty(propID=PropertyID.PID_ParticleNumberDensity,
objectID=objID.OBJ_CONE, time=1.5)
# print('%.10f' % pDens.value)
# Plot data to file
logger.info("Saving vtk")
v = fHeat.field2VTKData()
v.tofile('testHeat.vtk')
