[論文レビュー] 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

ひとこと要約

本稿では、ポisson＝ナーンスト＝プランク方程式を有限要素法と統合することで、軟らかくイオン性の材料およびイオントロニクスデバイスの電気化学的および機械的挙動を連成的にモデル化するマルチフィジックス計算フレームワークを提示する。主な貢献は、メッシュ解像度およびアスペクト比に対して収束性と感度の低さを示す、安定した次元なし有限要素法の定式化であり、実際の条件下での複雑なイオントロニクスデバイス挙動の正確なシミュレーションを可能にする。

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）を有する3次元有限要素メッシュを用い、対称性およびフラックスゼロの境界条件を適用する。
• ステディステート挙動の評価のため、τ = 8の次元なし時間単位にわたり平衡化ステップを適用する。

実験結果

リサーチクエスチョン

• RQ1 長さ対デバイー長比（L/λD）が、軟らかくイオン性の接合部における電位分布にどのように影響を与えるか？
• RQ2 イオン導体の幾何形状のアスペクト比（L/W）に起因する電気化学的応答への感度はどの程度か？
• RQ3 断面における有限要素数が、電位予測の精度にどの程度影響を及えるか？
• RQ4 本モデルは、大変形下における軟らかくイオン性材料の電気化学的・機械的連成応答を正確に捉えることができるか？
• RQ5 メッシュ密度および幾何形状の変化にかかわらず、モデルの解は収束性を示し、安定性を保っているか？

主な発見

• 電位はL/λDが80に達するまで増加し、以降は追加の電位増加が5％未満に抑えられ、飽和状態に達することが示された。したがって、L/λD = 100で十分な精度のシミュレーションが可能である。
• L/W比が0.1、1、10のいずれであっても、正規化された電位は一定であり、広く薄い幾何形状では1次元近似が成立することが示された。
• 断面に1つの有限要素を用いた場合と4つまたは9つを用いた場合の電位プロファイルは区別がつかず、メッシュ収束および簡略化モデルの妥当性が確認された。
• モデルの解は、メッシュの細分化および幾何的変化に対して頑健で、実世界のイオントロニクスデバイスのシミュレーションに信頼性を示す。
• 次元なしFEM定式化により、安定的かつ正確なニュートン・ラプソン収束が実現され、複雑なイオントロニクスデバイスのシミュレーションに応用可能であることが検証された。
• 本フレームワークは、イオン移動、電場、機械的変形の連成効果を正確に捉えることができ、イオントロニクスシステムの予測的シミュレーションを可能にした。