Suppose I have the parabolic system $$u_t=\nabla\cdot(k(x)\nabla u)+f,\quad (x,t)\in\Omega\times I$$ with Dirichlet boundary conditions $$u=g, \quad x\in\partial\Omega$$ and initial condition $$u(x,t)= h,\quad t=0.$$
Often times in engineering, we are more interested in the asymptotic (steady state) behavior of this PDE rather than the transient behavior. So, we sometimes neglect the time derivative term and solve the elliptic system $$-\nabla\cdot(k(x)\nabla u)=f,\quad (x,t)\in\Omega\times I$$ instead. The assumption is that over infinite time, $$\lim_{t\rightarrow \infty}u_{\mathrm{parabolic}}(x,t)= u_{\mathrm{elliptic}}(x,t)$$.
I've observed that when $f\equiv 0$, this limit is true, but I'm not sure if this is the case for arbitrary $f$ or if there are other necessary conditions to guarantee this limit to be true. Do the boundary conditions have to converge asymptotically to a constant value in order for the parabolic solution to converge to the elliptic solution?
Though I'm my question is formulated in the continuous case, I'm also curious if the same conditions are true for the discrete case. That is, assuming I use a stable and consistent finite difference scheme to approximate $u_{\mathrm{parabolic}}$ and $u_{\mathrm{elliptic}}$, should I expect $$\lim_{t\rightarrow \infty}u_{\mathrm{parabolic}}^\mathrm{fdm}(x,t)= u_{\mathrm{elliptic}}^\mathrm{fdm}(x,t)$$ if $u_{\mathrm{elliptic}}^\mathrm{fdm}$ and $u_{\mathrm{parabolic}}^\mathrm{fdm}$ are discretized on the same spatial grid and $\Delta t\rightarrow\infty$?