冲击波的计算机实践联系.docxVIP

冲击波物理实践报告实践项目一要求：基于计算和曲线对空气和氦气中的冲击波规律进行分析。自行设定初态条件（p0、v0、T0、u0等），编程或利用软件画出以下曲线：①p-v平面上的冲击压缩线；②?/?0随p/p0的变化曲线；③u随p/p0的变化曲线；④D随p/p0的变化曲线；⑤T/T0随p/p0的变化曲线。基于具体计算结果指出空气和氦气中强冲击波的具体条件。1、空气假定初始条件，（1）p-v平面上的冲击压缩线公式：当（2）随p/p0的变化曲线公式：分析图二，随着p/p0的增大，上升速度渐渐放缓，趋近于某一个值。计算p/p0+1相对于p/p0对应的的百分误差，即，得到下列图像：可以看出，当趋近于1000时，的相对变化值接近于0.计算时，，即空气强冲击波极限压缩比为6。（3）u随p/p0的变化曲线公式：由图像分析波后速度u随p/p0的变化趋势，正相关。扩大p/p0的范围，得到如下图像：求极限，，因此没有极限波后速度。（4）D随p/p0的变化曲线公式：D与p/p0成正相关，扩大p/p0的范围，得到如下图像：求极限，，因此理论上冲击波速度没有极限。（5）T/T0随p/p0的变化曲线公式：观察图五发现，p/p0 在一定的范围内，T/T0 随着p/p0的变化呈线性变化。拟合图像，得到黑色虚线拟合方程，进一步印证了线性趋势。增大p/p0到10000，得到下面图像，发现在很大范围内，T/T0 随着p/p0的变化近似于线性。（6）强冲击波的具体条件之前我们得到了压缩比的相对变化曲线：假定的相对变化率达到0.01%时，达到极强冲击波条件，计算得到此时的波后压强比约为273.96。p/p0=273.96即为强冲击波的具体条件。2、氦气假定初始条件：（1）p-v平面上的冲击压缩线公式：当（2）随p/p0的变化曲线公式：分析图二，随着p/p0的增大，上升速度渐渐放缓，趋近于某一个值。计算的百分误差，得到下列图像：可以看出，当趋近于1000时，的相对变化值接近于0.计算时，，即空气强冲击波极限压缩比为6。（3）u随p/p0的变化曲线公式：由图像分析波后速度u随p/p0的变化趋势，正相关。扩大p/p0的范围，得到如下图像：求极限，，因此没有极限波后速度。（4）D随p/p0的变化曲线公式：D与p/p0成正相关，扩大p/p0的范围，得到如下图像：求极限，，因此理论上冲击波速度没有极限。（5）T/T0随p/p0的变化曲线公式：观察图五发现，p/p0 在一定的范围内，T/T0 随着p/p0的变化呈线性变化。拟合图像，得到黑色虚线拟合方程，进一步印证了线性趋势。增大p/p0到1000000，得到下面图像，发现在很大范围内，T/T0 随着p/p0的变化近似于线性。（6）强冲击波的具体条件之前我们得到了压缩比的相对变化曲线：假定的相对变化率达到0.01%时，达到极强冲击波条件，计算得到此时的波后压强比约为191.04。p/p0=191.04即为强冲击波的具体条件。实践项目二利用计算机求解例题3.7和习题3.30。求解习题3.30：第1步：根据已知条件求出初态声速： p0=1 atm，u0=0，T0=300 K，第2步：根据冲击波关系求出1区和2区的状态首先利用式(3.3.21)求出温度：利用p1=2 atm，求得T1=369 K；利用p2=4 atm，求得T2=480 K然后利用式(3.3.19)，求声速：结果为c1=385 m/s，c2=439 m/s。接下来利用式(3.8.31)求波后质点速度由此求得u1=182 m/s，u2=394m/s。第3步：两个右行冲击波作用反射一个左行稀疏波、投射一个右行冲击波。对于右行冲击波有：对于左行稀疏波有：p3=p4u3=u4求解例题3.7：第1步：根据已知条件求出初态声速： p0=1 atm，u0=0，T0=300 K，第2步：根据冲击波关系求出1区和2区的状态首先利用式(3.3.21)求出温度：利用p1=2 atm，求得T1=369 K；利用p2=4 atm，求得T2=480 K然后利用式(3.3.19)，求声速：结果为c1=385 m/s，c2=439 m/s。接下来利用式(3.8.31)求波后质点速度由此求得u1=182 m/s，u2=-394m/s。第3步：根据冲击波的两个公式(3.8.32)和(3.8.33)画p-u曲线（计算压强随波后质点速度之间的关系），将数据点保存在一个数据文件中，然后观察两组数据，当速度与压强相等（近似）时便是解。这是阻抗匹配法。p3=p4=7.04 atm；u3=u4=-208 m/sT3=558K， T4=566Kc3=473 m/s ； c4=477 m/s结果：求交点得到（-208.695，7.0371），即结果为速度-208.695m/s,