Giulianot

25 Reputation

2 Badges

9 years, 25 days

MaplePrimes Activity


These are questions asked by Giulianot


 

dsolve({Q(0) = 0, Q(t) = (1.375*4190)*(80-T__1(t)), Q(t) = (1.375*4190)*(T__2(t)-38.2), diff(Q(t), t) = (0.1375e-1*(T__1(t)-T__1s(t)))*((T__1(t)+T__1s(t))*(1/2)), diff(Q(t), t) = (0.1375e-1*(T__2s(t)-T__2(t)))*((T__2s(t)+T__2(t))*(1/2)), diff(Q(t), t) = (240*0.1375e-1)*(T__1s(t)-T__2s(t))/(0.1e-2)}, numeric)

Error, (in dsolve/numeric/DAE/initial) missing initial conditions for the following: {T__1s}

 

``

i got 3 diff ecuations with two algebraic ones. a system of DAEs. there is only a derivative included on systems, for which it's necesary only one initial condition for solving the system, which is Q(0)=0. why maple wants to know initial conditions for T_1s. it's not supposed to calculate it itself?
 

Download ecuation_2.mwecuation_2.mw


 

diff(Q(t), t) = k*A*(T__1s(t)-T__2s(t))/d

diff(Q(t), t) = h__1(t)*A*(T__1(t)-T__1s(t))

diff(Q(t), t) = h__2(t)*A*(T__2s(t)-T__2(t))

Q(t) = m__1*c__p*(T__1i-T__1(t))

Q(t) = m__2*c__p*(T__2(t)-T__2i)

h__1(t) = k(T__1(t), T__1s(t))*(.825+.387*(g*h^3*c__p*beta(T__1(t), T__1s(t))*rho(T__1(t), T__1s(t))^2*(T__1(t)-T__1s(t))/(k(T__1(t), T__1s(t))*mu(T__1(t), T__1s(t))))^(1/6)/(1+(.492*k(T__1(t), T__1s(t))/(c__p*mu(T__1(t), T__1s(t))))^(9/16))^(8/27))^2/h

h__2(t) = k(T__2(t), T__2s(t))*(.825+.387*(g*h^3*c__p*beta(T__2(t), T__2s(t))*rho(T__2(t), T__2s(t))^2*(T__2s(t)-T__2(t))/(k(T__2(t), T__2s(t))*mu(T__2(t), T__2s(t))))^(1/6)/(1+(.492*k(T__2(t), T__2s(t))/(c__p*mu(T__2(t), T__2s(t))))^(9/16))^(8/27))^2/h

 

 

rho(T__1(t), T__1s(t)) = 999.9399+2.1082425*10^(-2)*(T__1(t)+T__1s(t))-1.77436275*10^(-3)*(T__1(t)+T__1s(t))^2+.438696375*10^(-5)*(T__1(t)+T__1s(t))^3-.6189861563*10^(-8)*(T__1(t)+T__1s(t))^4

beta(T__1(t), T__1s(t)) = -(4.216485*10^(-2)-7.097451*10^(-3)*(T__1(t)+T__1s(t))+2.63217825*10^(-5)*(T__1(t)+T__1s(t))^2-4.9518879*10^(-8)*(T__1(t)+T__1s(t))^3)/(999.9399+2.1082425*10^(-2)*(T__1(t)+T__1s(t))-1.77436275*10^(-3)*(T__1(t)+T__1s(t))^2+.438696375*10^(-5)*(T__1(t)+T__1s(t))^3-.6189861563*10^(-8)*(T__1(t)+T__1s(t))^4)
mu(T__1(t), T__1s(t)) = 2.414*10^(247.8/(.5*(T__1(t)+T__1s(t))+133)-5)

k(T__1(t), T__1s(t)) = -9.481411*10^(-6)*(T__1(t)+T__1s(t))^2+2.1356735*10^(-3)*(T__1(t)+T__1s(t))+.5599920949

 

 

rho(T__2(t), T__2s(t)) = 999.9399+2.1082425*10^(-2)*(T__2(t)+T__2s(t))-1.77436275*10^(-3)*(T__2(t)+T__2s(t))^2+.438696375*10^(-5)*(T__2(t)+T__2s(t))^3-.6189861563*10^(-8)*(T__2(t)+T__2s(t))^4

