[논문 리뷰] Modeling coupled electrochemical and mechanical behavior of soft ionic materials and ionotronic devices

Nikola Bosnjak, Max Tepermeister외 1명|arXiv (Cornell University)|2022. 04. 26.

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.

연구 동기 및 목표

• 전기화학적 및 기계적 반응 간의 복잡한 상호작용을 포괄적으로 반영하지 못하는 이온트로닉 장치용 시뮬레이션 도구의 부족을 해결한다.
• 등가 회로 모델의 한계를 극복하여 재료 조성, 장치 기하학적 형상 및 전기화학적 거동 간의 연관성을 확립한다.
• 큰 변형 하에서의 소프트 이온성 고분자에 대해 전기화학, 기계학, 이온 이동을 통합한 연속체 수준의 다물리적 모델을 개발한다.
• 재료 특성과 기하 형상에 기반하여 이온트로닉 장치(예: 센서, 액추에이터, 에너지 수확 장치)의 예측 시뮬레이션을 가능하게 한다.
• 기하 특성과 메쉬 해상도에 대한 체계적 민감도 분 析를 통해 모델의 강건성을 검증한다.

제안 방법

• 대규모 변형 기계학과 결합된 포아송-너른스트-플랭크(Poisson-Nernst-Planck, PNP) 방정식을 다루는 비차원화 유한요소법(FEM) 프레임워크를 수립한다.
• 전기적 잠재력, 이온 농도, 기계적 변위에 대한 제어 방정식의 약한 형태를 갈레르킨 방법을 사용해 구현한다.
• 비선형 시스템의 뉴턴-라프슨 해법을 위해 요소 수준의 잔차와 탄성 행렬(예: KΦΦ, KM(i)M(i))을 유도한다.
• 특성 길이(L), 데바이 길이(λD), 시간(τ)를 사용해 비차원화를 적용하여 수치적 조건 개선 및 수렴성 향상에 기여한다.
• 3차원 유한요소 메쉬를 사용하여 시뮬레이션을 수행하며, 요소 수는 1000개(1 nm × 1 µm × 1 µm)로 설정하고 대칭 및 불순도가 없는 경계 조건을 적용한다.
• 정적 상태 거동 평가를 위해 총 τ = 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 형식은 안정적이고 정확한 뉴턴-라프슨 수렴을 가능하게 하여 복잡한 이온트로닉 장치 시뮬레이션에의 적용 타당성을 검증한다.
• 프레임워크는 이온 이동, 전기장 및 기계적 변형의 연성 효과를 성공적으로 포착하여 이온트로닉 시스템의 예측 시뮬레이션을 가능하게 한다.