A Numerical Method for Incompressible Flow with Heat Transfer

A numerical method for the convective heat transfer problem is developed for low speed flow at mild temperatures. A simplified energy equation is added to the incompressible Navier-Stokes formulation by using Boussinesq approximation to account for the buoyancy force. A pseudocompressibility method is used to solve the resulting set of equations for steady-state solutions in conjunction with an approximate factorization scheme. A Neumann-type pressure boundary condition is devised to account for the interaction between pressure and temperature terms, especially near a heated or cooled solid boundary. It is shown that the present method is capable of predicting the temperature field in an incompressible flow