《常用的离散格式》PPT课件

《《常用的离散格式》PPT课件》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、2.4常用的离散格式在使用有限体积法建立离散方程时，很重要的一步是将控制体积界面上的物理量及其导数通过节点物理量插值求出。引入插值方式的目的就是为了建立离散方程，不同的插值方式对应于不同的离散结果。插值方式常称为离散格式。1离散格式假设速度场已知，则为求解离散方程，需计算广义未知量在边界e和w处的值。为完成这一任务，必须决定界面物理量如何通过节点物理量的插值表示。各种不同的插值方法就构成了不同的离散格式。中心差分格式一阶迎风格式混合格式指数格式乘方格式22.4.1术语的约定对离散格式的讨论以一维稳态对流扩散方程为例，不涉及瞬态项。定义两个新的物理量：F：通过界面上单位面积的对流质量通量，简称

2、对流质量流量；D：界面的扩散传导性在此基础上定义一个非常重要的参数Peclet数3Peclet数(佩克莱数)1976年Roache提出，网格或单元Peclet数可以用来度量某点处φ的对流和扩散的强度比例。从物理意义来说，Peclet数表示的是对流与扩散强度之比。而Re雷诺数表示作用于流体微团的惯性力与粘性力之比。Pr数是动量扩散厚度与热量扩散厚度之比。(粘性与导热性之比）佩克莱特数(Peclet)是普朗特数(Prandtl)与雷诺数(Reynolds)的乘积，实际上反映了惯性力与紊动强度的比值。4PeE(x)PeECentraldifferencingscheme中心差分格式Wede

3、terminethevalueofatthefacebylinearinterpolationbetweenthecellcenteredvalues.就是界面上的物理量采用线性插值公式来计算。具有二阶精度，扩散占优时计算具有较高精度。但流动为强对流时，计算的收敛性和精度较差。interpolatedvalue5中心差分格式的特点及适用性Pe小于2时，中心差分格式的计算结果与精确解基本吻合。Pe大于2时，数值解完全失去物理意义。原因：离散方程系数出现负值，负的系数会导致物理上不真实的解。解决办法：将计算网格加密，使网格间距变小，降低Pe数的大小。（计算量随之增大）基于此限制，中心差分格式不

4、能作为对于一般流动问题的离散格式，需创建其它更合适的格式（对纯扩散稳态，如热传导是适用的）。6对流扩散方程的精确解7精确解随Pe数的变化（Pe＝0纯扩散，Pe增大对流增强）8具体算例（不同计算工况意味着不同Pe数）9第一种工况Pe=0.2尽管网格粗糙，但数值解与精确解非常接近。10第二种工况Pe=5数值解在精确解周围振荡，计算精度不可接受。离散方程中系数有负值。11第三种工况（加密网格Pe=1.25）12Firstorderupwindscheme界面上的未知量恒取上游节点的值PeE(x)PeEFlowdirectionThisisthesimplestnumericalscheme

5、.Itisthemethodthatweusedearlierinthediscretizationexample.Weassumethatthevalueofatthefaceisthesameasthecellcenteredvalueinthecellupstreamoftheface.Themainadvantagesarethatitiseasytoimplementandthatitresultsinverystablecalculations,butitalsoverydiffusive.Gradientsintheflowfieldtendtobesmearedout,as

6、wewillshowlater.Thisisoftenthebestschemetostartcalculationswith.interpolatedvalue在确定界面物理量时考虑了流动方向13一阶迎风格式的数学描述及特点考虑了流动方向影响具有一阶截差unconditionallystable(nooscillations)（任何条件下都不会引起解的振荡，离散方程系数总是大于零）不足：夸大了扩散项的影响，numericallydiffusive14centraldifferencescheme:upwindscheme:1516171819混合格式Centraldifferencings

7、chemeismoreaccuratethanthefirstorderupwindscheme,butitleadstooscillationsinthesolutionordivergenceifthelocalPecletnumberislargerthan2.ItiscommontothenswitchtofirstorderupwindincellswherePe>2.Suchanapproachi