HYPRE（High Performance Preconditioners）是由劳伦斯利弗莫尔国家实验室（LLNL）开发的高性能并行线性代数求解器库，专为大规模科学计算设计，支持多种网格结构和求解器。其中，SStruct（Semi-Structured）接口是HYPRE中用于处理半结构化网格问题的核心模块之一。

---

一、SStruct 接口简介

1.1 什么是 SStruct？

SStruct 是 HYPRE 中介于 Struct（纯结构化网格） 和 IJ（完全非结构化） 之间的接口层，用于处理由多个结构化网格块（称为“parts”）组成的问题，这些块之间可以有不同的拓扑、网格尺寸、变量类型，甚至可以重叠或通过接口耦合。

SStruct 的设计目标是：

• 支持多部件（multi-part）网格系统；
• 允许不同部件之间通过“接口”（interface）耦合；
• 每个部件内部保持结构化网格的高效性；
• 支持并行分布式内存计算（MPI）；
• 可灵活处理复杂几何或多物理场耦合问题。

1.2 核心概念

• Part（部件）：每个 part 是一个独立的结构化网格块（如 StructGrid），可以有不同的维度、变量数、网格尺寸。
• Variable（变量）：每个网格点上可以定义多个变量（如压力、速度分量等），支持节点型、边型、面型、单元中心型等。
• Stencil（模板）：定义每个变量与其邻居的连接关系（类似有限差分模板）。
• Graph（图）：描述整个系统的稀疏矩阵结构，包括部件内部和部件之间的耦合。
• Matrix & Vector：基于 Graph 构建的矩阵和向量对象，用于求解 Ax=b。
• Interface coupling（接口耦合）：不同 parts 之间通过“接口”建立连接，例如共享边界点或插值关系。

---

二、SStruct 使用流程（简化版）

以下是使用 SStruct 接口的基本步骤：

#include "HYPRE_sstruct_mv.h"
#include "HYPRE_sstruct_ls.h"

// 1. 创建网格对象
HYPRE_SStructGrid grid;
HYPRE_SStructGridCreate(MPI_COMM_WORLD, ndim, nparts, &grid);

// 2. 设置每个 part 的变量类型和网格范围
HYPRE_SStructGridSetExtents(grid, part, ilower, iupper);
HYPRE_SStructVariable vars[1] = {HYPRE_SSTRUCT_VARIABLE_CELL};
HYPRE_SStructGridSetVariables(grid, part, nvars, vars);

// 3. 声明部件之间的接口耦合（可选）
HYPRE_SStructGridSetNeighborBox(...); // 用于周期性或重叠网格
HYPRE_SStructGridAddVariables(...);   // 添加接口点上的变量

// 4. 完成网格设置
HYPRE_SStructGridAssemble(grid);

// 5. 创建 stencil（模板）
HYPRE_SStructStencil stencil;
HYPRE_SStructStencilCreate(ndim, size, &stencil);
HYPRE_SStructStencilSetEntry(stencil, entry, offset, var);

// 6. 创建图（Graph）
HYPRE_SStructGraph graph;
HYPRE_SStructGraphCreate(MPI_COMM_WORLD, grid, &graph);
HYPRE_SStructGraphSetStencil(graph, part, var, stencil);

// 7. 添加接口耦合（跨部件连接）
HYPRE_SStructGraphAddEntries(graph, part, index, var, ...);

// 8. 完成图设置
HYPRE_SStructGraphAssemble(graph);

// 9. 创建矩阵和向量
HYPRE_SStructMatrix A;
HYPRE_SStructVector b, x;
HYPRE_SStructMatrixCreate(MPI_COMM_WORLD, graph, &A);
HYPRE_SStructVectorCreate(MPI_COMM_WORLD, grid, &b);
HYPRE_SStructVectorCreate(MPI_COMM_WORLD, grid, &x);

HYPRE_SStructMatrixInitialize(A);
HYPRE_SStructVectorInitialize(b);
HYPRE_SStructVectorInitialize(x);

// 10. 设置矩阵系数
HYPRE_SStructMatrixSetBoxValues(A, part, ilower, iupper, var, stencil, values);
HYPRE_SStructMatrixSetValues(A, part, index, var, nentries, entries, values); // 点对点
HYPRE_SStructMatrixAddEntries(A, ...); // 用于接口耦合项

HYPRE_SStructMatrixAssemble(A);
HYPRE_SStructVectorAssemble(b);
HYPRE_SStructVectorAssemble(x);

// 11. 设置求解器（如 PCG + PFMG）
HYPRE_SStructPCGCreate(MPI_COMM_WORLD, &solver);
HYPRE_SStructPCGSetTol(solver, 1e-6);
HYPRE_SStructPCGSetup(solver, A, b, x);