beta(T__2(t), T__2s(t)) = -(4.216485*10^(-2)-7.097451*10^(-3)*(T__2(t)+T__2s(t))+2.63217825*10^(-5)*(T__2(t)+T__2s(t))^2-4.9518879*10^(-8)*(T__2(t)+T__2s(t))^3)/(999.9399+2.1082425*10^(-2)*(T__2(t)+T__2s(t))-1.77436275*10^(-3)*(T__2(t)+T__2s(t))^2+.438696375*10^(-5)*(T__2(t)+T__2s(t))^3-.6189861563*10^(-8)*(T__2(t)+T__2s(t))^4)
mu(T__2(t), T__2s(t)) = 2.414*10^(247.8/(.5*(T__2(t)+T__2s(t))+133)-5)

k(T__2(t), T__2s(t)) = -9.481411*10^(-6)*(T__2(t)+T__2s(t))^2+2.1356735*10^(-3)*(T__2(t)+T__2s(t))+.5599920949

 

"`h__1`(t)=(-9.481411*10^(-6) (`T__1`(t)+`T__1s`(t))^(2)+2.1356735*10^(-3) (`T__1`(t)+`T__1s`(t))+0.5599920949)/(h) (0.825+(0.387 ((g h^(3) `c__p` (-(4.216485*10^(-2)-7.097451*10^(-3) (`T__1`(t)+`T__1s`(t))+2.63217825*10^(-5) (`T__1`(t)+`T__1s`(t))^(2)-4.9518879*10^(-8) (`T__1`(t)+`T__1s`(t))^(3))/(999.9399+2.1082425*10^(-2) (`T__1`(t)+`T__1s`(t))-1.77436275*10^(-3) (`T__1`(t)+`T__1s`(t))^(2)+0.438696375*10^(-5) (`T__1`(t)+`T__1s`(t))^(3)  -0.6189861563*10^(-8) (`T__1`(t)+`T__1s`(t))^(4))) (999.9399+2.1082425*10^(-2) (`T__1`(t)+`T__1s`(t))-1.77436275*10^(-3) (`T__1`(t)+`T__1s`(t))^(2)+0.438696375*10^(-5) (`T__1`(t)+`T__1s`(t))^(3)  -0.6189861563*10^(-8) (`T__1`(t)+`T__1s`(t))^(4))^(2) (`T__1`(t)-`T__1s`(t)))/((-9.481411*10^(-6) (`T__1`(t)+`T__1s`(t))^(2)+2.1356735*10^(-3) (`T__1`(t)+`T__1s`(t))+0.5599920949) 2.414*10^((247.8)/(0.5 (`T__1`(t)+`T__1s`(t))+133)-5)))^((1)/(6)))/((1+((0.492 (-9.481411*10^(-6) (`T__1`(t)+`T__1s`(t))^(2)+2.1356735*10^(-3) (`T__1`(t)+`T__1s`(t))+0.5599920949))/(`c__p` 2.414*10^((247.8)/(0.5 (`T__1`(t)+`T__1s`(t))+133)-5)))^((9)/(16)))^((8)/(27))))^(2)"

"`h__2`(t)=(-9.481411*10^(-6) (`T__2`(t)+`T__2s`(t))^(2)+2.1356735*10^(-3) (`T__2`(t)+`T__2s`(t))+0.5599920949)/(h) (0.825+(0.387 ((g h^(3) `c__p` (-(4.216485*10^(-2)-7.097451*10^(-3) (`T__2`(t)+`T__2s`(t))+2.63217825*10^(-5) (`T__2`(t)+`T__2s`(t))^(2)-4.9518879*10^(-8) (`T__2`(t)+`T__2s`(t))^(3))/(999.9399+2.1082425*10^(-2) (`T__2`(t)+`T__2s`(t))-1.77436275*10^(-3) (`T__2`(t)+`T__2s`(t))^(2)+0.438696375*10^(-5) (`T__2`(t)+`T__2s`(t))^(3)  -0.6189861563*10^(-8) (`T__2`(t)+`T__2s`(t))^(4))) (999.9399+2.1082425*10^(-2) (`T__2`(t)+`T__2s`(t))-1.77436275*10^(-3) (`T__2`(t)+`T__2s`(t))^(2)+0.438696375*10^(-5) (`T__2`(t)+`T__2s`(t))^(3)  -0.6189861563*10^(-8) (`T__2`(t)+`T__2s`(t))^(4))^(2) (`T__2s`(t)-`T__2`(t)))/((-9.481411*10^(-6) (`T__2`(t)+`T__2s`(t))^(2)+2.1356735*10^(-3) (`T__2`(t)+`T__2s`(t))+0.5599920949) 2.414*10^((247.8)/(0.5 (`T__2`(t)+`T__2s`(t))+133)-5)))^((1)/(6)))/((1+((0.492 (-9.481411*10^(-6) (`T__2`(t)+`T__2s`(t))^(2)+2.1356735*10^(-3) (`T__2`(t)+`T__2s`(t))+0.5599920949))/(`c__p` 2.414*10^((247.8)/(0.5 (`T__2`(t)+`T__2s`(t))+133)-5)))^((9)/(16)))^((8)/(27))))^(2)"

