T. Akao, T. Watanabe, K. Nagata
Vertical confinement effects on a fully developed turbulent shear layer
Physics of Fluids, 34 055129 2022

The PDF is available here.  
This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article may be found at https://doi.org/10.1063/5.0090686.

Abstract

The effects of vertical confinement on a turbulent shear layer are investigated with large-eddy simulations of a freely developing shear layer (FSL) and a wall-confined shear layer (WSL) that develops between two horizontal walls. In the case of the WSL, the growth of the shear layer is inhibited by the walls. Once the walls prevent the development of the shear layer, highly anisotropic velocity fluctuations become prominent in the flow. These anisotropic velocity fluctuations are recognized as elongated large-scale structures (ELSS), whose streamwise length is much larger than the length scales in the other directions. Spectral analysis confirms that the turbulent kinetic energy is dominated by the ELSS, whose streamwise length grows continuously. A proper orthogonal decomposition can effectively extract a velocity component associated with the ELSS. The isotropy of the Reynolds stress tensor is changed by the presence of the ELSS. These changes in flow characteristics due to the ELSS are not observed in the FSL, where the shear layer thickness increases continuously. These behaviors of the WSL are consistent with those of stably stratified shear layers (SSSLs), where flow structures similar to ELSS also develop when the vertical flow development is confined by the stable stratification. The vertical confinement by the walls or stable stratification strengthens mean shear effects. The flow behavior at large scales in the WSL and SSSL is consistent with rapid distortion theory for turbulence subject to mean shear, suggesting that the development of ELSS is caused by the mean shear.

日本語訳 (DeepL翻訳)

完全発達した乱流せん断層における垂直方向の拘束効果

乱流せん断層に対する鉛直拘束の効果を、自由発達せん断層(FSL)と二つの水平壁間に発達する壁拘束せん断層(WSL)のラージ・エディ・シミュレーションによって調べた。壁拘束せん断層の場合、せん断層の成長は壁によって阻害される。壁によってせん断層の発達が阻害されると、流れに非等方的な速度変動が発生するようになる。この非等方的な速度変動は、流れ方向の長さが他の方向の長さスケールよりはるかに大きい細長い大規模構造（ELSS）として認識される。スペクトル解析の結果、乱流の運動エネルギーは、長さが連続的に増加するELSSによって支配されていることが確認された。固有直交分解により、ELSSに関連する速度成分を効果的に抽出することができる。また、レイノルズ応力テンソルの等方性は、ELSSの存在によって変化することがわかった。このようなELSSによる流れの特性の変化は、せん断層の厚さが連続的に増加するFSLでは観測されない。これらのWSLの挙動は、安定成層によって鉛直方向の流れの発達が拘束されるとELSSと同様の流動構造が発達する安定成層中のせん断層（SSSL）の挙動と一致する。壁や安定成層による鉛直方向の拘束は、平均剪断効果を強める。WSLとSSSLにおける大きなスケールの流れの挙動は、平均せん断を受ける乱流に対するrapid distortion theoryと一致し、ELSSの発達が平均せん断に起因することを示唆するものであった。

K. Yamamoto, T. Ishida, T. Watanabe, K. Nagata
Experimental and numerical investigation of compressibility effects on velocity derivative flatness in turbulence
Physics of Fluids, 34 055101 2022

The PDF is available here.  
This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article may be found at https://doi.org/10.1063/5.0085423.

Abstract

Compressibility effects on the velocity derivative flatness 𝐹∂𝑢′/∂𝑥 are investigated by experiments with opposing arrays of piston-driven synthetic jet actuators (PSJAs) and direct numerical simulations (DNS) of statistically steady compressible isotropic turbulence and temporally evolving turbulent planar jets with subsonic or supersonic jet velocities. Experiments using particle image velocimetry show that nearly homogeneous isotropic turbulence is generated at the center of a closed box from interactions between supersonic synthetic jets. The dependencies of 𝐹∂𝑢′/∂𝑥 on the turbulent Reynolds number 𝑅𝑒𝜆 and the turbulent Mach number MT are examined both experimentally and using DNS. Previous studies of incompressible turbulence indicate a universal relationship between 𝐹∂𝑢′/∂𝑥 and 𝑅𝑒𝜆. However, both experiments and DNS confirm that 𝐹∂𝑢′/∂𝑥 increases relative to the incompressible turbulence via compressibility effects. Although 𝐹∂𝑢′/∂𝑥 tends to be larger with MT in each flow, the 𝐹∂𝑢′/∂𝑥 in the turbulent jets and the turbulence generated from PSJAs deviate from those in incompressible turbulence at lower MT compared with isotropic turbulence sustained by a solenoidal forcing. The PSJAs and supersonic planar jets generate strong pressure waves, and the wave propagation can cause an increased 𝐹∂𝑢′/∂𝑥, even at low MT. These results suggest that the compressibility effects on 𝐹∂𝑢′/∂𝑥 are not solely determined from a local value of MT and depend on the turbulence generation process.

