The decay constant should be modified to include $\nabla \cdot \mathbf{u}$. Woodcock tracking or similar method should be use for DNPs transport.
