Poiseuille flow (steady) in 3D pipe
Dear Dolfin developpers,
Following the example given by the "navier-stokes" testcase on the L-shape domain (from the tutorials), I'm trying to create a Poiseuille testcase on the simple 3D pipe geometry(radius = 1, length = 10).
The testcase is written in C++ and the files (mesh+main+
https:/
I'm prescribing the following boundary conditions:
* Wall = non-slip condition (sub-domain =1)
* Inlet = parabolic velocity (sub-domain = 3)
* Outlet = zero-pressure (sub-domain = 2)
Though the test case is stationnary, I'm using Chorin method (like in demo_navier_stokes) because I intend in the future to prescribe a time-varying profile at the inlet. Nevertheless I don't expect the solution to be different than a Poiseuille flow with this solver.
The time step increment is set to "dt = 0.05" and the final time to "T = 10".
Looking at the results of the simulation, it seems that the velocity profile at the inlet "expands" a bit like a fan (the velocity arrows diverge from the pipe axis), it stabilizes after a while but finally crashes for t = 6.1s (see pictures). Which I really don't understand. I was expecting the Poiseuille flow to develop all the way until the outlet and stay like this until the end of the simulation. Do you have any idea what I'm doing wrong in my files?
Thanks a lot for your help.
Best regards.
Emilie Sauvage.
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- DOLFIN Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask SAUVAGE Emilie for more information if necessary.