How to call the subroutine to compute two temperature energy conservation equations?
 

I am working on the simulation of two temperature arc model in which the energy conservation governing equation is in the form of electron temperature and heavy particle temperature ( see attached file). I can use “solve(TEM1)” to solve one electron energy conservation equation but the problem is how to call the internal subroutine to solve the other energy conservation equation in the form of heavy particle temperature.

I wonder whether some guy can give some suggestion on this problem.

Many thanks

WeizongWang

Since the continuity and momentum equations appear to be shared by the two phases, you can simply generate a second energy equation by activating the solution of a scalar variable, say C1 which you can name T2. You will then need to modify each term in this scalar transport equation ( transient, convection, diffusion ) so as to convert it into an energy equation for the solid phase. This can be done by means of the INFORM facility, which provide Q1 commands to modify such terms, i.e MODCON, MODDIF, etc. It will then be necessary to introduce any source terms and interphase source terms, with appropriate linearisation. Again, this can be done in the Q1 file using the INFORM facility, which is documented here:
http://www.cham.co.uk/phoenics/d_pol...nc/in-form.htm.