ELAS3D Xtal 고성능 결정탄성 해석기와 자동 미세구조 생성

초록

ELAS3D‑Xtal은 Fortran과 OpenMP 기반으로 재작성된 NIST ELAS3D의 확장판으로, 결정학적 이방성을 고려한 회전 강성 텐서를 사전 계산하고, 포인트‑블록 Jacobi 전처리를 적용한 행렬‑프리 PCG 솔버를 이용해 3차원 다결정 및 결함이 포함된 마이크로스트럭처의 탄성장을 고속으로 계산한다. 단일 멀티코어 PC에서 기본 CG 대비 10배, 블록 전처리 CG는 53‑61배 가속을 달성했으며, Eshelby 포함 해석으로 정확성을 검증하였다. 또한 공간 필터링 기반 통계적 다결정 생성, XCT 기반 기공 삽입, 텍스처 할당 등을 자동화하여 HDF5 형식으로 전장 데이터를 출력한다.

상세 분석

ELAS3D‑Xtal은 기존 NIST ELAS3D가 갖고 있던 한계, 즉 단일 스레드 실행과 단순 스칼라 Jacobi 전처리, 그리고 이방성 재료에 대한 비효율적인 처리 방식을 근본적으로 개선하였다. 첫 번째 핵심 개선은 OpenMP를 이용한 공유‑메모리 병렬화이다. 데이터 구조를 3‑D 정규 격자로 유지하면서, 이웃 인덱스 테이블(ib)을 미리 구축해 메모리 접근을 연속적으로 만들고, 루프 차단(blocking) 기법을 적용해 L1/L2 캐시 활용도를 높였다. 이는 메모리 대역폭이 제한적인 현대 CPU에서 행렬‑벡터 곱을 효율적으로 수행하게 한다.

두 번째 핵심은 포인트‑블록 Jacobi 전처리이다. 각 격자점에서 3개의 변위 자유도를 3×3 블록으로 묶어 전처리 행렬 M⁻¹을 미리 계산한다. 이는 강성 텐서의 비대각 성분이 큰 결정 이방성 재료에서 스칼라 전처리보다 훨씬 나은 수렴성을 제공한다. 전처리 행렬은 각 격자점 주변 8개의 요소(27‑점 스텐실)에서 얻은 로컬 강성 서브블록을 합산해 역행렬을 구함으로써 얻는다.

세 번째는 회전 강성 텐서의 사전 계산이다. 각 결정 입자마다 Rodrigues 벡터로 정의된 오리엔테이션을 3×3 회전 행렬 R로 변환하고, 이를 6×6 변환 행렬 Q로 확장해 C_global = Q C_local Qᵀ 형태로 전역 강성을 얻는다. 이 과정은 초기화 단계에서 한 번만 수행되므로 런타임 오버헤드가 거의 없다.

네 번째는 마이크로스트럭처 자동 생성 파이프라인이다. 공간 필터링을 이용해 통계적으로 캘리브레이션된 입도와 종횡비를 갖는 다결정을 생성하고, 병렬 타원형 Voronoi Tessellation을 통해 입자 할당을 수행한다. 또한 XCT에서 추출한 기공 중심·반경 정보를 직접 삽입해 실제 결함 형태를 재현한다. 이러한 전처리 단계는 모두 Fortran 모듈로 구현되어 HDF5 파일에 저장되며, 후속 ELAS3D‑Xtal 솔버와 원활히 연계된다.

성능 평가에서는 100³500³ 격자(10⁶1.25×10⁸ 셀)에서 블록 전처리 PCG가 순차 CG 대비 53‑61배 가속을 보였으며, 800³ 격자(5.12×10⁸ 셀)까지 확장 가능함을 시연했다. 정확도는 3‑D Eshelby 포함 해석과 비교해 오차가 1 % 이하로 유지되었다. 마지막으로, LPBF SS316L 합금에 대한 실제 기공 형태(가스, L‑of‑F, 키홀)를 적용한 사례에서, 결정 이방성에 따른 응력 집중 차이를 정량적으로 분석함으로써 AM 결함 메카니즘 연구에 바로 활용할 수 있음을 입증했다.

이러한 설계와 구현은 고해상도 다결정·결함 시뮬레이션을 단일 워크스테이션 수준에서 수행할 수 있게 하여, 연구자와 엔지니어가 복잡한 재료 거동을 빠르게 탐색하고 설계 최적화에 적용할 수 있는 기반을 제공한다.