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

Mass conservation in multiphase problem

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 1, 2010, 09:37
Default Mass conservation in multiphase problem
  #1
cna
New Member
 
candice
Join Date: Mar 2010
Posts: 4
Rep Power: 6
cna is on a distinguished road
Hi all,

I am working on a two-phase flow in a hydrocyclone using the eulerian multiphase model. A mixture air/water is entering at 4m/s the hydrocyclone of a dimension of approximately 1m high. Even when the solution seems to be stable, only approximately 2kg/s on a total of 4kg/s are leaving by the outlets. I tried to decrease the limit of the residues from 0.001 to 0.0001, to increase the discretization orders (momentum and turbulent quantities from 1st to scd order) but this deficit in mass doesn't change a lot.
I used the k-e RNG Swirl dominated flow and the Quick volume fraction discretization.
In the boundary conditions I retain the default turbulence specification method (k-e) instead of the Intensity and Hydraulic diameter and I worked with the 3d version and not the 3ddp.I don't know if these parameters can explain my lack of mass conservation...?
Is there someone having an idea for my problem?
Thanks a lot
cna is offline   Reply With Quote

Old   April 1, 2010, 18:15
Default
  #2
Senior Member
 
karine
Join Date: Nov 2009
Posts: 158
Rep Power: 6
thecfduser is on a distinguished road
Hi
1-Why are u enetring air in a hydroclone??
2-how are u choosing ure dispersed phase dimension????

For ure mass conservation problem:
1- are u using steady or unsteady solver????
2- did u monitor any data ????
3-3d or 3ddp will give u the same results.
4-u are not using QUICK for quantities other than volume fraction??
thecfduser is offline   Reply With Quote

Old   April 2, 2010, 08:42
Default
  #3
cna
New Member
 
candice
Join Date: Mar 2010
Posts: 4
Rep Power: 6
cna is on a distinguished road
The device is a separator (Stairmand), sorry.
I took 1mm of diameter for the water dropplets but without any specific reason.

I am using the unsteady solver.
I monitored the volume fraction and I could see from the distribution of the volume fraction that it didn't change a lot at the end of the simulation (during 8 secondes on a total of 15 secondes), as it is for the mass imbalance.
I used Quick only for the volume fraction.

Thank you for your answers
cna is offline   Reply With Quote

Old   April 2, 2010, 14:45
Default
  #4
Senior Member
 
karine
Join Date: Nov 2009
Posts: 158
Rep Power: 6
thecfduser is on a distinguished road
Hi
are u monitoring mass at the outlet ????
I mean that if u look to the mass at the outlet, it will not be the same as at the inlet because u are using unsteady solver, so u can allow accumulation. But as u say, in average it will be the same. I dont know if ure problem is here...
Anyway, i advice u to use QUICK for all quantities as in a highly swirling flow, a high order discretisation is very important.
Ure time step is how much??? do u have problems with turbulent viscosity ?
thecfduser is offline   Reply With Quote

Old   April 3, 2010, 08:30
Default
  #5
cna
New Member
 
candice
Join Date: Mar 2010
Posts: 4
Rep Power: 6
cna is on a distinguished road
At what it seems to me the end of the simulation I have checked the report fluxes for the mass several times and so I can see that the mass at the outlet doesn't change significantly in order to balance the inlet flux.
In the beginning my time step was of 0,0005s but in the end I could increase it to 0,01 without divergence, I think because the solution was nearly steady-state.
I had some convergence problems when I tried to use second order discretization but not with the 1st order discretization (for momentum and turbulent quantities).
cna is offline   Reply With Quote

Old   April 3, 2010, 18:56
Default
  #6
Senior Member
 
karine
Join Date: Nov 2009
Posts: 158
Rep Power: 6
thecfduser is on a distinguished road
Hi
in fact, as ure flow is unsteady, ure flow at the outlet will never be the same at the inlet. Are u sure that the inlet flow was always superior to ure outlet flow??? try to monitor it to be sure. Remember that in cyclones, u often have a reversed flow (a reversed core) that can extend to the outlet, so most of the time, u can have a small outlet flow, that abruptly get bigger for short moments...
As u say, with second order discretization u can get convergence problem. the way to avoid this is to refine time step. 0,0005 is very good . 0.01 is too big: even if u dont diverge (specially that u are using an implicit time stepping) ure solution of course will not be good. according to my experience, ure velocities can very enormous......
My advices: use QUIck, use a fixed time step. Verify that ure mesh is enough small. A very large mesh (specially with distorsed elements) or a non conformal grid, all this can cause big mass imbalance.
( 1 iteration by time step is sufficiant, but keep ure time step<0,0005)

If u have any of these problems, tell me and i can give more help
thecfduser is offline   Reply With Quote

Old   April 4, 2010, 10:16
Default
  #7
cna
New Member
 
candice
Join Date: Mar 2010
Posts: 4
Rep Power: 6
cna is on a distinguished road
Thank you very much for all your advices, I will apply them!
cna is offline   Reply With Quote

Reply

Tags
eulerian multiphase model, hydrocyclone, mass conservation

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
mass conservation problem with VOF model pat77 FLUENT 2 November 5, 2012 07:07
Mass conservation law violation in stationary case vitmalin CFX 3 September 3, 2007 13:32
Mass Conservation in LES Jaswant Main CFD Forum 8 July 4, 2005 22:32
mass conservation in diverging flow yonghong yan Main CFD Forum 4 July 27, 2002 02:06
problem in mass flow boundary bapi CFX 2 December 3, 2001 22:47


All times are GMT -4. The time now is 04:14.