我目前不清楚怎么用paraview处理这个
但我感觉paraview应该可以弄，可能你得外挂paraview的代码脚本，我不太熟悉这个

我也忘了..

这个site是什么意思

目前我也不清楚 知道了可以告诉我一下 :mihu:

@youmengtian 能请教一下你是怎么测试tut中的算例的吗？我计算算例的时候先用scalarTransportFoam计算出一个时间段的结果，再使用PODSolver就报错了。请问你是怎么计算的呀

@东岳 谢谢东岳老师

@东岳 谢谢东岳老师！

分离区大小的问题一般是和粘性有关。看一下用的是什么湍流壁面函数，个人经验，湍流壁面函数对分离区的预测通常是不准确的。我计算的时候一般y+小于1，不采用壁面函数（当然在openfoam里是所谓的低雷诺数模型）。

不太清楚这个，只能建议看看源文件了..

@cfdngu 这是哪篇文献，想拜读一下:chitang:

@zwl 这些都是类模板，你需要找到对这些类模板进行特化的源文件。

参考：
thermophysicalModels/basic/psiThermo/psiThermos.C
thermophysicalModels/basic/rhoThermo/rhoThermos.C

打开“设置”——“更新和安全”——“针对开发人员”——选“开发人员模式”—退出

找到 "设置 "-"应用和功能" -下翻打开 "相关设置"下的"程序和功能"--在左边点击 "启用或关闭windows功能" --找到 "适用于linux的windows子系统" -勾选后点确定——退出

打开应用商店搜索 ubuntu ， 会显示几个ubuntu系统，选择一个系统的（我选的是ubuntu 18那个），下载之后打开，等几秒钟，设置自己的用户名和密码。然后关闭。

然后去www.openfoam.com （如果用的是百度,打开前几个网页应该都是可以的，注意区分版本，有v1612，v172,v182）,打开之后选择顶头The open source CFD toolbox下的Download,
选择win10 native,找到大标题下Install OpenFOAM 下的OpenFOAM-v1806-windows10.tgz,点击下载，点进去大概等几秒钟就会自动下载（需要点确定下载就点确定）。在c盘的用户（Users）文件里新建一个Dwonload文件夹，找到你下载的文件，把他拷到
新建的Download的文件夹。

点桌面左下角的开始键，找到ubuntu,直接打开你的ubuntu子系统.
输入命令

cd/mnt/ ls cd c/ //(打开c盘) ls //(会显示你的c盘里的文件夹) cd Users/ //(打开用户文件夹) ls cd Download //(打开之前建的Download文件夹) cp OpenFOAM-v1806-windows10.tgz /home/(ubuntu的用户名)/ //(注意v1806是版本) ls cd ls ls tar xvf OpenFOAM-v1806-windows.tgz ls ls OpenFOAM ls OpenFOAM/OpenFOAM-v1806/etc/bashrc bashrc .nano

按方向键切到最后，再空两行，输入 source OpenFOAM/OpenFOAM-v1806/etc/bashrc
按ctrl o(不是0)，然后enter键保存，然后ctrl X退出。
然后关闭终端
打开一个新的终端，输入

cd OpenFOAM/OpenFOAM-v1806/tutorials 选择你想运行的算例测试一下，安装完毕。

参考视频链接https://v.qq.com/x/page/b0564p64f68.html

应该有一个WSGG的开源模型，你可以找一下，你可以以此为参考，改成你自己的模型

thermodynamics { Hf 0; Sf 0; CpCoeffs<8> ( 1000 -0.05 0.003 0 0 0 0 0 ); } \endverbatim The polynomial expression is evaluated as so: \f[ Cp = 1000 - 0.05 T + 0.003 T^2 \f]

所以应该是1

感谢！已更新 @D-Benjamin :huahua: http://www.cfd-china.com/topic/1441/openfoam用户指南-勘误/2

谢谢两位老师的关注，我会再多孔介质的资料先学学看，谢谢

@东岳 谢谢

inletOutlet边界条件针对流入和流出的情况需要作区分，比如流出的时候是零法向梯度，流入的时候需要指定某个值。判断流入和流出的信息是通过phi来识别的。OpenFOAM中边界处的phi在大于0的时候，全部流出。因此对于这一类条件，需要指定phi，你的案例中，指定的是phi.air

@东岳 i其实是第一个网格，n是最后一个网格
也就是
先求每一个网格中的alpha * x 再对整个计算域求和
然后求和所有计算域中alpha的值
两者做比值

