喷射器的模拟，求教

nickchen0625

我做的是一个喷射制冷系统中的一个喷射器模拟，喷射器尺寸是用索科洛夫的方法设计的，画完网格就像附件显示的那样，直接都是用size0.3的四边形网格，设置的是工作气体（3.2Mp，90℃）、引射气体进口（0.38Mp，2℃）都是压力进口，边界设置只用了总压和温度，出口是压力出口，也只设置了总压和温度（0.92Mp，32℃），流体是用R134制冷剂，物性参数用多项式输入，用了基于密度的轴对称隐式求解器，选了energy方程，k-ε方程（默认设置），
问题1：solve-control-solution里为啥没有energy方程选项，
问题2：计算老是发散，改了松弛因子，库朗数，只是把发散往后推了几步，而且发散前往往都有reversed flow in 4 faces on pressure-outlet 6.
time step reduced in 4 cells
absolute pressure limited to 1.000000e+000 in 4 cells on zone 2
temperature limited to 1.000000e+000 in 4 cells on zone 2
turbulent viscosity limited to viscosity ratio of 1.000000e+05 in 11 cells
对于初始化里的东西用的默认的
求解啊，大侠。。。最好能不能qq远程协助下。

浪子轻狂

网格是四边形网格，但计算过程，温度、压力、湍流粘度全部超限了，你这参数设置得多糟糕啊。另外，你这个温度和压力下，134a会不会发生相变？

nickchen0625

浪子大哥指的的参数设置是指哪个方面？温度，压力都是喷射器的工作参数应该是不会变的吧，是不是是湍流粘度的设置问题，我都是默认1的，
至于相变，，虽然进出口都是气相的，其实喷射器中间引射流体与工作流体的混合过程是存在相变的，

[ 本帖最后由 nickchen0625 于 2012-5-17 09:44 编辑 ]

nickchen0625

浪子大哥，我还想问一下，material里的物性参数设置不合理会不会引起那些问题，我是用自己拟合的一个多项式来设置的，这个多项式只适合气相的流体的。

浪子轻狂

默认肯定是不行的，湍流参数自己估算一下，适当修改一下温度压力以及湍流粘性的区间值。另外，既然存在相变，而你这么算，体现不了相变过程，我觉得你这么算不太合理，当然，如果仅有极少一部分气相发生了液化，倒也可以这么简化。

nickchen0625

浪子大哥，我把进出口边界的turbulence的条件都采用turbulence intensity和hydraulic diameter 的设置，水力直径是入口直径已知的，但是湍流强度我不知该如何确定，选的是10%，湍流模型采用了Realizablek −ε模型，收敛情况较之前有所改善，但是energy和continuty的收敛超不好，而且几乎相等，两条线重合的。并且整个计算过程始终有turbulent viscosity limited to viscosity ratio of 1.000000e+005 in 806 cells，reversed flow in 27 faces on pressure-inlet 6.这里的pressure inlet 6指的是引射流体的入口。不知道该改什么参数了，库朗数已经调成0.5了，算了9000步才这样。

[ 本帖最后由 nickchen0625 于 2012-5-17 16:03 编辑 ]

maaijie2009

谁有丙烷四步反应的方程式外加阿累尼乌斯反应速率的各个系数（指数前因子、温度指数、反应活化能等），可以发到我邮箱altanhorwa@yahoo.cn 谢谢。

nickchen0625

没有做过这方面的东西哦

通流

要查错，至少要把流场贴出来。

不过，我想，你的问题是超音速喷管的启动问题。在启动过程中，激波会通过管道（起码是管道中的一部分。用定常算法，是很难通过这一个过渡过程的。所以，我想你应该用非定常计算。等到流场差不多定下来以后，再用定常求解器。所有的参数，基本可以选用默认的值。我遇到过类似的问题。

nickchen0625

好的，谢谢版主，我试试看

nickchen0625

版主，用了基于压力的隐式非定常的求解器，时间步长1s，每步200次，计算了10s，流场的结果如图，喷嘴出口出现了激波，但是在最后的喷射器出口也出现了超音速，也是激波吗？
还有，我用了基于密度的隐式非定常求解器，设置与基于压力的一样的，但是算了3步就发散了，为什么啊。
而且计算过程中老是有reversed flow in 3 faces on pressure-outlet 5.

reversed flow in 1 faces on pressure-inlet 6.

[ 本帖最后由 nickchen0625 于 2012-5-18 11:16 编辑 ]

uesoft

一般都不用喷射制冷，能量效率很低。

nickchen0625

不要沉了呀。

通流

激波都到出口了。再算一会，激波可能就出去了，或者缩回来了。

你的时间步长取得挺有意思的。步长一秒钟。一秒的时间，压力波能够传播300多米。不知你的计算域有多长。就像你在房间大小的区域里找钱包。可是你的每一步都跨出去几百米的话，你怎么能够找到钱包？

nickchen0625

哦，是这样啊，谢谢版主，没做过激波的模拟，我的喷射器从喷嘴出口到扩压器出口的距离是0.2m左右，扩压器的长度是0.048m，那步长要怎么取啊，是不是要满足时间步长乘以步长数量在乘以音速等于传播距离？

[ 本帖最后由 nickchen0625 于 2012-5-19 11:15 编辑 ]