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

A new population balance equation solver for OpenFOAM

Register Blogs Community New Posts Updated Threads Search

Like Tree2Likes
  • 1 Post By sharonyue
  • 1 Post By sujata

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 9, 2020, 00:00
Default A new population balance equation solver for OpenFOAM
  #1
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 849
Rep Power: 18
sharonyue is on a distinguished road
Dear Foamers,

Multiphase system can be simulated by the Eulerian-Eulerian (E-E) method and Eulerian-Lagrangian (E-L) method. The drawback of the traditional two-fluid method is that only one velocity can be employed for both phases (although a multi-fluid method can be constructed). Another promising method is to use Eulerian method for the continuous phase and a generalized population balance equation (GPBE) for the dispersed phase. The GPBE should be solved numerically by existing method, such as the quadrature-based moments method (QBMM). Therefore, we call this method as E-QBMM, although strictly speaking it is based on Euler framework. Recently, we implemented such solver in OpenFOAM and make it open-sourced [1]. This solver was also validated by several gas-liquid test cases [2]. On the other hand, I am very interested into the algorithm. Therefore, I wrote it down and hope it can help others to understand the theory of QBMM [3]. All equations can be found there.

I also would like to share some other interesting works, because the QBMM algorithm was extensively investigated/developed by the research group in Iowa State University, Ecole Polytechnique and others. For example, Alberto Passalacqua also developed his own PBM solver: OpenQBMM [4]. He also started a forum for QBMM users, which looks quite professional..http://openqbmm.freeforums.net/


[1] twoWayGPBEFoam: An open-source Eulerian QBMM solver for monokinetic bubbly flows: https://www.sciencedirect.com/scienc...10465519303728
[2] Comparison of Eulerian QBMM and classical Eulerian–Eulerian method for the simulation of polydisperse bubbly flows: https://aiche.onlinelibrary.wiley.co...1002/aic.16732
[3] Quadrature-based moment methods for the population balance equation: an algorithm review https://www.jianguoyun.com/p/DQwBtyMQ9s3ZBhjDrKYC
[4] https://github.com/OpenQBMM
__________________
My OpenFOAM algorithm website: http://dyfluid.com
By far the largest Chinese CFD-based forum: http://www.cfd-china.com/category/6/openfoam
We provide lots of clusters to Chinese customers, and we are considering to do business overseas: http://dyfluid.com/DMCmodel.html

Last edited by sharonyue; February 9, 2020 at 00:18. Reason: typo
sharonyue is offline   Reply With Quote

Old   February 9, 2020, 12:34
Default
  #2
HPE
Senior Member
 
HPE's Avatar
 
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 931
Rep Power: 13
HPE is on a distinguished road
Thank you sharonyue,

I think OpenQBMM is now part of OpenFOAM-v1912 as a submodule.

Does your code come as a submodule of OpenFOAM as well?

Is it inside OpenQBMM? Or do we need to compile it separately?
HPE is offline   Reply With Quote

Old   February 9, 2020, 22:48
Default
  #3
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 849
Rep Power: 18
sharonyue is on a distinguished road
Hello Herpes,

Our code is not based on OpenQBMM. We implemented it from scratch and it should be compiled separately. I mentioned OpenQBMM is because 1) all of our algorithm belongs to QBMM. 2) I admire their effort to make it open source as well.
HPE likes this.
__________________
My OpenFOAM algorithm website: http://dyfluid.com
By far the largest Chinese CFD-based forum: http://www.cfd-china.com/category/6/openfoam
We provide lots of clusters to Chinese customers, and we are considering to do business overseas: http://dyfluid.com/DMCmodel.html
sharonyue is offline   Reply With Quote

Old   April 19, 2020, 07:23
Default
  #4
New Member
 
Masoud Kahnooji
Join Date: Sep 2017
Location: Iran
Posts: 5
Rep Power: 9
masoudbme90 is on a distinguished road
Hi Sharonyue
Thank you for your information.

What is the advantage of E-QBMM methods in comparison with two-phase euler-euler method (multi-fluid method) like twophaseeulerfoam?
masoudbme90 is offline   Reply With Quote

Old   August 28, 2020, 11:12
Default Population balance in Openfoam 4.1
  #5
New Member
 
sujata
Join Date: Dec 2019
Posts: 10
Rep Power: 7
sujata is on a distinguished road
I am a new user of Openfoam and currently using reactingmultiphaseEulerfoam 4.1. I want to add Population Balance. How can I do that? Any lead would be appreciated.
Ahyar likes this.
sujata is offline   Reply With Quote

