This circle boundary layer case involves a heat source being applied to a vortex flow around a circular annulus of radius r = [ 1 , 2 ]. The numerical heat flux at the inner boundary is checked with the analytical heat flux at the inner boundary to a specified tolerance. All boundaries are set to a thermal boundary layer condition which is a function of r and θ. The applied heat source is also a function of r and θ, and the initial condition is set to 1. | 
