Non-invasive imaging of solute redistribution below evaporating surfaces using 23Na-MRI

본 연구는 $^{23}\text{Na-MRI}$를 활용하여 고투과성 다공성 매질에서의 밀도 구동 핑거링 불안정성이 증발하는 용질의 상당한 하향 재분배를 유발함으로써 표면 포화를 방지하는 반면, 저투과성 매질은 극심한 표면 근처 농도 축적을 허용한다는 것을 실험적으로 입증한다.

핵심

개요: 소금, 모래, 그리고 "폭우" 효과

두 개의 모래 컵이 있다고 상상해 보세요. 한 컵은 거칠고 굵은 자갈(이것을 "빠른 모래"라고 부릅시다)로 채워져 있고, 다른 한 컵은 매우 미세하고 고운 가루(이것을 "느린 모래"라고 부릅시다)로 채워져 있습니다. 두 컵에 소금물을 부어 완전히 적신 다음, 물이 증발하도록 따뜻하고 건조한 바람 속에 둡니다.

연구진이 던진 질문은 이것입니다: 물이 사라질 때 소금은 어디로 가는가?

상식적으로는 물웅덩이에 막이 생기는 것처럼 소금이 맨 윗표면에만 쌓일 것이라고 생각할 수 있습니다. 하지만 이 연구에서는 특수한 "마법 카메라"(MRI)를 사용하여 모래 내부에서 실제로 어떤 일이 일 벌어지는지 관찰했고, 놀라운 사실을 발견했습니다: 그것은 전적으로 모래 구멍의 크기에 달려 있었습니다.

도구: "마법 카메라"

젖은 모래를 파헤치지 않고 내부를 보기 위해, 과학자들은 병원의 MRI 스캐너와 유사한 장치를 사용했습니다. 하지만 환자의 무릎 사진을 찍는 대신, 그들은 나트륨(소금의 'Na')을 볼 수 있도록 조율했습니다.

이것은 마치 젖은 모래 속에서 소금이 빛나는 모습을 볼 수 있는 카메라와 같습니다. 그들은 3D 스냅샷을 시간별로 찍어 소금이 어떻게 움직이는지, 마치 방 안을 이동하는 군중을 타임랩스 영상으로 보는 것처럼 관찰할 수 있었습니다.

실험: 두 가지 서로 다른 이야기

연구진은 투과성(물이 얼마나 쉽게 흐를 수 있는지)이 매우 다른 두 종류의 모래를 대상으로 실험을 진행했습니다.

1. "느린 모래" (미세한 입자)

• 무슨 일이 일어났나: 표면에서 물이 증발함에 따라, 소금은 위로 올라갈 곳이 없었습니다. 소금은 표면 근처의 아주 작고 좁은 틈새에 갇혀버렸습니다.
• 결과: 소금은 맨 윗부분에 두꺼운 농축 층을 형성하며 쌓였고, 거의 고체 암석(결정화)이 될 정도의 상태에 도달했습니다.
• 비유: 사람들이(물 분자) 나가려고 노력하지만 문이 너무 작은 붐비는 복도를 상상해 보세요. 사람들은 출구 바로 앞에서 갇히게 되고, "소금"(무거운 배낭)이 그곳에 쌓여 다른 사람들이 빠져나가는 것을 매우 어렵게 만듭니다. 소금이 윗부분을 막아버렸기 때문에 증발 속도가 현저히 느려졌습니다.

2. "빠른 모래" (굵은 입자)

• 무슨 일이 일어났나: 처음에는 느린 모래처럼 소금이 윗부분에 쌓이려고 했습니다. 하지만 이 모래의 구멍은 더 컸기 때문에 다른 일이 일어났습니다. 윗부분의 소금물이 매우 무거워졌습니다(밀도가 높아짐).
• 결과: 중력이 지배권을 잡았습니다. 무거운 소금물은 윗부분에 머물 수 없었기에, 마치 수영장에 떨어지는 무거운 돌처럼 모래 아래로 가라앉았습니다. 이는 소금을 기둥 깊숙이 운반하며 아래로 내려가는 "플룸(plume)" 또는 "손가락(finger)" 모양을 만들어냈습니다.
• 비유: 넓은 경기장에 있는 군중을 상상해 보세요. 사람들이 떠날 때, 무거운 배낭(소금)을 멘 한 그룹이 너무 무거워져서 앞줄에 머물 수 없게 됩니다. 대신, 그들은 군중 사이를 지나 뒤쪽으로 가라앉습니다. 윗부분은 비교적 깨끗하게 유지되며, "소금"은 모래 깊숙이 재분배됩니다.

이것이 중요한 이유 (논문에 따르면)

이 연구는 과학자들이 이전까지 추측만 해왔던 이론을 확인해 줍니다: 소금은 항상 위에 머무는 것이 아닙니다.

