Correct -1 in force calculation procedure (F = -gradU)
Roll back to normal flow
Average vlocity correction for energy conversation