新手求问冲击波管喷出的低马赫数supersonic的 边界条件问题

[Ryo]

用paraview显示的速度只在inlet附近变化 其他部分不变化。求指点下🙏

/Users/Air/Desktop/屏幕快照 2017-10-17 上午11.52.10.png

[Ryo]

[李东岳]

网格多少？小的话我给你试试。

[Ryo]

@李东岳 直径3mm长30mm的tube 用Gmesh 细化到0.1mm的网格 谢谢李老师！

图片地址)

[李东岳]

9000个网格？你上传上来我下载试试

[Ryo]

@李东岳 感谢回复🙏 这个是tube的mesh https://pan.baidu.com/s/1miBVoeO
这个tube只是为了简便练习下用的 实际想做有个小的3mm长的inlet然后测冲击波从tube里刚喷出那段时间的密度计测，这个是模型的mesh https://pan.baidu.com/s/1cLIYzO

谢谢李老师

[李东岳]

最好把整个算例传上来我直接运行。

[Ryo]

@李东岳 https://pan.baidu.com/s/1kVd1wl5
Thanks！

[李东岳]

国外打国内的网站特别慢。百度盘打不开。你把网格发送到沃邮箱li.dy@dyfluid.com

[李东岳]

你提供的算例无法运行，报错文件缺失：

dyfluid@dyfluid:~/桌面/tubecase/10$ ./Allrun 
Cloning resolved case from .
Running blockMesh on /home/dyfluid/桌面/tubecase/10/resolved
Running decomposePar on /home/dyfluid/桌面/tubecase/10/resolved


--> FOAM FATAL ERROR: 
Cannot open file "system/decomposeParDict"

    From function int main(int, char**)
    in file foamDictionary.C at line 427.

FOAM exiting

Running sonicFoam in parallel on /home/dyfluid/桌面/tubecase/10/resolved using  processes
Cloning modelled case from .
Running blockMesh on /home/dyfluid/桌面/tubecase/10/modelled
Running decomposePar on /home/dyfluid/桌面/tubecase/10/modelled


--> FOAM FATAL ERROR: 
Cannot open file "system/decomposeParDict"

    From function int main(int, char**)
    in file foamDictionary.C at line 427.

FOAM exiting

Running sonicFoam in parallel on /home/dyfluid/桌面/tubecase/10/modelled using  processes
         line 0: warning: Cannot find or open file "resolved/postProcessing/probes/0/p"
         line 0: warning: Cannot find or open file "modelled/postProcessing/probes/0/p"
         line 0: No data in plot

dyfluid@dyfluid:~/桌面/tubecase/10$

[Ryo]

@李东岳 在 新手求问冲击波管喷出的低马赫数supersonic的 边界条件问题 中说：

Allrun

李老师，

感谢回复🙏

发的文件夹里面的Allrun，Allclean文件是之前tutorial里复制粘贴过来的，忘了给删掉所以不能用，我用的是Gmesh进行网格划分已经划分过了不需要再用blockMesh了，

另外算这个例子的时候感觉运算量没那么大就没有分cpu去算，所以没有设置decomposeParDict。我现在已经给您加上decomposerParDict文件了如果您想用的话改一下里面的cpu个数就行了，如果不用的话直接在终端输入sonicFoam应该就ok。

另外，我用的版本是openFoam Plus的mac版的，如果用什么格式不对的话，按照tutorial里的复制一下报错的地方应该就行。

下面是case的dropbox链接，
这个是sonicFoam https://www.dropbox.com/s/2t2rc8k3vnl9fqp/sonicFoam.zip?dl=0

这个是rhoCentralFoam
https://www.dropbox.com/s/vzgr1r84fvxk1ts/rhoCentralFoam.zip?dl=0

这两个例子离散格式和用的求解器不一样，这些求解器和离散格式 的意思还没有完全搞懂。。不知道是不是这里出的问题。求老师给看一下，什么样的离散格式和求解器才适合这个例子。

Ryo

[李东岳]

我现在下载看看。

[李东岳]

把你的矩阵求解器这样改一下试试：

"(p|U|e).*"
    {
        solver          PBiCGStab;
	preconditioner	DILU;
        smoother        symGaussSeidel;
        tolerance       1e-15;
        relTol          0;
    }

后续还运行不成功的话，我可以帮你继续再看看，不过麻烦把网格做成六面体网格并适当减少网格数量（5万以下） :cheeky:

[Ryo]

@李东岳
谢谢李老师。

还是跟之前一样的错误。。tube里面的气体不往前面流动。
这个是把网格换成六面体的case
sonicFoam
https://www.dropbox.com/s/bb0a0lmzb67h8o1/sonicFoam 3.zip?dl=0

rhoCentralFoam
https://www.dropbox.com/s/m5ayc8xkdd9ffi4/rhoCentralFoam.zip?dl=0

Ryo

[李东岳]

你的压力是19600？是低压么。还是196000？

非常抱歉，你这个网格太多了我笔记本跑不动。为了去掉网格因素，最好使用全六面体网格，你的这个网格还是四面体的。

[Ryo]

@李东岳

感谢李老师回复！

用snappyHexMesh画六面体网格后发现问题解决了！！
但是还有个问题为什么画的网格感觉好多小的瑕疵。。。。而且这还是完全没有细化也没有添加layers的情况。。

还想请老师帮忙看下snappyHexMeshDic。
这个是snappyHexMeshDic：
https://www.dropbox.com/s/ehr4vknk6injnvr/snappyHexMeshDict?dl=0

这个是case
https://www.dropbox.com/s/xzhhygrywtj3g2r/1.zip?dl=0

Ryo

[李东岳]

如果可以评分的话，我给snappyHexMesh及格分。直接换网格生成软件做全六面体网格，snappyHexMesh只是为了开源框架搭建的。

[Ryo]

@李东岳
谢谢 李老师 有点激动真的快哭了。。

[李东岳]

[Ryo]

@李东岳 感谢老师之前的解答现在又有问题想请教一下。

现在把层流模型换成RAS KomegaSST乱流模型后结果也顺利出来了，研究目的是观察冲击波密度场变化，于是在controldict里加上了密度function，为了想观察冲击波和重播前面的涡流更清晰，想把初始场密度调成1kg/立方米，但是密度好像是根据粘度和初始压力等算出来的 不能调初始场的密度（修改0文件夹里的rho文件 结果还是没有变化）。想问下怎么调节初始密度场。