Froude数 — CAE用語解説
Froude数
先生、Froude数って船の設計で大事って聞いたんですけど、どういう数なんですか？
理論と物理
Froude数の定義と物理的意味
Froude数って、教科書には「慣性力と重力の比の平方根」と書いてありますが、具体的にどういう場面で重要なんですか？
自由表面を持つ流れ、例えば船の造波抵抗やダムの放流、河川流で本質的な役割を果たします。具体的には、長さスケール
「1に近づく」と言われましたが、具体的にFr=0.9とFr=1.1では何が違うんですか？
大きな違いは、下流の影響が上流に伝わるかどうかです。Fr=0.9の常流では、水面の小さな擾乱（例えば石を投げる）が上流にも下流にも波として伝播します。一方、Fr=1.1の射流では、その擾乱は下流側にしか伝わりません。これは情報伝播の速度が流れの速度を下回るためで、衝撃波（水力ジャンプ）が発生する条件になります。
船の抵抗でよく聞く「造波抵抗」とFroude数の関係は？
船体が作る波長
数値解法と実装
CFDにおけるFroude数の取り扱い
CFDで自由表面流をシミュレーションする時、Froude数をどう設定すればいいんですか？入力パラメータとして直接指定するのでしょうか。
ソルバーによって異なりますが、多くの場合、重力加速度
を個別に入力し、ソルバー内部でFroude数を計算します。例えば、Ansys FluentのVOF（Volume of Fluid）法でダムブレークを解く場合、「Operating Conditions」で重力ベクトルを設定し、入口境界条件で速度を定義します。Frは結果として得られる無次元数です。
モデル縮小（スケールモデル）を使った実験では、Froude数相似則が使われると聞きますが、CFDではその必要はないんですよね？
CFDでは実スケールで計算できるのが利点ですが、検証のためにはFroude数相似則が重要です。例えば、実機の1/25スケール模型で実験データを取得した場合、CFDでも同じ縮尺モデルと、相似則から導かれた速度（
非定常計算で、初期状態からFrが時間とともに変化するような現象（例えば、加速する船）をシミュレーションする時、何に気をつけるべきですか？
重要なのは、時間刻み
実践ガイド
自由表面流解析のワークフロー
船体周りの流れをVOF法で解析する具体的な手順で、Froude数を意識するポイントはどこですか？
主に3点です。1. **前処理**: 計算領域は、船体後方に少なくとも船長の3〜5倍、水面上下にも十分な領域を確保します。Frが高いほど、波が長く伸びるためです。2. **初期条件**: 静水中の船体から計算を始める「加速法」が一般的で、目標のFrに達するまでの過渡状態も含めて評価します。3. **後処理**: Frに応じて、造波抵抗の寄与を評価するため、船体表面の圧力積分と、運動量理論による抵抗算出を両方行い、クロスチェックします。
メッシュの切り方は？特に水面付近は細かくする必要がありますか？
必須です。自由表面が存在する予想領域には、少なくとも10〜20層の細かいメッシュを、界面キャプチャの精度を高めるために貼ります。Ansys Meshingの「Inflation」や、Siemens Star-CCM+の「Prism Layer Mesher」で作成します。メッシュサイズの目安は、予想される波長
計算が発散しやすい条件、例えばFrが1に近いような場合、安定化するためのコツはありますか？
はい。まず、疑似時間刻み法を用いる定常計算ではなく、非定常計算から始めることを勧めます。次に、流速と圧力の連成解法には「Coupled Algorithm」よりも「PISO」スキームの方が安定しやすい場合があります。また、界面の急峻さを緩和するために、界面シャープニングスキームの係数をデフォルト（通常1.0）から0.5程度に下げることも有効です。OpenFOAMの`interFoam`ソルバーでは、`nAlphaCorr`（界面再構成の反復回数）を増やすことで対応できます。
ソフトウェア比較
主要CAEソフトにおける自由表面流解析機能
Ansys Fluent、Star-CCM+、OpenFOAMで、Froude数が関わる自由表面流解析のアプローチに違いはありますか？
コアとなるVOF法の考え方は同じですが、実装とユーザビリティに差があります。**Ansys Fluent**は「Multiphase Model」内で明示的VOFと暗黙的VOFを選択でき、造波抵抗のモニタリング用に「Froude Number」を式で定義するカスタム関数が用意されています。**Siemens Star-CCM+**は、専用の「VOF Waves」物理モデルがあり、規則波や irregular wave を入射境界条件として直接設定でき、Fr相似則を使ったスケーリング機能が組み込まれています。**OpenFOAM**の`interFoam`はコードレベルで制御可能で、界面補正アルゴリズム（MULES）のパラメータを細かく調整できますが、初期設定はやや難易度が高いです。
専用の船舶解析ソフト（例えば、SHIPFLOWやFine/Marine）は、汎用CFDソフトと比べて何が優れているんですか？
最大の利点は、船体抵抗の成分分解が自動化されている点です。SHIPFLOW（Flowtech）は、粘性抵抗（RANS法）と造波抵抗（ポテンシャル流理論）を別々に計算する「分割法」を採用しており、Frを変数とした造波抵抗係数
構造解析との連成（流体構造連成：FSI）で、Frが影響する事例はありますか？
高速艇やヨットの帆のフラッター解析が典型例です。Ansys WorkbenchやSIMULIA（Abaqus/CFD）の連成プラットフォームを使用します。ここで重要なのは、流体側の荷重（圧力）の変動周波数が、構造物の固有振動数と一致しないかどうかです。Frが変化すると、剥離渦の周波数（ストローハル数に関連）や造波のパターンが変わり、加振力のスペクトルが変化します。特にFrが0.4〜0.6の範囲では、船首波と船尾波が干渉して強い波励起振動を起こす可能性があり、連成解析で評価が必要です。
トラブルシューティング
自由表面流解析のよくあるエラーと対策
計算中に自由表面がどんどん拡散して、はっきりした界面にならないのですが、原因は何ですか？
主に3つの原因が考えられます。1. **メッシュが粗すぎる**: 先述したように、波長に対して十分細かいメッシュが必要です。2. **界面シャープニングが不十分**: Fluentの「Sharpening/Damping Factors」、Star-CCM+の「Interface Compression」の値を調整します（通常は1.0を基準に増減）。3. **時間刻みが大きすぎる**: 界面の移動を正確に追跡するには、CFL数（Courant数）を1以下、理想的には0.5以下に保つ必要があります。特にFrが大きい（流速が速い）計算では、自動的に時間刻みを調整する機能を使うと良いでしょう。
造波抵抗の値が、実験値や経験式と比べて明らかに大きすぎる（または小さすぎる）場合、どこを疑えばいいですか？
まず、**計算領域の出口境界**を疑います。出口が船体に近すぎたり、静圧条件に設定していると、波が出口で反射して船体に戻り、抵抗が過大評価されます。出口には「波を吸収する」ためのダンピングゾーン（Wave Absorption）を設けるか、Orlanski型などの放射境界条件を使用します。次に、**乱流モデル**です。船体近傍の境界層を適切に解像できていないと、粘性抵抗の評価が狂い、造波抵抗との分離が不正確になります。最後に、**収束判定**。非定常計算では、少なくとも10波長分以上の計算時間を確保し、抵抗値の時間平均を取る必要があります。
Frをパラメータとして変化させた一連の計算で、ある一点（例えばFr=0.45）だけ計算が発散してしまいます。なぜですか？
それは、そのFr数で何らかの**共鳴現象**や**不安定な流れモード**が発生している可能性が高いです。例えば、船体から発生する波の干渉パターンが劇的に変化するFr数（「ホローの干渉条件」）付近では、流れ場が非常に敏感になります。対策としては、1) 発散するFrの計算だけ、より小さい時間刻みから始める、2) 初期条件を、収束した近傍のFr（例えばFr=0.4）の結果から継承する（restart計算）、3) より強力な数値拡散（高次精度スキームから1次風上差分へ一時的に変更）を導入して不安定性を抑える、などが挙げられます。このような特異点の存在自体が、物理的に重要な知見になることもあります。
関連トピック
なった
詳しく
報告