CFD有限体积3-网格

最新推荐文章于 2023-09-24 00:02:53 发布

觅°

最新推荐文章于 2023-09-24 00:02:53 发布

阅读量644

点赞数

分类专栏： CFD 笔记文章标签：算法

本文链接：https://blog.csdn.net/qq_41817886/article/details/128876954

版权

笔记同时被 2 个专栏收录

11 篇文章 1 订阅

订阅专栏

CFD

8 篇文章 0 订阅

订阅专栏

参考文献：

F. Moukalled, L. Mangani, M. Darwish. The finite volume method in computational fluid dynamics. An advanced introduction with Openfoam and Matlab[M]. Chapter 6.

本文仅介绍结构化网格。

1. 有限体积网格

以“在网格上计算参数梯度”问题说明FVM对网格的要求。

网格单元上参数C的平均梯度为

根据散度定理，将体积分转化为面积分

根据中值定理

上式中 $\bar{\phi _{f}}$ 为界面中心点处参数值，该值由与该界面相邻的网格单元中心点（cell-cenered）参数值确定，如下图所示，f处值由C处和F1处参数值插值获得。

因此，FVM需要网格生成模块提供的信息有：界面向量（ $\vec{S_{f}}$ ，包括面积大小和方向）、单元中心点和相邻单元中心点参数值（ $\phi _{C}, \phi_{F}$ ），以及计算界面上参数值（ $\bar{\phi_{f}}$ ）的插值方法（也可称为参数变化型线）。假设

$\bar{\phi_{f}}=g_{C}\phi_{C}+g_{F_{1}}\phi_{F_{1}}$

gc和gf1的给定方法代表不同插值方法，给定一种方法

$g_{F_{1}}=V_{C}/(V_{C}+V_{F_{1}}), g_{C}=V_{F_{1}}/(V_{C}+V_{F_{1}})$

2. 结构化网格

2.1 网格拓扑

网格不仅需要提供网格几何坐标值、界面参数等，还需要提供网格拓扑信息。网格拓扑信息代表了网格局部索引与全局索引之间的关系（需要网格拓扑是因为问题的全局求解矩阵是根据各网格单元的积分方程构造的）。对结构化网格，其拓扑是已知的，无需另外定义。

对二维结构化网格，全局坐标n与局部坐标(i,j)之间的关系为

对三维结构化网格，全局坐标n与局部坐标（i, j, k）之间的关系为

2.2 几何信息

对网格单元（i, j)

其四周界面为S1(i, j), S2(i, j), S1(i+1, j), S2(i, j+1) （为界面名，界面方向由全局索引较小单元指向全局索引较大单元）。当针对单元时，规定单元界面法向量向外，因此

存储信息：2D和3D结构化网格的相关信息存储可以用1D列表/矩阵，如对Ni*Nj的2D网格，可采用Ni*Nj长度的1D列表存储信息，单元（i, j）所存储列表位置为其对应的全局索引n。

2.3 单元参数

对单元（i, j）

其梯度根据公式

可写作

$\bar{\nabla\phi_{i,j}}=\frac{1}{V_{i,j}}(\bar{\phi}_{i+1/2, j}\vec{S}_{i+1/2,j}+\bar{\phi}_{i-1/2, j}\vec{S}_{i-1/2,j}+\bar{\phi}_{i, j+1/2}\vec{S}_{i,j+1/2}+\bar{\phi}_{i, j-1/2}\vec{S}_{i,j-1/2})$

根据界面方向，上式

$=\frac{1}{V_{i,j}}(\bar{\phi}_{i+1/2, j}\vec{S1}_{i+1,j}-\bar{\phi}_{i-1/2, j}\vec{S1}_{i,j}+\bar{\phi}_{i, j+1/2}\vec{S2}_{i,j+1}-\bar{\phi}_{i, j-1/2}\vec{S2}_{i,j})$

当采用全局索引时，上式可写作

$\bar{\nabla\phi_{n}}=\frac{1}{V_{n}}(\bar{\phi}_{n+1/2}\vec{S}_{n+1/2}+\bar{\phi}_{n-1/2}\vec{S}_{n-1/2}+\bar{\phi}_{n+Ni/2}\vec{S}_{n+Ni/2}+\bar{\phi}_{n-Ni/2}\vec{S}_{n-Ni/2})$