加速声学生成聚集体预测:一种融合声压梯度和斯托克斯拖曳傅里叶特征的物理信息神经网络

《Computer Physics Communications》:Accelerating acoustic agglomeration predictions: A physics-informed neural network with fourier featuring for sound pressure gradient and stokes drag

【字体: 大 中 小 】 时间:2026年10月23日 来源:Computer Physics Communications 3.9

编辑推荐:

   •利用PINN进行大规模声诱导粒子运动预测。•分离训练与积分以提升建模效率。•PINN预测结果与FEM仿真高度一致。•FF-PINN可准确预测粒子运动与碰撞行为。•FF-PINN模型仿真速度远快于FEM。引言颗粒物(PM)是对人类健康的主要危害之一,亟需有效的去除策略[1]。

  
  • •
    利用PINN进行大规模声诱导粒子运动预测。
  • •
    分离训练与积分以提升建模效率。
  • •
    PINN预测结果与FEM仿真高度一致。
  • •
    FF-PINN可准确预测粒子运动与碰撞行为。
  • •
    FF-PINN模型仿真速度远快于FEM。

引言

颗粒物(PM)是对人类健康的主要危害之一,亟需有效的去除策略[1]。特别是像PM2.5这样的细颗粒物,由于其粒径微小,传统过滤方法难以捕捉。为克服这一局限,研究人员探索了在去除前增大颗粒粒径的预处理技术。其中,声学聚结(AA)在烟气处理、柴油排气系统和环境空气过滤等领域引起了越来越多的关注[[2], [3], [4], [5]],因为它能促进颗粒聚结并便于后续去除[6]。
AA描述了声场激发下颗粒的相互作用、碰撞与聚结过程,在此过程中声波将机械能传递给悬浮颗粒[6,7]。在驻波场中,颗粒通常趋向于迁移到压力节点处,此处的振荡流体运动在颗粒与载流介质之间产生相对运动,从而产生Stokes阻力[8,9]。这种声学运动增强了碰撞概率,导致更大聚集体形成。随着过程持续,颗粒在节点区域积累,形成更容易被去除的团簇。与传统去除方法相比,AA具有可扩展性、相对较低的能耗以及在下游分离前通过增强聚结来改善细颗粒捕集效率等优势[10,11]。
AA实验性能取决于多种关键参数,这些参数通过调节底层物理机制共同决定净聚结结果。在声学方面,声强提供颗粒运动所需的能量[12,13],频率决定颗粒相对于惯性的携带动量效果[14,15,5,6],波型则定义了空间力场[9,16]。对于颗粒而言,其类型影响碰撞与黏附动力学[5,3,12],粒径分布通过异速效应驱动相对运动[17,18],数浓度则设定了基准碰撞概率[19,15,18]。就流体介质而言,温度同时影响粘度和潜在凝结现象[11],从而改变粘性阻力[20,21]及颗粒间结合力。这些因素的非线性相互作用影响了声压梯度力(SPGF)和Stokes阻力等初级力,以及声学尾流效应和相互辐射压力效应等次级机制,最终决定聚结效率。
鉴于上述参数的复杂性与敏感性,研究人员通常依赖物理试验前的数值仿真。AA研究领域可采用并广泛使用多种数值范式,每种方法都有其假设、优势和局限。离散元方法(DEM)将单个颗粒视为受压力梯度力、流体动力力和体力作用的显式实体[22,23]。DEM可表征碰撞解析、滚动/滑动摩擦及详细轨迹,因此适用于颗粒级相互作用至关重要的研究。然而,DEM的高保真度以高昂的计算成本为代价。使用小时间步长解析大量颗粒并评估详细的力模型往往变得极其昂贵,特别是在声学驱动和耦合流场需要在声学时间尺度上更新时。
另一方面,有限元方法(FEM)提供了一种基于连续介质的框架,其中控制声学与流体的方程在计算域上离散化,并通过基于网格的数值近似求解[24]。该方法特别擅长以高空间保真度解析耦合场,如声压、流体速度和颗粒输运。然而,实现这种保真度通常需要精细的空间网格和小的时间增量,这可能导致可观的计算开销和内存消耗,尤其对于高频声学激励和长时间瞬态仿真[25]。
传统数值方法以计算时间换取物理细节,一种创新替代方案是将机器学习引入仿真流程,以加速最耗时的组件,同时保留其余模型的物理结构。物理信息神经网络(PINNs)提供了一种利用物理控制方程训练模型的框架。神经网络被训练以近似从相关输入到目标量的映射,而非在每个时间步通过直接数值求解计算特定量。训练完成后,PINN通过单次前向评估即可输出该量,其成本比通过复杂计算和大量自由度对每个颗粒反复进行解析/数值求解低数个数量级。如今,PINN及相关深度学习求解器在非线性PDE和孤子承载系统方面进展迅速。Yu等人开发了定制库策略(PhyCuLi),从稀疏波数据中发现控制NLS型方程并重构非线性波动力学[26]。Wang等人提出了IRON单键Net,将有限差分数值损失与聚类重建相结合,以改善Manakov方程的多参数识别[27]。Yang等人采用神经算子网络预测(2+1)维玻色-爱因斯坦凝聚体中的亮孤子和暗孤子[28]。Jia等人进一步引入了参数驱动的Transformer-时间网络,用于从激光腔参数端到端预测超快孤子演化[29]。这些进展展示了PINN类和算子学习架构在高频和非线性场动力学方面的通用性;本工作将该范式专门应用于通过力分解的FF-PINN进行颗粒分辨声学驱动。
尽管闭式SPGF和Stokes阻力表达式对孤立颗粒状态评估来说开销不大,但AA仿真需要在每个时间步对每个颗粒反复评估。对于包含数百万颗粒并经过大量时间步积分的群体,这意味着每次仿真需要数十亿次力评估。因此,FF-PINN采用离线训练和在线推理策略:解析公式一次性生成训练标签,而运行时力预测通过批量神经前向传递完成,利用GPU/CPU张量并行。
本文介绍了傅里叶特征物理信息神经网络架构(FF-PINN),用于预测SPGF和Stokes力作用下的颗粒运动与碰撞,其中局部声压梯度与流体振荡耦合。FF-PINN在保持物理可解释性的同时,在高空间或时间分辨率仿真中实现了显著加速。如图1所示,本文后续章节记录了具体的FF-PINN架构、训练过程、与FEM/DEM数据的验证以及初步参数研究。通过将FF-PINN框架应用于从复杂动力学方程中分解的元素训练,可在保持预测精度的同时显著提高仿真效率。

