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

a problem of incompressible two phases(water and gas) with VOF Method

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 17, 2011, 09:52
Default a problem of incompressible two phases(water and gas) with VOF Method
  #1
New Member
 
vivien
Join Date: Jan 2011
Posts: 8
Rep Power: 5
zhangweisnoopy is on a distinguished road
Hi,all
I'm building a 2D simulation program for a two-phase incompressible viscous fluid simulation (water and air) with VOF method. The VOF method I use is the Hirt & Nichols SOLA-VOF scheme but the reconstruction method is PLIC (Young).
In a two-phase approach, when the air water interface is no longer considered as a free-surface, but the interface is still reconstructed using the VOF algorithm. The function Fs is used to determine the
aggregated density value inside a grid cell, using weighted
averaging .The incompressible Navier-Stokes equations without surface tension is solved with variable density and Viscosity considering water and gas.

The effective density and Viscosity are computed as
rho = rho_liquid*Fs+(1-Fs)*rho_gas

Fs is the fraction of the cell that is filled with fluid in VOF.


I use the MAC staggered meshes.So density and Viscosity at both cell center and face should be caculated. I ues the program to simulate the break dam over a dry bed ,but it doesn't work because the pressure possion equation isn't convergent after several hundreds steps...It seems that the interface is not accurate(I don't know it is the problem) ..I can figure it out .
Give some advise please,thanks very much!!
zhangweisnoopy is offline   Reply With Quote

Old   August 8, 2011, 12:56
Default
  #2
Senior Member
 
Join Date: Aug 2011
Posts: 251
Rep Power: 6
leflix is on a distinguished road
First check if your Navier-Stokes solver performs well in monophasic case (same rho, same mu).
Then for large density and/or viscosity ratios instead of using arithmetic mean as you did, use rather harmonic mean:
1/rho = Fs(1/rho_liquid) + (1-Fs)(1/rho_gas)
endly check the scheme you use to solve the hyperbolic equation for Fs. It matters a lot ! bad schemes can induce pb of mass conservation and thus pb to solve your pressure equation which is related to that.
leflix is offline   Reply With Quote

Old   August 8, 2011, 21:28
Smile
  #3
New Member
 
vivien
Join Date: Jan 2011
Posts: 8
Rep Power: 5
zhangweisnoopy is on a distinguished road
Quote:
Originally Posted by leflix View Post
First check if your Navier-Stokes solver performs well in monophasic case (same rho, same mu).
Then for large density and/or viscosity ratios instead of using arithmetic mean as you did, use rather harmonic mean:
1/rho = Fs(1/rho_liquid) + (1-Fs)(1/rho_gas)
endly check the scheme you use to solve the hyperbolic equation for Fs. It matters a lot ! bad can induce pb of mass conservation and thus pb to solve your pressure equation which is related to that.
Thank you very much leflix!!
I will try what you suggest.Thanks for your help very much!
zhangweisnoopy is offline   Reply With Quote

Old   August 10, 2011, 22:47
Default
  #4
New Member
 
vivien
Join Date: Jan 2011
Posts: 8
Rep Power: 5
zhangweisnoopy is on a distinguished road
Quote:
Originally Posted by leflix View Post
First check if your Navier-Stokes solver performs well in monophasic case (same rho, same mu).
Then for large density and/or viscosity ratios instead of using arithmetic mean as you did, use rather harmonic mean:
1/rho = Fs(1/rho_liquid) + (1-Fs)(1/rho_gas)
endly check the scheme you use to solve the hyperbolic equation for Fs. It matters a lot ! bad schemes can induce pb of mass conservation and thus pb to solve your pressure equation which is related to that.


Hi , leflix!
I have tried what you suggest.In monophasic case,the Navier-Stokes solver performs well .
I have used harmonic mean for large density and/or viscosity ratios, it doesn't work because the pressure possion equation isn't convergent after several hundreds steps.
It seems that there are some problems with VOF,of which the reconstruction method is PLIC (Young),while the tracking method is Lagrangian method(Gueyffier) instead of Eulerian method(Hirt and Nichols).
But I have used the scheme to handle one-phase problem such as generating wave and dam break, considering air is void and getting velocity and pressure of free surface via interpolation. So can you give some adivise?
Thank you very much.

zhangweisnoopy is offline   Reply With Quote

Old   August 11, 2011, 06:00
Default
  #5
Senior Member
 
Join Date: Aug 2011
Posts: 251
Rep Power: 6
leflix is on a distinguished road
Hi Vivien,

Visibly the problem comes from a bad mass conservation. So if your code performs well in the monophasic case, you have to search around two points:
1) is the way the volume fraction Fs is computed
2) is the implementation of your reconstruction method PLIC
I think the problem should be there.

