In this project we implement a 4th order Runge-Kutta integration scheme to solve the Blasius equation for laminar flow over a flat plate. One of the three boundary conditions is out at the freestream (η = ∞). Therefore, we implement a shooting method based on the half-interval method. Once the solution for the flow is in hand, we use those results in solving the temperature distribution for two cases. The first involves a fixed wall temperature and no viscous dissipation. The second includes viscous dissipation included and the wall kept adiabatic. Both cases result in a tridiagonal system for the temperature distribution. For the latter scenario, we present the final temperature distributions below:
Note that a much more versatile and highly interactive forced convection boundary layer virtual laboratory is available as part of the Heat Transfer Today collection as a Visual Studio executable.