カーマン渦巻き通り

空気中の円筒の後ろの渦列の可視化(左)。円筒付近の空気中にグリセロール蒸気を放出することで流れを可視化している。

流体力学において、カルマン渦列(カルマンうずれれ、フォン・カルマンうずれれ)とは、渦放出と呼ばれる過程によって引き起こされる渦の繰り返しパターンであり鈍い物体の周りの流体流れの不安定な剥離の原因となる。[ 1 ]

これは、技術者で流体力学者のセオドア・フォン・カルマンにちなんで名付けられ、[ 2 ]吊り下げられた電話線や電力線の「歌声」や、特定の速度での自動車アンテナの振動などの現象の原因です。

フォン・カルマン渦列の数学的モデリングは、k-イプシロンSSTk-オメガ、レイノルズ応力を含むナビエ・ストークス方程式の完全解法やラージエディシミュレーション(LES)乱流モデルなど、さまざまな手法を用いて行うことができる。 [ 3 ] [ 4 ]ギンツブルグ・ランダウ方程式などの動的方程式を数値的に解くことによる方法、[ 5 ] [ 6 ] [ 7 ]または双複素変数の使用による方法もある。[ 8 ]

分析

円筒形の物体によって作られた渦列のアニメーション。物体の反対側の流れは異なる色で表示され、渦が物体の交互の側から放出されていることを示しています。
太平洋からモハーベ砂漠の山々を越えて東へ空気が急速に流れ込む様子を地上から見たカーマン渦列効果。雲に関連するカーマン渦列の活動は宇宙から観測されることが多いため、地上から観測されるこの現象は非常に稀です。
ハードディスクの2次元液体における渦列

渦列は、レイノルズ数( Re ) の範囲で指定される特定の流速範囲でのみ形成されます。通常、 Reの限界値は約 90 です。流れの (グローバル) レイノルズ数は、物体の周りまたはチャネル内の流体の流れにおける慣性力粘性力の比率の尺度であり、流体の流れ全体のグローバル速度の無次元パラメータ として定義できます。 ここで、 ReLあなたLν0ν0μ0ρ0ReLあなたLρ0μ0{\displaystyle {\begin{aligned}\mathrm {Re} _{L}&={\frac {UL}{\nu _{0}}}\\\nu _{0}&={\frac {\mu _{0}}{\rho _{0}}}\\\mathrm {Re} _{L}&={\frac {UL\rho _{0}}{\mu _{0}}}\\\end{aligned}}}

  • U = 自由流の流速(つまり、流体の境界から離れた流速、 U 、静止した流体に対する物体の速度、またはベルヌーイ方程式で計算された非粘性流速など)、これは元のグローバル フロー パラメーター、つまり無次元化されるターゲットです。
  • L = 本体またはチャネルの特性長さパラメータ
  • ν 0 =流体の 自由流動粘性パラメータ、これは次の比率です:
    • ρ 0 = 基準流体の密度。
    • μ 0 = 自由流流体の動粘度

一般的な流れ (通常は非圧縮性または等温流れと見なすことができます) では、動粘性は流れ場全体にわたってどこでも均一で、時間的に一定であるため、粘性パラメータについて選択の余地はなく、これは当然、検討対象の温度における流体の動粘性になります。一方、基準長さは常に任意のパラメータであるため、異なる障害物の周りの流れや異なる形状のチャネル内の流れを比較する場合には特別な注意を払う必要があります。全体的なレイノルズ数は、同じ基準長さを基準にする必要があります。実際にこれが、翼型およびチャネルの流れデータの最も正確な情報源がレイノルズ数で基準長さを指定している理由です。基準長さは、実行する解析によって変わる可能性があります。円筒や球などの円形断面を持つ物体の場合、通常は直径を選択します。翼、一般的な非円筒、鈍角体、胴体や潜水艦のような回転体の場合、実際には安定した設計入力となるのは、通常、プロファイル弦長、プロファイル厚さ、またはその他の特定の幅です。流路の場合、通常は流体が流れる 水力直径です。

