CFD中文网

厉害厉害了老铁 你回国了么？

同理，下述网页中都需要处理$\sum_f p_f^{t}\bfS_f$为$\frac{1}{V_\rP}\sum_f p_f^{t}\bfS_f$，均已更正，谢谢

http://dyfluid.com/rhoPimpleFoam.html
http://dyfluid.com/rhoSimpleFoam.html
http://dyfluid.com/simplefoam.html

请问有大佬知道如何在Fluent中利用UDF添加Maxwell速度滑移和温度跳跃边界条件吗

在思考一个审稿人的意见，太难了

速度的物质倒数两端乘以密度，再化简，就是时间项和对流项加在一起。

密度就是单位体积的质量，速度的物质导数就是加速度。
所以时间项和对流项加在一起就是单位体积内的$F=ma$，即$F=\rho a$。这就是惯性力吧？

力分为体积力和表面力，本来表面力应该是产生相对于体心的力矩。但是无限小体积的情况下把力臂忽略了（我猜的），力矩造成的转动靠物理规则演化出来而不故意设计出来。所以左边是流体微元变化需要的力，右边是流体微元真实受到的力，两者相等，就推出来流体微元怎么变化的了。

物理意义有助于了解方程怎么来的，但是后来方程形式因为代换化简早没有了当初的模样，只能从大概上理解。像广义的源项，这个和运输方程放在一起好理解。NS方程化简到最后是三个运输方程，运输的物理量分别是密度、速度、熵还是啥。源项，就像一个无源封闭体内进来多少就出去多少，但是有了源，就兴许进来的少出去的多。就是对运输产生了扰动。

要是搁到动量方程上，比如加个源项，磁流体在磁场中受电磁力，或者离心机里受额外的等效重力。

另外粘性造成的应力是和速度有关系的吧？叫本构关系，用来区别牛顿和非牛顿流体的。系数小不打紧，只要湍流涡够多，速度够快。耗散能量的能力就很强。

如题。

好的，谢谢您！我找到问题所在了，模拟稳定之后就可以正常运算了。@ice_flow
我关于Q和Lambda2还有一个问题：
我的入口采用固定流速在一个粗网格上模拟了一下，发现模拟结果就能较为清晰的展现涡旋的结构。85327f57-9836-4fd0-b064-54b24323035f-4370c9d80aec1a40b3dcc5e1703a754.png
可是当我使用空渠道进行数值模拟，得到一组携带湍流特性的流速断面；然后作为添加了阻水建筑物后的模型的流速进口，计算出来Q和Lambda2总是展示出大量的稀碎的涡!
57baf78d-fe53-4675-9544-003d081397cc-e402149df84a5275c2af261e79309689.png
77fc94f6-a4da-444c-b4f6-dc6d68958b5c-54d7b6532cabd4e736203df5b16d0d41.png
我想问一下怎么才能得到清晰明了的涡旋结构图呢？希望能得到大佬们的解答！

@ljpgbl
这个是全貌
0_1539505430619_Clip_20181014_162251.png

@李东岳 主要是磁共振成像设计和基于磁共振图像的心脑血管力学分析

@cccrrryyy 感谢大佬答疑解惑，也谢谢您的分享:140:

@yuhxFoamer
不是老师哈，我最近也在做网格验证，所以可以一起讨论下。
1.最小网格尺寸看起来变化了，网格的质量好像也不一致，看起来不是结构网格。要考虑设置相同的计算参数。
2.几何尺寸发生变化会导致雷诺数、参考面积这些的变化。我看你的计算域在三个方向上都发生了改变，检查一下基于这些几何参数的其他量是不是也相对需要改动。

你可以按照这个方案来检查，首先在原本能够顺利计算的网格上直接加密，看看是否发散。然后使用缩放命令改变网格（不再加密）transformPoints -scale "( 2 2 1)"，设置对应尺寸下的物理参数计算验证。一般来说没有很密的边界层的话可以这样调整，这样不行的话就进一步考虑离散方法和数值精度的调整，先用一阶的，再反复调整这样。

