Se diseñó un método numérico para resolver las ecuaciones de Navier-Stokes tridimensionales, incompresibles y dependientes del tiempo en flujos turbulentos de canal térmico. Su originalidad radica en el uso de varios métodos conocidos para discretizar el problema y su naturaleza paralela. Se ha utilizado la formulación de vorticidad-laplaciano de la velocidad, por lo que la presión ha sido eliminada del sistema. El calor se modela como un escalar pasivo. Cualquier otra cantidad modelada como escalar pasivo puede ser estudiada muy fácilmente, incluyendo varias al mismo tiempo. Estos métodos se han utilizado con éxito para extensas simulaciones numéricas directas de flujo térmico pasivo para varias condiciones de contorno.