In this physically-based tornado simulation, the tornado-scale approach techniques are applied tosimulate the tornado formation environment. The three-dimensional Navier-Stokes equations for incompressible viscousfluid flows are used to model the tornado dynamics. The boundary conditions applied in this simulation lead torotating and uplifting flow movement as found in real tornadoes and tornado research literatures.Moreover, a particle system is incorporatedwith the model equation solutions to model the irregular tornado shapes.Also, together with appropriate boundary conditions, varied particle controlschemes produce tornadoes with different shapes. Furthermore,a modified metaball scheme is used to smooth the density distribution. Texturemapping, antialising, animation andvolume rendering are applied to produce realistic visual results. The rendering algorithm is implemented inOpenGL.