CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT > Fluent Multiphase

mass not balance in multiphase simulation

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

View Poll Results: Are you having the similar mass balance issue?
Yes 47 92.16%
No 4 7.84%
Voters: 51. You may not vote on this poll

Like Tree5Likes
  • 2 Post By CeesH
  • 1 Post By BlnPhoenix
  • 2 Post By CeesH

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 13, 2015, 13:11
Unhappy mass not balance in multiphase simulation
  #1
New Member
 
Join Date: Mar 2015
Posts: 3
Rep Power: 11
delphi_d is on a distinguished road
I'm setting a fluidized bed simulation, which has two phases: gas and solid. Eulerian model is used. Geometry is a 2D rectangle with around 5500 square-shaped cells. Gas comes in from the bottom velocity inlet and goes out from the top pressure outlet. Solid is initially patched within part of the bottom area with volume fraction of 0.55.

The problem is that, the total mass of solid phase keeps dropping, with the solid mass flow rates of all the boundaries to be 0. Since I use constant density, the total volume of solid phase is dropping in the same way, while volume of gas phase is increasing.

I'm checking for the total mass of solid phase because I had a chemical reaction in the case, and change of total mass may indicate the progress of reaction. Then the problem happened, and I reduced the case to be without chemical reaction to still find the problem.

Now I have no clue where could those masses go? Or is this just a normal situation?
delphi_d is offline   Reply With Quote

Old   March 19, 2015, 20:20
Default solved
  #2
New Member
 
Join Date: Mar 2015
Posts: 3
Rep Power: 11
delphi_d is on a distinguished road
Still don't know what is wrong.
Good thing is, after I use syamlal-obrien in stead of gidaspow for the drag law, all things turned normal. That's actually a UDF from Fluent tutarial 6.3, which is a modified version of syamlal-obrien.
delphi_d is offline   Reply With Quote

Old   March 24, 2015, 06:55
Default
  #3
Senior Member
 
Cees Haringa
Join Date: May 2013
Location: Delft
Posts: 607
Rep Power: 0
CeesH is on a distinguished road
Hi,

I actually have a similar problem for stirred tank (gas/liquid) - the mass of liquid is not conserved, even though there are no sources for liquid in the domain, and drainage at the top through the pressure outlet does not explain the imbalance (for one, liquid mass is often growing and the outlet-backflow is gas only)

Did you find any reason so far?
OuYangming and Darko like this.
CeesH is offline   Reply With Quote

Old   March 24, 2015, 11:12
Default
  #4
New Member
 
Join Date: Mar 2015
Posts: 3
Rep Power: 11
delphi_d is on a distinguished road
Hi, CeesH

I still don't know what specific reason caused this. But I would guess it's some computing issues rather than model setup problem. Once I tried mixture model rather than Eulerian model for the multiphase, then problem solved. But mixture model is only for extremely low fraction of solid.
The new thing is, after I refine my mesh to 4 times the original cell number, without decreasing time step, that imbalance mass occurs again. So I would suggest you to set a smaller number for you time step, or coarsen your mesh. Even try mixture model maybe, if you are using Eulerian.
delphi_d is offline   Reply With Quote

Old   April 3, 2015, 10:52
Default
  #5
New Member
 
Majid
Join Date: Jan 2014
Location: Canada
Posts: 23
Rep Power: 12
sharif88 is on a distinguished road
Hi guys,

I am also facing the same problem, my case is a static bed placed at the bottom of an annuls, water is coming into the domain from inlet and my goal is to see if the bed begin to erode at a given water flow rate,
the weird thing is that I tried Eulerian-Eulerian model and monitored the total mass of the secondary phase, it keeps dropping continually, I conducted a test with DDPM model I still see the total mass of the secondary phase dropping but when I check the fate of particles they all are still in the domain, how is that possible?

anyways I am going to try your suggestion of changing drag law to see if it helps,

Thanks
sharif88 is offline   Reply With Quote

