Density problem in fluent
 

Hello

I am modeling a 3d spillway in Fluent.(vof multiphase)
Depths and speeds and pressures are correct. (in free surface volume fraction is 0.5) But in Density or Air concentration (C) i have a problem that it is confused me.
In free surface(volume fraction=0.5) of my model density is about 500 kg/m3 (air concentration is C=50) !!!
I mean in upper than free surface it is going to 1.225 (air density) (C=90 or 95 or 99).
I mean in upper than free surface the density should be about 1.225 (air concentration C=100) .
I mean after free surface it should be just air with density=1.225

Please help
Thanks

I am sorry, but I don't understand what you are saying... What is your problem? Which thing is not what you thought it should be?

???? Your pic and profile are correct! What do you want? You are plotting the density of mixture and that is why you are getting averaged values

By the Way do you mean with C the void fraction

yes c is void fraction.
in free surface density should be 1.225 or upper than but in fluent it is 500.
please see attached image

You are using VOF have u read something about phase average of the navier stokes equations. You will never get your 998 at the free surface from any programm because you are solving the phasic average NS-Equations: Your free surface is of a finite thickness this thickness does depend on your mesh your discretisation order and the volume fraction recontrusction say tracking or caprturing. Your result are OK. At the free surface you have 50% probability to have air as same as water that is why you will get there the arithmetic average almost

Thanks for your replies.

I want to understand this problem for my thesis.
If in defence meeting, the judges ask me, why in free surface the air concentration is 0.5, (because in common and nature it should be 0.9 or upper) how can i explaine this problem? What justifications i should offer? If i say my modeling is OK and it is just related to Fluent VOF, they will not accept. :( (because in common and nature air concentration should be 0.9 or upper):confused:

You define free surface as (volume fraction=0.5), so that is what you get. If you would define free surface as (volume fraction=0.9), you would get a different result. This answer might be a bit too simple, because maybe I don't understand your question.

Is your question: "why is the density on the free surface, which I define as having a liquid volume fraction of 50%, around 500 kg/m3, while I would expect it to have a liquid volume fraction of 90% and therefore have a density of 900 kg/m3?"?

I am defining free surface in volume fraction=0.5 because it is usual in many papers. With defining free surface in volume fraction=0.9 the depth of flow is not same with experimental data.(it is larger than main depth)
Yes. I want with defining free surface in volume fraction=0.5, the density in free surface be about 1.225 (close to air). If it is not possible and my simulation is OK, so how can i doing with this problem? I need an explanation for this problem.

In theory guide (Fluent version 14.5), section 1.3.2, equation (1-15):



Here is the mixture density, is the volume fraction of phase l and is the density of phase l.
In your case, both alpha's (volume fractions) are 0.5, so:



And just put in the numbers, this value is 500 kg/m3 in your case. This can easily be calculated without simulation, so your simulation is giving the correct numbers.

The concept is not so difficult to understand: in your 'free surface', half of the cell is filled with water, the other half with vapor. The density will therefore be exactly in the middle. You can only get a different free surface density if you define your free surface in a different way.

You have first of all to be familiar with the theory behind what you are calculating say multiphase flows and how to get the underlying physics from the conservations equations of mass, momentum and energy. We helped there more than necessary it is now up to you to learn or not.

Keep it real

I am familiar with the theory behind what i am calculating multiphase flows.

Can i equivalent density=500 with density=1.225 ? I mean, moving upper density from free surface to down.

If you desparately want the free surface to have a density of 1.225, you should change the density of water to 1.225. It will not make any sense, but if it makes you happy, just do it.

I think you are asking the wrong question. You should be asking for yourself:
- What is a free surface exactly?
- Why do I think that the free surface should have a density of 1.225, and is that correct? (Hint: "because I read it in an article" is not sufficient, you should understand it)

Your questions show that you are not familiar with the theory behind calculating multiphase flows. That is not something to be ashamed of, I am not familiar with many theories.

My last pot here:

if you consider the surface of water: and you introduce a probe to capture the surface and the water is moving the probability to get the surface or water is 50% that is you in the phasic average NS you get an average density. In the realtiy you only have either water or air this leading to discontinuity on the surface but if you want to calculate that you have to smooth the gradient this yields to an average quantity..

Does the mesh effects in distance from contours?
I mean if my mesh be smaller than present, then the contour of density=1.225 is closer to density=500?

generally the free surface has to be fine meshed let say 10-15 cells mesh resolution would affect the solution especially if you want to resolve ther free-surface. Talking now about the density i still do not know why you think you have a problem. Your problem is not the simulation rather your understanding of what you expect from the simulation

Just tell me.
Can density of mesh effects on distance between contours?
I mean if i use denser mesh, the distance between contour of 0.5 and contour of 0.9 is fix or can change in relative to big mesh?

if you soliton is grid independent then you will get nor change otherweise some changes will occur. But you will still get the average density of 500 at the free surface..

If volume fraction=0.5 is free surface then can you see (without considering VOF and theories), what are the upper volume fractions? Are they wave or entrapped air or entrained air?

Try to select/deselect "node values" when plot contours or on 2d chart. I think you will get what you want.

Daniele

I tried but it is not effecting more in my results . What is "node values"?