[论文解读] Transport of Charged Particles: Entropy Production and Maximum Dissipation Principle

Chia-Yu Hsieh, YunKyong Hyon+3 more|arXiv (Cornell University)|Jul 31, 2014

Spectroscopy and Quantum Chemical Studies参考文献 18被引用 4

一句话总结

本文提出了一种改进的泊松-能斯特-普朗克（PNP）系统，通过能量变分法引入了离子种类间相对速度差异导致的额外耗散。该改进系统与原始PNP系统具有相同的平衡态，但由于额外的熵产生，动力学行为有所不同，通过Galerkin方法和Schauder不动点定理严格证明了局部解的存在性，有限元模拟结果表明尽管平衡态相同，瞬态行为仍表现出发散特征。

ABSTRACT

In order to describe the dynamics of crowded ions (charged particles), we use an energetic variation approach to derive a modified Poisson-Nernst-Planck (PNP) system which includes an extra dissipation due to the effective velocity differences between ion species. Such a system is more complicated than the original PNP system but with the same equilibrium states. Using Schauder's fixed-point theorem, we develop a local existence theorem of classical solutions for the modified PNP system. Different dynamics (but same equilibrium states) between the original and modified PNP systems can be represented by numerical simulations using finite element method techniques.

研究动机与目标

• 为标准PNP理论因排除体积效应和高离子密度而失效的拥挤、非理想环境中的离子输运建模。
• 利用热力学一致的变分原理，将离子种类间相对运动引起的额外耗散引入PNP框架。
• 为由此产生的改进PNP系统建立严格的局部经典解理论，尤其针对增强的非线性性。
• 通过数值比较改进与原始PNP系统的动力学行为，突出显示尽管平衡态相同，其瞬态行为存在显著差异。

提出的方法

• 采用能量变分法，结合最小作用量原理（LAP）处理保守力，以及最大耗散原理（MDP）处理耗散力。
• 通过添加与阳离子与阴离子间相对速度平方成正比的项，修改耗散泛函：$ \frac{k_B T}{D_{n,p}} c_n c_p |\vec{u}_n^* - \vec{u}_p^*|^2 $，以表征种类间的阻力。
• 将改进的PNP系统推导为一组耦合的、非线性的、退化的抛物型PDE，用于描述离子浓度与静电势。
• 在合适的函数空间框架下，利用Galerkin方法和Schauder不动点定理，严格证明经典解的局部存在性。
• 采用带有隐式时间离散化的有限元法（FEM）格式求解系统，确保数值稳定性与能量定律的近似。
• 在相同边界条件与初始数据下，比较改进与原始PNP系统的数值解，以隔离额外耗散的影响。

实验结果

研究问题

• RQ1 在拥挤离子系统中，引入与相对速度相关的耗散如何影响离子输运的动力学行为？
• RQ2 带有额外种间阻力项的改进PNP系统是否能保持与原始PNP系统相同的平衡态？
• RQ3 改进PNP系统的数学结构为何？尽管非线性性增强且可能存在退化，其是否仍允许局部经典解的存在？
• RQ4 在相同平衡态与边界条件下，改进与原始PNP系统的瞬态动力学行为有何差异？
• RQ5 额外耗散项在多大程度上改变了离子输运模拟中的熵产生率与能量耗散分布？

主要发现

• 改进的PNP系统引入了一项依赖于阳离子与阴离子间相对速度的新耗散项，能够捕捉原始PNP模型中缺失的种间阻力效应。
• 尽管增加了复杂性与非线性性，改进的PNP系统与原始PNP系统具有相同的平衡态，因为二者均由相同的自由能泛函决定。
• 通过Galerkin方法与Schauder不动点定理，严格证明了改进PNP系统经典解的局部存在性定理。
• 使用有限元法进行的数值模拟证实，在相同边界条件下，改进与原始PNP系统达到相同的平衡态。
• 两套系统表现出显著不同的瞬态动力学行为，改进系统因额外耗散项导致更高的熵产生率，尤其在时间相关的耗散分布中表现明显。
• 图2中对耗散率$\triangle^*$（改进）与$\triangle$（原始）的比较表明，额外耗散改变了系统趋向平衡的弛豫路径，但未改变最终状态。