关于《Mixed Finite Elements in IR³》的学术报告
本文档是一篇发表于《Numerische Mathematik》期刊第35卷（1980年）第315-341页的原创性研究论文。作者是法国巴黎综合理工学院应用数学中心的 J. C. Nédélec。论文的核心贡献在于，为解决三维空间中的Maxwell方程组和弹性力学方程组等矢量场问题，提出并系统构建了两类全新的非协调有限元族。
一、 研究背景与目标
学术领域： 本研究隶属于计算数学和工程计算领域，具体方向为有限元方法（Finite Element Method, FEM）的开发与应用，特别是针对混合有限元（Mixed Finite Elements）和矢量有限元（Vector Finite Elements） 的构造。
研究动因与背景： 在20世纪70年代末，有限元法已成功应用于许多标量场问题（如热传导、结构力学中的位移场）。然而，当问题涉及具有特定物理约束和连续性的矢量场时（如电磁场中的电场强度E和磁感应强度B，或弹性力学中的应力和应变），传统的基于节点标量插值的有限元方法遇到了根本性困难。这些困难主要体现在： 1. 物理定律的精确满足： 在电磁学中，Maxwell方程组要求磁场B是无散的（div B = 0），而电场E的旋度（curl E）与磁场变化相关。在界面上，物理定律要求法向分量（如B·n）或切向分量（如E×n）连续。传统有限元难以在离散层面天然地保持这些微分约束和界面连续性条件。 2. 函数空间的适配： 描述这类矢量场需要特定的Sobolev空间，即H(curl)空间（函数及其旋度平方可积）和H(div)空间（函数及其散度平方可积）。构造能“保形”（conforming）地嵌入这些空间的有限元（即有限元函数本身属于这些连续空间）是一个重大挑战。在此之前，P. A. Raviart和J. M. Thomas等人已在二维问题上取得了突破，但三维问题的构造更为复杂。
研究目标： 本研究的核心目标是系统地构造出适用于三维空间的两族新的有限元： 1. 在四面体网格上构造的、保形于H(curl)空间的有限元族。 用于近似电场E等场量，其切向分量在单元边界上连续。 2. 在四面体网格上构造的、保形于H(div)空间的有限元族。 用于近似磁场B等场量，其法向分量在单元边界上连续。 此外，论文也简要给出了这两族有限元在立方体（六面体）单元上的对应构造。最终目标是将这些新元素应用于Maxwell方程组和弹性力学方程组的数值近似。
二、 研究方法与详细流程
本研究是理论构造与数学证明性质的，不涉及实验或数值模拟。其工作流程可概括为以下几个核心步骤：
步骤一：构建理论基础与多项式空间 为了构造满足特定微分性质（与curl和div算子相关）的有限元，作者首先在n维空间（主要关注n=2,3）中定义并研究了一系列特殊的多项式矢量空间。这是整个构造的代数基础。 1. 定义对称化微分算子 e^k u： 对于一个k+1阶可微的矢量函数u，定义其k阶对称化微分形式e^k u（公式(1)）。当k=1时，e^1 u即为弹性力学中的应变张量。该算子的关键性质是，当e^k u = 0时，蕴含着u的特定高阶偏导数关系（引理1）。 2. 定义核心多项式空间 N_k 和 S_k： 基于算子e^k u，定义空间 N_k = {u ∈ (P_k)^n : e^k u = 0}，其中P_k是次数≤k的多项式空间。S_k是N_k中的齐次（k次）多项式子空间。论文详细计算了n=2,3时这些空间的维数和具体形式（例如，N1对应三维空间中的刚体运动模式：平移和旋转）。这些空间与curl算子的核（即无旋场）密切相关（引理3）。 3. 定义用于H(div)的多项式空间 ID_k： 定义 IDk = (P{k-1})^n ⊕ P_{k-1}^0 · r，其中r是位置向量(x1, x2, …)^T。这个空间被设计为使得其元素与位置向量r的点积具有特定性质（命题1后的备注）。该空间在适当的仿射变换下具有不变性（引理5），并且其散度为零的子空间ID_k^0与S_k空间通过curl算子同构（命题3），这揭示了两个目标有限元族之间的深刻内在联系。
步骤二：构造四面体上的H(curl)保形有限元 这是本研究的第一个主要成果。作者为任意阶次k定义了一个有限元三元组（单元K，形函数空间P，自由度Σ）。 1. 单元K： 三维四面体。 2. 形函数空间P： 采用步骤一中定义的 N_k 空间。 3. 自由度Σ（定义4）： 自由度是作用在形函数上的一组线性泛函，用于唯一确定形函数。它们被精心设计以保证切向连续性： * 类型1（边自由度）： 在每条边a上，对切向分量u·t沿边积分，权函数取为P_{k-1}中的多项式。共6k个自由度。 * 类型2（面自由度）： 在每个面f上，对切向分量uf（二维向量）沿面积分，与(P{k-2})^2中的矢量函数作内积。共4k(k-1)个自由度。 * 类型3（体自由度）： 在单元内部，对u进行体积分，与(P_{k-3})^3中的矢量函数作内积。共k(k-1)(k-2)/2个自由度。 自由度的总数等于N_k空间的维数nk = (k+3)(k+2)k/2。 4. 关键证明： 论文通过一系列引理和定理证明了该元素的唯一可解性（unisolvent）（定理1）和H(curl)保形性。证明的核心逻辑是：首先证明若一个形函数的所有自由度为零，则其旋度必为零（引理7）；进而由旋度为零推出该形函数是某个势函数的梯度；再利用边界自由度（类型1，2）为零证明该势函数在边界上为常数；最后利用内部自由度（类型3）为零证明该势函数在整个单元上为零，从而形函数为零。保形性的证明（命题5）则表明，只要两个相邻单元在公共面上的边和面自由度一致，则切向分量在该面上自然连续。 5. 误差估计（定理2）： 通过仿射变换到参考单元，并利用插值算子的性质，论文给出了插值误差的估计：|u - Πu|{H(curl)} ≤ C h^k |u|_{H^{k+1}}，其中h是单元尺寸，k是多项式阶次。这证明了方法的收敛性。
步骤三：构造四面体上的H(div)保形有限元 这是与第二步并列的第二个主要成果，构造思路相似但目标空间不同。 1. 单元K： 三维四面体。 2. 形函数空间P： 采用步骤一中定义的 ID_k 空间。 3. 自由度Σ（定义5）： 设计以保证法向连续性： * 类型1（面自由度）： 在每个面f上，对法向分量u·n沿面积分，权函数取为P{k-1}中的多项式。共4 * (k(k-1)/2) = 2k(k-1)个自由度。 * 类型2（体自由度）： 在单元内部，对u进行体积分，与(P{k-2})^3中的矢量函数作内积。共3 * ((k+1)k(k-1)/6) = k(k-1)(k+1)/2个自由度。 自由度总数等于ID_k空间的维数m_k = (k+3)(k+1)k/2。 4. 关键证明： 同样证明了唯一可解性和H(div)保形性（定理3）。证明逻辑为：由面自由度为零推出法向分量在面上为零；由体自由度为零和散度定理推出散度为零；结合ID_k空间的结构，最终证明形函数为零。 5. 误差估计（定理4）： 给出了类似的插值误差估计。
步骤四：构造立方体（六面体）上的对应有限元族 为了网格剖分的灵活性，论文简要地将上述构造推广到立方体单元上（第2节）。利用张量积形式的多项式空间（如Q{l,m,n}表示在x1, x2, x3方向最高次数分别为l, m, n的多项式），分别定义了保形于H(curl)（定义6，空间为Q{k-1,k,k} × Q{k,k-1,k} × Q{k,k,k-1}）和H(div)（定义7，空间为Q{k,k-1,k-1} × Q{k-1,k,k-1} × Q_{k-1,k-1,k}）的有限元，并给出了相应的自由度定义和唯一可解性证明（定理5，7）及误差估计（定理6，8）。这些构造是四面体构造的类比，但基于张量积空间。
步骤五：应用示例 论文在第三部分简要概述了新元素在几个重要问题中的应用框架，展示了其应用价值： 1. 电介质中的Maxwell方程组（3.1节）： 提出将电场E用H(curl)元离散，磁场B用H(div)元离散的混合变分形式。这样，离散解自动满足磁场无散（div B_h = 0）和界面法向/切向连续性条件，严格保持了物理定律的离散形式。 2. 导体中的涡流方程（3.2节）： 给出了基于新元素的另一种混合变分形式。 3. 弹性力学与Stokes方程（3.3节）： 针对具有特定边界条件（切向位移固定，法向应力自由）的线弹性问题和Stokes流问题，提出了一种基于H(curl)元离散位移、基于标准标量元离散压力（体积变化）的混合变分格式。论文甚至利用新元素的特性（引理11），证明了该格式满足离散的inf-sup条件（定理9），从而保证了数值方法的稳定性，并给出了误差估计。
三、 主要研究结果
本研究的主要成果是理论上的构造和证明： 1. 成功构造了两类新的三维有限元族： 明确给出了在四面体和立方体单元上、任意阶次k的、分别保形于H(curl)和H(div)空间的形函数空间和自由度集。 2. 完成了严格的理论证明： 对每一类元素，都证明了其唯一可解性、保形性以及插值误差估计。这些证明构成了完整的数学理论体系。 3. 揭示了深刻的数学联系： 通过多项式空间N_k, ID_k及其与微分算子curl, div的关系（如命题3指出curl是S_k到ID_k^0的同构），阐明了两类元素并非孤立存在，而是通过微分算子紧密关联的对偶体系。这种结构上的优美是方法成功的关键。 4. 搭建了应用于重要物理问题的桥梁： 通过第3节的应用示例，具体展示了如何将新元素应用于Maxwell方程、涡流方程和弹性力学的混合变分形式，并证明了在特定边界条件下弹性问题格式的稳定性。这些应用框架为后续大量的工程和科学计算研究奠定了基础。
四、 结论与价值
结论： J. C. Nédélec在本研究中系统性地解决了三维矢量有限元构造的核心难题，提出了以他名字命名的Nédélec棱边元（H(curl)元）和Nédélec面元（H(div)元）（尽管论文中未直接命名，但已被学术界广泛称为Nédélec elements）。这些元素能够精确描述矢量场在界面上的切向或法向连续性，并自然满足相应的微分约束。
科学价值与应用价值： 1. 理论价值： 丰富了混合有限元方法理论，为在H(curl)和H(div)空间中设计高阶、高精度离散格式提供了通用框架。论文中关于多项式空间和微分算子关系的研究具有基础性意义。 2. 应用价值： 极大地推动了计算电磁学（CEM）的发展。Nédélec元已成为时域和频域有限元法（FEM/FETD/FEFD）求解Maxwell方程组的标准工具，被广泛应用于天线设计、微波器件、光学、电磁兼容等领域。它们也被用于计算流体力学（CFS）中的某些问题，以及涉及应力/应变张量的固体力学计算。 3. 工程价值： 使得有限元法能够更物理、更稳健地处理复杂的多物理场耦合问题，其中不同场量需要满足不同的连续性条件。
五、 研究亮点
六、 其他有价值内容
论文中还包含一些有价值的细节，例如：讨论了自由度在仿射变换下的不变性问题，并给出了更具不变性的实用自由度选择建议（如用边矢量代替单位切向矢量）；在证明弹性问题格式稳定性时（定理9），引入了一个基于势函数分解的巧妙论证（引理11及后续推导），展示了新元素在分析复杂混合系统时的独特优势。