空力プロファイルの場合、基準長さは解析によって異なります。実際、翼セクションや薄型プロファイルの空力係数の基準長さとしてもプロファイル弦が選択されることが多く、その主な目標は揚力係数または揚抗比を最大化することです (つまり、薄型翼理論では通常、異なるプロファイルを比較するための流速パラメータとして弦レイノルズを使用します)。一方、フェアリングやストラットの場合、与えられたパラメータは通常、流線型にする内部構造の寸法であり (簡単にするために円形断面の梁と考えます)、主な目標は抗力係数または抗力/揚力比を最小化することです。したがって、自然に基準長さにもなる主な設計パラメータは、プロファイル弦ではなく、プロファイルの厚さ (流れ方向に垂直なプロファイル寸法または面積) です。

Re値の範囲は、渦が発生する物体のサイズと形状、および流体の動粘性によって異なります。基準長さが通常円筒の直径dである円筒の後流の場合、この範囲の下限はRe ≈ 47です。 [ 9 ] [ 10 ]渦は円境界の各側から連続的に発生し、後流に渦の列を形成します。この交互配置により、1 つの列の渦の中心が、他の列の 2 つの渦の中心の中間点の反対側に位置することになり、図に示すような特徴的なパターンが生じます。最終的に、渦が下流に移動するにつれてそのエネルギーは粘性によって消費され、規則的なパターンは消えます。Re 値が 188.5 を超えると流れは 3 次元になり、円筒に沿って周期的に変化します。[ 11 ] Reが10 5オーダーを超えると抗力危機となり、渦の放出が不規則になり乱流が発生します。

一つの渦が放出されると、物体の周囲に非対称な流れのパターンが形成され、圧力分布が変化します。つまり、渦が交互に放出されることで、物体に周期的な横方向の力が加わり、振動を引き起こします。渦放出の周波数が物体や構造物の固有振動数と近い場合、共鳴が発生します。この強制振動が適切な周波数で発生すると、吊り下げられた電話線や送電線が「鳴る」ように振動したり、車のアンテナが特定の速度でより強く振動したりします。

気象学では

フアン・フェルナンデス
ソコロ
チリ沖のファン・フェルナンデス諸島の周りを流れる風によって生じたフォン・カルマン渦列(左)と、ソコロ島付近のフォン・カルマン渦の衛星ループ()。

島や孤立した山などの障害物の上空を大気が流れると、フォン・カルマン渦列が発生することがあります。該当する高度に雲層が存在する場合、この渦列は視認可能になります。このような雲層の渦列は衛星から撮影されており[ 12 ]、障害物から400km(250マイル)以上も到達することがあります。渦の直径は通常20~40km(12~25マイル)です[ 13 ] 。

エンジニアリングの問題

滑り止め円筒状障害物の周りの渦列のシミュレーション
同じシリンダーにフィンが付き、横渦が相互作用する領域を減らすことで渦列を抑制します。

乱流の少ない環境では、高層ビルは高さ方向に均一であれば、カルマン街路を形成できる。都市部では、近くに他の高層ビルが多数あるため、これらのビルによって生じる乱流によってコヒーレントな渦の形成が妨げられることがある。[ 14 ]物体の側面に渦が発生する周期的な横風力は、渦誘起振動を引き起こし、構造物を損傷する恐れがあるため、非常に望ましくない。そのため、潜水艦の潜望鏡から産業用煙突、高層ビルまで、さまざまな構造物を設計する際に、技術者が渦放出の影響を考慮することが重要である。このような工学構造物を監視するために、圧縮センシングなどのスマートセンシングアルゴリズムを使用して、カルマン街路を効率的に測定することができる。[ 3 ]

