improvements in Navier-Stokes example
Hi all,
I started playing around with Dolfin's Navier--Stokes capabilities and used http://
Unfortunately, starting with the given example code, it was quite hard to achieve. For example, I tried to model a box with no inflow or outflow, just the force of gravity acting. You'd expect no flow, but replacing the respective boundary conditions in the demo example and setting f=(0,-1) does give something. This is because the boundary conditions on the tentative velocity step and the pressure update don't correctly incorporate f.
Other than that, I noticed that one doesn't really need to solve an equation system for the final step u1: Just setting
u1 -= k * grad(p1) should suffice.
I would give that example a revision unless that is something you rather want to keep it the way it is.
--Nico
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 Nico Schlömer for more information if necessary.