Abstract: | Finite-difference-like discretizations are developed for the time-dependent Navier–Stokes equations and the thermal energy equation on Delaunay triangulations of the flow domain. The flow is assumed to be thermally expandable; that is, the density varies only with temperature. A characteristic- like (CL) method is used to discretize the temporal and convective terms. The dual variable method reduces the size of the discrete system by a factor of five. |