Oded,
I tried to solve a slightly modified version of step-22 and obtained strange
results. The only difference in my case with respect to step-22 is the
definition of the boundary conditions. In my case, no Dirichlet boundary
conditions are applied on the velocity. Instead of that, a pressure boundary
condition is applied such that If x =- 2 then p(x) = -2, and if x = 2 then x = 2.
I expected to get a solution showing a pressure gradient along the x-axis, and
a velocity field with a non-positive x-component. However, for some reason,
the solution is different. Attached are the resulting VTK files from running
my code.
> [...]
It would be great to know if someone has an idea what is wrong in my
implementation.
You are asking too much of the people on this mailing list -- namely, to debug
your code. You will need to learn the skills to figure out what the issue is
yourself, because you will need to debug these problems many times over in
your career.
I will note that the Stokes equation does not allow prescribing a pressure. It
allows prescribing a boundary traction (the normal component of the stress). I
have no idea whether that's related to your problem, but it at least seems
like an issue.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
