ode45 for Thwaites Method

19 visualizaciones (últimos 30 días)
Zachary Lawless
Zachary Lawless el 24 de Feb. de 2020
Respondida: Pravin Jagtap el 27 de Feb. de 2020
I am trying to use the ode45 solver to solve boundar layer momentum thickness using Thwaites Method for a flat plate
The equation I am trying to put into ode45 is with the initial given condition for theta being which is from the analytical blasius solution. and are known for each point along the flat plate and I want to input them into ode45. I am however at a loss at how to translate this into Matlab.
My attempt:
[xode theta] = ode45(@(x,theta) thwaites(x,Ve,dVedx,v,theta) ,[0 1],Theta0)
function dthetadx = thwaites(x,Ve,dVedx,v,theta)
sym theta
dthetadx = .225 .* v./(Ve .* theta) - (3 .* theta)./Ve .* dVed
end

Respuestas (1)

Pravin Jagtap
Pravin Jagtap el 27 de Feb. de 2020
Hello Zachary,
Refer to the following documentation link of ‘ode45’ function. This will give your better insights into how to use ‘ode45’ function and its suitability for your problem along with examples:
Hope this will help you.

Categorías

Más información sobre Ordinary Differential Equations en Help Center y File Exchange.

Productos


Versión

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by