// 12. 求解
HYPRE_SStructPCGSolve(solver, A, b, x);

// 13. 清理
HYPRE_SStructPCGDestroy(solver);
HYPRE_SStructMatrixDestroy(A);
HYPRE_SStructVectorDestroy(b);
HYPRE_SStructVectorDestroy(x);
HYPRE_SStructGraphDestroy(graph);
HYPRE_SStructStencilDestroy(stencil);
HYPRE_SStructGridDestroy(grid);

---

三、应用场景举例

3.1 多块结构化网格（Multi-block Structured Grids）

场景：CFD 模拟中，复杂几何被划分为多个结构化块（如机翼、机身、尾翼），每个块独立网格化，块之间通过接口匹配。

SStruct 优势：

• 每个块内部使用 Struct 高效存储；
• 块间接口通过 HYPRE_SStructGraphAddEntries 或 SetNeighborBox 耦合；
• 支持非匹配网格接口（需用户手动插值）。

---

3.2 多物理场耦合（如流固耦合 FSI）

场景：流体区域和固体区域分别用不同网格离散，交界面上通过位移/应力耦合。

SStruct 实现：

• Part 0：流体网格（如速度、压力变量）；
• Part 1：固体网格（如位移变量）；
• 在交界面网格点上，通过 AddEntries 添加耦合项（如流体压力作用于固体边界）；
• 可构建整体系统矩阵并用 AMG 或 GMRES 求解。

---

3.3 自适应网格细化（AMR）中的层次网格

场景：在 AMR 中，不同层级的网格块（粗网格 + 细网格块）需要耦合。

SStruct 支持：

• 每个 AMR 层级作为一个 part；
• 通过插值/限制算子在接口处耦合（需用户显式设置矩阵项）；
• 可结合 HYPRE 的 BoomerAMG 或 SMG 求解器。

注：HYPRE 也有专门的 AMR 接口（如 SStruct + HYPRE_SStructMatrixSetInterpData），但需手动管理插值权重。

---

3.4 多材料区域或多区域问题

场景：材料属性在不同区域突变（如复合材料、多孔介质），每个区域用独立网格描述。

SStruct 优势：

• 每个材料区域作为一个 part；
• 区域交界处设置界面条件（如通量连续、跳跃条件）；
• 变量类型可不同（如一个区域用标量，另一个用向量）。

---

3.5 周期性边界或重叠网格（Chimera网格）

场景：周期性边界（如圆柱绕流）或重叠网格（如动网格）。

SStruct 支持：

• 使用 HYPRE_SStructGridSetPeriodic 设置周期性；
• 使用 HYPRE_SStructGridSetNeighborBox 设置重叠区域映射；
• 在重叠区通过插值耦合矩阵项。

---

四、SStruct 与其他接口对比

接口	适用网格类型	灵活性	性能	典型应用
Struct	单一块结构化网格	低	最高	简单规则区域
SStruct	多块结构化/半结构化	中高	高	多块CFD、FSI、AMR
IJ	完全非结构化	最高	中等	有限元、复杂几何
ParCSR	非结构化（底层格式）	高	高（优化后）	IJ/SStruct 的底层实现

SStruct 在灵活性和性能之间取得良好平衡，是处理“局部结构化、全局半结构化”问题的理想选择。

---

五、注意事项与技巧

1. 变量索引管理：确保每个 part 的变量编号和 stencil 中的 var 一致。
2. 接口耦合需显式设置：SStruct 不自动处理块间耦合，必须手动调用 AddEntries 或设置插值。
3. 并行分区：每个 part 可独立分区，但需注意跨进程接口通信。
4. 求解器选择：
   • 对 Poisson 类问题：PFMG、SMG；
   • 对非对称/多物理场：GMRES + BoomerAMG；
   • 接口耦合可能导致矩阵非对称，慎用 CG。
5. 调试建议：使用 HYPRE_SStructMatrixPrint 输出矩阵结构检查耦合是否正确。

---

六、官方资源推荐

• HYPRE 官网：https://hypre.readthedocs.io
• SStruct 用户手册：https://hypre.readthedocs.io/en/latest/sstruct.html
• 示例代码：src/test/sstruct 目录下有多个完整示例（如 sstruct.c, sstruct_ls.c）

---

总结

HYPRE SStruct 接口适用于由多个结构化网格块组成的半结构化问题，在多块CFD、流固耦合、AMR、多材料模拟等领域有广泛应用。它在保持结构化网格计算效率的同时，提供了灵活的块间耦合机制，是处理复杂工程问题的强大工具。合理使用 SStruct 可显著提升大规模并行求解的效率和可扩展性。

如需进一步深入，建议从官方示例 sstruct.c 入手，逐步构建多部件耦合系统。