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

[@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.

[@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 ?

[@foojanka3162]

Hi, Many thanks for your thorough explanation! I had a question regarding the units of the output data. For instance, for the single phase permeability example I input my own raw format image and the result of permeability was 0.05. The problem is that I have no idea of the physical units. Could you please help me out on this? I do not know how to do the unit conversion from lattice.

[@TuxRiders]

Hi, I believe it's a better idea if you contact the authors of the software for a more precise answer on your unit conversion question: doi.org/10.1016/j.softx.2022.101097.

[@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...

[@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.

[@kehu1549]

I think it is a lattice unit rather than a physical unit.

[@mazaherhayatdavoudi]

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.

[@mazaherhayatdavoudi]

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.

[@mazaherhayatdavoudi]

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