第三章-传热学数值计算方法技巧.ppt

本章任务：前两章已阐述了获得物理现象的理论预测有一些明显的好处，且一些有趣的现象是受一些微分方程支配的，这些方程已用通用方程归纳了，接下来的任务就是推导求解这个通用方程的方法—即所谓的离散化方法。;§3.1 数值方法的本质;eg. 设? 随 x 的变化遵循高次多项式的形式; ;2. 离散化方程：所取网格节点上未知因变量 ? 值的代数方程，此方程由支配 ? 的微分方程推导而得。;§3.1-3 离散化方程的结构;2. 离散化方法;§3.2 推导离散化方程的方法;取左端及右端的前三项，并进行相加或相减，便可得中心差分的近似式：;为了提高精度，可以得到截断误差更高阶的差分表达式。;§3.2-2 用多项式拟合法建立导数的差分表达式;2. 采用二次曲线拟合，可得到具有二阶精度的空间导数;eg: 如图所示，已知区域内部与边界节点的温度，物体的导热系数λ=const。试用多项式拟合法确定穿过壁面的热流密度。;由上面三式可解得：;§3.2-3 控制容积公式;有限差分法，仅当网格极其细密时，离散方程才满足积分守恒； 而有限体积法,即使在粗网格情况下，也显示出准确的积分守恒。 有限差分法,只考虑网格点上的数值而不考虑值在网格点之间如何变化。 有限体积法在寻求控制体积的积分时，必须假定值在网格点之间的分布，插值函数只用于计算控制体积的积分，得出离散方程之后，便可忘掉插值函数； 上述观点，可使得对微分方程中不同的项采取不同的插值函数有完全的自由。;3. 积分法实施的步骤;§3.3 一个说明性的例子;分布函数：通常有阶梯式分布和分段线性分布两种;分段线性分布:变量在网格节点间呈线性分布，可以用来计算变量的梯度，有时也用于计算变量在时域上的分布。;说明：;二、指导原则;x;②.总平衡的要求;应用线性化的源项表达式，离散化方程的形式一样，但系数有所改变。;四、不同离散方法的比较;③.控制容积法(平衡法) 优点：导出过程物理概念清晰，离散系数具有一定的物理意义，并可以保证所得方程具有守恒特性。 缺点：不便于对方程进行数学特性的分析。 ④.这两种方法分别展示了有限差分法与有限元法这两种数值解法的基本特点，有限容??法更具有吸引力。;§3.4 四项基本法则; 如图所示，在P、E两个控制容积的公共界面e上，离开P 控制容积穿过e界面的q、qm及动量通量? 应各自等于穿过e界面进入E 控制容积的相应的量。 若公共界面e上的型线选择不妥，可能导致界面上连续性受到破坏 如上图所示，界面上采用了二次曲线，由于从P控制容积及E 控制容积来确定的二次曲线拟合点不完全相同，在e界面上，无论是?值还是其导数从P、E两侧控制容积确定的值均不相等，使得格式失去守恒性。;①;若公共界面e上的导热系数选择不妥，也可能使界面上连续性被破坏;从物理意义上看，若从界面的两侧计算所得通量不能相互抵消，相当于在界面上存在一个由计算而造成的源或汇，使总体计算误差增加。;法则2：正系数 ;法则3：源项的负斜率线性化;法则4：相邻节点系数之和规则;法则4应用条件：