[论文解读] Modeling coupled electrochemical and mechanical behavior of soft ionic materials and ionotronic devices

Nikola Bosnjak, Max Tepermeister+1 more|arXiv (Cornell University)|Apr 26, 2022

Advanced Sensor and Energy Harvesting Materials参考文献 103被引用 15

一句话总结

本文提出了一种多物理场计算框架，通过将泊松-能斯特-普朗克方程与有限元法相结合，模拟软离子材料和离子电子器件中电化学与机械行为的耦合。关键贡献在于一种经过敏感性分析验证的鲁棒、无量纲化有限元公式，该公式表现出良好的收敛性且对网格分辨率和纵横比不敏感，从而能够在真实条件下准确模拟复杂离子电子器件的行为。

ABSTRACT

Recently there has been an increase in demand for soft and biocompatible electronic devices capable of withstanding large stretch. Ionically conductive polymers present a promising class of soft materials for these emerging applications due to their ability to realize charge transport across the polymer network, while preserving the desired mechanical and chemical features. As opposed to electron transfer in traditional electrical conductors, the charge transport across these polymers is achieved through ion migration. When such materials are used in combination with electrical systems, they are known as ionotronic devices. The ability to simulate device performance based on its material composition and geometry would accelerate and improve ionotronic device design. The main challenge in developing reliable simulation capabilities for ionically conductive polymers is the complex and coupled electro-chemo-mechanical behavior. In this work we address this challenge by introducing a multiphysics framework incorporating the coupled effects of ion transport, electric fields and large deformation. The utility of the developed multiphysics model is showcased by simulating representative ion transport problems and the operation of soft ionotronic devices.

研究动机与目标

• 解决缺乏能够捕捉电化学与机械响应复杂相互作用的综合性离子电子器件仿真工具的问题。
• 克服等效电路模型的局限性，后者无法将材料组成、器件几何结构与电化学行为关联起来。
• 开发一种连续介质层面的多物理场模型，整合电化学、力学与离子输运，用于大变形下的软离子聚合物。
• 基于材料参数和几何结构，实现对离子电子器件（如传感器、执行器和能量收集器）的预测性仿真。
• 通过系统性地分析几何特征与网格分辨率的敏感性，验证模型的鲁棒性。

提出的方法

• 构建一种将泊松-能斯特-普朗克（PNP）方程与大变形力学耦合的无量纲化有限元法（FEM）框架。
• 采用伽辽金法对电势、离子浓度和机械位移的控制方程推导弱形式。
• 推导单元级别的残差和切线矩阵（例如，KΦΦ、KM(i)M(i)），用于牛顿-拉夫森法求解非线性系统。
• 通过特征长度（L）、德拜长度（λD）和时间（τ）进行无量纲化，以改善数值条件和收敛性。
• 采用包含1000个单元（1 nm × 1 µm × 1 µm）的三维有限元网格进行仿真，施加对称性和零通量边界条件。
• 在τ = 8个无量纲时间单位内进行平衡化步骤，以评估稳态行为。

实验结果

研究问题

• RQ1 长径比与德拜长度之比（L/λD）如何影响软离子结中的电势分布？
• RQ2 离子导体几何形状的纵横比（L/W）对电化学响应的敏感性如何？
• RQ3 截面内有限元数量在多大程度上影响电势预测的准确性？
• RQ4 该模型能否准确捕捉软离子材料在大变形下的耦合电-化学-机械响应？
• RQ5 该模型的解在不同网格密度和几何构型下是否保持收敛与稳定？

主要发现

• 当L/λD增加至80时，电势达到饱和，此后进一步增加L/λD带来的电势增量不足5%，表明L/λD = 100已足够实现准确模拟。
• 在L/W为0.1、1和10时，归一化电势与长度比无关，表明宽而薄的几何结构可采用一维近似。
• 在截面中仅使用一个有限元所得到的电势分布与使用4个或9个单元的结果无法区分，证实了网格收敛性及简化建模的有效性。
• 模型的解对网格细化和几何变化不敏感，表现出鲁棒性，适用于真实离子电子器件的仿真。
• 无量纲化FEM公式可实现稳定且准确的牛顿-拉夫森收敛，验证了其在复杂离子电子器件仿真中的适用性。
• 该框架成功捕捉了离子迁移、电场与机械变形之间的耦合效应，实现了对离子电子系统的预测性仿真。