When using volume procedural (for example, for gridless advection), the resolution of the container becomes equal to 100 voxels along each axis. Therefore, volume quality must be increased:

volume quality =  (size along a larger dimension) / 100

This value is best set on the object.
If the rendering takes too long, you can create a multiplier for volume quality and reduce its value by a couple of times.

In volume procedural, you can use VDB, but you need to do vector split. The results of working with VDB and volume are almost the same.

Also in volume procedural, you can increase octree divisions. This eliminates the need to calculate empty areas. A value of 64 is sufficient; large values will take a long time to prepare before rendering.