Dirichlet Boundary Condition Bug

udcm-su / NTMpy

Simulation of multi temperature and multi layer heat diffusion

http://udcm.fysik.su.se/code

MIT License

18 stars 11 forks source link

Dirichlet Boundary Condition Bug #10

Closed VaSca92 closed 5 years ago

VaSca92 commented 5 years ago

Dirichlet Boundary condition do not work properly:

Right boundary condition seems to be fixed at zero and also the numerical value of the left boundary condition is not correct when the conductivity is not unitary

A bug fix is needed

luksen99 commented 5 years ago

The problem is fixed. There was an issue with the computation of the flux in the main loop. Now we are treating the right-hand side boundary correctly and distinguish between Neumann and Dirichlet type correctly.