Felix,
> Is the rate of convergence as expected? >> > > I believe so. I did a linear regression of the convergence graphs in > logscale and got a value of : > -7.6 for the speed : I believe it corresponds to the -8 from second order > finete element in 2D > -3.2 for the pressure : I beliebe it corresponds to the -4 from the first > order finite element in 2D > > Good! > >> If u = v = 0 in your test case, there might still be a problem with these >> equations. >> > > I'm not sur to understand what you are saying here. Can you elaborate why > you think there is a problem in my equations ? > If the reference solution is 0, you can't see if you are missing a constant factor in your equation even when the rate of convergence looks correct. > > >> Since the ansatz space on non-affine mapped geometries is not polynomial >> anymore, I would also try a cartesian mesh for the case that the solution >> is contained >> in your ansatz space. >> > > I don't really understand what you mean by ansatz space although I have > looked it up on google. Do you suggest that I try the X=0 solution on a > square mesh ? > Yes, but then you also have to be careful with the boundary conditions. If you have an analytical solution with homogeneous Dirichlet boundary conditions for \phi, I would definitely try that first. The ansatz space is the space spanned by the basis functions you are using for your solution variables. In case the mapping from the reference cell to the real cell is affine and you use polynomial base functions (on the reference element), the ansatz space is polynomial as well. Then, it is easier to manufacture a solution that is contained in the ansatz space. Best, Daniel -- 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 [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CAOYDWb%2BJz7DUmMGdN3YkbdubHsOEws91D_cT6REMfGL5LJRGXg%40mail.gmail.com.
