paraView-求展向平均问题

小考拉

各位前辈好！

我现在在计算三维水翼空化绕流，下图呈现的是alpha.water体积分数云图，现在想将三维流场进行展向平均（z方向），得到某一时刻的z方向平均流场，请问在paraView里面怎样处理呢？

哪位前辈可以指点一下思路呢？谢谢

benqing

@小考拉 paraview能不能做不清楚，但是可以写程序实现，下面是cfdonline上的程序，我用过，可以实现，你可以研究一下，源程序好像只有速度压力这些东西，alpha.water的话需要自己添加进去
https://www.cfd-online.com/Forums/openfoam-post-processing/106966-spanwise-averaging-three-dimensional-field.html

小考拉

谢谢前辈指点，我去试一试

小考拉

@benqing 在 paraView-求展向平均问题中说：

@小考拉 paraview能不能做不清楚，但是可以写程序实现，下面是cfdonline上的程序，我用过，可以实现，你可以研究一下，源程序好像只有速度压力这些东西，alpha.water的话需要自己添加进去
https://www.cfd-online.com/Forums/openfoam-post-processing/106966-spanwise-averaging-three-dimensional-field.html

前辈您好！

我最近在研究您回复我的cfd online上的求解展向时均的代码，cfd online上面有两个代码，分别如下和，请问前辈您之前用的是哪一套代码呢？

我现在在研究第一套代码的spanwiseAvg.C文件，但里面的一些遍历网格单元的代码自己不是很清楚，，请问前辈有哪些资料可以参考学习一下呢？

另外在第一套代码中（spanwiseAvg.C）的file文件中为什么是EXE = $(FOAM_USER_APPBIN)/spanwiseAvg呢？感觉看起来编译的是个求解器，

自己的of代码基础比较薄弱，还望前辈能够赐教，感激不尽

benqing

@小考拉我用的第二个
遍历网格单元的代码不知道你指的是什么意思
EXE = $(FOAM_USER_APPBIN)/spanwiseAvg这个就是编译一个程序，可以直接在终端调用，确实跟求解器的编译方法是一样的，其实求解器也是一个程序，类似于win里面的.exe的应用程序，只不过在这里需要通过终端调用

Mania

@小考拉您好，请问您成功编译这个代码了么，我现在也要做展向平均，有位前辈给我推荐了这个程序，但是我不知道怎么编译😂您能指教一下么

浪迹天大

推荐一个可用的代码，支持 OF-4 和 OF-7
https://github.com/ZmengXu/sampledPlaneAverage

小考拉

@Mania 我记得当时前辈推荐的这个代码里面有一个说明，直接放到有关的目录下面，编译链接即可。将postaverage文件加入到constant文件中，然后直接在终端输入命令 postaverage就可以生成对应的展向平均文件了。。还是比较简单的，建议你自己多试一试

vbcwl

@浪迹天大请问下这个怎么使用啊

coolhhh

@小考拉前面大佬提到cfd-online的程序，我下载第二个程序运行了下，源程序不做任何修改不会自动输出文件。我基于cfd-online第二个程序修改了一个在of8运行的程序，但有些疑问：

1. 修改的程序

程序下载链接：postSpanwiseAverage_of8.tar.gz
主要修改collapse.H文件，将internalField修改为primitiveFieldRef

collapse.H

  
    scalarField spanavgUMeanX
    (
        meshIndexing.collapse(UMean.component(vector::X)())
    );
    
    // avgUx.internalField() = spanavgUMeanX;
    // of8
    avgUx.primitiveFieldRef() = spanavgUMeanX;
    avgUx.write();

2. 问题1

虽然在postSpanwiseAverage.C增加了runTime.write()，但不知为何运行不会自动输出文件。需要在collapse.H增加avgUx.write()才可输出结果。这是什么原因造成的？

postSpanwiseAverage.C

     forAll(timeDirs, timeI)
    {
        runTime.setTime(timeDirs[timeI], timeI);
        Info<< "Collapsing fields for time " << runTime.timeName() << endl;
 
        #include "readFields.H"
        //#include "calculateFields.H"
 
        // Average fields over spanwise direction
        #include "collapse.H"
 
        runTime.write();
 
        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
            << "  ClockTime = " << runTime.elapsedClockTime() << " s"
            << nl << endl;
    }

3. 问题2

根据源程序avgUx.primitiveFieldRef() = spanavgUMeanX;，可以看出meshIndexing.collapse()计算的是内部场的结果，而没有计算边界场的结果（如下图速度场所示，入口的速度为0，没有求平均）。现在想对壁面的剪切应力求展向平均，但该程序无法计算，该程序遍历网格单元的代码没太看懂，想问各位大佬有无比较好的方法方便计算？

pengdi

您好您后来展向平均怎么做的呢


	scalarField spanavgUMeanX
	(
	meshIndexing.collapse(UMean.component(vector::X)())
	);

	// avgUx.internalField() = spanavgUMeanX;
	// of8
	avgUx.primitiveFieldRef() = spanavgUMeanX;
	avgUx.write();

	forAll(timeDirs, timeI)
	{
	runTime.setTime(timeDirs[timeI], timeI);
	Info<< "Collapsing fields for time " << runTime.timeName() << endl;

	#include "readFields.H"
	//#include "calculateFields.H"

	// Average fields over spanwise direction
	#include "collapse.H"

	runTime.write();

	Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
	<< " ClockTime = " << runTime.elapsedClockTime() << " s"
	<< nl << endl;
	}

CFD中文网

paraView-求展向平均问题

1. 修改的程序

2. 问题1

3. 问题2