@xpqiu 十分感谢！我会仔细根据您提的意见学习的，谢谢！

谢谢...

感谢管理员，

这个minimal example 不写文件确实是因为我忘了加上runTime.write():cold_sweat: 。它并不能重现我原来solver里面的问题。我再找找原因。

@李东岳 谢谢版主:laughing:

感谢，我看了下，UserGuide里面U28说的是把scale关掉，这会导致所有的箭头长度和箭头头部大小都一样。目前从ParaView的控制面板上还找不到保持箭头头部大小不变，而箭头长度随速度大小scale的选项。

一般的工具（Matlab, matplotlib, Paraview）里面， quiver图 里面的箭头大小都跟着箭头长度一起变，所以在矢量很小时候，效果很糟糕。如下图就是Matlab里面quiver画出来的，
0_1461126738148_upload-25983765-853b-4a0a-b72b-a0ffce9f781b

不过后来找到一个比较好的办法，在latex里面使用pgfplot宏包可以解决。
比如如下代码(来自这里)

\documentclass{standalone} \usepackage{tikz} \usepackage{pgfplots} \begin{document} \begin{tikzpicture} \begin{axis}[domain=-3:3, view={0}{90}] \addplot3[blue, quiver={u={1}, v={(x-y)}, scale arrows=0.15}, -stealth,samples=20] {0}; \end{axis} \end{tikzpicture} \end{document}

可以产生画出：
alt text
可以看到箭头大小不随箭头长度改变。

@李东岳 这一项和interpolationSchemes改为linear就可以运行，用MUSCL和venLeer就不可以

@聆听whatever

Hi,

可以看看这个：链接文本。欢迎补充交流。:baby:

我的并行出问题了，不并行倒是可以，我想具体学一学怎么用Debug模式来调试，能否举个例子教一下我呢？谢谢各位了~

@lhzhu 在 如何拷贝/创建 volVectorField 在每一个时间步文件夹中 中说：

只要0目录里面有一个场文件，decomposePar 都会在processorX里面生成子区域的场文件。是不是你写的0/HAG_BC 的file header 里面的object 不是 HAG_BC ?

decomposePar 后 在“processorX” 中的 "0" 目录里有这个HAG_BC, 问题是下一时间步的时候没法创建这个HAG_BC啊，
所以会有那个error.
不知道怎么解决

@李东岳 在 有关侵入边界法。 中说：

Hi，

OpenFOAM-extend里面有IBM方法，Extend和官方的OpenFOAM侧重不同，几个我了解到的不同主要是IBM，粘弹流体，动网格表面跟踪。

还有可压缩流动高分辨率格式。https://github.com/Unofficial-Extend-Project-Mirror/Unofficial-Extend-Project-Mirror.github.io/wiki

http://www.tfd.chalmers.se/~hani/kurser/OS_CFD/HrvojeJasak/DiscretisationBestPractice.pdf

看了这个，有点懂了，问题应该在于离散格式和梯度限制器啥的。
Finite Volume Discretisation in OpenFOAM Best Practice Guidelines

Hi,

并没有用过这个，不过cfd-online有关waves2Foam的讨论已经多达2000贴。从常识来看，附加湍流波高减小是可以理解的。

This questions is not clear. What is 液力变矩器?

@搬运工不好当 说：

求速度方程之后为什么没有求压力方程，而

Hi，

我觉得能量方程放在压力方程之后和之前结果是相同的。就像湍流求解，标量传输放在压力方程之前和压力方程之后也是相同的。

你可以自己测试一下并反馈？

@李东岳 谢谢东岳大神。祝论坛越来越好。

@李东岳 东岳兄，这个是一个基于OpenFOAM和LAMMPS的求解器，有一篇文章叫：Sun_2016_SediFoam A general-purpose, open-source CFD-DEM solver for particle-laden，可以一看。这个求解器其实功能是和CFDEM差不多，也就是求解密相流动的。其不属于CFDEM一部分。

Hi，

这和你的求解方程有关，例如：

fvm::ddt(rho) + fvc::div(U)
和
fvm::laplacian(nu, p) = fvc::div(phi)

是对称矩阵。更进一步的，第一个是对角阵。

fvm::ddt(U) + fvm::div(U)

则是非对称矩阵。

of用的是自己的容器

@xiaohsu2006 是根据温度和压力算出来的，　constant 参数后面的值是根据什么给的呢？