不好意思 数据在这里
0_1542843045909_48db69f6-3dde-40bc-9b8d-70c1b51caa28-image.png Dd [um] Number
0.00000220 17200.00000
0.00000439 71000.00000
0.00000659 97900.00000
0.00000878 84400.00000
0.00001100 53900.00000
0.00001320 35000.00000
0.00001540 23900.00000
0.00001760 15300.00000
0.00001980 12200.00000
0.00002200 9060.00000
0.00002410 6550.00000
0.00002630 5540.00000
0.00002850 4090.00000
0.00003070 3440.00000
0.00003290 2530.00000
0.00003510 2210.00000
0.00003730 1960.00000
0.00003950 1490.00000
0.00004170 1240.00000
0.00004390 938.00000
0.00004610 860.00000
0.00004830 743.00000
0.00005050 652.00000
0.00005270 550.00000
0.00005490 427.00000
0.00005710 339.00000
0.00005930 282.00000
0.00006150 268.00000
0.00006370 194.00000
0.00006590 208.00000
0.00006810 171.00000
0.00007030 131.00000
0.00007240 106.00000
0.00007460 125.00000
0.00007680 76.00000
0.00007900 66.10000
0.00008120 52.80000
0.00008340 63.80000
0.00008560 23.50000
0.00008780 42.40000
0.00009000 34.10000
0.00009220 33.50000
0.00009440 18.30000
0.00009660 21.60000
0.00009880 17.70000
0.00010100 14.00000
0.00010300 10.40000
0.00010500 10.30000
0.00010800 10.20000
0.00011000 6.71000
0.00011200 0.00000
0.00011400 3.29000
0.00011600 6.51000
0.00011900 0.00000
0.00012100 0.00000
0.00012300 0.00000
0.00012500 0.00000
0.00012700 0.00000
0.00013000 0.00000
0.00013200 0.00000
0.00013400 0.00000
0.00013600 2.98000
0.00013800 2.96000
0.00014100 0.00000
0.00014300 0.00000
0.00014500 0.00000
0.00014700 0.00000
0.00014900 0.00000
0.00015100 0.00000
0.00015400 0.00000
0.00015600 0.00000
0.00015800 0.00000
0.00016000 0.00000
0.00016200 0.00000
0.00016500 0.00000
0.00016700 0.00000
0.00016900 0.00000
0.00017100 0.00000
0.00017300 0.00000
0.00017600 0.00000
0.00017800 0.00000
0.00018000 0.00000
0.00018200 0.00000
0.00018400 0.00000
0.00018700 0.00000
0.00018900 0.00000
0.00019100 0.00000
0.00019300 0.00000
0.00019500 0.00000
0.00019800 0.00000
0.00020000 0.00000
0.00020200 0.00000
0.00020400 0.00000
0.00020600 0.00000
0.00020900 0.00000
0.00021100 0.00000
0.00021300 0.00000
0.00021500 0.00000
0.00021700 0.00000
0.00022000 0.00000
0.00022200 0.00000
0.00022400 0.00000
0.00022600 0.00000
0.00022800 0.00000
0.00023100 0.00000
0.00023300 0.00000
0.00023500 0.00000
0.00023700 0.00000
0.00023900 0.00000
0.00024100 0.00000
0.00024400 0.00000
0.00024600 0.00000
0.00024800 0.00000
0.00025000 0.00000
0.00025200 0.00000
0.00025500 0.00000
0.00025700 0.00000
0.00025900 0.00000
0.00026100 0.00000
0.00026300 0.00000
0.00026600 0.00000
0.00026800 0.00000
0.00027000 0.00000
0.00027200 0.00000
0.00027400 0.00000
0.00027700 0.00000
0.00027900 0.00000
0.00028100 0.00000
0.00028300 1.84000
0.00028500 0.00000
0.00028800 0.00000
0.00029000 0.00000
0.00029200 0.00000
0.00029400 0.00000
0.00029600 0.00000
0.00029900 0.00000
0.00030100 0.00000
0.00030300 0.00000
0.00030500 0.00000
0.00030700 0.00000
0.00031000 0.00000
0.00031200 0.00000
0.00031400 0.00000
0.00031600 0.00000
0.00031800 0.00000
0.00032100 0.00000
0.00032300 0.00000
0.00032500 0.00000
0.00032700 0.00000
0.00032900 0.00000
0.00033200 0.00000
0.00033400 0.00000
0.00033600 0.00000
0.00033800 0.00000
0.00034000 0.00000
0.00034200 0.00000

从代码上解释更好说：

scalarField U(-phi/patch().magSf()); //其中的phi你要给定（指定通量），然后即可计算出速度U的大小（模） vectorField n(patch().nf()); //面的法向 if (phiName.dimensions() == dimVelocity*dimArea) { // volumetric flow-rate operator==(n*U); //速度 } else if (phiName.dimensions() == dimDensity*dimVelocity*dimArea) { const fvPatchField<scalar>& rhop = patch().lookupPatchField<volScalarField, scalar>(rhoName_); // mass flow-rate operator==(n*U/rhop); //速度

你看下/combustion/fireFoam/LES/flameSpreadWaterSuppressionPanel的速度边界条件设置。

还没见过教材里面讨论low-Mach-number assumption，气动声学这块我还没研究过 :zoule: