《Computational Particle Mechanics》：Implementation of a single-phase SPH pneumatic power take-off model for small- and full-scale OWC devices within DualSPHysics
编辑推荐：
【编辑推荐】本研究针对固定式和浮动式振荡水柱(OWC)波能转换器的气动动力输出(PTO)建模难题，提出了一种在无网格SPH框架内融合解析腔室模型的单相模拟新方法。该方法在计算成本可控的前提下，有效捕捉了腔内空气的可压缩性效应，并针对不同尺度OWC装置进行了实验验证。其意义在于为OWC性能的精准评估及优化设计提供了一种高效可靠的新型模拟工具。
在寻求可持续清洁能源的时代背景下，海洋波浪能以其高能量密度和可预测性，成为了极具潜力的可再生能源选择。其中，振荡水柱(OWC, Oscillating Water Column)波能转换器(WEC, Wave Energy Converter)因其结构简单、研究深入且原型部署广泛，被认为是目前最成熟的波能技术之一。其核心工作原理是利用波浪驱动一个腔室内的水柱振荡，进而压缩或扩张其上方的空气，空气压力变化驱动涡轮机发电，这个系统被称为气动动力输出(PTO, Power Take-Off)。
然而，要精确评估OWC的能量转换性能，特别是预测其在真实海洋（全尺度）环境下的表现，准确模拟其内部空气与水的相互作用至关重要。传统的计算方法，如基于网格的计算流体动力学(CFD)，虽然可以实现高保真的两相（水-气）模拟，但其在模拟自由表面流动、浮体运动以及复杂的网格重构方面面临挑战。与此同时，无网格方法，特别是光滑粒子流体动力学(SPH, Smoothed Particle Hydrodynamics)，在处理这类问题时展现出独特优势，但其在模拟OWC的PTO时也面临两难：进行两相模拟会带来极高的计算成本，而传统的单相简化模型（如在自由表面添加浮动板）又可能无法准确反映空气压力对水柱的动态反馈效应，尤其是在考虑空气可压缩性时。那么，能否找到一种计算高效、物理准确的SPH模型，来模拟OWC的气动PTO，并揭示空气可压缩性在不同尺度下的影响呢？
为了回答这个问题，Beatrice Mina、Bonaventura Tagliafierro、Iván Martínez-Estévez、Corrado Altomare、Moncho Gómez-Gesteira、Giovanni Malara和Felice Arena等研究人员在《Computational Particle Mechanics》期刊上发表了一篇研究论文，提出并验证了一种新的单相SPH建模方法，用于模拟OWC的气动PTO，该方法能够同时兼顾计算效率和物理准确性，并首次在SPH框架下研究了浮动OWC中的空气可压缩性效应。他们的核心策略是：不直接模拟空气粒子，而是将水相SPH模型与一个描述腔室内空气压力变化的解析模型相耦合。当水柱在腔室内振荡时，通过追踪水粒子来计算腔内的水体积，从而反推空气体积的变化。基于空气体积变化率，可以计算出通过代表涡轮机的“孔口”的空气流量，再通过一个标准的孔口压力降方程，计算出施加在水粒子上的等效空气压力，并将其作为源项加入SPH的动量方程中，从而实现空气压力对水柱运动的动态影响。该研究还进一步将模型扩展，包含了空气的可压缩性效应，使其能够应用于对压缩性敏感的大尺度模拟。
为了开展研究，研究人员主要采用了以下几个关键技术方法：首先，在开源软件DualSPHysics中实现了所提出的单相SPH气动PTO模型。其次，利用该模型分别对固定式OWC、U型OWC以及浮式（垂荡）OWC装置进行了数值模拟。最后，将数值模拟结果与现有文献中的实验数据进行对比，以验证模型在不同结构（固定/浮动）和不同尺度（实验室/全尺度）下的准确性与可靠性。
研究结果：
- •
模型验证与固定OWC性能分析：研究人员首先将提出的单相SPH模型（采用不可压缩空气假设）应用于固定式OWC和U-OWC的模拟，并将结果与实验室实验数据进行对比。结果显示，模型能够准确地复现实验观测到的腔内水位振荡和压力变化。随后，他们利用该模型研究了空气可压缩性对固定OWC性能的影响。结果表明，在模型尺度（实验室尺度）下，采用不可压缩或可压缩假设得出的能量吸收功率预测值相差不大。然而，当模拟全尺度OWC时，空气可压缩性的影响变得显著。如果忽略空气可压缩性，会导致对装置吸收功率的高估。这一发现与基于网格的CFD方法的成熟结论一致，从而证明了所提出的SPH模型能够捕捉到这一关键的物理现象。
- •
浮式OWC的模型验证：该研究进一步将模型应用于一个结构更为复杂的浮式（垂荡）OWC装置。这个装置同时结合了气动和液压两种PTO。研究人员分别用不可压缩和可压缩两种腔室模型公式进行了模拟，并与实验室测试数据进行了对比。验证结果表明，即使在装置同时存在浮体运动和复杂PTO配置的情况下，所提出的SPH模型依然能够准确地预测装置的动力学响应和性能。这项工作是首次在考虑空气可压缩性的前提下，在SPH框架下对浮式OWC装置进行的研究，展示了该方法处理复杂浮式结构气动PTO问题的能力。
- •
算法效率与可扩展性：文中特别指出，所采用的通过追踪腔内水粒子来计算空气体积的算法，其计算复杂度为O(n)（n为粒子数）。这意味着无论是对于固定OWC还是浮式OWC（二维或三维），其计算成本都仅与粒子总数成线性关系，具有良好的计算效率和可扩展性，为未来三维及更复杂场景的应用奠定了基础。
结论与讨论：本研究成功地在DualSPHysics中开发并验证了一种新颖的单相SPH气动PTO模型，该模型通过耦合解析腔室模型，能够有效地模拟固定和浮式OWC装置在不同尺度下的性能。模型的核心优势在于，它在保持单相SPH计算效率的同时，通过间接但物理合理的方式纳入了空气压力的动态效应，特别是空气可压缩性的影响。研究表明，空气可压缩性在评估全尺度OWC性能时不可忽略，而所提出的模型能够准确地反映这一点。更重要的是，该模型在浮式OWC上的成功应用，证明了其在处理具有复杂运动和多物理场耦合的海洋能装置方面的潜力。这项工作不仅为OWC的PTO建模提供了一种高效、可靠的SPH替代方案，也为未来利用无网格方法进行波能转换器的优化设计和性能预测开辟了新的途径。
