Laser Powder Bed Fusion Melt Pool Dynamics for Different Geometric Variations and Powder Layer Heights: High-Fidelity Multiphysics Modeling vs 2025 NIST Experiments

이 논문은 오픈 소스 솔버 'LaserBeamFoam'을 활용한 고충실도 다물리 시뮬레이션이 2025 년 NIST 실험 결과와 정량적으로 일치하여, 금속 적층 제조의 용융 풀 역학 및 공정 최적화를 위한 예측 도구로서의 유효성을 입증했습니다.

핵심

이 논문은 **3D 프린팅 (특히 금속 3D 프린팅)**의 핵심인 '용융 풀 (녹아내리는 금속의 물방울)'이 어떻게 움직이는지, 그리고 그 모양이 어떻게 결정되는지를 컴퓨터로 아주 정밀하게 시뮬레이션한 연구입니다.

쉽게 비유하자면, **"고온의 레이저로 모래 (금속 가루) 위에 그림을 그릴 때, 그 그림이 어떻게 변형되는지 예측하는 '디지털 예언자'를 개발한 이야기"**라고 할 수 있습니다.

주요 내용을 일상적인 비유로 설명해 드릴게요.

---

1. 연구의 배경: 왜 이 연구가 필요할까요?

금속 3D 프린팅은 항공기 부품이나 인공 관절처럼 정밀하고 튼튼한 물건을 만듭니다. 하지만 레이저로 금속 가루를 녹일 때, 가루의 두께나 부품의 모양에 따라 녹는 금속의 모양 (용융 풀) 이 달라집니다.

• 비유: 마치 스프레이로 페인트를 칠할 때를 생각해 보세요. 벽이 평평하면 페인트가 고르게 퍼지지만, 벽에 모래알이 얹혀 있거나 벽의 크기가 작으면 페인트가 튀거나 모양이 달라집니다.
• 문제점: 지금까지는 레이저의 힘이나 속도만 조절하면 된다고 생각했지만, **'가루의 두께'**나 **'부품의 크기'**가 결과물에 얼마나 큰 영향을 미치는지는 잘 몰랐습니다.

2. 연구 방법: '가상의 실험실' 만들기

연구팀은 실제 실험을 하기 전에 컴퓨터 안에서 정밀한 시뮬레이션을 돌렸습니다. 이를 위해 OpenFOAM이라는 오픈소스 소프트웨어를 기반으로 한 **'LaserBeamFoam'**이라는 도구를 사용했습니다.

• 비유: 이 도구는 마치 가상 현실 (VR) 게임과 같습니다.
  • 실제 실험실에서는 레이저를 쏘고 가루를 뿌려서 실패하면 비용이 많이 들지만, 이 시뮬레이션은 컴퓨터 안에서 수천 번 실패해도 비용이 들지 않습니다.
  • 이 프로그램은 열, 액체 금속의 흐름, 증발, 레이저 반사 등 모든 물리 법칙을 동시에 계산합니다. 마치 날씨 예보 모델이 바람, 기온, 습도를 모두 고려하듯, 이 모델도 금속이 녹는 모든 과정을 고려합니다.

3. 핵심 발견: 가루 두께에 따른 '흡수율'의 비밀

가장 중요한 발견은 가루의 두께가 레이저를 얼마나 잘 흡수하는지에 영향을 준다는 점입니다.

• 비유:
  • 바닥이 평평한 경우 (가루 없음): 레이저가 바닥에 닿으면 반사되어 많이 튕겨 나갑니다. (햇빛이 거울에 반사되는 것처럼요.)
  • 가루가 얇게 깔린 경우: 가루 입자들이 레이저를 여러 번 반사시키며 에너지를 가둡니다. (숲속에서 햇빛이 나뭇잎 사이로 스며들며 여러 번 반사되는 것처럼요.)
  • 가루가 두껍게 깔린 경우: 더 이상 반사가 늘어나지 않고 일정한 수준에 도달합니다. (숲이 너무 깊어져서 더 이상 빛이 들어오지 않는 것처럼요.)

연구팀은 이 현상을 수학적으로 설명하는 **'유효 흡수율 공식'**을 개발했습니다. 가루가 두꺼워질수록 레이저가 더 잘 흡수된다는 사실을 공식에 넣어서, 시뮬레이션의 정확도를 극적으로 높였습니다.

4. 검증: NIST 의 '시험 문제' 풀기

이 연구의 가장 큰 특징은 **미국 국립표준기술연구소 (NIST)**에서 2025 년에 진행한 'AM-Bench'라는 정밀한 시험 문제를 풀었다는 점입니다.

