Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Velocity and Temp. profiles for stagnation point (m=1)

  1. Apr 10, 2007 #1

    Can anyone assist me in the findings of the velocity and temp. profiles for stagnation point flow (m=1) for Varoius values of Pr (.6,1,1.5,5,10). At this moment I am using Matlab's ode45 to solve this problem along with applying the Blasius Eqn (f''' +1/2f*f''=0) and Energy Eqn. (2*theta''+Pr*f*theta'=0). I have been reading on the use of ode45 and I think I know what I need to do, but I am not sure how to incorporate the above equations along with the B.C. that are f'=0, f=0, theta=o and phi=0.

    Can someone please suggest any info., that would be great.

    so far this is what I have in my M-file - this is the main part of the code:

    tintval=[ti tf];
    bcinit=[0.0 0.0 0.33206];

    [t,y]= ode45(@stagnation,tintval,bcinit);

    %t=the scalar time
    %y=the column vector
    %ode45 is the solver
    %@rigid is the function handle calling function
    %[0.0 15.0] is the time to be evaluated from ti to tf
    %[0.0....etc] is the initial conditions

    this next part is the function that is being called by ode45:

    function df=stagnation(t,f)

    dfdt = [

    Of course it is not much, but that is why I need some help. Thank you
  2. jcsd
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook

Can you offer guidance or do you also need help?
Draft saved Draft deleted