Old   July 4, 2021, 04:48
Default
  #6
Member
 
Muhammad Ahyar
Join Date: Mar 2020
Posts: 30
Rep Power: 6
Ahyar is on a distinguished road
Quote:
Originally Posted by sujata View Post
I am a new user of Openfoam and currently using reactingmultiphaseEulerfoam 4.1. I want to add Population Balance. How can I do that? Any lead would be appreciated.
I have the same plan for my project, have you found out how to set up a case with Population Balance into the openfoam case? thanks in advance
Ahyar is offline   Reply With Quote

Old   November 25, 2021, 21:10
Default
  #7
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 849
Rep Power: 18
sharonyue is on a distinguished road
Quote:
Originally Posted by sujata View Post
I am a new user of Openfoam and currently using reactingmultiphaseEulerfoam 4.1. I want to add Population Balance. How can I do that? Any lead would be appreciated.

Hello there, sorry for my late reply. The latest OpenFOAM provides a sectional method. U can have a look.
__________________
My OpenFOAM algorithm website: http://dyfluid.com
By far the largest Chinese CFD-based forum: http://www.cfd-china.com/category/6/openfoam
We provide lots of clusters to Chinese customers, and we are considering to do business overseas: http://dyfluid.com/DMCmodel.html
sharonyue is offline   Reply With Quote

Old   November 25, 2021, 21:13
Default
  #8
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 849
Rep Power: 18
sharonyue is on a distinguished road
We also implemented the QBMM for gas-particle flow, please see the following link:

https://authors.elsevier.com/c/1e8bP7QGmEkfF
[Implementation of CHyQMOM in OpenFOAM for the simulation of non-equilibrium gas-particle flows under one-way and two-way coupling]
__________________
My OpenFOAM algorithm website: http://dyfluid.com
By far the largest Chinese CFD-based forum: http://www.cfd-china.com/category/6/openfoam
We provide lots of clusters to Chinese customers, and we are considering to do business overseas: http://dyfluid.com/DMCmodel.html
sharonyue is offline   Reply With Quote

Old   November 26, 2021, 07:59
Default
  #9
Senior Member
 
Mark Olesen
Join Date: Mar 2009
Location: https://olesenm.github.io/
Posts: 1,715
Rep Power: 40
olesen has a spectacular aura aboutolesen has a spectacular aura about
Quote:
Originally Posted by sharonyue View Post
We also implemented the QBMM for gas-particle flow, please see the following link:

https://authors.elsevier.com/c/1e8bP7QGmEkfF
[Implementation of CHyQMOM in OpenFOAM for the simulation of non-equilibrium gas-particle flows under one-way and two-way coupling]

How much overlap would this have with OpenQBMM? Where possible it would make sense to combine efforts.
olesen is offline   Reply With Quote

Old   January 4, 2022, 01:50
Default plotting the population of bubbles
  #10
Member
 
Muhammad Ahyar
Join Date: Mar 2020
Posts: 30
Rep Power: 6
Ahyar is on a distinguished road
Quote:
Originally Posted by sharonyue View Post
Hello there, sorry for my late reply. The latest OpenFOAM provides a sectional method. U can have a look.
do you know how to plot the bubble population in the new openfoam version? i want to plot the bubble population from the population balance tutorial that mentioned in this article https://openfoam.org/guides/populati...ance-openfoam/

some sources said that to plot the population i can use the NDF plot but i have no idea to do it with openfoam output data
Ahyar is offline   Reply With Quote

Old   January 24, 2022, 07:53
Default
  #11
Senior Member
 
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 354
Rep Power: 16
Santiago is on a distinguished road
Quote:
Originally Posted by sharonyue View Post
Dear Foamers,

Multiphase system can be simulated by the Eulerian-Eulerian (E-E) method and Eulerian-Lagrangian (E-L) method. The drawback of the traditional two-fluid method is that only one velocity can be employed for both phases (although a multi-fluid method can be constructed). Another promising method is to use Eulerian method for the continuous phase and a generalized population balance equation (GPBE) for the dispersed phase. The GPBE should be solved numerically by existing method, such as the quadrature-based moments method (QBMM). Therefore, we call this method as E-QBMM, although strictly speaking it is based on Euler framework. Recently, we implemented such solver in OpenFOAM and make it open-sourced [1]. This solver was also validated by several gas-liquid test cases [2]. On the other hand, I am very interested into the algorithm. Therefore, I wrote it down and hope it can help others to understand the theory of QBMM [3]. All equations can be found there.

