On 10/11/25 10:36, Arthur Bawin wrote:
A first general question :
- The Lagrange multiplier should be defined only on the cylinder
boundary, a codimension 1 space. What is the simplest way of dealing
with a <dim-1> finite element space within a FESystem that contains the
other spaces of dimension dim ? Currently, and to use a single FESystem
and DOFHandler, I define lambda as a P2 field over the whole domain and
constrain all of its non-cylinder DOFs to zero, which works but is not
very elegant.
Perhaps not very elegant, but that's how we generally do it. You may
want to compare with how ASPECT does it:
https://github.com/geodynamics/aspect/blob/main/source/mesh_deformation/interface.cc
I believe that FENothing could be used to this end, but
when I tried it seems that it is not fully implemented for simplices, is
that correct ?
Perhaps. That's probably worth fixing. If you can create small test
cases that illustrate what doesn't work, it shouldn't be very difficult
to get those fixed.
The main, more niche question :
[...]
I think you'll have to debug that yourself :-) The usual suggestion
applies as always: make it simple. Try a test case with very low
Reynolds number where flow is stationary, or perhaps even prescribe the
flow. Then you know what the forces will be, and you can think about
what the displacement should be.
Best
W.
--
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 visit
https://groups.google.com/d/msgid/dealii/5db0a212-c591-4863-a48d-f1b236e7298a%40colostate.edu.
