Steady, Quasi-Periodic, and Chaotic Regimes of Natural Convection Around a Heated Cylinder in a Square Cavity

技術要約：正方形キャビティ内の加熱された円柱周囲における自然対流の定常、準周期、およびカオス的レジーム

問題提起
本研究は、中心に加熱された円形シリンダー（半径 $R = 0.2L$）を含む低温の正方形エンクロージャ内における空気（$Pr = 0.71$）の非定常自然対流を調査するものである。この差熱正方形キャビティは標準的なベンチマークであるが、内部に加熱された円柱を導入することで、複雑な流れの再編成が引き起こされる。先行文献の多くは、定常解または限定的な非定常範囲に焦点を当てており、しばしば人工的な対称性の制約（例：流線関数を円柱上でゼロに固定する）を課すことで、非対称な流れモードを抑制している。その結果、この構成における非定常性の開始となる臨界レイリー数（$Ra_c$）、カオスへの具体的な経路、および流れのレジームの全シーケンスは未解決のままであった。本研究は、対称性を強制せずに非対称モードを自然に成長させることで、定常層流からカオスへの遷移をシミュレーションし、これらのギャップを埋めることを目的としている。

手法
著者らは、ボディ適合メッシュを用いた圧力ベースの有限要素法（FEM）を採用した。主な手法の特徴は以下の通りである：

• 支配方程式： 原始変数下での、ブシネスク近似を用いた非定常2次元ナビエ・ストークス方程式およびエネルギー方程式を解いた。垂直中心線上には対称条件を課していない。
• メッシュ生成： レイキャスト投影に続くハーモニック・スムージングを用いて、デカルト背景格子からボディ適合メッシュを生成した。この手法により、座標変換やインマースド境界補間を用いることなく、曲面円柱境界を解像している。
• 数値スキーム： 対流拡散項に対して2次アダムス・バッシュフォース法を用いた、陽的な分数ステップ（投影）アルゴリズムを使用した。圧力ポアソン方程式は、非構造格子を効率的に扱うために、ウォームスタート付きの共役勾配法と要素別（EbE）行列フリー積を用いて解かれた。
• 検証および格子独立性： 解法はDe Vahl Davisのベンチマーク（$Ra = 10^6$）に対して検証された。格子収束指数（GCI）調査により、円柱のヌセルト数における細かい格子の不確かさが0.84%であることを確認し、格子独立性を確認した。
• シミュレーション範囲： $10^4$ から $5 \times 10^7$ までの8つのレイリー数についてシミュレーションを行った。円柱平均ヌセルト数（$Nu_{cyl}$）の時間履歴を記録し、キャビティ内に分散した11個のプローブ点を用いて、エネルギースペクトル密度（ESD）およびポアンカレ写像を生成した。

主要な結果
本研究は、3つの異なる流れのレジームを特定し、それらの間の遷移を特徴付けている：

1. 定常レジーム ($10^4 \le Ra \le 5 \times 10^6$):

   • 流れは定常かつ垂直中心線に対して対称であり、2つの鏡像反転循環セルと、円柱から上壁へと立ち上がる単一の熱プルームによって特徴付けられる。
   • $Ra = 5 \times 10^6$ において、最初の物理的な対称性の破れの兆候が現れる。負の流線関数コアの大きさ（$|\psi_{min}|$）が正のコア（$\psi_{max}$）を約3.6%上回るが、流れは視覚的には対称のままである。
   • 熱伝達は上壁によって支配されており（全壁熱伝達量の約55%を吸収）、下壁は下部キャビティの熱成層のため、2%未満となっている。

2. 非定常性の開始 ($5 \times 10^6 < Ra_c < 7.5 \times 10^6$):

   • 非定常性の開始となる臨界レイリー数は、$5 \times 10^6$ と $7.5 \times 10^6$ の間に挟まれている。これは、$Nu_{cyl}$ の標準偏差における急激な跳躍（3桁以上）によって決定された。
   • この臨界値は、対称性を強制した従来の定常解析の研究が示唆するものよりも高く、非対称モードの抑制が遷移を人工的に遅らせていることを裏付けている。

3. 非定常レジーム ($Ra \ge 7.5 \times 10^6$):

   • リミットサイクル ($Ra = 7.5 \times 10^6$): 流れは単一周波数の周期振動（リミットサイクル）を示す。無秩序は円柱上方のプルームに由来する。
   • 準周期 ($Ra = 10^7$): 第2の非整数周波数が現れ、流れを2次元トーラス上に置く。ポアンカレ断面は閉じた曲線を示し、下部の成層プールが二次的な周波数を感知し始める。
   • カオス ($Ra = 5 \times 10^7$): 流れは広帯域のストレンジ・アトラクタへと遷移する。ポアンカレ断面は散布された点の雲を示し、エネルギースペクトル密度は連続的になる。無秩序はプルームから側壁を下り、最終的に下部の成層プールに到達し、広範な時間スケールにわたって活動的になる。
   • 平均における対称性： 瞬間のスナップショットは強く非対称であるが、時間平均場は二方向の対称性を回復する。これは、システムが鏡像状態を時間的に等頻度で訪問するという統計的な効果である。

熱伝達相関

• 単一の冪乗則 $Nu_{cyl} = 0.764 Ra^{0.226}$ が、全範囲（$10^4 \le Ra \le 5 \times 10^7$）にわたる時間平均円柱ヌセルト数に適合し、偏差は6%未満である。
• 定常枝の指数は $n \approx 0.224$ であり、非定常枝の指数は $n \approx 0.243$ である。両方の値は、古典的な層流境界層の予測（$n=0.25$）とカスタングのソフト・タービュランスのスケーリング（$n \approx 0.286$）の間に位置しており、流れが完全な乱流ではなく、変動場の影響を受けていることを示唆している。
• 下壁のヌセルト数（$Nu_{bot}$）は $Ra$ の増加とともに減少する。これは、熱成層の強化と、下部コーナーのエディ内に冷たいジェットが捕捉されることに起因する、直感に反する傾向である。

意義と主張
本論文は、人工的な対称性の制約なしに、中心加熱円柱構成における定常からカオスへの完全な遷移シーケンスを、数値的に包括的に特徴付けた最初の研究であると主張している。主な貢献は以下の通りである：

• 臨界レイリー数： 非対称モードの発達を許容することで、より正確な $Ra_c$ の範囲を確立し、対称性の制約による従来の過大評価を修正した。
• カオスへの経路： 遷移が周期倍分岐カスケードではなく、ルエル・タケンス・ニューハウス経路（リミットサイクル $\to$ 2次元トーラス $\to$ ストレンジ・アトラクタ）に従うことを確認した。
• 無秩序の空間的伝播： 無秩序が熱プルームから発生し、下方へと伝播すること、および下部の成層プールが最も高いレイリー数に至るまで比較的静穏であることを示した。
• リファレンスデータ： 内部幾何学と非定常ダイナミクスが相互作用する構成における、臨界レイリー数、レジーム遷移、およびグローバルな熱伝達相関に関する信頼性の高いベンチマークデータを提供した。

著者らは、本研究が2次元、単一の円柱半径（$R=0.2L$）、および $Pr=0.71$ に限定されていることを指摘しており、3次元、オフセンター配置、および様々なプラントル数への拡張を今後の課題として示唆している。