这是在无粘性和不可压缩情况下作用于矩形横截面的振荡圆柱体的流体动力流体力的解析解的简单实现（即拉普拉斯方程的解）。 解决方案取自以下附录 B： Brumley、Douglas R.、Michelle Willcox 和 John E. Sader。 “浸入流体中的矩形截面圆柱体的振荡。” 流体物理学 22.5 (2010): 052001。 在纵横比大于 1000 时，结果由渐近解近似。

不可压计算流体苦学

用matlab的特征线法求解瞬变流，针对可压缩流体在管道中的流动。

欧拉 可压缩流体的欧拉方程

请注意，虽然每个子函数都可以独立运行，但我建议使用顶级函数 compressible.m 以避免格式混乱。 该函数求解与等熵可压缩流相关的关系、法向激波关系、具有热量增加和摩擦的等熵流、Prandtl-Meyer 函数和马赫角的解以及斜激波的 Theta-Beta-Mach 关系。 它对任何输入值和任何 gamma 值都这样做，其中 gamma 是比热的比率。 还要注意，这个函数可以处理输入向量，如果可以的话，它会以与输入相同的形状返回它们。 该功能可以通过三种方式使用： -如果未指定输入或输出，它将运行GUI，提示用户确定哪种表和输入类型以及将使用的伽玛。 然后 GUI 将结果显示在集成到图中的表中。 - 如果包含输入但不包含输出，则该函数将在工作区中打印结果。 如果这些函数用于查找参考值或家庭作业分配，则上述方法很有用。 如果要在函数内使用计算，则第三个选项是： - 如果包含输入和输出

函数 flow2d 为由复势给出的不可压缩流体的二维势流生成流线、速度场和动态压力场的等高线图。 假设流动是静止的并且流体没有自由表面。 重力的影响被忽略。 例子： flow2d('z + 1/z') 产生流线的等高线图，速度场和压力场对应于势 f(z) = z + 1/z。

工程流体力学第七章理想不可压缩流体的有旋流动