Modeling lamps using commercial packages

1
Modeling lamps using commercial packages
COST Model Inventory Workshop, Funchal, April
2005
ROUFFET Jean-Baptiste CPAT Université de
Toulouse, France
2
CONTENTS

• Our modeling of HID lamps
• Use of FEMLAB
• Why have we choose commercials packages?
• Conclusion

COST Model Inventory Workshop, Funchal, April
2005
3
OUR MODELING OF HID LAMPS
LTE
FLUID MODEL
We have to calculate transport coefficients
versus T and p

• ASUMPTIONS
• LTE
• No demixion
• No erosion
• Racte
• DC

We have to calculate discharge composition versus
T and p
COST Model Inventory Workshop, Funchal, April
2005
4
OUR MODELING OF HID LAMPS
We calculate the pressure that we have in the
burner versus wall temperature for a given mixture
MTDaTa
We calculate discharge composition versus T and
p(Tw)
We calculate and tabulate transport coefficients
versus T and p

• We solve equations
• Energy balance
• Momentum balance
• Current conservation

FEMLAB
COST Model Inventory Workshop, Funchal, April
2005
5
OUR MODELING OF HID LAMPS
DISCHARGE PCA
Anode, we choose V0
Boundaries conditions given by the M. Benilovs
team model
g
COST Model Inventory Workshop, Funchal, April
2005
6
OUR MODELING OF HID LAMPS
Current density
Temperature
COST Model Inventory Workshop, Funchal, April
2005
7
FEMLAB
COST Model Inventory Workshop, Funchal, April
2005
8
PLAN
COST Model Inventory Workshop, Funchal, April
2005
9
FEMLAB
COST Model Inventory Workshop, Funchal, April
2005
10
FEMLAB
COST Model Inventory Workshop, Funchal, April
2005
11
FEMLAB
COST Model Inventory Workshop, Funchal, April
2005
12
FEMLAB
COST Model Inventory Workshop, Funchal, April
2005
13
FEMLAB
COST Model Inventory Workshop, Funchal, April
2005
14
FEMLAB
COST Model Inventory Workshop, Funchal, April
2005
15
FEMLAB
TEMPERATURE
VELOCITY
CURENT (Log)
COST Model Inventory Workshop, Funchal, April
2005
16
Why have we choose commercials packages?
Lamps can be considered like a converter of
electric energy in radiative energy
COST Model Inventory Workshop, Funchal, April
2005
17
Why have we choose commercials packages?
T
V (or the lamps conductivity)
j
POWER SUPPLY
Electric network
LAMP
CANDELABRA
ELECTRODES
T
I
T
Cathodique fall
OUT PUT OF THE SYSTEM
COST Model Inventory Workshop, Funchal, April
2005
18
Why have we choose commercials packages?

• LAMP MODULE
• BALLAST MODULE
• ELECTRODE MODULE

Simulink is a Matlab toolbox that enable to make
model with module
FEMLAB CAN GENERATE SIMULINK BLOCK AND IS
MULTIPHYSICS
T
V (or the lamps conductivity)
j
POWER SUPPLY
Electric network
LAMP
CANDELABRA
ELECTRODES
T
I
T
Cathodique fall
OUT PUT OF THE SYSTEM
COST Model Inventory Workshop, Funchal, April
2005
19
Why have we choose commercials packages?
Algorithms development are made ? optimization is
made (rapidity) Conviviality (user
interfaces) User community Support Cost no
access to all parameters

-
COST Model Inventory Workshop, Funchal, April
2005
20
CONCLUSION
The use of FEMLAB enable us to have quickly
results without to develop algorithms and
programs The possibility of convert a FEMLAB
model in Simulink is a good things when we want
to work with other peoples (collaboration) We
have to add demixion and go toward the 3D
LIMITATIONS ? We cant model very high
compressible flow (shock tube) ? 3D (memory
problems)
COST Model Inventory Workshop, Funchal, April
2005
21
THANK YOU
COST Model Inventory Workshop, Funchal, April
2005
22
COST Model Inventory Workshop, Funchal, April
2005
23
COST Model Inventory Workshop, Funchal, April
2005
24
COST Model Inventory Workshop, Funchal, April
2005
25
COST Model Inventory Workshop, Funchal, April
2005
26
COST Model Inventory Workshop, Funchal, April
2005
27
Les équations utilisées pour modéliser le plasma
28
Les conditions aux limites
29
Le projet NumeLiTe
For Twgt1090K, we have in gas phase
62,5 of Hg 18.75 of Tl
18,75 of I
Mercury vaporisation
Thallium iodine vaporisation
Only gas phase
Evolution of the relative proportion of elements
in gas phase versus wall temperature.
Evolution of system pressure versus wall
temperature.
Wall temperature changes system pressure and
modifies proportions of each element in gas phase
30
Calcul des densités des espèces
31
Les coefficients de transport
Tw Psys Hg 700 K 1,83 atm
99,90 800 K 2,10 atm 98,88 900 K
2,43 atm 93,89 1000 K 2,93 atm
80,05 1100 K 3,73 atm 62,50 1400 K 4,76
atm 62,50
32
Le coefficient démission net
Pour un plasma à léquilibre thermodynamique
local, on peut relier le coefficient démission
nette à la divergence du flux radiatif
En considérant un plasma homogène et isotherme de
rayon Rp, nous pouvons écrire
étant le coefficient dabsorption monochromatique
corrigé de lémission stimulée
33
Le coefficient démission nette
34
Les équations du modèle
35
Résultats
Temperature profile versus radius for different z
r (m)
z
r