
[Sponsors] 
December 16, 2009, 16:57 
Rewriting twoPhaseEulerFoam in conservative form

#1 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,895
Rep Power: 26 
Hello,
I started to work on an improved version of a twofluid solver with the goal of dealing with gasparticle flows with the traditional kinetic theory closures for the particle phase and, eventually, in future, more advanced techniques. At first I thought to work on a twophase code, because there still are a few questions to be addressed before adding complexity to the problem. So the steps I plan to follow are:
Currently I have a rough implementation with a fictitious equation of state (I'm assuming the granular temperature to be constant), without the fix that causes the velocity peak at the interface. P.S. Of course the code will be released, as soon as it works
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image. OpenQBMM  An opensource implementation of quadraturebased moment methods Last edited by alberto; December 16, 2009 at 17:44. 

April 22, 2010, 21:30 
Track this work

#2 
Member
Roro Wang
Join Date: Mar 2010
Location: Cambridge, MA, USA
Posts: 30
Rep Power: 7 
Hi Alberto,
This work sounds very attracting. How is it going on? I was flooded by the code. twoPhaseEulerFoam is not easy to be adapted to new turbulence lib structure, especially for a roockie as me. I look forward to your advance and wish all the best. Thanks. Roro 

April 23, 2010, 06:29 

#3 
Senior Member
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 14 
It is interesting!
I did some work to create a RANS multiphase turbulence library in the style of the single phase turbulence libraries. They generalised reasonably well, at least it makes writing new models easier. It would be interesting to compare how we've both done it.
__________________
Laurence R. McGlashan :: Website 

April 23, 2010, 21:14 

#4 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,895
Rep Power: 26 
Hi Laurence and Roro :)
@Laurence: how is it going with your work? This work is going on, a bit slowly due to other obligations, but it's not stuck. I can give a short status update, even if there still are some open questions. Please, keep in mind I'm interested in a code to compute gasparticle flows, so some of the developments cannot be directly generalized to a gasliquid system. In particular, I dropped the pressure equation based on total continuity, and I implemented two pressure equations, one for the fluid and the other for the particle phase. Proper interpolation practises removed the peak at the interface (it's really all already in OpenFOAM, in some of the solvers) The algorithm to enforce the particle packing limit is what gave me more troubles, and there still are some open question on this point to improve the stability of the solver, given the level of stiffness introduced by the equation of state. Best, 

April 27, 2010, 07:12 

#5 
Senior Member
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 14 
Very slowly...
I've got a multifluid solver alongside DQMoM solving a coupled PBE implemented now (multiple dispersed phase momentum equations so that each node is transported with a different velocity) , and it seems to work ok. There are very slight differences between the DQMoM and the QMoM solution, but I suppose that is to be expected, and they should vanish with sufficient refinement. I'm still having problems when solving the countercurrent contactor (due to the phase inversion). I can get a solution which looks reasonable, but I have a lot of validation to do and it's not very stable. There was a big problem with some of the breakup models (at least, with the way I'd implemented them!), so I hope to have that cleared up soon.
__________________
Laurence R. McGlashan :: Website 

April 27, 2010, 12:06 

#6 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,895
Rep Power: 26 
The difference between QMOM and DQMOM is, I would say, OK. Btw, an improved flavour of multivariate QMOM will come out soon
About twoPhaseEulerFoam stability, how does it crash? Alpha going above 1 or drag not converging? Best, Alberto 

April 27, 2010, 12:11 

#7 
Senior Member
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 14 
It's when the coupling between the CFD and the PBE is activated. The CFD obviously does not like having too big/small diameters being returned to the drag force. I can get around it by imposing cut offs but it's rather arbitrary.
__________________
Laurence R. McGlashan :: Website 

April 27, 2010, 12:33 

#8 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,895
Rep Power: 26 
Probably a more implicit treatment of the drag would help. The semiimplicit treatment is known to work in some situation, but in general the partial elimination algorithm (PEA) or a coupled solution of the momentum equations is required for stability purposes.
The PEA is relatively easy to put into OF, the coupled solution requires much more work, since no coupled solver is available in OpenFOAM at the moment. MFIX uses PEA for nphases, and it is very robust though. Best, 

April 27, 2010, 16:27 

#10  
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,895
Rep Power: 26 
Quote:
It is surely an effective solution, but probably also one of the reasons why FLUENT multifluid simulations are not exactly fast, at least compared to what MFIX can do. It is surely possible to use the partial elimination with twophases, and you can derive an incomplete version of the PEA for n phases (See MFIX  www.mfix.org documentation and Mathiesen's work (Prof. Hjertager group). I did something similar in MFIX to couple QMOM for kinetic equations to the fluid solver, and there I've 8 velocities (quadrature abscissas) per each particle phase ). Drop me an email if you want to talk about it. I have some notes and references on the topic. Best, Alberto 

April 28, 2010, 13:15 

#11 
Senior Member
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 14 
I've just written the PEA for a twofluid case, it seems to give much better results, I'll post a derivation and comparison up soon after I've checked it over.
__________________
Laurence R. McGlashan :: Website 

April 28, 2010, 14:36 

#12 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,895
Rep Power: 26 
Yes, if you have small bubbles, with strong coupling, the difference is very clear


September 18, 2010, 06:02 