Old   August 21, 2017, 04:42
Default
  #6
Senior Member
 
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 22
vasava will become famous soon enough
Has any one found solution to this mysteriously disappearing mass?
vasava is offline   Reply With Quote

Old   August 21, 2017, 05:12
Default
  #7
Senior Member
 
Join Date: Aug 2014
Location: Germany
Posts: 292
Rep Power: 13
BlnPhoenix is on a distinguished road
In general in Euler-Euler you need to find a suitable time step size for a given mesh size and make sure the calculation converges in each time step. Also higher order discretization shemes should be used. Once these factors can be ruled out and mass imbalance still persists, deeper investigation into the problem is needed.
soheil_r7 likes this.
BlnPhoenix is offline   Reply With Quote

Old   August 22, 2017, 07:08
Default
  #8
Senior Member
 
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 22
vasava will become famous soon enough
Thanks for the advise.

I am trying a simple case where initially particles are suspended in the air and they are allowed to settle down with time. I am interested in packing pattern once everything is settled down.
vasava is offline   Reply With Quote

Old   August 23, 2017, 07:25
Default
  #9
Senior Member
 
Cees Haringa
Join Date: May 2013
Location: Delft
Posts: 607
Rep Power: 0
CeesH is on a distinguished road
I suppose the issue in my simulations arose from the fixed flowfield approach during mixing. What we did is first solve the flowfield (residuals ~E-4) in steady state, after which mixing was conducted with the transient solver. since no flow periodicities were expecte, the flowfield was frozen (i.e. not solved every timestep), and only the scalar transport equation was solved each step.

Although the residuals were decently low, the equations are of course not exactly satisfied, and by freezing the equations there may be a consistent net. imbalance in momentum locally. Since this is frozen, it consistently exists every timestep, and may act as a "source" of tracer material via the convective terms in the scalar transport equation. Indeed, unfreezing the equations and conducting a fully transient simultation solves the issues - but at a huge time penalty.
soheil_r7 and Darko like this.
CeesH is offline   Reply With Quote

Old   August 25, 2017, 06:12
Default
  #10
Senior Member
 
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 22
vasava will become famous soon enough
Thanks Cees Haringa for sharing your experience.

The case we were trying seems to be working alright now. The continuity residual is still high (around 0.1 with smallest time-step I can afford), but the good thing is that the mass is not disappearing as in earlier cases. The higher continuity is perhaps due to open surface boundary condition where air is allowed to move in-n-out freely or due to coarse mesh.

Anyways, its moving in a positive direction.

Also, when we consulted Ansys about it, they recommended us to switch to DDPM. But since tracking 20 millions of particles could make simulation very long, we are sticking to Eulerian for now.
vasava is offline   Reply With Quote

Old   August 25, 2017, 08:50
Default
  #11
Senior Member
 
Cees Haringa
Join Date: May 2013
Location: Delft
Posts: 607
Rep Power: 0
CeesH is on a distinguished road
Hehe, did they make a quote for a new supercomputer as well to run those ddpm simulations
CeesH is offline   Reply With Quote

Old   November 10, 2017, 23:21
Default
  #12
New Member
 
Hussein
Join Date: Jul 2017
Posts: 7
Rep Power: 8
H.Almrayatee is on a distinguished road
Hi vasava

I'm facing the same problem the mass is continuously decreasing, can you please explain how you are able to solve the problem using Eulerian model

Thank
H.Almrayatee is offline   Reply With Quote

Reply

Tags
fluidized bed, mass balance error, multiphase unsteady

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Using VOF model and Lagrangian Multiphase in the same simulation ganesh_krishnan STAR-CCM+ 2 March 17, 2017 02:34
Multiphase simulation of bubble rising Niru CFX 5 November 25, 2014 14:57
Low Mixing time Problem Mavier CFX 5 April 29, 2013 01:00
mass balance novice FLUENT 0 March 21, 2006 03:33
mass balance of single species antonietta FLUENT 0 November 11, 2005 12:11


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