Maybe there are something wrong on the boundary condition on phi

#77 · open · 1 comments

View on GitHub ↗

candy29333

Hello, I’ve been using this software to write programs recently. On a whim, I printed out the values of the electric potential on the upper and lower surfaces. To my surprise, I found that they differ significantly from the prescribed fixed values (0 and 5). Could you please advise whether there might be an issue somewhere in the code? <img width="271" height="378" alt="Image" src="https://github.com/user-attachments/assets/fe836282-2fab-4fe9-81f1-f007c97379e8" />

Comments

ajnonaka

Are those values from the cell-centers adjacent to the boundary? If so this could make sense as the Dirichlet values are imposed exactly on the boundary (the face). So the cell-centered values will move off of the prescribed values by a little bit.