[论文解读] Phase-field theory of brine entrapment in sea ice: Short-time frozen microstructures

Silke Thoms, Bernd Kutschan+1 more|arXiv (Cornell University)|May 1, 2014

Arctic and Antarctic ice dynamics参考文献 2被引用 3

一句话总结

本文提出一种相场模型，用于模拟海冰中卤水包裹的早期微结构形成过程，将非守恒序参量动力学与质量守恒下的盐分扩散相耦合。该模型预测的主导波长约为198 μm，与观测到的海冰板状间距高度一致，表明过冷度和凝固点等热力学参数可直接决定微结构长度尺度。

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.

研究动机与目标

• 理解在短时冻结条件下，海冰初始冻结阶段卤水包裹的形成机制。
• 开发一种基于物理原理的相场模型，实现盐分质量守恒，并捕捉卤水网络的微结构演化过程。
• 确定过冷度和凝固点等热力学参数如何影响冰-卤水结构的尺寸与稳定性。
• 通过实验观测的卤水层纹理和自然海冰微结构（特别是板状间距）对模型进行验证。
• 建立朗道参数与水及盐水体系可测量物理性质之间的关联。

提出的方法

• 基于朗道-金兹堡梯度动力学构建相场模型，描述一阶相变过程，将序参量演化与盐分扩散相耦合。
• 从生成泛函推导出耦合演化方程，确保盐分质量守恒，从而与非守恒模型（如图灵模型）相区别。
• 采用指数时间差分法求解一维和二维时间依赖偏微分方程。
• 进行线性稳定性分析，以两个朗道参数（凝固参数α₁和结构参数α₃）表示相图。
• 利用物理参数校准模型：过冷度、过热度及水的凝固点（T₀ = 233.15 K）。
• 将临界畴尺寸λc与主导波长关联，关系为λc = 2π/κc，其中κc由稳定性分析导出。

实验结果

研究问题

• RQ1 在海冰卤水包裹的早期形成阶段，何种微结构长度尺度出现？其如何由热力学条件决定？
• RQ2 相场模型中的凝固参数与结构参数如何调控稳定冰-卤水图案的形成？
• RQ3 该模型在多大程度上能再现实验观测到的卤水层纹理及自然海冰中的板状间距？
• RQ4 与非守恒模型相比，引入盐分质量守恒如何提升相场模型的物理真实性？
• RQ5 该模型的临界畴尺寸与真实海冰中观测到的200–500 μm冰片层间距之间存在何种关系？

主要发现

• 在真实条件下（α₁ = 0.111482，α₃ = 1.99，T₀ = 233.15 K，ΔT = 0.032 K过冷度），模型预测的临界畴尺寸为λc = 198 μm，与观测到的海冰板状间距一致。
• 线性稳定性分析显示，主导波长约为198 μm，对应系统中最不稳定的模态。
• 模型参数与物理性质直接关联：过冷度、过热度及水的凝固点，支持物理校准。
• 在高过冷度（ΔT = 6.3 K）下，获得临界畴尺寸λc = 0.8 μm，与盐水凝固的形态稳定性理论一致。
• 由稳定性分析导出的相图识别出空间结构形成的区域，仅由α₁和α₃决定，与扩散率无关。
• 二维数值模拟再现了垂直取向、片状结构的卤水包裹，与实验观测到的卤水层形态一致。