KTP 결정 내 연속파 레이저 가열 해석을 위한 오픈소스 열 모델링 툴킷

초록

본 연구는 포타슘 티타닐 포스페이트(KTP) 결정 내에서 연속파(CW) 가우시안 빔에 의한 2차 고조파 발생(SHG) 시 열 분포를 계산하는 오픈소스 유한 차분 툴킷을 공개한다. 이 툴킷은 결정 기하학 및 물성 정의, 경계 및 냉각 조건 모델, 과도 및 정상 상태 솔버를 포함하며, 사용자가 제공한 파라미터를 기반으로 시공간적 온도장을 계산해 출력한다. 기존 연구를 통합한 이 구현체는 온도 의존적 열전도율, 복사를 포함한 대류 냉각, 다양한 열전달 계수 등 여러 시나리오를 재현 가능한 파이프라인으로 검증하였고, MIT 라이선스로 공개되었다.

상세 분석

본 논문에서 소개된 툴킷의 기술적 핵심은 KTP 결정의 비등방성 열방정식을 원통 좌표계(r, z)에서 유한 차분법(FDM)으로 수치 해석하는 것이다. 가장 중요한 기여는 단순화된 상수 열전도율(K) 모델을 넘어, 실제 KTP의 열전도율이 온도에 따라 변하는 K(T) 모델을 구현했다는 점이다. 검증 결과에 따르면 80W 펌프 조건에서 정상 상태 중심 온도는 상수 K 모델 대비 K(T) 모델에서 약 70K 높게 계산되어, 고출력 SHG 설계 시 정확한 열 해석의 중요성을 보여준다.

두 번째 주요 통찰은 경계 조건 모델링의 정교함에 있다. 툴킷은 단순한 고정 온도 경계 조건이 아닌, 로빈(Robin) 경계 조건(대류 + 선택적 복사)을 구현한다. 복사 열전달은 표면적이 큰 경우(예: 반경 5mm, 스폿 3mm)에만 유의미하게 작용하여 입출력면 온도를 각각 약 15K, 10K 낮추는 것으로 분석되었다. 또한, 냉각 계수(h)에 대한 민감도 분석을 통해 공기 중 자연 대류(~6.5–50 W/m²K) 영역에서는 h 변화에 무관하지만, 고성능 마운트 또는 물 냉각(~15,000–20,000 W/m²K)을 적용하면 표면 온도가 약 20K 낮아지는 것을 확인하였다.

아키텍처 측면에서, 이 툴킷은 포트란(Fortran)으로 작성된 고성능 계산 커널을 기반으로 하며, 기하학/물성 정의, 경계 조건, 솔버 모듈로 구성되어 있다. 사용자는 펌프 출력, 빔 스폿 크기, 흡수 계수, 결정 크기, 냉각 조건 등을 입력하면, 솔버가 시간에 따른 온도 필드 T(r, z, t)와 정상 상태의 방사형/축방향 프로파일을 텍스트 파일로 출력한다. 이는 후속 SHG 효율 분석이나 열 렌즈 효과 계산을 위한 입력 데이터로 직접 활용될 수 있다.

현재 툴킷의 주요 한계는 CW 가우시안 펌핑에 국한되어 있으며, 펄스형 excitation이나 비가우시안 빔 형상은 지원하지 않는다. 또한 KTP 물성에 특화되어 있어 LBO, BBO 같은 다른 비선형 결정으로의 확장은 미래 로드맥에 포함되어 있다. 이러한 명확한 범위 설정과 검증된 베이스라인 제공은 연구자들이 복잡한 열-광학 결합 현상을 분석할 때 신뢰할 수 있는 출발점을 마련해 준다는 점에서 실용적 가치가 크다.