diff(Q(t), t) = k*A*(T__1s(t)-T__2s(t))/d, diff(Q(t), t) = A*(T__1(t)-T__1s(t))*(-9.481411*10^(-6)*(T__1(t)+T__1s(t))^2+2.1356735*10^(-3)*(T__1(t)+T__1s(t))+.5599920949)*(.825+.387*(((-g*h^3*c__p*(4.216485*10^(-2)-7.097451*10^(-3)*(T__1(t)+T__1s(t))+2.63217825*10^(-5)*(T__1(t)+T__1s(t))^2-4.9518879*10^(-8)*(T__1(t)+T__1s(t))^3)/(999.9399+2.1082425*10^(-2)*(T__1(t)+T__1s(t))-1.77436275*10^(-3)*(T__1(t)+T__1s(t))^2+.438696375*10^(-5)*(T__1(t)+T__1s(t))^3-.6189861563*10^(-8)*(T__1(t)+T__1s(t))^4))*((999.9399+2.1082425*10^(-2)*(T__1(t)+T__1s(t))-1.77436275*10^(-3)*(T__1(t)+T__1s(t))^2+.438696375*10^(-5)*(T__1(t)+T__1s(t))^3-.6189861563*10^(-8)*(T__1(t)+T__1s(t))^4)^2))*(T__1(t)-T__1s(t))/(2.414*10^(247.8/(.5*(T__1(t)+T__1s(t))+133)-5)*(-9.481411*10^(-6)*(T__1(t)+T__1s(t))^2+2.1356735*10^(-3)*(T__1(t)+T__1s(t))+.5599920949)))^(1/6)/(1+((.492*(-9.481411*10^(-6)*(T__1(t)+T__1s(t))^2+2.1356735*10^(-3)*(T__1(t)+T__1s(t))+.5599920949))/(2.414*c__p*10^(247.8/(.5*(T__1(t)+T__1s(t))+133)-5)))^(9/16))^(8/27))^2/h, diff(Q(t), t) = A*(T__2s(t)-T__2(t))*(-9.481411*10^(-6)*(T__2(t)+T__2s(t))^2+2.1356735*10^(-3)*(T__2(t)+T__2s(t))+.5599920949)*(.825+.387*(((-g*h^3*c__p*(4.216485*10^(-2)-7.097451*10^(-3)*(T__2(t)+T__2s(t))+2.63217825*10^(-5)*(T__2(t)+T__2s(t))^2-4.9518879*10^(-8)*(T__2(t)+T__2s(t))^3)/(999.9399+2.1082425*10^(-2)*(T__2(t)+T__2s(t))-1.77436275*10^(-3)*(T__2(t)+T__2s(t))^2+.438696375*10^(-5)*(T__2(t)+T__2s(t))^3-.6189861563*10^(-8)*(T__2(t)+T__2s(t))^4))*((999.9399+2.1082425*10^(-2)*(T__2(t)+T__2s(t))-1.77436275*10^(-3)*(T__2(t)+T__2s(t))^2+.438696375*10^(-5)*(T__2(t)+T__2s(t))^3-.6189861563*10^(-8)*(T__2(t)+T__2s(t))^4)^2))*(T__2s(t)-T__2(t))/(2.414*10^(247.8/(.5*(T__2(t)+T__2s(t))+133)-5)*(-9.481411*10^(-6)*(T__2(t)+T__2s(t))^2+2.1356735*10^(-3)*(T__2(t)+T__2s(t))+.5599920949)))^(1/6)/(1+((.492*(-9.481411*10^(-6)*(T__2(t)+T__2s(t))^2+2.1356735*10^(-3)*(T__2(t)+T__2s(t))+.5599920949))/(2.414*c__p*10^(247.8/(.5*(T__2(t)+T__2s(t))+133)-5)))^(9/16))^(8/27))^2/h, Q(t) = m__1*c__p*(T__1i-T__1(t)), Q(t) = m__2*c__p*(T__2(t)-T__2i)

