《Energies》：Valorization of Sargassum via Hydrothermal Co-Liquefaction with Food Waste: Effects of Reaction Temperature and Feedstock Ratio on Biocrude Yield and Fuel Properties
Md Mostafizur Rahman and
Toufiq Reza
编辑推荐：
本文提出了一种名为径向基函数-有限差分矩阵算子（RBF–FDMO）的混合数值方法，旨在提高传统的有限差分矩阵算子（FDMO）技术在求解由拉普拉斯-泊松（Laplace–Poisson）方程支配的三维（3D）电磁场分析时的准确性和灵活性。传统的数值方法在处理非均
本文提出了一种名为径向基函数-有限差分矩阵算子（RBF–FDMO）的混合数值方法，旨在提高传统的有限差分矩阵算子（FDMO）技术在求解由拉普拉斯-泊松（Laplace–Poisson）方程支配的三维（3D）电磁场分析时的准确性和灵活性。传统的数值方法在处理非均匀离散化和复杂几何形状时，常面临计算复杂度和灵活性有限的问题。在所提出的方法中，空间导数的近似采用基于径向基函数（RBF）的插值而非源自泰勒（Taylor）级数展开的有限差分格式。该公式使得能够在均匀和非均匀有限差分（FD）网格上构建高精度的导数算子，从而提高了数值鲁棒性以及对复杂几何形状的适应性。研究人员首先将所提方法的性能与FDMO在三维基准问题中进行了比较，其均方根（RMS）误差和最大误差均降低了两个数量级以上。此外，研究人员开发并首次将RBF-FDMO方法应用于IEEE Std. 80?中规定的接地系统（GS）配置分析，以及越南一个实际的110 kV变电站接地系统分析。所获得的电位分布、接地电阻以及接触电压（touch voltage）和跨步电压（step voltage）证实了该方法的有效性和可靠性。结果表明，所提方法具有公式简单和计算效率高的特点，使其成为接地系统分析和设计中有限元法（FEM）和边界元法（BEM）等传统方法的实用替代方案。
论文解读：基于RBF–FDMO方法的接地系统三维数值模拟研究
研究背景与意义
接地系统（Grounding System, GS）是电力系统设计与运行中的关键组成部分，主要承担系统工作接地、设备保护接地和防雷接地三大功能。当电力系统发生绝缘击穿、短路接地或雷击放电等故障时，接地系统为故障和浪涌电流提供低阻抗入地路径，从而保障保护装置的有效动作，维护人员安全以及整个电力系统的电磁兼容性与稳定性。准确求解描述该问题的三维拉普拉斯（Laplace）方程，对于评估接地电阻、接触电压（Touch Voltage）和跨步电压（Step Voltage）等关键安全参数至关重要，这也是确保符合IEEE Std. 80?等国内外标准的基础。
目前，分析具有不规则几何形状和非均匀土壤条件的复杂接地系统时，通常采用有限元法（FEM）、边界元法（BEM）或矩量法（MoM）等高级数值方法。然而，这些方法存在固有缺陷：FEM虽灵活准确，但涉及复杂的网格划分，会产生大规模稀疏矩阵，消耗大量计算资源和内存，且人工边界的引入会带来建模不确定性；BEM虽减少了维度，但在处理非均匀或多层土壤时数学复杂，且密集矩阵会导致内存需求和计算时间增加；MoM则受限于密集阻抗矩阵的增长及对简化土壤模型的依赖。这些传统方法普遍面临精细结构化网格要求高、矩阵系统庞大、对非均匀土壤适应性差等挑战。因此，研究人员开展了本研究，旨在开发一种更高效、灵活且保持数值精度的新型数值框架——RBF–FDMO方法，用于实际接地系统分析。该论文发表于《Energies》期刊。
主要关键技术方法
研究人员采用了一种混合数值方法，即将径向基函数（RBF）插值技术集成到有限差分矩阵算子（FDMO）框架中，形成RBF–FDMO方法。该方法通过RBF插值近似空间导数，替代了传统的泰勒级数有限差分格式，从而能在均匀和非均匀三维有限差分网格上构建高精度导数算子，并利用克罗内克积（Kronecker product）将一维算子扩展至三维。为了优化RBF的形状参数，研究人员采用了Victor Bayona等人提出的算法以确定最优形状参数。研究通过解决一个三维泊松（Poisson）方程基准问题验证精度，并将方法应用于四种接地系统案例：IEEE Std. 80?中的方形、矩形、L形接地网以及越南一个110 kV漂浮太阳能变电站的实际接地系统，使用非均匀FD网格进行离散化，并与BEM、FEM及IEEE Std. 80?的计算结果进行对比。
研究结果
4.1. 应用于三维基准问题的RBF-FDMO方法
研究人员通过求解一个三维泊松方程基准问题来评估RBF-FDMO的性能。结果表明，均方根（RMS）误差和最大误差在形状参数约为特定值时达到最小。在测试的多种RBF（GA, MQ, IQ, IMQ）中，逆二次（IQ）函数表现出更优的精度和数值稳定性平衡，因此被选用于后续接地系统模拟。与传统的3D FDMO方法相比，RBF-FDMO的RMS误差和最大误差均降低两个数量级以上，且能以显著更少的离散点达到同等精度，证明了其计算经济性。
4.2. 用于接地系统（GS）分析的RBF-FDMO方法
研究人员将IQ-FDMO方法应用于不同接地系统配置的电位分布及安全参数计算，土壤假设为均匀介质，远场边界设为狄利克雷（Dirichlet）边界，地表设为齐次诺伊曼（Neumann）边界。
4.2.1. 案例1：方形接地系统
针对IEEE Std. 80?中的R04P16方形接地网（含垂直接地极），研究人员使用非均匀FD网格和IQ-FDMO进行模拟。结果显示，该方法计算的接地电阻为0.402 Ω，地电位升（GPR）为20.1 kV，与BEM结果（0.4 Ω）高度吻合，较IEEE Std. 80?（0.426 Ω）降低5.6%。最大接触电压为28%，与BEM一致；最大跨步电压仅为6.0%，较BEM（16%）和IEEE Std. 80?（11.7%）显著降低，表明该方法在提升安全评估准确性及降低跨步电压危害方面的有效性。
4.2.2. 案例2：矩形接地系统
针对IEEE Std. 80?中63 m × 84 m的矩形接地网，研究人员利用非均匀FD网格计算。IQ-FDMO得到的最大电位1.014 p.u.和最小电位0.978 p.u.与FEM结果（1.013 p.u.和0.979 p.u.）接近，峰值电位变化0.036 p.u.略高于FEM的0.034 p.u.，源于网格构造差异，但证实了其准确性。安全参数方面，接地电阻0.886 Ω较IEEE Std. 80?（0.876 Ω）高1.14%，GPR为5.03 kV，接触电压14.9%较标准（14.8%）高0.68%，而最大跨步电压6.0%较标准（7.76%）降低22.7%，显示了其在安全评估中的可靠性与优势。
4.2.3. 案例3：L形接地系统
针对IEEE Std. 80?中带垂直接地极的L形接地网（主区105 m × 70 m，附加35 m × 35 m），研究人员采用非均匀FD网格模拟。IQ-FDMO与FEM得到的电位分布几乎一致，最大、最小电位及变化量差异极小，确认了该方法适用于任意几何形状接地系统的可靠性。安全参数对比显示，接地电阻（0.647 Ω）和GPR（5.23 kV）与IEEE Std. 80?一致；最大接触电压14.9%和跨步电压6.0%略低于标准的15.4%和8.1%，提供了更精细且保守的电压分布评估。
4.2.4. 案例4：越南实际接地系统
研究人员分析了越南Da Mi水电站水库上110 kV–63 MVA漂浮太阳能变电站的接地系统（含镀锌钢导体、垂直接地极及30 m深GEM井）。IQ-FDMO得到的电位分布与FEM高度吻合，倾向于预测稍大的电位变化，利于安全合规。计算的接地电阻0.5 Ω与现场测量值0.49 Ω仅差0.01 Ω。在安全评估中，归一化最大接触电压（10.3%）远低于IEEE Std. 80?限值25.1%（安全裕度59.8%），归一化最大跨步电压（3.0%）远低于12.98%限值（裕度76.7%），系统完全满足安全标准，无需额外修正。
讨论与结论
研究人员在讨论中指出，与FEM、BEM和MoM相比，有限差分法（FDM）因其概念和计算的简洁性而突出，RBF–FDMO框架通过克罗内克积高效构建三维算子，适用于大规模域和复杂几何形状。在形状参数较小时，RBF近似取代了经典的泰勒级数有限差分近似，显著提高了精度；结合最优形状参数选择算法，可进一步提升性能。该方法在离散域中使用仅几千个节点即可达到可比或更优精度，而传统FDMO需数十万个点，极大减轻了计算负担，适合大规模或几何复杂问题。此外，该框架可适应多种接地系统配置，具有低计算复杂度和简单实现的特点，是高压变电站接地设计和优化的有效工具，可按IEEE、IEC等标准评估安全参数。
结论部分总结道，本研究开发了RBF–FDMO方法，通过克罗内克积将此前提出的FDMO公式扩展至全三维问题，无需FEM求解器所需的复杂网格生成或区域离散化。通过整合RBF，显著增强了导数近似的准确性、平滑性和数值稳定性，在精度和计算成本间取得更好平衡。需注意本研究基于理想条件（如均匀土壤和标准配置），计算效率主要体现在算子公式的简洁及较少离散点达到高精度，未来将扩展至瞬态、频变分析及更复杂土壤模型，并集成至商业软件。该研究的创新在于构建了统一均匀/非均匀节点网格上的三维RBF增强有限差分算子方案，兼具FDM简洁性和RBF处理不规则节点的灵活性，通过纯代数算子易于实现，适用于多样接地配置。经基准问题验证及实际110 kV变电站应用，结果与FEM、BEM高度一致，且公式更简单、计算效率具竞争力。总之，该研究引入了一种新型三维RBF–FDMO公式，结合了FDM、RBF插值及克罗内克积算子构建的优势，具备数学简洁、高精度和高效率的特点，弥合了传统基于网格的求解器与新兴无网格技术的差距，在高压变电站早期设计和优化任务中潜力显著。
