有限差分法解雷诺方程.docxVIP

有限差分法解雷诺方程 有限差分法解雷诺方程 PAGE 有限差分法解雷诺方程 雷诺方程的数值解法——有限差分法 雷诺方程的一般形式为 ? 1求解雷诺方程的方法有很多，有限差分法是最为常用的一种方法 2有限差分法的解题步骤： 将所求的方程量纲一化。目的是减少自变量和应变量的数目，同时，用量纲一话的解具有通性。 将求解域划分成等距或者不等距的网格，网格的划分根据精度要求来定。 将方程写成线性形式。 这一步是有限差分法的关键，主要是将二次导数和一次导数展开成一次项 例如变量?在（i，j）点的偏导数为 ?? 类似于f(x)= ?? ? ? 当然，在边界上时，采用单向差分： ?? ?? 或者向后差分公式： ?? ?? 所以，先将雷诺方程写成标准形式： A 其中A、B、C、D、E都是常数，将一次、二次导数的展开式带入到上式中可得： ? 其中 C C C C G=- 这样，我们得到了有限差分法的计算方程，对于每一个节点都可以写出这样一个方程，再利用边界上的点满足的边界条件 本次matlab算法基于稳态、等温、层流和不可压缩牛顿流体假设的Reynolds 方程为 ? 展开得到 h 将h3和?h 求出各点压强之后承载量 W= 摩擦力为流体的剪切力 F= 其中 τ=η 此时，我们取z=h/2，这样得到τ的平均值