A hybrid sharp-diffuse interface approach to accurately model melt pool dynamics with rapid evaporation in laser-based processing of metals

기술 요약: 용융 풀 역학을 위한 하이브리드 샤프-디퓨즈 인터페이스 접근법

문제 정의
레이저 기반 금속 가공(예: 레이저 빔 용접, 분말 베드 융합)에서 용융 풀 역학의 예측 시뮬레이션을 위해서는 용융체-가스 계면에서의 복잡한 열-유체 역학적 상호작용을 해결해야 한다. 증발 유도 반동 압력(recoil pressure) 및 온도 의존적 표면 장력과 같은 주요 구동력은 계면 온도에 매우 민감하게—종종 지수 함수적으로—의존한다. 표준 디퓨즈 인터페이스(DI) 정식화는 파쇄 및 합체와 같은 복잡한 위상 변화를 처리하는 데는 견고하지만, 유한한 전이 영역 전반에 걸쳐 극심한 열 구배를 뭉개뜨리는(smear) 경향이 있다. 이러한 뭉개짐 현상은 계면 온도의 예측에 결정적인 오류를 도입하며, 이는 결과적으로 계면 힘의 예측과 그에 따른 용융 풀 역학의 정확도를 저하시킨다. 반대로, 샤프 인터페이스(SI) 방법은 높은 열적 충실도를 제공하지만, 큰 밀도 비와 고도로 동적인 다상 유동을 포함하는 복잡한 수치적 복잡성 및 견고성을 확보하는 데 어려움을 겪는 경우가 많다. 기존 모델들은 대개 물리적 충실도와 계산 가능성 사이의 절충안을 강요하며, 열적 부정확성을 보상하기 위해 지나치게 미세한 격자 해상도나 광범위한 파라미터 보정을 요구하기도 한다.

방법론
저자들은 고해상도의 샤프 인터페이스 열 모델과 디퓨즈 인터페이스 유동 모델을 결로하여, 완전히 샤프 인터페이스 방식의 다상 시뮬레이션이 요구하는 계산 비용을 들이지 않고도 높은 충실도의 모델링을 달야내는 하이브리드 샤프-디퓨즈 인터페이스(HSDI) 접근법을 제안한다.

