
[Sponsors] 
phase diameter and phase compressibility in multiphaseEulerFoam solver 

LinkBack  Thread Tools  Display Modes 
March 9, 2016, 07:06 
phase diameter and phase compressibility in multiphaseEulerFoam solver

#1 
Member
Michael Frank
Join Date: Aug 2012
Posts: 45
Rep Power: 6 
Hello,
I have recently started using OpenFOAM. Specifically, I am trying to model a multiphase/multicomponent bubble columnlike problem. It seems to me that the multiphaseEulerFoam solver is perfectly suitable for this task. However, I have a couple of questions. First, for the calculation of the interphase interactions, the Eulerian model assumes that each phase forms bubbles of a constant diameter. In the absence of any experimental data, is there any empirical model that we can use to calculate this diameter? My second question regards the compressibility of each phase. The documentation states that the solver is for compressible fluid phases (http://cfd.direct/openfoam/userguide/standardsolvers/). Glancing at the code however, I can see that the energy (or temperature) equation is not solved. Am I missing something? Finally, commercial software packages often consider a multiphase model called the mixture model (or at least that is what ANSYS is calling it). This is much more efficient than its Eulerian counterpart and might be suitable for my case. Does OpenFOAM have something equivalent? Thanks in advance for any help Mike 

March 9, 2016, 17:02 

#2 
Senior Member
Kent Wardle
Join Date: Mar 2009
Location: Illinois, USA
Posts: 208
Rep Power: 13 
Mike,
1. You are correct that the solver as released can handle a constant diameter only. There is a diameterModel class to add in other models and there is one that can handle size changes of bubbles due to external pressure changes. There are numerous equations in literature for estimating bubble size from properties and flow conditions. Not sure on accuracy. Probably more reliable is taking a case which the bubble size is 'known' to start. 2. multiphaseEulerFoam does not include compressibility. Depending on the version of OF you are using, twoPhaseEulerFoam might. 3. One of the main features of multiphaseEulerFoam is to enable eulereuler multiphase (momentum equation per phase with interphase drag) along with a VOFstyle sharpened interface where this is needed. Fluent's 'mixture model' is a single shared momentum equation with a slip velocity to account for drag between phases. This is only appropriate for certain kinds of problems as the Fluent manual can tell you. Hope that helps. Kent 

March 14, 2016, 06:13 

#3 
Member
Michael Frank
Join Date: Aug 2012
Posts: 45
Rep Power: 6 
Hi Kent,
Thanks for your reply. It is indeed very useful. Regarding the diameter model, I see that I can use the isothermalDiameter model. To use this do I simply specify it in the transportProperties dictionary? Looking at the code I understand that there is no need for programming and recompiling the solver. As for the compressibility, is it possible to adjust the code to be compressible (i.e. following a similar approach to the twoPhaseEuler solver) or is this a much more complicated task? Thanks again for your help Mike 

March 14, 2016, 10:03 

#4 
Senior Member
Kent Wardle
Join Date: Mar 2009
Location: Illinois, USA
Posts: 208
Rep Power: 13 
Addition of compressibility would be much more involved as it requires insertion of density as a variable and rearranging how the PDE is written and solved. You can take a look at an earlier version of twoPhaseEulerFoam and compare with recent to get an idea.
It is reasonably straightforward to implement a Boussinesq approximation density variation into multiphaseEulerFoam to account for natural convection flows where density variation throughout the system is not small but the difference is critical to the physics. 

March 15, 2016, 19:18 

#5 
Senior Member
Join Date: Jan 2013
Posts: 310
Rep Power: 6 
My understanding for the compressibility included in multiphaseEulerFoam is because they include "p = p_rgh + rho*gh;". the hydrodynamics and the thermodynamics are not decoupled. So the compression would also lead to the change of the density. However, since the temperature is solved in that solver, so the change of T would not.


March 16, 2016, 04:27 

#6 
Member
Juho Peltola
Join Date: Mar 2009
Location: Finland
Posts: 83
Rep Power: 9 
You may also want to check out the new reactingTwo/MultiphaseEulerFoam solvers:
https://github.com/OpenFOAM/OpenFOAM...ctingEulerFoam 

March 16, 2016, 19:47 

#7 
Member
Michael Frank
Join Date: Aug 2012
Posts: 45
Rep Power: 6 
Thank you all for your help.
openfoammaofnepo, so the multiphaseEulerFoam solver is in fact compressible? Thinking about it, the isothermal diameter model must also include compressibility effects since the bubbles expand and contract with pressure. Juho the solver you suggested looks very interesting. Does it have any documentation? 

Tags 
compressible flow, multiphase, multiphaseeulerfoam, phase diameter 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
[multiphase] C_UDMI in DEFINE_PROPERTY for diameter of secondary phase  SJSW  Fluent UDF and Scheme Programming  1  December 12, 2015 13:06 
Eulerian phase gas diameter  hwet  Fluent Multiphase  3  August 14, 2015 05:59 
phase names in multiphaseEulerFoam: strange behaviour  kobelak  OpenFOAM Running, Solving & CFD  0  May 11, 2015 18:35 
Derive fixed number of phase from clone () in multiphaseEulerFoam Solver  Aj Nair  OpenFOAM Programming & Development  0  December 16, 2013 21:41 
compressible two phase flow in CFX4.4  youngan  CFX  0  July 1, 2003 23:32 