Taylor-Maccoll equations for axis-symmetric conical flow using the ode45 solver

25 views (last 30 days)
SRIKESHAV VENKATARAMAN NIKHILESH
Answered: Rohan Kokate on 18 Dec 2020
I wrote the differential equations a system of equations. I have reproduced the function as follows - function dvrdtheta= TMequation(theta,vr) gamma=1.4; dvrdtheta=zeros(2,1); dvrdtheta(1) = vr(2); dvrdtheta(2) =(((gamma-1)/2)*(1-vr(1).^2-vr(2).^2)*(2*vr(1)+cotd(theta).*vr(2))-vr(1).*vr(2).^2)./(vr(2).^2-((gamma-1)/2)*(1-vr(1).^2-vr(2).^2)); I just don't know where I got this wrong. Please help me fellas!
  1 Comment
Rajan Bhandari
Rajan Bhandari on 4 May 2018
I also have been working to solve this equation but I got the problem with integration tolerance. What is your problem actually?

Sign in to comment.

Answers (1)

Rohan Kokate
Rohan Kokate on 18 Dec 2020
Attaching the link to the code which solves the Taylor Maccoll equation numerically. Hope this helps!
https://www.mathworks.com/matlabcentral/fileexchange/84490-taylor-maccoll-solution-for-supersonic-flow-over-a-cone

Community Treasure Hunt

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

Start Hunting!

Translated by