可以用newcommand

\newcommand\xx{\begin{equation}} 
\newcommand\xxx{\end{equation}} 

\xx
1+1=2
\xxx

输出：

是这种吗

回复自己。。。用bindata
https://www.openfoam.com/documentation/guides/latest/doc/guide-fos-forces-forces.html

最近在研究圆柱绕流。 用的算例是李老师提供的算例：http://www.dyfluid.com/icofoam.html

现在我把圆柱从二维（Z方向一层网格）拉伸到三维计算圆柱绕流。 但是现在的forceCoeffs 和Forces function只能生成整个三维圆柱的升力和阻力。 如果我想生成沿管长方向的升力分布，例如均分成100个分段的升力。 Openfoam能够见到做到这点吗

隔壁有人提供了这种后处理方法，但是看起来需要每个时间步单独处理，而且三维的数据量太大了，，这样好像不能每个时间步自动处理。

These are my paraviewsteps in short:

1)Extracted Surface
2)Generate Surface Normals
3)Calculated:
cp = normals*(p2-p1); ->Pressure Drag
cf=tangentialsmuwallGradU_y -> Friction Drag (or calculated directly by foamCalc components wallGradU
4) Build Integral Drag=S(cp_X+cf_X)dA2/rhovvA

最近也在研究这个函数。 看了半天没有明白。看了一下2.X的pisofoam中用的是turbulence->divDevReff(U). 这个比较好找。。 请问一下老师这两个可以当做是一个东西吗

@五好青年 tutorials里面的wingMotion

@dzw05 您好~ 那我可以理解成图上的差距就是SST模型模化的雷诺应力和真实的雷诺应力的差距吗？

@cccrrryyy 您好~

我将URANS的结果缩比成与Quasi-DNS相同范围则不能显示了。

关于收敛的问题。使用的都是非定常计算，fvSolution里面的值都是默认的。时间步长收敛性上，库朗数基本上在0.5左右。网格收敛性都经过验证了，而且用的网格都非常密了，这样一个流场用了15W个网格计算的。

所以这个结果是不合理的是吗。。。

算圆柱绕流，利用icofoam分别不用湍流模型(laminar)和SSTkomega湍流模型算出了两种工况

层流模型利用 fieldAverage生成雷诺应力

 fieldAverage
    {
        type            fieldAverage;
        functionObjectLibs ("libfieldFunctionObjects.so");
        enabled         true;
        outputControl   outputTime;
        fields
        (
            U
            {
                mean        on;
                prime2Mean  on;
                base        time;
            }
        );
    }

sstkomega模型利用Postprocss -func R生成雷诺应力

但是两种方法生成的雷诺应力完全不同，quasiDNS的雷诺应力XX场呈现下图且随时间基本不变

SSTkomega的雷诺应力XX场为

且像涡街一样随着时间而变化。

两者算出来的雷诺应力差别有这么大吗？ 还是我哪里的计算出了问题？

gradU输出格式是