什么“delta_vel_pos_AI”管(G)的意思吗?

2视图(30天)
PY”class=
PY 2021年7月24日
评论道: PY2021年8月10日
的源代码simscape组件管(G), 237 - 241
流速变化百分比
delta_vel_AI = (1 / rho_I - 1 / rho_AI) * mdot_A /区域;
delta_vel_BI = (1 / rho_I - 1 / rho_BI) * mdot_B /区域;
delta_vel_pos_AI =√delta_vel_AI ^ 2 + (0.001 * a_AI) ^ 2);
delta_vel_pos_BI =√delta_vel_BI ^ 2 + (0.001 * a_BI) ^ 2);
和线326 - 328
%动量平衡
p_A - p_I = = + pressure_loss_AI delta_vel_pos_AI * mdot_A /区域;
p_B - p_I = = + pressure_loss_BI delta_vel_pos_BI * mdot_B /区域;
delta_vel_pos_AI和delta_vel_pos_BI是什么意思?
根据文档 刚性管道气体流- MATLAB MathWorks德国 动量平衡,不考虑任何有关疏通的声速流。这个词的意思是什么“(0.001 * a_AI) ^ 2”吗?我试图找到它的两个引用documentaion但是找不到任何相关的。
不应该只是吗
p_A - p_I = = + pressure_loss_AI delta_vel_AI * mdot_A /区域;
所以它匹配动量平衡方程表示的文档?
提前谢谢你。
2的评论
PY”class=
PY 2021年7月27日
谢谢你的回复。
我试着编辑的代码和使用用户定义的块,和相同的输入参数,导致压降非常不同,如果我删除这一项,所以我猜它应该有一个物理意义,而不只是为了更好的数值鲁棒性。
在参考相关方程中的声速对可压缩流的章,但是我找不到这个方程和方程可演绎成这一项。
我期待着你进一步的解释。谢谢你!

登录置评。

接受的答案

亿丰唐”class=
亿丰唐 2021年7月30日
嗨PY,
这个词使用技巧来避免使用 abs () 函数,它可能引发频繁的零交点当mdot很小。使用 sqrt () 制定和价值总是> 0,没有零交点会解决方程时被触发。
后,需要计算abs动量平衡方程 (delta_vel_AI) * mdot_A /区域 。在这里 delta_vel_pos_AI 大约是腹肌 (delta_vel_AI) delta_vel_AI 不是很小。如果 delta_vel_AI 很小或接近0.001 * speed_of_sound, mdot_A 词也会小,对总动量的影响也应该但在平衡。
我很好奇你发现后修改。你能分享.ssc文件和测试模型?
3评论
PY”class=
PY 2021年8月10日
谢谢你的建议。
如果有任何结果与开发团队讨论后,我应该感激,如果你能提供给我更新的信息。

登录置评。

更多的答案(0)

类别

找到更多的在气体模型帮助中心文件交换

s manbetx 845


释放

R2020a

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!