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

PLIC-VOF with a Conservative formulation at high density ratio

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 9, 2018, 05:54
Post PLIC-VOF with a Conservative formulation at high density ratio
  #1
New Member
 
Abhishek
Join Date: Feb 2018
Posts: 4
Rep Power: 8
abhishek13 is on a distinguished road
Greetings!
I am trying to develop a Two Phase VOF solver with a conservative formulation. I have successfully implemented the following:
  • PLIC-VOF in 2D (for coupling with a Naviers Stokes Solver)
  • Non-conservative formulation for a dam-break problem with density ratio of 1000
  • Conservative formulation for a dam-break problem with denstiy ratio of 10
At higher density ratios, while using the conservative formulation, there is an unphysical increase in the velocity near the interface. I know that this is caused by a sudden jump in density which causes an ill-posed matrix while solving for pressure. I would like to know how to rectify the issue to make the conservative solver work at higher density ratios.
abhishek13 is offline   Reply With Quote

Old   February 9, 2018, 06:36
Default
  #2
Member
 
Join Date: Dec 2012
Posts: 92
Rep Power: 13
beer is on a distinguished road
Basically what you have to do is not to use the average density for the interface density in the mass flux. Rather use a mass flux, which is calculated from the mass fluxes of the two phases:
if the face densities for the two phases are \rho_{f1} and \rho_{f2}, the face normal velocity is v_f, the face are is A, the overall volume flux on the face is

\dot{V}_{Sum}=v_F*A.

Now, the volume flux for the second phase \dot{V}_2 is known from the PLIC reconstruction, because that is what you need PLIC for. This means you can calculate the volume flux for the first phase from these numbers:

\dot{V}_1=\dot{V}_{Sum} - \dot{V}_2

Now just use these two volume fluxes to calculate the mass flux using the density for the two phases on the face seperately:

\dot{m}_f=\dot{V}_{f1}*\rho_1 + \dot{V}_2*\rho_{f2}

This procedure worked very well for me. I can't give you a reference though, because I kind of came up with it myself. However, I am fairly sure someone has done and published it before.

Last edited by beer; February 9, 2018 at 11:27.
beer is offline   Reply With Quote

Old   February 9, 2018, 09:46
Default
  #3
New Member
 
Abhishek
Join Date: Feb 2018
Posts: 4
Rep Power: 8
abhishek13 is on a distinguished road
Hello! Thanks for the reply. Are you suggesting that the density at the faces should be calculated using the amount of both fluids present in the volume (face velocity*dt)*A?
abhishek13 is offline   Reply With Quote

Old   February 9, 2018, 11:28
Default
  #4
Member
 
Join Date: Dec 2012
Posts: 92
Rep Power: 13
beer is on a distinguished road
No, sorry I should have stated that. You should calculate the face density for each phase seperately. I edited my answer above.
beer is offline   Reply With Quote

Old   February 11, 2018, 23:44
Default
  #5
New Member
 
Abhishek
Join Date: Feb 2018
Posts: 4
Rep Power: 8
abhishek13 is on a distinguished road
Sorry but i still dont get what you are trying to say. Where do i need to use the mass flux?

Last edited by abhishek13; February 13, 2018 at 03:46.
abhishek13 is offline   Reply With Quote

Old   February 12, 2018, 04:38
Default
  #6
New Member
 
nkalkote's Avatar
 
Nikhil
Join Date: Mar 2016
Location: India
Posts: 15
Rep Power: 10
nkalkote is on a distinguished road
Dear Abhishek,

This is caused due to disparity in Gamma (specific heat ratio calculation) Please refer the following paper,

https://www.sciencedirect.com/scienc...21999115005021

titled "Preventing numerical oscillations in the flux-split based finite difference method for compressible flows with discontinuities".

As this paper is highly condensed you might have to cross refer the references for more details.


Regards,
nkalkote is offline   Reply With Quote

Old   February 15, 2018, 05:07
Default
  #7
Member
 
Join Date: Dec 2012
Posts: 92
Rep Power: 13
beer is on a distinguished road
Quote:
Originally Posted by abhishek13 View Post
Sorry but i still dont get what you are trying to say. Where do i need to use the mass flux?
For the discretization of the convection term and for the source term in the pressure correction equation.
beer is offline   Reply With Quote

Old   February 15, 2018, 11:23
Default
  #8
Senior Member
 
Saideep
Join Date: Apr 2015
Location: INDIA
Posts: 203
Rep Power: 12
Saideep is on a distinguished road
Might be of use. A simple way is to filter the volume fraction several times so that you dont see a sharp jump in the fluid properties but a gradual jump. This is something that the software "Gerris Flow Solver" using PLIC-VOF does.

However, I did find that for some density ratios (air water) keeping the viscosity ratio as 1, the solver really struggles.
Saideep is offline   Reply With Quote

Reply

Tags
conservative form, twophase fluid, volume of fluid


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
pressure eq. "converges" after few time steps maddalena OpenFOAM Running, Solving & CFD 69 July 21, 2011 08:42
Problem on high density ratio in Level Set method Kai Yan Main CFD Forum 10 December 25, 2007 07:12
conservative total pressure and density Atit Koonsrisuk CFX 0 June 24, 2007 13:45
VOF formulation mvee FLUENT 3 March 12, 2007 07:35
Specified Operating Density for Compressible VOF Andrew Wick FLUENT 1 November 14, 2005 05:58


All times are GMT -4. The time now is 09:07.