Euler方程组到UTSD方程的近似.pdfVIP

科学时代 Euler方程组到UTSD方程硇近似 王 丽 ／上海电机学院 [摘 要]利用流体速度在X方向和y方向的分量阶数的不同，给出了Euler方程组近似为UTSD方程的详细推导。 [关键词 ]Euler方程组 UTSD方程 在研究激波反射的问题中，人们曾采用了各种不同的模型方 它由入射激波i，反射激波 ，，马赫激波m和稀疏波 构成， 程组来描述流体运动，其中有完全欧拉方程组，位势流方程，非 反射点为 。在本节中我们仅在三叉点 的邻域内考虑如何将 线性波动方程 ，非定常跨音速小扰动方程 (UnsteadyTransonicSmall Euler方程组近似为USTD方程组。因为在图一中，流体是平行于 DisturbanceEqu~ion，简称UTSDE)等 [1112131141o虽然Euler方程组最 坡面的，故我们旋转坐标系使新坐标系(量，)的横坐标平行于坡面， 完整地描述了无粘流体的运动，但形式也是最复杂的。由于UTSD 即 方程组的非线性程度较低，又抓住流场变化中不同方向流动参 fli1一fCOS8sin]f] 数变化速率不同的特点，故我们可以将Euler方程组加以简化。 歹j—l—sin0c0sjlJ (2) 当vonNeumannparadox被提出以后，也有人试图在UTSD方程组 假设流体速度在 方向和．方向的分量分别是s阶和s％阶的， 的框架下来解释这一现象。下面我们将给出Euler方程组近似为 且设反射点 的坐标为((sec)，O)，其中c0是入射激波前的音速， UTSD方程组的详细推导。 =t。下面我们求具有如下形式的Euler方程组 (1)的近似解 我们考虑如下的欧拉方程组： p0 p1 f户 +( )+(y)Y：0， 0 0 ：：lc量，)，，+％l： p t p()，+：D()+(+p)嚣=0， ㈩、 ． p。 lp’J p 其中P。和 分别是入射激波iS的密度和压强，且(，Y，)由 (pe)，+(( +p))+(( +p)V)=0， 下式定义 其中 是密度，“和v分别是x和Y方向的速度分量，P是压强， = 专 ㈩ 8是内能，且满足理想气体的状态方程p=(一1)pe(一≠)， 将 (2)一 (4)代入 Euler方程组 (1)，得 其中y是气体的绝热指数。 A。+A1￡ +A2e+O(￡)=0， ‘ 假设一束激波i冲向一倾斜角为0的斜坡，发生如图一所示的 马赫反射结构。 CBoo++Bc1-s￡ ++C：2ze￡++O(s￡)：0。，， L(’) 1／ ／ D0+Dl￡ +D2￡+O(8 )=0． 因 为 ￡可 以是 任 意 小 的正 数，故 (5) 中 的系 数