Solving second order inhomogenous PDE by separation of variables requires homogenization of the boundary conditions. Let's say we are looking at 1D heat equation. From intuition, if we have fixed temperature on both sides(inhomogeneous Dirichlet-Dirichlet boundary conditions), there is no heat coming in or out of the 1D bar, meaning as time goes to infinity, the bar will reach an equilibrium state where the temperature would no longer depend on time meaning. So when times go to infinity the solution would be a function u(x) (so-called homogenization function), meaning the heat equation is: $$d^2u/dx^2=0$$ with the Dirichlet boundary conditions. The solution to this is $$u=c1*x+c2$$ and by applying the the conditions we can find c1 and c2. Same intuition works to generate the homogenization function if we have boundary condition Dirichlet-Neumann/Neumann-Dirichlet conditions ($u(0,t)=f(t), du(l,t)/dx=g(t) / du(o,t)/dx=f(t), u(l,t)=g(t)$). Unfortunately, when we have Neumann-Neumann conditions ($du(0,t)/dx=f(t),du(l,t)/dx=g(t)$), this intuition is not applicable and I cannot get the homogenization function in it's proper form which is $$u(x,t)= (f(t)+g(t))*x^2/l - f(t)x$$. My guess is that when there is prescribed flux on both sides there cannot be equilibrium temperature as time goes to infinity. I can see that if $d^2u/dx^2=constant$ I will reach the answer but what is the physical intuition to this statement?
P.S. Sorry if some of the terminologies are not correct in English. Also sorry for the unedited equations, I do not know how to write them out.
Edit: Basically I am trying to solve the following problem:
Having the heat equation $$\partial u(x,t)/\partial t=-\epsilon \partial^2u(x,t)/\partial x^2$$
$$ \partial u(0,t)/\partial x=f(t)$$ $$ \partial u(l,t)/\partial x=g(t)$$ $$ u(x,0)=v(x) $$
So to apply separation of variables we need to make the BC homogenous. So we make the assumption that the solution is in the form:
$$u(x,t)=u_0(x,t)+\tilde u(x,t) $$ which basically converts the equation to:
$$\partial u_0(x,t)/\partial t-\epsilon\partial^2u_0(x,t)/\partial x^2=-\partial \tilde u(x,t)/\partial t+\epsilon\partial^2\tilde u(x,t)/\partial x^2=\tilde s(x,t)$$
and the BC transform to: $$\partial u_0(0,t)/\partial x=0$$ $$ \partial u_0(l,t)/\partial x=0$$ $$ u_0(x,0)=v(x)-\tilde u(x,0)$$ which is a homogenous problem with a source and modified initial conditions. The only thing we need is to find the explicit form of $\tilde u$ which I know is in the form of $\tilde u(x,t)= (f(t)+g(t))*x^2/l - f(t)x$ but I do not know how to get to it.