Solving Blasius Equation with the Shooting Method

This code solves the Blasius equation (third-order ordinary differential equation) for boundary layer flow over a flat plate.
2.7K Downloads
Updated Sat, 03 Nov 2018 16:39:01 +0000

View License

The equation we wish to solve is f''' + (1/2)*f*f'' with f(0) = 0, f'(0) = 0, f'(inf) = 1. This equation arises in the theory of fluid boundary layers, and must be solved numerically. We recast this problem as a system of first-order ODEs: y = [f; f'; f''] = [y(1); y(2); y(3)] so that dy/dEta = y' = [f'; f''; f'''] = [y(2); y(3); -(1/2)*y(1)*y(3)] with y(1)(0) = 0, y(2)(0) = 0, y(2)(inf) = 1. This new system of equations may then be solved numerically using the shooting method. For a description of this numerical method, please follow the given link: https://en.wikipedia.org/wiki/Shooting_method.

Cite As

Mohammad Alkhadra (2024). Solving Blasius Equation with the Shooting Method (https://www.mathworks.com/matlabcentral/fileexchange/69310-solving-blasius-equation-with-the-shooting-method), MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2018b
Compatible with any release
Platform Compatibility
Windows macOS Linux
Tags Add Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
Version Published Release Notes
1.0.0