概述
RT 空间（Raviart–Thomas space）是一类用于有限元方法的向量场函数空间，尤其适合处理混合有限元与守恒型（conservation law）偏微分方程。其关键构造思想是在网格单元上用多项式函数刻画向量场，同时在单元边界上施加特定连续性条件。最常见的情形是要求向量场的法向分量在相邻单元面之间具有合适的连续性，从而使离散解在“通量”层面更符合守恒律。
在数值分析实践中，RT 空间通常与“标量/压力”等辅助变量配对使用，形成混合离散格式。通过选择多项式次数、网格几何映射与自由度组织方式，RT 空间在稳定性（离散守恒与适定性）、收敛性以及工程实现之间取得平衡。此外，RT 家族还存在多种变体，例如不同阶次或低阶/高阶版本，用于在计算成本与精度之间做权衡。
1. 概念与动机
1.1 守恒律与通量连续的数值需求
许多守恒型模型的核心物理量不是直接的“量本身”，而是“量的流动或传输”，即通量。连续性条件常常体现在：在单元面上，通量从一侧进入与从另一侧离开应当匹配（在合适意义下）。若离散空间不保证这种边界层面的连续性，容易出现局部守恒误差，导致数值解出现非物理解、振荡或误差累积。
因此，构造能自然表达通量并在界面上维持正确连续性的函数空间，是保证守恒离散质量的一条主线。
1.2 为什么选“向量场”而非仅标量场
在许多问题中，未知量或关键中间量是向量形式：例如电磁学中的电/磁通量、流体与多孔介质中的达西速度（速度向量）、混合形式中对梯度的再变量化等。仅使用标量场通常只能在弱形式下“间接”表示通量，往往需要额外处理才能获得良好的守恒性质。
RT 空间直接把通量变量以向量场方式嵌入有限元空间，使“法向通量”的几何与物理含义更清晰，进而更容易控制界面上的一致性。
1.3 Raviart–Thomas 家族的基本定位
Raviart–Thomas 家族是一类满足特定相容性要求的向量有限元空间，典型特点是与 \(H(\mathrm{div})\) 相关的结构相匹配：它们强调向量场的散度（divergence）与法向分量在网格层面的可控性。借助这种相容性，RT 空间经常用于混合有限元法中，把“流量/通量”与“压力/势函数”等变量组合成稳定的离散体系。
2. 函数空间的定义
2.1 在参考单元上的构造思路
RT 空间的构造通常遵循“从参考单元到实际单元”的范式。首先在参考几何（如参考三角形、参考四面体或参考四边形/六面体的变体）上定义局部向量多项式空间；随后通过映射（仿射或更一般的几何映射）将基函数迁移到物理网格单元。这样做的好处是：局部多项式结构可以统一管理，而全局网格几何的影响通过映射进入。
2.2 典型的局部多项式结构
RT 家族的局部空间通常由“固定次数的多项式向量场”与“与位置变量相关的附加项”构成。直观理解是：向量场不仅允许有常规多项式变化，还允许带有与坐标方向相关的结构，从而使其法向分量落在合适阶次的面多项式空间中。通过这种精心选取，RT 基函数能够在面上形成合适的自由度，从而实现边界连续性与良好近似性。
2.3 法向分量连续性（H(div) 相关）要求
RT 空间最显著的连续性特征是对法向分量的控制。在二维/三维网格中，单元共享的公共面（边/面）上，向量场的法向分量在离散意义下保持连续或满足配套的“拼接条件”。这种性质与 \(H(\mathrm{div})\) 空间的相容性直接相关：离散散度算子与法向通量的界面耦合更自然，从而提升混合格式的稳定性表现。
2.4 与网格划分（单元、面、边）的关系
RT 空间的构造强依赖网格的拓扑划分。局部函数的自由度往往与网格的面（在二维即边）相关，而不是仅与单元内部点相关。这样一来，法向通量的连续性可以通过面上的自由度来实现，确保在单元之间“对得上”。在实现层面，这意味着：组装过程需要对共享面做一致的编号与约束/拼接。
3. 离散化与有限元形式
3.1 典型混合有限元框架
混合有限元通常把原问题改写为包含两个（或更多）变量的系统：一个变量更接近通量（向量），另一个变量更接近势函数或标量未知量。RT 空间通常被用于通量变量的离散空间；而压力或标量变量则由另一类标量有限元空间提供。
在弱形式中，会出现散度算子与梯度/对偶算子的耦合项。RT 空间的相容性使这些耦合在离散层面保持良好结构，避免不适定或严重锁定现象。
3.2 变分形式中的角色：通量变量与压力/标量变量
在典型混合格式中，通量变量被放在向量空间（RT 或其变体），其散度对应到标量变量的源项或守恒关系；标量变量通常通过拉格朗日乘子或直接的变分项与通量互相作用。结果是：守恒律可从离散方程组中体现为更直接的代数平衡，而非仅依靠数值后处理实现“近似守恒”。
3.3 代数结构：自由度与基函数的选择
将局部基函数拼接到全局时，会把自由度组织为“面相关”的模式。由于法向分量连续性被严格编码到空间定义中，离散系统的矩阵结构也往往呈现与守恒算子对应的特性。合理选取基函数与自由度，可以减少“塌缩”（空间退化到不合适的子空间）的风险，并提升数值求解稳定性。
3.4 与其他空间的耦合（如标量空间的选择）
RT 空间的稳定表现依赖于与标量空间的配对方式。一般而言，标量空间需要与向量空间的散度像（或相近结构）在阶次与维数上匹配，才能满足常见的适定性条件（如对应的离散交叉条件）。因此，实践中通常会选择与 RT 阶的阶次相容的标量有限元空间，从而确保混合系统不仅能收敛，而且不会产生不合理的自由度失衡。
4. 性能与理论性质
4.1 稳定性：离散守恒与适定性（概念层面）
在理论层面，RT 空间常用于构造满足离散适定性（well-posedness）的混合有限元离散。稳定性的关键直观来源包括：
- \(H(\mathrm{div})\) 相容性让散度与法向通量之间的结构在离散层面保持一致；
- 面相关的自由度安排使通量在界面上更容易满足局部平衡；
- 与标量空间的配对阶次控制了自由度之间的匹配关系。
这些因素共同改善离散系统的条件数表现与误差传播行为。
4.2 近似性质与收敛阶（概念层面）
RT 空间的近似能力由其多项式次数与几何映射决定。一般而言，随着阶次提高，能够更好逼近光滑解中的向量场及其散度结构；对混合问题而言，误差不仅体现在通量变量，也体现在压力/标量变量。收敛阶通常与所选 RT 阶次以及配套标量空间阶次有关，并在足够规则的网格与解光滑性条件下表现为稳定的代数增长趋势。
4.3 保守性与局部通量平衡
局部守恒是 RT 空间的亮点之一。由于法向通量的拼接方式与 \(H(\mathrm{div})\) 结构相一致，在满足混合方程组的前提下，离散解往往能在单元或控制体意义下实现通量平衡。对工程计算而言，这通常意味着：局部误差不易“凭空产生”，守恒相关的量更可控。
4.4 与网格畸变及高阶基函数的关系（概念层面）
数值性能会受到网格质量与基函数阶次的共同影响。网格畸变可能改变几何映射的稳定性，进而影响误差常数；而高阶基函数虽然提高逼近能力，但也可能带来更严格的实现要求，例如数值积分精度、映射误差控制与求解器对高维自由度的适配。因此在实际计算中，往往需要在网格质量、阶次选择与计算成本之间做综合平衡。
5. 自由度（Degrees of Freedom）与实现
5.1 法向通量相关自由度的含义
RT 空间的自由度常与面上的法向分量相关。直观上，可以把某个单元面上的若干“通量矩”（例如法向分量与面多项式的加权积分）作为自由度。这样定义的优势在于：当相邻单元共享同一面时，只要它们在该面上的法向自由度匹配，便能自然实现法向连续性。
5.2 基函数构造的常见策略（概念层面）
基函数构造常遵循“确保自由度可识别”的原则：基函数应满足在各面自由度上呈现类似 Kronecker delta 的性质，从而让系数可以直接由自由度求得。实际实现中，常用参考单元上解析构造，再通过映射与适配规则生成物理单元基函数。对于高阶情形，需要更多面多项式自由度与对应的基函数组织方式来覆盖更高的变化模式。
5.3 组装流程：局部到全局的映射
组装通常包含以下关键环节：
- 在每个单元上计算局部矩阵与局部载荷；
- 将局部自由度与全局编号对应；
- 处理共享面的连续性要求（通过全局自由度的唯一编号实现拼接，或在离散框架中通过约束结构实现）；
- 汇总得到全局线性系统。
由于 RT 自由度主要与面相关，面共享关系是组装正确性的基础。
5.4 质量检查：数值守恒与误差评估
实现完成后，常见的质量检查包括：
- 检查离散通量的界面连续性是否符合预期（例如面上法向分量的拼接误差）；
- 在可验证基准算例上评估收敛曲线；
- 若问题具有天然守恒量，则检查离散解是否体现局部或全局守恒趋势。
这些检查有助于发现自由度定义、映射、数值积分或配对空间阶次不匹配等问题。
6. 常见阶次与变体
6.1 RT0、RT1 等低阶空间的直观理解（概念层面）
低阶 RT 空间提供了最直观的通量离散方式：法向分量在每个面上呈现相对简单的多项式结构，因此自由度数量较少、实现相对直接。以 RT0 为例，其通量在面上的表达能力较基础，但常常足以用于许多工程需求或作为更高阶方法的基准。
6.2 高阶 RT 空间的递推/升级思路（概念层面）
高阶 RT 空间可通过提高局部多项式次数来增强近似能力。升级通常伴随自由度的增加：面上的自由度数量随着阶次增长而增长。实现时往往需要更高精度的数值积分以及更细致的基函数与自由度管理，以确保理论阶次在离散层面真正体现。
6.3 RT 与其他 H(div)-相容空间的比较
在保持 \(H(\mathrm{div})\) 相容性的框架下，除了 RT 空间外还存在其他相关构造（例如不同家族的向量有限元）。比较时通常看：
- 离散稳定性与适定性表现；
- 近似阶次与常数；
- 自由度规模与实现复杂度；
- 与特定物理模型的耦合便利性。
在一些情形下，不同家族在性能或实现上可能各有优势，但 RT 空间因其结构清晰与经典应用而仍非常常见。
6.4 特殊网格与几何映射下的注意事项
当网格包含畸变单元、非仿射映射或更一般的网格拓扑时，需要注意映射规则对基函数性质的影响。若映射导致离散空间偏离预期的等价结构，可能影响法向连续性或近似性质。工程上通常通过提高几何处理质量、选择合适的映射精度或采用兼容的离散策略来降低此类问题。
7. 应用领域
7.1 泊松方程混合形式与通量-势函数表述
泊松方程常可写成“势函数-通量”耦合的混合形式，例如把通量视为势函数梯度的再变量化。RT 空间作为通量离散空间，能更自然地处理由梯度与散度之间关系形成的守恒结构，从而在数值上兼顾稳定性与通量合理性。
7.2 流体与多孔介质：Darcy 型问题（概念层面）
Darcy 型问题中的速度通常体现为某种与压力梯度相关的向量场。采用 RT 空间离散通量变量时，法向通量自由度可以直接对齐界面输运的物理含义，使局部守恒与流量平衡更易实现。在多材料或不规则边界条件下，RT 离散常用于保证通量计算的可靠性。
7.3 电磁与波动问题中的通量离散（概念层面）
在电磁学的某些表述中，相关的守恒或边界传输量可以被视为向量通量。RT 类空间强调法向通量的连续控制，因此在与电磁模型的特定变量选取相匹配时，可用于实现界面上更合理的离散传输行为。实际选择往往取决于模型采用的变量形式与边界条件类型。
7.4 工程守恒问题中的数值实践
工程场景中，常见需求包括局部守恒、通量核算准确以及混合方程组的稳定求解。RT 空间因其面相关自由度与 \(H(\mathrm{div})\) 相容性，适合将这些需求转化为可验证的离散性质。结合合适的标量配对空间与求解器策略，能在保持物理一致性的同时获得较好的数值效果。
8. 相关概念与延伸阅读
8.1 H(div) 空间与匹配连续性
RT 空间与 \(H(\mathrm{div})\) 的关系是其理论基础之一。理解 \(H(\mathrm{div})\) 的法向追踪（界面连续性）与散度定义，有助于把 RT 空间的连续性条件理解为对守恒结构的离散对应。
8.2 与 Nédélec（H(curl)）等体系的互补关系（概念层面）
在守恒型离散中，除了强调 \(H(\mathrm{div})\) 的空间外，电磁相关问题还常涉及 \(H(\mathrm{curl})\) 的相容性。Nédélec 类型空间通常用于对偶的“旋度/切向”结构进行相容离散。两者在许多混合或耦合框架中具有互补性：一个侧重法向通量，一个侧重切向场或旋度结构。
8.3 混合有限元的条件与配对原则
混合有限元并非只选某个空间即可获得稳定性，还需要合适的空间配对。常见原则包括：向量空间的散度像与标量空间的阶次匹配、自由度维数在离散约束下不失衡等。理解这些配对原则有助于正确选择 RT 空间的阶次及其标量伴随空间。
8.4 进一步资料：经典文献脉络（方向性）
深入学习通常可从以下方向入手：有限元函数空间的构造思想、\(H(\mathrm{div})\) 相容与守恒离散理论、混合有限元的适定性分析，以及 RT 家族与相关家族在不同 PDE 变量选取下的应用对照。阅读经典文献时，可以重点关注从理论构造到数值实现的桥梁环节，例如自由度定义、映射规则与误差估计的对应关系。