Hello Jochen,
if the local error indicators are rotationally symmetric, then also the
refined grid should be rotationally symmetric. So there must something
go wrong in your code. Perhabs you could provide us with some more
information: Where does it go wrong? Check the error indicators and do
refinement by hand. Are there differences between your refined grid and
the automatically generated one?
The square grid geometry is not a problem. If the boundary conditions
are fine, then this also works for rotationally symmetric solutions.
Best Regards,
Markus
Am 21.04.11 09:36, schrieb Jochen Abhau:
Hi,
I am studying a time dependent PDE on a square domain and want to
compare exactness of the solution to the number of dofs. I use the Kelly
estimator for local refinement. For testing the program, I use a
parameter setting allowing for a rotational symmetric analytic solution
at every time. Now I observe that the locally refined mesh is by far NOT
rotational symmetric, but rather symmetric along the coordinate axes,
although the (analytic) Kelly estimator is rotational symmetric for
rotational symmetric functions.
Now I wonder if anyone has experienced similar problems? Is there an
easy way to get around this problem which seemingly comes from the
square grid geometry?
Right now I am getting correct solutions, but local adaptation works so
poorly that there is nearly no advantage compared to uniform meshes.
Thanks,
Jochen
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
