The expression "-y^2+1.0" had calculated as "(-y)^2+1.0". The "1-y^2" and the "-(y^2)+1.0" are OK. I have attached an IncNavierStokesSolver example. [flow.xml](/uploads/98d0ddd3a44fdad74342aa6c8d2e6906/flow.xml)