利用计算机进行数值模拟计算,数值模拟法

差分法

差分法是把原来求解物体内随空间、时间连续分布的温度问题,转化为求在时间领域和空间领域内有限个离散点的温度值问题，再用这些离散点上的温度值去逼近连续的温度分布。差分法的解题基础是用差商来代替微商，.这样就将热传导微分方程转换为以节点温度为未知量的线性代数方程组，得到各节点的数值解。用不同方法定义差商可得到不同的差分格式:向前差分、向后差分、平均差分、中心差分、加列金格式等。不同的差分格式其误差和稳定性各不相同。如向前差分计算过程简便，然而是有条件稳定的;向后差分则是无条件稳定的;而平均差分虽然精度较高但容易发生振荡等。因此，使用差分法时要选择合理的差分格式、合理的网格划分和计算步长，以尽可能减少误差，保证解的精度和稳定性。差分法的长处是对于具有规则的工件外形和均质材料的温度场求解，它的程序设计和计算过程比较简单，收敛性也较好。

2.有限元法

有限元法是根据变分原理来求解热传导问题微分方程的一种数值计算方法。有限元法的解题步骤是先将连续求解域分割为有限个单元组成的离散化模型，再用变分原理将各单元内的热传导方程转化为等价的线性方程组，最后求解全域内的总体合成矩阵。由于有限元法的单元形状可以比较任意，因此更能适合于具有复杂形状的物体。对于由几种不同材料组成的物体，可以利用不同材料的界面进行单元分割。特别是可以根据实际问题需要设置单元和节点分布的稀疏，这样就可以在不增加节点和计算量的前提下提高计算精度。此外，由于有限元法是用统一的方法对区域内节点和边界节点列出计算格式，因此适用于任意的边界条件，且在计算精度上比较协调。

有限元法是随着计算机技术发展起来的一种新颖、有效的数值计算方法，它不仅可以用于传热分析，而且用于电磁场、流体力学、结构力学等所有连续介质和场的问题的分析采用有限元法计算温度场时便于与后续相关计算统一计算方法，如铸造残余应力分析，凝固裂纹倾向与凝固组织分析等。

20世纪80年代以来，基于有限元法的计算机数值模拟技术在世界各国先后得到突飞猛进的发展，陆续研究出用于解决各种工程问题的多种有限元分析商业软件，包括：用于液态金属凝固过程的温度场、应力应变场的数值模拟;凝固组织、凝固缺陷分析与性能预测等。

总的来说，差分法和有限元法已成为解决温度场计算问题的两种主要计算方法，它们各有所长。在实际计算过程中,经常将这两种方法结合起来进行计算，例如在空间域采用有限元法进行离散,时间域采用差分法进行离散。