• 조밀한 모래에서는: 소금이 표면에 머물며 매우 농축되고, 이는 증발을 늦춥니다. 이는 건축 자재(풍화 작서를 유발함)나 토양(염류화를 유발함)과 같은 분야에 좋지 않은 영향을 미칩니다.
• 느슨한 모래에서는: 소금이 아래로 가라앉습니다. 이는 표면의 염도를 낮게 유지하고, 증발이 꾸준한 속도로 계속되게 하며, 소금이 표면에 딱딱한 껍질을 만드는 대신 땅속 깊은 곳으로 이동하게 함을 의미합니다.

"손가락" vs "껍질"

연구진은 실제 관찰 결과와 컴퓨터 시뮬레이션을 비교했습니다. 컴퓨터 모델은 MRI에서 본 것과 정확히 일치하는 결과를 예측했습니다:

• 느린 모래는 윗부분에 소금 "껍질(crust)"을 형성했습니다.
• 빠른 모래는 아래로 뚝뚝 떨어지는 소금 "손가락(fingers)"을 만들었습니다.

또한 그들은 "레이리 수(Rayleigh number)"라는 개념(유체가 충분히 무거워 가라앉을 수 있는지를 측정하는 세련된 방식)을 사용하여 수학적 검증을 했습니다. 수학적 계산은 빠른 모래가 불안정하여 가라앉을 것이라고 예측했고, 느린 모래는 그대로 머물 것이라고 예측했습니다. MRI 카메라는 이 수학적 예측이 옳다는 것을 증명했습니다.

요약

이 논문은 과학자들이 특별한 카메라를 사용하여 소금과 모래에 관한 미스터리를 해결한 탐정 이야기와 같습니다. 그들은 모래 알갱이의 크기가 교통 통제관 역할을 한다는 것을 발견했습니다:

• 작은 알갱이는 소금을 문 앞에 가두어 무거운 교통 체증(껍질)을 만듭니다.
• 큰 알갱이는 무거운 소금이 지하로 내려가도록 허용하여, 더 많은 물이 증발할 수 있도록 길을 터줍니다.

이는 호수가 마르거나, 건물을 손상시키거나, 정원의 토양에 영향을 미치는 등 자연계에서 소금이 어떻게 이동하는지 이해하는 데 도움을 줍니다.

기술 요약

기술 요약: $^{23}\text{Na-MRI}$를 이용한 증발 표면 하부의 용질 재분배에 대한 비침습적 이미징

문제 제기
다공성 매질에서의 염수 증발은 토양 염류화 및 건축 자재 풍화와 관련된 중요한 과정이다. 기존 연구들은 증발 단계(SS1–SS3)를 확인하였으나, 용질이 표면에 축적되는 역학은 여전히 복잡하다. 최근의 이론적 모델(예: Bringedal et al., 2022)은 완전히 포화된 고투과성 다공성 매질에서 증발이 밀도 구동 불안정성(density-driven instabilities)을 유발할 수 있음을 시사한다. 염분이 표면에 축적됨에 따라 발생하는 유체 밀도의 증가는 하향 핑거링 흐름(downward fingering flow)을 유발하여, 증발 전면으로부터 용질을 하부로 재분배할 수 있다. 그러나 이러한 현상은 주로 이론적으로만 관찰되었거나 2차원 Hele-Shaw 셀에서만 관찰되었다. 자연적인 윅킹(wicking) 조건의 다공성 매질 내에서 밀도 구동 하향 흐름이 발생하는지, 특히 고유 투과율이 표면 농축에서 하향 재분배로의 전이에 어떻게 영향을 미치는지에 대한 3차원 실험적 증거가 부족한 실정이다.

연구 방법론
이러한 공백을 메우기 위해, 저자들은 대조적인 고유 투과율을 가진 두 가지 다공성 매질 샘플, F36(중사, $k \approx 7.83 \times 10^{-12} \, \text{m}^2$)과 W3(매우 고운 모래, $k \approx 3.2 \times 10^{-14} \, \text{m}^2$)에 대해 증발 실험을 수행하였다. 두 샘플 모두 초기에는 1.0 mol L⁻¹ NaCl 용액으로 포화되었으며, 증발을 유도하는 난류 N₂ 가스 스트림이 존재하는 동안 포화 상태를 유지하기 위해 하부 저장소로부터 동일한 용액을 공급받았다(윅킹 조건).

본 연구의 핵심 혁신은 NaCl의 시공간적 재분배를 모니터링하기 위해 비침습적인 $^{23}\text{Na-MRI}$(나트륨-23 자기공명영상)를 적용한 것이다.

