第2章冲击波导论课件.ppt
《第2章冲击波导论课件.ppt》由会员分享,可在线阅读,更多相关《第2章冲击波导论课件.ppt(114页珍藏版)》请在三一文库上搜索。
1、2.1 波的基本概念 1、波(Wave)波通常可以分为两大类:一类是电磁波,另一类是机械力学波。当介质(Medium)受到外界作用(如振动、冲击等)时,介质的局部状态参量就会发生变化,这就是扰动(Disturbance)。波就是扰动的传播。2.1 波的基本概念 n如果活塞突然向右移动,便有波向右传播。n在扰动传播过程中,扰动介质与未扰动介质之间存在一个界面,这个界面就叫波阵面波阵面(Wave front)。n扰动在介质中的传播速度叫做波速波速(Wave velocity)。(要与介质的质点速度区分)2.1 波的基本概念 n如果扰动前后介质的状态参数变化量与原来的参数量相比是很微小的,则称这种扰
2、动为弱弱扰扰动动(Weak disturbance)或小小扰扰动动。弱扰动的特点是各种参数的变化量是微小的、逐渐的和连续的。如果扰动前后介质的状态参数发生突跃变化,则称这种扰动为强扰动强扰动(Strong disturbance)。2.1 波的基本概念 2、声波声波(sound wave)n声波是一种弱扰动波。弱扰动在介质中的传播速度就叫声声速速。它是气体动力学中一个非常重要的参数。n下面以活塞在直管中移动 所引起的气体扰动的传播 来建立声速c与其它参数 的关系式。如图所示。(1)式中x为t1时刻扰动传播的距离,x=ct1 x1为时刻活塞运动的距离,x1=ut1 A0为活塞的截面积。代入(1)
3、式可得:消去t1后可得:(2)2.1 波的基本概念 n 质量守恒质量守恒(Conservation of Mass):2.1 波的基本概念 n动量守恒动量守恒(Conservation of Momentum):气体受到扰动后的动量等于作用在其上面的冲量。化简后得:(3)(2)式代入(3)式得:(4)由(2)式可得:(5)2.1 波的基本概念 把(5)式代入(4)式得:(6)n由于声波为弱扰动波,波阵面过后介质状态变化为一微小量,故有 ,因此,(6)式变为:(7)n看作等熵过程:(8)2.1 波的基本概念 对于理想(多方)气体,其等熵方程为:(9)则 (10)所以理想气体的声速为:(11)又由
4、 可得:(12)2.1 波的基本概念 n对于地表面上的空气,可近似地视为理想气体,将 ,代入上式可得:(13)将 代入(13)式可得 340m/s。2.1 波的基本概念 n需要指出的是,只有对于小扰动,才成立,扰动才以声速传播。对于 的扰动,其传播速度大于声速,扰动越强,传播速度将越高。2.1 波的基本概念 3、压缩波和稀疏波压缩波和稀疏波 压缩波压缩波(Compression Wave):扰动传过后,介质的压力、密度、温度等状态参数增加的波称为压缩波。其特点是波传播的方向与介质质点运动方向相同。稀疏波稀疏波(Rarefaction Wave):扰动传过后,介质的压力、密度、温度等状态参数下降
5、的波称为稀疏波稀疏波,其特点是波传播的方向与介质质点运动方向相反。2.1 波的基本概念 n在一个连续的,缓慢的压缩过程中,每一小步的压缩都是一种等熵变化,但由于每经一步压缩后气体的温度都要上升,气体的声速必将上升,这样下一步的压缩波的波速逐渐增加,一旦集中起来,状态参数的变化将不再连续,就会发生突跃,弱扰动变成强扰动强扰动。2.1 波的基本概念 n由于稀疏波的膨胀飞散是按顺序连续进行的,所以稀疏波传播中介质的状态变化是连续的,如图24中的压力变化。图24稀疏波现象 2.1 波的基本概念n在稀疏波扰动过的区域中,任意两相邻端面的参数都只差一个无穷小量,因此稀疏波的传播过程属于等熵过程,它的波速等
6、于介质当地的声速或音速(Local sound speed)。2.2 气体的平面一维流动 n 所谓一维流动,是指在某一空间坐标x等于常数的平面上流体参数都是均匀分布的,并且在给定坐标x处的流体参数都只与时间t变化相关的流动。平面一维流动规律的求解目标是确定一维流场中介质参数随时间t和空间x的变化规律,如p=p(x,t),T=T(x,t),u=u(x,t),=(x,t)。n一维流动又可分为一维定常流动和一维不定常流动。2.2.1气体一维流动的基本方程组 n气体在平面一维流动下,满足质量守恒、动量守恒和能量守恒,其对应的方程分别叫质量方程(连续方程)、动量方程(欧拉方程)、能量方程。1、连续方程(
7、质量方程)(1)该式为一维不定常流动的连续方程。2、欧拉方程(动量方程)(2)2.2.1气体一维流动的基本方程组 3 3、能量方程能量方程n在不考虑气体的粘性和热传导的情况下,气体的流动是等熵的。(3)4 4、状态方程状态方程n由于S可表示为p和 的函数,故等熵流动条件可表示为:n对于理想气体,其等熵方程为:(4)这样,便可由连续方程、欧拉方程、能量方程和状态方程求解气体一维等熵流动的四个未知量 。2.2.2 以u、c为求解参量的方程组 n为使前面建立起来的气体一维等熵流动的方程组的物理意义更容易理解,将它们稍加变换。引入声速c代替p和 。由声速公式及等熵方程可得:(5)将(5)式两边微分并同
8、时除以 ,得 (6)2.2.2 以u、c为求解参量的方程组 由(5)式知,(7)把(6)式代入(7)式,可得:(8)2.2.2 以u、c为求解参量的方程组 将(6)式代入连续方程(1)式,可得 (9)将(8)式代入欧拉方程(2)式,可得 (10)2.2.2 以u、c为求解参量的方程组 将(9)、(10)两式相加和相减,整理可得 (11)这个方程组即是以u、c为变量描述气体一维等熵不定常流动规律的方程组。确定气体一维等熵流动过程中气体各参数时的时间、空间变化规律,归结为解此偏微分方程组。2.2.2 以u、c为求解参量的方程组n小扰动波在静止介质中是以音速进行传播的,在一维情况下,静止气体中小扰动
9、波的传播速度为c。在流动介质中,小扰动波的传播速度为介质流动速度u与当地音速c的叠加,即 。顺介质流动方向传播的扰动取正号,逆介质流动方向传播的扰动取负号。2.2.2 以u、c为求解参量的方程组n在 条件下,(11)式可表示为 对t的全导数形式,并且该导数为零,即 (12)即 (13)2.2.2 以u、c为求解参量的方程组 由此可以看出,方程(11)在 条件下描述的是两个量的推进规律:n即由 所确定的状态(或扰动)以速度 顺气体流动方向(即x轴的正方向)传播;n而由 所确定的状态(或扰动)以速度 逆气体流动方向传播。2.2.3方程组的特征线及一般解2.2.3方程组的特征线及一般解 dx/dt=
10、u+c和dx/dt=u-c分别代表一维等熵流动介质中扰动沿x轴的正向和反向传播的速度,我们称它们为(11)式的特特征征或特特征征方方程程。它们的积分各自代表xt平面上的一簇曲线,叫做特特征征线线。其中在xt平面上由dx/dt=u+c所确定的特征线称为第第一一簇簇特特征征线线,用C表示;而由dx/dt=u-c所确定的特征线称为第二簇特征线第二簇特征线,用C表示。2.2.3方程组的特征线及一般解 n这两簇特征线分别描述的是物理状态量 ,即扰动波以速度 沿x轴的正向或负向传播的轨迹。n因此,对于一维等熵不定常流动方程组(11)式,有沿着C特征线 (14)2.2.3方程组的特征线及一般解 沿着C特征线
11、 (15)式中,I+,I-称为黎曼(黎曼(RiemannRiemann)不变量不变量。它们在u,c平面上可用两簇相互平行的直线来描述,称为方程组(11)在速度平面上的特征线。它们在沿着各自的特征线(C和C)传播时保持不变。如图25所示。2.2.3方程组的特征线及一般解 图25 特征线2.2.3方程组的特征线及一般解 n方程(14)和(15)为方程组(11)的一般解。n在k3的最普通的情况下,由于(u+c)和(u-c)都是x和t的函数,即右传波的传播速度受反方向波的影响,因此(14)式和(15)式无法得到精确的解析解。一般采用数值积分法或特征线法近似求解。2.2.3方程组的特征线及一般解 n在x
12、t平面上,假设曲线AB上的各点状态参数已知。C和C分别表示AB线上各点发出的不同簇的特征线,求解流场D内各点的状态参数。图26 特征线法解流场参数 2.2.3方程组的特征线及一般解【解】近似认为 ,并近似把x-t平面上特征线的一小段视为直线。在曲线AB上选取一系列的点M1,M2,Mi等。由于ui和ci已知,过Mi点作特征线C,其特征方程为:2.2.3方程组的特征线及一般解 而过Mi+1点作特征线C,其特征方程为:其中x和t为C和C相交点Mi空间坐标和时间,u和c为该相交点Mi的状态参量。同样可以求出AB上各点的参量,依次可求出任意位置点的状态参量。2.2.4 方程组的特殊解简单波流动2.2.
13、4 方程组的特殊解简单波流动 前面讨论的(14)和(15)式是方程组(11)式的通解,流场中可同时存在左传波和右传波。如果流场中只有一个方向传播的扰动波,即波未进入的区域介质处于静止状态或稳定流动状态,这种波就称为简简单单波波,其解称之为方程组的特殊解。简简单单波波:它的某一族特征线上的黎曼不变量是同一个常数,即该族各条特征线上的黎曼不变量彼此相等。2.2.4 方程组的特殊解简单波流动 当给定如下条件,即 (16)对(16)式分别对t和x求偏导,得 将这两式代入(11)式可得 (17)2.2.4 方程组的特殊解简单波流动 n该式表明,沿特征线dx/dt=u+c,有du/dt=0。即u常数。n由
14、16)式知,c亦为常数。因此dx/dt=u+c就可以积分了。因此 (18)2.2.4 方程组的特殊解简单波流动 同理,当 时,有 (19)式中,是u的任意函数,由边界条件确定。由(18)和(19)式即可确定简单波的向前波(右传波)和向后波(左传波)流动区内任一点的参数u和c。2.2.4 方程组的特殊解简单波流动 n为了阐明简单波的性质,我们来考察下面两种情况。(1)活塞向左加速运动,如图所示。图2-7 右传系数波 2.2.4 方程组的特殊解简单波流动 n当活塞向左加速拉动时,便形成一系列的简单稀疏波向右传播,并以当地声速传播,因此,活塞向左拉动时发出的第一道稀疏波是以静止气体当地的音速u0+
15、c0=c0的速度向右传播的,特征线如图2-7所示。n该特征线的右边为静止气体区域,故该区域内的特征线也都是平行的。活塞向左加速拉动而发出的各个后续右传稀疏波,是在扰动过的气体中传播的,因此第n道波的传播速度(un+cn)总是比其前面的的(n-1)道波的传播速度(un1+cn1)要慢。因此后面的各道波的特征线C是发散的。2.2.4 方程组的特殊解简单波流动(2)活塞向右渐渐加速运动,如图2-8所示。图2-8 压缩波随t的变化 2.3 平面正冲击波2.3 平面正冲击波 n冲冲击击波波(Shock wave),又称激激波波,是一种强烈的压缩波,其波阵面通过的前后参数变化很大,它是一种状态突跃变化的传
16、播。冲冲击击波波阵阵面面(Shock front)实际上有一定的厚度,其厚度约为几个分子平均自由程,在这个厚度上各物理量发生迅速的、但却是连续的变化,这是由于物质具有粘性和热传导的原因。但在工程计算上可以不考虑粘性和热传导等耗散效应,而将冲击波视为一个没有厚度的间断面。因此,可以说冲击波阵面是一种强间断面。2.3.1 基本关系式2.3.1 基本关系式 n 设有一冲击波以恒定的速度向右传播,如图2-9所示。图2-9 平面正冲击波阵面 2.3.1 基本关系式n波的右边,尚未扰动的介质,参数为:。n波的左边,扰动的介质,参数为:。n为方便起见,把坐标系建立在波阵面上。则未扰动的介质以D-u0的速度向
17、左流入冲击波阵面,扰动的介质以D-u的速度从波阵面流出。1、质质量量守守恒恒(Conservation Conservation of of MassMass):单位时间内流入波阵面的质量等于流出的质量。即:(1a)将 ,上式变为:(1b)2.3.1 基本关系式 2、动动量量守守恒恒(Conservation of Momentum):单位时间内作用介质上的冲量等于其动量的改变。冲量:动量变化:因此 (2a)即 (2b)2.3.1 基本关系式(3)能能量量守守恒恒(Conservation of Energy):冲击波传播视为绝热过程,忽略介质的粘性和热传导效应等能量耗散。单位时间内从波阵面右
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 冲击波 导论 课件
