首页/文章/ 详情

【经典教材翻译】17-激波的相互作用和分离

1月前浏览1050

1

激波相互作用和分离激波

许多超音速应用涉及斜激波。这些包括从实体墙上反射的激波、激波的交叉、高速车辆产生的激波等。


2

激波从墙面反射

第一个应用是斜激波以角度β冲击墙面,如图37所示。这个斜激波可能是由于在马赫数M1的超音速流动中放置了一个楔形物而产生的。流动被偏转了一个角度θ。底部的墙面将流动重新导向水平方向。这产生了一个反射激波,它与墙面成β'的角度。虽然流动偏转角度相等但方向相反,但反射波的角度不一定等于入射波的角度。它们只有在波是马赫波并且没有偏转流动时才会相同。每波前的马赫数不同,因此对于类似的转向角度,斜波角度也不同。入射激波是在马赫数为M1的流动中产生的,而反射激波是在马赫数为M2的地方产生的。

图37:激波从墙面反射

可以通过依次将前面的斜波结果应用于两个波,轻松计算出图示的流动特性和压力分布。


3

两个激波的交叉

当两个激波以图38所示的角度相互交叉时,就会发生两个激波的交叉。如果激波的强度相等,相互作用就像每个激波都从流动的中心线反射一样,实际上这是一条流线。流线总是可以被当作墙面来对待,以便计算过程与单个激波反射相同。

图38:两个对称激波的交叉

当激波的强度不相等时,流场就会失去对称性。在相互作用的下游出现了一个新的特征,如图39所示。这被称为滑流流Slip Stream。这将流动分为两部分 - (1) 顶部激波处理过的流动和 (2) 底部激波处理过的流动。滑流流是这样的,即压力P3和流动角度δ在其中连续。密度、温度和其他属性将是不同的。要解决这种情况,需要一个迭代过程,以便滑流两侧的压力相等。

图39:两个非对称激波的交叉


4

强解 - 分离激波

另一个情况是流动冲击一个物体,其前端是这样的,使得θ > θmax(其中θmax是前一节中计算的斜激波允许的最大流动偏转角度)。在这种情况下,激波不附着在前端,而是远离它 - 即,它是分离的Detached,如图40所示。激波不再是一条直线,而是弯曲的,其形状和强度取决于M1和物体的几何形状。在中心线上,激波是一个正常激波在a处。远离中心线,激波减弱,并在d处接近斜波。从a到d存在一系列强激波解。

图40:楔形物前的分离激波

这种激波波解的组合使得激波下游的流场有些复杂。正常激波下游的流动是亚音速的。因此,在激波下游的中心线附近有一个亚音速流动区域。这个区域的范围取决于物体的几何形状和自由流马赫数。远离中心线,激波最终达到一个带有后方超音速流动的弱斜波解。一个声速线将超音速流动与亚音速区域分开。

如果物体是钝的,如图41所示,那么在所有马赫数下激波都是分离的。

图41:钝体前的分离激波

物体和激波之间的距离称为激波立距Shock Stand Off Distance,它随着马赫数的减小而减小。可以使用CFD技术或实验测量来计算这个距离。


5

马赫反射

在许多应用中,激波在实体表面上的简单反射并不发生。当下游的马赫数M2是这样的,即不可能有弱斜激波解时,就会发生这种情况。在这些情况下,激波的反射不是在实体墙上发生的,而是在一定距离之外。如图42所示。现在流动中有一个三重点,然后是滑流流。这种现象称为马赫反射Mach Reflection。同样需要通过迭代计算属性,直到三重点背后的压力在两个下游区域相等。

图42:马赫反射




来源:CFD饭圈
通用
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-09-08
最近编辑:1月前
CFD饭圈
硕士 分享CFD文章,感谢关注
获赞 22粉丝 23文章 378课程 0
点赞
收藏
作者推荐

Fluent软件中的Y+值

无量纲变量Y-Plus代表从壁面到第一个网格单元中心的距离,同时,Y-Plus是一种无量纲距离。换句话说,Y-Plus使得单元中心与壁面之间的垂直空间无量纲化,这个值用于评估网格单元大小(高度)在壁面上的比例。另一方面,Y-Plus表示单元内湍流和层流效应的比率。如果Y-Plus值较低,该单元的流动是层流的。如果Y-Plus值较大,该单元中的流动将是湍流的。因此,如果Y-Plus太小,那么单元中的流动是层流慢速的,壁面函数将不可用。1、Fluent中Y-Plus范围在使用AnsysFluent中的壁面函数时,Y-Plus应该在以下范围内指定:根据上表,如果Y-Plus的值等于300,我们将使用K-Epsilon(StandardWallFunctions,Non-equilibriumWallFunctions,ScalableWallFunctions)湍流模型。如果我们有Y-Plus值为1,我们将使用K-Epsilon(EnhancedWallTreatment),K-OmegaSST,和Spalart-Allmaras。2、举个实际例子为了更好地理解上述内容,我们尝试以在ANSYSFluent软件上完成的实际例子来解释naca0018翼型的情况。在零攻角下模拟围绕naca0018翼型的湍流。我们将使用K-Epsilon(StandardWallFunctions,Non-equilibriumWallFunctions,ScalableWallFunctions),K-Epsilon(EnhancedWallTreatment),K-OmegaSST,Spalart-Allmaras湍流模型来解决这个问题。2.1壁面定律粘度和壁面剪切应力是靠近壁面的重要特性。我们根据这些变量定义了粘性尺度,即近壁区域合适的速度和长度尺度。使用方程1到4,我们可以指定单元到壁面的距离或y的值。当前问题中使用的初始和边界条件如下:2.2Y-Plus等于300时如何计算Y-Plus等于300时,与壁面相邻的单元的高度将这个值放入上述关系中,y的值为0.006889。使用相同的网格、材料参数、求解器和边界条件运行模拟,以说明基于单元的方法的影响。Y-Plus的区域平均最大值和最低值对应于上述数据。使用的网格,使得y-plus=300,可能如下所示。这个案例中,在fluent中打印了Y-Plus的值,显示壁面附近的第一层高度已经被正确分布:从上面的图中可以看出,Y-Plus的值在30到300之间。2.2Y-Plus等于1时同样的计算步骤,Y-Plus为1时,y的值为0.0000229,这表明第一层的高度值变得小得多。为了使网格具有高质量,我们必须增加边界层的数量,或者通过使用纵横比,我们可以测量最后一个边界层与其邻近网格。当Y-Plus1时产生的网格如下;根据图示,与之前的状态相比,网格单元的数量增加了。这个案例中,在fluent中打印了Y-Plus的值,显示壁面附近的第一层高度已经被正确分布:如上图所示,Y-Plus的值低于1。2.3结论下表显示了1和300之间Y-Plus的差异:来源:CFD饭圈

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