- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
特征线理论基础 气体力学、空气动力学及计算流体力学 主干与重点基础内容 重要性 数学物理方法 Method of Mathematical Physics 运输方程特性 考虑运输方程 对于给定的初始条件 ,其任一时刻的解为 解在信息传播线上的导数 在信息传播线上取线段( ),看是否根据该线段两端的解的微分 确定解的导数。如果能,那么信息传播线附近区域的解可以通过泰勒展开求得。导数满足的方程为 以上方程如果有解,那么导数 便能由信息传播线上的解唯一确定。 导数表达式 直接求解得 信息传播线上的导数 以上代数方程组有解的充要条件是 而在信息传播线上 ,因此 。从而代数方程组无解,即导数不存在,亦即无法通过泰勒展开由一条独立的信息传播线上的解获得附近的解。 导数表达式 直接求解得 结论 信息沿信息线传播 由一条独立的信息线无法确定解的(尤其是法向的)导数,从而不能通过泰勒展开由信息线上的解确定附近区域的解 定义:信息传播线称为特征线,即特征线是流场中任一点上信息沿之传播的曲线。特征线的斜率称为特征值(运输方程和伯格斯方程的特征值分别为a和u) 特征线周边区域的解不能由特征线上的解求出 知道特征线上一点(如 )的解,在特征线上其它点 上的解便可以确定(例如,对于运输方程 )。 系统方程的情形 设系统方程为 在任一点,假设特征线的斜率为 ,即特征线方程 为 特征线上的微分关系式 构造方程组 特征值 类似前面分析,特征线斜率(特征值)满足方程 对于含有m个分量的方程组,有m个特征值,即m条特征线。 相容关系式 沿特征线的导数定义为 对每个 ,左特征向量由下式确定 注意,左特征向量各分量对应书中 相容关系式 将左特征向量乘以 得 以上即为相容关系式,也称为特征线形式的方程,对于每个特征值,都有一个(这样的标量)相容关系式,属于常微分方程。因此,沿特征线方向的变化可以通过求常微分方程得出 相容关系式的作用 沿每条特征线 ,都有一个相容关系式。每个相容关系式定义了一个组合量 沿特征线的变化。往往该组合量沿特征线为常数(称为不变量,对于气体动力学方程,称为黎曼不变量)。 思考题 对每个特征值,由 知 讨论 与相容关系式有何联系?书上说,二者是等价的,是否真的等价? 等价性? 对每个 ,只有一个相容关系式 一般系统方程 对于m阶系统方程,总共只有m个相容关系式,但有 个分子行列式(因为有m个特征值,而对每个特征值,有m个分量),由m个分母行列式为0使得独立的分子行列式为0的个数为 。 因此,如果m2, 即 ,二者就不应该等价。但我们能证明,针对有三个方程的气体动力学方程是等价的。 特征线法的作用:工程问题 在非特征线曲线L上给定初始条件,求曲线外任一点P上的值。对于有m个分量的系统,必有过P点的m条特征线与曲线L上的m个点 相连。由于沿每条特征线都有一个相容关系式将P点和曲线L上的一个不变量相关联,所以可以确定m个这样的不变量 ,从而唯一确定P点所有的流动参数。 高维问题特征线理论 在一维情况下,特征线为为曲线;在二维情况下应该为曲面即特征面;三维情况下为特征体。 有时考虑一些特征面或特征体上的特殊曲线,作为特征线。 有时将问题投影到各方向(如x,y,或z),把被投影方向作为一维问题考虑。 二维问题特征面 基本方程 在特征面上,由相距(dx,dy)的两点的解所定义的微分 dw定义偏导数,因此偏导数所满足的方程为 二维问题特征面续 考虑特征面上的相互正交的曲线,满足方程 于是有 二维问题特征面续 下面看在什么情况下导数不存在
文档评论(0)