Wolfgang, thanks, I'm inclined towards using FE_Nothing. Since this becomes quite technical, I'm thinking of a following test. Can you please confirm if this is a sensible way of validation.
(i) Construct a single cell quad mesh. Possibly refine n times. (ii) At the top (arbitrary choice) boundary set no-flux BCs. Set (smooth) Dirichlet BCs on all other sides. (iii) Solve laplace problem. (iv) Construct a two cells quad mesh such that the top boundary of the mesh from (i) became an internal face shared by the two cells. (v) Set cell #1 to FE_Q and cell #2 to FE_Nothing. Refine n times. (vi) Impose no-flux face integral at the shared faces. Set same Dirichlet BCs as in (ii) on all other sides of cell #1. (vii) Solve laplace problem. (viii) Compare solution in cell #1 to the one in step (iii). They should be equal. > I'll add another thought: If you have non-zero Neumann boundary > conditions, > and you approximate the original domain by one that approximates the > boundary > by a stair-step function as you seem to do, then the numerical solution > does > not converge to the exact one. That's because the stair-step boundary has > a > different length than the exact boundary (and doesn't converge to the > exact > length), and the total flux into the domain equals the average flux on the > boundary times the length of the boundary. > An interesting point, although I guess in a row of all approximations I have already done in this problem, this one won't be the first to hit me back. -- 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/8d94125b-501d-421f-874d-0a3700e9062c%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.
