【用于高效、高精度热力学计算的计算包括（ H2、paraH2、orthoH2、N2、O2、Ar、H2O 和 CO2）（Matlab代码实现）

 💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现及文章

---

💥1 概述

冲击管被用于研究爆轰波，或者短暂地获得高压和高温，以便研究燃烧或其他反应。Glass 等人（Glass I.I. 1953）描述了一个简化的一维模型，忽略了气体粘度，适用于完美气体（即具有恒定热容的理想气体）的情况。在下文中，该模型被推广至非理想气体。

最初，冲击管由两个部分组成：高压（或驱动）部分和低压（或被动）部分。一个隔膜将这两个部分分隔开。在时间零时，隔膜破裂，气体从驱动侧流向被动侧。图1显示了波的轨迹和几毫秒后的一个样本压力曲线。

这些工具可用于进行高效、高精度的热力学计算。可用的种类包括 H2、paraH2、orthoH2、N2、O2、Ar、H2O 和 CO2。

左侧的稀疏波以局部声速向驱动气体传播。图中称之为“膨胀扇”。与此同时，一道冲击波向右侧的被动气体传播。由于我们忽略了粘度和因此摩擦，所以在稀疏波的右端和冲击波之间，压力和流速是恒定的。我们忽略了驱动气体和被动气体之间的任何混合，因此这个接触面会随着气体流速移动。在图1中，有四个区域：
1) 未受干扰的被动气体
2) 冲击波后的被动气体
3) 完全膨胀的驱动气体
4) 未受干扰的驱动气体
如果管道右侧有一个刚性的端壁，冲击波将撞击这面墙并反射，定义另一个区域：
5) 反射波后（右侧）的气体。
膨胀波将撞击驱动部分的末端，形成一个反射波向右传播。随着时间的增加，一系列相互作用的波模式会出现，变得越来越复杂。这些都没有包含在模型中。

详细文章见第4部分。

📚2 运行结果

部分代码：

% Tank parameters:
  tank.wallmass = 500;    % kg
  tank.cp_wall = 100;     % J/kg/K 
  tank.V = 5;             % m3
  tank.Area = 2;          % m2
  tank.heatcoeff  = 0;    % W/m2/K  (Adiabatic tank)
  
  Ta = 273.15+15;   % ambient temoerature
  T0 = Ta;
  p0 = 20e5;
  T_up = Ta;
  p_up = 350e5;
  w = 0.2;  % Filling rate (kg/s)
  p_final = 300e5;
  [t,T,N,p] = fill_gas_tank(th,tank,T0,p0,T_up,p_up,Ta,w,p_final);

  figure;
  subplot(311);
  plot(t,T-274.15)
  ylabel \circC
  title('Filling of H_2 tank with constant rate')
  subplot(312);
  plot(t,p*1e-5);
  ylabel bar
  subplot(313);
  plot(t,N*th.Mw);
  ylabel kg  
  xlabel seconds  
  
  fprintf('Final temperature %6.1f\n',T(end)-273.15)
  fprintf('Temperature increase %6.1f\n',T(end)-T(1))

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

🌈4 Matlab代码实现及文章