|
[Sponsors] |
PLIC-VOF with a Conservative formulation at high density ratio |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 9, 2018, 05:54 |
PLIC-VOF with a Conservative formulation at high density ratio
|
#1 |
New Member
Abhishek
Join Date: Feb 2018
Posts: 4
Rep Power: 8 |
Greetings!
I am trying to develop a Two Phase VOF solver with a conservative formulation. I have successfully implemented the following:
|
|
February 9, 2018, 06:36 |
|
#2 |
Member
Join Date: Dec 2012
Posts: 92
Rep Power: 13 |
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 and , the face normal velocity is , the face are is , the overall volume flux on the face is . Now, the volume flux for the second phase 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: Now just use these two volume fluxes to calculate the mass flux using the density for the two phases on the face seperately: 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. |
|
February 9, 2018, 09:46 |
|
#3 |
New Member
Abhishek
Join Date: Feb 2018
Posts: 4
Rep Power: 8 |
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?
|
|
February 9, 2018, 11:28 |
|
#4 |
Member
Join Date: Dec 2012
Posts: 92
Rep Power: 13 |
No, sorry I should have stated that. You should calculate the face density for each phase seperately. I edited my answer above.
|
|
February 11, 2018, 23:44 |
|
#5 |
New Member
Abhishek
Join Date: Feb 2018
Posts: 4
Rep Power: 8 |
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. |
|
February 12, 2018, 04:38 |
|
#6 |
New Member
Nikhil
Join Date: Mar 2016
Location: India
Posts: 15
Rep Power: 10 |
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, |
|
February 15, 2018, 05:07 |
|
#7 |
Member
Join Date: Dec 2012
Posts: 92
Rep Power: 13 |
||
February 15, 2018, 11:23 |
|
#8 |
Senior Member
Saideep
Join Date: Apr 2015
Location: INDIA
Posts: 203
Rep Power: 12 |
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. |
|
Tags |
conservative form, twophase fluid, volume of fluid |
|
|
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 |