ode45 for Thwaites Method
조회 수: 6 (최근 30일)
이전 댓글 표시
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.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/273407/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/273408/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/273409/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/273410/image.png)
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
댓글 수: 0
답변 (1개)
Pravin Jagtap
2020년 2월 27일
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.
댓글 수: 0
참고 항목
카테고리
Help Center 및 File Exchange에서 Ordinary Differential Equations에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!