第三章偏微分方程的数学性质对CFD的影响－.pptVIP

* 确定偏微分方程的一般方法（特征值法） 4 举例说明特征值法 B 特征值法 将方程改写成： 其中： 其中： 为求[K]-1，先写出[K]的代数余子 NEXT * 确定偏微分方程的一般方法（特征值法） 4 举例说明特征值法 B 特征值法 其转置矩阵还是 推出 * 确定偏微分方程的一般方法（特征值法） 4 举例说明特征值法 B 特征值法 求矩阵[N]的特征值λ，令 展开 和特征线法结果完全相同。 * 确定偏微分方程的一般方法（特征值法） 4 举例说明特征值法 C 小结 这里我们要指出： 有些方程组，特征值可能既有虚数又有实数。这些情况下，方程组既不是椭圆型的又不是双曲型的，他们表现出双曲型－椭圆型的混合特性。 * 不同类型偏微分方程的一般性质 双曲型方程一般性质 抛物型方程一般性质 椭圆型方程一般性质 * 双曲型方程一般性质 用CFD术语讲，由双曲线方程决定的流程可以“推进”求解。 * 流体力学中，下列类型的流动由双曲型偏微分方程决定的： 定常无粘超声速流动 非定常无粘流动 双曲型方程一般性质 * 抛物型方程一般性质 和双曲线方程一样，抛物型方程也适用于“推进”求解。 * 抛物型方程一般性质 流体力学中，下列流动模型的控制方程是抛物线方程： 定常边界层流动 “抛物化”粘性流动 非定常热传导 * 涉及椭圆方程的问题常常被称为“陪审团”问题，因为区域内的解依赖于所有的边界。 必须在整个边界上给定边界条件： 在边界上指定未知函数 u和v－Dirichlet条件 在边界上指定未知函数的导数， 椭圆型方程一般性质 * 椭圆型方程一般性质 流体力学中，下列流动模型的控制方程是椭圆型方程： 定常亚声速无粘流动 不可压无粘流动 * Nanjing University of Technology 南京工业大学机械与动力工程学院 凌祥 * 第三章 偏微分方程的数学性质对CFD的影响 拟线性偏微分方程的分类 确定偏微分方程的一般方法（特征值法） 不同类型偏微分方程的一般性质 定解问题的适应性 * 基本概念 记号 函数 偏微分方程 拟线性偏微分方程的分类 * 偏微分方程阶数 一阶 二阶 四阶 拟线性偏微分方程的分类 * 偏微分方程线性性质（一阶为例） 拟线性偏微分方程的分类 线性 拟线性 非线性 * 当一个n价偏微分方程的系数依赖于n阶导数时，此方程是非线性的 当系数依赖于m阶导数，而m＜n时，它就是拟线性的 方程的线性性质极为重要，因为对于线性和拟线性偏微分方程，它们的许多解析性质已被了解，而对于非线性偏微分方程则必须逐个地去研究它。 拟线性偏微分方程的分类 偏微分方程阶数线性性质 * 拟线性偏微分方程的分类 * 拟线性偏微分方程的分类 抛物型 椭圆型 双曲型 ? 如何得出的 * 拟线性偏微分方程的分类 考虑如下拟线性方程组 u和v是未知数，都是x、y的函数。 我们可以把u和v想象成xy平面的连续速度场 * 拟线性偏微分方程的分类 考虑xy平面的任意一点，如图中P点 写出u和v的全微分 系数矩阵 （1）、（2）、（3）、（4）组合成一个方程组： * 拟线性偏微分方程的分类 用克莱默法则求解方程（5）中的 一般情况下， 在P点值是不变的，与通过P点的方向无关。 然而，有个重要的例外： 从P点开始移动，此时 特征线 * 拟线性偏微分方程的分类 特征线 展开得 令 方程演化为 这样我们可以积分求得特征曲线y=y(x) * 拟线性偏微分方程的分类 通过P点得斜率： 判别式： 二次曲线是双曲线 二次曲线是抛物线 二次曲线是椭圆 * 确定偏微分方程的一般方法（特征值法） 上一节借助克莱默法则得到偏微分方程分类，本节讲另外一种判别方法－特征值法： （8） （9） 先假设上节方程（1）（2）中的 f1、f2为0： 定义列向量： 其中： * 确定偏微分方程的一般方法（特征值法） 利用 的特征值就可以确定方程组的类型 如果 的特征值均为实数—— 双曲型 如果 的特征值均为复数—— 椭圆型 * 二阶波动方程 1 双曲型方程 确定偏微分方程的一般方法（特征值法） 拟线性形式 特征方程 I为单位矩阵 * 1 双曲型方程 确定偏微分方程的一般方法（特征值法） 二阶波动方程 特征方程有两个实根 波动方程因此是双曲型，容易求出两个左特征向量为 * 1 双曲型方程 确定偏微分方程的一般方法（特征值法） 二阶波动方程 对应的特征相容关系为： 沿特征线 沿特征线 * 1 双曲型方程 确定偏微分方程的一般方法（特征值法） 二阶波动方程 引入变量置换： 其中f，g为任意可微函数，因此，波动方程的通解为： * 1 双曲型方程 确定偏微分方程的一般方法（特征值法） 考虑初值问题，初始条件为： 则波动方程的解为：