I'm having a problem where i want to solve a 2D parabolic heat equation on a rectangular grid with a hole in the middle (circle with curved boundaries).
I have a scheme that applies to the solution of every point on the grid,(see attached file scheme.jpg) but I'm unsure how to detect the distance to the neighboring problematic grid points (For example, when a grid point's neightbor is inside the circle, or on the circle boundary, the distance to the circle boundary would be less than the usual distance of 1).
Please look at my attached picture to see an illustration of the problem. Problematic points are denoted in "red". The "coordinates" for each point is the matrix index for each point of alpha and beta.
So in MATLAB I want to make a code that somehow detects the value of alpha and beta for these problematic grid points next to the circle boundary. Can anybody give me a hint of where to start? I'm unsure how to approach and code this.
Best Jonathan Gundorph
Best Answer