[論文レビュー] Homogenization of the Stefan problem, with application to maple sap exudation

Isabell Graf, John M. Stockie|arXiv (Cornell University)|Nov 12, 2014

Advanced Mathematical Modeling in Engineering参考文献 22被引用数 3

ひとこと要約

本稿では、周期的微細構造における熱伝達および相変化をモデル化するため、2スケール収束を用いた周期的均質化フレームワークを構築する。まず、融解するアイスバーを含む簡略化されたステファン問題に適用し、次に複雑なシロナガラスの樹液放出現象に適用する。主な貢献は、すべての微細構造物理学を参照セル問題に埋め込んだ、温度変化の有効なマクロスケールモデルの確立であり、これにより春の融雪期における木の樹液圧力生成の効率的数値シミュレーションが可能になる。

ABSTRACT

The technique of periodic homogenization with two-scale convergence is applied to the analysis of a two-phase Stefan-type problem that arises in the study of a periodic array of melting ice bars. For this "reduced model" we prove results on existence, uniqueness and convergence of the two-scale limit solution in the weak form, which requires solving a macroscale problem for the global temperature field and a reference cell problem at each point in space which captures the underlying phase change process occurring on the microscale. We state a corresponding strong formulation of the limit problem and use it to design an efficient numerical solution algorithm. The same homogenized temperature equations are then applied to solve a much more complicated problem involving multi-phase flow and heat transport in trees, where the sap is present in both frozen and liquid forms and a third gas phase is also present. Our homogenization approach has the advantage that the global temperature field is a solution of the same reduced model equations, while all the remaining physics are relegated to the reference cell problem. Numerical simulations are performed to validate our results and draw conclusions regarding the phenomenon known as sap exudation, which is of great importance in sugar maple trees and few other related species.

研究の動機と目的

• 周期的媒体における相変化を伴う2相ステファン問題の厳密な数学的フレームワークの開発。
• 液体、凍結、気体の樹液相を含むシロナガラスの木部における複雑な多相流および熱輸送のモデル化。
• 参照セル問題を通じてマクロスケール温度ダイナミクスを微細構造物理学から分離し、効率的な数値シミュレーションを可能にする。
• シミュレーションによるモデルの検証と、シロナガラスの木において樹液が放出される物理的メカニズムの解明。

提案手法

• 周期的均質化と2スケール収束を用い、水中に周期的配置のアイスバーを含む簡略化されたステファン問題に適用する。
• マクロスケール温度場と周期的単位胞 $Y$ 内の微細スケール参照セル問題を有する2スケール定式化を導入する。
• マクロスケール問題がグローバル温度を支配し、セル問題が局所的な相変化ダイナミクスを捉える、連成系を解く。
• 弱いおよび強い定式化を用いて、均質化方程式を数値的に取り扱えるアルゴリズムを導出する。
• エネルギー法およびグローヴァルの不等式を用いて、2スケール極限解の事前推定、存在、一意性、収束性を証明する。
• 同じ均質化モデルを、すべての微細構造物理学（多相流、浸透圧、気体効果）を参照セル問題に統合して、完全なシロナガラスの樹液放出問題に適用する。

実験結果

リサーチクエスチョン

• RQ1 温度依存の拡散係数を伴う2相ステファン問題に、周期的均質化を厳密に適用する方法は何か？
• RQ2 相転移を伴う周期的媒体における熱伝達および相変化の2スケール極限系の構造は何か？
• RQ3 均質化モデルは、計算複雑性を低減しつつ、シロナガラスの木における樹液放出の本質的物理をどのように捉えているか？
• RQ4 同一の均質化方程式を、簡略化モデルと完全な生物学的問題の両方に適用可能であり、物理的解釈が一貫しているか？
• RQ5 微細構造プロセスの物理的忠実性を保ちながら、均質化系を効率的に解く数値アルゴリズムは何か？

主な発見

• エネルギー推定およびグローヴァルの不等式を用いた証明により、2スケール極限解は存在し、一意的かつ適切な関数空間において強く収束することが示された。
• 均質化マクロスケール問題は有効な温度場を支配するが、すべての微細構造的相変化ダイナミクスは周期的単位胞 $Y$ 上の参照セル問題に符号化されている。
• モデルは、木部における完全な多相流問題の複雑さを、グローバル温度の1組の式に簡略化し、微細構造効果はセル問題に統合している。
• 数値的シミュレーションにより、モデルが春の融雪期における圧力上昇といった樹液放出の主要な特徴を再現できることを確認した。
• 細かい細胞構造を解像しないで済むため、このアプローチは、全生物レベルでの生物学的現象の研究に適している。
• 均質化方程式は、簡略化モデルと完全なシロナガラスの樹液放出問題の両方で一貫性があり、ロバストであることが確認され、この手法の一般性が裏付けられた。