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

wangweizong · March 10, 2012, 15:38

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

phoenics_cfd · March 21, 2012, 10:27

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.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
strong conservation form of eulers equations	Dahyalal Bhate	Main CFD Forum	0	September 5, 2002 19:22
2D CFD code using SIMPLE algorithm	bfan	Main CFD Forum	3	June 22, 2002 23:01
Crank-Nicolson technique for Conservation form of equations	Yogesh Talekar	Main CFD Forum	2	July 20, 1999 02:45
Why FVM for high-Re flows?	Zhong Lei	Main CFD Forum	23	May 14, 1999 14:22
CONSERVATION EQUATIONS IN NONINERTIAL SYSTEM!!!	zhang zhanwang	Main CFD Forum	0	August 17, 1998 08:38

March 21, 2012, 10:27		#2
phoenics_cfd Member michael malin Join Date: Mar 2009 Location: london Posts: 35 Rep Power: 18	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.