Hello everyone !
I haven’t been able to find a topic that deals with this weird issue, so here it is:
I’m trying to solve a 3D Navier-Stokes problem using IPCS ( Incremental Pressure Correction Scheme ). It is the same core principle as in Fenics tutorial ft07.
Here is my method:
-First, I create a STEP file with Onshape (here, a cylinder)
-I mesh the STEP file and tag all the boundaries in GMSH and then convert it to XDMF.
-Finally I import the XDMF file in Fenics, compute Navier-Stokes solution and plot the results with Matplotlib and Paraview.
Here is my issue:
I can see no arrows to represent the velocity field in matplotlib, only points. When plotting in Paraview arrows seems to be unit vectors on the x direction, which is very unrealistic.
However, the norm values seems quite relevant: zero at the cylinder’s wall and maximum velocity near the center.
Any idea on why this is happening ?
(dowload the github repository)
Minimal working example:
Script named that is needed to import the create_mesh function:
.msh file needed: