Question: How do I solve diffusion non-gomogenous PDE in bipolar coordinates in Maple?

Hi! I need a help for solve PDE problem (diffusion problem) in bipolar system coordinates. Actually, to begin with, I wanted to try to solve the homogeneous diffusion equation in a bipolar coordinate system with given boundary conditions. The process is not stationary, for simplicity, the diffusion coefficient is a constant (or linearly depends on η). The source function is also specified on the right side. The problem is that I cannot construct a correct exact solution, and the PDESolve function gives a very strange solution. The question is whether it is possible to solve a homogeneous and then a non-homogeneous solution automatically, or will it be necessary to manually specify parts of the general solution (homogeneous and particular solution)?

A document with a brief description and code for solving the equation is provided in the file. I ask for help with solving this problem, any valuable comments will be useful!

PDE.mw

====UPDATE===
I try to set Neumann conditions for solve this problem, but get error:
Error, (in PDEtools:-Library:-NormalizeBoundaryConditions) unable to isolate the functions {u(0,eta), u(xi,1.1), u(2*Pi,eta)} in the given boundary conditions {u(0,eta) = u(2*Pi,eta), u(xi,1.1) = 0}

PDE_upd.mw

Please Wait...