@东岳 东岳老师，问题已解决。非定常计算一开始要求时间步长要比较小，先得到一个较稳定的流场，在别的求解器里可能会先用稳态求解器算一个初场。而我上面的算例一开始时间步长较大，到后面越来越不稳定最终发散。开始用变时间步长让对应的库郎数从1以下慢慢增加即可。

我看了下，确实不是连续的。在ManualInjection.C文件中：

template<class CloudType> Foam::scalar Foam::ManualInjection<CloudType>::timeEnd() const { // Injection is instantaneous - but allow for a finite interval to // avoid numerical issues when interval is zero return ROOTVSMALL; }

我看了一下，可以用cellZoneInject：

Description Injection positions specified by a particle number density within a cell set. User specifies: - Number density of particles in cell set (effective) - Total mass to inject - Initial parcel velocity Properties: - Parcel diameters obtained by PDF model - All parcels introduced at SOI

你试试这个？

PS. 你需要用topoSet做一个cellSet出来toturial里面有

@zhangxc0223 可以用cyclic边界

@东岳 成功了，谢谢东岳老师了

@东岳
在interfaceProperties.C中编程如下
fvc::interpolate(fvc::grad(sigma) - fvc::grad(alpha)*(fvc::grad(alpha)&fvc::grad (sigma))*mag(fvc::grad (sigma))) & mesh.Sf()
将此项添加到速度和压力方程。其中当计算压力项时发散了。

@东岳 好的，明白了。谢谢李老师！

//- Breakage kernel SBre_[kth].ref() += breakagePtr_->g()[i]()*weis_[i] *( DSDPtr_->beta(abs_[i], kth)() - pow(abs_[i], kth) ); //- DSD TO DO implemented //- Coalescence kernel if (coalescence_) { for (int j

上周五自己做的评论，DSD函数还没植入，周一来的时候完全忘了。吭哧吭哧算了好几天，感觉没有破碎的效果呢？

检查代码发现函数还没植入，返回的0..

我目前不清楚怎么用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