Lattice Boltzmann simulation of multiphase flow in porous media using Palabos [On Demand 13]

Multiphase flow in porous media is a complex phenomenon with numerous applications in various fields. Lattice Boltzmann method has proven to be a valuable tool for studying multiphase flow phenomena, enabling insights into transport processes and aiding in the design and optimization of porous media applications. In the previous video, we had a look at the Palabos lattice Boltzmann code, which includes numerical schemes to simulate multiphase flow relatively easily. In this video, we have a look at these features and the MPLBM-UT tool built on top of them to simulate multiphase flow in porous media.
Video series:
Please refer to the first video of the “On Demand” series to find out more about the purpose of these videos. You can find the list of videos belonging to this series at tuxriders.com/videos/on-demand/ and • On Demand
Materials:
You can find the materials developed in the videos in this repository on GitHub: github.com/TuxRiders/on-demand
Topics covered:
🎯 Examples of porous media flow in Palabos
🎯 MPLBM-UT software as a helper for multiphase porous media flow in Palabos
🎯 Simulating relative permeability in porous media
Lecturer: Mojtaba Barzegari mbarzegary.github.io/
To learn more about the goals of the TuxRiders project, please visit our website at tuxriders.com.
Chapters in this video!
################
00:00 - Introduction
01:48 - Porous media flow in Palabos
03:52 - Introducing MPLBM-UT software
05:25 - Installing MPLBM-UT
09:39 - Exploring porous media examples
10:55 - Running steady-state relative permeability example
15:07 - Visualizing multiphase flow results in ParaView

[@TuxRiders]

I'm terribly sorry. there is a mistake in the video at 7:47. THE PYTHON VERSION SHOULD BE 3.8 and not 3.7.

[@cellian8488]

Thanks a lot. It is really helpful for my research

[@TuxRiders]

Glad it was helpful :)

[@cellian8488]

By the way, is there any way to simulate flow with a large density ratio (like a water-air system)? I tried a density less than 1 for wetting fluid (fluid 2), but it didn't work...

[@md.irfankhan9139]

I am new in this field and was trying hard to understand this. You explained really well. Thank you so much.

[@TuxRiders]

Glad it was helpful :)

[@md.irfankhan9139]

Hi, I am trying to measure permeability of a 3D porous media (200x200x200) using single phase example. How can i change the fluid properties like density, viscosity of the simulation ? and i I take the domain size 100x100x100 will it give me accurate result ? can you help me in this regard ?

[@zq88]

there is also an open source package called LBPM, if you are interested and convenient, can you explain how to install and use it ?

[@TuxRiders]

good suggestion. yes, I used it a bit so will make a video on it soon.

[@zq88]

@@TuxRiders Thank you sir! Looking forward to it!!!

[@khaledahmed2070]

I really thank you for the explanation of this helpful tool. I still could not figure out the units of the Single-phase permeability obtained from this model. Is it in mD or D?

[@TuxRiders]

I'm sorry, I totally forgot to answer your question. I think it should be Darcy as the code uses SI units. I hope I'm not wrong here.

[@user-cb8lj4nn4f]

is it possible to run this simulation in colaboratory? I tried but not being able to execute it.

[@TuxRiders]

I'm not a collab user, so my answer may not be correct. I guess it should not be possible to run it there because this Python codes rely on C++ codes running behind the scene, the compilation of which can be an issue as well.

[@user-cb8lj4nn4f]

I tried to implement all steps. everything went correct except while running "python 2_phase_sim.py" command, I come to this error: "Traceback (most recent call last):
File "/content/drive/MyDrive/Palabos/MPLBM-UT/examples/steady_state_relperm/2_phase_sim.py", line 241, in
inputs = organize_2_phase_outputs_for_relperm(inputs, sim_counter)
File "/content/drive/MyDrive/Palabos/MPLBM-UT/examples/steady_state_relperm/2_phase_sim.py", line 178, in organize_2_phase_outputs_for_relperm
last_rho_config = rho_files[-1]
IndexError: list index out of range". would you help me solve this error? thanks!

[@TuxRiders]

it seems that the simulations cannot run, which is maybe due to some MPI problems. I specified 6 as the number of cores (processes) in the config file in this video, but you may need to decrease it according to the the available resources in your system. try to run it with 2 or 1 first.

[@user-cb8lj4nn4f]

I am trying to run the code in colaboratory. I have done all steps. even the the code is executed and shows "Running sim 8 of 8
PoreSpy Snw = 0.9999889559101527". but all of a sudden the error "Traceback (most recent call last):
File "/content/drive/MyDrive/Palabos/MPLBM-UT/examples/steady_state_relperm/2_phase_sim.py", line 241, in
inputs = organize_2_phase_outputs_for_relperm(inputs, sim_counter)
File "/content/drive/MyDrive/Palabos/MPLBM-UT/examples/steady_state_relperm/2_phase_sim.py", line 178, in organize_2_phase_outputs_for_relperm
last_rho_config = rho_files[-1]
IndexError: list index out of range" comes to me. I do not know what the error is. I would appreciate if you could help me. Thanks@@TuxRiders