理论力学电子教案.pptx

《理论力学电子教案.pptx》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、二维气体流动1有一个平面无旋超音速流动，它被限制在夹角为15。的直固壁面中，沿圆弧线1-2-3-4气流是径向的，并且具有相同的马赫数M=1.348，试用逐点的特征线法确定流场中各点的M和压力p/pe。题目解：图1(a)为物理平面网格点，(b)为速度平面，(a)图上沿初值曲线选择均匀分布节点l，2，3和4。各点的M和为已知值，见表1。速度平面上相应的节点为l´，2´，3´和4´。在速度平面上，由经过点1´的第I族特征线和经过2’点的第Ⅱ族特征线的交点定出5´。2图1(a)为物理平面网格点3图1(b)为速度平面4根据1´-5´的平均角(-μ)和2´-5´的平均角(+μ)，见图2在物理平面

2、上绘制特征线段1-5和2-5，这样就定出了点5的位置，同样可以绘制出其他相应的点6，7，8，……。其具体做法如下：解题步骤表1平面节点的已知数据节点1234M1.3481.3481.3481.348º7.5º2.5º-2.5º-7.5º5图2流场节点的气流方向角6在节点1处，M1＝1.348，查普朗特-迈耶表得到ν(M1)＝ν1＝7.5º，1＝7.5º，μ＝47.9º，速度平面上点1´的坐标号码为解题步骤在节点2处，M2＝1.348，ν(M2)＝ν2＝7.5º，2＝2.5º，速度平面上点2´的坐标号码为7在速度平面上的点l´的第I族特征线与点2´第II族特征线的交点5´的I＝392.

3、5，II=597.5，所以按ν(M5)查表得M5=1.435，μ＝44.2º。其他各点可按同样办法处理。解题步骤8在确定点8和8´时，注意到点8在物体的边界面上，即点8´处的速度方向是已知的，因此可由经过点5´的第II族特征线来确定点8´的位置，然后再以点5至点8的平均角+μ值，在物理平面上绘出相应的特征线，它与壁面的交点就是节点8，而后用类似办法定出其他壁面点。具体做法如下:因为已知8=7.5º，II＝597.5，故题目9按ν(M8)=ν8查表得M8=1.52。以此办法就可以在物理平面和速度平面上绘出整个特征线网。表2是各节点的计算结果数据，下划线数据为已知数据。解题步骤10表2各节

4、点计算结果节点号码I号码II(º)ν(M)(º)Mμ(º)(+μ)(º)(-μ)(º)p/p01392.56007.57.51.34847.955.4-40.40.3382395597.52.57.51.34847.950.4-45.40.3383397.5595-2.57.51.34847.945.4-50.40.3384400592.5-7.57.51.34847.940.4-55.40.3385392.5597.55.0101.43544.249.2-39.20.29963955950101.43544.244.2-44.20.2997397.5592.5-5.0101.435

5、39.239.2-49.20.2998390597.57.512.51.5248.648.6-33.60.26511已知喷管超音速段的进口端马赫数M=1，出口端马赫数ME=1。试用特征线法设计此缩放喷管超音速段长度为最短的型线(假定管内为二维流动，＝1.4)。解：如图1(a)为物理平面上喷管超音速段的气流流向及波系图，喷管中心线为一条流线，故可用一直固壁替代，图中虚线a为音速线，点0为凸曲壁扰动点，用4条入射膨胀波来代替整个膨胀波束，因此相应的有4条反射膨胀波，气流穿过8条膨胀波后为直匀的超音速气流，ME=2。题目12图1(a)物理平面13图1(b)速度平面14解题步骤当M=1时，ν(M

6、)=0º；ME=2时，ν(ME)=26.38º,总的普朗特-迈耶角变化为26.38º，则其平均增量为:其相应的方向角的变量为：Δ=±3.298°，其正负号分别用于特征线I和II上。亚音速a的上游与反射膨胀波10-10下游，各为均熵均焓，气流方向相同的流场通过全部膨胀波系后=0º，并必须与收敛段壁面轮廓相一致。由于扩张段集中在点0，由此可知壁面型线的最大倾角max为：15解题步骤点4、7、9、10为消波段壁面型线上的点，各点有相同的转折角3.298º。从区域a至e的膨胀过程，气流跨越4根第I特征线01，02，03和04，其第II特征线为常数。在起始线上，M=l，νa=a=max，和

7、e区域νe=e=max，所以第II特征线II为：16解题步骤从区域e至o的膨胀过程，气流经过4根第II特征线线4-4，7-7，9-9和10-10，故其第II特征线保持不变。对于e区，νe=13.19º，e=13.19º，对于o区，νe=26.38º，o=0º，所以第I特征线I为：表3为特征线法的计算结果，表中数字下面有划线的，为两个求解用的初始参数。17表3各区间计算结果节点号码I号码II(º)νa(M)(º)