日本語訳 (DeepL翻訳)

乱流中の速度微分平坦度に対する圧縮性効果の実験的および数値的検討

速度微分平坦度𝐹∂𝑢′/∂𝑥に対する圧縮性の影響を、対向するピストン駆動合成ジェットアクチュエータ(PSJA)を用いた実験と、統計的に定常な圧縮性等方性乱流と時間発展する乱流平面噴流の直接数値シミュレーション(DNS)によって調査した。粒子画像流速測定法を用いた実験により、閉じたボックスの中心で超音速合成噴流の相互作用からほぼ一様な等方性乱流が生成されることが示された。乱流レイノルズ数𝑅𝑒𝜆と乱流マッハ数MTに対する𝐹∂𝑢′/∂𝑥の依存性を実験とDNSにより調べた。非圧縮性乱流に関するこれまでの研究では、𝐹∂𝑢′/∂𝑥と𝑅𝑒𝜆の間に普遍的な関係があることが示されている。しかし、実験とDNSの両方で、非圧縮性乱流に対して圧縮性効果により𝐹∂𝑢′/∂𝑥が増加することが確認された。それぞれの流れにおいて、MTが大きくなると𝐹∂𝑢′/∂𝑥が大きくなる傾向があるが、乱流噴流とPSJAの乱流はソレノイド加振による等方性乱流と比較して、低MTで無圧縮乱流の𝐹∂𝑢′/∂𝑥から外れた値となった。PSJAと超音速平面噴流は強い圧力波を発生させ、その波動伝播によって低MTでも𝐹∂𝑢′/∂𝑥が増加することがある。これらの結果は、𝐹∂𝑢′/∂𝑥に対する圧縮性の効果はMTの局所的な値のみから決まるものではなく、乱流生成過程に依存することを示唆している。

Y. Zheng, K. Nagata, T. Watanabe
Turbulent characteristics and energy transfer in the far field of active-grid turbulence
Physics of Fluids, 33 115119 2021

The PDF is available here.  
This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article may be found at https://doi.org/10.1063/5.0070912.

Abstract

Turbulent characteristics in the far field of active-grid turbulence have been investigated through wind tunnel experiments using hot-wire anemometry. Two forcing protocols are employed following previous studies: one is the double-random mode and the other is the open mode with the grid remaining static with minimum blockage. The integral length scale L for the double-random modes slightly decreases with streamwise distance in the far field as observed in the near field of the active-grid turbulence. The nondimensional dissipation rate 𝐶𝜀 for the double-random modes is around 0.5. This asymptotic value is different from those reported in previous active-grid turbulence experiments and could be nonuniversal. The equilibrium scaling 𝐿/𝜆=𝐶𝜀𝑅𝑒𝜆/15 (λ is the Taylor microscale and 𝑅𝑒𝜆 is the turbulent Reynolds number) with a constant 𝐶𝜀 is established in the far field of the double-random modes regardless of active-grid motions. The sum of production and destruction terms in the enstrophy budget equation for homogeneous and isotropic turbulence 𝑆+2𝐺/𝑅𝑒𝜆 (S is the skewness of the longitudinal velocity derivative and G is the destruction coefficient) is proportional to 𝑅𝑒𝜆^−1 and close to zero in the present active-grid turbulence, suggesting that the equilibrium scaling is possibly related to the balance between the production and destruction of the enstrophy.

日本語訳 (DeepL翻訳)

アクティブ格子乱流の遠方場における乱流特性とエネルギー輸送

アクティブ格子乱流の遠方場における乱流特性を、熱線流速計を用いた風洞実験により調べた。一つはダブルランダムモードで，もう一つはグリッドが最小限の閉塞で静止しているオープンモードである．ダブルランダムモードの積分長さスケールLは、アクティブ格子乱流の近傍場で観測されたように、遠方場では流向距離とともにわずかに減少することがわかった。ダブルランダムモードの無次元散逸率𝐶𝜀は 0.5 程度である. この漸近的な値は、過去の活動格子乱流実験で報告された値とは異なっており、普遍的でない可能性がある。ダブルランダムモードの遠距離場では、アクティブ格子の運動によらず、 𝐿/𝜆=𝐶𝜀𝑅𝑒𝜆/15 (λ はテイラー・マイクロスケール、𝑅𝑒𝜆 は乱流レイノルズ数) の平衡スケーリングが成立している。均質・等方乱流のエンス トロフィー方程式における生成項と破壊項の和𝑆+2𝐺 (S は縦速度微分の歪度、G は破壊係数) は𝑅𝜆-1 に比例し、現在のアクティブ格子乱流では 0 に近い値になっている。は、平衡スケーリングがエンストロフィーの生成と破壊のバランスに関係している可能性を示唆する。

