非稳态流体网络方法在发动机空气冷却系统中的应用探索.pdfVIP

对方程(10)积分得到 流路的简化网络图，包括14个节点(4个边界节 r=r(石，L碍，i，r2，K) (11) 若有M个元件汇流到同一节点i，节点总温可由下 式得到 ／●、 巩∽ 、J 彳 = ，L舰N 砭、J (12) 埘∑艄 踟 肘∑心 ‰ 其中M为汇聚到节点f的元件总数，N表示元件。 图3计算模型空气流路网络图 Fluidnetworktheair model Fig．3 of system 气体状态方程 对气体．应用理想气体状态方程 ，，。 点，lO个内部节点)和14个元件。表2为进出口 ●-’l P=pRT 、’ 已知条件，表3为初始时刻部分节点的设计压力。 其中R为气体常数。 需要指出的是由于模拟的主要目的是确定沿程压 1．2压力修正方程组 力的非稳态特性。所以未列出太多温度条件，但 是能量方程是参与了迭代计算的。 对每一元件由两时间步间流量变化关系有 Am2】『，lf+矗f一，％ (14) ‰：硒罴之黜燃量气m。既 5 将公式(5)代入(14)得到 节点号码 ／盘；b／x10Pa 温度，l( m什Af(pin，Pout，，in，reut，f，KI)一 (15) mf(pm，pout，‘n，rout，f，K1)=Am 另一方面，由全微分关系可以得到 (16) 址∑》n+∑缸t 表3部分节点压力初始设计值 ．!垫芝i坐虫雹型P!!：!唑：型：!!!!磐竖12丝! o．5 Pa 节点编号压力，×l Pa节点编号压力／x105 厶A，(‰，‰，‰，‰，tKt)一 五(pi。，‰，‰，，椰t，f，蜀) (17)4结果及分析 对所有元件建立上述关系式，得到线性压力修正 方程组，用CROUT方法求解方程组得到各节点的 图4、图5是分别是节点3和节点4的非稳态 压力修正值却f，各节点新的压力值为 压力特性模拟结果，从图中可以看到这两个节点 Pi=臃十衄 (18) 压力的变化趋势非常相似，先是在开始的小段时 方程组(17)中的系数用到压力，所以要进行迭间内保持稳定(后面所节点都有这个现象)，这可 代求解，直到满足如下的总压残差百分比和流量 能是由于非稳态模拟的初始参数还没进入正常的 迭代造成的，而后压力迅速增大，大约在0．002秒 残差百分比计算精度要求 达到最大值，而后开始变小，变小的速度比之前 IAmI I△口I