@东岳 成功解决了，在头文件里加了些东西，好有成就感 。准备复盘，写个教程
开心

今年九月份博士入学，研究生做了2年半的OpenFOAM，做的课题不难，在算法上对求解器也没有变动。博士阶段，自身非常想往CFD这方面靠一靠。感觉自己现在就属于那种在OpenFOAM门口站的那种人（不知道有没有人跟我有一样的感觉，就是OF的水平没多高，就是会调调参数，在控制方程里面做点小改动什么的，都能做到，但是再就不行了），十分想想听听老师前辈们的故事 目的就是能打破我目前的状态，正在看那本 The Finite Volume Method in Computational Fluid Dynamics,那本书反正已经拿起来放下好几次了，这回心里立了flag了，我就是啥也没看懂，我也一页一页把那本书翻一遍。

@Calf-Z-DNS 之前硕士的时候在做示踪剂在水中传输过程中，密度随浓度变化的问题，所处理的问题和你遇到的差不多，我的压力也是收敛的不太好，大概在0.01，然后硕士毕业时间紧迫问题就放下了，看见有人也在做这个问题，突然感觉我以前做的工作还挺有意义的（以前觉得没什么意义 ），现在突然想捡起来做一做了

OpenFOAM紧靠壁面网格结点上的物理量（比如湍动能）计算值怎么用代码提取出来？

@李东岳 是个稳态计算过程，不随时间变化，但是公式中r值更新会影响边界的U的值

在此更新一下，我知道问题出现在哪了，问题在于我计算的方程是输运方程，速度U是我给定的速度场，在计算的时候这个速度场的边界条件一直没更新。我现在想在求解输运方程的程序里更新一下速度的边界条件应该加点啥代码呢？
我查了几个，不知道对不对，大家有想法的捧个场，没想法的，捧个人场也行，嘻嘻嘻嘻

Uc.correctBoundaryConditions();
           fvOptions.correct(Uc);

或者上面这两行代码啥意思，大家有知道解释一下也行哈，谢谢~

@李东岳 就是我最后的计算结果文件中，速度那个文件里，使用这个边界条件的那个边界上没有速度，就还是下图这样的

我想象中是这样的：

我现在想在壁面处写一个速度的边界条件，具体表达式如下：

其中，k为湍动能，算好的值，被我放在0文件夹下；r为半径，不是算好的值，是目标算例计算过程中计算得到的值；n为表示速度的方向，就是壁面的法相方向向外。
我打算用codeFixedValue边界条件实现，然后下面是我写的代码，已经编译成功。

INLET_WALL
    {
        type            codedFixedValue; 
        value           uniform (0 0 0);
        name            Fdc1;
        code
        #{
          const fvPatch& boundaryPatch = patch();
            const fvBoundaryMesh& boundaryMesh = boundaryPatch.boundaryMesh();
            const fvMesh& mesh = boundaryMesh.mesh();
             scalarField r = mesh.lookupObject<volScalarField>("r");
            const volScalarField& k = mesh.lookupObject<volScalarField>("k");
            const vectorField Cf=patch().Cf();
  const vectorField a = patch().Sf()/mag(patch().Sf());
             vectorField Fdc1(Cf.size());
              forAll(Cf, faceI)
              {
                Fdc1[faceI]=(a[faceI]*0.01*0.3*4*k[faceI]*r[faceI])/(8e-7*pow(6,1.0/3));
               }
               operator==(Fdc1);
          #};
    }

问题是当我运行算例的时候，显示了使用了这个边界条件，但是计算结束后，我去找边界处速度的值，结果发现没有值，所以我感觉好像上面程序里面的运算好像没有进行一样，想让大家帮我看看我这个程序还有什么不对的地方吗？导致出现这种情况。

@李东岳 好像不是这个 又显示了mesh没有声明....

@李东岳 在 OpenFOAM小代码 中说：

//const fvBoundaryMesh& boundaryMesh = boundaryPatch.boundaryMesh();
//const fvMesh& mesh = boundaryMesh.mesh();

我想用codedFixedValue,写一个壁面的边界条件，壁面都是wall，而不是patch，所以这里想用这两个命令，但是编译过程中总是提示我是不是应该声明点啥呢？我也查了网上的，发现全是在inlet使用这个功能，然后入口是patch

想请教一下大家，我想用OpenFOAM算连续介质中固体颗粒的碰撞聚合问题，固体颗粒的量还比较大，我想问一下，就拉格朗日方法来说，用MPPICFoam可以实现嘛？我主要关心的点在于两个：
(1)MPPICFoam能不能算连续介质中固体颗粒的碰撞聚合问题；
(2)我看了一些我这个课题前人的文献，基本上都是减少了颗粒的数量，比如原本工业反应器中可能有上亿的颗粒，他们给简化到了几十万或者几万，MPPICFoam计算这种数量很多的颗粒问题，计算量是不是过于大了？如果过于大了，我也简化数量的话，计算的话是否有研究的意义。

@dingcy 自问自答了，在tecplot里输出非压缩的视频文件，然后放到ps里，增加分辨率，然后再渲染视频，最后的视频就清楚啦。

为啥我用tecplot输出的视频都好模糊啊，大佬们有方法吗？