Asking for help about Taichi-LBM3D #7
Comments
|
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 it looks like there is some spontaneous
imbibition for blue phase. I remember that I implement a forced-recoloring
at the inlet/outlet to force the color of fluid become blue/red.
you can try to set these two lines like this to force the inlet as blue and
outlet as red: (also probably you need to lower surface tension and/or
increase body force to allow the fluid going to the geometry a little bit
faster) 😀
bc_psi_x_left, psi_x_left = 1, -1.0 # boundary condition for
phase-field: 0 = periodic,
bc_psi_x_right, psi_x_right = 1, 1.0 # 1 = constant value on the
boundary, value = -1.0 phase1 or 1.0 = phase 2
…On Tue, 21 Mar 2023 at 05:15, Ffffffffffly ***@***.***> wrote:
Dear Dr.Yang,
I am a LBM learner about multiphase flow in porous media. Your work is so
excellent, and we can spend more time in the physicis problem itself.
When I run the 'lbm_solver_3d_2phase.py' in the folder '2phase', something
strange happens. The invading phase is appeared at the outlet after 10000
steps, but the invading phase doesnot breakthrough at that moment. Is this
caused by the small body force when using periodic boudary condition?
[image: image]
<https://user-images.githubusercontent.com/92718998/226521430-a994a64a-26ad-4cd2-995b-47330932cfb1.png>
—
Reply to this email directly, view it on GitHub
<#7>, or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQCVBKWQFZGFEIRWGZDW5E2QNANCNFSM6AAAAAAWB6R5N4>
.
You are receiving this because you are subscribed to this thread.Message
ID: ***@***.***>
|
|
Wow, thanks for your reply! It really works! The strange phenomenon seems like counter-current imbibition due to poor driving force along the injecting direction. This forced recoloring strategy is very useful, like the modified periodic condition proposed by Liu. |
您好，当运行文件夹“2phase”中的“lbm_solver_3d_2phase. py”后，参数都是默认的吗？结果的展示是在paraview里面看的，但是只有8个不同形态，请问参数这些都是默认的吗还是有什么别的参数设置或者是我查看的方式有问题啊？我这边没有形成这样的流动形态。谢谢您！ |
|
你好，程序里面的参数都是可以修改的，需要根据你想要的流体情况修改，还有可视化输出频率也是可以修改的。
程序这部分定义了
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<0
psi_solid = 0.7
CapA = 0.005
这部分定义 程序输入geometry 大小，bodyforce, 两种流体的黏度，接触角大小（接触角的cos值）和表面张力
#Boundary condition mode: 0=periodic, 1= fix pressure, 2=fix velocity;
boundary pressure value (rho); boundary velocity value for vx,vy,vz
bc_x_left, rho_bcxl, vx_bcxl, vy_bcxl, vz_bcxl = 0, 1.0, 0.0e-5, 0.0,
0.0 #Boundary
x-axis left side
bc_x_right, rho_bcxr, vx_bcxr, vy_bcxr, vz_bcxr = 0, 0.995, 0.0, 0.0,
0.0 #Boundary
x-axis left side
bc_y_left, rho_bcyl, vx_bcyl, vy_bcyl, vz_bcyl = 0, 1.0, 0.0, 0.0, 0.0
#Boundary
x-axis left side
bc_y_right, rho_bcyr, vx_bcyr, vy_bcyr, vz_bcyr = 0, 1.0, 0.0, 0.0,
0.0 #Boundary
x-axis left side
bc_z_left, rho_bczl, vx_bczl, vy_bczl, vz_bczl = 0, 1.0, 0.0, 0.0, 0.0
#Boundary
x-axis left side
bc_z_right, rho_bczr, vx_bczr, vy_bczr, vz_bczr = 0, 1.0, 0.0, 0.0,
0.0 #Boundary
x-axis left side
bc_psi_x_left, psi_x_left = 1, -1.0 # boundary condition for
phase-field: 0 = periodic,
bc_psi_x_right, psi_x_right = 0, 1.0 # 1 = constant value on the
boundary, value = -1.0 phase1 or 1.0 = phase 2
bc_psi_y_left, psi_y_left = 0, 1.0
bc_psi_y_right, psi_y_right = 0, 1.0
bc_psi_z_left, psi_z_left = 0, 1.0
bc_psi_z_right, psi_z_right = 0, 1.0
这部分是边界条件，定义了是不是使用周期边界，还是压力边界还是速度边界，建议使用周期边界搭配bodyforce这样最稳定
下面是定义是不是固定颜色在边界处，比如这个是X左边固定为 non-wetting phase, 所以你看到的右边的有non-wetting phase
出现应该就是周期边界产生的
可以通过把bc_psi_x_right, psi_x_right = 1, 1.0
来让出口也给固定为 wetting phase
for iter in range(80000+1):
这行定义了整个程序要运行多少步
if (iter%500==0):
这行定义了屏幕输出频率，这个是每 500步在屏幕输出一个summary
if (iter%10000==0):
这行定义了用于paraview 可视化的vtr文件输出频率，这个例子是每10000步输出一个
上面的数据都是可以根据你的需要更改的
…On Sat, 7 Oct 2023 at 06:43, renxiaosa00 ***@***.***> wrote:
亲爱的杨博士， 我是一名多孔介质多相流LBM学习者。你的工作非常出色，我们可以花更多的时间来研究物理问题本身。
当我运行文件夹“2phase”中的“lbm_solver_3d_2phase.py”时，发生了一些奇怪的事情。10000步后，在出口处出现侵入相，但此时侵入相并未突破。这是由于使用周期性边界条件时体力较小造成的吗？
[image: 图像]
<https://user-images.githubusercontent.com/92718998/226521430-a994a64a-26ad-4cd2-995b-47330932cfb1.png>
您好，当运行文件夹“2phase”中的“lbm_solver_3d_2phase.
py”后，参数都是默认的吗？结果的展示是在paraview里面看的，但是只有8个不同形态，请问参数这些都是默认的吗还是有什么别的参数设置或者是我查看的方式有问题啊？我这边没有形成这样的流动形态。谢谢您！
—
Reply to this email directly, view it on GitHub
<#7 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQFALBFG6NBYAK6EUILX6DTZXAVCNFSM6AAAAAAWB6R5N6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONJRGYYTAMJWHE>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
|
Prediction of capillary hysteresis in a porous material using
lattice-Boltzmann methods and comparison to experimental data and a
morphological pore network model
<https://scholar.google.com/scholar?oi=bibs&cluster=10617077284738717967&btnI=1&hl=en>
B Ahrenholz, J Tölke, P Lehmann, A Peters, A Kaestner… - Advances in Water
Resources, 2008
…On Thu, 21 Mar 2024 at 08:05, Ffffffffffly ***@***.***> wrote:
杨老师您好，请问程序中转换矩阵M的参考文献是哪一篇，方便分享一下吗？感谢您！
—
Reply to this email directly, view it on GitHub
<#7 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQD3OLANZVNSSPACJQLYZKIC5AVCNFSM6AAAAAAWB6R5N6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDAMJRGUYTQNBWGU>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
|
杨老师，您好。我最近在做一些LBM计算稳态相渗的工作，阅读了您的博士论文，有些困惑想要请教您。 按照您提出的方法，首先开展排驱和渗吸模拟，获得一定饱和度下的初始相分布，再施加体积力和周期边界，待稳定后根据两相达西计算相渗。（博士论文P116：“After the permeability of both phases converged, the relative permeability for a desired saturation was obtained.”） 我的疑问是如何判断稳定？比如非湿相nw处于低饱和度时，nw不能连续横跨多孔介质两端，以液滴方式流动，那么就会经历卡断、界面破碎与重合，势必会对流场造成影响。我觉得最理想的情况是流场会根据nw流动状态呈现出周期性波动。那么这种情况如何判断稳态条件，如何计算稳定后的Q？ 谢谢杨老师，期待您的回复！ @yjhp1016 |
|
判断稳定这个事确实有一定的主观性，在我处理的样本里面，均质性比较好的砂岩，基本上跑一段时间你如果画出湿润相和非湿润相的速度随时间变化曲线就能看出来基本变化很小了，就是稳定了，但是也有些尺度跨度很大的碳酸盐，跑了很久还是能看到速度曲线一直波动很大，如果能看出来事周期性变化，那也可以认为事稳定了，如果不是，可能就还得继续让流动继续一段时间，或者可视化下油水分布，看看是不是分布合理，人工判断下。这个和样本的复杂度是高度相关的
…On Wed, 26 Jun 2024 at 15:19, Ffffffffffly ***@***.***> wrote:
杨老师，您好。我最近在做一些LBM计算稳态相渗的工作，阅读了您的博士论文，有些困惑想要请教您。
按照您提出的方法，首先开展排驱和渗吸模拟，获得一定饱和度下的初始相分布，再施加体积力和周期边界，待稳定后根据两相达西计算相渗。（博士论文P116：“After
the permeability of both phases converged, the relative permeability for a
desired saturation was obtained.”）
我的疑问是如何判断稳定？比如非湿相nw处于低饱和度时，nw不能连续横跨多孔介质两端，以液滴方式流动，那么就会经历卡断、界面破碎与重合，势必会对流场造成影响。我觉得最理想的情况是流场会根据nw流动状态呈现出周期性波动。那么这种情况如何判断稳态条件，如何计算稳定后的Q？
谢谢杨老师，期待您的回复！ @yjhp1016 <https://github.com/yjhp1016>
—
Reply to this email directly, view it on GitHub
<#7 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQAO6CFOSCSY6YHJBE3ZJLEXHAVCNFSM6AAAAAAWB6R5N6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCOJRHAZTKOJXGY>
.
You are receiving this because you were mentioned.Message ID:
***@***.***>
|
谢谢您的解答！还有一个问题，既然LBM计算稳态相渗十分耗时，为什么不直接模拟实验室非稳态相渗测试呢（随驱替饱和度变化计算）？是JBN方法忽略毛细力的假设不适合孔尺度流动吗？我只看到过一篇边驱边计算的文章“Direct simulations of two-phase flow on micro-CT images of porous media and upscaling of pore-scale forces”，是通过能量计算出达西力再用两相达西方程算RP，也不是通过JBN |
|
JBN
是实际生成种经常使用的方法，好处就是省时省力，但是我们认为稳态更准确更能反应由于是油气运移初期的正确分布等等。确定就是测量非常费时费力。可根据需要选择
…On Wed, 26 Jun 2024 at 15:49, Ffffffffffly ***@***.***> wrote:
判断稳定这个事确实有一定的主观性，在我处理的样本里面，均质性比较好的砂岩，基本上跑一段时间你如果画出湿润相和非湿润相的速度随时间变化曲线就能看出来基本变化很小了，就是稳定了，但是也有些尺度跨度很大的碳酸盐，跑了很久还是能看到速度曲线一直波动很大，如果能看出来事周期性变化，那也可以认为事稳定了，如果不是，可能就还得继续让流动继续一段时间，或者可视化下油水分布，看看是不是分布合理，人工判断下。这个和样本的复杂度是高度相关的
… <#m_7876738697947640361_>
On Wed, 26 Jun 2024 at 15:19, Ffffffffffly *@*.*> wrote:
杨老师，您好。我最近在做一些LBM计算稳态相渗的工作，阅读了您的博士论文，有些困惑想要请教您。
按照您提出的方法，首先开展排驱和渗吸模拟，获得一定饱和度下的初始相分布，再施加体积力和周期边界，待稳定后根据两相达西计算相渗。（博士论文P116：“After
the permeability of both phases converged, the relative permeability for a
desired saturation was obtained.”）
我的疑问是如何判断稳定？比如非湿相nw处于低饱和度时，nw不能连续横跨多孔介质两端，以液滴方式流动，那么就会经历卡断、界面破碎与重合，势必会对流场造成影响。我觉得最理想的情况是流场会根据nw流动状态呈现出周期性波动。那么这种情况如何判断稳态条件，如何计算稳定后的Q？
谢谢杨老师，期待您的回复！ @yjhp1016 <https://github.com/yjhp1016>
https://github.com/yjhp1016 <https://github.com/yjhp1016> — Reply to this
email directly, view it on GitHub <#7 (comment)
<#7 (comment)>>,
or unsubscribe
https://github.com/notifications/unsubscribe-auth/AJEDKQAO6CFOSCSY6YHJBE3ZJLEXHAVCNFSM6AAAAAAWB6R5N6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCOJRHAZTKOJXGY
<https://github.com/notifications/unsubscribe-auth/AJEDKQAO6CFOSCSY6YHJBE3ZJLEXHAVCNFSM6AAAAAAWB6R5N6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCOJRHAZTKOJXGY>
. You are receiving this because you were mentioned.Message ID: @.*>
谢谢您的解答！还有一个问题，既然LBM计算稳态相渗十分耗时，为什么不直接模拟实验室非稳态相渗测试呢（随驱替饱和度变化计算）？是JBN方法忽略毛细力的假设不适合孔尺度流动吗？我只看到过一篇边驱边计算的文章“Direct
simulations of two-phase flow on micro-CT images of porous media and
upscaling of pore-scale forces”，是通过能量计算出达西力再用两相达西方程算RP，也不是通过JBN
—
Reply to this email directly, view it on GitHub
<#7 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQGXB6IOIXR4OG5FTVDZJLIIZAVCNFSM6AAAAAAWB6R5N6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCOJRHEYDIMRUHE>
.
You are receiving this because you were mentioned.Message ID:
***@***.***>
|
老师，您好，我主要从事含水合物沉积物气-水两相流动方面研究，很感谢您开发了这么优秀的程序，看了您的文章“A Single-Phase and Multiphase Lattice Boltzmann Solver on Cross-Platform Multicore CPU/GPUs”之后，想学习更多内容，请问在哪去查找您的博士论文，我想系统学习一下您开发这套程序的过程 |
|
… On Wed, 3 Jul 2024 at 13:58, Sam Ultraman ***@***.***> wrote:
JBN
是实际生成种经常使用的方法，好处就是省时省力，但是我们认为稳态更准确更能反应由于是油气运移初期的正确分布等等。确定就是测量非常费时费力。可根据需要选择
… <#m_-2608658186854581103_>
On Wed, 26 Jun 2024 at 15:49, Ffffffffffly *@*.*> wrote:
判断稳定这个事确实有一定的主观性，在我处理的样本里面，均质性比较好的砂岩，基本上跑一段时间你如果画出湿润相和非湿润相的速度随时间变化曲线就能看出来基本变化很小了，就是稳定了，但是也有些尺度跨度很大的碳酸盐，跑了很久还是能看到速度曲线一直波动很大，如果能看出来事周期性变化，那也可以认为事稳定了，如果不是，可能就还得继续让流动继续一段时间，或者可视化下油水分布，看看是不是分布合理，人工判断下。这个和样本的复杂度是高度相关的
… <#m_7876738697947640361_> On Wed, 26 Jun 2024 at 15:19, Ffffffffffly @.>
wrote: 杨老师，您好。我最近在做一些LBM计算稳态相渗的工作，阅读了您的博士论文，有些困惑想要请教您。
按照您提出的方法，首先开展排驱和渗吸模拟，获得一定饱和度下的初始相分布，再施加体积力和周期边界，待稳定后根据两相达西计算相渗。（博士论文P116：“After
the permeability of both phases converged, the relative permeability for a
desired saturation was obtained.”）
我的疑问是如何判断稳定？比如非湿相nw处于低饱和度时，nw不能连续横跨多孔介质两端，以液滴方式流动，那么就会经历卡断、界面破碎与重合，势必会对流场造成影响。我觉得最理想的情况是流场会根据nw流动状态呈现出周期性波动。那么这种情况如何判断稳态条件，如何计算稳定后的Q？
谢谢杨老师，期待您的回复！ @yjhp1016 <https://github.com/yjhp1016>
https://github.com/yjhp1016 <https://github.com/yjhp1016>
https://github.com/yjhp1016 <https://github.com/yjhp1016>
https://github.com/yjhp1016 <https://github.com/yjhp1016> — Reply to this
email directly, view it on GitHub <#7
<#7> (comment) <#7 (comment)
<#7 (comment)>>>,
or unsubscribe
https://github.com/notifications/unsubscribe-auth/AJEDKQAO6CFOSCSY6YHJBE3ZJLEXHAVCNFSM6AAAAAAWB6R5N6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCOJRHAZTKOJXGY
<https://github.com/notifications/unsubscribe-auth/AJEDKQAO6CFOSCSY6YHJBE3ZJLEXHAVCNFSM6AAAAAAWB6R5N6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCOJRHAZTKOJXGY>
https://github.com/notifications/unsubscribe-auth/AJEDKQAO6CFOSCSY6YHJBE3ZJLEXHAVCNFSM6AAAAAAWB6R5N6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCOJRHAZTKOJXGY
<https://github.com/notifications/unsubscribe-auth/AJEDKQAO6CFOSCSY6YHJBE3ZJLEXHAVCNFSM6AAAAAAWB6R5N6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCOJRHAZTKOJXGY>
. You are receiving this because you were mentioned.Message ID: @.>
谢谢您的解答！还有一个问题，既然LBM计算稳态相渗十分耗时，为什么不直接模拟实验室非稳态相渗测试呢（随驱替饱和度变化计算）？是JBN方法忽略毛细力的假设不适合孔尺度流动吗？我只看到过一篇边驱边计算的文章“Direct
simulations of two-phase flow on micro-CT images of porous media and
upscaling of pore-scale forces”，是通过能量计算出达西力再用两相达西方程算RP，也不是通过JBN — Reply to
this email directly, view it on GitHub <#7 (comment)
<#7 (comment)>>,
or unsubscribe
https://github.com/notifications/unsubscribe-auth/AJEDKQGXB6IOIXR4OG5FTVDZJLIIZAVCNFSM6AAAAAAWB6R5N6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCOJRHEYDIMRUHE
<https://github.com/notifications/unsubscribe-auth/AJEDKQGXB6IOIXR4OG5FTVDZJLIIZAVCNFSM6AAAAAAWB6R5N6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCOJRHEYDIMRUHE>
. You are receiving this because you were mentioned.Message ID: @.***>
老师，您好，我主要从事含水合物沉积物气-水两相流动方面研究，很感谢您开发了这么优秀的程序，看了您的文章“A Single-Phase and
Multiphase Lattice Boltzmann Solver on Cross-Platform Multicore
CPU/GPUs”之后，想学习更多内容，请问在哪去查找您的博士论文，我想系统学习一下您开发这套程序的过程
—
Reply to this email directly, view it on GitHub
<#7 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQFSHYS32VMA6IFAKBLZKPYRFAVCNFSM6AAAAAAWB6R5N6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDEMBWGAYTSOJQG4>
.
You are receiving this because you were mentioned.Message ID:
***@***.***>
|
谢谢老师！ |
Dear Dr.Yang,
I am a LBM learner about multiphase flow in porous media. Your work is so excellent, and we can spend more time in the physicis problem itself.
When I run the 'lbm_solver_3d_2phase.py' in the folder '2phase', something strange happens. The invading phase is appeared at the outlet after 10000 steps, but the invading phase doesnot breakthrough at that moment. Is this caused by the small body force when using periodic boudary condition?
The text was updated successfully, but these errors were encountered: