首页/文章/ 详情

LES 中的 commutation error

1年前浏览1033
大涡模拟不仅对网格的解析度有要求,而且对网格质量也有更高的要求。这里主要介绍一种由网格带来的误差,commutation error,并简单分析其对结果的影响。

在 LES 模拟的程序中,常用的一种过滤方法是隐式过滤,即,用当地网格的尺度作为该处的过滤尺度。这种方法在各种类型的网格上实现起来都相对简单。只是,对于稍微复杂一点的几何构体,生成的网格的尺度总是不可能一样,这也就意味着过滤尺度是变化的。如果相邻网格的尺度变化很大,这也将引起相邻网格的过滤尺度相差很大,这时就会带来严重的 commutation error。

Commutation error 产生的根本原因是,当相邻网格过滤尺度不一样时:

下面举例说明。
假设有一个场量,解析值为:

其导数为:

在如下图的网格中来进行过滤:

第一次,对点 P 和 N处进行过滤,过滤直径都是 ,使用 top hat 过滤函数,即

则:

类似地,  

另一方面,对梯度使用过滤,得


这说明,当相邻网格的过滤尺度一致时,commutation error 为零。但是,当过滤尺度不一致时,考虑右边的网格中心为N的情形,此时:

这时就产生了 commutation error。

所以,网格尺度不均匀现象严重的区域就会出现显著的 commutation error。边界层是容易出现网格尺度不一致的区域,在壁面上,严格来说,过滤尺度必须是零,否则将破坏无滑移条件。这种情况还可以通过 DES 这样的方法来处理。但是,核心区也可能出现相邻网格不一致的情形。局部网格加密是一种实用的节省计算量的方法,这种方法只需要在重要区域进行加密,不必全局加密。但是局部加密就导致了相邻网格尺度不一致,由此就带来了显著的 commutation error。所以,使用局部加密时,要注意不要让相邻两种等级的网格的边界落在重要区域

局部加密效应可以从下图看出:

当从细网格向粗网格过渡时,会出现亚格子粘度的突降;  
当从粗网格过渡到细网格时,亚格子粘度会突然上升。  
有一种办法是对过滤尺度进行光滑处理  

C△S的取值约为 1.5。

用这种方法可让局部加密带来的问题得到部分缓解。亚格子粘度在粗细网格交接附近的过渡也更光滑。

误差总是无可避免的,降低误差的技术也有很多。从“计算性价比”(即为了得到一定准确度的结果所耗费的计算量) 的角度看,没有哪个结果就一定是严格地错误或者正确,只是“计算性价比”不同。有时候为了降低计算量,在一些区域即便产生了严重的误差,也是可以接受的。关键是要对误差差生的原因有清楚的认识,尽量不要让很大的误差出现在重要的区域。

本篇图/文/公式完全取材整理自 Eugene de Villiers 的博士论文 “The Potential of Large Eddy Simulation for the Modeling of Wall Bounded Flows”,特此声明。  
来源:  
http://xiaopingqiu.github.io/2016/08/21/commutationError/  
来源:多相流在线
OpenFOAM多相流湍流核能电力材料控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2023-06-23
最近编辑:1年前
积鼎科技
联系我们13162025768
获赞 134粉丝 119文章 319课程 0
点赞
收藏
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