I am trying to solve the Heat Equation in 2D for a circular domain and I used the example attached, however, for some reason I do not get any answer from it, and in principle, it seems that I am following the same steps as in the original document from wolfram tutorials. Any help will be much appreciated. I am using version 11.1.1
thank you so much, it worked now, however raised then another question, let say I use Plot3D[Evaluate[u[t,x,y] /. bvpdisk], …] it should give a similar answer also, and for some reason it does not , could you please help me to understand the reason?
NDSolve would return a rule which could be used like this. But in this case you are using NDSolveValue, which returns the solution as an Interpolating Function. So this is for direct use. This has not provided a rule for u, or even a value for u. All it has done is assign a value to bvpdisk, which is the interpolating function.
thank you once again, got it.
I am trying to solve the heat equation inside and outside a region, boundary conditions are based on exchange of heat across the boundary separating both regions. I see NDSolveValue accepts BVP for a single function as well as a description of the region for that specific function, how could I implement the solution with NDSolveValue inside a circle and outside and to match both solutions at the boundary with Mathematica?