• 이미징 프로토콜: 4.7 T Bruker 스캐너를 사용하였다. W3의 더 빠른 $T_2$ 이완 시간을 수용하기 위해 서로 다른 펄스 시퀀스(F36에는 3D RARE, W3에는 MSSE)를 채택하였다.
• 교정: 알려진 NaCl 농도(0 ~ 5 mol L⁻¹)를 가진 샘플로부터 도출된 교정 곡선을 사용하여 정량적 Na 농도 지도를 생성하였다. F36에 대해서는 선형 관계가 확립된 반면, 표면 이완 효과로 인해 W3에 대해서는 2차 다항식이 필요하였다.
• 검증: 샘플이 실험 전반에 걸쳐 완전히 포화된 상태를 유지하는지 확인하기 위해 $^{1}\text{H-MRI}$를 병행 사용하였다. MRI로 도출된 Na 질량을 증발에 의한 질량 손실 데이터와 비교하여 질량 수지를 계산하였다.
• 수치 모델링: 실험 결과는 밀도 구동 흐름을 포함하고 불안정성을 유발하기 위한 작은 무작위 섭동으로 초기화된 DuMux 시뮬레이터를 사용한 2D 수치 시뮬레이션과 비교되었다.

주요 결과

1. 용질 거동의 차이: 실험 결과, 투과율에 따라 근본적으로 다른 용질 재분배 패턴이 나타났다.

   • W3 (저투과성): 표면 근처에서 점진적이고 얕은 Na 농축을 보였다. 상부 1 mm 내의 농도는 약 6.7 mol L⁻¹에 달했다(용해도 한계인 6.14 mol L⁻¹에 근접하거나 초과함). 하향 핑거링은 관찰되지 않았다. 결과적으로 표면 농도가 증가함에 따라 증기압이 감소하여 증발률이 약 9 mm d⁻¹에서 ~4 mm d⁻¹로 급격히 떨어졌다.
   • F36 (고투과성): 초기 표면 농축을 보였으나, 1시간 이내에 밀도 구동 하향 플룸(핑거링)이 발달하였다. 이 플룸은 NaCl을 기둥 깊숙이 재분배하여 극단적인 표면 축적을 방지하였다. 표면 농도는 비교적 완만하게 유지되었으며(최대 2.5 mol L⁻¹), 실험 전반에 걸쳐 증발률은 상대적으로 일정하게 유지되었다(6 mm d⁻¹).

2. 안정성 분석: 관찰된 거동은 Wooding 등(1997a)이 정의한 안정성 기준과 일치한다. F36의 레이리 수($R_\delta$)는 33.5(불안정성 지시)로 계산된 반면, W3는 0.59(안정)였다. Bringedal 등(2022)의 선형 안정성 분석은 F36의 불안정성 발생 시점을 0.61시간으로 예측하였으며, 이는 실험에서 핑거링이 약 1시간 만에 나타난 관찰 결과와 일치하였다.

3. 모델 비교: DuMux를 이용한 수치 시뮬레이션은 실험 데이터와 질적인 일치를 보였다. 모델은 W3의 표면 농축과 F36의 하향 핑거 형성을 성공적으로 재현하였다. 그러나 후기 단계의 농도 프로파일에서는 경계 효과(유한한 샘플 높이) 및 모델의 단순화된 경계 조건으로 인한 불일치가 존재하였다.

의의 및 주장
본 연구는 증발하는 포화 다공성 매질 내에서 밀도 구동 용질 재분배가 일어남을 3차원 규모에서 최초로 실험적으로 확인하였다. 저자들은 본 연구가 투과율이 표면 크러스트 형성으로 이어질지 아니면 심층 용질 재분배로 이어질지를 결정하는 결정적 요소라는 이론적 예측을 검증한다고 주장한다.

구체적으로, 논문은 다음과 같이 단언한다:

• 고투과성 매질은 증발 플럭스에 대응하는 밀도 구동 하향 흐름을 허용하여, 표면 포화 및 크러스트 형성을 지연시킨다.
• 저투과성 매질은 용질을 표면 근처에 가두어, 농도의 급격한 증가, 증발률 감소, 그리고 조기 크러스트 형성을 초래한다.
• $^{23}\text{Na-MRI}$는 신호 대 잡음비(SNR) 문제에도 불구하고, 밀도 구동 불안정성을 해상할 수 있는, 다공성 매질 내 용질 재분배를 정량화할 수 있는 유효한 비침습적 도구이다.

저자들은 이러한 발견이 플레이아(playas) 및 소금 평원과 같은 염수 지하수 시스템의 증발률, 염 결정(salt crust) 형성 시기, 그리고 수문학을 예측하는 데 시사점을 준다고 결론짓는다. 그들은 향후 연구에서 다양한 종류의 염과 불균질한 매질을 탐구해야 한다고 완곡하게 제안하지만, 주요 기여는 이전에 이론으로만 존재했던 밀도 구동 흐름 메커니즘을 실험적으로 검증한 데 있다.