CFD中文网

学

@yangom 无滑移边界条件，下面就是固壁面，无滑移就是壁面上速度为0，那下表面就需要从入口处开始发展发展出一个边界层。在这个边界层中，沿壁面法向向外，速度当然是从0逐渐变化成外流速度的
你直接看也看不出来有没有速度，画几条曲线出来。
1 靠近近壁面某个高度上，沿流动方向的一条线，画这条线上沿流动方向的速度
2 选几个位置，沿垂直于壁面方向的一条线，画这条线上沿流动方向的速度

学

小球这里没速度是和你速度边界条件的设置有关系，你应该是把下面那个边界设置成了noSlip。
需要更多的信息。

学

加源项的时候，这个源项的表达式好像有点问题，比如下面添加vectorSemiImplicitSource的时候源项是( (299.4600 0 0) 0)而不是(299.4600 0 0)。原因是加入的源项表达式为：
$$
S=S_C+S_Pϕ_P
$$
标量源项的代码也是这么写的：

        injectionRateSuSp
        {
            variable_name     (Sc Sp);
        }

momentumSource
{
   type vectorSemiImplicitSource;
   active on;
   

   vectorSemiImplicitSourceCoeffs
   {
      selectionMode all;
      //volumeMode      absolute; // specific
      volumeMode        specific;
      injectionRateSuSp
      {
         U           ( (299.4600 0 0) 0); //partial p / partial x
      }
   }
}

两个参考：
https://caefn.com/openfoam/fvoptions-acousticdampingsource
https://xiaopingqiu.github.io/2016/03/20/fvOptions2/

刚又找了一下之前的笔记，不是上面这个原因，因为你用的是vectorCodedSource。我之前也刚好写过scalarCodedSource，也发现量纲不对导致的问题。之前记录的笔记是：
有一个问题是感觉并不需要乘以cell的体积。官方给的文档里面是有乘以体积的，但这个网页就很迷惑，https://www.openfoam.com/documentation/guides/v2012/doc/guide-fvoptions-sources-coded.html ，上面的表达式是除以体积，下面的代码又是乘以体积了。找到一个人的注释 https://xiaopingqiu.github.io/2016/03/20/fvOptions2/

// fvMatrix<Type> 类中对“+=”操作符进行了重载，所以，eqn与Su的相加，相当于eqn+Su*mesh.V()，要不然eqn与Su的量纲不一致。
eqn += Su + fvm::SuSp(Sp, psi);

所以还是要乘以体积。

coolhhh

@Afr1yne 如果模型什么都没变，只是网格加密， u k e正常也是一样的。把fvSchemes换成最稳定的格式计算看看，再逐个改二阶格式试算，看哪个设置导致的发散

coolhhh

@Afr1yne 如果普通的均匀网格都发散，有点难搞，得再仔细查一查什么原因。fvSchemes和fvSolution改为更稳定的格式试试

coolhhh

基本不用修改，发散大概率是网格问题。看看fvSchemes和fvSolution怎么设置的，有的格式在好点的网格能算，复杂的网格就发散

coolhhh

对于不可压缩湍流，of输出的wallShearStress还要乘以空气密度。另外是否可能实验的数据无量纲化处理了？找下类似的研究，别人的文献数量级是多少，判断哪个计算错了

学

有一个在区域当中加源项以实现消波的模块
https://caefn.com/openfoam/fvoptions-acousticdampingsource

学

function object surfaces
实时输出就用

writeControl timeStep;
writeInterval 1;

coolhhh

我没具体操作过，正常画网格时候就要像你那个图下面挖空一样，额外挖空一个小长方形

CFD中文网

讲师

帖子

成员列表