I also would like to share some other interesting works, because the QBMM algorithm was extensively investigated/developed by the research group in Iowa State University, Ecole Polytechnique and others. For example, Alberto Passalacqua also developed his own PBM solver: OpenQBMM [4]. He also started a forum for QBMM users, which looks quite professional..http://openqbmm.freeforums.net/


[1] twoWayGPBEFoam: An open-source Eulerian QBMM solver for monokinetic bubbly flows: https://www.sciencedirect.com/scienc...10465519303728
[2] Comparison of Eulerian QBMM and classical Eulerian–Eulerian method for the simulation of polydisperse bubbly flows: https://aiche.onlinelibrary.wiley.co...1002/aic.16732
[3] Quadrature-based moment methods for the population balance equation: an algorithm review https://www.jianguoyun.com/p/DQwBtyMQ9s3ZBhjDrKYC
[4] https://github.com/OpenQBMM
Do you have a repo for your code (github, other)? Would you mind sharing it?
Santiago is offline   Reply With Quote

Old   January 26, 2022, 03:52
Default Regarding the solver OF version compatibility
  #12
New Member
 
Akshay
Join Date: Jan 2020
Posts: 28
Rep Power: 6
enthusiast is on a distinguished road
Quote:
Originally Posted by sharonyue View Post
Dear Foamers,

Multiphase system can be simulated by the Eulerian-Eulerian (E-E) method and Eulerian-Lagrangian (E-L) method. The drawback of the traditional two-fluid method is that only one velocity can be employed for both phases (although a multi-fluid method can be constructed). Another promising method is to use Eulerian method for the continuous phase and a generalized population balance equation (GPBE) for the dispersed phase. The GPBE should be solved numerically by existing method, such as the quadrature-based moments method (QBMM). Therefore, we call this method as E-QBMM, although strictly speaking it is based on Euler framework. Recently, we implemented such solver in OpenFOAM and make it open-sourced [1]. This solver was also validated by several gas-liquid test cases [2]. On the other hand, I am very interested into the algorithm. Therefore, I wrote it down and hope it can help others to understand the theory of QBMM [3]. All equations can be found there.

I also would like to share some other interesting works, because the QBMM algorithm was extensively investigated/developed by the research group in Iowa State University, Ecole Polytechnique and others. For example, Alberto Passalacqua also developed his own PBM solver: OpenQBMM [4]. He also started a forum for QBMM users, which looks quite professional..http://openqbmm.freeforums.net/


[1] twoWayGPBEFoam: An open-source Eulerian QBMM solver for monokinetic bubbly flows: https://www.sciencedirect.com/scienc...10465519303728
[2] Comparison of Eulerian QBMM and classical Eulerian–Eulerian method for the simulation of polydisperse bubbly flows: https://aiche.onlinelibrary.wiley.co...1002/aic.16732
[3] Quadrature-based moment methods for the population balance equation: an algorithm review https://www.jianguoyun.com/p/DQwBtyMQ9s3ZBhjDrKYC
[4] https://github.com/OpenQBMM

Can the twoWayGPBEfoam be compiled in open foam version 9?

is it for openfoam.org or openfoam.com software?
enthusiast is offline   Reply With Quote

Old   February 2, 2022, 08:23
Default
  #13
New Member
 
Akshay
Join Date: Jan 2020
Posts: 28
Rep Power: 6
enthusiast is on a distinguished road
Quote:
Originally Posted by sharonyue View Post
Dear Foamers,

Multiphase system can be simulated by the Eulerian-Eulerian (E-E) method and Eulerian-Lagrangian (E-L) method. The drawback of the traditional two-fluid method is that only one velocity can be employed for both phases (although a multi-fluid method can be constructed). Another promising method is to use Eulerian method for the continuous phase and a generalized population balance equation (GPBE) for the dispersed phase. The GPBE should be solved numerically by existing method, such as the quadrature-based moments method (QBMM). Therefore, we call this method as E-QBMM, although strictly speaking it is based on Euler framework. Recently, we implemented such solver in OpenFOAM and make it open-sourced [1]. This solver was also validated by several gas-liquid test cases [2]. On the other hand, I am very interested into the algorithm. Therefore, I wrote it down and hope it can help others to understand the theory of QBMM [3]. All equations can be found there.

