MATLAB 기반 MCPlas 툴박스로 COMSOL 플라즈마 모델링 재현성 확보

초록

MCPlas는 MATLAB 함수 모음으로, JSON 기반 입력 데이터를 이용해 COMSOL에서 비열 플라즈마의 방정식 기반 유체‑포아송 모델을 자동으로 생성한다. 전자 수송을 세 가지 방식으로 제공하고 1D·2D(직교·극좌표·원통) 기하를 지원한다. Argon 저압 방전(DC·RF) 사례를 통해 COMSOL Plasma Module과 결과를 비교 검증했으며, 복잡한 반응 메커니즘을 JSON으로 관리해 다른 시뮬레이션 툴에서도 재사용 가능하도록 설계하였다.

상세 분석

본 논문은 저온 비열 플라즈마 모델링에서 재현성과 투명성을 확보하기 위해 MCPlas라는 MATLAB‑COMSOL 연동 툴박스를 제안한다. 핵심은 모든 모델 정의 정보를 JSON 스키마에 맞춰 구조화함으로써 데이터 검증과 교환을 자동화한다는 점이다. 기존 상용 툴은 방정식과 경계조건이 블랙박스로 숨겨져 있어 사용자 맞춤형 수정이 어렵고, 복잡한 반응 메커니즘(RKM)을 수작업으로 입력해야 하는 비효율성이 있었다. MCPlas는 이러한 문제를 다음과 같이 해결한다.

1. JSON 기반 입력: 종족·상태·전하·반응계수·전송계수 등을 JSON 파일에 기술하고, LXCat 및 Plasma‑MDS 스키마와 호환한다. 이는 FAIR 원칙을 충족시키며, 다른 시뮬레이션 플랫폼(예: PLASIMO, CFD‑ACE)에서도 동일 파일을 재사용할 수 있게 한다.

2. 자동 모델 생성: LiveLink™ for MATLAB 모듈을 이용해 JSON을 파싱하고, COMSOL General Equation Module에 직접 방정식, 경계조건, 초기·수치 파라미터를 삽입한다. 사용자는 MATLAB 스크립트만 수정하면 복잡한 2D·극좌표·원통 좌표 모델을 즉시 구축할 수 있다.

3. 전자 수송 모델의 다양성: 기존의 전형적인 Drift‑Diffusion Approximation (DDAc) 외에 DDA53(간소화된 에너지 플럭스)와 새롭게 도입한 DDAn(Legendre 다항식 기반 고차 모멘트) 세 가지 옵션을 제공한다. DDAn은 전자 충돌 주파수와 고차 전송 계수를 포함해 저압·대기압 모두에서 정확도가 향상된다는 기존 연구(Refs. 23, 38, 40)를 구현한 최초의 툴이다.

4. 경계조건 구현: Hagelaar et al. (2015)의 전자 및 중성 입자 경계조건을 그대로 적용하고, 전자 반사계수·2차 전자 방출·이온 흐름 등을 사용자 정의 가능하게 했다. 이는 플라즈마‑전극 상호작용을 정밀히 모델링하는 데 필수적이다.

5. 검증 사례: Argon 저압 방전(DC 및 RF) 두 가지 기준 사례를 선정해 MCPlas와 COMSOL Plasma Module(CPM)의 전자 밀도, 전위, 전류 파형 등을 비교하였다. 결과는 전자 수송 모델 선택에 따라 차이가 크게 나타났으며, 특히 DDAn을 적용했을 때 실험 데이터와 가장 근접한 전압‑전류 특성을 보였다. 이는 전송 모델의 물리적 정확도가 시뮬레이션 결과에 미치는 영향을 명확히 보여준다.

6. 재사용성 시연: 동일 JSON 파일을 PLASIMO와 같은 다른 플라즈마 시뮬레이션 툴에 입력해 동일한 방전 조건을 재현함으로써 데이터 교환의 실용성을 입증했다.

전반적으로 MCPlas는 모델 정의의 표준화, 자동화된 방정식 생성, 고급 전자 수송 옵션 제공을 통해 플라즈마 시뮬레이션의 투명성·재현성을 크게 향상시킨다. 또한 오픈소스 MIT 라이선스로 배포돼 연구 커뮤니티가 자유롭게 확장·수정할 수 있다.