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

Non-uniform Dirichlet BCs appear to not allow time dependence #406

Open
wband opened this issue Jan 16, 2025 · 4 comments
Open

Non-uniform Dirichlet BCs appear to not allow time dependence #406

wband opened this issue Jan 16, 2025 · 4 comments

Comments

@wband
Copy link
Contributor

wband commented Jan 16, 2025

I've been developing a phase-field fracture app that uses a time-dependent non-uniform Dirichlet boundary condition to steadily propagate a crack across the domain. Everything seems to be working except for the time dependent boundary condition, which stays identically the same at all iterations.

Here is the app in question: https://github.com/wband/phaseField/tree/fracture-app/applications/fracture

@landinjm
Copy link
Contributor

Hmm this might be because it has a TIME_INDEPENDENT field. I can double check constraint distribution.

@david-montiel-t
Copy link
Contributor

david-montiel-t commented Jan 20, 2025

The time-dependent nonuniform Dirichlet BC issue which seems to be the same as what you are describing was originally solved in commit 2ade5c2c62ea6b7af0fbab0d7dba7b690d88ff1e.

@david-montiel-t
Copy link
Contributor

david-montiel-t commented Jan 20, 2025

I'm pretty sure I tested this fix with problem 6d of the training set which has a TIME_INDEPENDENT field.

@wband
Copy link
Contributor Author

wband commented Jan 20, 2025

In the current version, it looks like the code block starting around line 181 of init.cc will overwrite all of the BC flags to match those of the last field.

Other issue is that the 'DIRICHLET' and 'NON_UNIFORM_DIRICHLET' flags appear to be mutually exclusive, which meant that the conditional triggering application of Dirichlet BCs needed to be 'DIRICHLET or NON_UNIFORM_DIRICHLET' instead of just 'DIRICHLET'.

I'm making a pull request to fix this.

@landinjm landinjm mentioned this issue Feb 4, 2025
17 tasks
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

3 participants