CFD中文网

@coolhhh 谢谢您！我尝试做一下。

@coolhhh 我的速度分布是从文件中读取的，射流入口速度是按照下图给出：

@coolhhh 您提到的这个turbulentDFSEMInlet应该跟InflowGenerator差不多，我看是在V版本中的。但是我的速度分布跟时间相关，我想在这个基础上人为的添加随机扰动。

@coolhhh 您好！这里的decayingTurbulenceInflowGenerator是生成的随机湍流入口边界，与时间有关。比如它其中一个版本的源码链接文本。

如题，我想实现一个依时类边界与人工湍流的耦合，比如以下两种。请问各位大佬该如何实现呢？给个提示也行
依时类边界：

type            timeVaryingUniformFixedValue;
fileName        "$FOAM_CASE/150MPa-U";
outOfBounds     clamp;

人工湍流：

value             uniform (600 0 0);
type              decayingTurbulenceInflowGenerator;
direction         1;
LField            uniform 0.00011;//integral scale = width
refField          uniform (600 0 0);
RField            uniform (0.0729 0 0 0.0729 0 0.0729);
Umean             (600 0 0);

@soulx7 按照模版正常使用就行，或者你搜索一下这个边界，有很多相关说明。

@soulx7 就是OpenFOAM自带的无反射边界条件waveTransmissive。但后来我好像把上边界定义为了壁面计算的。

@李东岳 是的李老师，都是涡量输运方程，但不是一个实验室哈哈哈

@尚善若水 后处理可以直接导出涡度，但是我的计算中涡度是个中间量，所以我写出了计算公式，再把它加入到求解器中了

@李东岳 哈哈是的，李老师。我自己提前写了一下，担心不准确，来请教一下各位。
（1）体积膨胀效应中的涡度omega是不是应该定义为fvc::curl(U)呀？
（2）斜压效应中的原式是▽ρ×▽p(密度梯度和压力梯度)，fvc::grad(rho) ^ fvc::grad(p)中间是叉乘吧？
（3）粘性效应中的粘性应力张量tau是不是应该定义为tau=-muEff()*dev(twoSymm(fvc::grad(U_)))呢？
关于涡量和粘性应力张量的定义是要在求解器中重新声明吧？关于上述三个效应的量化，我不确定能否在计算结果中直接导出呢？还请李老师指点。

如下图红框所示，我想在OpenFOAM中输出这三项所代表的场数据，请问该如何实现呢？跪求

如题，请教各位老师同学，我在LES模拟后，对数据进行后处理时，发现一个问题：定压比热的场值经过平均后被缩小了，但是其他场的平均值没有这个问题，请问各位这是怎么回事呢？平均后的峰值(8000)明显低于真实值(120000)，如图所示：

求问各位，以往的OpenFOAM Workshop 的论文集在哪可以找到？在这个网站里没有呢？https://openfoamworkshop.org/previous-workshops/
我想找一个OpenFOAM Workshop 2016年的一些报告，请问有没有知晓如何获取的？

@cccrrryyy 感谢感谢！给我指了一个好的修改方向！我做过RANS模型的对比，确实影响不大，但我看文献里有的是Smagorinsky算的很好，比如文献1，有的就是算的误差比较大，比如文献2，一头雾水 还是得自己测试呀！

@cccrrryyy 您说的这个方法我在文献也看过，不过跟我这个方法完全不同，就没太关注。
我也按照您的方式调整一下特征长度，看看效果怎么样。特征长度我现在给的是管径的10%。您调整后与实验对的怎么样？
我试了of自带的turbulentInlet，效果不好，别的我也没再找过。我看我这个方向的文献大多基于inflowGenerator。我现在觉得问题有两个：一个是湍流强度，就是雷诺应力影响比较大；再一个就是LES模型，我想试试动态Smagorinsky呢，无奈没找到extend版本的，也没再做修改。

@李东岳 李老师，这个是超临界压力下的射流，工质是氮气。具体工况如下：

@cccrrryyy 您好！我算的是可压缩射流，物性变化比较剧烈，具有高密度梯度。射流与环境存在传热过程。
蓝色的线，射流的液核长度比较短，跟实验差别比较大，其他区域还算理想。我正在测试一个Prt=1.0的Case，看看结果。
您说的对！是给的雷诺应力，雷诺应力计算的时候是(UI)^2，这个I应该指的湍流强度吧？

最近做了一些算例测试，更新一下结果，欢迎各位提出宝贵意见！
入口脉动通过InflowGenerator产生，主要修改了入口的湍流强度和湍流普朗特数。有关残差问题，我尝试缩小至1e-9，结果并没有变化。
另外，图中的红色曲线是两篇不同文献的结果，同样使用Smag模型，看起来差距很大。这里附上参考文献2的一些模拟结果，感觉是Smagorinsky模型不太合适呢。
240万网格模拟结果：

330万网格模拟结果：

参考文献2的模拟结果（SM：Smagorinsky；WL：Wale；VR：Verman；SMD：Dynamic Smagorinsky）：

东岳老师，各位同行，我近期也一直在修改参数重新计算，我抽空整理一下，放到论坛，还望各位大佬多指点。

关于OpenFOAM中在采用K-Omega SST模型时，使用/不使用壁面函数的Y+及边界设置问题，可以参考这个链接链接文本，解释的很详细。