有限体积法详解.ppt

输运性 网格Peclet数大，上游节点变量值对下游影响大，下游对上游影响小； Pe =0，上下游影响一样。 中心差分格式使节点P处场变量对所有相邻节点一样，没有反映出扩散和对流的差别，不能体现输运方程的方向性，在高Peclet数数时，中心差分格式不具有输运特征。 中心差分格式离散方程特点 （1）守恒性：满足； （2）有界性：Pe小于2时满足，不然不满足； （3）输运性：没有； （4）计算精度：二阶精度， Pe小于2时精度较高， Pe大于2或流动为强对流时，收敛性和精度均较差。 双曲型方程特征值体现了微分方程的解（或者说扰动、信息）的传播方向，表达了波动、能量等的传播方向，应该充分考虑这一物理特性，与之相适应。 恒取上游节点处的值建立差分格式。 一阶上风差分格式 对无源扩散问题： 只对对流项的差分格式进行改变 离散方程为： 解: 2，3，4点： 边界： 整理： 工况1 工况2 显示上风差分格式考虑了流动的方向性，在有较强对流输运状况时具有计算优势。 上风差风格式性质 （1）守恒性：满足； （2）有界性：满足，计算结果不会出现振荡或不收敛； （3）输运性：考虑了流动方向，满足； （4）计算精度：一阶精度；同时扩散项采用中心差分格式，随着Pe数增大，对流输运增强，扩散输运减弱，所以，当Pe足够大，仍然保持不变的扩散输运强度，必然带来误差，称为假扩散。 混合差分格式 Spalding提出：当网格Pe数小于2时采用中心差分格式计算控制容积界面值，具有二阶精度；当Pe数大于2时采用上风差分格式计算控制容积界面对流输运量，同时忽略扩散输运量，计算精度一阶，但能反映流动的输运特征。 工况2 上风格式有扩散 混合格式无扩散 QUICK格式 上风格式和中心差分格式的优缺点； 1979年Lenoard提出用于计算控制容积界面的二次插值计算格式。“对流项的二次上风插值”的英文缩写（Quadratic Upwind Interpolation for Convective Kinematics）. Quick格式利用控制容积界面两侧3个节点值进行插值计算。其中两个节点位于界面的紧邻两侧，另一个节点位于上风侧远邻点。WW， W， P的拟合曲线计算 ，通过W，P，E计算 QUICK格式特性 1）守恒性：满足； 2）有界性：有条件满足， Pe大于8/3时有可能不满足； 3）输运性：总采用上风两个节点，下风一个节点进行插值，所以满足； 4）精度：具有3阶截差，精度较高，要保证Pe小于8/3。 在用有限体积法计算时，和有限差分法一样，方程的解是用单元节点上离散点值构成的，而不关心单元间的状态变量是怎么变化的，也就是不关心解的分布。 有限单元法中一旦选定了分布曲线，就确定了状态变量的分布函数是不同的。 尽管我们在离散方程的时候要假定单元分布曲线，但是这只是为了推导公式时计算积分近似而采用的一些辅助关系式。一旦离散化方程推导出来了，就可以不用再管这些分布曲线近似。 在积分离散时，根据数值模拟的需要，对控制方程中的每一项都可以采用不同的分布曲线来近似单元界面上的状态变量，而不必要追求近似假设的一致性。 方程离散1 一维稳态扩散问题 按节点整理后 得到: 最后解方程组得到节点的值 例题 解: 1点控制容积积分 方程离散1 二维扩散 一维稳态对流扩散问题 方程离散2 稳态对流扩散方程 连续性方程 对控制体积分 定义: 解: 对中间节点2,3,4: 边界节点1: 整理得到: 边界节点5: 整理得到: 工况1 改进办法:需要增加网格数 工况2 工况3 差分格式问题 中心差分格式? 在扩散方程中,未引起数值解和分析解之间较大的差别,但在包含对流项的对流扩散方程中,某些计算条件下需要加密网格结果才合理. 边界的计算形式? 差分格式问题 控制容积界面处变量的近似计算格式是否对数值计算结果有影响？ 采用其他差分格式是否能提高计算精度？ 差分格式近似计算式在流场计算中的物理意义是什么？ 差分格式应该满足三个特性：守恒性、有界性和输运性。 守恒性 如果对一个离散方程在定义域的任一有限空间内作求和的运算（相当于连续问题中对微分方程作积分），所得的表达式满足该区域上物理量守恒的关系时，则称该离散格式具有守恒特征。 有限体积法正是从物理量守恒这一基本要求出发提出的。有限体积法的离散化方程满足了单个控制体积的平衡，当然在整个计算区域内，诸如质量、动量等物理量的积分守恒也就都能精确得到满足。无论在数值计算中采用巨大数目的细网格和少数的粗网格，数值解也照样显示准确的积分平衡。有限体积法的离散思想自动满足守恒定律，如质量守恒，动量守恒，能量守恒等等。所以有限体积法是守恒定律的一种最自然的表现形式。 有界性 离散方程为代数方程组，求解时需要用迭代方法，获得收敛解。 扩散问题