Info
이 질문은 마감되었습니다. 편집하거나 답변을 올리려면 질문을 다시 여십시오.
How to solve a set of coupled ODEs (4th Order)?
조회 수: 2 (최근 30일)
이전 댓글 표시
Could somebody help me solve this set of four coupled ODEs?
tpcl = @(xi,y) [y(2); ((1 + (y(2)*y(2)))^(1.5))*(y(3) - (A/(y(1)^3)))/sigma; -3*mu_l*y(4)/(rho_l*h_lv*(y(1)^3)); (Twall - Tsat - (Tsat*y(3)/(rho_l*h_lv)))/((y(1)/k_l) + R_int)];
xi_step = 0:1e-8:xi_end;
[XI,Y] = ode15s(tpcl,xi_step,[delta_ad deltaprime_ad delta_p_ad Q_ad]);
I tried both ode45 and ode15s, but the results are unrealistic/ The constants used here are as follows:
A = 1e-20;
R_FC = 24.597;
c = 1;
sigma = 0.01;
rho_l = 1593.84;
k_l = 0.0544;
mu_l = 0.0004683;
rho_v = 11.61;
h_lv = 88000;
Twall = 334;
Tsat = 329;
R_int = ((2 - c)/(2*c))*(Tsat^(1.5))*((2*pi*R_FC)^(0.5))/(rho_v*(h_lv^2));
delta_ad = (A/(rho_l*h_lv*((Twall/Tsat) - 1)))^(1/3);
deltaprime_ad = 0;
delta_p_ad = A/(delta_ad^3);
Q_ad = 0;
댓글 수: 0
답변 (0개)
이 질문은 마감되었습니다.
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!