Heat Transfer Today

Educational Software for
Heat and Mass Transfer

Student Projects

Self-similar Boundary Layer Calculations

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:

Temperature arising in self-similar, laminar boundary layer flow on a flat plate due to viscous dissipation. Three representative fluids, Mercury ( Pr << 1), Water ( Pr >1 )and Air ( Pr ~ 1 ) are included for comparison. As expected, the computed wall temperatures (adiabatic wall temperatures) are found to be about the square root of the Prandtl number.

 

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.

Click for pdf writeup

© 2025 Heat Transfer Today – All Rights Reserved.