Which velocity-pressure coupling algorithm do you use ?
Have you tried to under-relax the pressure ?
leflix is offline   Reply With Quote

Old   August 12, 2011, 03:58
Default
  #6
New Member
 
vivien
Join Date: Jan 2011
Posts: 8
Rep Power: 5
zhangweisnoopy is on a distinguished road
Quote:
Originally Posted by leflix View Post
Hi Vivien,

Visibly the problem comes from a bad mass conservation. So if your code performs well in the monophasic case, you have to search around two points:
1) is the way the volume fraction Fs is computed
2) is the implementation of your reconstruction method PLIC
I think the problem should be there.

Which velocity-pressure coupling algorithm do you use ?
Have you tried to under-relax the pressure ?
Thank you very much !leflix
I want to show you some results for more advise~


I use the program to simulate the break dam over a dry bed. In this test(Numerical simulation of breaking waves using a two-phase flow model), a square computation domain with the length and height equal to 22.8 m is set. A water column with the width L and height 2L (L = 5.7 m) is assumed at the left side of the computation domain. And the spatial steps in the horizontal and vertical directions are 0.285while Δt is 0.0025.


When T=200*0.0025=0.5s, the figure of the flow field is

200.jpg

When T=269*0.0025s, the pressure Possion equation isn't convergent .The figure of the flow field is


269.jpg


I have used the VOF method to handle one-phase problem such as generating wave and dam break, considering air is void and getting velocity and pressure of free surface via interpolation. And there are no problem with mass conservation....


So I donít know what the problem is. Can you give me more advice?
zhangweisnoopy is offline   Reply With Quote

Old   August 12, 2011, 07:49
Default
  #7
Senior Member
 
Join Date: Aug 2011
Posts: 251
Rep Power: 6
leflix is on a distinguished road
Quote:
Originally Posted by zhangweisnoopy View Post

When T=269*0.0025s, the pressure Possion equation isn't convergent .The figure of the flow field is


Attachment 8789



It is clear on this figure that the interface is not well captured.


Quote:

I have used the VOF method to handle one-phase problem such as generating wave and dam break, considering air is void and getting velocity and pressure of free surface via interpolation. And there are no problem with mass conservation....



In fact I don't understand well what you are doing when you say that you handle one-phase problem and getting velocity and pressure of free surface via interpolation....
Could you explain?

What is the convection scheme you use to compute Fs field?
leflix is offline   Reply With Quote

Old   August 12, 2011, 07:58
Default
  #8
Senior Member
 
Join Date: Aug 2011
Posts: 251
Rep Power: 6
leflix is on a distinguished road
Another question what is the velocity-pressure algorithm that you use?

Do you under-relax the pressure?

Can you display the isocontour 0.5 of the Fs field? It should coincides with the position of the free surface
leflix is offline   Reply With Quote

Old   August 12, 2011, 08:08
Default
  #9
New Member
 
vivien
Join Date: Jan 2011
Posts: 8
Rep Power: 5
zhangweisnoopy is on a distinguished road
Quote:
Originally Posted by leflix View Post

It is clear on this figure that the interface is not well captured.





In fact I don't understand well what you are doing when you say that you handle one-phase problem and getting velocity and pressure of free surface via interpolation....
Could you explain?

What is the convection scheme you use to compute Fs field?
hi!leflix
I have used the VOF to handle one-phase problem considering air is void and getting velocity and pressure of free surface via interpolation. The interopolation is carried out to satisfy the freesurface boundary conditions.That is to say ,the N-S equations are solved in the domain which is full of fluid.
In the VOF method, the reconstruction method is PLIC (Young), while the tracking method is Lagrangian method ( D.Gueyffier, J.Li, A.Nadim, R.Scardovelli, S.Zaleski, Volume-of-fluid interface tracking with smoothed surface stress methods for three-dimensional flows ) instead of Eulerian method (Hirt and Nichols) .
Thank you!
zhangweisnoopy is offline   Reply With Quote

Old   August 12, 2011, 08:25
Default
  #10
New Member
 
vivien
Join Date: Jan 2011
Posts: 8
Rep Power: 5
zhangweisnoopy is on a distinguished road
Quote:
Originally Posted by leflix View Post
Another question what is the velocity-pressure algorithm that you use?

Do you under-relax the pressure?

Can you display the isocontour 0.5 of the Fs field? It should coincides with the position of the free surface
The N-S equations are solved through projection method.
T=200*0.0025s
200t.jpg
T=269*0.0025s
269t.jpg
Thank you!
can you give me you email address so that I can discuss with you ?Thank you
zhangweisnoopy is offline   Reply With Quote

Reply

Tags
incompressible two phases, vof

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



All times are GMT -4. The time now is 13:31.