有无大佬帮忙解决一下问题，我现在在模拟石蜡的凝固融化，一开始石蜡的热物性参数用的是常数，模拟的结果就很收敛，趋势也正确，但是因为固液的热物性是不一样的，一把石蜡的密度、导热率和比热，只要有一个改成分段线性定义，模拟很快就发散；后面密度用Boussinesq，导热和比热用分段线性定义也不发散了，但是只有在石蜡热膨胀系数设为0的情况下才行，不设为0也发散，但是设为0的话密度不变就看不到自然对流对融化的影响了，不知道是什么原因导致的

你好 我现在和你一样的问题 怎么解决的啊