• 비유: 마치 수학 경시대회에서 출제된 매우 어려운 문제를 풀고, 정답지 (실제 실험 데이터) 와 비교해 본 것과 같습니다.
  • 연구팀은 가루 두께가 0μm(평평한 바닥), 80μm, 160μm 인 세 가지 상황과 부품 크기 (1cm x 5cm, 5cm x 5cm) 를 모두 시뮬레이션했습니다.
  • 결과: 컴퓨터가 예측한 금속의 깊이, 너비, 높이 등이 실제 실험 결과와 놀라울 정도로 일치했습니다. 오차가 매우 적었습니다.

5. 의의: '디지털 트윈'의 완성

이 연구는 단순히 실험을 대체하는 것을 넘어, 미래의 3D 프린팅을 설계하는 핵심 도구가 되었습니다.

• 비유: 이제 우리는 **디지털 트윈 (Digital Twin)**이라는 개념을 사용할 수 있게 되었습니다.
  • 실제 비행기를 만들기 전에, 컴퓨터 안에 완벽한 '가상의 비행기'를 만들어서 레이저를 쏘고 녹는 과정을 미리 확인합니다.
  • 만약 가루 두께를 80μm 로 하면 부품이 약해질 수 있다는 것을 컴퓨터에서 미리 알아내면, 실제 실험을 하기 전에 설계도를 수정할 수 있습니다.

요약

이 논문은 **"금속 3D 프린팅에서 가루의 두께와 부품 크기가 녹는 금속 모양에 미치는 영향을, 컴퓨터 시뮬레이션으로 아주 정확하게 예측하는 방법을 찾아냈다"**는 내용입니다.

이는 마치 요리사가 "재료를 얼마나 두껍게 깔아야 맛이 잘 우러나는지"를 실험 없이도 정확히 계산해내는 것과 같습니다. 덕분에 앞으로 더 정밀하고 결함이 없는 금속 부품을 저렴하고 빠르게 만들 수 있는 길이 열렸습니다.

기술 요약

논문 요약: 레이저 분말 베드 융합 (LPBF) 의 용융 풀 역학 및 NIST 2025 실험 검증

1. 연구 배경 및 문제 제기 (Problem)

• 배경: 금속 레이저 분말 베드 융합 (PBF-LB/M) 은 항공우주, 의료, 에너지 분야에서 복잡한 금속 부품을 제조하는 핵심 적층 제조 기술입니다. 그러나 기공 (porosity), 균열 (cracking), 잔류 응력 등 결함 발생이 빈번하여 부품의 신뢰성과 품질 보증이 주요 과제입니다.
• 문제점: 기존 연구들은 주로 레이저 출력, 스캔 속도, 스포트 직경 등의 공정 변수가 용융 풀 (melt pool) 역학에 미치는 영향을 다루었습니다. 그러나 **분말 층 두께 (powder layer height)**와 **부품 기하학적 형상 (geometric variations)**이 용융 풀 거동에 미치는 영향은 상대적으로 덜 이해되어 있습니다.
• 한계: 기존 고충실도 (High-fidelity) 시뮬레이션은 계산 비용 문제로 인해 단일 트랙 (single-track) 또는 소수의 트랙 (최대 8 개) 에 국한되는 경우가 많았습니다. 또한, 분말 베드 특성과 멀티 트랙 스캐닝을 모두 고려한 통합 모델링과 이를 검증할 수 있는 잘 제어된 벤치마크 데이터의 부재가 한계였습니다.

2. 방법론 (Methodology)

• 시뮬레이션 프레임워크: 오픈소스 유한체적법 (FVM) 솔버인 **'LaserBeamFoam' (OpenFOAM 기반)**을 사용하여 고충실도 다물리 (multiphysics) 시뮬레이션을 수행했습니다.
• 물리 모델링:
  • 연속체 모델: 열전달, 유체 흐름, 증발, 반동 압력 (recoil pressure), 마랑고니 대류, 자유 표면 진화를 통합하여 모델링했습니다.
  • 레이저 상호작용: 가우시안 열원과 프레넬 (Fresnel) 기반 흡수율, 그리고 키홀 (keyhole) 내부의 다중 반사를 고려한 레이 트레이싱 (Ray-tracing) 알고리즘을 적용하여 레이저 에너지 흡수를 정밀하게 모사했습니다.
  • 상변화 및 분말 처리: 엔탈피 - 다공성 (enthalpy-porosity) 기법을 사용하여 용융/응고 과정을 처리했습니다. 특히, 분말 입자가 완전히 용융되기 전까지 유동성을 갖지 않도록 하는 '마스크된 액체 분율 (masked liquid fraction)' 개념을 도입하여 분말 베드의 물리적 거동을 정확히 표현했습니다.
