PDE solution positivity constraint

I have a 2d PDE and NDSolve returns a solution that takes negative values but I want the function to be positive

     f[x,y]>=0  

everywhere. How do I include this simple constraint?