zh-shen

thanks for your quick reply! I think the result has reached equilibrium, you can see the screenshot 

do i still run longer?

To reduce the computing time, I reduced fy. and this is the code  the formular I used to calculate the max velovity is $u_{max}={gH^2\over 8\nu}$, then the analytical one...

:-)

hi, i recalculated two case of poiseuille flow with your new code the formular that i used to compute the max velocity is $u_{max}={fH^2\over 8 \rho \nu}$ the code of...

but i get the same results evenif if write ti.init(default_ip=ti.i64,default_fp=ti.f64,arch=ti.cpu,kernel_profiler=False, print_ir=False) in the beginning of the code :disappointed_relieved:

:sob:

i want to calculate the permeability of a porous medium, so i need to validate the code first

we measure the permeability of a porous medium in the expriment is 0.4mD, at present, we can't get satisfying simulation results from your code. I don't want to deep into...

thanks for your advice :smile: