OpenFOAM:并行区域划分理解(Domain Decomposition)

OpenFOAM Domain Decomposition理解

OpenFOAM使用Domain Decomposition方法对计算域进行并行划分(命令: decomposePar), 对应的字典配置文件system/decomposeParDict。调用decomposePar后，会根据字典内的设定方法和方法对应的配置参数，按指定数目对计算区域进行划分，每个划分对应一个子目录，该子目录下的划分数据(polyMesh下)会记录本进程(processor)需要处理的(points,faces,cells)信息，以及和其他进程关联的划分边界信息(faces)，通过MPI启动多个processors并行运行求解器程序(加上 -parallel选项)

mpirun --hostfile <machines> -np <nProcs> <foamExec> <otherArgs> -parallel >log &
1

---

并行划分方法

---

源码位置

~/OpenFOAM-v2112/OpenFOAM-v2112/src/parallel/decompose
$ tree
.
├── Allwclean
├── Allwclean-mpi
├── Allwmake
├── AllwmakeLnInclude
├── Allwmake-mpi
├── decompose
│   ├── decompositionInformation.C
│   ├── decompositionInformation.H
│   ├── decompositionModel.C
│   ├── decompositionModel.H
│   ├── dimFieldDecomposer.C
│   ├── dimFieldDecomposerFields.C
│   ├── dimFieldDecomposer.H
│   ├── fvFieldDecomposer.C
│   ├── fvFieldDecomposerFields.C
│   ├── fvFieldDecomposer.H
│   ├── Make
│   │   ├── files
│   │   └── options
│   ├── pointFieldDecomposer.C
│   ├── pointFieldDecomposerFields.C
│   └── pointFieldDecomposer.H
├── decomposition.dox
├── decompositionMethods
│   ├── decompositionConstraints
│   │   ├── decompositionConstraint
│   │   │   ├── decompositionConstraint.C
│   │   │   └── decompositionConstraint.H
│   │   ├── geometric
│   │   │   ├── geometricConstraint.C
│   │   │   └── geometricConstraint.H
│   │   ├── preserveBaffles
│   │   │   ├── preserveBafflesConstraint.C
│   │   │   └── preserveBafflesConstraint.H
│   │   ├── preserveFaceZones
│   │   │   ├── preserveFaceZonesConstraint.C
│   │   │   └── preserveFaceZonesConstraint.H
│   │   ├── preservePatches
│   │   │   ├── preservePatchesConstraint.C
│   │   │   └── preservePatchesConstraint.H
│   │   ├── refinementHistory
│   │   │   ├── refinementHistoryConstraint.C
│   │   │   └── refinementHistoryConstraint.H
│   │   └── singleProcessorFaceSets
│   │       ├── singleProcessorFaceSetsConstraint.C
│   │       └── singleProcessorFaceSetsConstraint.H
│   ├── decompositionMethod
│   │   ├── decompositionMethod.C
│   │   └── decompositionMethod.H
│   ├── geomDecomp
│   │   ├── geomDecomp.C
│   │   └── geomDecomp.H
│   ├── hierarchGeomDecomp
│   │   ├── hierarchGeomDecomp.C
│   │   └── hierarchGeomDecomp.H
│   │   ├── Make
│   │   ├── files
│   │   └── options
│   ├── manualDecomp
│   │   ├── manualDecomp.C
│   │   └── manualDecomp.H
│   ├── metisLikeDecomp
│   │   ├── metisLikeDecomp.C
│   │   └── metisLikeDecomp.H
│   ├── multiLevelDecomp
│   │   ├── multiLevelDecomp.C
│   │   └── multiLevelDecomp.H
│   ├── noDecomp
│   │   ├── noDecomp.C
│   │   └── noDecomp.H
│   ├── randomDecomp
│   │   ├── randomDecomp.C
│   │   └── randomDecomp.H
│   ├── simpleGeomDecomp
│   │   ├── simpleGeomDecomp.C
│   │   └── simpleGeomDecomp.H
│   └── structuredDecomp
│       ├── structuredDecomp.C
│       └── structuredDecomp.H
├── faDecompose
│   ├── faFieldDecomposer.C
│   ├── faFieldDecomposerFields.C
│   ├── faFieldDecomposer.H
│   ├── faFieldDecomposerReadFields.C
│   ├── faMeshDecomposition.C
│   ├── faMeshDecomposition.H
│   └── Make
│       ├── files
│       └── options
├── kahipDecomp
│   ├── kahipDecomp.C
│   ├── kahipDecomp.H
│   └── Make
│       ├── files
│       └── options
├── metisDecomp
│   ├── lnInclude
│   ├── Make
│   │   ├── files
│   │   └── options
│   ├── metisDecomp.C
│   └── metisDecomp.H
├── ptscotchDecomp
│   ├── Make
│   │   ├── files
│   │   └── options
│   ├── ptscotchDecomp.C
│   └── ptscotchDecomp.H
└── scotchDecomp
    ├── Make
    │   ├── files
    │   └── options
    ├── scotchDecomp.C
    └── scotchDecomp.H

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118