计算流体力学的基础P2偏微分方程的性质.pptVIP

Slide * 4.4.2 抛物型方程-例 3. 非定常热传导 假设流体的温度梯度是速度的函数、无附加的体积热，且内能e=cvT: K=const 一维情况： 热扩散率 Slide * 4.4.3 椭圆型方程 1.特征线为虚数，故与特征线有关的解法不适用； 2.无有限影响区域和依赖区域，流场参数信息可以向任何方向传播； 3.图中P点参数影响整个区域的信息，同时区域内任意点的参数也影响P点的参数。 边界条件的类型： Slide * 4.4.3 椭圆型方程-例 1. 定常亚音速无粘流动 2. 不可压缩无粘流动 典型的低马赫数流动 超音速钝体绕流问题的解决 Slide * 与尖头体相比，钝头体的气动加热明显降低。 20世纪50-60年代，最让人困惑的空气动力学问题之一。Why？同一个问题，出现了两个明显不同的区域，Ma1和1，此时偏微分方程的性质不同（抛物线和椭圆型），如何发展一致有效的方法呢？后来，布鲁克林工业大学的Moretti和Abbett发展了时间推进方法，解决了此难题。 * * 与尖头体相比，钝头体的气动加热明显降低。 20世纪50-60年代，最让人困惑的空气动力学问题之一。Why？同一个问题，出现了两个明显不同的区域，Ma1和1，此时偏微分方程的性质不同（抛物线和椭圆型），如何发展一致有效的方法呢？后来，布鲁克林工业大学的Moretti和Abbett发展了时间推进方法，解决了此难题。 * 第四章 偏微分方程的性质Behavior of Partial Differential Equations Slide * 超音速钝体绕流问题的解决 Slide * 方程的精确解： 含义： 以常速度c向右传播。 波形，振幅保持不变 * （常用）特例：常系数线性单波方程 偏微方程的分类及特征 基本概念：椭圆型、双曲型、抛物型方程 1. 一阶偏微分方程 初值： u x t=0 u x t=t0 t=0时刻与t=t0时刻物理量的分布 t x t=t1 t=t2 t=t3 x-ct=const 重要概念： 特征线 自变量空间的一条曲线，该曲线上物理量的方程可简化 A B c0 扰动波向右传播： 左端(A)需要给定边界条件； 右端(B)只能被动接受，无法给定边界条件 （即使给定，对计算域也无任何影响, 且造成B端的非适定性）。 c0 扰动波向左传播： 右端(B)需要给定边界条件； 左端(A)无需给定 线性单波方程的边界条件: 对于初值问题，如果微分方程解的定解域中存在、唯一、且连续依赖于初始值，则称数学问题的提法是适定的。 * 有限空间 重要基本概念，需掌握 初值： 问题： 如何给定边界条件？ （一般形式）一阶线性偏微方程 采用特征线法，可转化为常微分方程 考虑曲线G: 显然， 沿着该曲线G有： 如果该曲线G满足： 则有： 偏微方程在特征线上变成了常微分方程 特征线 特征相容关系 （特征线上物理量的简化方程） * x y 特征线简化了方程，在空气动力学领域应用广泛 Copyright by Li Xinliang * 演示： 如何利用特征线计算物理量 特征线 特征相容关系 计算域 x y 步骤： 1）设定积分步长 （根据精度需求设定，例如0.1） 2） 在边界上选取初始点 ，由边界条件确定该点的物理量值 3） 根据特征线及特征相容关系数值积分，求出特征线下一个点的坐标 和函数值 。递推下去，计算出整条特征线的（离散）坐标及物理量的（离散）值。 4）在边界上选取新的点，重复步骤3），计算出整个计算域物理量的分布 特征线法是空气动力学重要的计算方法。早期（计算机出现之前），是主要的CFD手工计算方法之一。 2. 一阶常系数偏微方程组 如果矩阵A 可以被对角化： 令： 有 即： m个方程完全解耦， 可独立求解 有m 条特征线： m个特征相容关系式： 如果矩阵A能够（相似变换）对角化，则原方程是双曲型的 * 如果矩阵A 具有m个实特征值， 这些特征值共具有m个线性无关的特征向量， 则称为双曲型方程 一阶拟线性偏微分方程组和m条特征线上的m个特征相容关系（常微分方程）等价。 如果A的特征值为m重根，而且对应的独立特征向量数小于m，则称为抛物型方程。 如果其A的特征值均为复数，则称为椭圆型方程 组合情况： 双曲-椭圆型 双曲-抛物型 * 3. 高阶偏微方程—— 可转化为一阶方程组 原方程化为一阶方程组： 转化为一阶偏微方程组 矩阵 特征方程（3）有两个互异