计算传热学 对流扩散方程的离散格式精选文档.ppt

计算传热学 对流扩散方程的离散格式本讲稿第一页，共三十三页主要内容n n原始变量法原始变量法n n两个问题两个问题n n交错网格交错网格n n交错网格的压力修正方法交错网格的压力修正方法n n同位网格的压力修正方法同位网格的压力修正方法本讲稿第二页，共三十三页6.1 原始变量法n原始变量非原始变量n基于密度基于压力n耦合求解分层求解本讲稿第三页，共三十三页6.1 原始变量法本讲稿第四页，共三十三页6.1 两个问题n 压差n如何求解压力本讲稿第五页，共三十三页6.1 两个问题n 压差本讲稿第六页，共三十三页6.1 两个问题n 压差本讲稿第七页，共三十三页6.1 两个问题n 压差棋盘型压力场棋盘型压力场本讲稿第八页，共三十三页6.1 两个问题n压力如何求解二维问题：二维问题：u,v,p,T,四个未知量四个未知量四个方程：连续，动量，能量四个方程：连续，动量，能量压力没有独立方程，如何求解？压力没有独立方程，如何求解？本讲稿第九页，共三十三页6.1 两个问题n压力如何求解二维问题：二维问题：u,v,p,T,四个未知量四个未知量四个方程：连续，动量，能量四个方程：连续，动量，能量压力如何求解？压力如何求解？压力作用？压力作用？本讲稿第十页，共三十三页6.1 交错网格n 压差的处理交错网格引入交错网格引入1-压差压差本讲稿第十一页，共三十三页6.1 交错网格n 压差的处理u，v的控制容积P的控制容积本讲稿第十二页，共三十三页6.1 交错网格n动量方程的离散本讲稿第十三页，共三十三页6.1 交错网格n界面参数插值本讲稿第十四页，共三十三页6.1 压力修正方法n压力修正法的基本思想压力没有独立方程迭代修正的思想利用连续方程u假定压力场，假定压力场，P*u利用压力场求解动量离散方程，得出利用压力场求解动量离散方程，得出u*，v*u利用连续方程构造相应的方法，改进压力场，要求改进的压力利用连续方程构造相应的方法，改进压力场，要求改进的压力场能够满足连续性方程，修正值场能够满足连续性方程，修正值 p,u,v本讲稿第十五页，共三十三页6.1 压力修正方法n压力修正法的基本思想u假定压力场，假定压力场，P*u利用压力场求解动量离散方程，得出利用压力场求解动量离散方程，得出u*，v*u利用连续方程构造相应的方法，改进压力场，要求改进的压利用连续方程构造相应的方法，改进压力场，要求改进的压力场能够满足连续性方程，修正值力场能够满足连续性方程，修正值 p,u,vuu*+u,v*+v,p*+p作为新的解，开始下一层次的计算作为新的解，开始下一层次的计算如何获得如何获得p？如何在如何在p的基础上获得的基础上获得u，v本讲稿第十六页，共三十三页6.1 压力修正方法n压力修正法的基本思想如何获得如何获得p？如何在如何在p的基础上获得的基础上获得u，v只有一个连续方程，因此，只能求解一个变量，需要把只有一个连续方程，因此，只能求解一个变量，需要把u，v转换成转换成p的形式的形式本讲稿第十七页，共三十三页6.1 压力修正方法n速度修正值的获得p的基础上获得的基础上获得u，v本讲稿第十八页，共三十三页6.1 压力修正方法n速度修正值的获得速度修正值速度修正值u以上方程可以求解，但是比较复杂以上方程可以求解，但是比较复杂忽略相邻节点的修正值影响，只保留压力修正值的影响忽略相邻节点的修正值影响，只保留压力修正值的影响本讲稿第十九页，共三十三页6.1 压力修正方法n速度修正值的获得速度修正值速度修正值u，v改进后的速度为改进后的速度为P如何获得呢？如何获得呢？本讲稿第二十页，共三十三页6.1 压力修正方法n压力修正值方程从连续方程出发，构造压力修正方程从连续方程出发，构造压力修正方程控制容积积分法离散得到控制容积积分法离散得到本讲稿第二十一页，共三十三页6.1 压力修正方法n压力修正值方程压力修正方程压力修正方程本讲稿第二十二页，共三十三页6.1 压力修正方法n压力修正值方程压力修正方程的源项压力修正方程的源项b=0意味着什么？意味着什么？本讲稿第二十三页，共三十三页6.1 压力修正方法n压力修正值方程压力修正方程的边界条件压力修正方程的边界条件不必引入边界上的压力修正值信息，设置相应的系数为不必引入边界上的压力修正值信息，设置相应的系数为0即可即可本讲稿第二十四页，共三十三页6.1 压力修正方法nSIMPLE算法的步骤SIMPLE：压力耦合方程的半隐方法：压力耦合方程的半隐方法本讲稿第二十五页，共三十三页6.1 压力修正方法nSIMPLE算法的发展SIMPLER（构造了一个压力方程，求解作为（构造了一个压力方程，求解作为p*）SIMPLECSIMPLEXSIMPLEST本讲稿第二十六页，共三十三页6.1 同位网格上的SIMPLE算法n同位网格交错网格数据存储不方便交错网格数据存储不方便同位网格变量均存储在节点上同位网格变量均存储在节点上如何解决前面所述的两个问题？如何解决前面所述的两个问题？本讲稿第二十七页，共三十三页6.1 同位网格上的SIMPLE算法n同位网格动量方程的离散形式动量方程的离散形式本讲稿第二十八页，共三十三页6.1 同位网格上的SIMPLE算法n同位网格连续方程的离散形式连续方程的离散形式界面上的速度如何处理呢？界面上的速度如何处理呢？在这里引入了在这里引入了1-deta压差压差本讲稿第二十九页，共三十三页6.1 同位网格上的SIMPLE算法n同位网格同位网格同位网格的同位网格的SIMPLE算法思想与交错网格相同算法思想与交错网格相同u，v，p的关系类似于交错网格的关系类似于交错网格ue，uw等的表达式中的系数需要插值等的表达式中的系数需要插值本讲稿第三十页，共三十三页6.1 同位网格上的SIMPLE算法n同位网格同位网格u，v，p的关系的关系采用略去邻点影响的方法采用略去邻点影响的方法本讲稿第三十一页，共三十三页6.1 同位网格上的SIMPLE算法n同位网格同位网格压力修正方程的获得压力修正方程的获得代入代入本讲稿第三十二页，共三十三页6.1 同位网格上的SIMPLE算法n同位网格同位网格计算步骤计算步骤与交错网格类似，注意动量插值的使用与交错网格类似，注意动量插值的使用本讲稿第三十三页，共三十三页