diff(Q(t), t) = k*A*(T__1s(t)-T__2s(t))/d, diff(Q(t), t) = A*(T__1(t)-T__1s(t))*(-0.9481411000e-5*(T__1(t)+T__1s(t))^2+0.2135673500e-2*T__1(t)+0.2135673500e-2*T__1s(t)+.5599920949)*(.825+.387*(-.4142502071*g*h^3*c__p*(0.4216485000e-1-0.7097451000e-2*T__1(t)-0.7097451000e-2*T__1s(t)+0.2632178250e-4*(T__1(t)+T__1s(t))^2-0.4951887900e-7*(T__1(t)+T__1s(t))^3)*(999.9399+0.2108242500e-1*T__1(t)+0.2108242500e-1*T__1s(t)-0.1774362750e-2*(T__1(t)+T__1s(t))^2+0.4386963750e-5*(T__1(t)+T__1s(t))^3-0.6189861563e-8*(T__1(t)+T__1s(t))^4)*(T__1(t)-T__1s(t))/(10^(247.8/(.5*T__1(t)+.5*T__1s(t)+133)-5)*(-0.9481411000e-5*(T__1(t)+T__1s(t))^2+0.2135673500e-2*T__1(t)+0.2135673500e-2*T__1s(t)+.5599920949)))^(1/6)/(1+.4087338992*((-0.9481411000e-5*(T__1(t)+T__1s(t))^2+0.2135673500e-2*T__1(t)+0.2135673500e-2*T__1s(t)+.5599920949)/(c__p*10^(247.8/(.5*T__1(t)+.5*T__1s(t)+133)-5)))^(9/16))^(8/27))^2/h, diff(Q(t), t) = A*(T__2s(t)-T__2(t))*(-0.9481411000e-5*(T__2(t)+T__2s(t))^2+0.2135673500e-2*T__2(t)+0.2135673500e-2*T__2s(t)+.5599920949)*(.825+.387*(-.4142502071*g*h^3*c__p*(0.4216485000e-1-0.7097451000e-2*T__2(t)-0.7097451000e-2*T__2s(t)+0.2632178250e-4*(T__2(t)+T__2s(t))^2-0.4951887900e-7*(T__2(t)+T__2s(t))^3)*(999.9399+0.2108242500e-1*T__2(t)+0.2108242500e-1*T__2s(t)-0.1774362750e-2*(T__2(t)+T__2s(t))^2+0.4386963750e-5*(T__2(t)+T__2s(t))^3-0.6189861563e-8*(T__2(t)+T__2s(t))^4)*(T__2s(t)-T__2(t))/(10^(247.8/(.5*T__2(t)+.5*T__2s(t)+133)-5)*(-0.9481411000e-5*(T__2(t)+T__2s(t))^2+0.2135673500e-2*T__2(t)+0.2135673500e-2*T__2s(t)+.5599920949)))^(1/6)/(1+.4087338992*((-0.9481411000e-5*(T__2(t)+T__2s(t))^2+0.2135673500e-2*T__2(t)+0.2135673500e-2*T__2s(t)+.5599920949)/(c__p*10^(247.8/(.5*T__2(t)+.5*T__2s(t)+133)-5)))^(9/16))^(8/27))^2/h, Q(t) = m__1*c__p*(T__1i-T__1(t)), Q(t) = m__2*c__p*(T__2(t)-T__2i)

(1)

"(->)"

``

``

``

i have a system with 5 dif equations and five unknows. i have told to maple to solve it numerically with interactively solve comand (right cilck button). the window open like it normally does and i put values to my parameters, with an initial condition for the system (Q(0)=0). then i press numerically solve and that's all, the program just keep evaluating with no answer. i wait for 15 min, which i think is too much time, and got any answer yet.

hope you can help with this

thanks.. 
 

Download propuesta_transfer.mw


 

T := proc (t) options operator, arrow; 80+(-1)*1.4375*53*(1-exp((-1)*0.13775e-1*2.875*k*h__1*h__2*t/(1.4375^2*4190*(k*(h__1+h__2)+0.1e-2*h__1*h__2))))/2.875 end proc