コンクリート製の冷却塔は、特に密集して建設された場合、さらに深刻な不安定性を引き起こす可能性があります。1965年、フェリーブリッジ発電所Cでは、強風時に渦 放出により3基の冷却塔が倒壊しました。

元のタコマナローズ橋の崩壊は、当初は渦放出による過度の振動が原因とされたが、実際には空力弾性フラッターによって引き起こされた。

カルマン乱気流は飛行機にとっても問題であり、特に着陸時に問題となる。[ 15 ]

ソリューション

チューンドマスダンパー(TMD)

円筒形物体(高い煙突やマストなど)の渦放出を防止し、不要な振動を軽減するために、同調マスダンパー(TMD)がよく用いられます。同調マスダンパーは、渦放出によって引き起こされる振動を打ち消すために特別に設計・調整されたマス・スプリングシステムで構成される装置です。多くの場合、スプリングまたはダンパーを介して構造物に取り付けられます。

渦放出による振動を緩和する同調質量ダンパーの有効性は、ダンパーの質量、構造物への配置、システムの調整といった要因に依存します。エンジニアは、構造力学と渦放出現象の特性を綿密に解析し、同調質量ダンパーの最適なパラメータを決定します。多くの場合、バネは質量をケーブルで吊り下げることで表現され、同じ共振周波数を持つ振り子システムを形成します。質量は、渦放出の卓越周波数と一致する固有振動数を持つように慎重に調整されます。

構造物が渦放出による振動を受けると、同調質量ダンパーは構造物と逆位相で振動します。これにより振動が抑制され、振幅が低減され、共振や構造損傷の可能性が最小限に抑えられます。

フィンとストレーキ

渦流を分散させるためのストレークが取り付けられた煙突

このような円筒形の物体の不要な振動を防ぐ別の解決策は、下流側に取り付けることができる縦方向のフィンです。このフィンが円筒の直径よりも長い場合、渦の相互作用を防ぎ、結果として渦が付着したままになります。

当然のことながら、高層ビルやマストの場合、相対的な風はあらゆる方向から吹く可能性があります。そのため、大きなネジ山のような螺旋状の突起(ストレーキ)が上部に配置されることがあります。これにより、非対称な三次元の流れが効果的に形成され、渦の交互放出が抑制されます。これは一部の自動車アンテナにも見られます。[ 16 ] [ 17 ]

その他の対策

高層ビルに対するもう一つの対策は、高さに応じて直径を変化させること、例えばテーパー状にすることである。これにより、建物全体が同じ周波数で駆動されるのを防ぐことができる。[ 18 ]

この式は、250 < Red < 200000 の 範囲で一般的に当てはまります 。 St=0.198(119.7Red) {\displaystyle {\text{St}}=0.198\left(1-{\frac {19.7}{{\text{Re}}_{d}}}\right)\ }St=fdU{\displaystyle {\text{St}}={\frac {fd}{U}}}

  • f = 渦放出周波数。
  • d = シリンダーの直径
  • U = 流速。

この無次元パラメータ St はストローハル数として知られており、1878 年に電信線の一定のハミング音または歌声を初めて調査したチェコの物理学者、ヴィンセンツ ストローハル(1850 年 - 1922 年)にちなんで名付けられました。

歴史

渦列はセオドア・フォン・カルマンにちなんで名付けられたが[ 19 ] [ 20 ]、彼は[ 21 ] 、この渦列が以前にアーヌルフ・マロック[ 22 ]アンリ・ベナール[ 23 ]によって研究されていたことを認めている。カルマンは著書『空気力学』の中でこの話を次のように語っている。[ 24 ]

