Numerical homogenization of a heat transfer problem with phase transitions in heterogeneous frozen soils

광활하고 얼어붙은 북부 지역에서 지표면은 단순한 흙이 아닙니다. 그것은 세계 최대의 도시, 파이프라인, 그리고 도로를 지탱하는 구조적 토대입니다. 수 세기 동안 이 영구 동토층은 단단하게 유지되어 왔으며, 그 결빙 상태는 인간의 건설을 위한 견고한 기반암 역할을 해왔습니다. 그러나 기후가 온난해짐에 따라 이 얼어붙은 땅이 녹기 시작했습니다. 토양 내부의 얼음이 녹으면 지표면은 강도를 잃게 되고, 이는 불균일한 침하와 건물의 균열로 이어집니다. 이 해빙이 정확히 얼마나 깊게 진행될지 예측하는 것은 안전과 공학의 문제입니다. 하지만 토양 자체는 거대한 난관을 제시합니다. 실제 토양은 균일한 얼음 덩어리가 아니라, 서로 다른 재료들이 층, 렌즈, 포켓 형태로 뒤섞인 혼돈스러운 혼합물이며, 불과 몇 센티미터의 거리에서도 그 특성이 변합니다. 이러한 복잡한 지표면의 녹는 현상을 완벽하게 정확하게 시뮬레이션하려면 컴퓨터가 모든 층을 추적해야 하는데, 이는 현재의 기술로는 대규모 장기 예측을 수행하기에 부족한 계산 능력을 요구합니다.

이를 해결하기 위해 연구자들은 수치적 균질화(numerical homogenization)라고 불리는 기법을 사용합니다. 스펀지를 통과하는 물의 흐름을 예측한다고 상상해 보십시오. 스펀지의 모든 작은 구멍과 섬유를 모델링하는 대신, 스펀지를 물이 통과할 수 있는 '평균적인' 능력을 가진 하나의 재료로 취급하는 것입니다. 이것이 균질화의 본질입니다. 즉, 무질서하고 상세한 현실을 컴퓨터에서 빠르게 실행 가능한 더 단순하고 평균화된 모델로 대체하는 것입니다. 그러나 지표면이 녹을 때, 이 평균화 작업은 함정이 됩니다. 얼어붙은 토양과 녹은 토양 사이의 경계인 '전선(front)'은 컴퓨터 계산 블록의 가장자리에서 깔끔하게 멈추지 않습니다. 그것은 블록의 한복판을 가로질러 지나갑니다. 여기서 딜레마가 발생합니다. 만약 토양 블록이 절반은 얼어 있고 절반은 녹아 있다면, 컴퓨터는 그 성질을 어떻게 결정해야 할까요? 얼음처럼 행동해야 할까요, 물처럼 행동해야 할까요, 아니면 그 중간 어딘가여야 할까요? 여기서 내리는 선택이 모델이 건물의 침하를 예측할지 아니면 건물이 잘 서 있을지를 결정합니다.

최근 연구에서 북동부 연방 대학교의 티무르 자카로프(Timur Zakov)는 이 세 가지 '혼합된' 토양 블록을 처리하는 방법을 테스트함으로써 이 특정 딜레마를 다루었습니다. 연구는 얼어붙은 지면 위에 놓인 가열된 건물을 대상으로 진행되었습니다. 이는 구조물의 온기가 아래의 토양을 천천히 녹이는 시나리오입니다. 연구팀은 자연적인 지표면의 혼돈을 모사하기 위해 실제적인 무작위 패턴의 토양 전도성을 사용하여 이 환경의 디지털 트윈을 만들었습니다. 그런 다음 그들은 각 토양 블록이 실제 흙과 얼음의 층보다 훨씬 큰 조밀한 격자(coarse grid)에서 시뮬레이션을 실행하여, 상변화 전선이 블록의 중앙을 가로지르게 만들었습니다. 그들은 이 잘린 블록 내부에서 어떤 일이 일어나는지 결정하기 위해 세 가지 뚜렷한 방법을 비교했습니다.

