Electric potential scaling in case of two temperature electric arc simulation
Posted: Sun Mar 10, 2019 11:40 pm
Hello evereone,
I am developing a two temperature model of electric arc in Code_Saturne. In this model I have two thermal scalars: electron enthalpy and heavy species enthalpy. These two scalars have diffferent diffusivity and different source terms. The source term for electrons is Joule power and radiation loss while for heavy species it is the heat transfer from electrons. Hence I have to reach the balance between them, The transport properties (thermal and electrical conductivities, viscosity) are defined by the electrons temperature Te and the ratio of temperatures Te/Th.
When I switched my electric arc model from the LTE simulation with just one thermal scalar to two thermal scalars I found that the arc fades away every time. The total Joule power goes down as well as the voltage. I started looking for a reason and found that in the arc voltage scaling algorithm we have the following formulas for the correction coefficient:
Kcorrection = coepoa = Iimposed/Icomputed = Iimposed*Vimposed/Qtotal_computed
but if coepoa >1.05 then
Kcorrection = coepot = SQRT(20*MIN{(Enthalpy per kg in cell i)/(Joule power per kg in cell i)})
So I corrected it in order to take into account both enthalpies, but the arc still fades away.
Could you please tell me what am I doing wrong with this scaling algorithm and what is the origin of the second formula for the correction coefficient? Probably it must be changed to implement the 2T model? I am a little concerned especially by the number 20 in the formula.
Thank you all in advance,
Rodion
I am developing a two temperature model of electric arc in Code_Saturne. In this model I have two thermal scalars: electron enthalpy and heavy species enthalpy. These two scalars have diffferent diffusivity and different source terms. The source term for electrons is Joule power and radiation loss while for heavy species it is the heat transfer from electrons. Hence I have to reach the balance between them, The transport properties (thermal and electrical conductivities, viscosity) are defined by the electrons temperature Te and the ratio of temperatures Te/Th.
When I switched my electric arc model from the LTE simulation with just one thermal scalar to two thermal scalars I found that the arc fades away every time. The total Joule power goes down as well as the voltage. I started looking for a reason and found that in the arc voltage scaling algorithm we have the following formulas for the correction coefficient:
Kcorrection = coepoa = Iimposed/Icomputed = Iimposed*Vimposed/Qtotal_computed
but if coepoa >1.05 then
Kcorrection = coepot = SQRT(20*MIN{(Enthalpy per kg in cell i)/(Joule power per kg in cell i)})
So I corrected it in order to take into account both enthalpies, but the arc still fades away.
Could you please tell me what am I doing wrong with this scaling algorithm and what is the origin of the second formula for the correction coefficient? Probably it must be changed to implement the 2T model? I am a little concerned especially by the number 20 in the formula.
Thank you all in advance,
Rodion