CFD中文网

@dzw05 湍流模型中有一个运动黏度nu，正常情况下在后面赋值即可。我当时在网上查询如何提取升力系数和阻力系数时有篇帖子提示要在nu后面先再输入一个nu字符再赋值，由于之后的计算我都用的是这个错误的湍流模型文件，导致阻力系数求解一直出现问题。后来一个一个文件排查的时候发现是这里的问题。

DEM加速比很高
CFD好像一般

@两月三年 good ，好主意。我试试哈。

各位专家好：
在多孔介质模拟的过程中，我发现了一个问题，虽然在发表文章中省略了。但是确实是一个直观的问题，我的模型，RANS方法中，不管采用哪一种，在多孔介质的近壁面处的拟合的速度与实测值总是具备一定的差距。这种差距我分析认为是在较高阻力降的多孔介质环境下，气流在流经多孔介质表面时，一部分流体会因为高阻降而产生流向变化，形成比较复杂的近壁面流场。而在这个近壁面流场中，因为Fluent默认是采用的多孔介质域或者是多孔跳跃面，应该是对这层近壁面没有处理的，依然按照充分发展层处理，造成这层近壁面的流速（甚至是流向）与实测值都有较大区别。
虽然能够较好的拟合实时压降（Endo based equations），但是如果在近壁面流场不符合实际的话，那么很显然对于最后的表面覆尘的预算难以符合实际情况。不知道大家对于这种多孔介质的近壁面流是否有过研究，或者能够有比较好的解决思路？欢迎大家提出宝贵的想法和建议，感谢，再次感谢！
下图为提出新模型对于滤管外壁的模拟和实测结果对比：
0_1544421734755_compare.png

可以用输入边界法造波，比如IHFoam，waves2Foam等，现在的of也自带边界输入法造波了，只是功能还比较单一。

之前一直帮老师干活，好久没算东西了~接下来看看mfix代码吧，想加点东西进去。

二维的vorticity只有一个方向 实际上从涡的形成来看 二维似乎并不适合一些涡量的模拟
但是对于一些简单的流场 实际上2D具备非常大的优势 他的处理简单 很多算法很好自编程实现
并且网格可以完全mapped 最大的优势 还是在于后处理的一目了然
这就是为什么大多数理论研究还是集中2D 但是2D可能不是太适用Fluent Fluent主要接轨工程实际问题

@东岳 对，后来查看了一下userguide 中的solver部分，里面PCG的全称是preconditioned conjugate gradient，就一下明白了

压差是分析很多工程设备的关键参数，在之前，主要都是针对多孔介质的Ergun模式来处理这一块。
也跟试验对比很好。但是，最近新开发的一个模型，虽然与试验的结果对比比较好，但感觉总是有问题。
于是不知道CFD中静压值到底是怎么得出来的？
对于一块多孔介质好像比较好处理，但是如果这个时候多孔介质表面沉积了灰，这个时候这块的整个压降会升高，然后灰越堆越多整个表面的压降都会升高。然而，这个局部压降和表面灰的分布十分相关，造成整个表面的这个灰分的压降相差一千多Pa。而最后一般的检测压降点，都是压差变松器的两个接口，一般都是设置在容器外壁的。但是如果我模拟检测点和试验一样的时候，这个时候的压差比我的表面计算压差平均值要小得多。（压差是1200多Pa，但是我的多孔介质有800多Pa，加上平均粉尘层的压降有640多Pa，合起来应该有1400多Pa，而不是1200多Pa.）
这就产生了一个问题，这个通过N-S方程求解的压力项到底是怎么处理的?
当他通过一层多孔介质的时候，多孔介质内部有压力源项，而在外部没有，那么当他气体从多孔介质外向内流通的时候，他的这层压力和速度是怎么计算的？如果这个时候我再添加一个源项在多孔介质层，而这个源项在多孔介质外是没有的，那么这层额外的源项是不是也会通过增加流体的静压来表达呢？（通过伯努利方程似乎很好理解，但是N-S是怎么实现的，他的这个\Delta P 是怎么获得的？）

这个问题我解决了，就是在周期性关联的时候，在edit periodic/periodic veteric 时要在中心轴上把所有的veterivic都要关联。即一个点要点两次

在使用Fluent 6dof+overset+隐式更新时，发现写出的motion history中的重心坐标与实际明显不符，差的很大，而在udf中通过DT_CG得到的却是正确的重心坐标，不知道有没有人遇到过这样的问题，是什么原因导致的呢？

@冬天大太阳 请问博主找到设置方法了吗?找到了，能不能向您请教一下，有偿答疑，谢谢您，打扰了。

adept

如果我没记错的话，流体中文网上的ICEM版块里有人画过4119桨的结构网格，也传了文件在上面，你可以去上面找一下看还在不在。（耐心往后翻一下）

多谢，已更新

@东岳 内网格，尝试着接着划吧

@东岳 谢谢！

0_1543372305704_捕获.JPG

https://www.ansys.com/-/media/ansys/corporate/resourcelibrary/article/aa-v2-i2-coupling-momentum-and-continuity.pdf

刚看了下，Fluent那面也只是速度和压力耦合求解。能量抛在了外面。

但是另一个文档速度压力能量放在了一起 ： http://www.bakker.org/dartmouth06/engs150/05-solv.pdf

速度压力耦合求解在高超那面用的多，如果分析N-S方程的数学特征（双曲抛物），你需要把u, v, w, p当成矢向量，然后求特征值特征向量。接下来求解这个矩阵就是耦合算法。我目前只是了解到思路，没亲自植入过。可

@gemini 谢谢您，服务fluent版本跟调试用的是一致的，我去试试修改成您说的这样

想解决IV型激波，入射激波如何通过边界条件设置产生