[論文レビュー] 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の主な波長を予測し、過冷却度や 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 といった熱力学的パラメータが、氷-塩水構造のサイズと安定性にどのように影響するかを特定すること。
• 実験的塩水層のテクスチャーや観察された海氷の微細組織、特にプレートレット間隔と比較してモデルを検証すること。
• ランダウパラメータと水および塩水系の測定可能な物理的性質との間の関係を確立すること。

提案手法

• 一次相転移を記述するランダウ＝ジンブルグ勾配ダイナミクスに基づくフェーズフィールドモデルを構築し、秩序パラメータの進化と塩の拡散を結合した。
• 生成汎関数から導出した連立進化方程式を用いて塩の質量保存を保証し、チューリングモデルのような非保存的モデルとは区別した。
• 1次元および2次元の時間依存偏微分方程式を数値的に解くために指数型時間差分法を用いた。
• 線形安定性解析により、2つのランダウパラメータ（凍結パラメータ α₁ と構造パラメータ α₃）を用いた相図を導出した。
• 過冷却度、過熱度、および水の freezing point（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として出現し、これは系における最も不安定なモードに対応した。
• モデルのパラメータは過冷却度、過熱度、および水の freezing point といった物理的性質と直接関連しており、物理的キャリブレーションが可能である。
• 高い過冷却度（ΔT = 6.3 K）下では臨界ドメインサイズ λc = 0.8 μm が得られ、塩水の凝固における形態的安定性理論と整合的であった。
• 安定性解析から導かれた相図により、空間的構造が形成される領域が、拡散度とは無関係に α₁ と α₃ のみに依存して特定された。
• 2次元の数値シミュレーションにより、垂直に配向したシート状の塩水インクルージョンが再現され、実験的に観察された塩水層の形態と一致した。