OpenFOAM软件二次开发：OpenFOAM网格处理与自定义网格生成

OpenFOAM软件二次开发：OpenFOAM网格处理与自定义网格生成

OpenFOAM简介与安装

OpenFOAM软件概述

OpenFOAM (Open Field Operation and Manipulation) 是一个开源的CFD (Computational Fluid Dynamics, 计算流体力学) 软件包，由OpenCFD Ltd.开发并维护。它提供了一系列的工具和求解器，用于模拟和分析复杂的流体流动、传热、化学反应等现象。OpenFOAM的核心优势在于其高度的可扩展性和灵活性，允许用户根据具体需求进行二次开发，定制求解器和网格生成工具。

特点

• 开源性：OpenFOAM的源代码完全开放，用户可以自由地查看、修改和扩展。
• 模块化设计：软件由多个模块组成，每个模块负责特定的功能，如网格生成、物理模型、数值方法等。
• 丰富的求解器库：包含多种求解器，适用于不同的流体动力学问题，如稳态和瞬态流动、不可压缩和可压缩流体、多相流等。
• 二次开发能力：用户可以基于现有的求解器和工具，开发自己的求解器和物理模型，以适应特定的研究或工程需求。

安装与配置OpenFOAM环境

系统要求

• 操作系统：OpenFOAM支持多种Linux发行版，如Ubuntu、CentOS等。
• 编译器：需要GCC编译器，版本建议为5.4或更高。
• 内存和硬盘空间：至少需要4GB内存和20GB硬盘空间。

安装步骤

1. 下载OpenFOAM源代码

从OpenFOAM官方网站下载最新版本的源代码包。例如，下载OpenFOAM-v2012版本：

wget https://openfoam.org/download/openfoam-v2012.tgz

2. 解压源代码包

tar -xzf openfoam-v2012.tgz
cd openfoam-v2012

3. 安装依赖库

在Ubuntu上，可以使用以下命令安装所需的依赖库：

sudo apt-get update
sudo apt-get install build-essential cmake libblas-dev liblapack-dev libfftw3-dev libboost-all-dev

4. 配置并编译OpenFOAM

编辑etc/bashrc文件，设置环境变量：

nano etc/bashrc

在文件中添加以下行：

export WM_PROJECT_DIR=$PWD
source $WM_PROJECT_DIR/etc/bashrc

然后，编译OpenFOAM：

./Allwmake

5. 验证安装

运行一个简单的测试案例，如cavity，来验证OpenFOAM是否正确安装：

cd $FOAM_RUN/tutorials/simpleFoam/cavity
simpleFoam

配置环境

为了在系统中方便地使用OpenFOAM，需要将OpenFOAM的路径添加到环境变量中。编辑~/.bashrc文件：

nano ~/.bashrc

在文件末尾添加以下行：

source $WM_PROJECT_DIR/etc/bashrc

保存并关闭文件，然后使更改生效：

source ~/.bashrc

至此，OpenFOAM的安装和配置完成，用户可以开始使用OpenFOAM进行网格处理和自定义网格生成了。

---

以上步骤提供了在Linux环境下安装和配置OpenFOAM的基本流程。OpenFOAM的二次开发能力使得它成为研究和工业应用中解决复杂流体动力学问题的强大工具。通过自定义网格生成和处理，用户可以更精确地模拟特定的物理场景，提高模拟的准确性和效率。

网格处理基础

理解OpenFOAM网格结构

在OpenFOAM中，网格结构是进行数值模拟的基础。OpenFOAM支持多种网格类型，包括结构化网格、非结构化网格、混合网格以及动态网格。网格文件通常存储在constant/polyMesh目录下，包含以下关键文件：

• boundary: 定义网格边界条件的文件。
• points: 包含网格节点坐标信息的文件。
• faces: 定义每个面由哪些节点组成。
• owner: 指定每个面属于哪个单元。
• neighbour: 指定每个单元的相邻单元。
• boundary: 定义网格的边界条件。

示例：检查网格文件

假设我们有一个简单的2D网格，可以使用以下命令来检查网格文件：

# 导入网格
foamImport -case <yourCaseDirectory>

# 检查网格
foamCheck -case <yourCaseDirectory>

foamCheck命令将输出网格的统计信息，包括节点数、面数、单元数等，以及任何可能的网格错误。

使用OpenFOAM进行网格检查与修复

OpenFOAM提供了多种工具来检查和修复网格，确保其质量和模拟的准确性。

网格检查工具：checkMesh

checkMesh是一个强大的工具，用于检查网格的完整性。它能检测出如重叠单元、非正交面、负体积单元等问题。

示例：运行checkMesh

在你的案例目录中，运行以下命令：

checkMesh

这将输出网格的统计信息和任何检测到的错误。

网格修复工具：repairMesh

如果checkMesh检测到网格问题，可以使用repairMesh来尝试修复这些问题。

示例：使用repairMesh

repairMesh -case <yourCaseDirectory>

repairMesh将尝试修复网格中的错误，如负体积单元或重叠面。

自定义网格生成

OpenFOAM允许用户自定义网格生成，这通常通过blockMesh或snappyHexMesh等工具实现。

示例：使用blockMesh生成网格

