스크립트 언어 기반 연속체 다중물리 시뮬레이션 컴파일러 Nsim 프로토타입

초록

본 논문은 고전장 이론으로 기술되는 연속체 물리 시스템을 대상으로, 물리 방정식을 런타임에 자동으로 병렬화된 유한요소 코드로 변환하는 프레임워크를 제시한다. 파이썬 등 스크립트 언어에 다중물리 확장 모듈을 구현함으로써 사용자는 순수 스크립트 수준에서 시스템 정의, 경계조건, 물성치 등을 기술하고, 컴파일러가 내부적으로 효율적인 수치 해석 코드를 생성한다. 마그네틱 시뮬레이션 패키지 Nmag과 반응‑확산 형태소 발생 모델 두 사례를 통해 구현 가능성과 성능을 검증한다.

상세 분석

이 연구는 연속체 물리학에서 흔히 사용되는 편미분 방정식(PDE)을 고수준 스크립트 언어로 기술하고, 이를 자동으로 병렬화된 유한요소(FEM) 코드로 변환하는 ‘시뮬레이션 컴파일러’ 개념을 제시한다. 핵심 아이디어는 물리 방정식의 기호적 표현을 추출하고, 이를 변분 원리 혹은 강형(weak form) 형태로 변환한 뒤, 메쉬와 자유도(DOF) 정보를 결합해 행렬 연산을 생성하는 과정이다. Nsim은 이러한 파이프라인을 파이썬 바인딩을 통해 제공함으로써, 사용자는 파이썬 스크립트 안에서 필드 변수, 물성 파라미터, 경계조건, 초기조건 등을 선언하고, ‘solve()’ 호출만으로 시뮬레이션을 실행할 수 있다.

컴파일러는 런타임에 C++ 코드 조각을 생성하고, MPI 기반의 PETSc 라이브러리와 연동해 분산 메모리 환경에서 선형·비선형 시스템을 해결한다. 이때 자동 미분(AD) 기법을 활용해 비선형 항의 야코비안을 효율적으로 구성하고, 사용자 정의 함수도 C++ 레벨에서 JIT(Just‑In‑Time) 컴파일이 가능하도록 설계하였다. 또한, 메쉬 생성·관리, 영역별 물성치 할당, 적응형 시간 적분 등 전형적인 FEM 워크플로우를 추상화해 스크립트 수준에서 일관된 API를 제공한다.

두 사례 연구는 프레임워크의 범용성을 입증한다. 첫 번째는 마이크로자기학 시뮬레이션 Nmag으로, LLG(라르민-길버트-고든) 방정식과 마이크로스케일 전자기 상호작용을 동시에 다루며, 복잡한 비선형 경계조건과 비등방성 항을 자동으로 처리한다. 두 번째는 Turing 패턴을 생성하는 반응‑확산 방정식으로, 다중 스칼라 필드와 비선형 반응항을 포함한다. 두 경우 모두 기존 전용 코드와 비교했을 때 구현 편의성은 크게 향상되었으며, 병렬 효율성도 손실 없이 유지되었다.

이 논문의 주요 기여는 (1) 물리 방정식의 기호적 정의를 런타임에 수치 코드로 변환하는 자동화 파이프라인, (2) 파이썬과 같은 고수준 언어와 고성능 병렬 라이브러리(PETSc, MPI)의 원활한 연동, (3) 다중물리 확장성을 고려한 모듈식 설계, (4) 실제 과학·공학 문제에 적용 가능한 검증 사례 제공이다. 한계점으로는 현재 지원되는 PDE 형태가 변분 원리 기반에 국한되고, 복잡한 비선형 경계조건이나 비정형 메쉬에 대한 최적화가 추가 연구가 필요하다는 점을 들 수 있다. 향후 연구에서는 자동 미분을 통한 고차 미분 연산 지원, GPU 가속 통합, 그리고 더 일반적인 미분 연산자(예: 비보존형 흐름)까지 확장하는 것이 목표이다.