Hi Alex,
Yes, at first glance it does look to me like the integration is a problem. You
are using the same FEFaceValues object to perform the integration in both the
material and spatial configurations, and this is not quite correct with the way
that you have things set up. The following
Wolfgang,
Thanks very much for taking time to walk through the code. I was hoping
this would work as it looks like a neat way, at least to me! When I was
posting this question, I also thought why dealii still doesn't have some
functionality like this which makes few things easy. But I was not
Hi Jean-Paul,
Thanks for your response. The spatial normal vectors are now correct. It
seems my lack of grounding in differential geometry was the culprit here.
However, the integrated forces are still not as I expect, but I think this
is a separate issue, so I will create a new thread.