调查的压力和速度测量系统不稳定现象也见图
1。压力传感器(Kulite xcq - 062 - 25 - a)的误差在±0.1%是安装充裕的裹尸布表面压缩机。压力探头的直径是1.7毫米。压力测量沿叶轮流道以下几点:叶轮入口(I.I.),叶轮上腹部(即时),和扩压器进口(扩散指数一个和D.I.b)。以上三个点的直径是0.757
D
2,0.758
D
2和1.024
D
2,分别。此外,impeller-discharge流测量的知识分布扩散圆周方向(D.I.1-D.I.4),使用由裂膜探针(DANTEC 55 r57)恒温的风速计(90 n10 DANTEC,简化)。裂膜探针的长度是8毫米胶片长度为1.25毫米。测量的位置
D
=
339年mm位于叶轮出口和扩散之间的前缘。四个计量点放置在0度,90度,180度,240度。impeller-discharge流量测量是遍历从中心到裹尸布边每扩散宽度的5%(20分)。不稳定压力和速度跟踪测量25.6 kHz的采样频率。时间分辨率是大约20个样本/叶片传递。
3所示。计算过程3.1。控制方程
目前使用内部进行了模拟计算流体动力学(CFD)代码,验证了各种涡轮机械流。代码解决了控制方程的连续性,三维可压缩n - s方程,能量方程,理想气体的状态方程。对流通量进行了通量差分分裂(FDS) [
18),这是扩展到使用单调的三阶逆风保守的法律(MUSCL)插值方案
19]。粘性通量测定的二阶中心差分使用高斯定理。矩阵免费高斯-赛德尔(生产厂)隐式算法用于时间集成(
20.]。DES方法,这是一种混合的方案涉及大涡模拟(LES)和跑
21),采用湍流模型。DES是基于剪切应力传输(SST)
k -
ω湍流模型(
22]。这个模型包含依赖当地的紊流长度尺度。模型常数
C
D
E
年代设置为推荐值Strelets et al。科氏力和离心力被认为是惯性力项的相对坐标系统。跑和莱斯地区叶轮和扩压器通道如图
2。这个DES仿真代码可以变换RANS-mode墙内进入LES-mode附近地区流体体积。
叶轮和扩压器通道内跑/ LES地区。
叶轮流道
扩散器通道
3.2。计算域
在数值模拟计算域采用如图
3。网格系统包括所有的流动通道(14叶轮和15扩散通道)。计算域分为三个区域:移动叶轮,静止的扩散器,叶尖间隙区域移动。蜗壳地区是不包括在这个仿真。整个域包括叶尖间隙区域,多次拉丝的网格生成结构化网格使用AutoGrid5 ver9.1 (NUMECA国际)。在整个环分析,叶轮(包括间隙区域)和扩散器系统26.2和2480万个细胞,分别。计算网格包含七种主要和分流叶片和15扩散叶片和有5100万个细胞。细胞墙的宽度为0.1
μ米,这对应于一个
y +参数是约等于一个沿固体表面。从电网依赖测试部分环分析,网格密度足够为了捕捉的不稳定行为扩散叶片附近的漩涡。
压缩机性能的实验和数值结果表现出图
4。数值DES的总压力上升特征得到不稳定分析。和总压力上升流动系数定义如下:
(1)
ϕ
=
60
问
π
2
D
2
2
B
2
N
,
Ψ
t
=
60
2
Δ
p
t
ρ
π
2
D
2
2
N
2
/
2
。
压缩机的性能。
ODV实验结果用图中的黑色圆圈。的七个操作点进行了模拟设计点(
ϕ
o
p
t
=
0.24)到深失速点(
ϕ= 0.10)。定结果不稳定的DES分析获得的测量结果具有很好的一致性,而且使用ODV,除了结果失速点(
ϕ= 0.14和0.10)。原因之一是,额外的损失从蜗壳数值模拟没有得出结论。特别是在深失速点(
ϕ= 0.10),压缩机性能数据获得的数值,利用CFD,比实验,获得更高的价值。结果关于性能非常不匹配的问题
ϕ= 0.10详细考虑后。
此外,瞬时套管压力波动在无因次时间在无叶片的空间
t
∗见图-31.5 = 0.0
7。纵轴是每个扩散叶片,标签与数字从1到15。套管静压较低地区,由一个黑色虚线圆,传播到下一个随着时间的推移扩散通道。这个地区的旋转速度大约是叶轮转速的25%。因此,它可能会得出结论,这种低静压地区是由细胞扩散器旋转失速。此外,扩散器入口的径向速度波动测量如图
8。红色和蓝色线代表了扩散通道数字1 - 2和7 - 8,分别。红色虚线,如图
8表明,径向速度的平均值计算扩散通道。如这个图所示,扩散通道数字7 - 8的径向速度低于平均径向速度,因为扩散器旋转失速细胞接近扩散通道数字从无量纲时间7 - 8,
t
∗
一个
=
0.0,如图
7。在无量纲时间,
t
∗
c
=
4
。
5的径向速度扩散通道数字7 - 8明显减少。因此,这些结果表明扩散通道内的流动是停滞不前的堵塞效应扩散器失速细胞形成。
套管在扩散器入口静压波动从无量纲时间
t
∗
=
0.0来
t
∗
=
31.5(CFD)。
径向速度波动以扩散器入口(CFD)。
然后,瞬时旋转的结构低静压地区的无量纲时间,
t
∗
c
=
4
。
5,见图
9。涡结构在几个停滞扩散通道数字7 - 11,如图
7,由Q-definition可视化
24通过无量纲螺旋性)和彩色。在流场无量纲时间,
t
∗
c
=
4.5,tornado-type涡连接裹尸布墙和扩散器吸力面观察扩散叶片的吸力面号码7和9。还开发了纵向涡的扩压器叶片吸力面数字8和10。这个漩涡叫做列弗。列弗是开发和扩散器内形成大量堵塞通道在非设计工况点。此外,另一个旋转的结构被发现扩散叶片压力面附近的数字8和10。这漩涡形成堵塞咽喉区,诱导强势逆转流附近的裹尸布。因此,它是确定发达停滞细胞内形成旋涡的结构。
瞬时旋转的扩散器结构在扩散段7号到11号摊位。
理解的不稳定行为扩散器失速细胞,瞬时旋转的结构内扩散通道数字在无量纲时间7 - 8
t
∗
一个,
t
∗
b,
t
∗
c是可视化。无量纲的漩涡核心的螺旋性扩散通道内数字7 - 8如图
10。首先,tornado-type漩涡的吸力面扩散叶片是在无量纲时间生成的
t
∗
一个
=
0.0。tornado-type分离是由积极的迎角附近的扩散叶片裹尸布,这导致了停滞细胞接近扩散通道数字7 - 8。然后,tornado-type涡是发达国家和列弗是观察到裹尸布/吸力面角无量纲时间
t
∗
b
=
2.1。发达列弗是扩展到下一个扩散器叶片(8)因为不良高扩散通道内的压力梯度。此外,观察LEV腿的吸力面扩散叶片数量7。列弗是由于二次流的影响扩散叶片吸力面。吸力面附近的流卷起的裹尸布边,因为垂直压力梯度,形成了纵向涡在裹尸布/吸力面角。在无量纲时间,
t
∗
c
=
4.5,列弗与下一个扩散器前缘,形成了喉咙堵塞,引起强烈的反转流附近的裹尸布。因此,分离涡生成被认为是在扩压器叶片吸力面8号随着时间的推移,由于入射角扩散叶片8号是非常积极的。新成立的tornado-type涡和列弗也观察到扩散叶片前缘附近的7号。
接下来,停滞阶段细胞的产生机制是通过实验分析调查
ϕ= 0.10。速度和压力进行了测量,揭示了详细的流场,舞台失速发生时。图
14显示了径向速度的测量结果和汽缸壁压力痕迹(30 Hz低通滤波器)在扩散器入口。低通滤波器由傅里叶反变换的测量光谱减少高频范围。顶部左边和右边的数字代表中心端的径向速度跟踪和裹尸布,分别。左边两个痕迹,同时正确的测量。这些数据表明,径向逆流基本上发生在裹尸布边的强度压力波动小的无叶片的空间(蓝色背景颜色:条件(a))。然而,裹尸布上的径向速度突然加速和减速轮毂一侧(红色背景颜色:条件(b))。然后,还有一个大强度的压力和速度波动(绿色背景颜色:条件(c))。截止频率(30 Hz)接近观察旋转摊位频率(25 Hz和22赫兹)。然而,截止频率是合适的,因为我们专注于整个间歇过程(条件(a) - (c))。此外,高空impeller-discharge流的径向速度分布的条件(a)和(b)在D.I.1-D.I测量。4使用由两个裂膜探头,如图
15。探测器设置在中心端为了捕捉失速条件(a) (c)。另一个速度探头是测量速度分布在每个条件。流动条件(a),裹尸布的一面都是逆转点。此外,在D.I.1-D.I逆流。4moved from the shroud side to the hub side on condition (b). Thus, it was assumed that the flow condition in vanless space at
ϕ= 0.10分为三个行为和交换使不稳定。
径向速度和压力波动痕迹扩压器进口(
ϕ= 0.10,30 Hz低通滤波器)。
高空impeller-discharge流的径向速度分布情况(a)和(b)。
的功率谱速度波动在每个条件无叶片的空间如图
16。速度波动(以25 Hz)引起的扩散器失速附近发现了裹尸布边条件(a)。然后,扩散转移到中心旋转波动方面条件(b)。最后,速度的大强度波动造成的停滞阶段被发现在整个跨度,特别是在中心附近的一面。此外,低通滤波速度跟踪分析小波变换来估计旋转不稳定行为的摊位出现在压缩机。速度信号的小波变换,
u
(
t
),定义如下:
(2)
W
一个
,
b
=
1
一个
∫
ψ
∗
T
u
t
d
t
,在哪里
t是时间,
ψ
∗
(
T
)是一个母小波,
一个是一个比例因子,
b是一个时间变换因子,然后呢
T
=
(
t
- - - - - -
b
)
/
一个。在这部作品中,伽柏母小波
(3)
ψ
∗
T
=
1
2
π
一个
e
- - - - - -
T
2
/
2
σ
2
e
- - - - - -
我
T是申请获取旋转失速扰动。在这里,恒定的参数,
σ,被定义为
σ= 32这项工作。图
17给出了小波变换的等高线
W
(
一个
,
b
)低通滤波的径向速度在中心端痕迹。纵轴表示频率,
f。径向速度减速时,扩散器失速波动出现。25 Hz的非定常扰动消失了,和停滞阶段波动条件(c)被发现。最后,停滞阶段波动不存在条件(a)。因此,它被认为是扩散器失速旋转在裹尸布边突然搬到中心,演变成停滞阶段。
速度在高空波动方向分布在三个条件(一)~ (c)。
小波变换的等高线图
W
(
一个
,
b
)低通滤波的径向速度在中心端痕迹。
图
18显示了同步测量的结果无叶片的径向速度的空间和汽缸壁压力痕迹叶轮通道内(30 Hz low-ass过滤器)。上面的图是径向速度跟踪中心一侧。底部的三个数字代表了套管压力子午方向的痕迹,I.I.,即时消息,D。我。一个。From the results at I.I. and I.M., the scale of pressure fluctuation in condition (b) was larger than condition (a), because the diffuser stall was considered to be expanded to impeller passages. In addition, the stage stall fluctuation was also found in the impeller passages. Figure
19显示每个流动状态的功率谱测量叶轮子午方向的裹尸布墙。红、绿、蓝线表示I.I.的压力波动水平,即时消息,D。我。一个。From these results, the diffuser stall cell was not found within the impeller passages on the flow condition (a). However, the pressure fluctuation of the diffuser stall occurred within the impeller passages on condition (b). Thus, the diffuser stall was considered to expand into the impeller passages before the stage stall occurred. Finally, the stage stall cell was formed by the expansion of the diffuser stall on condition (c). Furthermore, the pressure fluctuations measured at 55 Hz were found from the power spectra of I.I. on conditions (a) and (b). Based on the previous report [
12,
13),这55 Hz压力波动是由于叶轮旋转失速。轴流式压缩机,据报道,阻碍细胞,由多个漩涡,细胞合并成一个单一的漩涡。单涡细胞阻塞流在转子和定子的段落的质量流率降低(
14]。建议停滞细胞结构在轴向压缩器的早期阶段是类似于上面的扩散阻碍细胞观察。CFD分析能够捕捉叶轮和扩压器失速细胞;然而扩散器的发展停滞阶段停滞不存在流量系数
ϕ= 0.10。因此,实验和CFD结果发生在性能不匹配
ϕ= 0.10,因为停滞阶段压缩机性能恶化。详细的实验和CFD分析将在未来的研究进行检查叶轮和扩压器/阶段停滞细胞之间的关系。