以三维速度矢量为例
滑移边界条件的完整定义是：
A. 面法向的速度分量为零
B. 面切向的应力分量都为零
条件 A 控制 1 个自由度, 条件 B 控制 2 个自由度
A 和 B 一起控制速度矢量的 3 个自由度
对称边界条件的完整定义是:
先把域内速度场外插到边界面上，记为 Ub
然后边界面上的速度设为 Ub 沿边界面切向的分量
也就是说，无滑移边界条件需要指定边界面上切应力为零
而对称边界条件，是直接指定边界面上的速度矢量
一般来说，对于无粘流动，任意曲面的对称边界条件就等价于无滑移边界条件
但对于有粘流动，如果指定对称边界条件
那么边界面上依然有可能存在切向的应力
@cccrrryyy 我搞明白了，直接把收敛的速度场复制到 0 文件夹下，续算就行了；收敛场的结果相当于用 non-uniform field initialization ，里面也包含边界条件设置，可以直接放到 0 文件下当作另一个求解器的初始场；
@cccrrryyy 是在算RANS，网格也很规整，过渡也很平滑；似乎是 reactingFoam 直接算冷态超声速场比较容易发散；所以打算先用可压的求解器整一个初始场
我想试一下用 reactingFoam 算一个弱超声速场（Ma=2左右）的燃烧流动
dt=1.e-10s 还是容易发散，无法收敛
所以想试一下，先用 rhoCentralFoam 算一个收敛的冷态流场
然后把冷态流场的结果设为 reactingFoam 的初始场，然后续算
现在的问题是，rhoCentralFoam 没有组分场的信息
想请教各位大神：
如何把 rhoCentralFoam的结果（如：速度，压力，密度，湍流 k-espilon）
赋值给 reactingFoam 的初始场
而组分场（如：O2, N2, H2）还是按照 reactingFoam 本身的 0 文件夹的方式来初始化
比如运行 reactingFOAM 的时候
屏幕上不停的输出：
DILUPBiCG: Solving for h, Initial residual = 2.1859e-07, Final residual = 2.1859e-07, No Iterations 0
ExecutionTime = 8545.97 s ClockTime = 8756 s
Courant Number mean: 0.00199512 max: 0.401162
deltaT = 2.09015e-06
Time = 0.00879189
min/max(T) = 291, 1880
time step continuity errors : sum local = 9.41463e-08, global = 1.84917e-08, cumulative = 1.25902e-07
这些个没营养的信息
而且屏幕输出还占用大量时间
该怎么关闭呢？不让它们从屏幕上输出来
@星星星星晴 雀食祖传的，而且还得把 Fortran 的转成 OpenFOAM 。。。
@星星星星晴 啊我懂了~非常感谢！就是得用正儿八经的 case 跑 。。
@星星星星晴 你是说需要随便找一套网格，然后按照CFD模拟的设置，先做 decomposePar，这个 Test-parallel 才能正常运行吗？ 可这个程序不是跟网格一点关系都没呀。。。难道说因为版本太低了无法自动识别。。。OTZ
@星星星星晴 可以告诉我哪个教程的哪一页有讲吗？我没找到啊 :)
我想学习一下怎么用 OpenFOAM：applications/test/parallel 里面的并行函数
遇见如下几个问题：
1.直接编译运行显示需要读 system/controlDict 文件，但是 test 中没有给
2.随便找了个 controlDict 文件可以读，且可以运行 Test-parallel
3.mpirun -np 4 Test-parallel 运行之后，并没有真的并行，而是把 Test-parallel 独立运行了四次
4.mpirun -np 4 Test-parallel -parallel 运行报错：
[0]
[0]
[0] --> FOAM FATAL ERROR:
[0] Cannot read "/home/xq/OpenFOAM/OpenFOAM-2.3.1/applications/test/parallel/system/decomposeParDict"
[0]
FOAM parallel run exiting
[0]
[XQ:04726] 1 more process has sent help message help-btl-vader.txt / cma-permission-denied
[XQ:04726] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages
5.应该如何运行 Test-parallel ？实现 Test-paralle.C 中代码正常运行的结果？
先谢过各位大佬！
@李东岳 多谢！
我想在 OpenFOAM 中构建一个依托于网格格心的 整数 场
该怎么做呢？有没有现成的东西能用？
这个场在每个时间步都会重复更新并使用
@Hungryandfool
Operator 翻译过来叫做 ”算子“，也就是说，微分方程的空间项可以看作是不同的算子，比如 “对流”（散度），“扩散”（梯度）；
Splitting 的意义在于不同项可以相对独立的去求解，甚至采用完全不同的求解器，这样有利于提高计算效率，减小计算复杂性
$\partial T/\partial t = 1/r^2\partial[r^2\lambda(\partial T/\partial x)]/\partial x$
在原点 (r=0) 处，粘性边界条件应该怎么设置？
@lsprxd 是这么回事`
@程迪 追求的不是效果。。。是希望学习一门技术。。。。。
啊。。。。看到了这篇文献。。。它真的不能算是正儿八经的 【预处理】。。。。。
Von Neumann Stability Analysis 在普通的CFD教材里面讲的太简单了
但是到复杂非线性问题（比如Euler方程）中分析时间推进格式的阻尼特性 (damping characteritics) 的时候就比较复杂
而且一般期刊文章中只给应用结果而且缺乏细节
它背后的逻辑是什么，我一直没有搞清楚。。。
希望有大神能推荐一些详细一点的教材，能把思路讲的比较透彻的