Advertisement
Guest User

Untitled

a guest
Jun 16th, 2019
71
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.46 KB | None | 0 0
  1. Clear["Global`*"]
  2. Ac = 1.98 e - 6;
  3. k = 6.7;
  4. P = 0.005;
  5. U = 4;
  6. Cp = 0.5263;
  7. q = 300;
  8. [Sigma] = 5.6703*10^-8;
  9. [Epsilon] = 0.5;
  10. To = 295;
  11. Ti = 295;
  12. Tl = 2273;
  13. h = 0.0025;
  14. l = 0.005;
  15. Ta = 500;
  16.  
  17. eqns = {(Ac*k)*x''[t] - (P*Ac*U*Cp)*x'[t] + P*q ==
  18. 0, (Ac*k)*y''[t] - (P*Ac*U*Cp)*y'[t] - [Sigma]*[Epsilon]*
  19. P*(Ta - Ti)*(Ta^2 - Ti^2) y[t] == 0, x[0] == To, y[l] == Tl,
  20. x[h] == y[h], x'[h] == y'[h]};
  21.  
  22. s = DSolve[eqns, {x[t], y[t]}, t];
  23.  
  24. y[t] /. s;
  25. N[y[t] /. s /. t -> 30]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement