다공성 매체 열전달 방정식에 대한 적합도 기반 방사형 기저함수 적분·미분 콜로케이션 비교

초록

본 논문은 다공성 매체 내 자연 대류 열전달 방정식에 두 가지 방사형 기저함수(RBF) 콜로케이션 방법을 적용한다. 하나는 적분 기반 IRBF, 다른 하나는 미분 기반 DRBF이다. 수치 실험 결과 IRBF가 DRBF에 비해 수렴 속도가 빠르고 안정성이 높으며, 고정밀 해를 얻는 데 유리함을 확인하였다.

상세 분석

본 연구는 방사형 기저함수(RBF)를 이용한 콜로케이션 기법을 두 가지 형태로 구현한다. 첫 번째는 미분 연산을 직접 적용하는 전통적인 DRBF(Differentiation‑RBF)이며, 두 번째는 미분 대신 적분 연산을 활용하는 IRBF(Integration‑RBF)이다. DRBF는 기저함수의 미분값을 계산해 미분 방정식에 바로 대입하는 방식으로, 기저함수의 고차 미분이 필요할 경우 수치적 불안정성과 큰 조건수를 초래한다. 반면 IRBF는 미분 연산을 적분 형태로 변환한 뒤, 적분된 RBF를 이용해 시스템을 구성한다. 이 과정에서 미분 연산에 비해 노이즈에 대한 민감도가 낮아지고, 행렬의 스펙트럼이 더 균일해져 수치 해석의 안정성이 크게 향상된다.

연구 대상은 다공성 매체 내 자연 대류에 의해 발생하는 비선형 열전달 방정식이다. 이러한 방정식은 비등방성 열전도, 다공성 매체의 투과성, 그리고 방사형 대류 계수 등을 포함해 복잡한 비선형성을 띤다. 특히 원자력 폐기물 저장용 캔스터 설계에서는 온도 분포와 열전달 효율이 안전성 평가의 핵심이므로, 고정밀 수치 해법이 필수적이다. 논문은 차원 축소와 비선형 변환을 통해 2차원 축축한 방정식으로 정리한 뒤, 경계 조건(디리클레·노이만 혼합)과 초기 조건을 RBF 콜로케이션 노드에 적용한다.

실험에서는 다중 스케일의 RBF(가우시안, 멀티쿼드리틱)와 다양한 형태의 노드 배치를 검토하였다. 결과는 IRBF가 동일한 노드 수에서 DRBF보다 평균 절대 오차가 1~2 차수 감소하고, 수렴률이 2차 이상 향상됨을 보여준다. 또한, IRBF는 행렬의 조건수가 DRBF에 비해 평균 10배 이상 낮아, 선형 시스템 해석 시 발생하는 수치 오차와 연산 비용을 크게 절감한다. 이러한 장점은 특히 고해상도 시뮬레이션이 요구되는 원자력 폐기물 저장 설계 단계에서 실용적이다.

요약하면, IRBF는 DRBF가 갖는 미분 연산의 수치적 취약점을 적분 변환을 통해 극복하고, 다공성 매체 열전달 문제에 대해 높은 정확도와 빠른 수렴을 제공한다. 이는 향후 복합 물리 현상 모델링에서 RBF 기반 콜로케이션 방법을 선택할 때 중요한 지침이 될 것이다.