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

Ryo

@李东岳感谢老师回复 🙏

Mac上还是没有找到什么原因wmake出错，但是在虚拟机（Linux Mint）上在按照openfoam wiki上添加温度的教程编译成功了。
发现在createfields.H里加Oobject::MUST_READ,IObject::AUTO_WRITE，和在/systym/controlDict里面加
0_1513060706790_屏幕快照 2017-12-12 15.37.48.png
效果一样都只是根据初始p，U，粘度等来计算密度，不能直接读取0文件夹里面的rho文件。。。
这是0文件夹里面的rho文件虽然设置的是1但是计算后在paraview上显示是0.23
0_1513060816206_屏幕快照 2017-12-12 15.31.18.png
怎么才能跟fluent一样直接调节初始场密度呢。
麻烦老师在看下🙏

Ryo

@李东岳
老师好！发现openfoam是不能直接进入密度场的因为openfoam用的是kinematic viscasity.fluent用的是动力粘度。
0_1513326954486_屏幕快照 2017-12-14 16.26.35.png
所以感觉还是得调节p，之前的压力边界条件我设置的是outlet为0.2大气压，inlet为zerogradient。现在调换后inlet为0.2大气压，outlet为zerogradient，initial为1个大气压。然后按照
0_1513327734658_屏幕快照 2017-12-15 17.46.38.png
更改fvscheme后达到预期的冲击波外的区域密度场都1左右，但是发现冲击波前面的涡很微弱。
这是实验得到的图像。冲击波前面有很明显的蓝色的涡。
0_1513328028553_屏幕快照 2017-11-20 13.59.32.png
这是openfoam得到的图像，冲击波前面基本没有涡
0_1513328098516_屏幕快照 2017-12-15 17.52.52.png

李东岳

volScalarField rho
    (
        IOobject
        (
            "rho",
            runTime.timeName(),
            mesh
        ),
        thermo.rho()
    );

OpenFOAM里面的密度是通过thermo.rho()更新的，因此你不能直接给密度。

更改fvscheme后达到预期的冲击波外的区域密度场都1左右，但是发现冲击波前面的涡很微弱。这是实验得到的图像。冲击波前面有很明显的蓝色的涡。

云图看起来不是很客观，最好出plot图看看

Ryo

@东岳 Hi!
又要来麻烦老师了🙏

之前的涡的存在问题已经解决。原因是因为压力边界条件设置问题，还有发现k和epsilon之前按照公式算的是10的5次方级别的数字但是我把他调节成0.5也能算而且精度更准，这个还不知道什么原因。。

现在（其实从之前）就有个关于paraview的问题，把data放到paraview里面显示rho跑第一遍的colour legend和先显示p跑一遍在切换到rho显示的colour legend有些许的不同。

这是显示rho跑第一遍的图 colour legend只显示到0.23
0_1542768671839_屏幕快照 2018-11-21 11.46.20.png

这是先跑完一遍p在切换到rho， colour legend就可以显示到0.06
0_1542768794615_屏幕快照 2018-11-21 11.46.33.png

求老师给点指示

李东岳

什么是“先跑一遍rho和先跑完一遍p”...

Ryo

@东岳
sorry 😂
先跑一遍rho：把data放到paraview里面选择rho然后点击play
0_1542872176845_屏幕快照 2018-11-22 16.10.30.png

先跑一遍p在调节到rho：把data放到paraview里面先选择p然后点击play在从p调节到rho
0_1542872217819_屏幕快照 2018-11-22 16.09.03.png

想说的意思可能就是第一遍play显示的数据图像和play完第一遍后切换到其他性质（比如从密度切换到压强，压强切换到密度）后在切换回来显示的不一样。一样的data为什么会显示的不一样呢？问题解释的可能有点傻，请见谅。。

李东岳

0_1542931053530_捕获.JPG

你点击一下这个试试

Ryo

@东岳
懂了，多谢老师。

Ryo

@东岳
又来麻烦老师了。之前用snappyHexMesh做的网格感觉不好现在想换用icem cfd来做，弄了好几天发现喷口旁边座位阻碍物的圆柱的网格有大瑕疵（表面的都很好），找了很多视频里面都没有找到该怎么去掉圆柱里面这些不需要的网格。
另外icem cfd做好的网格怎么导入openfoam和paraview现在也没找出来方法。请老师给点指示。
0_1543480873745_Inked捕获_LI飞.jpg