CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > FLUENT

VOF mass loss

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   November 27, 2012, 05:31
Post VOF mass loss
  #1
Member
 
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 4
Smaras is on a distinguished road
Again one more time as i havent got any repsonse

Hello all,

I am doing VOF simulation with 2ddp. Now using two phase air and water. Taking water as first phase and air as second. With small inlet with 96 m/s air impingement over water surface. Lance to water bath distance is 154 mm with 100 mm depth of water bath and having 100 mm radius of cylinder. Outlet is at top defined as pressure outlet. Now the problem is that when i am doing calculation it's diverging but the results are totally disappointing as the penetration depth of the air jet should have reached 14 mm while it simulation its stuck at 7.5 and further there is a mass loss of around 0.1 Kg out of 3 Kg.
Have used coarser mesh and its showing better result. Now according to the theory the finer mesh should have given better result but here it's 180 degrees.

And further i am using URF as
pressure 0.2
momentum 0.3
Tur. dissipation rate 0.8
Tur. viscosity 0.8
Step size 0.0001

Non iterative with Factional Scheme
Volume Fraction = geo reconstruct
Turb. KE = quick

Would be thankful if some would give suggestions.

Regards,
Smaras



Smaras is offline   Reply With Quote

Old   November 27, 2012, 07:37
Default
  #2
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Quote:
Originally Posted by Smaras View Post
Now the problem is that when i am doing calculation it's diverging but the results are totally disappointing...
You should not look at results before you fix the numerics.
Get a converging simulation!
RodriguezFatz is online now   Reply With Quote

Old   November 27, 2012, 07:50
Post
  #3
Member
 
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 4
Smaras is on a distinguished road
Quote:
Originally Posted by RodriguezFatz View Post
You should not look at results before you fix the numerics.
Get a converging simulation!
Thanks Rodriguez but i have a converging solution that's the problem. And the numerical calculation and also some research paper i am consulting also giving the same depth of penetration which is not confirmed with fluent.

Therefore i am doing some mistake. I need help in this regard.

Regards,
Smaras

Last edited by Smaras; November 27, 2012 at 07:51. Reason: mistakes
Smaras is offline   Reply With Quote

Old   November 27, 2012, 07:55
Default
  #4
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Do you say that in your originial post you meant "converging" instead of "diverging"?
RodriguezFatz is online now   Reply With Quote

Old   November 27, 2012, 07:57
Default
  #5
Member
 
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 4
Smaras is on a distinguished road
Quote:
Originally Posted by RodriguezFatz View Post
Do you say that in your originial post you meant "converging" instead of "diverging"?
Mistake,

I meant to say diverging from the numerical result and the results from the research paper

Thanks for correction.

Regards
Smaras is offline   Reply With Quote

Old   November 27, 2012, 08:01
Default
  #6
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Alright, now I get you.

Can you post some pics of your setup and mesh? And also the residuals...
What I don't understand: How can your simulation converge, when you have a mass loss?
RodriguezFatz is online now   Reply With Quote

Old   November 27, 2012, 09:14
Post
  #7
Member
 
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 4
Smaras is on a distinguished road
Quote:
Originally Posted by RodriguezFatz View Post
Alright, now I get you.

Can you post some pics of your setup and mesh? And also the residuals...
What I don't understand: How can your simulation converge, when you have a mass loss?
That i cannot understand as well.

Mass-Loss:






Mesh:





Penetration Depth:



Residuals:



Phases:




If anything else is left please i can post that as well. Vielen Dank Rodriguez

Regards,

Smaras
Smaras is offline   Reply With Quote

Old   November 27, 2012, 09:39
Default
  #8
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Just to be sure, I am getting your setup right:

1) The box you are simulating is a cylinder, the axis is on the left of your pictures?
2) On the bottom you got the water and a needle blasts air vertically on your water surface?
3) The needle is on the top left of your domain, where I can see the little gap in the green grid?
4) "Penetration depth" is the (average) difference of the height of the water surface at your symmetry axis, compared to the steady state without the needle fan?
5) Air from the needle escapes through the pressure outlet at the top?
6) On the right and on the bottom are walls?
RodriguezFatz is online now   Reply With Quote

Old   November 27, 2012, 09:59
Default
  #9
Member
 
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 4
Smaras is on a distinguished road
Quote:
Originally Posted by RodriguezFatz View Post
Just to be sure, I am getting your setup right:

1) The box you are simulating is a cylinder, the axis is on the left of your pictures?
2) On the bottom you got the water and a needle blasts air vertically on your water surface?
3) The needle is on the top left of your domain, where I can see the little gap in the green grid?
4) "Penetration depth" is the (average) difference of the height of the water surface at your symmetry axis, compared to the steady state without the needle fan?
5) Air from the needle escapes through the pressure outlet at the top?
6) On the right and on the bottom are walls?


