采用焓-多孔介质方法模拟相变材料熔化过程交流讨论
-
各位Foamer大家好，很高兴获得中文网的邀请码和大家进行交流。 首先说下我的情况：我的课题是关于电场作用下相变材料熔化过程研究，前期课题在学校主要以实验为主。但受疫情影响，目前在家采用openfoam模拟电场力作用下流动传热问题。前期在看了有限体积法基本原理、压力速度耦合求解等策略后，目前能够在动量方程中植入电场力来模拟单相液体流动传热问题，处于入门进阶阶段。这个过程必须要感谢东岳老师在“东岳流体”官网上提供的基本求解器解析及其他素材，使我入门在openfoam入门过程少走了很多弯路，再次对东岳老师及提供的平台表示感谢。
言归正传，我目前在进行采用焓-多孔介质模型模拟恒定密度下纯物质熔化过程。目前我所知道的相变材料熔化过程能量方程可分为求解温度方程或求解焓方程。在求解温度的能量方程中，由于温度和液相分数存在非线性的耦合关系，Voller等人提出了对温度和液相分数进行迭代的方法（文献1：A General Enthalpy Method for Modeling Solidification Processes），目前openfoam中的fvoption/solidificationMeltingSource 应该基于这个原理；Fabian等人通过引入误差函数来构建液相分数和温度的关系（文献2：Shell-and-tube type latent heat thermal energy storage: numerical analysis and comparison with experiments），同时基于误差函数的求解方法及模拟相变材料熔化的求解器在cfd-online上也引起了广泛的讨论（https://www.cfd-online.com/Forums/openfoam-solving/93620-melting-problem-looking-appropriate-solvers-13.html#post652310）。
在通过求解焓方程方面，受到文献3（Physics of humping formation in laser powder bed fusion）的启发，我在openfoam中植入了求解焓的方程，其中基本控制方程见附图1
通过求解焓获得温度和液相分数。并采用开发的求解器模拟金属镓的熔化过程，其中镓选取和文献1相同的物性，得到的模拟结果中的总液相分数和文献1中的结果相差不大，其中文献1的数值模拟结果和实验进行过对比，其中文献1的结果总体相对可靠，经常当做benchmark，。但是附图2中的固液相分界面位置和文献1结果存在偏差，
也就是我得到的相界面上部较文献1落后，底部熔化快于文献1，对于出现该问题的原因我至今还没找到，因此特开此帖和大家进行交流讨论，最后附上我的solver和case，
热切希望大家能给出一些建议。PS:本帖文风略显正式，希望大家不要见怪哦~
-
我觉得你这个结果已经非常接近了，可能不是算法上有大的遗漏，可能是一些边边角角，格式啊网格啊边界条件啊什么的问题？
-
首先感谢东岳老师给出的建议。经过一番比较后，发现数值离散格式（upwind VS quick）、液相与固相温差（Tl-Ts）果然会对固液相分界面的位置相对明显的影响。其中，对于纯金属镓的熔化模拟，尽可能选择较小的固液相温差，如0.025k。同时这也充分说明了搞CFD真的靠经验~ @东岳
-
同时对于大多数物质在熔化和凝固过程体积均会发生变化，即固相和液相密度不同，现在一方面参考东岳老师的buoyantPimpleFoam的解析，即需要植入密度求解方程，一方面参考文献：An Optimum Enthalpy Approach for Melting and Solidification with Volume Change。我在想如何在openfoam中植入这一过程，期待以后能跟东岳老师多学习交流~
-
挺好挺好，
buoyantPimpleFoam的解析那个我还没更新 https://www.cfd-china.com/topic/3665/关于buoyantpimplefoam的几个问题/12 rhoPimpleFoam目前更完善。
我这面带了一个学生目前打算做这方面工作，我还没找好具体的方向，不过也是类似的。熔化和凝固。我得先看看这方向是不是可以做的，然后再让他做
-
@沐风 您好! 我最近刚开始学习OF。关于您的 文献1的validation，我也想试一下，请问在什么标准求解器上修改比较好？想要应用boussinesq 假设，但看buoyantBoussinesqPimpleFoam 求解器描述为湍流的。但该算例的自然对流Ra=6*10^5,即为层流。希望您指教一下，非常感谢！！
