Physics
Lava lamp
Watch warm wax rise and cool wax sink.

Try this
Follow the wax as heat changes its density. Lower heater power to reverse the buoyancy. CUDA, OpenCL and Vulkan execute the same grid equations; CPU is the fallback.
The model and its limits
Three-dimensional thermal flow with buoyancy, viscosity, pressure projection, transported wax fraction and approximate surface tension. Wax volume is corrected for numerical transport drift. This is a bounded teaching model, not a validated multiphase CFD solver; phase change and chemistry are omitted.
Equations used in the module
Equations F_b = (ρ_l - ρ_w) V g ρ_w = ρ_0 / (1 + α (T - T_0))


