CFD中文网

@Leon-Zhao 感谢大兄弟

谢谢！有效！

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

请问有前辈知道gmesh如何能划分如图所示得网格吗，图1是图2的block

block-tuya.jpg model-tuya.jpg

目前算例移植到计算机集群上计算，模型是水面上漂浮的平板。大约计算到0.8s时出现这个报错，查找原因是内存分配问题；
8fda3931-8566-4c1d-96e3-08c9fb3e932f-image.png
有查看了相关的log文件，文件中给出的报错是如下所示
d5c39409-328d-471d-be7b-801f479d1faa-image.png
为什么算例运行中会出现nu不被定义的报错呢？是因为计算过程中，网格出现问题嘛？

@李东岳
好，谢谢李老师。。。
等到了李老师回复

这在以前的openfoam里面是出现过这种问题的。我详细的跟踪过twoPhaseEulerFoam的算法演变。你这个求解器应该是ESI那面的吧？可以问问他们什么时候能把alpha、rho搞进去

挺好的,ratio很小了

@硫酸亚铜 大神，又看到你了，我要膜拜你了。

@cfd-china 这是我的毕业课题:tongue_out:
我采用了全六面体网格，起初本来打算学习OpenFOAM，但是编程基础不是很好，所以最后利用Fluent来计算。

我现在做的模拟是根据一篇外文文献，里面有详细的几何模型数据，还有一些实验数据，文章实验数据中有8个值左右，但是实验值均在0.72浮动，没有多大变化，所以我先比较了其中一种情况。

因为是毕业课题，所以担心被答辩老师质疑。比如网格不多时模拟的结果好，为什么要用网格多的来模拟。我就懵了:confused: 。（我的工科生，流体力学知识仅仅学过工程流体力学）

0_1462072825790_14.jpg
调大此值，如果还出现警告，建议改为标准初始化

诸位大佬，XiFoam中所使用的燃烧模型，是基于反应变量的湍流燃烧模型。
b=1:未燃混合气b=0:已燃混合气
是否可以考虑在其中增加详细的化学反应机理，进行计算。

@bestucan :mihu: 确实是直接删除block就行，俺发现好多操作都跟人家一样。得到得结果却不一样，可能因为我用的不是正版

@李东岳

常见的问题是求$\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.有的连矩阵都不需要组装，何苦去求逆...

@lrl3512 要不您加一下我qq 2559727310

@Rachel0096 接口部分内容如下：5252f072-1693-4abd-b204-3e6712b25d77-image.png 674fe43f-b3e0-49ad-a2de-daf068772bc0-image.png

@赵一铭 不是用软件，是编程计算，现在的问题是一个不收敛，和出现负值

很普适性的问题。如果你说的是如何判定非稳态的收敛性的话，可以监控一些参数，如果稳定波动一定周期，就可以认为收敛了。

@李东岳 哈哈，其实是从nodebb的帮助中抄来的:cheeky: