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

Add function to compute gaps for Mortar #19

Open
ebchin opened this issue Aug 18, 2023 · 1 comment
Open

Add function to compute gaps for Mortar #19

ebchin opened this issue Aug 18, 2023 · 1 comment

Comments

@ebchin
Copy link
Member

ebchin commented Aug 18, 2023

It would be convenient to have a way to compute gaps without updating forces and Jacobian contributions in the mortar method. See this code in Serac. There, when we are using a penalty enforcement of the constraint we have to call update() twice for the code to work properly: once to update the gaps for computing the updated pressures, then again to compute the forces and Jacobian with the updated pressures.

@srwopschall
Copy link
Collaborator

Yeah right now you'd have to modify the implicit evaluation option to compute gaps only and then call update(). I don't think the gap, residual, jacobian calculations in the mortar physics routine are guarded correctly, but that is an easy fix. We could easily write an API function that computes the gaps for a given coupling scheme and currently registered mesh data, but if you wanted to update the geometry of the overlaps you may as well call update. If that feels awkward, maybe we have an API routine that actually wraps update(), but to the user feels like they are getting what they are asking for? Let's discuss more.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants