MATLAB and Simulink resources for Arduino, LEGO, and Raspberry Pi

Learn moreOpportunities for recent engineering grads.

Apply Today
Asked by Shashank on 28 Jun 2013

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;

*No products are associated with this question.*

## 0 Comments