On Fri, Jun 21, 2019 at 9:36 AM 서승진 <tjtmdwls1...@gmail.com> wrote:

>  I tried to use the smoother after reading your email, and  it seems that
> the LaplaceMeshSmoother is quitely similar with smooth function in
> MeshTools::Modification namespace.
>
> After trying this, I found a new problem.
>
>  As you said, boundary nodes are not moved.
>
>  However, internal nodes escape from its domain, and as you said, it seems
> that they ingnore the boundary nodes...
>

Are you still running the smoother for 1000 iterations as in the original
code snippet you posted? That may be causing an issue, as I don't think
I've ever tested Laplace smoothing for more than about 2 iterations. Note
that even though boundary nodes are not moved, they should still be
accounted for in the "averaging" process of the smoother, and therefore
nodes should not be able to go outside the domain. Theoretically I think
it's possible for Laplace mesh smoothing to send points outside the domain
for non-convex geometries, though off the top of my head I can't really
construct a case where this would happen for Quad4s... and I don't think
you have a non-convex domain (?) so I'm not really clear on how your
internal nodes are escaping.

--
John

_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to