Jianhui Yang

Thanks for your feedback, no problem, take your time, I can use CUDA for now :-) Kind regards, Jianhui On Fri, 26 Aug 2022 at 08:16, PENGUINLIONG ***@***.***> wrote: >...

Hello, thanks for your interest, yes, I feel the same, it's due to periodic boundary conditions for velocity field and also the capillary force is bigger than inertia force, so...

你好，程序里面的参数都是可以修改的，需要根据你想要的流体情况修改，还有可视化输出频率也是可以修改的。 程序这部分定义了 nx,ny,nz = 131,131,131 #nx,ny,nz = 131,131,131 fx,fy,fz = 5.0e-5,0e-5,0.0 #niu = 0.1 niu_l = 0.1 #psi>0 niu_g = 0.1 #psi wrote: > 亲爱的杨博士， 我是一名多孔介质多相流LBM学习者。你的工作非常出色，我们可以花更多的时间来研究物理问题本身。 > 当我运行文件夹“2phase”中的“lbm_solver_3d_2phase.py”时，发生了一些奇怪的事情。10000步后，在出口处出现侵入相，但此时侵入相并未突破。这是由于使用周期性边界条件时体力较小造成的吗？ > >...

Prediction of capillary hysteresis in a porous material using lattice-Boltzmann methods and comparison to experimental data and a morphological pore network model B Ahrenholz, J Tölke, P Lehmann, A Peters,...

Hi, yes, you can use your own image, you first need to convert your 3D CT image into a text based data file (0=void, 1=solid), with for z in 1..nz:...

Hi, They should be same, do you have more information what are different? porosity? The visualization you showed is a counter-surface plot? On Thu, 21 Sept 2023 at 09:02, hangqqq...

Ah, I see, 输入图像（一般都是8bit的灰度图像）都必须要进行二值化才可以，应为作为LBM只需要有0=孔，1=固体， 所以程序把灰度图像转换成了0 1 矩阵 On Thu, 21 Sept 2023 at 09:39, hangqqq ***@***.***> wrote: > Hi， > txt文档里面几乎数值几乎都是255，但模拟之后显示的contour比例远小于txt中所显示的。 > Almost all values in the txt document are...

porosity around 0.16 looks right, I'm a little confused by the number 10, proportion is a value between 0.0-1.0? why it is 10? On Thu, 21 Sept 2023 at 09:54,...

the counter plot is the iso-surface between 0 and 1. As you can see, due to the openeess at the inlet outlet boundary(not a closed manifold), probably paraview don't know...

You can load all vti files together when you use open files then you can use threshold to select phase field for water or oil only [image: image.png] At the...