The Square 2D case involves thermal diffusion in a horizontal flow on a unit square domain. The numerical heat flux at the bottom boundary (y=0) is checked with the analytical heat flux to a specified tolerance. All boundaries are set to a thermal condition which is a function of x and y, and the initial condition is set to 1. |