节选摘要

基本声波方程

在驻波条件下,一维传播的声波位移可表示为关于位置与时间的函数[10]χ(x,t)=?Acos(kx)sin(ωt)其中χ是声场中流体分子的位移,x是坐标,A是位移振幅,k是波数,ω是角频率,t是时间。
位移振幅和波数定义为A=p0ρcωk=2πλ其中p0是声压幅值,ρ是流体密度

结果与讨论

为评估FF-PINN在声学场中预测颗粒输运的可靠性,将颗粒数浓度(PNC)的相对变化与传统DEM仿真和FEM结果进行比较。除非另有说明,本节结果均使用表1中所示配置的驻波场获得。

结论

本工作开发了一种傅里叶特征物理信息神经网络框架,通过替代反复评估声压梯度力和Stokes阻力来加速颗粒分辨AA仿真。通过将声学驱动分解为空间、时间和速度依赖分量,所提出的FF-PINN在保留控制物理规律的同时显著简化了学习任务。引入傅里叶特征编码进一步实现了精确的

代码可用性

本工作提出的FF-PINN框架使用Python和PyTorch深度学习库实现。源代码和代表性算例公开托管于FF-PINN的GitHub仓库,可通过以下链接访问: https://github.com/YinHang-jpg/FF-PINN-for-AA。该仓库包含神经网络实现、训练脚本、推理流程以及重现本文报告的主要数值结果所需的示例配置。

CRediT作者贡献声明

Hang Yin:概念化、方法论、数据整理、正式分析、初稿撰写。Pengzhan Liu:审阅与编辑。Bing Feng Ng:概念化、指导、资金获取、审阅与编辑。
Hang Yin|Pengzhan Liu|Bing Feng Ng
相关新闻
生物通微信公众号
微信
新浪微博
  • 搜索
  • 国际
  • 国内
  • 人物
  • 产业
  • 热点
  • 科普

热点排行

    今日动态 | 人才市场 | 新技术专栏 | 中国科学人 | 云展台 | BioHot | 云讲堂直播 | 会展中心 | 特价专栏 | 技术快讯 | 免费试用

    版权所有 生物通

    Copyright© eBiotrade.com, All Rights Reserved

    联系信箱:

    粤ICP备09063491号