Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Selection of relaxation time tau_f in LBM_3D_SinglePhase_Solver.py #28

Open
gxw-ll opened this issue Jul 5, 2024 · 5 comments
Open

Comments

@gxw-ll
Copy link

gxw-ll commented Jul 5, 2024

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
Copy link
Author

gxw-ll commented Jul 5, 2024

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

@yjhp1016
Copy link
Owner

yjhp1016 commented Jul 5, 2024 via email

@gxw-ll
Copy link
Author

gxw-ll commented Jul 5, 2024

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

@yjhp1016
Copy link
Owner

yjhp1016 commented Jul 5, 2024 via email

@gxw-ll
Copy link
Author

gxw-ll commented Jul 5, 2024

这个需要和你使用的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: @.*>

谢谢老师

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants