@Wayne 对，电磁场里面会用到

我们都知道，速度可以做Helmholtz分解，用速度势的梯度和流函数旋度叠加的形式，求解二维粘性可压缩翼型绕流中的速度势和流函数相当于求两个泊松方程，但是需要给出边界条件，一般情况下，在速度场已知的情况下，即涡量场和胀量场已知，该如何分别给速度势和流函数的物面和远场边界条件？

@Elibathe康 感谢！

更精准一些的话，D为全导。d后为单变量的函数。:cheeky:

@wwzhao :thanks!

@cfd-china 写文章，做实验。。

谢谢，那个我也看了（但没仔细看，只是特意找了扩散项）。应该也是和Menon一样的结论。因为Menon的文章截图比较方便，才引用他的

@pippowong
可以简单说说这种非常大尺寸的大气模拟和普通的流场模拟的区别么？以及相关的侧重点？

@hangsz Anderson 的中文版啊
0_1464596732984_351878237561634601.jpg 0_1464596737400_157033197177427341.jpg 0_1464596744921_63913991367426126.jpg

@cfd-china 我觉得不是，并不是大于$\Delta$的窦唯大涡。

为何无法插入图片？我试着把公式截图发上来，但是好像不行:sad:

Jacobi是并行性最好的算法，LU-SGS并行性较差。

@搬运工

高斯迭代法处理20阶矩阵就已经非常慢了。通常和多重网格+高斯赛德尔求解器一起处理。并且OpenFOAM这种成熟的商软里面矩阵求解器算法非常精炼。通常由主求解器+预处理+光顺+重组构成。多重网格更是复杂。另外，在编程的时候使用动态、寻址的方法寄存系数。

所以如果你要写一个快速的求解器：

使用外挂包； 潜心研究稀疏线性系统求解理论；我建议花大量时间研读：YousefSaad, 2003。看看大型稀疏矩阵怎么求解的。而不是3/4阶矩阵的高斯塞德尔求解器。:cheeky:

@Elibathe康

直接和金钱相挂钩啊！飞机公司确实需要好好算计算计。哈哈。:cheeky:

@Elibathe康

没用过，不过在CFD中源项计算用过比较简单的5点高斯勒让德积分。精度没度量过。不过足够了...

:sunglasses:

SRF和MRF只是坐标系变换，将旋转问题变成stationary问题，所以湍流模型还是适用的。

但是对于旋转问题或者流线弯曲问题，湍流模型是要进行修正的。

参看： 《Some improvements in Menter's k-omega SST turbulence model》

Navier-Stokes 方程中的对流项写作 $\nabla \cdot (\textbf{u}\phi)$，扩散项写作 $\nabla \cdot (\nu \nabla \phi)$。

对流表示流体微团从空间某点运动到另一点的过程，而扩散则表示流体物理性质浓度变化的过程。

对动量方程而言 $\phi=\textbf{u}$，因此对流项变为 $\nabla \cdot (\textbf{u} \textbf{u})$，为非线性项。一些算法将其中一个速度作为已知量，另一个作为未知量离散求解。可参考icoFoam解析。

@Wayne

显性方法中要求 $\frac{\Gamma \Delta t}{\rho (\Delta x)^2}<0.5$，以及苦郎树小于1. 如果你的例子库朗数很低。那满足第一个条件即可。

参见：M. Peric，143页。

ｏｐｅｎｆｏａｍ　边界条件引入到离散控制方程求解中，实际是把边界节点的值引入到了源项来看待的吧？最想不明白的就是再靠近壁面的体积单元Ｐ，再其通过惯性力/压力差/剪切力形成动量方程的时候，是不是还要把壁面剪切应力考虑进来。

@ywj123450
《数值传热学（第2版）》是我唯一力荐的一本中文CFD理论书。