Visualizing3Dcontoursoftotalvelocity:
In[]:=
Show[bmesh["Wireframe"],ContourPlot3D[Evaluate[Sqrt[xVel[3.65,x,y,z]^2+yVel[3.65,x,y,z]^2+zVel[3.65,x,y,z]^2]],{x,y,z}∈mesh],MeshNone]
Out[]=