如何使用UDF修改LES模型中的滤波尺度？

[西湖冷月]

定义了一个新的LES亚格子应力模型，其中有部分改进涉及到滤波尺度，主要的：在原本的Smagorinsky模型中滤波尺度（网格在x y z方向间距乘积的1/3次方）的基础上，修改为（1）计算单位涡矢量与网格单元每一个顶点与网格体中心点之间的矢量的叉乘；（2）查询并确定上一步中得到的叉乘矢量之间相减的模的最大值；（3）对上一步中得到的最大值分别与网格在x y z方向的间距进行加权平均（有表达式），引入局部速度梯度表达滤波尺度；
想请教各位前辈，在fluent中用UDF进行修改的话，这三步分别需要哪些宏命令，UDF基本思路大概是什么，实施的时候需要对网格节点编码有什么特殊要求吗（例如，是否需要用到查询结构化网格节点编码结果的操作，fluent内部网格节点编码应该是非结构化的）？
谢谢大家，恳请高手不吝赐教！！

[李东岳]

虽然在OpenFOAM里面实现很简单，但是我不太熟悉Fluent爱莫能助了

问一下，这是你提出的滤波尺度？还是？有相关文献没有？
感觉可以发JFM啊。
小伙子我看好你

[西湖冷月]

@东岳 别人提了一部分，自己改了一部分，推导还没完全搞定，大致思路是有了，就是万水千山走不动啊

[hhj0219]

有一篇，中国农业大学的博士论文《离心泵三维流动大涡模拟方法研究》就是通过UDF修改SGS模型，你可以看看，不知道对你有帮助没。
同时请教一个问题，在分析LES结果分析时，往往需要统计长时间（假如40s）特征面上的物理量，而时间步长又特别小，每个时间步保存一个.dat文件多半不现实。。。有其他什么办法吗？

[西湖冷月]

@hhj0219 在开始计算之前，在fluent的计算界面下面有时均化统计的选项，可以设定想要的时均统计方法和结果（不同于自动保存的设定方案），计算结束后在后处理的时候就会有时均量，不知道这样行不行？

[hhj0219]

@西湖冷月 我不仅仅是为了统计时均值。主要目的是：提取瞬态流场，用于后续相干结构的研究，如拟序结构的检测，涡的识别，模态分析等。。

[hhj0219]

@西湖冷月 而且还存在一个问题，fluent提供的时均统计，好像是基于迭代步（sampling interval），每个时间步最初的迭代，残差肯定相对较大，这样全部统计，是否不合理。。也不知道我理解是否有问题，还望共同探讨一下。

[西湖冷月]

@hhj0219 这个sampling interval是采样间隔，在ANSYS Fluent Users Guide中的LES-Specific Solution Strategies部分有介绍，采用也是按照时间步来的。
另外，如果要是分析相干结构，貌似没啥好方法能减少数据存储量，我见过一个算燃气轮机用LES的，总时长3s，数据3个TB多，就这还嫌算的糙。确实如果哪位同仁有经验心得分享分享

[hhj0219]

@西湖冷月 :谢谢指正。我去翻来看看。以前在用户手册里，一直没有看见。