All depends on what do you want to do with your velocity field later. General ways is:
1. Split velocity and density into different caches.
2. Save velocity with volume compress as 16bit float.
3 . Convert your volumes to vdb.
4. Clamp your velocity with vdb clip using density as mask.