Laminar flow over a flat plate is a problem that has been studied extensively, both analytically and experimentally. A solution to a simplified form of this problem can be obtained when several assumptions are made. These assumptions were first presented by Blasius in 1908. Comparisons of experiments with this solution have proved to be accurate.
Consider a flow at speed U past a thin plate of length L. The time a fluid particle spends near the body is approximately L/U, while the time required for viscous effects to spread across the streamlines is of order √vL/U³. If the diffusion time is much shorter than the residence time, the viscous region will be thin. This implies that the Reynolds number should be large enough. This condition is known to hold true in the region 1000< Re < 106 outside of which the flow becomes turbulent.
The Blasius solution assumes that there is an inner boundary layer region below an outer region at which the flow remains the same as free stream conditions. A numerical simulation is set up to simulate the same conditions.
A Mach number of 0.1 and a Reynolds number of 1.0 x 105 were chosen for the simulation. The total computational domain extended two length units in the streamwise direction with a leading edge at x = 0. The extent of the domain in the normal direction is four length units. The computational grid employed 200 x 50 x 3 grid points in the i, j, and k directions, which also align with the x, y, and z axes respectively. The simulation type is simple CFD with compact scheme used for spatial differencing and the Beam-Warming scheme used for time differencing. The problem is solved as though two-dimensional in i and j directions. (Note that, although no spatial differencing is done in the k-direction, three grid points are still required in this direction as this is the minimum number of grid points allowed in any direction.)
The computational mesh is set up with graded elements close to the wall in the normal direction. In addition, the mesh is graded slightly in the streamwise direction, favoring more grid points close to the inlet. The computational mesh is shown in Figure P2.1(a). The initial conditions are as follows:
Change the directory to the subfolder with the set up file. Start the simulation by
mpirun –np 1 mpiaeroflo.exe < flateplate.afl
White, F. M., “Viscous Fluid Flow”, 2nd Ed., McGraw-Hill, 1991. pp 233-237.
© 2025 All Rights Reserved.