proc (t) options operator, arrow; 80+(-1)*1.4375*53*(1-exp((-1)*0.13775e-1*2.875*k*h__1*h__2*t/((1.4375^2*4190*(k*(h__1+h__2)+0.1e-2*h__1*h__2)))))/2.875 end proc

(1)

T(204.39)

53.50000000+26.50000000*exp(-8.094482719*k*h__1*h__2/(8658.242188*k*(h__1+h__2)+8.658242188*h__1*h__2))

(2)

T(429.63)

53.50000000+26.50000000*exp(-17.01469059*k*h__1*h__2/(8658.242188*k*(h__1+h__2)+8.658242188*h__1*h__2))

(3)

T(872.5)

53.50000000+26.50000000*exp(-34.55372656*k*h__1*h__2/(8658.242188*k*(h__1+h__2)+8.658242188*h__1*h__2))

(4)

fsolve({53.50000000+26.50000000*exp(-34.55372656*k*h__1*h__2/(8658.242188*k*(h__1+h__2)+8.658242188*h__1*h__2)) = 56.15, 53.50000000+26.50000000*exp(-17.01469059*k*h__1*h__2/(8658.242188*k*(h__1+h__2)+8.658242188*h__1*h__2)) = 64.1, 53.50000000+26.50000000*exp(-8.094482719*k*h__1*h__2/(8658.242188*k*(h__1+h__2)+8.658242188*h__1*h__2)) = 72.05}, {h__1 = 0 .. 1000, h__2 = 0 .. 1000, k = 0 .. 1000})

fsolve({53.50000000+26.50000000*exp(-34.55372656*k*h__1*h__2/(8658.242188*k*(h__1+h__2)+8.658242188*h__1*h__2)) = 56.15, 53.50000000+26.50000000*exp(-17.01469059*k*h__1*h__2/(8658.242188*k*(h__1+h__2)+8.658242188*h__1*h__2)) = 64.1, 53.50000000+26.50000000*exp(-8.094482719*k*h__1*h__2/(8658.242188*k*(h__1+h__2)+8.658242188*h__1*h__2)) = 72.05}, {h__1, h__2, k}, {h__1 = 0 .. 1000, h__2 = 0 .. 1000, k = 0 .. 1000})

(5)

``

it would be great if you could explain me why maple cant numerically solve this set of 3 ecuatios for h1, h2, k being my unknows terms, it works with an exponential function. there is a way to solve it?

thanks..
 

Download heat_ecuation.mw

Also always i try to solve that strange ecuation maple give me like an aswer in terms of Z, with an integral on it, i don't get nothing. I evaluate a point with right click, evaluate a point. Choose the variables value, and the time at wich i want to evaluate the ecuation. Later it give me the same strange solution but with the variables replaced for the values i gave to the program. Then i choose right click, aproximate, and don't matter if i select 5, 10 or more digits, it launchs to me the same strange ecuation, not the number i'm waiting for.

So that i want to know now is how can i manipulate a numerical solution to a diff ecuation.

For example if i got this system of diff ecuation

dx/dt = (-k/m) (x srqt(x^2 + y^2))

dy/dt = (-k/m) (y srqt(x^2 + y^2)) - g

Which is a projectile trajectory taking care of air drag force, and i want to know at which initial speed i have to fire the body to reach 300m if it is fired from the initial point at (0,0), and it's initial velocity on x is 3 times it's initial velocity on y (initial launch angle condition)? 

How can i input that conditions to make maple give me a solution for something like that?

Hi, i got this diff ecuation to answer. It is an ecuation of a proyectile, it give the angle of speed with respect to time. I typed it on wolfram; got an aswer. When i tried on maple without any initial condition, got a good answer. With an initial condition, didn't get nothing. I even tried with the interactive solver, got a 'can't solve'. So this is my ecuation:

dx/dt = -[g cos(x)] / [vo - g t sin(x)]

I had to find the initial coindition for myself after that, i supect what the problem is. Maple gave an aswer in terms of arctan (y,x), and i think, based on my own try, it doesn't know how to find a term if it's inside the arctan function, typed in that way. So i rewrited it like arctan (y/x), anf found my constant C1.

 

Another problem i got was about this system of diff ecuations:

dx/dt = (-k/m) (x srqt(x^2 + y^2))

dy/dt = (-k/m) (y srqt(x^2 + y^2))

Again was wolfram who gave me a good answer to the problem. So you are going to understand im thinking wolfram is better on symbolic calculation. Thanks anyway for your help and effort.

 

 

 

 

1 2 Page 2 of 2