在Matlab中使用Newton求解Colebrook方程

时间:2017-10-02 01:31:09

标签: matlab equation-solving fluid-dynamics

我制作了这个剧本:

这部分用于创建de变量Re和Ru所以,我可以替换使用另一个。接下来,创建表达式,以便MATLAB可以推导出它。

spark-submit -v --master local ...

这是Newton-Rapson的脚本,它适用于其他表达式。

Re = 152788; % Reynolds
Ru = 4e-05; % roughnesss
syms x a b
dy = subs(1/sqrt(x)+2*log10((a/3.7)+2.51/(b*sqrt(x))), [a, b], [Ru,Re]);
df1 = diff(dy);

当我运行脚本时,迭代永远不会结束,继续无限。获得的数据摩擦系数为0.016

0 个答案:

没有答案