What is the more "fenics-y" approach to this FSI problem?
3 months ago by
1) Extract the BoundaryMesh of the Mesh, and then define a subdomain on the BoundaryMesh using a meshfunction (defined on the bmesh) to mark the desired Mesh facets. Then the FSI problem is solved using a partitioned approach: The solid equations are solved on a function space defined on the subdomain of the bmesh, and then mapped back to the function space of the 2D fluid mesh. I can't find a built-in mapping for this so I adapted the approach from the answer to this post (https://www.allanswered.com/post/wrogg/#gnrlj ). However, it seems like this mapping breaks down for finer meshes, since I am getting non-zero values on my function defined on the mesh (instead of only on the boundary. Perhaps there is a better way to map from bmesh subdomains back to the parent mesh?
2) Define the subdomain on the boundary of the 2D mesh, using a facet function, marking all other boundaries and the interior. Then integrate along the boundary and set the other boundaries and interior dofs == 0 with a DirchletBC. However, I keep getting the "integral of type cell, cannot contain a reference normal" error.
Thanks for your input!
Community: FEniCS Project
Please login to add an answer/comment or follow this question.