time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(em)
# Excitation spectrum
ex = Property.Property(value=ex_em_import.getEx(),
propID=PropertyID.PID_ExcitationSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(ex)
# 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)
# Solve
logger.info('Entering mieApp solver')
mieApp.solveStep(0)
logger.info('mieApp solved')
# Connect properties
logger.info('Connecting Mie properties...')
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(em)
# Excitation spectrum
ex = Property.Property(value=ex_em_import.getEx(),
propID=PropertyID.PID_ExcitationSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(ex)
# 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)
def xsim(z, *params):
#create the apps:
mieApp = MMPMie('localhost')
tracerApp = MMPRaytracer('localhost')
comsolApp = MMPComsolDummy('localhost')
# Point data conversion to false. Speeds up testing
tracerApp._convertPointData = False
# Set default LED json
tracerApp.setDefaultInputFile('./DefaultLED3.json')
#Set refractive indexes to mieApp:
pri = Property.Property(1.83,
PropertyID.PID_RefractiveIndex,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1)
mieApp.setProperty(pri)
# same for particle_type_2
pri2 = Property.Property(1.84,
PropertyID.PID_RefractiveIndex,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_2)
mieApp.setProperty(pri2)
hmri = Property.Property(1.55, PropertyID.PID_RefractiveIndex,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_CONE)
mieApp.setProperty(hmri)
#set also another particle type to MieApp:
pScat2 = Property.Property(0,
PropertyID.PID_ScatteringCrossSections,
valueType=ValueType.Vector,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_2)
pPhase2 = Property.Property(0,
propID=PropertyID.PID_InverseCumulativeDist,
valueType=ValueType.Vector,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_2)
mieApp.setProperty(pScat2)
mieApp.setProperty(pPhase2)
# Connect fields
fTemp = comsolApp.getField(FieldID.FID_Temperature, 0)
fHeat = comsolApp.getField(FieldID.FID_Thermal_absorption_volume, 0)
tracerApp.setField(fTemp)
tracerApp.setField(fHeat)
#set the constant properties from *params:
p_max = 35.0
p_min = 1.0
p_num = 50
w_max = 1100.0
w_min = 100.0
w_num = 1001
# Particle density
dens_p = 5.0 # g/cm3
# Silicone density
dens_host = 1.1 # g/cm3
# log mean in microns
mu1 = 2.4849
# log standard deviation in microns
sigma1 = 0.3878
mu2 = 2.35137
sigma2 = 0.627218
mu = [mu1, mu2]
sigma = [sigma1, sigma2]
mup = Property.Property(value=mu1, propID=PropertyID.PID_ParticleMu, valueType=ValueType.Scalar, time=0.0, units=None, objectID=objID.OBJ_PARTICLE_TYPE_1)
sigmap = Property.Property(value=sigma1, propID=PropertyID.PID_ParticleSigma, valueType=ValueType.Scalar, time=0.0, units=None, objectID=objID.OBJ_PARTICLE_TYPE_1)
mup2 = Property.Property(value=mu2, propID=PropertyID.PID_ParticleMu, valueType=ValueType.Scalar, time=0.0, units=None, objectID=objID.OBJ_PARTICLE_TYPE_2)
sigmap2 = Property.Property(value=sigma2, propID=PropertyID.PID_ParticleSigma, valueType=ValueType.Scalar, time=0.0, units=None, objectID=objID.OBJ_PARTICLE_TYPE_2)
mieApp.setProperty(mup)
mieApp.setProperty(sigmap)
mieApp.setProperty(mup2)
mieApp.setProperty(sigmap2)
# Particles diameters in micrometers
d = np.linspace(p_min, p_max, p_num)
#pdf = lognorm(sigma, scale=np.exp(mu))
#set the changing properties from z:
# Weight fractions
weight_frac = z / 100.0
print("z, weight_frac = ", z, weight_frac)
#run the simulation:
# Number of rays to trace
pRays = Property.Property(value=1000000,
propID=PropertyID.PID_NumberOfRays,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_CONE)
tracerApp.setProperty(pRays)
n_particles = Property.Property(value=2,
propID=PropertyID.PID_NumberOfFluorescentParticles,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_CONE)
tracerApp.setProperty(n_particles)
#Phosphor efficiencies:
p_eff1 = Property.Property(value=0.8, propID=PropertyID.PID_PhosphorEfficiency, valueType=ValueType.Scalar, time=0.0, units=None, objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(p_eff1)
p_eff2 = Property.Property(value=0.7, propID=PropertyID.PID_PhosphorEfficiency, valueType=ValueType.Scalar, time=0.0, units=None, objectID=objID.OBJ_PARTICLE_TYPE_2)
tracerApp.setProperty(p_eff2)
# Emission spectrum
em = Property.Property(value=ex_em_import.getEm(),
propID=PropertyID.PID_EmissionSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(em)
# Excitation spectrum
ex = Property.Property(value=ex_em_import.getEx(),
propID=PropertyID.PID_ExcitationSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(ex)
# 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)
# Emission spectrum 2
em2 = Property.Property(value=ex_em_import.getEmRed(),
propID=PropertyID.PID_EmissionSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_2)
tracerApp.setProperty(em2)
# Excitation spectrum 2
ex2 = Property.Property(value=ex_em_import.getExRed(),
propID=PropertyID.PID_ExcitationSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_2)
tracerApp.setProperty(ex2)
# Absorption spectrum 2
aabs2 = Property.Property(value=ex_em_import.getAbsRed(),
propID=PropertyID.PID_AsorptionSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_2)
tracerApp.setProperty(aabs2)
# Solve Mie
mieApp.solveStep(0)
# Connect functions
pScat = mieApp.getProperty(PropertyID.PID_ScatteringCrossSections, 0,
objectID=objID.OBJ_PARTICLE_TYPE_1)
pPhase = mieApp.getProperty(PropertyID.PID_InverseCumulativeDist, 0,
objectID=objID.OBJ_PARTICLE_TYPE_1)
pS2 = mieApp.getProperty(PropertyID.PID_ScatteringCrossSections, 0.0,
objectID=objID.OBJ_PARTICLE_TYPE_2)
pP2 = mieApp.getProperty(PropertyID.PID_InverseCumulativeDist, 0.0,
objectID=objID.OBJ_PARTICLE_TYPE_2)
tracerApp.setProperty(pScat)
tracerApp.setProperty(pPhase)
tracerApp.setProperty(pS2)
tracerApp.setProperty(pP2)
particles_in_um3 = []
j=0
for w_frac in weight_frac:
particles_in_um3.extend([st.particlesInVolumeLogNormWeightTotal(w_frac, dens_p, dens_host, mu[j], sigma[j], particle_diameters=d)])
j = j+1
# Particle density
#NOTE: If more than 1 particle type, use ValueType.Vector !
vDens = particles_in_um3
pDens = Property.Property(value=vDens,
propID=PropertyID.PID_ParticleNumberDensity,
valueType=ValueType.Vector,
time=0.0,
units=None,
objectID=objID.OBJ_CONE)
tracerApp.setProperty(pDens)
# Solve
tracerApp.solveStep(0, runInBackground=False)
#get the simulation result value and return it:
p = tracerApp.getProperty(PropertyID.PID_LEDCCT, 0,
objectID=objID.OBJ_LED)
print('CCT: %f' % p.getValue())
p = tracerApp.getProperty(PropertyID.PID_LEDRadiantPower, 0,
objectID=objID.OBJ_LED)
print('RadiantPower: %f' % p.getValue())
p = tracerApp.getProperty(PropertyID.PID_LEDColor_x, 0,
objectID=objID.OBJ_LED)
print('Color x: %f' % p.getValue())
p = tracerApp.getProperty(PropertyID.PID_LEDColor_y, 0,
objectID=objID.OBJ_LED)
print('Color y: %f' % p.getValue())
p = tracerApp.getProperty(PropertyID.PID_LEDSpectrum, 0,
objectID=objID.OBJ_LED)
print(p.getValue()['wavelengths'][0:10])
print(p.getValue()['intensities'][0:10])
#return CCT:
p = tracerApp.getProperty(PropertyID.PID_LEDCCT, 0,
objectID=objID.OBJ_LED)
return (p.getValue())
def xsim(z, *params):
#create the apps:
mieApp = MMPMie('localhost')
tracerApp = MMPRaytracer('localhost')
comsolApp = MMPComsolDummy('localhost')
# Point data conversion to false. Speeds up testing
tracerApp._convertPointData = False
# Set default LED json
tracerApp.setDefaultInputFile('./DefaultLED.json')
# Connect functions
pScat = mieApp.getProperty(PropertyID.PID_ScatteringCrossSections, 0,
objectID=objID.OBJ_PARTICLE_TYPE_1)
pPhase = mieApp.getProperty(PropertyID.PID_InverseCumulativeDist, 0,
objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(pScat)
tracerApp.setProperty(pPhase)
# Connect fields
fTemp = comsolApp.getField(FieldID.FID_Temperature, 0)
fHeat = comsolApp.getField(FieldID.FID_Thermal_absorption_volume, 0)
tracerApp.setField(fTemp)
tracerApp.setField(fHeat)
#set the constant properties from *params:
p_max = 35.0
p_min = 3.0
p_num = 50
w_max = 1100.0
w_min = 100.0
w_num = 1000
# Weight fractions
#weight_frac = np.array([24]) / 100.0 #CHANGING PARAMETER!See below
# Particle density
dens_p = 5.0 # g/cm3
# Silicone density
dens_host = 1.1 # g/cm3
# log mean in microns
mu = 2.4849
# log standard deviation in microns
sigma = 0.3878
# Particles diameters in micrometers
d = np.linspace(p_min, p_max, p_num)
pdf = lognorm(sigma, scale=np.exp(mu))
#set the changing properties from z:
# Weight fractions
weight_frac = np.array([z]) / 100.0
print("z, weight_frac = ", z, weight_frac)
#run the simulation:
# Number of rays to trace
pRays = Property.Property(value=100000,
propID=PropertyID.PID_NumberOfRays,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_CONE)
tracerApp.setProperty(pRays)
# Emission spectrum
em = Property.Property(value=ex_em_import.getEm(),
propID=PropertyID.PID_EmissionSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(em)
# Excitation spectrum
ex = Property.Property(value=ex_em_import.getEx(),
propID=PropertyID.PID_ExcitationSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1)
tracerApp.setProperty(ex)
# 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)
#logger.info('Properties set!')
# Solve Mie
mieApp.solveStep(0)
for w_frac in weight_frac:
particles_in_um3 =\
st.particlesInVolumeLogNormWeightTotal(w_frac, dens_p,
dens_host, mu, sigma,
particle_diameters=d)
# Particle density
vDens = particles_in_um3
pDens = Property.Property(value=vDens,
propID=PropertyID.PID_ParticleNumberDensity,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_CONE)
tracerApp.setProperty(pDens)
# Solve
tracerApp.solveStep(0, runInBackground=False)
#get the simulation result value and return it:
p = tracerApp.getProperty(PropertyID.PID_LEDCCT, 0,
objectID=objID.OBJ_LED)
print('CCT: %f' % p.value)
p = tracerApp.getProperty(PropertyID.PID_LEDRadiantPower, 0,
objectID=objID.OBJ_LED)
print('RadiantPower: %f' % p.value)
p = tracerApp.getProperty(PropertyID.PID_LEDColor_x, 0,
objectID=objID.OBJ_LED)
print('Color x: %f' % p.value)
p = tracerApp.getProperty(PropertyID.PID_LEDColor_y, 0,
objectID=objID.OBJ_LED)
print('Color y: %f' % p.value)
p = tracerApp.getProperty(PropertyID.PID_LEDSpectrum, 0,
objectID=objID.OBJ_LED)
print(p.value['wavelengths'][0:10])
print(p.value['intensities'][0:10])
#return CCT:
p = tracerApp.getProperty(PropertyID.PID_LEDCCT, 0,
objectID=objID.OBJ_LED)
return (p.value)
tracerApp.setProperty(em)
# Excitation spectrum
ex = Property.Property(
value=ex_em_import.getEx(),
propID=PropertyID.PID_ExcitationSpectrum,
valueType=ValueType.Scalar,
time=0.0,
units=None,
objectID=objID.OBJ_PARTICLE_TYPE_1,
)
tracerApp.setProperty(ex)
# 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)
# Emission spectrum 2
em2 = Property.Property(
value=ex_em_import.getEmRed(),
propID=PropertyID.PID_EmissionSpectrum,
valueType=ValueType.Scalar,
time=0.0,