I am now updating my 2D navier stokes equation with moving boundary to 3D, here is the code Parallel_3d_NS_moving.edp (4.4 KB) Th.mesh (734.3 KB).
However, I meet a error report when I run it for some time as
I do not know the the exact issue, but I am not sure you can use convect when running FreeFEM with PETSc in parallel. Maybe looking at the solution can also provide the insight.
Yes, after reading this post Convect( ) in parallel - #2 by prj. I know that using ‘convect’ inside a ‘varf’ will generate wrong results for high timesteps. But it seems work at small time step as this post Convect( ) in parallel - #6 by mojtaba.barzegari shows. Anyway, I will change ‘convect’ in ‘varf’ with ‘convectParallel’ to see what will happen.
Many thanks for your reply.