Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Facetset writer #1132

Open
termi-official opened this issue Jan 17, 2025 · 2 comments
Open

Facetset writer #1132

termi-official opened this issue Jan 17, 2025 · 2 comments
Labels
feature good first issue Something to work on as a new contributor!

Comments

@termi-official
Copy link
Member

Since I have been asked for it and for debug purposes, it will be helpful to have a separate Ferrite.write_facetsets function, analogue to Ferrite.write_cellsets.

It should be noted that for quadratic geometric elements we are missing functions to extract all nodes on facets, as Ferrite.faces is reserved for the extraction of the vertices of the face only (describing the uniquely identifiable part of a face).

@termi-official termi-official added feature good first issue Something to work on as a new contributor! labels Jan 17, 2025
@KnutAM
Copy link
Member

KnutAM commented Jan 19, 2025

Just in case someone reads this before its been resolved, a current workaround is (given the dofhandler, dh with a field :u)

ch = ConstraintHandler(dh)
add!(ch, Dirichlet(:u, facetset, Returns(1.0), [1]))
close!(ch)
VTKGridFile(filename, dh) do vtk
    Ferrite.write_constraints(vtk, ch)
end

to export facetset.

@KnutAM
Copy link
Member

KnutAM commented Jan 19, 2025

It should be noted that for quadratic geometric elements we are missing functions to extract all nodes on facets, as Ferrite.faces is reserved for the extraction of the vertices of the face only (describing the uniquely identifiable part of a face).

I believe we can use

local_node_nrs = Ferrite.facetdof_indices(geometric_interpolation(CellType))[facet_nr]

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
feature good first issue Something to work on as a new contributor!
Projects
None yet
Development

No branches or pull requests

2 participants