evstigneevnm

Stokes example Rhie-Chow is correct?

5

Using stokesCavity.py added verification of `div(V)` as: ```python velocity_d = FaceVariable(mesh=mesh, rank=1) velocity_d[0] = xVelocity.arithmeticFaceValue; velocity_d[1] = yVelocity.arithmeticFaceValue; velocity_d[..., mesh.exteriorFaces.value] = 0. velocity_d[0, mesh.facesTop.value] = U div_value = max(abs(velocity_d.divergence.value)) ```...