350 likes | 503 Views
Modeling lamps using commercial packages. COST – Model Inventory Workshop, Funchal, April 2005. ROUFFET Jean-Baptiste CPAT – Université de Toulouse, France. CONTENTS. Our modeling of HID lamps Use of FEMLAB Why have we choose commercials packages? Conclusion.
E N D
Modeling lamps using commercial packages COST – Model Inventory Workshop, Funchal, April 2005 ROUFFET Jean-Baptiste CPAT – Université de Toulouse, France
CONTENTS • Our modeling of HID lamps • Use of FEMLAB • Why have we choose commercials packages? • Conclusion COST – Model Inventory Workshop, Funchal, April 2005
OUR MODELING OF HID LAMPS LTE FLUID MODEL We have to calculate transport coefficients versus T and p • ASUMPTIONS: • LTE • No demixion • No erosion • Ra=cte • DC We have to calculate discharge composition versus T and p COST – Model Inventory Workshop, Funchal, April 2005
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
OUR MODELING OF HID LAMPS DISCHARGE PCA Anode, we choose V=0 Boundaries conditions given by the M. Benilov’s team model g COST – Model Inventory Workshop, Funchal, April 2005
OUR MODELING OF HID LAMPS Current density Temperature COST – Model Inventory Workshop, Funchal, April 2005
FEMLAB COST – Model Inventory Workshop, Funchal, April 2005
PLAN COST – Model Inventory Workshop, Funchal, April 2005
FEMLAB COST – Model Inventory Workshop, Funchal, April 2005
FEMLAB COST – Model Inventory Workshop, Funchal, April 2005
FEMLAB COST – Model Inventory Workshop, Funchal, April 2005
FEMLAB COST – Model Inventory Workshop, Funchal, April 2005
FEMLAB COST – Model Inventory Workshop, Funchal, April 2005
FEMLAB COST – Model Inventory Workshop, Funchal, April 2005
FEMLAB VELOCITY TEMPERATURE CURENT (Log) COST – Model Inventory Workshop, Funchal, April 2005
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
Why have we choose commercials packages? T V (or the lamps conductivity) j POWER SUPPLY Electric network LAMP ELECTRODES CANDELABRA T I T Cathodique fall OUT PUT OF THE SYSTEM COST – Model Inventory Workshop, Funchal, April 2005
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 ELECTRODES CANDELABRA T I T Cathodique fall OUT PUT OF THE SYSTEM COST – Model Inventory Workshop, Funchal, April 2005
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
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 can’t model very high compressible flow (shock tube) 3D (memory problems) COST – Model Inventory Workshop, Funchal, April 2005
THANK YOU COST – Model Inventory Workshop, Funchal, April 2005
Le projet NumeLiTe For Tw>1090K, 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
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
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 d’absorption monochromatique corrigé de l’émission stimulée
Résultats Temperature profile versus radius for different z r (m) z r