1 概念介绍

文章二十一、文章二十二分别介绍了动网格的Layering方法和弹簧(Spring)光顺方法，这篇文章我们来介绍一下扩散(Diffusion)光顺。

2 基本原理

原理：扩散光顺方法与弹簧光顺原理上不同，它使用扩散方程进行网格的拉伸或压缩。如下图，其中γ为扩散系数，会对网格运动产生影响。

对于边界距离

对于网格体积

当不勾选remeshing选项时，不会发生网格重构，因此节点个数及连接关系并不会改变，只是改变了节点的相对位置。

使用条件：可适用于任意的网格类型；光顺方法一般会和remeshing方法共同使用。扩散光顺的方法相较于弹簧光顺适用于大变形。

3 模型描述

本例的模型与文章二十二相同，如图所示。模型有四个面，分别为底面(bottom)、上面(top)、左侧面(left side)和右侧面(right side)。其中top面保持不动，bottom面以-0.01m/s的速度沿z轴的负方向运动，内部网格随着底部边界运动而被压缩，但网格节点并没有合并或分离。

4 导入网格

使用Workbench打开工程文件，文件在本文末尾链接资源内。

5 Scale网格尺寸

Scale修改网格尺寸。如图所示。

确保计算域尺寸是我们所需要的。

6 设置求解器

由于网格运动是随时间进行的，因此动网格技术必须为瞬态模拟。将计算工况设置为瞬态。其他求解器设置可根据具体工况进行选择。

7 设置动网格模型

在模型树下双击Dynamic mesh，此时Task Page下页面如图

勾选Dynamic mesh，勾选Smoothing后，点击Settings，出现Smoothing方法的设置页面

点击Advanced，出现Mesh Smoothing Parameters参数设置选项。

Diffusion Function：

用于定于扩散函数，有两个选项Boundary distance和Cell volume，当选择某一选项时，在Diffusion Parameter输入参数，从而来定义扩散函数。其中Boundary distance和Cell volume表示的是边界距离和网格体积，其扩散函数如下。

对于边界距离

对于网格体积

两者区别主要是边界运动扩散到内部网格的方式不同。无特殊要求，可任选一种进行扩散。

Diffusion Parameter：

用于指定扩散参数，所指定的参数就是上述两个公式中的α，通过改变α的值可改变扩散系数γ，从而影响扩散。

当α取为0时，γ=1，表示均匀扩散。边界的运动能够均匀扩散到其他网格上。当α增大时，γ增大，使得扩散更远。一般α取值在0-3之间。

一般只需设置Diffusion Function和Diffusion Parameter，其他设置保持默认。

AMG Stabilization：

指定光顺网格的具体算法，保持默认即可

Maximum Number of Iterations：

最大迭代步数，默认值为50。保持默认即可

Relative Convergence Tolerance：

收敛相对残差，默认为10-10

Verbosity：

用于指定是否在控制台输出光顺方法的残差，1表示输出，0表示不输出。

Generalized Boundary Distance Method：

只在扩散函数选择边界距离时才可用，上述公式中的d为标准距离，当勾选此选项时，d发生改变，被称为“generalized” boundary distance。保持默认即可。

Smoothing From Reference Position：

从参考位置进行网格光顺，对于周期性网格效果较好

Dynamic Mesh Zones：

单击Dynamic Mesh Zones下的Create/Edit，弹出如下窗口

bottom边界：刚体边界；Zone Names选择bottom；Type选择Rigid Body；Motion Attributes下只需要定义Motion UDF，导入文章二十二的UDF，选中此UDF，其他均保持默认，单击Create，则在Dynamic Mesh Zones栏下出现bottom。

Right side：变形边界，Zone Names选择Right side；Type选择Deforming；Geometry Definition下Definition选择Faceted，定义变形边界；

Left side：与Right side保持相同的设置。

8 查看动网格

在Task Page页面单击Preview Mesh Motion，设置时间步数为100，时间步长为0.1s，显示频率为1。

Diffusion Parameter=0

Diffusion Parameter=3

以上案例的cas和dat文件均可以免费获取，需要的朋友，只需要点赞关注收藏一键三连后私信我即可哦

二十三、动网格Smoothing Diffusion方法及实例相关推荐

二十二、动网格Smoothing Spring方法及实例
1 概念介绍最近一直忙着其他的事情,好久都没有更新,让各位家人们久等了. 今天我们接着聊一聊动网格问题,我们在文章二十一中介绍了动网格的Layering方法,但是Layering方法的局限性很高,只 ...
【Python学习系列二十三】Scikit_Learn库降维方法(矩阵分解)-PCAFA
1主成分分析PCA 1.1 精确PCA和似然估计 PCA基于最大方差的正交变量分解多维数据集.在scikit-learn库中,PCA的实现是先通过fit方法计算n维的特征值和特征向量,然后通过tran ...
Helm 3 完整教程（二十三）：使用 Files 方法在模板中读取文件内容
推荐阅读 Helm3(K8S 资源对象管理工具)视频教程:https://edu.csdn.net/course/detail/32506 Helm3(K8S 资源对象管理工具)博客专栏:https: ...
二十四、动网格Remeshing方法及实例
1 概念介绍这次我们做一个比较复杂,同时也是比较实用的案例-网格重构.网格重构(Remeshing)方法是动网格常用的方法,一般会和弹簧光顺方法及扩散光顺方法相结合使用.当网格变形达到一定的程度后, ...
COMSOL三维动网格步骤！！！
鄙人经过两天的潜心研究终于学会了,啊啊啊!目前网上大多数资料都是二维动网格步骤,对于三维的甚少,我必须要给大家分享一下我的成果! 这个是一些关于二维的步骤,亲测可用! COMSOL中的动网格_midi ...
二十一、动网格Layering方法及实例
1 概念介绍我们之前给出的案例计算域边界都是固定不变的,如文章二十的多孔介质.文章十八的DPM模型等,但实际上一些工况下计算域的边界可能发生运动或者变形,如风扇扇叶的运动.齿轮的运动等. 对于计算 ...
OpenCV学习笔记（二十一）——绘图函数core OpenCV学习笔记（二十二）——粒子滤波跟踪方法 OpenCV学习笔记（二十三）——OpenCV的GUI之凤凰涅槃Qt OpenCV学习笔记（二十
OpenCV学习笔记(二十一)--绘图函数core 在图像中,我们经常想要在图像中做一些标识记号,这就需要绘图函数.OpenCV虽然没有太优秀的GUI,但在绘图方面还是做得很完整的.这里就介绍一下相关 ...
一个Fluent动网格问题及解决方法的记录
一个Fluent动网格问题及解决方法的记录 case中有固体和液体两种区域.由于固体模仿的是推进剂非均匀燃烧,因此燃面上下两侧边界层外有一部分使用了三角形网格.在计算过程中,推进剂的下部和右侧均是流场 ...
第二部分：二十三种设计模式解读——什么是工厂方法模式
二十三种设计模式解读--什么是工厂方法模式 author:陈镇坤27 日期:2022年2月10日修改日期:2022年6月23日文章目录二十三种设计模式解读--什么是工厂方法模式一.工厂方法模式 ...

二十三、动网格Smoothing Diffusion方法及实例