yjhp1016 / taichi_LBM3D

A 3D sparse LBM solver implemented using Taichi
MIT License
267 stars 36 forks source link

Selection of relaxation time tau_f in LBM_3D_SinglePhase_Solver.py #28

Open gxw-ll opened 2 months ago

gxw-ll commented 2 months ago

Hello, Professor! I have encountered a problem with the example you provided in img_ftb1131.txt. The file LBM_3D_SinglePhase_Solver.py provides two tau_f. image when I use the self.tau_f=3.0*self.niu+0.5 The permeability results are as follows: image when I use the self.tau_f=self.niu/3.0+0.5 The permeability results are as follows: image

gxw-ll commented 2 months ago

How should tau_f be selected in permeability modelling of porous media?

yjhp1016 commented 2 months ago

Comment out code means it is out of date, you can delete them

Sam Ultraman @.***> 于 2024年7月5日周五 上午4:44写道:

How should tau_f be selected in permeability modelling of porous media?

— Reply to this email directly, view it on GitHub https://github.com/yjhp1016/taichi_LBM3D/issues/28#issuecomment-2210019804, or unsubscribe https://github.com/notifications/unsubscribe-auth/AJEDKQBJ4DND2M62E52I5RDZKYJCLAVCNFSM6AAAAABKMOCGIKVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDEMJQGAYTSOBQGQ . You are receiving this because you are subscribed to this thread.Message ID: @.***>

gxw-ll commented 2 months ago

谢谢老师!我看见您在lbm_solver_3d_sparse.py使用的tau_f=3.0niu+0.5 image 有些文章在三维情况也使用的tau_f=3.0niu+0.5,我们应该在什么条件下选择tau_f=niu/3+0.5,什么条件下选择tau_f=3.0*niu+0.5呢

yjhp1016 commented 2 months ago

这个需要和你使用的MRT的M矩阵匹配,所以不同文章的会有些不同

On Fri, 5 Jul 2024 at 09:45, xiwei gao @.***> wrote:

谢谢老师!我看见您在lbm_solver_3d_sparse.py使用的tau_f=3.0

niu+0.5 image.png (view on web) https://github.com/yjhp1016/taichi_LBM3D/assets/45071471/41056e03-fe2d-41e4-9f79-6c4d97de6430 有些文章在三维情况也使用的tau_f=3.0 niu+0.5,我们应该在什么条件下选择tau_f=niu/3+0.5,什么条件下选择tau_f=3.0*niu+0.5呢

— Reply to this email directly, view it on GitHub https://github.com/yjhp1016/taichi_LBM3D/issues/28#issuecomment-2210460489, or unsubscribe https://github.com/notifications/unsubscribe-auth/AJEDKQA7R4VYQRQEMCHT2B3ZKZMKBAVCNFSM6AAAAABKMOCGIKVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDEMJQGQ3DANBYHE . You are receiving this because you commented.Message ID: @.***>

gxw-ll commented 2 months ago

这个需要和你使用的MRT的M矩阵匹配,所以不同文章的会有些不同 On Fri, 5 Jul 2024 at 09:45, xiwei gao @.> wrote: 谢谢老师!我看见您在lbm_solver_3d_sparse.py使用的tau_f=3.0 niu+0.5 image.png (view on web) https://github.com/yjhp1016/taichi_LBM3D/assets/45071471/41056e03-fe2d-41e4-9f79-6c4d97de6430 有些文章在三维情况也使用的tau_f=3.0 niu+0.5,我们应该在什么条件下选择tau_f=niu/3+0.5,什么条件下选择tau_f=3.0niu+0.5呢 — Reply to this email directly, view it on GitHub <#28 (comment)>, or unsubscribe https://github.com/notifications/unsubscribe-auth/AJEDKQA7R4VYQRQEMCHT2B3ZKZMKBAVCNFSM6AAAAABKMOCGIKVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDEMJQGQ3DANBYHE . You are receiving this because you commented.Message ID: **@.***>

谢谢老师