Consider $$u_t - \Delta u + u = 0$$ $$\frac{\partial u}{\partial \nu} = 0$$ $$u(T) = u(0)$$ on a domain $\Omega$ (the BC is obviously on $\partial\Omega$. If $u$ solves this PDE, clearly as does $\lambda u $ for any constant $\lambda.$ So uniqueness is not there.
Suppose we have two solutions $u$ and $v$. Their difference $d = u-v$ satisfies $$d_t - \Delta d + d = 0$$ $$\frac{\partial d}{\partial \nu} = 0$$ $$d(T) = d(0)$$
Multiplying the equation by $d$ and integrating by parts we get $$\frac{1}{2}\frac{d}{dt}\int_{\Omega} d^2(t) + \int_{\Omega} |\nabla d(t)|^2 + \int_{\Omega} d^2(t) = 0$$ Now integrating by time $$\frac{1}{2}|d(T)|_{L^2} - \frac{1}{2}|d(0)|_{L^2} + \int_{0}^T |\nabla d(t)|_{L^2}^2 + \int_0^2 |d(t)|_{L^2}^2 = 0$$ but the first two terms cancel each other out, and we get $d=0$ in $H^1$. So this shows that there is a unique solution.
What am I doing wrong???!?!
Edit: maybe $0$ is the only solution to this problem? How to prove this if so? What happens in nonhomogenous case?
A slightly different approach: writing $y(x,t)=e^tu(x,t)$, the system becomes equivalent to the following homogeneous heat equation: $\displaystyle \frac{dy}{dt}-\Delta y=0$, $\displaystyle\frac{\partial y}{\partial\nu}=0,\, y(x,T)=e^Ty(x,0)$. Now multiplying by $y$ and integrating by parts, we get $\displaystyle \frac{d}{dt}\int_{\Omega}\frac{y^2}{2}=-\int_{\Omega}|\nabla y|^2$, so that $\displaystyle\int_{\Omega}\frac{y^2}{2}$ is decreasing in time. Hence $\displaystyle \int_{\Omega}y^2(x,T)=e^{2T}\int_{\Omega}y^2(x,0)\ge e^{2T}\int_{\Omega}y^2(x,T),$ and this implies $y(x,0)=y(x,T)=0$, and also since $\displaystyle\int_{\Omega}\frac{y^2}{2}$ is decreasing in time but always positive and is zero initially, $y(x,t)=0$.