From 14e79470fc5cec3a7ec91544b4093e0893851c22 Mon Sep 17 00:00:00 2001 From: kronbichler Date: Sat, 24 Aug 2013 02:02:52 +0000 Subject: [PATCH] Update comment about no normal flux boundary conditions around reentrant corners git-svn-id: https://svn.dealii.org/trunk@30471 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/include/deal.II/numerics/vector_tools.h | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/deal.II/include/deal.II/numerics/vector_tools.h b/deal.II/include/deal.II/numerics/vector_tools.h index f83d28fc42..5cb317c497 100644 --- a/deal.II/include/deal.II/numerics/vector_tools.h +++ b/deal.II/include/deal.II/numerics/vector_tools.h @@ -1529,7 +1529,13 @@ namespace VectorTools * different boundary indicators and call this function twice, once for each * boundary indicators; doing so will yield only one normal vector at this * point per invocation (because we consider only one boundary part at a - * time), with the result that the normal vectors will not be averaged. + * time), with the result that the normal vectors will not be averaged. This + * situation also needs to be taken into account when using this function + * around reentrant corners on Cartesian meshes. If no-normal-flux boundary + * conditions are to be enforced on non-Cartesian meshes around reentrant + * corners, one may even get cycles in the constraints as one will in + * general constrain different components from the two sides. In that case, + * set a no-slip constraint on the reentrant vertex first. * * *

Computing constraints in 3d

-- 2.39.5