I also would like to share some other interesting works, because the QBMM algorithm was extensively investigated/developed by the research group in Iowa State University, Ecole Polytechnique and others. For example, Alberto Passalacqua also developed his own PBM solver: OpenQBMM [4]. He also started a forum for QBMM users, which looks quite professional..http://openqbmm.freeforums.net/


[1] twoWayGPBEFoam: An open-source Eulerian QBMM solver for monokinetic bubbly flows: https://www.sciencedirect.com/scienc...10465519303728
[2] Comparison of Eulerian QBMM and classical Eulerian–Eulerian method for the simulation of polydisperse bubbly flows: https://aiche.onlinelibrary.wiley.co...1002/aic.16732
[3] Quadrature-based moment methods for the population balance equation: an algorithm review https://www.jianguoyun.com/p/DQwBtyMQ9s3ZBhjDrKYC
[4] https://github.com/OpenQBMM

I am getting an error when I tried to complile the code in OpenFoam v7.

In file included from oneWayGPBEFoam.C:64:0:
pEqn.H: In function ‘int main(int, char**)’:
pEqn.H:29:21: error: ‘class Foam::fvMesh’ has no member named ‘solver’; did you mean ‘solvers_’?
pEqn.solve(mesh.solver(p.select(piso.finalInnerIte r())));
^~~~~~
solvers_
/opt/openfoam7/wmake/rules/General/transform:25: recipe for target 'Make/linux64GccDPInt32Opt/oneWayGPBEFoam.o' failed
make[1]: *** [Make/linux64GccDPInt32Opt/oneWayGPBEFoam.o] Error 1
make[1]: Leaving directory '/mnt/d/solver/twoWayGPBEFoam/solvers/oneWayGPBEFoam'
/opt/openfoam7/wmake/makefiles/apps:39: recipe for target 'oneWayGPBEFoam' failed
make: *** [oneWayGPBEFoam] Error 2

I tried to compile this in /mnt/d/customsolver

Was this the issue? Should I take this folder in the solver directory and then run ./Allwmake
enthusiast is offline   Reply With Quote

Old   February 2, 2022, 12:36
Default Problem Sorted
  #14
New Member
 
Akshay
Join Date: Jan 2020
Posts: 28
Rep Power: 6
enthusiast is on a distinguished road
Quote:
Originally Posted by enthusiast View Post
I am getting an error when I tried to complile the code in OpenFoam v7.

In file included from oneWayGPBEFoam.C:64:0:
pEqn.H: In function ‘int main(int, char**)’:
pEqn.H:29:21: error: ‘class Foam::fvMesh’ has no member named ‘solver’; did you mean ‘solvers_’?
pEqn.solve(mesh.solver(p.select(piso.finalInnerIte r())));
^~~~~~
solvers_
/opt/openfoam7/wmake/rules/General/transform:25: recipe for target 'Make/linux64GccDPInt32Opt/oneWayGPBEFoam.o' failed
make[1]: *** [Make/linux64GccDPInt32Opt/oneWayGPBEFoam.o] Error 1
make[1]: Leaving directory '/mnt/d/solver/twoWayGPBEFoam/solvers/oneWayGPBEFoam'
/opt/openfoam7/wmake/makefiles/apps:39: recipe for target 'oneWayGPBEFoam' failed
make: *** [oneWayGPBEFoam] Error 2

I tried to compile this in /mnt/d/customsolver

Was this the issue? Should I take this folder in the solver directory and then run ./Allwmake
I was trying to compile this in the Windows WSL Linux extension (Ubuntu 18.04 LTS). But it gave me an error, so I switch to OF v6. It compiled successfully.
enthusiast is offline   Reply With Quote

Reply


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
suGWFoam: Richards equation solver for porous media flows liu OpenFOAM Announcements from Other Sources 4 February 10, 2021 16:14
Poisson equation in custom solver diverges for large boundary values babakflame OpenFOAM Running, Solving & CFD 2 May 23, 2017 21:11
Population Balance Equation in OpenFOAM Amjad Asad OpenFOAM Programming & Development 0 January 19, 2016 06:14
Star cd es-ice solver error ernarasimman STAR-CD 2 September 12, 2014 01:01
Population Balance Theory student Main CFD Forum 2 August 10, 2008 05:12


All times are GMT -4. The time now is 15:50.