首先感谢李老师的热情帮助~~
这个PPT中的第9页给出公式（1.5.9）。如果这个公式是正确的，则我的疑问就没有了，因为我原图中（4）式的右端项，第二项与第三项抵消了。但可疑的，这个公式（1.5.9）似乎并不正确。贴上来再麻烦李老师帮我看看：
0_1502119393399_7fd74579-5b34-4c15-ac3c-1cc3aae56236-image.png

0_1501229627502_figure1.PNG

Von Neumann Stability Analysis 在普通的CFD教材里面讲的太简单了
但是到复杂非线性问题（比如Euler方程）中分析时间推进格式的阻尼特性 (damping characteritics) 的时候就比较复杂
而且一般期刊文章中只给应用结果而且缺乏细节
它背后的逻辑是什么，我一直没有搞清楚。。。
希望有大神能推荐一些详细一点的教材，能把思路讲的比较透彻的

@Kingkong 搜如下教材

高性能计算并行编程技术_MPI并行程序设计_[都志辉]：非常适合入门
Using MPI ：结合上一个教材一起看
Iterative Methods for Sparse Linear Systems [ Yousef Saad ] ：介绍 deomain decomposition 的数学理论

@yuan_neu 建议看看《张量分析》的相关课程
二阶 【各向同性】的 【张量函数】 （ Dij(r) ）
只可能是上述 (6.25) 的表达形式

所以，你首先要明白“张量”，其次再搞清楚“张量函数”的定义，然后需要理解“各向同性”的概念，最后才能明白 6.25 式是怎么推导出来的

实际上，NS 方程的 Newton 本构关系，也是 “可以” 基于 【线性】和【各向同性】的假设直接写出来的 ！

推荐看看清华黄克智的《张量分析》教材

@李东岳 coupled solver比segregated solver能更快收敛，但在处理复杂几何形状物体时不太容易收敛[1]。

[1] http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2012/KlasJareteg/KlasJareteg-OF2012-Project-3-Version-2.0.pdf

@浪迹天大 是这样的，我最近也在读OpenFOAM中LES模型，发现Lund提到
The nearly universal approach is to simply write down the filtered Navier-Stokes equations together with an assumed model for the subgrid-scale stresses and then apply the desired spatial discretization to this “filtered” system. Although it is rarely mentioned, what one is doing by adopting this procedure is to imagine that the finite support of the computational mesh together with the low-pass characteristics of the discrete differentiating operators act as an effective filter. One then directly associates the computed velocity field with the filtered velocity. This procedure will be referred to as implicit filtering since an explicit filtering operation never appears in the solution procedure.
也就是说OpenFOAM程序除了动态模型外都没有对N-S方程进行显式的滤波操作，那么我们在方程植入中只要修改$\nu_{sgs}$就OK了，我在后面会做几个滤波函数的选择的对比，看看最后结果一样不一样

不是很清楚你说的那种涡脱落现象。不过用户指南中我记得有一章曾经表示，大体意思为：“如果用稳态求解器模拟比较强的非稳态问题会出现残差的震荡”。

能否附稳态求解涡街的文献？谢谢。

无意间发现这个帖子，里面提到了normal turbulent model,
https://www.cfd-online.com/Forums/openfoam/125473-low-reynolds-turbulence-models.html

0_1497945059409_2017-06-20 15-50-42屏幕截图.png

震荡已经非常非常小了，

能否提供一个关于气液两相流关于y+的文献？可能研究侧重点不同，目前我关注的领域对y+要求不是很严格。

一般情况下，多相流动湍流要依据连续相而定，因此我的回答是应该根据液体来算。在某些极端情况下，气体的湍流场可以认为和液相湍流场相同。

@赵一铭 润滑理论的雷诺方程就是NS方程简化了，不考虑厚度方向的压力变化。因为在润滑中动静环之间的液膜很薄，可以认为在厚度方向没有压力变化，是雷诺提出的，就叫雷诺方程了。不知道和那个雷诺是不是一个雷诺。

@赵一铭 目前想通过现有商业软件做一做···openfoam有点难，项目要求又很急····

@arctan 非常感谢！

@cfd-china link text!![alt text]0_1492654992588_QQ截图20170420102158.png 这是这个文献上说的在固壁上的条件；还有一篇文献直接在壁面设置为0了link text0_1492655090308_QQ截图20170420101901.png

@李东岳已经非常感谢了，我再继续努努力吧！:expressionless:

@赵一铭 整个程序我用的SIMPLE算法，所以连续性方程改成了压力修正方程来求解压力，不能直接求解出密度，怎么办？

谢谢，我刚把这个导入到程序中，但结果仍然不理想，解出得P'偏大，在修正速度时会导致速度变成负的，而且是进口速度得负几十倍，所以总感觉压力方程和速度方程不匹配，肿么破？？:crying:

不是很了解这个技术。

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