You always face times that you want to add new physical models to your solver. In two-phase simulations, drag model and turbulence model are perhaps two of the most important ones that one may want to change. In this tutorial you will learn how to add these physical models to the standard solver twoPhaseEulerFoam. Almost the same procedure may be applied for other solvers. Before you start, knowing C++ basics and programming can help you understand the whole procedure better.
- Solver: twoPhaseEulerFoam
- Compatible with: OpenFOAM 7, OpenFOAM 6, OpenFOAM v1912
- Model: New physical models such as RAS turbulence model and drag model
- Note: You must be familiar with C++ programming
- Description: The solver twoPhaseEulerFoam solves a complete set of compressible balance equations including continuity, momentum and energy equations for the two phase system: with one phase as dispersed phase. The model equations are closed by using closure relations for drag force, heat transfer coefficients, turbulence models, and other sub-models. A set of conventional sub-models are already implemented in the solver (see other tutorials in this series). You always face times when we want to test the solver outputs with new sub-models. Here, you will learn how to add these new sub-models to the existing ones: a new drag model and a turbulence model.
That was awesome.
Thanks to CEMF group.
Thanks. Glad it helped.
Hello Prof. Thank you for the tutorial.
However, I think Re should not be multiplied by solid volume fraction (alpha2) in the Difelice.C file. I also think that the return section in DiFelice.C file should not be multiplied by void fraction (max(pair_.continuous(), pair_.continuous().residualAlpha()).
Kindly check and let me know what you think
For DiFelice Model, Re is defined in these papers:
R. Di Felice, The voidage function for fluid-particle interaction systems, Int. J. Multiph. Flow. 20 (1994) 153–159, https://doi.org/10.1016/0301-9322(94)90011-6.
and
https://www.sciencedirect.com/science/article/abs/pii/S0921883118303212
(Effect of drag models on CFD–DEM predictions of bubbling fluidized beds with Geldart D particles)
sorry, alpha2 is void fraction (1-solid void fraction)
It is fluid volume fraction and this multiplication converts local velocity to somehow superficial velocity, since the original Re equation is formulated based on superficial velocity.
Thank for this tutorial. please on implementation of new drag model, i have problem about Make sub-folder. i don’t no why it is empty when a have create it.