Replies: 1 comment 1 reply
-
This is difficult. It might be easier to use soft constraint. |
Beta Was this translation helpful? Give feedback.
1 reply
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
-
Hello, everyone!
I am trying to perform a 2D amperometric electrochemical model by utilizing hard-boundary conditions
I have defined only one chemical species and its concentration has neumann boundary conditions on left, right and upper boundaries. On the other hand, the lower boundary is partitioned. While r<0.1, a dirichlet boundary condition with the concentration value of 0.2529 is applied upon the electrode boundary, and neumann boundary condition is valid for the rest of the lower boundary. Moreover, I have an initial condition as C=1.
I have defined the hard boundary conditions as below:
When I executed the code, all collocation and boundary points are always calculated as unity; thus, 0.2529 concentration value would never be predicted on the electrode boundary
I have schematically summarized the boundary/initial conditions with the picture below
Can you help me to solve this problem?
Thank you in advance!
Beta Was this translation helpful? Give feedback.
All reactions