研究反应流中电磁相互作用的计算方法

《Computer Physics Communications》：A computational approach for the study of electromagnetic interactions in reacting flows

　　

研究人员提出了一种用于模拟可压缩反应流中电磁相互作用的计算流体力学方法。开发的代码，命名为电磁积分器（EMI），基于SENGA直接数值模拟（DNS）软件。静电场和静磁场使用麦克斯韦方程组的高斯定律求解。电磁波传播通过使用显式时域有限差分法（FDTD）离散化安培定律和法拉第定律求解。电磁场方程与纳维-斯托克斯方程完全耦合，因此电磁场与流体之间的相互作用包含在公式中。相互作用项包括洛伦兹力、极化力和磁化力。这些力决定了影响动量输运、物种输运中的扩散速度以及能量守恒方程的体积力。此外，介质的性质通过电介电常数和电导率、电荷密度以及磁导率影响电磁波的传播。电磁场的解通过解析解和数值解进行了验证。静电场与物种、能量和动量守恒方程之间耦合的实现通过与已建立的代码获得的层流反应流数值解进行比较进行了验证。通过一系列层流甲烷-空气在静电场、静磁场和高频电磁场下的计算，研究了所提出公式的关键能力。还评估了在存在由带电物种运动产生的电流时静电公式的有效性。结果表明，EMI-SENGA能够捕捉电磁场对反应流的基本影响，以及带电物种的动力学及其对火焰形状和反应活性的影响。本研究提出的公式为反应流中的多物理场相互作用建模提供了一个全面的框架，这对利用电磁场操纵物种扩散和反应活性的新兴研发领域具有意义。

**论文解读：基于EMI-SENGA的反应流电磁相互作用多物理场数值框架**

**研究背景与问题**
在工程、物理和化学领域，多相反应流中的多物理场现象，特别是不同物种与外部电磁场之间的相互作用，是一个长期的研究课题。数值模拟是研究这些复杂多物理场现象的重要工具，尤其适用于难以或无法通过实验捕捉的小尺度过程。大量实验表明，外部电场和磁场能够影响反应流的行为，但现有数值研究通常做出简化假设，主要为了降低计算成本，却忽略了若干关键物理过程。这些问题包括：第一，研究通常假设静电场或静磁场，忽略了带电物种运动产生的电磁场；第二，通常忽略物种电矩或磁矩对电磁场分布的影响，且静磁场仅被简单施加而非基于边界条件和流动演化求解；第三，据作者所知，尚无数值研究探讨电磁波对反应流的影响，而实验表明交流电流和脉冲放电对反应流具有独特效应；第四，现有研究主要集中于一维或二维构型，三维模拟因计算成本高昂而极少开展。因此，开发一个全面包含所有物种电磁性质及其与电磁场相互作用的综合数值框架，并能够模拟电磁波在空间和时间上演化的计算模型，具有重要意义。本研究旨在提出一个统一的理论和数值框架，描述电磁场与反应流之间的相互作用，并开发高效代码以研究层流和湍流反应流中的电磁场效应，同时扩展框架以描述空间和时间变化的电磁波传播，并考察多种多物理现象对系统整体电动力学的贡献。

**主要关键技术方法**
研究人员基于SENGA直接数值模拟（DNS）框架，开发了电磁积分器（EMI）代码。静电场和静磁场通过求解麦克斯韦方程组的高斯定律获得，电磁波传播则采用显式时域有限差分法（FDTD）离散安培定律和法拉第定律。电磁场方程与纳维-斯托克斯方程完全耦合，纳入洛伦兹力、极化力和磁化力作为体积力，影响动量输运、物种扩散速度和能量守恒。介质性质（电介电常数、电导率、磁导率、电荷密度）通过本构关系与电磁场耦合。代码采用结构化网格，并行策略基于区域分解，并通过解析解及与Cantera、gprMax等代码的比较进行验证。

**研究结果**
**3.1 静电场求解器验证**：通过模拟化学计量甲烷-空气层流预混火焰，将EMI-SENGA结果与Cantera对比，验证了在考虑离子物种和静电相互作用时，代码对温度、密度、主要物种及离子物种质量分数的预测准确，静电场和力的实现正确。

**3.2 一维层流反应流在外部静电场下**：施加不同强度的外部电势，发现电荷密度分布受电场方向强烈影响。正电势下，电子向阳极（入口）运动并与氧气反应生成氧离子（O₂^?），导致入口附近电荷密度增大；负电势下则阳离子在反应区附近积累。洛伦兹力对带电物种产生显著加速，而极化力对中性物种的作用约小24个数量级，但因其质量分数高，可能仍影响反应流。

**3.3 一维层流反应流在均匀静磁场下**：使用不含离子物种的化学机制，施加不同磁势，发现磁场仅在火焰前锋附近因磁化率变化有微小改变（0.6 mA/cm），磁化力在二维构型中可能为零，均匀磁场对火焰形状和温度影响可忽略。但磁化力的两个贡献项（磁化率梯度和磁场梯度项）均不可忽略，且二者相互抵消导致合力为零。

**3.4 二维层流反应流在非均匀静磁场下**：模拟由圆柱磁体产生的非均匀磁场，发现磁化力对顺磁物种（如O₂）和抗磁物种（如CH₄）作用方向相反，顺磁力约强两个数量级。磁化力使反应区弯曲，温度在y方向变化2–26 K，表明边界磁势的分布直接影响火焰形状，准确边界条件至关重要。

**3.5 EMI-FDTD求解器验证**：对赫兹偶极子源进行三维FDTD模拟，与解析解及gprMax代码对比，EMI-FDTD在多种CPML参数和边界厚度下均获得一致结果，并行计算精度与单处理器一致。

**3.6 反应流中静电假设的评估**：利用EMI-FDTD求解无外场下的火焰，考虑感应电流，发现静电场解与FDTD解存在差异。感应电流导致电场分量E_x在反应区附近减小，电荷密度变化，且该变化主要由电流密度引起而非物种质量分数变化。电子迁移率（0.2 m²/(V·s) vs 1 cm²/(V·s)）对电导率影响显著，高迁移率下电流不可忽略，静电假设失效。

**3.7 电磁波源对反应流的影响**：施加100 GHz正弦波源，发现电磁波在域内产生波峰和波谷，电子所受洛伦兹力可达10¹⁰–10¹² N/kg，与静电场下相当，表明电磁波可局部影响带电物种分布，使电磁效应具有时空依赖性。

**讨论与结论**
本研究开发了EMI-SENGA代码，可模拟静、动态电磁场与反应流的多物理场耦合。验证表明代码准确。静电场下，电荷密度和电场在入口处受化学-扩散耦合影响显著；静磁场下，磁化力正确描述需包含磁化率梯度项和磁场梯度项，强非均匀磁场可改变火焰形状；FDTD与静电场的对比显示，高迁移率带电物种产生的感应电流不可忽略，静电假设在此条件下失效；电磁波源可在短时间尺度内产生局部力，影响电荷分布。EMI-SENGA为研究电磁场对反应流的基本效应及多物理场相互作用提供了全面框架，对利用电磁场调控物种扩散和反应活性的新兴领域具有重要意义。