A question about lbm_ solver_ 3d.py #18
Comments
|
Hi,
For the first question, in the snapshot you shared with me, I can see the
code is printing the result every 1000 time steps?
For the second question, lbm_solver_3d.py was old code to do single phase
simulation, you can ignor it. LBM_3D_SinglePhase_Solver.py is a script
with all single phase kernels (wrapped into classes). You don't need to
study this script if you just need to use LBM to do simulation. If you want
to make a simulation, you only need to write ~20 lines of code. all you
need to do is open example_cavity.py, make a copy and rename it, change
parameters inside (geometry file name, sizes, fluid property etc). Then run
it. We have a documentation with instructions, you can have a look:
https://yjhp1016.github.io/taichi_LBM3D/
…On Sat, 26 Aug 2023 at 10:35, Vivia93 ***@***.***> wrote:
Hello，
When I run lbm_ solver_3d.py,why only output the first step and not the
subsequent steps ?
[image: ***@***.******@***.***%1M]
<https://user-images.githubusercontent.com/130535694/263455946-e0661574-dd3f-4c11-ad18-1c0e6d80a8f3.png>
Also, may I ask what is the difference between lbm_ solver_ 3d.py file and
LBM_3D_SinglePhase_Solver.py file?
—
Reply to this email directly, view it on GitHub
<#18>, or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQB4SKPBSNEXPRYJXKTXXG7PZANCNFSM6AAAAAA37NW6E4>
.
You are receiving this because you are subscribed to this thread.Message
ID: ***@***.***>
|
|
杨教授您好，
我在进行单相流模拟的时候遇到了一点问题，我将自己CT扫描得到的几何模型代入代码，几何模型尺寸为200x100x100，但是输出的图像貌似不对，您可以帮忙看一下是哪里出了问题吗？
输出的图像：
叶落何翩翩
***@***.***
…------------------ 原始邮件 ------------------
发件人: "yjhp1016/taichi_LBM3D" ***@***.***>;
发送时间: 2023年8月26日(星期六) 晚上10:27
***@***.***>;
***@***.******@***.***>;
主题: Re: [yjhp1016/taichi_LBM3D] A question about lbm_ solver_ 3d.py (Issue #18)
Hi,
For the first question, in the snapshot you shared with me, I can see the
code is printing the result every 1000 time steps?
For the second question, lbm_solver_3d.py was old code to do single phase
simulation, you can ignor it. LBM_3D_SinglePhase_Solver.py is a script
with all single phase kernels (wrapped into classes). You don't need to
study this script if you just need to use LBM to do simulation. If you want
to make a simulation, you only need to write ~20 lines of code. all you
need to do is open example_cavity.py, make a copy and rename it, change
parameters inside (geometry file name, sizes, fluid property etc). Then run
it. We have a documentation with instructions, you can have a look:
https://yjhp1016.github.io/taichi_LBM3D/
On Sat, 26 Aug 2023 at 10:35, Vivia93 ***@***.***> wrote:
> Hello，
>
> When I run lbm_ solver_3d.py,why only output the first step and not the
> subsequent steps ?
> [image: ***@***.******@***.***%1M]
> <https://user-images.githubusercontent.com/130535694/263455946-e0661574-dd3f-4c11-ad18-1c0e6d80a8f3.png>
>
> Also, may I ask what is the difference between lbm_ solver_ 3d.py file and
> LBM_3D_SinglePhase_Solver.py file?
>
> —
> Reply to this email directly, view it on GitHub
> <#18>, or unsubscribe
> <https://github.com/notifications/unsubscribe-auth/AJEDKQB4SKPBSNEXPRYJXKTXXG7PZANCNFSM6AAAAAA37NW6E4>
> .
> You are receiving this because you are subscribed to this thread.Message
> ID: ***@***.***>
>
—
Reply to this email directly, view it on GitHub, or unsubscribe.
You are receiving this because you authored the thread.Message ID: ***@***.***>
|
我看不到你的图片，你能再发送一次吗，有可能是你的输入geometry 文件的顺序不太对，如果一个图片转换成0和1， 每张图是X,Y 方向，输入的文件里面的0 和1 的顺序应该是 这个顺序的 |
|
|
|
你用paraview 打开以后好像初始尺寸不太对，paraview 里面显示你的geometry 大小是多少？
另外程序是用一个三层循环去读取的
for k in range (nz):
for j in range (ny):
for i in range (nx)::
read solid[i][j][k]
…On Thu, 9 Nov 2023 at 11:30, Vivia93 ***@***.***> wrote:
[image: W2BVM5~2S1S {F7QLG)$R 3]
<https://user-images.githubusercontent.com/130535694/281715964-46e33c9d-ea20-4980-b24f-6aa5acb2ea13.png>
[image: @2{LWFA7AC_54W35XM}$AG]
<https://user-images.githubusercontent.com/130535694/281715994-90ce9171-6a2a-417d-b005-00a77a41054a.png>
[image: 0(5DH_5J3AW ***@***.***
<https://user-images.githubusercontent.com/130535694/281716027-2ccdce6d-8939-467f-b44f-4697099f3e8e.png>
—
Reply to this email directly, view it on GitHub
<#18 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQAPYT53FDFYWOADQULYDS5FZAVCNFSM6AAAAAA37NW6E6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTQMBTGY2TOOJXGM>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
|
杨教授， |
|
单相流一般就是压力和速度了，你可以画streamline 流线，可能效果更好些
…On Mon, 20 Nov 2023 at 12:25, Vivia93 ***@***.***> wrote:
杨教授，
您好。我看vtk文件输出的可视化只有密度和速度的变化，请问单相流该如何去看流体流动过程呢？
—
Reply to this email directly, view it on GitHub
<#18 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQFINQLH5LL4VT3NNNDYFND3NAVCNFSM6AAAAAA37NW6E6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTQMJYHE3DSOJWGY>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
|
杨教授，您好！我运行了单相流的程序，为什么我运行到24步以后，突然就不行了呢？运行过程中也没有报错。 |
|
那应该是设置上有些问题让simulation crash了。。。
粘度不要太小，尽量使用bodyforce来驱动流体，压力边界使用上需要特别小心，数值上不是很稳定
…On Wed, 6 Dec 2023 at 12:45, Vivia93 ***@***.***> wrote:
杨教授，您好！我运行了单相流的程序，为什么我运行到24步以后，突然就不行了呢？运行过程中也没有报错。
VTV.IZ6SRMLBR.FPN.J7Z5G.png (view on web)
<https://github.com/yjhp1016/taichi_LBM3D/assets/130535694/126a95e9-e265-4178-9d3e-62023f059155>
***@***.*** (view on web)
<https://github.com/yjhp1016/taichi_LBM3D/assets/130535694/adfc463d-592a-4351-9db3-a7f0bfcd210c>
—
Reply to this email directly, view it on GitHub
<#18 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQAEMMX7YCSAIUN6B23YIBSHXAVCNFSM6AAAAAA37NW6E6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTQNBSG44TMOJQGU>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
Hello，
When I run lbm_ solver_3d.py,why only output the first step and not the subsequent steps ?
Also, may I ask what is the difference between lbm_ solver_ 3d.py file and LBM_3D_SinglePhase_Solver.py file?
The text was updated successfully, but these errors were encountered: