CFD中文网

最近在做过滤器的仿真计算，计算后对得到的压降损失的来源有些疑惑，请问大家我只打开了粘性模型Realizable k -ε，并用了增强壁面函数，计算得到的压降损失都来自哪些？目前我学习到的应该是包括了流动过程中流体微团的变形所损耗的压降，但不知道还有没有其他压降损失来源，主要要的是壁面的摩擦阻力造成的压降损失有没有考虑进去，为此求教下大家

http://c.biancheng.net/view/2226.html

this就是当前对象的别称，有哪些成员，就看它所在的函数是哪个对象的成员函数了。
去定义这个对象的头文件查，与之同名的.c文件有各成员函数的定义。

不好判断 但感觉就是设置的问题 如果可以上传算例并且网格数量小的话 我可以给你看看

头晕眼花

sad
:crying:

有个方向，叫混沌之类的。详细可以了解下蝴蝶效应和非线性的一些内容。

均匀各向同性湍流，其雷诺应力矩阵非对角项等于0 竖向非均匀各向异性湍流，比如槽道流（水平向均匀），只有$ \left \langle u_{1} u_{2} \right \rangle $有值，其他非对角项均为0。Pope的书中也给出了解释。 现在问题是，均匀各向异性湍流，其雷诺应力矩阵非对角项是否为0？如果参照Pope书中类似推导，在空间3个方向都是均匀的，理论上$ \left \langle u_{1} u_{2} \right \rangle $，$ \left \langle u_{1} u_{3} \right \rangle $，$ \left \langle u_{2} u_{3} \right \rangle $都为0吗？

《Pope, S.B., Turbulent Flows. Vol. 12. 2000: Cambridge University Press.》89页
e0bada21-7ac0-4ae4-acc3-3319a012dfc7-image.png

$\mathrm{TV}=\sum_{i}^{}|\phi_{i+1}-\phi_{i}|$

可能是超算上有不只一个节点，所以节点间通过IP地址交流信息。
而linux系统里有一套“权限”和“环境变量”的继承规则。你加载 openfoam 环境变量的的那句应该是放在家目录的.bashrc了（就是source ~/OpenFOAM/etc/bashrc，大概是这句）

而 miprun 再通过IP调用其他节点的时候，这个网络通信的用户权限和环境变量，应该没有继承你的用户名的环境变量；而是系统的环境变量。所以不包含你的用户目录下.bashrc里的东西。

把你家目录那句 source 放到 /etc/profile 试试 。这是系统范围的环境变量。得是绝对路径：/home/<用户名>/OpenFOAM/……

当然也可能得在各个节点上都配上环境变量、装上 openfoam，这细致活得慢慢调……

部分求解器出现上述状况（编译成功，但是不能算），还有一部分求解器完全正常。

没用过多节点的，但问题是出在环境变量上没跑了。
到这个路径下
~/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin
这都是可运行的二进制文件。

直接./<文件名>，比如
./icoFoam
就知道是没编译好还是环境变量没配置好了。

十分感谢

就是下面这个
$$
\frac{\p T_{ij}}{\p x_i\p x_j}
$$
其实$T_{ij}=u_iu_j$，所以演变为
$$
\frac{\p u_iu_j}{\p x_i\p x_j}
$$
这个公式是可以写成$\nabla\cdot(\nabla\cdot(\bfU\bfU))$这种形式，$\frac{\p T_{ij}}{\p x_i\p x_j}$好像写不出来 :134:

大佬们看看？

@CarmeloSun 谢谢前辈

对于URANS，多了个瞬态项，那每个时间步计算的时均速度（时均动量）都是从第一个时间步开始累加后平均后获得么？

对于$\int_t^{t+\Delta t}$的$\Delta t$，从理论上URANS这个$\Delta t$是任意的，如果不考虑附加的湍流封闭模型如ke之类，$\Delta t$足够的小，就变成了LES，$\Delta t$足够的大，就变成了RANS。然而$\Delta t$并不需要显性的给定。所以对URANS取决定性的还是$\mu_t$

@星星星星晴 明白了.我这个方法确实直接套用不上.并且碰撞模型从同一个cell扩展到周围cells之间的碰撞,还会有其他的问题.我再想想,谢谢了.以后常联络!

有一节叫Advanced Multiphase Macro，也许有帮助。

For most standard UDFs written for multiphase models (for example, source term, material property, profile functions), variables that your function needs (domain pointers, thread pointers, and so on) are passed directly to your UDF as arguments by the solver in the solution process. All you need to do is hook the UDF to your model and everything is taken care of. For example, if your multiphase UDF defines a custom profile for a particular boundary zone (using DEFINE_PROFILE) and is hooked to the appropriate phase or mixture in Ansys Fluent in the relevant boundary condition dialog box, then appropriate phase or mixture variables will be passed to your function by the solver at run time.

There may, however, be more complex functions you want to write that require a variable that is not directly passed through its arguments. DEFINE_ADJUST and DEFINE_INIT functions, for example, are passed mixture domain variables only. If a UDF requires a phase domain pointer, instead, then it will need to use macros presented in this section to retrieve it. ON_DEMAND UDFS are not directly passed any variables through their arguments. Consequently, any on demand function that requires access to phase or domain variables will also need to use macros presented in this section to retrieve them.

JCP
让大家膜拜一下。:sunglasses:

@thegame 感谢老铁，将在下一版本更新

计算的是气液混输泵的气液两相流计算，air-water(25℃)的混合介质，air为离散相给定体积分数和固定的气泡直径，water为连续项，湍流模型为SST k-w，欧拉双流体模型。计算出来之后总感觉残差是不是波动有点大了。
有没有这方面的大佬教教我，是不是设置有问题还是其他的问题。下面是计算得出的残差
47a7ac58-69fe-4db7-a71d-e8ba6fd07ebc-image.png
800b8c63-4b22-47a5-9131-cc03980b1291-image.png
05ca6207-15dd-45ab-bb5e-9628346b8a9e-image.png
e252f689-2fd6-4d82-b46a-77c7453a624a-image.png

跟出口速度有关系，延长出口计算域试试？

udf有问题，你试试这个

#include "udf.h" DEFINE_PROFILE(inlet_velocity, t, i) { real x[ND_ND],y,vel; face_t f; begin_f_loop(f,t) { F_CENTROID(x, f, t); y = x[1]; if(y>0.0393) vel = 29.5; else vel = 29.5 * pow((y / 0.0393), 1/7); F_PROFILE(f, t, i) = vel; } end_f_loop(f,t) }