첫 번째 방법인 테이블 접근법(table approach)은 모든 가능한 토양 블록 상태에 대한 속성 목록을 미리 계산해 두고, 블록의 평균 온도에 따라 답을 찾아내는 방식입니다. 두 번째 방법인 노달 스킴(nodal scheme)은 블록의 모서리에서 토양의 상태를 결정하고 그 결정들을 평균화합니다. 세 번째 방법인 h-스무딩(h-smoothing)은 선형 온도 구배를 바탕으로 녹은 블록의 정확한 비율을 계산하며, 본질적으로 물로 변한 블록의 조각을 측정합니다. 연구진은 이 시뮬레이션을 2차원과 3차원에서 모두 실행하여, 모든 층을 추적하는 매우 상세한 '참조(reference)' 솔루션과 비교했습니다.

연구 결과는 엔지니어가 도구를 선택할 때 직면하게 될 놀라운 트레이드오프(trade-off)를 드러냈습니다. 목표가 토양의 온도를 예측하는 것이라면, h-스무딩 방법이 가장 정확한 것으로 나타났습니다. 이 방법은 다른 방법들이 어려움을 겪었던 조밀하지 않은 격자에서도 열 흐름을 가장 정밀하게 추적했습니다. 그러나 연구진이 녹는 전선, 즉 얼어붙은 지면과 녹은 지면 사이의 실제 경계를 조사했을 때, 이야기는 완전히 달라졌습니다. h-스무딩 방법은 온도 측면에서는 우수했지만, 해빙 전선의 위치를 예측할 때는 지속적으로 실제보다 더 깊게 녹는 것으로 과대평가했습니다. 즉, 지표면이 실제보다 더 깊게 녹을 것이라고 예측했습니다.

이러한 불일치의 원인은 에너지 보존 법칙에 있습니다. h-스무딩 방법은 녹는 구역의 형태를 정밀하게 표현하려는 노력 속에서, 매 시간 단계마다 시스템에 미세한 추가 열을 의도치 않게 더하게 되었습니다. 1년 동안 이 미세한 잉여 에너지는 축적되었고, 결과적으로 시뮬레이션된 지표면이 실제 지표면보다 더 많이 녹게 만들었습니다. 반면, 다른 두 방법, 특히 테이블 접근법과 노달 스킴은 '보수적'이었습니다. 이들은 에너지를 추가하거나 손실시키지 않았으며, 그 결과 온도 측정값은 약간 덜 정밀할지라도 해빙의 깊이는 훨씬 더 정확하게 예측했습니다. 이 연구는 만약 연구자가 토양의 정확한 온도를 알아야 한다면 h-스무딩이 더 우수한 선택이라는 것을 보여주었습니다. 하지만 목표가 해빙으로 인해 기초가 얼마나 깊게 가라앉을지를 아는 것이라면, 보수적인 방법들이 훨씬 더 신뢰할 수 있다는 것을 보여주었습니다.

또한 연구는 토양이 녹는 온도가 곳곳에서 동일하지 않고 서로 다른 층에 따라 변하는 경우를 탐구했습니다. 이 더 복und한 시나리오에서, 상세한 속성 조회를 사용하는 테이블 기반 방법은 두 텐서(two-tensor) 방법보다 해빙 전선의 위치를 예측하는 데 있어 더 뛰어난 성능을 보였습니다. 이는 매우 가변적인 지표면 조건에서는 더 상세한 사전 계산된 토양 상태 지도가 필요함을 시사합니다. 궁극적으로 이 연구는 엔지니어와 과학자들에게 명확한 가이드를 제공합니다. 즉, 단 하나의 '최선'의 방법은 없다는 것입니다. 선택은 전적으로 어떤 질문을 던지느냐에 달려 있습니다. 열적 정확도가 우선이라면 한 방법이 승리하고, 해빙의 물리적 경계가 우선이라면 다른 방법이 앞서 나갑니다. 이러한 차이점을 명확히 함으로써, 이 연구는 북극의 건물들을 보호하기 위해 사용되는 모델들이 단순히 빠른 것을 넘어 신뢰할 수 있도록 돕습니다.