[ziemowitzima]

Dear Foamers,
I am modeling cold gas injection (5K) into a tunnel filled with warm other gas (300K).
I noticed that cold gas is not expanding properly, there is to little of it in the tunnel.

did any of you encountered similar problem ?

regards
mm

[mturcios777]

What are the boundary conditions of your injection? When you say too little gas, do you mean compared with experiment?

[ziemowitzima]

1. inlet/ejection BC is:
a) for U - time dependent mass flow of the cold gas1.
b) for p - zeroGradient
c) for gas1 (cold) - constant = 1
d) for bulk gas2 (warm) - constant = 0
e) for temperature T = 5K

2. yes it is too little comparing to experiment and comparing to fluent results

I am using reactinFOAM with chemistry off.

Thanks for your interest !

mm

[mturcios777]

Originally Posted by ziemowitzima

1. inlet/ejection BC is:
a) for U - time dependent mass flow of the cold gas1.
b) for p - zeroGradient
c) for gas1 (cold) - constant = 1
d) for bulk gas2 (warm) - constant = 0
e) for temperature T = 5K

2. yes it is too little comparing to experiment and comparing to fluent results

I am using reactinFOAM with chemistry off.

Thanks for your interest !

mm

No problem. I'm not sure about your conditions for c/d, is that the mass fraction? I've found slightly better success with injection when the inlet pressure condition is inletOutlet. Also, what gases are you injecting, and what is the thermo model you are using?

[ziemowitzima]

yes,
c) and d) are just mass fraction.

how would you define inletOutlet condition for pressure ?
I do not know the value of the pressure at the eject/inlet.

I am injecting Helium, modelled as a perfect gas, but this assumption should be fine, because the same is used in Fluent.

for thermo model I use:

thermoType
{
type hePsiThermo;
mixture reactingMixture;
transport sutherland;
thermo janaf;
energy sensibleEnthalpy;
equationOfState perfectGas;
specie specie;
}

[mturcios777]

Well, inletOutlet switches the boundary between being constant value or zeroGradient depending on if the flow is into or out of the domain. Its more a stability helper for me. The extra entry is inletValue, which is the constant value the pressure will be if the flow is into the domain.

Try that with some pressure that is close to what you think it should be. Is this gas coming from a resevoir? Also, is the mesh grading similar to what the Fluent simulation used?

[ziemowitzima]

Hi,
The gas in forcedly ejected in to the tunnel with the mas flow approx 1kg/s.
So the flow is always into the domain (because it is inlet).

As far as I know, inletOutle BC should be used at the outlet from the domain ?