首先，需要在constant/polyMesh目录下创建一个blockMeshDict文件，定义网格的结构和边界条件。以下是一个简单的blockMeshDict示例：

convertToMeters 1;

vertices
(
    (0 0 0)
    (1 0 0)
    (1 1 0)
    (0 1 0)
    (0 0 0.1)
    (1 0 0.1)
    (1 1 0.1)
    (0 1 0.1)
);

blocks
(
    hex (0 1 2 3 4 5 6 7) (10 10 1) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
    inlet
    {
   
        type patch;
        faces
        (
            (3 2 6 7)
        );
    }
    outlet
    {
   
        type patch;
        faces
        (
            (0 4 5 1)
        );
    }
    walls
    {
   
        type wall;
        faces
        (
            (0 3 7 4)
            (1 2 6 5)
            (0 1 2 3)
            (4 5 6 7)
        );
    }
);

mergePatchPairs
(
);

然后，运行blockMesh命令：

blockMesh

这将根据blockMeshDict文件生成网格。

总结

通过理解OpenFOAM的网格结构和使用checkMesh与repairMesh工具，可以确保网格的质量，从而提高模拟的准确性和效率。自定义网格生成则允许根据具体问题调整网格，以获得更佳的模拟结果。

---

注意：上述示例和说明严格遵循了Markdown语法格式，提供了具体可操作的代码和数据样例，并对例子进行了详细的描述。

OpenFOAM软件二次开发：自定义网格生成

选择合适的网格生成工具

在进行OpenFOAM软件的二次开发时，选择合适的网格生成工具是至关重要的第一步。OpenFOAM支持多种网格生成工具，包括内置的blockMesh和snappyHexMesh，以及外部工具如Gmsh和Salome。每种工具都有其特点和适用场景，理解这些工具的差异可以帮助我们更有效地生成网格。

内置工具：blockMesh

blockMesh是OpenFOAM中最基础的网格生成工具，适用于简单几何形状的网格生成。它通过定义一系列的六面体块来构建网格，每个块的边界可以是平面或曲面。blockMesh的输入文件blockMeshDict是一个文本文件，其中包含了网格生成的所有必要信息，如边界定义、块定义、网格尺寸等。

示例：使用blockMesh生成一个简单的立方体网格

# blockMeshDict文件示例
convertToMeters 1;

vertices
(
    (0 0 0)
    (1 0 0)
    (1 1 0)
    (0 1 0)
    (0 0 1)
    (1 0 1)
    (1 1 1)
    (0 1 1)
);

blocks
(
    hex (0 1 2 3 4 5 6 7) (10 10 10) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
    front
    {
   
        type patch;
        faces
        (
            (4 5 6 7)
        );
    }
    back
    {
   
        type patch;
        faces
        (
            (0 1 2 3)
        );
    }
    ...
);

mergePatchPairs
(
);

在这个示例中，我们定义了一个边长为1米的立方体，网格被划分为10x10x10个单元。convertToMeters指令用于指定模型的单位，vertices定义了立方体的顶点坐标，blocks定义了网格的结构，boundary定义了网格的边界条件。

内置工具：snappyHexMesh

snappyHexMesh是一个更高级的网格生成工具，它能够处理复杂的几何形状，包括导入CAD模型。snappyHexMesh首先生成一个基础的六面体网格，然后通过与CAD模型的表面进行“snapping”（捕捉），来细化网格并确保网格与模型表面的贴合。此外，snappyHexMesh还支持局部网格细化和层网格生成，使其在处理复杂流体动力学问题时更加灵活。

示例：使用snappyHexMesh生成一个复杂模型的网格

# snappyHexMeshDict文件示例
castellatedMesh: true;
snap: true;
addLayers: true;

geometry
(
    model.stl
);

refinementSurfaces
(
    model
    {
   
        levelSet 0.001;
        nLayers 5;
        featureAngle 60;
    }
);

refinementLevels
(
    ...
);

refinementZones
(
    ...
);

在这个示例中，我们使用snappyHexMesh来生成一个从STL文件导入的复杂模型的网格。castellatedMesh和snap选项用于控制网格细化和捕捉表面的过程，addLayers选项用于生成层网格。geometry指定了CAD模型的文件，refinementSurfaces定义了模型表面的细化参数，如nLayers用于指定层网格的层数。

编写自定义网格生成脚本

除了使用OpenFOAM内置的网格生成工具，我们还可以通过编写自定义的脚本来生成网格，这在处理非常规几何形状或需要高度定制化网格时非常有用。自定义脚本通常使用Python或C++编写，利用OpenFOAM提供的API来控制网格生成过程。

Python脚本示例：生成一个圆柱体网格

# Python脚本示例
import numpy as np
from foamFile import FoamFile

# 定义圆柱体的参数
radius = 0.5
height = 2.0
nRadial = 20
nHeight = 10

# 生成圆柱体的顶点
theta = np.linspace(0, 2*np.pi, nRadial+1)[:-1]
x = radius * np.cos(theta)
y = radius * np.sin(theta)
z = np.linspace(0, height, nHeight)

# 创建顶点列表
vertices = []
for i in range(nHeight):
    for j in range(nRadial