#13 
Member
Bjorn H. Hjertager
Join Date: Mar 2009
Posts: 69
Rep Power: 8 
Alberto & Laurence
I am very interested to know about the progress in this area. As you migth know, my group has had success in using the multiphase PEA algorithmn in both bubble column and fluidised bed applications also using the DQMOM. We have our own fortran code but would be interested in using OF for thes problems. regards Bjorn 

September 18, 2010, 16:12 

#14  
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,895
Rep Power: 26 
Quote:
you have been one of my teachers at a summer school in Udine in 2006, when I was doing my PhD in Torino. I really appreciated the insights in the multiphase numerics you gave. They were of great help for my work. My focus is on gasparticle flows, and what I did up to now is to modify twoPhaseEulerFoam so that:
About QMOM, a colleague of mine in Dr Fox group, who will graduate next week has been working on a new flavour of QMOM for multivariate distributions (Conditional QMOM  CQMOM) applied to nanoprecipitation. She coupled it with a singlephase solver in OpenFOAM. At the moment, the implementation is not done "the OpenFOAM way", but the code is functional. Some extra updates. QMOM has also been extended to deal with kinetic equation (Boltzmannlike, to be clear), with the goal of directly approximating the kinetic equation to describe polydisperse flows and overcome the limitations of multifluid models. In this method QMOMlike techniques are used to solve for the dispersed phase velocity, without making restrictive assumptions to derive closures for the phase viscosity as in KTGF. This development is not done in OpenFOAM yet, since the main initial goal was to develop the approach without too many complications. However we have a functional code based on MFIX, and the procedure is colocated by nature, so porting it should not be too painful (we are thinking about it, since now the base multiphase code seems to be quite robust). Best, 

September 21, 2010, 04:35 

#15 
Senior Member
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 14 
Dear Bjorn,
I see you are at Stavanger. I lived there and in Bergen until I was 7. I can't remember much apart from the snow though. Yes, I've read the work of Stefano Bove in particular and enjoyed reading Neils Deen's work. I've been working on multiphase DQMoM for gasliquid/liquidliquid applications, where each quadrature node has its own velocity and hence a separate momentum equation is solved for. I also want to increase the number of internal coordinates in the population balance to simulate mass transfer in a real application. I've extended bubbleFoam to solve for N dispersed phases and coded multiphase PEA (hope it's correct  it does give equivalent results to the semiImplicit method with better residuals when the drag force is large). I also have some multiphase turbulence models and some surface forces implemented in a runtime selectable framework. I have problems with solution stability when solving a momentum equation for each quadrature node. The residuals are very oscillatory and some cases diverge completely when coalescence and breakup in the population balances are on. I am not sure how to proceed from here: I may look at using a coupled approach. Any suggestions would be appreciated. Alberto: Regarding the cut off for solving the momentum equations. We had convergence problems when imposing arbitrary cutoffs for reactive flows using DQMoM, and found it better to filter/smooth the source terms towards zero. I'm not sure if it would be relevant (some limited discretisation schemes probably do a similar thing?) but the paper is here.
__________________
Laurence R. McGlashan :: Website 

September 21, 2010, 05:29 

#16 
Member
Bjorn H. Hjertager
Join Date: Mar 2009
Posts: 69
Rep Power: 8 
Dear Laurence,
We also had some problems in solving the PBM using DQMOM in Fluent. You migth look at the PhD thesis of Jesper Madsen where he explains what he did: http://hugin.aue.auc.dk/publ.php?id=2007 Regards Bjorn 

September 21, 2010, 05:35 

#17 
Senior Member
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 14 
Hello,
So you found an operator splitting was necessary? I used that for the reacting flow 'version' of DQMoM, where the chemical time scales are much smaller than that of the flow, we needed to have an operator splitting. Interesting, this should be quick to try. Did you use RADAU5 or CVODE for the time integration? Regards,
__________________
Laurence R. McGlashan :: Website 

September 21, 2010, 10:14 

#18  
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,895
Rep Power: 26 
Quote:
About my problem with the cuttoff, I do not have instabilities related to it, but a certain sensitivity to the cutoff value. In other words the "small" value used as cutoff value cannot be too small (it is typically between 10^4 and 10^6). In some cases 10^6, which is what I would like to use, is too small and gives some trouble at interfaces (wrinkles), while 10^4 works fine, but it might be too high in some case, depending on the application. Thank you for your reference. I'll take a look and let you know 

September 30, 2011, 08:31 

#19 
New Member
Join Date: Sep 2011
Posts: 12
Rep Power: 6 
Hello there.
I've been working with multiphase simulations for a while now. One of the new goals of the taskforce I work with is the evaluate OpenFOAM as a possible substitute ofr other simualtions packages. Would one og guys mind telling me how does the code used in the tutorial example "DamBreak4phase" ensures that the sum of all components fractions will be equal to 1??? After reading your discussion I thought one of you could help me. Thanks a lot 

September 30, 2011, 10:31 

#20  
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,895
Rep Power: 26 
Quote:
You can find the implementation of this into multiphaseInterFoam in the solveAlphas method in multiphaseMixture.C. Best, 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Wind turbine simulation  Saturn  CFX  34  October 16, 2014 05:27 
RPM in Wind Turbine  Pankaj  CFX  9  November 23, 2009 05:05 
conservative, non conservative form????  vijesh joshi  Main CFD Forum  4  March 16, 2006 00:10 
Replace periodic by inletoutlet pair  lego  CFX  3  November 5, 2002 21:09 
Equations in full form  Mark Render  Main CFD Forum  3  June 20, 2000 15:05 