- #1
carpediem85
- 2
- 0
hi I am trying to write a Fortran 77 code to solve the Blasius equation numerically:
blasius equation: F''' + (1/2)*F*F'' = 0
boundary conditions: f(0)=0.0, F'(0)=0, and limit of F'(eta) as eta approaches infinity is 1.0.
This differential equation represents the velocity profile for an incompressible and laminar flow over a flat plate.
So far I know this about how to do it:
1. I need to take this 3rd order ode and build asystem of of 1st order odes. I do not how many, probably 3 or 2, I am not sure at this points.
2. i need to use the shooting method to approximate the last boudnary condition and find f''(0) that I can use to continue to implement Runge-Kutta march.
If anyone please knows or has already solved this problem , i would really appreciate that you write to me or post here your advice. i am confused and cannt get tarted..once i know the differential equations and the initial conditions properly, i can write the runge kutta algorithm in no time.
thanks.
blasius equation: F''' + (1/2)*F*F'' = 0
boundary conditions: f(0)=0.0, F'(0)=0, and limit of F'(eta) as eta approaches infinity is 1.0.
This differential equation represents the velocity profile for an incompressible and laminar flow over a flat plate.
So far I know this about how to do it:
1. I need to take this 3rd order ode and build asystem of of 1st order odes. I do not how many, probably 3 or 2, I am not sure at this points.
2. i need to use the shooting method to approximate the last boudnary condition and find f''(0) that I can use to continue to implement Runge-Kutta march.
If anyone please knows or has already solved this problem , i would really appreciate that you write to me or post here your advice. i am confused and cannt get tarted..once i know the differential equations and the initial conditions properly, i can write the runge kutta algorithm in no time.
thanks.