CFD中文网

@李东岳 是单塔式电脑

可以用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输出格式是

我是用的这个

https://blog.csdn.net/liuxiabing150/article/details/51920182/

没错吧。。

《Fluid Mechanics》 Frank M.White P239

《工程流体力学》 李玉柱

弹性力也是这么定义的

@bestucan

最后找到了解决办法，是边界条件的原因。本来上下圆柱连接的上下两个边界用的symmetry边界条件，改成cylic边界之后就能算了。。。具体原因还没研究出来

@fireztw

https://openfoamwiki.net/index.php/Parameter_Definitions_-_dynamicMotionSolverFvMesh

https://forum.cfdac.com/t/topic/660/2

这两个网站可以参考一下~~

@赵海盛 看起来没有问题的。。如果是自激振荡就是caculated，如果是强迫就要oscillatingdisplacement了。。。看边界条件是一个三维的圆形流场吗？。。

二维还是三维呢。。边界条件问题吧一般是

CPU型号： Intel(R) Xeon(R) Gold 6248R
系统：ubuntu 18.04
版本：openfoam-8

48 86.08
32 93.05
16 125.13
8 153.34
4 287.14
2 524.41
1 1013.26

正常了。。。

CPU型号： Intel(R) Xeon(R) Gold 6248R
系统：ubuntu 18.04
版本：openfoam-8

48 98.12
32 96.14
16 128.01
8 155.18
4 304.95
2 542.82
1 1068.49

没关超线程，物理核48 感觉不合理啊哈哈 明天再弄一次