Modeling efflux time from a tank. Equation 11 in the attached document

Question

Eric Gray 2017-1-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/322203-modeling-efflux-time-from-a-tank-equation-11-in-the-attached-document

回答： Sergey Kasyanov 2017-1-27

CHE 535-Computer Modeling in the UG Unit Ops Lab (1).pdf

The equation is first-order nonlinear ordinary differential equation. Says to combine Newton-Rhapson and Runge-Kutta methods to solve numerically.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Sergey Kasyanov 2017-1-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/322203-modeling-efflux-time-from-a-tank-equation-11-in-the-attached-document#answer_252291

在 MATLAB Online 中打开

I think it can be helpful.

You have this function with A, B, C, D are known:

 fun = A * (dHdt)^2 + B * (dHdt)^1.75 + C * (dHdt) + D* H;
 A=number;
 B=number;
 C=number;
 D=number;

Then you define initial condition:

 dHdt0=dHdt0;
 H0=H0;

Then you combine solving of function fun(dHdt) and integrate ODE fun(dHdt,H):

 %create new function
 fun1 = A * (dHdt)^2 + B * (dHdt)^1.75 + C * (dHdt);
 %define interval to integrate
 T=number;
 %define step
 dt = 1e-3;
 dH = [dHdt0,zeros(1,T/dt)];
 H = [H0,zeros(1,T/dt)];
 i=2;
 %integrate
 while i<T/dt+1
     %solving of fun(dHdt)
     dH(i)=vpasolve(fun1 + H(i-1) * D,dHdt);
     %integrate with euler method
     H(i)=H(i-1) + dH(i) * dt;
     i=i+1;
 end
 t=[1:i]*dt;
 plot(t,H);