Coupled Infrared Imaging and Multiphysics Modeling to Predict Three-Dimensional Thermal Characteristics during Selective Laser Melting

기술 요약: 선택적 레이저 용융(SLM)에서의 3D 열적 특성화를 위한 결합형 적외선 이미징 및 다물리 모델링

문제 정의
적층 제조(AM), 특히 선택적 레이저 용융(SLM) 과정 중 발생하는 레이저 가열은 미세조직 진화와 기계적 성질을 결정짓는 극심하고 일시적인 열적 조건을 유도한다. 이 분야의 핵심 과제는 충분한 시공간적 정확도로 이러한 조건들을 해상하는 것이다. 피로미터(pyrometry)나 이파장 온도계(two-wavelength thermometry)와 같은 기존 실험 기법들은 정확한 열 구배를 포착하기 위한 공간 해상도가 부족하거나 표면 측정에 국한되는 경우가 많다. 고속 적외선(IR) 이미징은 더 나은 해상도를 제공하지만, 고융점 합금의 경우 방사 신호를 실제 온도로 변환할 때 방사율 불확실성과 좁은 스펙트럼 대역에 적용된 슈테판-볼츠만 법칙의 한계로 인해 어려움이 따른다. 반대로, 고충실도 수치 시뮬레이션은 3차원(3D) 데이터를 제공할 수 있지만, 계산 비용이 많이 들고 시간이 오래 걸리며, 알려지지 않은 열물리적 특성으로 인해 실험 조건과 일치시키기 위한 인위적인 파라미터 조정이 필요한 경우가 많다. 또한, X선 투과 촬영과 같은 기술은 단일 단면 평면에 제한된다. 결과적으로, 고체-액체(S-L) 계면에서의 실험적으로 검증된 3D 일시적 응고 조건(열 구배 및 응고 속도)을 얻는 데 공백이 존재한다.

방법론
저자들은 동적인 하부 표면 온도 분포를 재구성하기 위해 통합된 실험-계산 프레임워크를 개발하였다. 방법론은 다음과 같다:

1. 고속 IR 이미징: 고속 IR 카메라(>10,000 fps)를 사용하여 워롯 니켈 기반 초합금(MAR-M247)에 대한 실험을 수행하여 과도적인 용융지 표면 복사를 포착하였다.
2. 커스텀 캘리브레이션: 광학 투과율, 센서 노출 시간, 대역 내 복사 강도를 포함한 특정 실험 파라미터를 플랑크 법칙(3–5.5 µm 범위에 대해 적분됨)과 통합한 커스텀 캘리브레이션을 생성하여, 복사 에너지를 표면 온도 맵으로 변환하였다.
3. 결합형 다물리 모델링: 실험적으로 도출된 시계열 2D 표면 온도 맵을 3D 열전도 및 고체-액체 상변화 모델(COMSOL)의 일시적 경계 조건으로 도입하였다. 이 접근 방식은 레이저 흡수, 마랑고니 흐름(Marangoni flow), 표면 복사와 같은 복잡한 현상을 명시적으로 모델링할 필요를 없애는데, 이는 해당 현상들이 측정된 표면 온도에 의해 본질적으로 포착되기 때문이다.
4. 검증 및 분석: 모델의 용융지 기하학적 구조 및 예측값은 ex situ 주사 전자 현미경(SEM) 및 전자 후방 산란 회절(EBSD) 맵을 통해 검증되었다. S-L 계면에서의 응고 속도($R$)와 열 구배($\nabla T$)를 계산하였으며, 이를 Trivedi 모델과 연관시켜 셀 간격을 예측함으로써 미세조직 특징들과 상관관계를 분석하였다.

주요 결과

