等效比热容法计算相变 该用哪个平滑函数？？？

The_Dawn

大家帮忙看下，在COMSOL中用等效比热容法计算相变该用哪个平滑函数？
1、C=Cp1+L*fldc2hs(T-Tm,dT), 貌似《CM多物理场建模与分析》这本书及天乐树的blog上都是这么说的。但我用这个表达式计算，不收敛啊。
2、C=Cp1+L/dT*flc2hs(T-Tm,dT), 利用这个表达式定义比热容就能收敛，其它设置和 上面相同。但是观察温度随时间变化曲线发现，从峰值温度下降呈分段的情况。这跟实际冷却过程不符啊，应该是连续冷却的过程啊。
各位大侠帮帮忙，感激不尽

cobra

这个是可行的~楼上可以参考网页http://www.cntech.com.cn/faq/comsol/2009/1231923737d3255.html

The_Dawn

cobra 发表于 2013-3-2 12:21
这个是可行的~楼上可以参考网页http://www.cntech.com.cn/faq/comsol/2009/1231923737d3255.html ...

这个我看过了。
中方论坛上貌似是第二种函数表式，但看过《多物理场建模与分析》和相关搜索结果，都是使用第一种表式。
观察这两个函数fldc2hs和flc2hs。是有明显区别的，前者是突变后回归到原始值，后者是阶跃的
考虑至此，等效热熔函数因为前者才是，对吗？

cobra

The_Dawn 发表于 2013-3-5 07:55
这个我看过了。
中方论坛上貌似是第二种函数表式，但看过《多物理场建模与分析》和相关搜索结果，都是使 ...

我的理解是，flc2hs能保证二次导连续，一般比较推荐这个。另一个我没有使用，不清楚，不过，应该区别就是不保证二次导连续。这两个函数都是阶跃的处理方法吧。。。

cobra

The_Dawn 发表于 2013-3-5 07:55
这个我看过了。
中方论坛上貌似是第二种函数表式，但看过《多物理场建模与分析》和相关搜索结果，都是使 ...

楼上可以加我q，注明名字、专业和学校/公司哈~一起探讨哦~q消息给你拉~虽然最近不一定有空，在忙毕业论文

Dreamcatcher927

你好，C=Cp1+L*fldc2hs(T-Tm,dT)和C=Cp1+L/dT*flc2hs(T-Tm,dT),，为什么潜热一个除了dT,一个没有除呢？C=Cp1+L*fldc2hs(T-Tm,dT)这里没有除dT的话，前后单位应该不一致吧？L*fldc2hs(T-Tm,dT)的单位是J/kg,和比热容的单位不一致啊？很困扰啊？谢谢您的解答~~

The_Dawn

Dreamcatcher927 发表于 2013-3-29 18:47
你好，C=Cp1+L*fldc2hs(T-Tm,dT)和C=Cp1+L/dT*flc2hs(T-Tm,dT),，为什么潜热一个除了dT,一个没有除呢？C=Cp ...

应该要除dT的

ghostlet

Dreamcatcher927 发表于 2013-3-29 18:47
你好，C=Cp1+L*fldc2hs(T-Tm,dT)和C=Cp1+L/dT*flc2hs(T-Tm,dT),，为什么潜热一个除了dT,一个没有除呢？C=Cp ...

fldc2hs是flc2hs的导数，所以不需要除dT。
纯粹从潜热概念来说，第一个式子更合适。因为相变之后热容会变为液态的值（此处默认假设固液两相相等而已）。