CFD中文网

@东岳 十分感谢，看了一下我的设置还是以物理核心为准吧，而且测试我感觉，多几个核心的速度几乎没有什么差距

@李东岳

常见的问题是求$\mathbf A \cdot \mathbf x=\mathbf b$中的$\mathbf x$，可以有直接解法（如LU），也可以用迭代法（如Gauss-Seidal），求逆矩阵肯定是下下之选，原因有几点： 操作太多，$~O(N^3)$; 稀疏矩阵的逆矩阵不是稀疏的； 你要求得$\mathbf x$还要再乘一次，$~O(N^2)$ 逆矩阵求了之后没有用，因为CFD方程组是非线性的，每次求的矩阵$\mathbf A$都不一样，所以你求得的逆矩阵只用一次，下次迭代又得重新求解。如果是针对线性问题，每次针对不同的$\mathbf b$求对应的$\mathbf x$，用逆矩阵还勉强说得过去。 现代的CFD程序似乎迭代法用得更多，JFNK, IDR(s)，最不济还有CG和BiCGSTAB.有的连矩阵都不需要组装，何苦去求逆...

自己写的代码终于可以出结果了，但是用cgns存坐标数据的时候发现只能存一半，很奇怪，代码举例如下
三维的面网格，节点数2x11，因此是1x10个网格, 坐标储存在x,y,z数组中

real, dimesion(2,11) :: x,y,z CALL CG_BASE_WRITE_F(INDEX_FILE, BASENAME, 2, 3, INDEX_BASE, IER) ISIZE(1,1) = 2 ISIZE(2,1) = 11 ISIZE(1,2) = 1 ISIZE(2,2) = 10 ISIZE(:,3) = 0 CALL CG_ZONE_WRITE_F(INDEX_FILE, INDEX_BASE, ZONENAME, ISIZE, Structured, INDEX_ZONE, IER) CALL CG_COORD_WRITE_F(INDEX_FILE,INDEX_BASE,INDEX_ZONE, RealSingle, & 'CoordinateX', x, INDEX_COORD, IER) YZ坐标同上

读取时代码为

real :: node(2,11) node = 0.0 call CG_COORD_READ_F(INDEX_FILE,INDEX_BASE,INDEX_ZONE, 'CoordinateX', RealSingle, (/1,1/),(/2,11/),node,IER ) write(*,*) node

node的读取结果为

x(1,1) x(2,1) ... x(1,6) 0.0 0.0 ... 0.0

请问有没有大佬能指点一下

@bart 这部分可以从湍流模型的假设和封闭性开始入手分析，我也不是特别懂，而且内容很多。0_1544074389927_第三章，湍流模型.pdf 附件是中科大fluent 讲稿，和fluent 帮助文件类似，可以参考。其他部分可以参考陶文铨数值传热学第九章湍流模型部分。

用有限差分法得到的二阶精度流场变量和用有限体积法得到的二阶精度流场变量，再分别求涡量场，用两种方法得到的涡量场是否都具有二阶精度或者一阶精度？我不太清楚这其中的原理。你说不会降阶，能否具体说一下原因？

@李东岳 安装ansys 时就可以同时安装cfx么？按照我经验，目前的ansys学术版里面可以直接安装cfx，不需要任何额外的操作。

在SCI里面用云图和矢量图解释不能做正统的解释，并且云图很容易造成视觉假象。一般是用数据图。

0_1540861155976_捕获.JPG

比如这个图，就是分析不同的阻力模型系数预测的Cd系数，可以看出他们有一些不同。接下来在文章中就可以说

由于Ma and Ahmadi模型没有考虑液体粘度的作用，其预测的阻力系数相对于其他模型预测的要低

单纯给一个云图，比如下面这个图，可以分析下涡的存在、涡的特征、相分数分布的均一程度等。

重要的是要能从理论上解释CFD预测的现象

如果不能从CFD理论解释plot以及各种云图的形状以及差异。那应该是CFD模型理解的不够深。比如有些传输模型，湍流耗散率比较大的时候，会导致产生一个较大的源项。这样在看云图的时候，就可以说

从图中可以看出，A区域湍流耗散率比较大，进而导致比较大的源项，从而变量A在此区域浓度较大。和xxx模型是一致的

比如一些模型会增加一个径向的力，那么在考虑这些模型的时候，你的浓度分布在颈向上应该分布不均：

从图中可以看出，由于考虑了径向的xxx力，其将变量A从壁面向管壁中心push，导致一个径向的浓度分布，和实验值更加温和。相比较没有考虑xxx力的情况，浓度分布差异较大。因此建议在通过CFD分析的时候开启xxx力模型获得更好的结果

@astremon 诶，还没有朝这个方向想过，以前一直udf里面只有只有函数和常量......

@李东岳 谢谢:threaten: :qinqin:

时间离散格式应该是计算时间导数的格式，模板没有限制；而推进格式应该要求模板是单侧偏心的。http://www.scholarpedia.org/article/General_linear_methods是一般的时间推进格式，偏心模板最多也就是隐式格式。

我推了一下，好像是这个$\nabla\cdot(\nabla\cdot(\bfU\bfU'))$，你验证下？

@kasitefy
非常对。

另外一个问题是这种加和的数值积分忒慢了

@禾斗石开 fluent

@bestucan 感谢回复

李博继续多出好文！

你两相间的传热项选的什么？两相的温度是怎样的？

@hustzzw 微信图片_20221101093137.png
就像这个图中里面圆柱区域和外面圆环柱区域

@星星星星晴 好的 我搞了一个

@bestucan 嗯嗯，感谢大佬。检查过，check和Quallity其实都没有问题，再次感谢~

@XIN 请问您解决了吗，我也遇到这个问题了