As I plan to allow more control than those physically available for
standards simulators here I show a small video with two near real smokes
The first candle like, use an interpolated velocity field of only 32x32x32
to calculate the density in a 64x64x64 field (faster to calculate than
full simulation) and the second more turbulent smoke use an iterated
vorticity confinement similar therm to boost turbulence.
edit from the blog feeder: the pseudo-Navier-Stokes solver is built inside blender