CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

multiphaseEulerFoam for packed bed reactor

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 20, 2012, 16:44
Default multiphaseEulerFoam for packed bed reactor
  #1
New Member
 
Yanxiang Shi
Join Date: Mar 2012
Location: Cambridge, MA, US
Posts: 29
Rep Power: 5
yanxiang is on a distinguished road
Foamers,

I am working on a two phase flow in packed bed reactor. Since the solid phase is not moving, I thought I can basically force it fixed by setting corresponding velocity to zero at each time step. So I did. In the UEqn.C, I added after line 9 so it became

Code:
    volVectorField& U = phase.U();

    if (phase.name() == "water")
    {
        U = U*0;
    }

    volScalarField nuEff(sgsModel->nut() + iter().nu());
It compiled OK. I took the bubbleColumn as the test case (tutorials/multiphase/multiphaseEulerFoam/bubbleColumn) and changed the property of the water phase to that of a solid (taken from tutorials/multiphase/twoPhaseEulerFoam/bed) but left the other phase (gas phase) unchanged. Initial conditions are shown in the attached picture t0.png. alphawater is really alphasolid in this case. Since the solid phase velocity is forced to be zero, I would expect it not moving. However, the solution shows that the solid phase is compressed till the volume fraction is 0.9. The sum of the two phases, indicated by alphasum, is, however, 1 (with slight deviation at some points due to numerical errors). I played around with the diamterModel parameters, it stayed the same.

So why would the solid phase compress itself if the velocity is zero and the alpha equations depend on that? What do I miss?

Thanks,
yanxiang
Attached Images
File Type: jpg t0.jpg (30.1 KB, 46 views)
File Type: jpg t24.jpg (32.6 KB, 33 views)
yanxiang is offline   Reply With Quote

Old   November 20, 2012, 18:45
Default
  #2
New Member
 
Yanxiang Shi
Join Date: Mar 2012
Location: Cambridge, MA, US
Posts: 29
Rep Power: 5
yanxiang is on a distinguished road
I think I figured out what went wrong myself.

So the solution procedures of the alpha equations in multiphaseEulerFoam is similar to the ones presented in Rusche's thesis (Eqn 3.58), only adding more equations. So with this equation, even Ua = 0, alpha_a still changes. Change this to the straightforward solution of the equations should solve the problem, but also would possibly introduce unboundedness.
yanxiang is offline   Reply With Quote

Old   November 21, 2012, 17:11
Default
  #3
New Member
 
Yanxiang Shi
Join Date: Mar 2012
Location: Cambridge, MA, US
Posts: 29
Rep Power: 5
yanxiang is on a distinguished road
Ok. Seems like nobody is really interested in this topic.

So anyway, my last post was wrong. The fact that the solid phase is compressed has nothing to do with how the alpha equations are solved. I implemented the straightforward solution of those, but what I got is simply worse. Compression still, and unboundedness. So for this simple two phase case, the solution is basically not to solve the alpha equations. alpha_solid = alpha_solid0, and alpha_gas = 1- alpha_solid.
yanxiang is offline   Reply With Quote

Old   December 14, 2012, 13:51
Default
  #4
Senior Member
 
Kent Wardle
Join Date: Mar 2009
Location: Illinois, USA
Posts: 195
Rep Power: 10
kwardle is on a distinguished road
Maybe it is not that no one is interested, just not fast enough to respond in 1 day...
If you are only looking at two phases, why not try with twoPhaseEulerFoam? Did you try to change the order that your phases are defined in transportProperties? Perhaps that way you would have the solid phase be the dependent as you have suggested should be the case.
I guess, fundamentally I don't see why multiphaseEulerFoam is the correct approach for you. Is there a reason why a porous approach does not work. Just some thoughts.
kwardle is offline   Reply With Quote

Old   December 14, 2012, 15:36
Default
  #5
New Member
 
Yanxiang Shi
Join Date: Mar 2012
Location: Cambridge, MA, US
Posts: 29
Rep Power: 5
yanxiang is on a distinguished road
Dear Kent,

Thanks for the reply, and for comforting. I wasn't patient enough.

You are right about the twoPhaseEulerFoam. I have discontinued any effort of trying to make the multiphaseEulerFoam work. It is not such a trivial task to do. So instead, I switch to twoPhaseEulerFoam and so far I have made some progress. However, just to clarify, there really is three phases, i.e., a static solid phase and flowing gas and liquid phases.

The porous media model isn't ideal because what I am interested in is the liquid distribution within the bed, whereas the porous media model implemented in OF only considers uniform porosity. I have implemented my own variable porous media model, but I still favor the multi-fluid approach as it provides more flexibility and to some degree is closer to the real physics because it accounts for the interactions between phases.

Now I am facing another issue. There is a brief discussion here:

basic question with 'ForAll' loop

So the issue is basically that I can't seem to limit the volume fraction within 0 and 1. I would appreciate it if you could help me with that.

yanxiang
yanxiang is offline   Reply With Quote

Reply

Tags
multiphase flow

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
Packed Bed Reactor shankar reddy FLUENT 7 March 2, 2014 10:57
urgent... multiphase cfd model of packed bed reactor.( VOF MODEL using) balu@gold6 FLUENT 4 July 26, 2012 10:37
Fixed bed reactor.....you can help...I am sure!! Amrrefay FLUENT 0 June 6, 2009 06:58
Simulating trickle-bed reactor (packed bed) brajesh FLUENT 0 October 17, 2005 06:54
Bubbling Fluidized Bed reactor burning biomass Mattias FLUENT 6 June 8, 2004 02:53


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