[...]プラントルは博士課程の学生カール・ヒーメンツに、円筒背後の流れの剥離を観察できる水路を作るという課題を与えました。境界層理論によって計算された剥離点を実験的に検証することが目的でした。そのためには、まず定常流における円筒周囲の圧力分布を知る必要がありました。ところが、ヒーメンツは驚いたことに、水路内の流れが激しく振動していることに気付きました。彼がこのことをプラントルに報告すると、プラントルは「明らかに円筒は円形ではありません」と言いました。しかし、円筒を非常に慎重に加工した後でも、流れは振動し続けました。そこでヒーメンツは、水路が対称ではない可能性があると指摘され、調整を始めました。私はこの問題を気にしていませんでしたが、毎朝研究室に来るとヒーメンツに尋ねました。「ヒーメンツさん、流れは安定していますか?」彼は非常に悲しそうに「常に振動しています」と答えました。

フォン・カルマンは自伝の中で、この発見が、幼子イエスを背負って水の中を歩いている聖クリストファースを描いたイタリアの絵画に触発されたと述べています。水中には渦が見られ、フォン・カルマンは「歴史家にとっての問題は、なぜクリストファースがイエスを背負って水の中を歩いていたのかということだったかもしれません。私にとっては、なぜ渦があるのか​​ということでした」と述べています。研究者たちは、この絵画はボローニャのサン・ドメニコ教会博物館に所蔵されている14世紀の絵画ではないかと推測しています。[ 25 ]

参照