• 검증 데이터: NIST (미국 국립표준기술연구소) 의 AM-Bench 2025 챌린지 (AMB2025-06-PMPG) 데이터를 벤치마크로 사용했습니다.
  • 조건: IN718 합금을 사용하며, 0µm(베어 플레이트), 80µm, 160µm 의 분말 층 두께와 1mm×5mm, 5mm×5mm 의 두 가지 기하학적 형상을 비교했습니다.
  • 공정 파라미터: 레이저 출력 285W, 스캔 속도 960mm/s, 해치 간격 0.11mm, 스포트 직경 72µm 로 고정.
• 모델 개선 (Enhanced Model): 초기 'NIST 제출 모델'은 분말 베드에서의 예측 정확도가 낮았습니다. 이를 개선하기 위해 **분말 층 두께에 따른 유효 레이저 흡수율 (effective laser absorptivity, $\eta_{eff}$)**을 지수 함수 형태로 모델링하여 도입했습니다. 이는 분말 베드 내부의 다중 반사로 인한 에너지 결합 증가를 계수화한 것입니다.

3. 주요 기여 (Key Contributions)

1. 고충실도 멀티 트랙 시뮬레이션: 계산 비용의 한계에도 불구하고, 45 개의 병렬 레이저 스캔 트랙을 시뮬레이션하여 NIST 벤치마크와 직접 비교할 수 있는 데이터를 생성했습니다.
2. 층 두께 의존적 흡수율 모델 개발: 분말 층 두께가 변함에 따라 레이저 흡수율이 어떻게 변화하는지 (베어 플레이트에서 분말 베드로 전환 시의 점근적 행동) 를 수학적으로 모델링하여, 다양한 층 두께 조건에서 일관된 예측 정확도를 달성했습니다.
3. NIST 2025 데이터와의 정량적 검증: 베어 플레이트뿐만 아니라 80µm 및 160µm 분말 층 조건에서도 용융 풀 깊이, 너비, 비드 높이, 중첩 깊이, 응고 면적, 희석 면적 등 모든 기하학적 메트릭에서 실험 데이터와 높은 일치도를 보였습니다.
4. 확장성 있는 예측 기법: 15 트랙까지의 고해상도 시뮬레이션 결과를 기반으로 지수 함수 및 선형 피팅을 통해 45 트랙까지의 거동을 외삽 (extrapolation) 하는 방법을 제시했습니다.

4. 결과 (Results)

• 정량적 일치: 개선된 모델 (Enhanced Model) 은 NIST 실험 데이터와 비교하여 용융 풀 깊이, 너비, 비드 높이 등 주요 지표에서 20% 이내의 오차를 보였으며, 특히 베어 플레이트 및 80µm/160µm 분말 층 조건에서 매우 우수한 예측 능력을 입증했습니다.
• 기하학적 영향:
  • 5mm×5mm 패드: 스캔 위치에 따라 잔류 열의 영향이 다르게 나타나며, 시뮬레이션이 이러한 열적 이력 효과를 정확히 포착했습니다.
  • 1mm×5mm 패드: 좁은 폭으로 인해 이전 트랙의 잔류 열이 다음 트랙의 용융 역학에 큰 영향을 미쳐 재용융 (remelting) 현상이 두드러졌으며, 모델이 이를 잘 재현했습니다.
• 흡수율 분석: 분말 층 두께가 증가함에 따라 레이저 흡수율이 비선형적으로 증가하여 포화되는 경향을 보였으며, 제안된 지수 포화 함수 ($\eta_{eff}(L)$) 가 이를 잘 설명했습니다.
• 오차 원인: 비드 높이와 응고 면적에서 실험 대비 약간 과대 예측된 것은 잔류 열 효과의 정밀한 해석 부족과 45 트랙 전체를 직접 시뮬레이션하지 않고 외삽을 사용했기 때문으로 분석되었습니다.

5. 의의 및 결론 (Significance)

• 공정 최적화 및 결함 완화: 분말 층 두께와 부품 형상 변화에 따른 용융 풀 거동을 정확히 예측할 수 있게 되어, LPBF 공정의 최적화와 결함 (lack-of-fusion, balling 등) 방지에 기여합니다.
• 디지털 트윈 통합: 본 연구에서 개발된 검증된 다물리 시뮬레이션 프레임워크는 적층 제조를 위한 디지털 트윈 (Digital Twin) 시스템 구축의 핵심 기초 기술로 활용될 수 있습니다.
• 학문적 혁신: 다양한 분말 층 두께와 부품 치수에 대해 NIST 2025 실험 데이터와 정량적으로 일치하는 첫 번째 계산적 접근법으로, 고충실도 시뮬레이션이 실험적 검증과 어떻게 결합될 수 있는지 보여주는 중요한 사례입니다.

이 연구는 계산 역학과 실험적 검증의 긴밀한 협력을 통해 적층 제조 공정의 예측 능력을 획기적으로 향상시켰으며, 향후 산업적 적용을 위한 신뢰할 수 있는 시뮬레이션 도구로서의 가능성을 입증했습니다.