我碰巧最近在看这方面的工作，打算在OKS课进行植入，中文还对不上，但是应该就是楼上说的 转捩模型。我参考的是Menter的文章A One-Equation Local Correlation-Based Transition Model，他们的参考案例也很有意思，T3A系列，单独用层流、湍流都不行，只能用附加转变的模型。下图：

捕获.JPG

考虑最简单的动量交换，液相：
\begin{equation}
\frac{\rd \bfU_\rc}{\rd t}=-\bfA
\end{equation}
气相：
\begin{equation}
\frac{\rd \bfU_\rd}{\rd t}=\bfA
\end{equation}
乍一看，如果$\bfA$是负的，那么会导致$\bfU_\rd$向下走（为负），$\bfU_\rc$（为正）。但是这并不符合物理，考虑一个管子的颗粒，如果颗粒向下走，必然会带动空气同时向下走。但是方程缺不是这种体现，为什么呢？

实际上，上述方程并没有写完整，完整形式应该是这样：
液相：
\begin{equation}
\frac{\rd \bfU_\rc}{\rd t}=-\bfA=-\frac{1}{\tau}(\bfU_\rc-\bfU_\rd)
\end{equation}
气相：
\begin{equation}
\frac{\rd \bfU_\rd}{\rd t}=\bfA=\frac{1}{\tau}(\bfU_\rc-\bfU_\rd)
\end{equation}
如果开始的时候$\bfU_\rc=0$，则为

见过一些，目前也不好解释。:141:

@chengan-wang 你好，请问你说少了一个定义是啥？我好像遇到了和你一样的问题，先谢谢啦。

你看看反应速率有没有 是太低了 还是模型设置有问题

@东岳 那请问东岳老师，在t=0，和t=1时刻，这个函数返回的值是什么意义呢？

@李东岳 前处理只能在Open FOAM里面建模和画网格吗

@huoyanxiadebing 在 说实话，根本不知道选择什么湍流模型 中说：

看到你的帖子，发现你做的都是外流场，对于外部流动，沿表面位置分布的雷诺数大于500000时，通常认为流动状态为湍流；沿障碍物的雷诺数大于20000时认为流动状态为湍流，如果你的模型中雷诺数不是很大的话，你用层流模型计算即可，这就是有时候你用层流模型算的比较准确的原因；对于雷诺数的计算，最重要的一点就是你要选好特征长度，对于没有结构物时，你要以水槽的长度或高度来算雷诺数，加上了结构物，就要判断是否产生了湍流现象，以结构物的某个长度及其附近的速度来计算雷诺数。

前辈您好，请问这里提到的两个临界值是否有相应文献的支撑？

如果要计算防波堤（例如沉箱防波堤和水下防波堤）附近的雷诺数来选取相应的湍流模型（适用于低雷诺数或高雷诺数的），这个特征长度应该取防波堤的高度吗？对于二维和三维情况，特征长度的选取是否有不同呢？

另外，结构物附近的速度是否有较为通用的估算方式？对于防波堤附近的速度，能否用基于浅水波理论的公式Uh=cη进行估算呢（水深h，波速c，自由面高度η）？

7f0a7c02-1beb-40d4-8130-9ba6e44bbb4a-image.png d45554ce-7a55-414b-a667-8ccc7afa71c1-image.png
从上至下依次为网格1，2，3。
如图，翼型不同攻角的Cl验证，使用fluent2022r1在Re=5e5, SST k-w模型进行计算，稳态使用simple求解器，在攻角为5°的情况下，网格1（单一网格）算得的cl为0.479，网格2（使用重叠网格计算）算得的cl为0.531，网格3（单一网格，在网格2的前景网格为基础上拓展）算得的cl为0.475,文献上的cl在0.459-0.483之间，上述网格的正交性都比较好,y+=1。目前遇到的问题是不知道为什么用重叠网格后升力系数会差这么大。
已使用过的方法：
①换不同的湍流模型计算，虽然结果有差异，但都在0.52-0.55区间波动
②改变计算域大小，从12c拓展到20c均无明显变化；改变前景网格区域大小，从d=5c拓展到d=20c（此时background为30c）均无明显变化
③使用不同的fluent版本（19，20，23），无明显差异
实在找不出原因，有没有大神给小菜鸟解惑_/|o