Implementing phase change in a heat transfer solver

ThibautDer · May 6, 2016, 09:32

Hi,

I'm working on a heat transfer solver for incompressible flows. In particular, this solver must be able to take account of phase change. To do this I would like to compute the dryness by implementing this:

If $h < h_{sat}$

$T = \frac{h}{c_p}$

else

$x = \left( h - \int^{T_{sat}}_{T_0} c_p dT \right) \bigg/ h_{fg}$

and

$T = T_{sat}$

I have no idea how to implement the if condition in a way to go through each mesh cells. Does anyone know how to do so in OpenFOAM ?

Thanks in advance

klilla · May 9, 2016, 08:38

You can go through every volScalar/Vector/Tensor field as:

forAll(h, cellI)
{
if(h[cellI]<hSat)
{
T[cellI] = ...;
}
....

}

This will loop over all the cells where h is defined, practically all the cells.
Best,
Lilla

ThibautDer · May 9, 2016, 08:49

Thank you !

Zhiheng Wang · May 9, 2016, 15:18

Hi kapa Lilla

I could not get you

"
You can go through every volScalar/Vector/Tensor field as:

forAll(h, cellI)
{
if(h[cellI]<hSat)
{
T[cellI] = ...;
}
....

}
"
Can you please elaborate a little ????

May 6, 2016, 09:32	Implementing phase change in a heat transfer solver	#1
ThibautDer New Member Join Date: Apr 2016 Posts: 2 Rep Power: 0	Hi, I'm working on a heat transfer solver for incompressible flows. In particular, this solver must be able to take account of phase change. To do this I would like to compute the dryness by implementing this: If $h < h_{sat}$ $T = \frac{h}{c_p}$ else $x = \left( h - \int^{T_{sat}}_{T_0} c_p dT \right) \bigg/ h_{fg}$ and $T = T_{sat}$ I have no idea how to implement the if condition in a way to go through each mesh cells. Does anyone know how to do so in OpenFOAM ? Thanks in advance Kummi likes this. Last edited by ThibautDer; May 9, 2016 at 06:21.

May 9, 2016, 08:38		#2
klilla Member Kapa Lilla Join Date: Mar 2009 Location: Bruxelles, Belgium Posts: 57 Rep Power: 17	You can go through every volScalar/Vector/Tensor field as: forAll(h, cellI) { if(h[cellI]<hSat) { T[cellI] = ...; } .... } This will loop over all the cells where h is defined, practically all the cells. Best, Lilla Kummi likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
How to simulate the flow and heat transfer of phase change materials?	rosol	FLUENT	9	April 10, 2016 14:22
Looking for a solver (Mach 0.4, Turbulent, Heat Transfer, Second Order)	fredo490	OpenFOAM Running, Solving & CFD	3	March 27, 2013 04:07
Water subcooled boiling	Attesz	CFX	7	January 5, 2013 03:32
what approach for direct mixing fluids with heat transfer and phase change?	hennas	Main CFD Forum	0	June 9, 2009 15:44
How can I increase Heat Transfer at Domain Interf?	B.Simon	CFX	3	October 28, 2008 18:53

May 9, 2016, 08:49		#3
ThibautDer New Member Join Date: Apr 2016 Posts: 2 Rep Power: 0	Thank you !

May 9, 2016, 15:18		#4
Zhiheng Wang Member Zhiheng Wang Join Date: Mar 2016 Posts: 72 Rep Power: 10	Hi kapa Lilla I could not get you " You can go through every volScalar/Vector/Tensor field as: forAll(h, cellI) { if(h[cellI]<hSat) { T[cellI] = ...; } .... } " Can you please elaborate a little ????