1. 샤프 인터페이스 열 모델 (CutFEM): 열 전달 문제는 샤프 인터페이스 언피티드(unfitted) 유한 요소법인 CutFEM을 사용하여 해결한다. 이 정식화는 금속-가스 계면에서의 온도 및 열 유속의 불연속성을 정확하게 해결한다. 저자들은 두 가지 변형을 제시한다:

   • SI-2P (결합된 2상 모델): 금속과 주변 가스 상 모두에서 열 방정식을 풀며, 니체 방법(Nitsche's method)을 통해 계면에서의 온도 연속성을 약하게 강제한다.
   • SI-1P (단상 모델): 금속 도메인에서만 열 방정식을 풀며, 가스 상은 부과된 열 유속을 제외하고는 단열 상태로 취급한다. 이 변형은 미해결 증발 모델에서 주변 가스의 열용량이 무시할 만하다는 점에 착안하였으며, 더 높은 계산 효율성을 제공한다.
     두 변형 모두 작은 컷 엘리먼트(cut elements)를 처리하기 위해 면 기반 고스트 페널티(face-based ghost-penalty) 안정화 기법을 사용한다.

2. 디퓨즈 인터페이스 유동 모델: 큰 밀도 비와 복잡한 계면 역학을 특징으로 하는 다상 유동 문제는 견고한 레벨 셋(level-set) 기반의 원-플루이드(one-fluid) DI 유한 요소 정식화를 사용하여 모델링된다. 이 접근법은 표면 장력과 반동 압력을 위해 연속체 표면력(CSF) 모델을 활용하여, 계면의 파쇄 및 합체를 자연스럽게 처리할 수 있게 한다.

3. 일관된 결합 전략: 두 프레임워크를 연결하기 위해, 샤프 인터페이스 온도장을 최근접 점 투영(closest-point projection) 알고리즘을 사용하여 금속-가스 계면 주변의 좁은 밴드 영역으로 확장한다. 이렇게 확장된 온도장($\tilde{T}_\Gamma$)은 디퓨즈 인터페이스 유동 프레임워크 내에서 온도 의존적 계면 힘(반동 압력 및 표면 장력)을 평가하는 데 사용된다. 이를 통해 SI 모델로부터 얻은 고정밀 열 데이터가 DI 유동 솔버를 일관되게 구동하도록 보장한다.

4. 수치 프레임워크: 구현에는 행렬 프리(matrix-free) 연산자 평가, 적응형 격자 세분화(adaptive mesh refinement), 그리고 MPI 기반 병렬화를 지원하는 deal.II 라이브러리가 사용된다. 시간 적분을 위해 분할된 연산자 분리(partitioned operator-splitting) 기법이 채택되었으며, 결합 항들을 명시적(explicitly)으로 처리하면서 개별 필드들을 암시적(implicitly)으로 해결한다.

주요 기여

• 새로운 하이브리드 정식화: 본 연구는 미해결 증발 모델을 기반으로 하는 SI CutFEM 온도 모델링과 DI FEM 유동 모델링을 결합한 최초의 2상 열-유체 역학용 HSDI 접근법을 제시한다.
• 열적 정확도: 표준 DI 모델이 실패하는 지점, 즉 용융 풀 표면에서의 온도장을 정확하게 예측할 수 있게 하는 CutFEM 기반 SI 정식화(SI-2P 및 SI-1P)를 개발하였다.
• 결합 알고리즘: 샤프 인터페이스 온도를 디퓨즈 전이 영역으로 확장하여, DI 유동 솔버 내에서 온도 의존적 소스 항을 정확하게 평가할 수 있는 일관된 전략을 제공한다.
• 벤치마킹: 결합된 열-유체 역학 모델의 정확도를 종합적으로 평가하기 위해 레이저 유도 증기 함몰(vapor depression)을 나타내는 새로운 2D 벤치마크 예시를 제안한다.

결과

• 열적 수렴성: 레이저-금속 상호작용을 모사한 순수 열 벤치마크에서, SI 모델은 2차 공간 수렴($O(h^2)$)을 보인 반면, 표준 DI 모델은 1차 수렴($O(h^1)$)만을 보였다.
• 격자 해상도 대비 정확도: 계면 온도 및 반동 압력에서 1%의 목표 정확도를 달성하기 위해, SI 모델은 DI 모델이 요구하는 것보다 두 자릿수 더 큰 유한 요소 크기를 허용했다. 구체적으로, SI 접근법은 DI 접근법보다 약 50배 더 성긴(coarser) 격자로 벤치마크 허용 오차를 달성했다.
• 결합 열-유체 역학: 새로운 2D 증기 함몰 벤치마크에서, HSDI 접근법은 동일한 격자 상에서 순수 DI 모델보다 한 자릿수 더 높은 정확도를 보였다. HSDI 접근법은 DI 모델에 비해 계면 기하 구조 및 운동 에너지 측면에서 훨씬 낮은 오차를 보이며 증기 함몰의 안정적 및 불안정적(keyhole) 역학을 성공적으로 예측했다.
• 3D 견고성: 노출된 금속판의 정적 레이저 유도 용융에 대한 3D 개념 증명 시뮬레이션은, HSDH 접근법이 수치적 불안정성 없이 keyhole 형성, 파쇄 및 합체를 포함한 복잡하고 고도로 동적인 2상 유동을 처리하는 데 있어 견고함을 입증했다.
• SI-1P의 효율성: 단상 열 변형 모델(SI-1P)이 완전 결합된 SI-2P 변형과 대등한 정확도를 유지하면서도 계산 효율성을 제공함을 보여줌으로써, 이 맥락에서 가스 상의 열 전달이 무시 가능하다는 가설을 검증하였다.

의의
본 논문은 HSDI 접근법이 용융 풀 시뮬레이션에 있어 예측 정확도와 계산 효율성 사이의 전례 없는 균형을 제공한다고 주장한다. 열 문제와 유동 문제를 각각의 최적 정식화(열은 샤프, 유동은 디퓨즈)로 분리함으로써, 이 방법은 지나치게 미세한 격자를 요구하는 순수 DI 모델과 복잡한 위상 변화를 처리하는 데 어려움을 겪는 순수 SI 모델의 한계를 극복한다. 이 방법론은 용융 풀 역학에 대한 고충실도 예측을 산업적으로 유의미한 길이 및 시간 척도에서 가능한 격자 해상도로 수행할 수 있게 한다. 저자들은 이 접근법이 레이저 빔 용접 및 레이저 분말 베드 융합과 같이 강력하게 국부화된 열원을 포함하는 광범위한 온도 구동 유동 문제에 적용 가능하며, 레이 트레이싱(ray-tracing) 및 온도 의존적 흡수율과 같은 더 상세한 물리학을 통합하기 위한 토대가 된다고 강조한다.