CFD中文网

是不是我的问题没有问到点子上？请各位赐教。

麻烦各位了，刚开始接触openfoam，困难重重呀。

要求解的问题是个一维轴对称问题（无限长的均匀圆柱），请问如何制定坐标系，使div等算符使在圆柱坐标系下运算？谢谢各位

@李东岳 您好，抱歉，很久没来了。当时在学习FVM，最近在尝试利用openfoam解决NS方程。

@李东岳 明白了，感谢您的回复。后续要学习openfoam了~

@李东岳 非常感谢，我确实表达存在问题。对于瞬态计算，是用前一时间步的值确定边界处值的形式，再请教一下，请问对于稳态计算呢？如何处理？

@李东岳 您说的我不是很理解，如下图，以MUSCL格式为例，我理解限制器的函数形式可能是0,2rf，(rf+1)/2, 或2，具体要通过rf的值来判断，但rf的值取决于节点的具体取 值，节点是未知量（或者说是待求量，FVM最终要形成求解节点值得线性方程组），我的问题是如何确定限制器的函数形式（是2rf 还是 （rf+1）/2）?

@李东岳 感谢您的回复，对于你所提到的几种格式面上值与节点值无关我是理解的，但对TVD等High resolution 我就不理解了，比如TVD（如下图），觉得通量限制器与节点值有关呀？

@huangsx 各位，是否是由于我对这个概念理解错误仍然没有问到点子上。能否麻烦各位推荐一些high resolution scheme（包括NVF和TVD）的基础资料（书籍或其他入门资料）？感谢各位

@李东岳 李老师，感谢您的回复。由于是我对一些基础概念理解不深，造成了表达上的问题（这可能也是我困惑的主要原因），我重新组织了我的问题如下，请您多多帮助和指教。

在利用HR scheme对对流项进行离散时，控制单元边界处的值通过一个分段函数确定，函数形式根据上流控制单元中心值确定，在进行稳态计算时，所以控制单元中心值都是待求量，请问各位此时如何利用HR scheme，如何确定求每个控制单元边界值函数形式？对于时间域计算，是否是根据上一时刻迎风控制单元中心值确定求边界值得函数形式？

@李东岳 非常乐意分享，只是我刚接触流体力学和计算流体力学，理解可能并不到位。稍后我会把我的理解写下来，请大家批评指正！

昨晚通过学习《The finite method in computational fluid dynamics-an advanced introduction with open-foam and matlab》找到了问题的答案，现在的问题是怎么删除发表的帖子呢？

初学有限体积法，最早看的是李人宪的《有限体积法基础》，原因是这书简单，容易入门。其中对一阶迎风格式的数值扩散原因解释为扩散项仍采用中央差分格式所以造成假扩散。对此很是疑惑。因为

1. 几乎所有演示一阶迎风格式的数值扩散的算例都只有对流项而没有扩散项；

2. 扩散项采用中央差分格式从物理上看似乎更合理。

请问各位两个问题：

1.一阶迎风格式数值扩散的原因？

2.稳态扩散对流方程应用一阶迎风格式是否同意存在数值扩散？