《An introduction to COMPUTATIONAL FLUID DYNAMICS, The Finite Volume Method》
2.7 The role of characteristics in hyperbolic equations
考虑上一节提到的双曲线型方程,令,。波动方程转化为
可以求得。其中和为任意函数,为问题的简单波动函数解,波以速度c和-c传播且波型和幅值不变。
考虑无限长细绳,初始条件为,代入求解可以得到:
从式中可以看出在点的解仅仅受到初始条件中范围内的值的影响,通过下图来说明。
在图中画出特征线和两条交于点。该点上的值受到依赖域上的值影响。可以通过波速确定,信息沿两簇平行线随时间增大传播,因此越过交于这点的两条特征线信息就不会传播到该点,只有依赖域内的值有可能影响到该点值。同理也可以推得,该点的信息在之后的时间里也只会传播到由这两条特征线所包围的影响域内的点。
2.8 Classification method of for simple PDEs
考虑一般二维二阶偏微分方程:
其特征方程(?)为
通过判断特征方程实根的个数判断偏微分方程类型。即通过判别式如下表判断。
其中可以为常数或函数。
N个独立变量的二阶偏微分方程,先将他们重写为如下标准形式:
其中,即将系数分为相等的两部分,凑成矩阵形式。求矩阵的特征值,根据特征值的情况判断偏微分方程所属类型。
前面所提到的拉普拉斯方程,波动方程和扩散方程都可以通过此两种方法判定所属类型。
2.9 Classification of fluid flow equations
表中可得,稳态的NS方程为椭圆型,瞬态的为抛物线型。无粘流体和粘性流体相比,由于缺少粘性的高阶项而大不一样,其分类和马赫数相关。稳态问题中,由于压力是以声速在流动中传播,如果马赫数小于1即流动速度小于声速,压力便可以传播到所有区域,方程特征即为椭圆型方程。而如果马赫数大于1,流动速度大于声速,压力只能传播到部分区域,方程特征即变为双曲线型特征。边界层流动仅含一项二阶扩散项,所以为抛物线型方程。
在稳态问题中出现了信息不能传播到整个区域的双曲线型和抛物线型方程,其求解是类似时间前进的方式,联想到COOLFluiD接触过的用瞬态逼近稳态的方法设置。