[論文レビュー] Transition from classical to ultimate melting

Edoardo Bellincioni, Kevin Zhong+5 more|arXiv (Cornell University)|Jan 10, 2026

Lattice Boltzmann Simulation Studies被引用数 0

ひとこと要約

論文は乱流中の融解転移を示す：小スケール（古典的）融解は Nu ~ Re_D0^{1/2}、大スケール（究極的）融解は Nu ~ Re_D0^{0.8–1}、物体サイズの4桁分のオーダーにわたる実験とDNSで観測。

ABSTRACT

Melting is omnipresent in nature and technology, with applications ranging from metallurgy, biology, food science, and latent thermal energy storage to oceanography, geophysics, and climate science, and occurring on all scales from sub-millimeter to global scales. The key objective is to understand the rate at which an object melts as a function of its size and of the ambient conditions. To achieve this it is important to be able to extrapolate from small scale experiments and observations to large or even global scales. This is done by scaling laws. However, these are only meaningful if there is no transition from one scaling relation to another one. Here we show, however, that for both fixed and freely-advected melting objects immersed in a turbulent flow a melting transition does exist, namely from slow melting at the small scales to fast melting at the large scales. We do so by controlled melting experiments and corresponding direct numerical simulations, covering four orders of magnitude in scale. The transition corresponds to the transition from a laminar-type boundary layer around the melting object to a turbulent-type boundary layer, i.e., from so-called classical turbulence to ultimate turbulence, with its enhanced transport properties. Our results thus provide a quantitative understanding of the flow physics of the melting process and thereby enable a better extrapolation and prediction of melt rates on large scales such as relevant in geophysics, oceanography, and climate science.

研究の動機と目的

• 融解速度が物体サイズと周囲の乱れ（固定・自由輸送物体）にどうスケールするかを解明する。
• 乱流境界層の異なるレジーム間のすきゃんの転移が融解速度に影響するかを特定する。
• Reynolds 数および流れパラメータの依存性と転移点を定量化する。
• 地球物理学的・エネルギー貯蔵への適用のため、ラボ規模から地球規模へ融解速度を外挿する枠組みを提供する。

提案手法

• 統計的に定常な均質等方乱流中で、対象サイズを変えて制御された氷球融解実験を行う。
• immersed-boundary 法を用いた直接数値シミュレーションでHITにおける融解を追跡し、ラグランジュ融解を含む。
• ステファン条件を用いて融解率を解析し、水側対流と関連づける。
• Eulerian（固定）と Lagrangian（自由に輸送）融解ケースを比較し、境界層効果を評価する。
• 先行研究のデータを再整理・再分析して、Re_D0 の4桁オーダーにわたるスケーリング傾向を確立する。

Figure 1: Present experimental and numerical setups (a) Experimental setup: twenty independently-controlled electric motors mounted on the vertices of a $\approx$ \qty 210 dodecahedral water tank drive propellers to generate HIT in the tank center. (b) Volume contours of instantaneous temperature in

実験結果

リサーチクエスチョン

• RQ1 小規模（古典的）融解と大規模（究極的）融解における融解率（Nu）とサイズ基準のレイノルズ数（Re_D0）とのスケーリングはどうなるか。
• RQ2 乱流中の融解速度は層流型と乱流型境界層の転移によって決まるのか。
• RQ3 転移点は Re_D0 または Re_λ の観点で、プラントル（Pr）・ステファン数・流れ条件にどう依存するのか。
• RQ4 Eulerian と Lagrangian の融解配置は融解率の挙動にどの程度差をもたらすか。
• RQ5 地球物理学的・産業的文脈で、非常に異なるスケールにわたる融解速度を統一的に予測できるモデルは存在するか。

主な発見

• 融解は二つのレジームを示す：小規模Re_D0では Nu ∝ Re_D0^{1/2}（層流様境界層）、大規模Re_D0では Nu ∝ Re_D0^{γ}、0.8 ≤ γ ≤ 1（乱流境界層）。
• レジーム間の遷移は Re_D0 ≈ 4000 付近で起き、境界層せん断レイノルズ数は約60に相当。
• 究極的レジームでは Nu ∝ Re × L(Re) に従い、対数補正を取り入れることで、調査範囲で実質的に Nu ∝ Re_D0^{0.8} を与える。
• 融解は Re_D0（物体スケール渦）に強く依存するが、Re_λ（全体流れの渦）にはあまり敏感でなく、局所渦相互作用が融解を支配することを示す。
• Eulerian 融解は Eulerian ケースで表面のせん断率が高いため、Lagrangian よりやや速い融解速度をもたらす。
• 本研究は Re_D0 の4桁オーダーにわたるデータを統合し、実験・DNS・文献データを組み合わせて転移とスケーリング法則を確立した。

Figure 2: Reynolds number controls melting rate of ice balls in both experiments and DNS Snapshots of a freely-advected melting ice ball at two different times at fixed $\textit{Re}_{\lambda}$ and fixed $\textit{Re}_{D_{0}}$ (as noted) and temporal evolution of the volume-equivalent radius for vario