I hope this image will clear the queries. While for the 4th
Penetration depth is the depth difference between the stable condition and the impingement of air onto water surface. the distance calculated is from the top i.e. Outlet.

And all dimension are in mm.

Regards.
Smaras is offline   Reply With Quote

Old   November 27, 2012, 10:04
Default
  #10
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Correct me if I am wrong, but as I understand it, the needle (lower left in your last picture) shouldn't be an axis, but a wall boundary. There can be just one cylinder axis.
RodriguezFatz is online now   Reply With Quote

Old   November 27, 2012, 10:09
Default
  #11
Member
 
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 4
Smaras is on a distinguished road
Quote:
Originally Posted by RodriguezFatz View Post
Correct me if I am wrong, but as I understand it, the needle (lower left in your last picture) shouldn't be an axis, but a wall boundary. There can be just one cylinder axis.
Ohh it was by mistake typing error. The top left is nozzle wall. While the lower left having length 254 is axis.
Smaras is offline   Reply With Quote

Old   November 27, 2012, 11:33
Default
  #12
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Ok, is your simulation time dependent? The residuals screenshot looks a bit strange... it looks like your first 2000 iterations have a residual in continuity equation of about 1.0
RodriguezFatz is online now   Reply With Quote

Old   November 27, 2012, 11:47
Default
  #13
Member
 
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 4
Smaras is on a distinguished road
Quote:
Originally Posted by RodriguezFatz View Post
Ok, is your simulation time dependent? The residuals screenshot looks a bit strange... it looks like your first 2000 iterations have a residual in continuity equation of about 1.0
Yup it is time dependent since i am using Fractional Scheme with NITA.
Smaras is offline   Reply With Quote

Old   November 27, 2012, 11:59
Default
  #14
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
But look at your residuals. Something ugly happens during the first 2000 iterations.

You can also do this: Under "Monitors" add a surface monitor and integrate the water density over your domain. Plot it each timestep to a window. I guess when your residuals start to converge (after about 2000 iterations), your integrals won't change any more over time, thus the total mass will be conserved.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is online now   Reply With Quote

Old   November 28, 2012, 04:33
Default
  #15
Member
 
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 4
Smaras is on a distinguished road
Quote:
Originally Posted by RodriguezFatz View Post
But look at your residuals. Something ugly happens during the first 2000 iterations.

You can also do this: Under "Monitors" add a surface monitor and integrate the water density over your domain. Plot it each timestep to a window. I guess when your residuals start to converge (after about 2000 iterations), your integrals won't change any more over time, thus the total mass will be conserved.
Rodriguez i decrease the time step size to 1e-3 from 1e-5.
Smaras is offline   Reply With Quote

Old   November 28, 2012, 04:35
Default
  #16
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Quote:
Originally Posted by Smaras View Post
Rodriguez i decrease the time step size to 1e-3 from 1e-5.
You could also try to keep the current stepsize and get convergence here! Maybe increase the number of iterations of each time step or change the solution algorithm.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is online now   Reply With Quote

Old   November 28, 2012, 04:41
Default
  #17
Member
 
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 4
Smaras is on a distinguished road
Quote:
Originally Posted by RodriguezFatz View Post
You could also try to keep the current stepsize and get convergence here! Maybe increase the number of iterations of each time step or change the solution algorithm.
By currewnt step size you mean 1e-3 or 1e-5
Smaras is offline   Reply With Quote

Old   November 28, 2012, 04:43
Default
  #18
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
I mean the old one (1e-3). I guess in your post above you meant "i decrease the time step size from 1e-3 to 1e-5" and not vice versa...
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is online now   Reply With Quote

Old   November 28, 2012, 04:45
Default
  #19
Member
 
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 4
Smaras is on a distinguished road
Quote:
Originally Posted by RodriguezFatz View Post
I mean the old one (1e-3). I guess in your post above you meant "i decrease the time step size from 1e-3 to 1e-5" and not vice versa...
other way i.e. from 1e-5 to 1e-3 that why its such a drastic change
Smaras is offline   Reply With Quote

Old   November 28, 2012, 05:04
Default
  #20
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
I am confused... you can not decrease to 1e-3 from 1e-5. This is an increase.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is online now   Reply With Quote

Reply

Tags
penetration, phase change, vof modeling

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
VOF mass loss Smaras FLUENT 2 November 20, 2012 07:36
Source terms due to mass transfer in VOF model ssamton FLUENT 0 March 5, 2012 01:03
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
Mass loss! So what? jinwon park Main CFD Forum 13 May 22, 2008 09:29
HELP! adding a mass source to VOF eqn. by UDF??? ROOZBEH FLUENT 4 May 29, 2003 09:54


All times are GMT -4. The time now is 08:59.