T. Watanabe, K. Tanaka, K. Nagata
Solenoidal linear forcing for compressible, statistically steady, homogeneous isotropic turbulence with reduced turbulent Mach number oscillation
Physics of Fluids, 33 095108 2021

The PDF is available here.  
This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. 

This article may be found at https://doi.org/10.1063/5.0062596.

Abstract

This study investigates a solenoidal linear forcing scheme with reduced oscillation of a turbulent Mach number MT for direct numerical simulations (DNS) of statistically steady, homogeneous isotropic turbulence. A conventional linear forcing scheme results in a large temporal oscillation of MT, where the maximum MT reaches about 1.1 times the time-averaged MT. Therefore, strong shocklets are generated when MT becomes large although such strong shocklets hardly appear when MT is close to the time-averaged value. DNS with the proposed forcing scheme confirms that the temporal oscillation of MT is effectively reduced by adjusting a forcing coefficient with a ratio between velocity variance and its steady state value prescribed as a parameter. The time-dependent forcing coefficient results in the variation of the power input to kinetic energy. Therefore, the temporal oscillation of the Reynolds number for this forcing scheme is as large as that for the conventional linear forcing. The ratio between the solenoidal and dilatational kinetic energy dissipation rates increases with MT, and the MT dependence is consistent between the present solenoidal linear forcing and the low-wavenumber solenoidal forcing in wavenumber space. The skewness and flatness of the velocity derivative become large compared with incompressible turbulence when MT exceeds 0.6. Both average and root-mean-squared fluctuation of the shock Mach number of shocklets increase with MT. The most typical thickness of shocklets decreases with MT and asymptotically approaches about 1.5 times the Kolmogorov scale. The shocklet thickness normalized by the Kolmogorov scale hardly depends on the Reynolds number.

日本語訳 (DeepL翻訳)

統計的に定常な圧縮性一様等方性乱流に対する乱流マッハ数変動を低減してソレノイド線形加振スキーム

本研究では，統計的に定常な等方性乱流の直接数値シミュレーション（DNS）において，乱流マッハ数MTの振動を低減したソレノイド線形加振スキームを検討する．従来の線形加振スキームでは、MTの時間振動が大きく、最大MTは時間平均MTの1.1倍程度に達する。そのため、MTが時間平均に近いときにはほとんど発生しないが、MTが大きくなると強い衝撃波が発生する。提案した加振スキームを用いたDNSでは、速度変動と定常値の比をパラメータとする加振係数を調整することで、MTの時間的振動を効果的に低減できることが確認された。時間依存の加振係数は、運動エネルギーへの入力パワーを変化させることになる。そのため、この強制方式ではレイノルズ数の時間的振動が従来の線形強制と同程度に大きくなる。運動エネルギー散逸率の発散・非発散成分の比はMTとともに増加し、MT依存性は波数空間において現在のソレノイド線形加振と低波数ソレノイド加振の間で一致することが示された。MTが0.6を超えると、非圧縮性乱流と比較して速度微分の歪み度と平坦度が大きくなることがわかった。衝撃波のマッハ数の平均値と二乗平均の揺らぎはMTとともに増加する。最も典型的な衝撃波の厚さはMTとともに減少し、漸近的にKolmogorovスケールの約1.5倍に近づく。Kolmogorovスケールで規格化した衝撃波の厚さはレイノルズ数にはほとんど依存しない。

T. Katagiri, T. Watanabe, K. Nagata
Statistical properties of a model of a turbulent patch arising from a breaking internal wave
Physics of Fluids, 33 055107 2021

The PDF is available here.  
This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article may be found at https://doi.org/10.1063/5.0046832.

Abstract

The turbulent patch arising from internal gravity wave breaking is investigated with direct numerical simulation of a stably stratified flow over a two-dimensional hill. The turbulent patch is distinguished from the non-turbulent wave region with potential vorticity. The turbulent patch is highly intermittent, and its location fluctuates with space and time. The buoyancy Reynolds number slowly decays with time in the turbulent patch and the mixing efficiency stays around 0.2. The turbulent patch is separated from the non-turbulent wave region by a turbulent/non-turbulent interfacial (TNTI) layer, whose thickness is about five times the Kolmogorov scale. The kinetic energy dissipation rate also sharply decreases from the turbulent to the wave region while the potential energy dissipation rate has a large peak within the TNTI layer. Both shear and stable stratification are strong in the upper area of the turbulent patch. On the other hand, the lower area has a small mean density gradient, i.e., weak stratification, which is related to the strong intermittency of the turbulent patch in the lower area. Furthermore, weak stratification in the lower area results in a low gradient Richardson number, which is below the critical value for the shear instability, and the roller vortex appears. The outer edge of the turbulent patch aligns with the perimeter of the roller vortex, and the vortex affects the spatial distribution of the turbulent patch.

