基于焓变守恒的相变过程变比热问题温度解析修正.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号： 基于焓变守恒的相变过程变比热问题温度解析修正 穆磊，贺志宏，马宇，董士奎 （哈尔滨工业大学航空航天热物理研究所，黑龙江哈尔滨市150001） TEL: 0451 E-mail: dongsk@hit.edu.cn 摘要：物质相变过程中，比热随相态转变和潜热释放而改变，而比热的变化又对温度产生影响，使得 求解相变过程正确的温度分布成为一件困难的事情。为解决此问题，许多学者提出了不同的解决办法， 但基本上都是基于相变过程温度和比热相互迭代的思想。本文在研究前人工作的基础上提出了基于焓 变守恒的相变过程变比热问题的温度解析修正法，给出了详细推导过程，并采用SPH方法模拟计算了 角域凝固问题，通过和解析解的比较得出，该修正方法物理意义清晰、计算精度高、计算速度快，能 够较好的模拟相变过程温度分布的变化。 关键词：相变过程；变物性；焓变守恒；温度修正方法；SPH模拟 0 引言 固液相变过程中比热随相态转变而变化，比热的改变将影响温度的分布，同时在固 液相变过程中潜热的吸收和释放也会影响温度的分布，因此准确计算固液相变过程中的 温度分布是一件很困难的事情。尤其是在采用等效热容概念计算固液相变过程中的温度 分布时，每一个时间步内所使用的比热通常是该时间步初始时刻的数值，且在时间步内 将比热处理为常数，但由于相态转换和潜热释放使得时间步内比热变化是较为剧烈的， 处理为常数会导致时间步终了时的计算温度偏离真实值，需要对该温度进行修正。已有 文献中采取时间步内迭代、温度回升[1] [2] 和时间后退法 进行处理，以上方法在求解过程 中利用了时间步内温度和比热相互迭代的思想来准确的求解温度分布，这会耗费大量的 计算时间。本文为提高计算效率，避免时间步内迭代，提出基于焓变守恒的相变过程变 比热问题的温度解析修正方法，给出了详细的推导过程，并结合SPH法对角域凝固问题 进行了模拟计算，验证了本文方法的正确性。 1 相变传热主控方程 根据傅里叶定律可以得到多维非稳态无内热源导热方程一般式为： dT ρc =∇⋅ λ∇( T ) (1) dτ 采用(1)式描述固液相变过程时，通常将固相和液相的物性参数处理为常数，而相变 区的物性参数采用线性化处理[3] 。若如此，固相区的物性参数用角标 表示，液相区的 s 物性参数用角标l 表示，则可将(1)式改写为以下两个方程： ∂T 在固相区有： c =∇⋅ ∇T (2) ρs s (λs ) ∂τ ∂T 在液相区有： c =∇⋅ ∇T (3) ρl l (λl ) ∂τ 在实际计算时为了表达式的统一，无论是固相区还是液相区采用统一的角标 表