3次元翼周りの計算-fluidlab.1.k.u.pptVIP

2003/02/19（水） プロペラまわり流れの数値計算法に関する研究 研究背景 研究目的 従来の研究との比較 プロペラ周り流れへのCFDの適用例として以下が有る。 宇都ら（１９９２） 舩野（２００１） しかし、構造格子や部分的構造格子を用いているため、 格子生成が困難であった。 そこで、本研究では、複雑形状まわりの格子生成が 容易である完全四面体格子を採用した。 さらに、解像度が必要なところに格子を集めるために、 「解適合格子法」を導入した。 解適合格子とは 計算結果を利用して、解像度向上の必要な場所に充分細かい格子が配置されるように、分割し直す格子。 翼端渦などは、細かく分割すべき場所を事前に予測しにくいので、解適合格子は特に有効。 当発表の流れ 検証実験 キャビテーションタンネルにて プロペラ単独性能試験 キャビテーションの観察 プロペラ模型の主要目 　 実験の条件 一様流中でプロペラを回転させ、圧力一定の下、流速を変化させて推力?トルクを計測 プロペラ回転数n=25[rps] プロペラレイノルズ数Rn= =1.1×10６ タンネル減圧により、 σn= =1.0まで7通り それぞれのにσnについて、流速を系統的に変化 流速は、前進係数Jが0.05ずつ増すように調節 実験結果(プロペラ単独性能) sn=1.5のとき キャビテーションが発生しているが、推力低下はわずかである。 実験結果のまとめ 数値計算 使用ソフトウェア 3次元CADによるモデリング 半径rにおける螺旋座標系(s,n) 座標変換 螺旋座標系　?　円筒座標系　?　直交座標系 CAD内で完成した三次元形状 計算領域と境界条件の定義 格子分割(翼面) 格子分割(全体) Fluent推奨の計算条件(計算1) 四面体、非構造格子に分割 計算結果(計算1) 誤差大 （J=0.6で推力12%、トルク47%の誤差） 圧力分布(計算1) 計算条件の改良(計算2) 圧力勾配による解適合格子を生成 計算条件を以下のように変更 計算結果(計算2) 誤差の低減 に成功 （J=0.6で 　　推力の誤差12%?3% 　トルクの誤差47%?19%） 翼端渦(J＝0.6) 計算1 計算2 まとめ① CFD計算精度の検証?向上の為の検証データを実験により取得 オフセットデータで与えられるプロペラ形状を3次元ソリッドモデルに変換する手順を明確化 まとめ② 　 Fluentの推奨する条件を用い、Gambitで格子生成、Fluentにより数値計算 推力係数 12%の誤差 トルク係数 47%の誤差 まとめ③ 　次に、以下のように改良して計算 格子の解適合分割 乱流モデルをRNGk-εに変更 運動量方程式の対流項のスキームをQUICKに変更 推力係数 誤差12%?3% トルク係数 誤差47%?19% 今後の課題 誤差を半減できたものの、トルク係数の誤差がまだ19%。格子解像度の不足が原因と考えられる。 格子解像度の不足している場所は、 翼面近傍の境界層（薄い三角柱の格子が必要） 下流側に放出された渦 解適合を解が収束するまで行う まとめです。まず、CFD計算精度の検証?向上の為の検証データを、実験により取得しました。 つぎに、オフセットデータで与えられるプロペラ形状を3次元ソリッドモデルに変換する手順を明確化しました。そして、 船舶海洋工学科　船舶コースの八巻です。プロペラ周りの流れの数値計算法に関する研究 について発表致します。 この発表の流れは、以下の様になります。 　数値計算結果の評価の基準とするための舶用プロペラまわり流れの実験、 　検証実験と同じ舶用プロペラのまわり流れを、標準的条件と改良した条件とで数値計算、 それぞれの結果と実験値と比較し、その改良法を検討、そして、まとめです。 計算領域は、半径方向に3R,上流下流方向とも5Rずつの円筒形を取りました。Rはプロペラの半径です。流れパターンの同じ領域が軸中心に90度ずつ４つ並ぶと考えられますので、そのうちの一つだけをとって4分の1円筒を計算領域としました。また、実際はプロペラが回るんですけれども、計算机上では座標が回転するという設定にしました。この領域を格子分割します。 舶用プロペラは従来、ポテンシャル理論に基づく理論計算と模型試験により設計されてきました。しかし、粘性や非線形性といった問題点が有ります。そこで、ナビエ?ストークス方程式に基づく数値計算による設計が必要になってくるのですが、精度の問題と計算を行う際の前処理の困難さからまだCFDを用いての計算は一般的ではありません。 そこで今度の研究の目的です。 シャフトやボスを含んだプロペラ幾何形状は複雑で、CADで完全に三次元形状を再