凍結ロータ法 — ボリュート付き遠心ポンプへの適用
遠心ポンプのモデル構成
遠心ポンプの全体モデルはどう組むんですか?
典型的な構成は以下だ。
- 吸込管: 静止域、非回転
- インペラ: 回転域、MRFまたはSliding Mesh
- ボリュート: 静止域、非軸対称
- インペラ-ボリュート界面: Frozen Rotor(定常)or Sliding Mesh(非定常)
ボリュートが非軸対称だからMixing Planeは適用できない。これがFrozen Rotorが重宝される理由だ。
ボリュートのメッシュはどうやって作りますか?
ボリュートは断面形状がスパイラル状に変化するから、TurboGridでは作れない。Ansys MeshingやFluent Meshingの非構造メッシュ、またはSTAR-CCM+の自動メッシュを使う。断面形状をスイープ方向に並べてヘキサドミナントメッシュにするのが品質面で有利だ。
ボリュート舌部の処理
カットオフ(舌部)付近は何が難しいんですか?
舌部はインペラ出口流と再循環流が衝突する領域で、圧力勾配が急峻だ。メッシュを特に細かくする必要がある。また、Frozen Rotorでは翼と舌部の相対位置で流れ場が大きく変わるため、圧力脈動評価にはSliding Meshが不可欠だ。
性能マップの作成
Frozen Rotorで遠心ポンプのH-Q曲線は作れますか?
作れるが、各流量点で複数位相の平均をとることを推奨する。手順は以下だ。
1. 設計流量でFrozen Rotor計算を収束
2. 翼を10度間隔で3~5位相回転させて再計算
3. 各位相の揚程を算術平均
4. 流量を変えて2-3を繰り返し
結構手間ですね…
だから実務ではFrozen Rotorで設計点付近の1~2点を素早く評価し、最終的な性能マップはSliding Meshの非定常計算で確定させるケースが多い。
レイノルズの実験(1883年)——乱流発見の瞬間
オズボーン・レイノルズは、管内の水にインクを流す実験で「層流から乱流への遷移」を発見しました。流速を上げていくと、インクの線がある瞬間にグチャグチャに乱れる。この劇的な瞬間を、レイノルズは数学的に $Re = \rho uD/\mu$ という無次元数で表現した。100年以上経った今も、CFDエンジニアが最初に確認するのはこのレイノルズ数です。
実務者のための直感的理解
この解析分野のイメージ
CFDって、要は「デジタル風洞」です。自動車メーカーが巨大な風洞実験設備に何億円もかけるところを、PCの中で再現できる。でも1つ注意——風洞実験なら「風を当てれば結果が出る」けど、CFDでは「メッシュの品質」と「乱流モデルの選択」という見えない品質要因がある。ここを手抜きすると、きれいなコンター図が出ても中身はデタラメ…なんてことになりかねません。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め(計算領域)、水の入り口と出口を設計し(境界条件)、ポンプの強さを設定する(流量条件)。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか? 意外と直感的ではありませんか?
初心者が陥りやすい落とし穴
「y+って何ですか?」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない!」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た!」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか(低速)、全開にするか(高速)。でもCFDではもう一つ——「どのくらい暴れた水を出すか」(乱流強度)も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね? CFDでも入口条件のミスは下流全体に波及します。
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
Project NovaSolver — CAE実務の課題に向き合う研究開発
「凍結ロータ法をもっと効率的に解析できないか?」——私たちは実務者の声に耳を傾け、既存ワークフローの改善を目指す次世代CAEプロジェクトに取り組んでいます。具体的な機能はまだ公開前ですが、開発の進捗をお届けします。
進捗通知を受け取る →