日本語訳 (DeepL翻訳)

砕波した内部波から発生する乱流パッチモデルの統計的性質

内部重力波の砕波によって生じる乱流パッチを，2次元の丘陵上の安定成層流の直接数値シミュレーションによって調べた．乱流パッチはポテンシャル渦度を持つ非乱流波動領域と区別される。乱流パッチは非常に間欠的であり、その位置は空間的、時間的に変動する。乱流パッチでは浮力レイノルズ数は時間とともにゆっくりと減少し、混合効率は0.2程度にとどまる。乱流パッチは非乱流波浪領域から厚さコルモゴロフスケールの5倍程度の乱流・非乱流界面（TNTI）層で隔てられている。運動エネルギー散逸率も乱流領域から波浪領域へ向かって急激に減少するが、位置エネルギー散逸率はTNTI層内で大きなピークを持つ。乱流パッチの上部ではシアーと安定成層がともに強い。一方、下部領域では平均密度勾配が小さく、すなわち成層が弱い。これは、下部領域での乱流パッチの強い間欠性に関連している。さらに、下部領域の成層が弱いため、勾配リチャードソン数が低く、シアー不安定性の臨界値以下となり、ローラー渦が出現する。乱流パッチの外縁はローラー渦の外周に一致し、渦は乱流パッチの空間分布に影響を与える。

M. Hayashi, T. Watanabe, K. Nagata
The relation between shearing motions and the turbulent/non-turbulent interface in a turbulent planar jet
Physics of Fluids, 33 055126 2021

The PDF is available here.  
This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. 

This article may be found at https://doi.org/10.1063/5.0045376.

Abstract

The relation between shearing motions and the turbulent/non-turbulent interfacial (TNTI) layer is studied with direct numerical simulation of a temporally evolving planar jet. Small-scale shear layers are detected with the triple decomposition of the velocity gradient tensor, which is decomposed into shear, rotation, and elongation tensors. The shear layers are found in the turbulent sublayer more frequently than in the turbulent core region although they hardly appear in the viscous superlayer. The shear layers undergo a biaxial strain with stretching in the shear vorticity direction and compression in the interface normal direction. This compressive strain is related to the non-turbulent fluid, which is relatively advected toward the shear layer. The shear layer thickness in the TNTI layer is well predicted by Burgers vortex layer. The velocity jump of the shear layer is about seven times the Kolmogorov velocity both in the turbulent core region and the TNTI layer. However, the layer thickness normalized by the Kolmogorov scale is about 6 in the turbulent core region and decreases in the TNTI layer, where consequently, the shear Reynolds number becomes small. The shear layers have significant contributions to the enstrophy production in the turbulent sublayer and the viscous enstrophy-diffusion toward the viscous superlayer. The shear layer and the outer edge of the TNTI layer have a curvature radius of about 50 times the Kolmogorov scale. The alignment between the shear layer orientation and the interface normal direction confirms that the shear layers near the interface are mostly parallel to the TNTI layer.

日本語訳 (DeepL翻訳)

乱流平面噴流におけるせん断運動と乱流・非乱流界面の関係

時間発展する平面噴流の直接数値シミュレーションにより、剪断運動と乱流・非乱流界面（TNTI）層の関係を研究した。速度勾配テンソルを剪断、回転、伸長テンソルに分解する三成分分解により、小スケールの剪断層が検出される。せん断層は、粘性超層にはほとんど現れないが、乱流コア領域よりも乱流副層に多く存在することがわかった。せん断層は、せん断渦度方向に伸び、界面法線方向に圧縮される二軸性の歪みを受ける。この圧縮ひずみは、非乱流が相対的にシア層に向かって移流されることに関係している。TNTI層におけるせん断層の厚さは、Burgers渦層によってよく予測される。せん断層の速度ジャンプは、乱流コア領域とTNTI層の両方で、コルモゴロフ速度の約7倍であった。しかし、Kolmogorovスケールで正規化した層厚は乱流コア領域で約6、TNTI層で減少し、その結果、せん断レイノルズ数は小さくなった。せん断層は、乱流下層でのエンストロフィー生成と粘性超層に向かう粘性エンストロフィー拡散に大きく寄与している。せん断層とTNTI層外縁の曲率半径はKolmogorovスケールの約50倍である。せん断層の配向と界面法線方向の一致から、界面付近のせん断層はほとんどTNTI層と平行であることが確認された。