バルブ流れ解析

バルブのCFD解析は、流量係数（Cv/Kv値）の予測、キャビテーション特性の評価、バルブ後方の圧力回復と騒音源の特定を目的とする。バタフライバルブ、ボールバルブ、ゲートバルブ、グローブバルブなど、バルブタイプごとに固有の流れ現象がある。

バルブの流量特性は流量係数 $C_v$（米国式）または $K_v$（欧州式）で表す。

$Q$ は流量 [US GPM]、$SG$ は比重（水=1）、$\Delta p$ は差圧 [psi] だ。SI単位系の $K_v$ との関係は:

$K_v$ の定義: $\Delta p = 1$ bar、水温15℃で $K_v$ [m³/h] の水が流れる。

CFDで得られた入口-出口間の圧力損失 $\Delta p$ と流量 $Q$ から上式で逆算する。ISA/IEC 60534規格に準拠した評価方法がある。

キャビテーション指数 $\sigma$ で評価する。

$p_1$ が上流圧力、$p_2$ が下流圧力、$p_v$ が飽和蒸気圧だ。$\sigma$ が臨界値 $\sigma_i$（Incipient Cavitation Index）を下回るとキャビテーションが始まる。

その通り。グローブバルブは流路が屈曲するため圧力回復が小さく、Vena Contracta（縮流部）の最低圧力と下流圧力の差が小さい。結果として、同じ$\Delta p$でもキャビテーションが起きにくい（$\sigma_i$が大きい）。

圧力回復係数 $F_L$ はIEC 60534で定義されるバルブ固有の係数だ。

そう。実験では下流の圧力タップで測定するから、Vena Contractaの位置を正確に特定するのが難しい。CFDなら流路内の最低圧力点を直接可視化できる。

バルブ内部は複雑な3D形状で、流路断面が急激に変化する。メッシュ品質が結果に大きく影響する。

そう。10%開度のバタフライバルブでは、ディスクと管壁の隙間が数mmしかない。この隙間に最低5層のセルを確保する必要がある。Inflation Layer（プリズム層）で対応する。

Cv値の算出には、一定流量でのΔpを求める方が精度が高い。Mass Flow Inlet + Pressure Outletの組み合わせが推奨だ。逆にΔpを固定して流量を求めるなら、Pressure Inlet + Pressure Outletを使う。

バルブ流れでは剥離、再付着、強い曲率効果があるため、SST k-omega が最も信頼性が高い。

Schnerr-Sauer模型またはZwart-Gerber-Belamri模型が広く使われる。VOF（Volume of Fluid）法と組み合わせて、気泡の生成（蒸発）と消滅（凝縮）を計算する。

Fluentのデフォルト値（Zwart: $C_{prod}=50$, $C_{dest}=0.01$, $R_B=10^{-6}$ m）は多くの場合で妥当な結果を出す。ただし運転圧力が非常に高い場合や、特殊な流体の場合はキャリブレーションが必要だ。

バタフライバルブの開度0°〜90°に対するCv値をCFDで算出し、メーカーカタログと比較する。

手順:

そう。バタフライバルブは全開でもディスクの厚さ分の圧損が残る。全開時の損失係数Kは0.2〜0.5程度だ。

安全弁（Relief Valve）の吹き出し圧力、排出量、反力をCFDで評価する。圧縮性ガスの場合、超音速流れ（チョーク流れ）が発生する。

チョーク流量は次式で計算できる:

そう。安全弁のAPIやASME認証にはCFDベースの$C_d$予測が認められつつある。ただし、CFDの結果は実験データとの比較による検証（V&V）が必須だ。

IEC 60534-8-3に基づく空力騒音の予測手法がある。CFDで得られるバルブ下流の乱流エネルギーと散逸率から、騒音のパワーレベルを概算できる。

より正確にはLES + FW-Hで音圧レベルを直接計算するが、RANS + Broadband Noise Source Modelでスクリーニングしてから、問題のあるケースだけLESで精査するのが実用的だ。

Simcenter FloEFDはCADソフトウェア（Solid Edge, NX, Creo, CATIA）に直接埋め込まれたCFDツールだ。バルブの3D CADモデルから直接CFD解析ができるので、バルブメーカーの設計者に人気がある。メッシュは直交格子ベースで自動生成される。

Ansys Workbenchでの自動化手順:

STAR-CCM+ではDesign Managerで同様の自動化が可能だ。Morpherを使えば、メッシュの再生成なしにディスク角度を変更できる場合もある。

Morphing（メッシュ変形）は小さな形状変更には有効だが、大きな開度変化（例: 10°→90°）ではメッシュ品質が劣化するので、再生成が必要になる場合が多い。

バルブの振動や疲労寿命を評価するために、CFDで求めた流体力を構造解析に渡すFSI（Fluid-Structure Interaction）が重要になっている。

特に問題となるケース:

そう。Ansys System Couplingを使えば、Fluent（流体）とMechanical（構造）を時間ステップごとに連成させて、弁体の振動を直接シミュレーションできる。

スラリーバルブ（固液二相）やフラッシングバルブ（気液二相）のCFD解析が活発だ。

Finnie modelやOka modelがFluent/STAR-CCM+に実装されている。バルブメーカーはCFDエロージョン予測をトリム材質の選定に活用している。

グローブバルブのトリム（絞り部品）を3Dプリンティングで製造し、CFD最適化された複雑な多段減圧構造を実現する研究が進んでいる。

1段の絞りで大きな圧力降下を起こすとキャビテーションと騒音が発生する。これを複数段の小さな絞りに分割して、各段の圧力降下を飽和蒸気圧以上に保つ設計だ。Fisherの WhisperFlo、MetsのMaX-Trimなどの商品がある。CFDで各段の圧力分布を最適化する。

プラントのバルブにセンサー（圧力、振動、アコースティックエミッション）を設置し、CFDモデルと照合してバルブの状態を推定するデジタルツインが実用化されつつある。

Valmet、Emersonなどの大手バルブメーカーがこの方向で製品開発を進めている。

チェックポイント:

そう。下流のRecovery Zone（圧力回復域）では位置によって圧力が大きく変わるから、評価位置を揃えることが非常に重要だ。

対策:

原因: バルブ後方のKarman渦列は本質的に非定常現象だ。定常RANS計算では収束しない場合がある。

バルブのアクチュエータ（駆動装置）の選定にはディスクに作用する流体トルクの予測が必要だ。

注意点:

そう。オフセット量をCFDのパラメトリック変数にして、最適なオフセット設計を探索することもある。

対策:

その通り。Total PressureとStatic Pressureを取り違えるとマッハ数が全く異なる値になる。出口にBack Pressureとして設定するのはStatic Pressureだ。