Great, this is already showing how much better lighting looks without caching.
Phase function - the downside of a realistic phase function is that scattering takes much longer to resolve. I've observed this with other path tracers (Corona renderer), volumes with values around 0.5-0.6 (forward scattering) will resolve much quicker than realistic values around 0.8-0.9. Visually, the higher values look much more realistic, especially when the sun is behind the camera. The tricky part is that the phase function changes across the cloud and to get realistic values/balance.
Since you mention backlit clouds looking wrong - would you still be willing to post it so we can have a peek even if it's not right yet?
Anyway, again great to see progress here and thanks for keeping us updated.