구형 입자용 삼상 스테판 문제의 고정격자 날카로운 경계 수치 해법

초록

본 논문은 밀도 차이와 운동에너지 항을 포함한 삼상(고체‑액체‑기체) 스테판 문제를 구형 입자에 적용한다. 유도된 소시간 해를 초기값으로 사용하고, 고정격자 날카로운 경계(immersed‑boundary) 방법으로 2차 정확도의 수치 해를 얻는다. 두‑상 나노입자 용융 결과와 비교해 검증했으며, 나노 규모에서는 운동에너지 항이 용융·비등 시간을 크게 좌우함을 확인했다.

상세 분석

이 연구는 기존 2상 스테판 모델이 갖는 “밀도 동일”·“잠열만 고려”라는 제한을 넘어, 고체‑액체‑기체 3상 전이를 동시에 기술한다는 점에서 혁신적이다. 구형 입자를 대상으로 삼상 경계가 각각 고체‑액체(용융면), 액체‑기체(비등면), 외부 대기(전체 경계)로 존재하며, 각 경계에서 온도·열플럭스·밀도·운동에너지의 불연속성을 Rankine‑Hugoniot 조건으로 정확히 구현한다. 특히, 운동에너지 항 δₘ 은 ρ 비율이 큰 경우(예: 금속‑증기)에도 유의미한 영향을 미치며, 이는 기존 연구에서 무시된 부분이다.

수치 해법은 고정격자 방식으로, 인터페이스를 Dirac δ 함수 형태의 표면 장력·곡률 항과 함께 immersed‑boundary(IB) 기법으로 추적한다. 공간·시간 모두 2차 중심 차분/볼륨 스킴을 사용해 수렴성을 확보했으며, 인터페이스 이동을 명시적으로 업데이트한다. 초기 조건은 “소시간 해석”을 통해 얻는다. ε ≪ 1 스케일링을 도입해 r = 1‑ε r̂, t = ε t̂ 형태로 방정식을 전개하고, 밀도비(ρ_SL, ρ_LV 등)와 열전도도비(κ_SL) 등을 포함한 비차원 파라미터를 정의한다. 이 과정에서 고체‑액체 경계의 Gibbs‑Thomson 관계와, 비등면의 온도 T_im = T_v − σ_VL T_v/(ρ_LV L_v R₂) 등도 반영한다.

검증 단계에서는 기존 문헌(Font et al.)의 2상 구형 나노입자 용융 해를 재현함으로써, 고정격자 IB 방법이 좌표 변환 방식(고정 도메인 매핑)과 동등하거나 우수함을 입증한다. 이후 3상 시뮬레이션을 수행해, 입자 반경이 수십 나노미터 수준일 때 비등면의 운동에너지 항이 용융·비등 속도를 약 30 % 이상 가속시킴을 확인했다. 반면, 마이크론 규모에서는 δₘ → 0이 되어 기존 2상 모델과 차이가 사라진다.

이 논문의 핵심 기여는 (1) 구형 좌표계에서 3상 스테판 문제를 전반적으로 수식화하고, (2) 밀도·운동에너지 불연속을 포함한 정확한 인터페이스 조건을 도출한 것, (3) 고정격자 IB 방법으로 2차 정확도를 유지하면서 복잡한 이동 경계를 직접 해결한 것, (4) 소시간 해를 이용한 초기화 기법을 제시해 수치적 안정성을 확보한 점이다. 이러한 접근은 금속 3D 프린팅·와이어링 등 고에너지 레이저 가공 공정에서 입자‑플라즈마 상호작용을 다중상 CFD 모델에 통합하는 기반을 제공한다.