This code accompanies the paper "Modeling of bubble coalescence and break-up using the Log-normal Method of Moments" by V. Habiyaremye, E.M.J. Komen, J.G.M. Kuerten and E.M.A. Frederix. For a description on what LogMoM does, please refer to that paper.
LogMoM is developed by Edo Frederix and Victor Habiyaremye as part of the PIONEER programme at the Nuclear Research and Consultancy Group (NRG), Westerduinweg 3, 1755 LE Petten, the Netherlands. PIONEER is financed by the Dutch Ministry of Economic Affairs and Climate Policy.
LogMoM is published under the GNU GPL Version 3 license.
LogMoM is distributed under the European Dual Use Codification N: EU DuC=N. Goods labeled with an EU DuC (European Dual-Use Codification) not equal to 'N' are subject to European and national export authorization when exported from the EU and may be subject to national export authorization when exported to another EU country as well. Even without an EU DuC, or with EU DuC 'N', authorization may be required due to the final destination and purpose for which the goods are to be used. No rights may be derived from the specified EU DuC or absence of an EU DuC.
- OpenFOAM-10 Foundation version. While it may compile against other versions, this is not tested and currently not supported.
- Python with numpy, scipy and matplotlib
- Make sure that OpenFOAM-10 is loaded into your environment
- Compile all libraries and apps with
./Allwmake
- All three cases can be run using
./prep.sh multiphaseEulerFoam
- The prep shell scripts set up the case by configuring the required parameters. Some prep scripts require arguments to be specified.
- For example, the TOPFLOW case can also be run in parallel, with
./prep.sh A 16 logmom false mpirun -np 8 multiphaseEulerFoam -parallel
- This runs TOPFLOW case A on a mesh with 16 cells across the radius using logmom without poly-celerity.
- After running the uniformCoalescence and uniformBreakup cases, pdf plots can be generated with
python plot.py
The LogMoM model comes as an additional diameterModel which can be
selected for each phase separately. To use LogMoM in your own cases, you
must first compile LogMoM as described above. This will put a shared
object called libLogMoM.so in your $FOAM_USER_LIBBIN path. Your case
must load this shared object by adding
libs
(
"libLogMoM.so"
);
to your case's system/controlDict file. Next, you can select the
threeMomentLogNormal model as a diameterModel for each phase, in the
constant/phaseProperties file. Note that this works only with the
multiphaseEulerFoam solver (in OpenFOAM-10). Additional parameters, such as
the coalescence and break-up models to use, must be provided as well. An example
is given in cases/TOPFLOW/constant/phaseProperties.LogMoM.m4. Finally, your
case should provide initial and boundary conditions for the N.<phase> and
A.<phase> fields, which are the zeroth and second diameter-based moments of
the size distribution of phase <phase>. Refer to the cases in cases/ for
examples.
- Frederix, E.M.A., Cox, T.L.W., Kuerten, J.G.M., & Komen, E.M.J. (2019). Poly-dispersed modeling of bubbly flow using the log-normal size distribution. Chemical Engineering Science, 201, 237-246.
- Habiyaremye, V., Komen, E.M.J., Kuerten, J.G.M., & Frederix, E.M.A. (2022). Modeling of bubble coalescence and break-up using the Log-normal Method of Moments. Chemical Engineering Science, 253, 117577.
- Habiyaremye, V., Kuerten, J.G.M., & Frederix, E.M.A. (2023). Comparison of population balance models for polydisperse bubbly flow. Chemical Engineering Science, 278, 118932.
- Frederix, E.M.A., Habiyaremye, V., Tajfirooz, S., & Kuerten, J.G.M. (2024). Extension of the two-fluid model to bubble size distribution moment velocities. Experimental and Computational Multiphase Flow, submitted.
For bug reports or support, feel free to contact Edo Frederix at [email protected]. Please note that this code is not maintained nor regularly updated, and is only tested with OpenFOAM-10. Questions related to other versions will thus not be answered.
LogMoM is provided by the copyright holders and contributors "as-is" and any express or implied warranties, including, but not limited to, the implied warranties of merchantability and fitness for a particular purpose are disclaimed. In no event shall the copyright owner or contributors be liable for any direct, indirect, incidental, special, exemplary, or consequential damages (including, but not limited to, procurement of substitute goods or services; loss of use, data, or profits; or business interruption) however caused and on any theory of liability, whether in contract, strict liability, or tort (including negligence or otherwise) arising in any way out of the use of this software, even if advised of the possibility of such damage.