[论文解读] 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

一句话总结

本文利用两尺度收敛的周期均质化框架，对周期性微结构中的传热与相变进行建模，首先将该方法应用于熔化的冰棒简化版Stefan问题，随后扩展至糖枫树汁渗出这一复杂现象。主要贡献在于构建了一个用于温度演化的简化宏观模型，能够捕捉有效热行为，所有微观物理过程均编码于参考元胞问题中，从而实现对春季解冻期间树木中树汁压力生成的高效数值模拟。

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.

研究动机与目标

• 开发一个严格的数学框架，用于对具有相变的周期性介质中的Stefan问题进行均质化。
• 对糖枫树木质部中的复杂多相流与热传输进行建模，包括液态、固态和气态树汁相。
• 通过参考元胞问题，将宏观温度动力学与微观物理过程分离，从而实现高效数值模拟。
• 通过模拟验证模型，并揭示糖枫树中树汁渗出的物理机制。

提出的方法

• 将周期性均质化与两尺度收敛方法应用于一个简化的Stefan问题，涉及水中的周期性冰棒阵列。
• 引入两尺度公式，包含宏观温度场与周期单位胞 $Y$ 中的微观参考元胞问题。
• 求解一个耦合系统，其中宏观问题控制全局温度，而元胞问题捕捉局部相变动力学。
• 采用弱形式与强形式的极限问题，推导出求解均质化方程的数值可处理算法。
• 通过能量方法与Gronwall不等式，推导先验估计，并证明两尺度极限解的存在性、唯一性与收敛性。
• 将相同的均质化模型应用于完整的糖枫树汁渗出问题，通过参考元胞问题处理所有微观物理过程（多相流、渗透、气体效应）。

实验结果

研究问题

• RQ1 如何严格地将周期性均质化方法应用于具有温度依赖扩散系数的两相Stefan问题？
• RQ2 在具有相变的周期性介质中，传热与相变的两尺度极限系统具有何种结构？
• RQ3 该均质化模型如何在降低计算复杂度的同时，捕捉糖枫树中树汁渗出的本质物理机制？
• RQ4 相同的均质化方程是否可同时应用于简化模型与完整生物问题，并保持一致的物理解释？
• RQ5 何种数值算法可高效求解均质化系统，同时保持微观过程的物理保真度？

主要发现

• 通过能量估计与Gronwall不等式证明，两尺度极限解在适当的函数空间中存在、唯一且强收敛。
• 均质化后的宏观问题控制有效温度场，而所有微观相变动力学均编码于周期单位胞 $Y$ 上的参考元胞问题中。
• 该模型成功地将树木质部中全尺度多相流问题的复杂性，简化为一组控制全局温度的方程，而微观效应则由元胞问题捕捉。
• 数值模拟证实，该模型能够再现树汁渗出的关键特征，包括春季解冻期间的压力累积。
• 该方法可避免解析细观细胞结构，从而实现高效模拟，适用于在整体生物体水平研究生物现象。
• 均质化方程在简化模型与完整糖枫树汁渗出问题中均表现出鲁棒性与一致性，验证了该方法的普适性。