参考文献

  1. ^ JE Cooper (2001). 「空力弾性応答」. S. Braun (編). 『振動百科事典』 . エルゼビア. pp.  87– 97. doi : 10.1006/rwvb.2001.0125 . ISBN 978-0-12-227085-7
  2. ^セオドア・フォン・カルマン、空気力学。マグロウヒル (1963): ISBN 978-0-07-067602-2ドーバー(1994年):ISBN 978-0-486-43485-8
  3. ^ a bバユンドゥル、チハン;ナムル、バリシュ (2021)。 「圧縮センシングを使用したフォン・カルマン渦の効率的なセンシング」。コンピューターと流体226 104975.arXiv : 2005.08325土井10.1016/j.compfluid.2021.104975S2CID 234828962 
  4. ^ Amalia, E.; Moelyadi, MA; Ihsan, M. (2018). 「乱流モデルと数値時間ステップが円柱のフォン・カルマン流挙動と抗力精度に与える影響」 . Journal of Physics: Conference Series . 1005 (1) 012012. Bibcode : 2018JPhCS1005a2012A . doi : 10.1088/1742-6596/1005/1/012012 . S2CID 126372504 . 
  5. ^ Albarède, P., & Provansal, M. 準周期的円筒後流とギンツブルグ・ランダウモデル. Journal of Fluid Mechanics, 291, 191-222, 1995.
  6. ^ Farazande, S. および Bayindir, C., フォン・カルマン渦と複素ギンツブルグ・ランダウ方程式のソリトンとの相互作用。国際工学応用数学会議 (ICAME) 2021年9月1日~3日 - トルコ、バルケシル
  7. ^ Monkewitz, PA, Williamson, CHK, Miller, GD, 円筒後流におけるカルマン渦の位相ダイナミクス. 流体物理学, 8, 1, 1996.
  8. ^ Kleine, Vitor G.; Hanifi, Ardeshir; Henningson, Dan S. (2022). 「双複素数を用いた2次元ポテンシャルフローの安定性」 . Proc. R. Soc. A. 478 ( 20220165). arXiv : 2203.05857 . Bibcode : 2022RSPSA.47820165K . doi : 10.1098/ rspa.2022.0165 . PMC 9185835. PMID 35702595 .  
  9. ^ Jackson, CP (1987). 「様々な形状の物体を通過する流れにおける渦放出の開始に関する有限要素法による研究」. Journal of Fluid Mechanics . 182 : 23–45 . Bibcode : 1987JFM...182...23J . doi : 10.1017/S0022112087002234 . S2CID 123071463 . 
  10. ^ Provansal, M.; Mathis, C.; Boyer, L. (1987). 「ベナール・フォン・カルマン不安定性:過渡状態と強制状態」 . Journal of Fluid Mechanics . 182 : 1– 22. doi : 10.1017/S002211208700223 .
  11. ^ Barkley, D.; Henderson, RD (1996). 「円筒後流の3次元フロケ安定性解析」. Journal of Fluid Mechanics . 322 : 215–241 . Bibcode : 1996JFM...322..215B . doi : 10.1017/S0022112096002777 . S2CID 53610776 . 
  12. ^ 「Rapid Response - LANCE - Terra/MODIS 2010/226 14:55 UTC」 Rapidfire.sci.gsfc.nasa.gov . 2013年12月20日閲覧
  13. ^ Etling, D. (1990-03-01). 「大島嶼からのメソスケール渦放出:回転成層流の実験室実験との比較」.気象・大気物理学. 43 (1): 145– 151. Bibcode : 1990MAP....43..145E . doi : 10.1007/BF01028117 . ISSN 1436-5065 . S2CID 122276209 .  
  14. ^アーウィン、ピーター・A. (2010年9月). 「渦と高層ビル:共鳴のレシピ」. Physics Today . 63 (9). アメリカ物理学会誌: 68–69 . Bibcode : 2010PhT....63i..68I . doi : 10.1063/1.3490510 . ISSN 0031-9228 . 
  15. ^ 「空港開館式典延期」2016年7月26日時点のオリジナルよりアーカイブ2016年10月18日閲覧。
  16. ^ Ahmed, Syed R. (2007-09-19), Crocker, Malcolm J. (ed.), "Aerodynamic Sound Sources in Vehicles—Prediction and Control" , Handbook of Noise and Vibration Control , Hoboken, NJ, USA: John Wiley & Sons, Inc., p. 1076, doi : 10.1002/9780470209707.ch87 , ISBN 978-0-470-20970-7、 2023年8月11日取得{{citation}}: CS1 maint: work parameter with ISBN (link)
  17. ^ 「風切り音」自動車車体の騒音と振動制御、英国チチェスター:John Wiley&Sons、Ltd、p. 362、2018年10月12日、doi10.1002/9781119515500.ch6ISBN 978-1-119-51550-0, S2CID  240055532 , 2023年8月11日取得{{citation}}: CS1 maint: work parameter with ISBN (link)
  18. ^ Crocker, Malcolm J. (2007-09-19), Crocker, Malcolm J. (ed.), "Vibration Response of Structures to Fluid Flow and Wind" , Handbook of Noise and Vibration Control , Hoboken, NJ, USA: John Wiley & Sons, Inc., pp.  1375– 1392, doi : 10.1002/9780470209707.ch116 , ISBN 978-0-470-20970-7、 2023年8月11日取得{{citation}}: CS1 maint: work parameter with ISBN (link)
  19. ^ T. フォン・カルマン: Nachr.ゲス。ヴィッセンシャフト。ゲッティンゲンの数学。物理学。 Klasse、509 ~ 517 ページ (1911 年) および 547 ~ 556 ページ (1912 年)。
  20. ^ T. フォン カルマン: および H. Rubach、1912: Phys. Z.、vol. 13、49–59ページ。
  21. ^ T. Kármán, 1954.「空気力学:歴史的発展の観点から見た選択されたトピックス」(コーネル大学出版局、イサカ)、pp.68-69。
  22. ^ A. Mallock, 1907: 空気抵抗についてProc. Royal Soc. , A79, pp. 262–265.
  23. ^ H. ベナール、1908 年: Comptes Rendus de l'Académie des Sciences (パリ)、vol. 147、839–842、970–972ページ。
  24. ^フォン・カルマン、T. (1954)。空気力学 (Vol. 203)。コロンバス: マグロウヒル。
  25. ^溝田健人;ズドラフコビッチ、ミッキー。グロウ、カイウ。アルフレッド、リーダー(2000 年 3 月)。「文化の中の科学」自然404 (6775): 226.土井: 10.1038/35005158ISSN 1476-4687