<p>This investigation examines instability mechanisms and turbulent transition in the stationary disk boundary layer of a rotor-stator cavity through an integrated approach combining two-dimensional direct numerical simulations (2D DNS), global linear stability analysis, and three-dimensional direct numerical simulations (3D DNS). The base flow characteristics, the global linear stability properties, and the nonlinear evolution processes leading to turbulence are explored. The research examines an annular, enclosed rotor-stator cavity with curvature parameter <InlineEquation ID="IEq1"> <EquationSource Format="TEX">\(Rm = (b^*+a^*)/(b^*-a^*)=1.8\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>R</mi> <mi>m</mi> <mo>=</mo> <mrow> <mo stretchy="false">(</mo> <msup> <mi>b</mi> <mo>∗</mo> </msup> <mo>+</mo> <msup> <mi>a</mi> <mo>∗</mo> </msup> <mo stretchy="false">)</mo> </mrow> <mo stretchy="false">/</mo> <mrow> <mo stretchy="false">(</mo> <msup> <mi>b</mi> <mo>∗</mo> </msup> <mo>-</mo> <msup> <mi>a</mi> <mo>∗</mo> </msup> <mo stretchy="false">)</mo> </mrow> <mo>=</mo> <mn>1.8</mn> </mrow> </math></EquationSource> </InlineEquation>, and an aspect ratio <InlineEquation ID="IEq2"> <EquationSource Format="TEX">\(L=(b^*-a^*)/(2h^*) = 5\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>L</mi> <mo>=</mo> <mrow> <mo stretchy="false">(</mo> <msup> <mi>b</mi> <mo>∗</mo> </msup> <mo>-</mo> <msup> <mi>a</mi> <mo>∗</mo> </msup> <mo stretchy="false">)</mo> </mrow> <mo stretchy="false">/</mo> <mrow> <mo stretchy="false">(</mo> <mn>2</mn> <msup> <mi>h</mi> <mo>∗</mo> </msup> <mo stretchy="false">)</mo> </mrow> <mo>=</mo> <mn>5</mn> </mrow> </math></EquationSource> </InlineEquation>, where <InlineEquation ID="IEq3"> <EquationSource Format="TEX">\(a^*\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>a</mi> <mo>∗</mo> </msup> </math></EquationSource> </InlineEquation> and <InlineEquation ID="IEq4"> <EquationSource Format="TEX">\(b^*\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>b</mi> <mo>∗</mo> </msup> </math></EquationSource> </InlineEquation> are the minimum and maximum radius, respectively, and <InlineEquation ID="IEq5"> <EquationSource Format="TEX">\(2h^*\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>2</mn> <msup> <mi>h</mi> <mo>∗</mo> </msup> </mrow> </math></EquationSource> </InlineEquation> denotes the inter-disk spacing. The current study examines Reynolds numbers <InlineEquation ID="IEq6"> <EquationSource Format="TEX">\(Re_\varphi = \Omega ^*_db^{*2}/\nu ^*\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>R</mi> <msub> <mi>e</mi> <mi>φ</mi> </msub> <mo>=</mo> <msubsup> <mi mathvariant="normal">Ω</mi> <mi>d</mi> <mo>∗</mo> </msubsup> <mmultiscripts> <mi>b</mi> <mrow /> <mrow> <mrow /> <mo>∗</mo> <mn>2</mn> </mrow> </mmultiscripts> <mo stretchy="false">/</mo> <msup> <mi>ν</mi> <mo>∗</mo> </msup> </mrow> </math></EquationSource> </InlineEquation> ranging from <InlineEquation ID="IEq7"> <EquationSource Format="TEX">\(10^4\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mn>10</mn> <mn>4</mn> </msup> </math></EquationSource> </InlineEquation> to <InlineEquation ID="IEq8"> <EquationSource Format="TEX">\(10^5\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mn>10</mn> <mn>5</mn> </msup> </math></EquationSource> </InlineEquation>, where <InlineEquation ID="IEq9"> <EquationSource Format="TEX">\(\Omega ^*_d\)</EquationSource> <EquationSource Format="MATHML"><math> <msubsup> <mi mathvariant="normal">Ω</mi> <mi>d</mi> <mo>∗</mo> </msubsup> </math></EquationSource> </InlineEquation> is the angular velocity of the rotor, <InlineEquation ID="IEq10"> <EquationSource Format="TEX">\(\nu ^*\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>ν</mi> <mo>∗</mo> </msup> </math></EquationSource> </InlineEquation> is the fluid’s kinematic viscosity. The results reveal that for <InlineEquation ID="IEq11"> <EquationSource Format="TEX">\(Re_\varphi \le 3.6 \times 10^4\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>R</mi> <msub> <mi>e</mi> <mi>φ</mi> </msub> <mo>≤</mo> <mn>3.6</mn> <mo>×</mo> <msup> <mn>10</mn> <mn>4</mn> </msup> </mrow> </math></EquationSource> </InlineEquation>, the base flow converges to a stable equilibrium state. However, when <InlineEquation ID="IEq12"> <EquationSource Format="TEX">\(Re_\varphi \)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>R</mi> <msub> <mi>e</mi> <mi>φ</mi> </msub> </mrow> </math></EquationSource> </InlineEquation> exceeds <InlineEquation ID="IEq13"> <EquationSource Format="TEX">\(3.7 \times 10^4\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>3.7</mn> <mo>×</mo> <msup> <mn>10</mn> <mn>4</mn> </msup> </mrow> </math></EquationSource> </InlineEquation>, the base flow destabilizes, giving rise to persistent velocity oscillations within the stationary disk boundary layer. Global linear stability analysis confirms these observations. For axisymmetric disturbances (azimuthal wavenumber <InlineEquation ID="IEq14"> <EquationSource Format="TEX">\(\beta = 0\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>β</mi> <mo>=</mo> <mn>0</mn> </mrow> </math></EquationSource> </InlineEquation>), the critical Reynolds number for global linear instability is <InlineEquation ID="IEq15"> <EquationSource Format="TEX">\(Re_\varphi = 3.67 \times 10^4\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>R</mi> <msub> <mi>e</mi> <mi>φ</mi> </msub> <mo>=</mo> <mn>3.67</mn> <mo>×</mo> <msup> <mn>10</mn> <mn>4</mn> </msup> </mrow> </math></EquationSource> </InlineEquation>, precisely demarcating the threshold for self-sustained oscillations in the base flow. For non-axisymmetric disturbances with <InlineEquation ID="IEq16"> <EquationSource Format="TEX">\(\beta = 15\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>β</mi> <mo>=</mo> <mn>15</mn> </mrow> </math></EquationSource> </InlineEquation>, a more complex dynamics is identified: the maximum linear temporal growth rate occurs at <InlineEquation ID="IEq17"> <EquationSource Format="TEX">\(Re_\varphi = 2.5 \times 10^4\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>R</mi> <msub> <mi>e</mi> <mi>φ</mi> </msub> <mo>=</mo> <mn>2.5</mn> <mo>×</mo> <msup> <mn>10</mn> <mn>4</mn> </msup> </mrow> </math></EquationSource> </InlineEquation>, while the critical threshold for global linear instability emerges at a lower <InlineEquation ID="IEq18"> <EquationSource Format="TEX">\(Re_\varphi = 1.58 \times 10^4\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>R</mi> <msub> <mi>e</mi> <mi>φ</mi> </msub> <mo>=</mo> <mn>1.58</mn> <mo>×</mo> <msup> <mn>10</mn> <mn>4</mn> </msup> </mrow> </math></EquationSource> </InlineEquation>, marking the onset conditions for a self-sustained spiral wave mode. The 3D DNS results reveal sophisticated nonlinear dynamics where spiral wave energy is transferred from higher to lower radial positions through complex mode interactions. When circular waves dominate in the flow field, their interaction with spiral waves across various radial positions generates a rich spectral composition with multiple distinct dominant frequencies within the boundary layer. At <InlineEquation ID="IEq19"> <EquationSource Format="TEX">\(Re_\varphi \ge 7.0 \times 10^4\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>R</mi> <msub> <mi>e</mi> <mi>φ</mi> </msub> <mo>≥</mo> <mn>7.0</mn> <mo>×</mo> <msup> <mn>10</mn> <mn>4</mn> </msup> </mrow> </math></EquationSource> </InlineEquation>, the flow transitions to a chaotic state with turbulent characteristics. Remarkably, even in this strongly chaotic regime, coherent circular wave structures remain observable within the chaotic flow field, although distinct dominant frequencies are no longer present in the spectral analysis.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Instability and transition to turbulence of the stationary-disk boundary layer of a rotor-stator cavity

  • Lei Xie,
  • Siyi Li,
  • Yaguang Xie,
  • Qiang Du,
  • Junqiang Zhu,
  • Ruonan Wang,
  • Qingzong Xu,
  • Zihao Zhu

摘要

This investigation examines instability mechanisms and turbulent transition in the stationary disk boundary layer of a rotor-stator cavity through an integrated approach combining two-dimensional direct numerical simulations (2D DNS), global linear stability analysis, and three-dimensional direct numerical simulations (3D DNS). The base flow characteristics, the global linear stability properties, and the nonlinear evolution processes leading to turbulence are explored. The research examines an annular, enclosed rotor-stator cavity with curvature parameter \(Rm = (b^*+a^*)/(b^*-a^*)=1.8\) R m = ( b + a ) / ( b - a ) = 1.8 , and an aspect ratio \(L=(b^*-a^*)/(2h^*) = 5\) L = ( b - a ) / ( 2 h ) = 5 , where \(a^*\) a and \(b^*\) b are the minimum and maximum radius, respectively, and \(2h^*\) 2 h denotes the inter-disk spacing. The current study examines Reynolds numbers \(Re_\varphi = \Omega ^*_db^{*2}/\nu ^*\) R e φ = Ω d b 2 / ν ranging from \(10^4\) 10 4 to \(10^5\) 10 5 , where \(\Omega ^*_d\) Ω d is the angular velocity of the rotor, \(\nu ^*\) ν is the fluid’s kinematic viscosity. The results reveal that for \(Re_\varphi \le 3.6 \times 10^4\) R e φ 3.6 × 10 4 , the base flow converges to a stable equilibrium state. However, when \(Re_\varphi \) R e φ exceeds \(3.7 \times 10^4\) 3.7 × 10 4 , the base flow destabilizes, giving rise to persistent velocity oscillations within the stationary disk boundary layer. Global linear stability analysis confirms these observations. For axisymmetric disturbances (azimuthal wavenumber \(\beta = 0\) β = 0 ), the critical Reynolds number for global linear instability is \(Re_\varphi = 3.67 \times 10^4\) R e φ = 3.67 × 10 4 , precisely demarcating the threshold for self-sustained oscillations in the base flow. For non-axisymmetric disturbances with \(\beta = 15\) β = 15 , a more complex dynamics is identified: the maximum linear temporal growth rate occurs at \(Re_\varphi = 2.5 \times 10^4\) R e φ = 2.5 × 10 4 , while the critical threshold for global linear instability emerges at a lower \(Re_\varphi = 1.58 \times 10^4\) R e φ = 1.58 × 10 4 , marking the onset conditions for a self-sustained spiral wave mode. The 3D DNS results reveal sophisticated nonlinear dynamics where spiral wave energy is transferred from higher to lower radial positions through complex mode interactions. When circular waves dominate in the flow field, their interaction with spiral waves across various radial positions generates a rich spectral composition with multiple distinct dominant frequencies within the boundary layer. At \(Re_\varphi \ge 7.0 \times 10^4\) R e φ 7.0 × 10 4 , the flow transitions to a chaotic state with turbulent characteristics. Remarkably, even in this strongly chaotic regime, coherent circular wave structures remain observable within the chaotic flow field, although distinct dominant frequencies are no longer present in the spectral analysis.