There is plenty of differential equations such as drift-diffusion ones that have a non-negative solution (positive initial and boundary conditions, positive forcing term). However the pod-based model reduction cannot assure that the positiveness is preserved. I would really like to know whether there is any derived model reduction approach that can guarantee the positiveness.