• 3D 용융지 재구성: 결합된 접근 방식은 3D 용융지 프로파일과 하부 표면 온도 분포를 성공적으로 재구성하였다. 모델이 예측한 용융지 깊이와 단면 프로파일은 다양한 레이저 공정 파라미터(예: 200 W @ 500 mm/s)에 대해 ex situ SEM 측정값과 높은 일치도를 보였다.
• 물성 변화에 대한 강건성: 민감도 분석 결과, 입력 열물리적 특성이 10% 변동하더라도 예측된 용융지 깊이는 1.5 µm 미만으로 변동하여, 특성이 잘 규명되지 않은 재료에 대한 이 방법의 유용성을 입증하였다.
• 유체 흐름 및 마랑고니 효과: 열적 마랑고니 흐름이 내부 온도 분포를 약간 균질화하는 것으로 나타났으나, 본 연구에서는 전도(conduction)가 용융지 내부에서 지배적인 열 전달 메커니즘임을 확인하였다(평균 페클레 수 < 1). 결과적으로, 실험적 표면 데이터에 의해 구동되는 단순 전도 모델만으로도 전체 유체 역학 모델링의 계산 비용 없이 필수적인 응고 조건을 포착하기에 충분하였다.
• 응고 역학: 연구 결과, 응고 속도($R$)의 상당한 공간적, 시간적 변화가 밝혀졌다. $R$은 용융지 가장자리에서의 $10^{-3}$ m/s부터 꼬리 부분에서의 $\sim 100$ m/s까지 3 자릿수 차이를 보였다. 특히, 표면 근처의 최대 응고 속도는 레이저 스캔 속도보다 높게 나타났는데($R_{max} \approx 1.5 \times V_{laser}$), 이는 단순한 기하학적 근사치가 $R$을 스캔 속도로 제한한다는 점과 상반되는 현상으로, 강한 대류 효과와 국부적인 열 구배에 기인한다.
• 미세조직 예측: 계산된 열 구배와 응고 속도는 Trivedi 모델을 통해 셀 간격을 예측하는 데 사용되었다. 예측된 셀 간격 분포(500 mm/s에서 중앙값 $\sim 0.7$ µm, 1000 mm/s에서 $\sim 0.52$ µm)는 실험적인 EBSD 및 SEM 측정값과 밀접하게 일치하였다. 모델은 스캔 속도 증가에 따라 셀 간격이 감소하는 경향을 성공적으로 포착하였다.
• 시간적 가변성: 열 구배와 셀 방향의 표준 편차 분석을 통해, 열적 조건이 응고 시작 단계에서는 가장 안정적이지만 용융지의 꼬리 부분으로 갈수록 상당한 가변성을 보인다는 것을 확인하였다. 이 높은 가변성 영역은 미세조직 전이(예: 주상정-등축정 전이) 및 결함 형성이 발생하기 쉬운 영역과 일치한다.

의의 및 주장
본 논문은 이 통합된 실험-계산 접근 방식이 상업적 AM에서 미세조직의 in situ 예측 및 최적화를 실현하는 데 필수적이라고 주장한다. 고속 IR 이미징을 경계 조건으로 활용함으로써, 이 방법은 순수 수치 시뮬레이션보다 훨씬 적은 계산 비용과 적은 가정을 사용하여 고충실도 3D 열적 특성화를 달anche한다. 저자들은 이 프레임워크가 다음과 같은 데이터 풍부한 플랫폼을 제공한다고 확언한다:

1. 미세조직 모델(예: 셀 간격에 대한 재료 의존적 상수 결정)의 검증 및 정교화.
2. 공정 파라미터에 의해 구동되는 미세조직 크기 및 방향 예측.
3. 열적 조건의 시간적 변화에 기반한 결함 취약 영역 식의 식별.
4. 열물리적 특성이 불확실한 새로운 재료 및 제조 시나리오(예: 분말 원료)로의 확장성(이 방법은 가정된 재료 파라미터가 아닌 실험적 경계 조건에 의존하기 때문).

본 연구는 현재의 최첨단 기술들이 2D 정보를 제공하거나 값비싼 3D 시뮬레이션을 요구하는 반면, 이 결합된 접근 방식이 실험적으로 검증된 시공간적 해상도의 3D 응고 조건을 제공함으로써 그 간극을 메워, AM 공정 제어에 필수적인 역할을 한다고 결론짓는다.