[논문 리뷰] Phase-field theory of brine entrapment in sea ice: Short-time frozen microstructures

Silke Thoms, Bernd Kutschan외 1명|arXiv (Cornell University)|2014. 05. 01.

Arctic and Antarctic ice dynamics참고 문헌 2인용 수 3

한 줄 요약

이 논문은 질량 보존 조건 하에서 비보존적 순서 매개변수 역학과 염소 확산을 결합하여 해빙에서 염수 포획의 초기 미세구조 형성 과정을 시뮬레이션하기 위한 단계장 모델을 개발한다. 모델은 관측된 해빙 플레이트릿 간격과 잘 일치하는 약 ~198 μm의 주요 파장(파장)을 예측하며, 이는 초과냉각과 freezing point와 같은 열역학적 매개변수들이 직접적으로 미세구조 길이 척도를 결정함을 보여준다.

ABSTRACT

We analyze the early phase of brine entrapment in sea ice, using a phase field model. This model for a first-order phase transition couples non-conserved order parameter kinetics to salt diffusion. The evolution equations are derived from a Landau-Ginzburg order parameter gradient dynamics together with salinity conservation. The numerical solution of model equations by an exponential time differencing scheme describes the time evolution of phase separation between liquid water with high salinity and the ice phase with low salinity. The numerical solution in one and two dimensions indicates the formation of one dominant wavelength which sets the length scale of short-time frozen structures. A stability analysis provides the phase diagram in terms of two Landau parameters. It is distinguished an uniform ice phase, a homogeneous liquid saline water solution and a phase where solidification structures can be formed. The Landau parameters are extracted from the supercooling and superheating as well as the freezing point temperature of water. With the help of realistic parameters the distribution of brine inclusions is calculated and found in agreement with the measured samples. The size of the ice domains separating regions of concentrated seawater depends on salinity and temperature and corresponds to the size of sea ice platelets obtained from a morphological stability analysis for the solidification of salt water.

연구 동기 및 목표

• 해빙의 초기 고체화 단계에서, 특히 짧은 시간 동안 고체화되는 조건에서 염수 포함물의 형성 메커니즘을 이해하기 위해.
• 염소 질량을 보존하는 물리적으로 타당한 단계장 모델을 개발하여 빙수 네트워크의 미세구조 변화를 포착하기 위해.
• 초과냉각과 freezing point와 같은 열역학적 매개변수들이 얼음-염수 구조의 크기와 안정성에 미치는 영향을 규명하기 위해.
• 실험적 염수 층의 미세구조와 관측된 해빙의 미세구조, 특히 플레이트릿 간격에 기반하여 모델을 검증하기 위해.
• 랜도 매개변수와 물 및 염수 시스템의 측정 가능한 물리적 성질 간의 관계를 설정하기 위해.

제안 방법

• 일阶 전이에 대한 랜도-긴즈버그 기울기 역학을 기반으로 한 단계장 모델을 수립하여 순서 매개변수의 진화와 염소 확산을 결합한다.
• 생성 기능에서 유도된 연립 진화 방정식을 사용하여 염소 질량 보존을 보장하며, 터링 모델과 같은 비보존적 모델과 구별된다.
• 일차원 및 이차원에서 시간에 따라 변화하는 편미분 방정식의 수치적 해를 구하기 위해 지수 시간 차분 기법을 사용한다.
• 선형 안정성 분석을 수행하여 두 랜도 매개변수인 freezing 매개변수(α₁)와 구조 매개변수(α₃)로 표현된 단계도를 유도한다.
• 초과냉각, 초과가열, 물의 freezing point(T₀ = 233.15 K) 등의 물리적 성질을 사용하여 모델 매개변수를 校정한다.
• 안정성 분석에서 유도된 κc를 바탕으로 임계 영역 크기 λc는 λc = 2π/κc로 정의된다.

실험 결과

연구 질문

• RQ1 해빙에서 염수 포획의 초기 단계에 어떤 미세구조 길이 척도가 나타나며, 이는 열역학적 조건에 의해 어떻게 결정되는가?
• RQ2 단계장 모델의 freezing 매개변수와 구조 매개변수가 안정적인 얼음-염수 패턴 형성에 미치는 영향는 어떠한가?
• RQ3 모델이 자연 해빙에서 관측된 실험적 염수 층의 미세구조와 플레이트릿 간격을 어느 정도 정확하게 재현할 수 있는가?
• RQ4 염소 질량 보존을 포함함으로써 비보존적 모델 대비 물리적 현실성은 어떻게 향상되는가?
• RQ5 모델의 임계 영역 크기와 실제 해빙에서 관측된 200–500 μm의 얼음 램셀 간격 간의 관계는 어떠한가?

주요 결과

• 실제 조건(α₁ = 0.111482, α₃ = 1.99, T₀ = 233.15 K, ΔT = 0.032 K 초과냉각)에서 모델은 관측된 해빙 플레이트릿 간격과 일치하는 임계 영역 크기 λc = 198 μm를 예측한다.
• 선형 안정성 분석을 통해 약 ~198 μm의 주요 파장이 유도되며, 이는 시스템에서 가장 불안정한 모드에 해당한다.
• 모델의 매개변수는 초과냉각, 초과가열, 물의 freezing point 등 물리적 성질과 직접 연결되어 있어 물리적 校정이 가능하다.
• 높은 초과냉각 조건(ΔT = 6.3 K)에서 임계 영역 크기 λc = 0.8 μm가 도출되었으며, 이는 염수의 응고에 대한 형태 안정성 이론과 일치한다.
• 안정성 분석에서 유도된 단계도는 α₁과 α₃에 의해 결정되는 영역에서 공간적 구조가 형성됨을 보여주며, 이는 확산도와 무관하다.
• 이차원 수치 시뮬레이션은 수직으로 배열된 시트 형태의 염수 포함물을 재현하며, 실험적으로 관측된 염수 층의 형태와 일치한다.