Simon,
visit0000.pngThank you for your respone. Sorry, I did not think about the
meaning of the modification in step.21. It was naive. But in my work, I solve
only the Darcy equation (v velocity, p pressure) on a more complicated mesh,
whichis similar to a Hyper_Ball. I have the same problems. So I tried to find
a “simple case” with the same behavior to understandthe problem.
For example I solve the Darcy equation on a "HyperCube" with constant
Dirichlet condition for the pressure (p=p_D on partial Omega). Then I get the
trivial solution v = 0, p = p_D (=2000).
But if I use a “Hyper_Shell” or “Hyper_Ball” instead of a “Hyper_Cube”, then
there are problems(see attachments) and I don't get the trivial solution. Why?
The question is: How exactly do you prescribe the boundary values? There are
several steps associated with getting boundary values right:
* You need to generate a geometry
* If you have curved boundaries, you also need to make sure that the different
parts of the boundary have the appropriate manifold objects associated with
them to ensure that the mesh is refined in the correct way
* You need to associate different boundary indicators to different parts of
the boundary
* You need to write a function that for a given points returns the appropriate
boundary value.
I don't know how you implement all of these steps, but you should double check.
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.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/9672f2da-afca-5552-7f2d-e55c363b8fe4%40colostate.edu.