CFD Online Discussion Forums

CFD Online Discussion Forums (
-   Main CFD Forum (
-   -   Pohlhausen energy equation MATLAB help (

abe_cooldude May 3, 2010 13:16

Pohlhausen energy equation MATLAB help
Hi all,

I am completely new to this forum, and am impressed with the amount of knowledge posted here!

I am taking few heat transfer and fluid courses this term, and one of them involves deriving solution for the Blasius equation and the Pohlhausen equation in MATLAB. I am not very familiar with MATLAB and don't know how to go about doing so. Both solutions are going to be in tabular form.

Blasius equation: 2f'''+ff''=0 with B.C. f'(0)=df/d(eta) @ eta=0 = 0, f(0)=0, and f'(inf)=df/d(eta) @ eta=inf = 1.

Pohlhausen equation: Tstar''+0.5fTstar'Pr=0 with B.C. Tstar(0)=0 and Tstar(inf)=1

I am more concerned with the second equation. Can anyone guide me in how to implement this in MATLAB?

I really appreciate all your help!


DoHander May 3, 2010 13:38


if you know how to implement Blasius, just apply the same procedure for Polhausen, use a Runge-Kutta method combined with shooting method to tune Tstar'(0)=T_tuned, where T_tuned must satisfy Tstar(inf)=1.


abe_cooldude May 3, 2010 16:58

That's the thing though. I don't know how to implement the Blasius equation in matlab. I have been reading about it ODE45 function in MATLAB and information on runge-kutta and shooting method, but still lost. :confused:

All times are GMT -4. The time now is 10:58.