流体の圧力と浮力 — パスカルの原理・ベルヌーイからCFD解析まで
1. 流体圧力とCAE — なぜエンジニアは圧力を気にするのか
水、空気、油などの流体はあらゆる工学システムに関わっている。圧力はその流体が構造物に与える力の根源であり、設計の可否を左右する。配管の破裂、ダムの決壊、航空機の与圧破壊、深海探査機の座屈——これらはすべて圧力の見積もりミスや想定外の圧力負荷に起因する事故だ。
CAEエンジニアにとって「圧力」は単なる数値ではなく、構造物の境界条件として直接入力される荷重である。CFD（計算流体力学）では圧力場を求めること自体が解析の目的であり、流体-構造連成解析（FSI）では流体側の圧力を構造側に転送してFEM解析を行う。圧力の物理を深く理解することは、解析モデルの境界条件設定を正確に行う基礎力となる。
2. 圧力の定義と単位
圧力は単位面積あたりに作用する力として定義される。
$$p = \frac{F}{A} \quad \text{[Pa = N/m}^2\text{]}$$圧力はスカラー量であり、向きを持たない。ただし圧力が面に作用するとき、その力は面の法線方向（外向き）に働く。
絶対圧・ゲージ圧・真空圧
圧力には3種類の表現方法があり、CAEの境界条件設定で混同しやすい。
- 絶対圧（Absolute Pressure）：真空（$p = 0$）を基準とした圧力。熱力学計算で使う。
- ゲージ圧（Gauge Pressure）：大気圧を基準（$p_{gauge} = 0$）とした圧力。圧力計が示す値。
$p_{abs} = p_{gauge} + p_{atm}$ - 真空圧（Vacuum Pressure）：大気圧より低い場合の差分。$p_{vacuum} = p_{atm} - p_{abs}$
標準大気圧は $p_{atm} = 101.325\,\text{kPa} = 1\,\text{atm} = 1013.25\,\text{hPa}$
身近な圧力スケール
| 対象 | 圧力（絶対圧） | 備考 |
|---|---|---|
| 標準大気圧 | 101.3 kPa | 海面高度での大気の重さ |
| 自動車タイヤ（ゲージ圧） | ~230 kPa | 約2.3気圧（ゲージ） |
| 血圧（収縮期） | ~16 kPa above atm | 120 mmHg |
| 水深100 m | ~1.1 MPa | 大気圧 + $\rho g h$ |
| 深海10,000 m（マリアナ海溝） | ~100 MPa | 約1000気圧 |
| 原子炉冷却水 | ~15 MPa | 水が300℃でも沸騰しない |
| ガスタービン燃焼器 | 1〜4 MPa | 高圧下で燃焼効率向上 |
3. 流体静力学（静水圧）
静止した流体中の圧力分布を記述するのが流体静力学の基本式だ。
$$p = p_0 + \rho g h$$ここで $p_0$ は液面の圧力（通常は大気圧）、$\rho$ は流体の密度 [kg/m³]、$g$ は重力加速度 [m/s²]、$h$ は液面からの深さ [m]。
重要な性質：圧力は水平面内では均一で、深さ（高さ）のみに依存する。容器の形状（広いか細いかなど）によらない。これを「パスカルのパラドックス」の基礎と呼ぶこともある。
ダムの水圧分布と合力
高さ $H$ の鉛直面（ダム壁面など）に作用する静水圧の合力 $F$ は、圧力分布を積分して求める。
$$F = \int_0^H \rho g h \cdot w \, dh = \frac{1}{2}\rho g H^2 w$$ここで $w$ はダム幅。合力の作用点（圧力中心）は底面から $H/3$ の高さにある（三角形分布の重心）。
*DLOAD, TYPE=HYDRO で水面位置と流体密度を指定するだけで、深さに比例した圧力を自動計算して各面要素に分配してくれる。
圧力水頭・速度水頭・位置水頭
静水圧の式を圧力水頭（m）として表すこともある。
$$\frac{p}{\rho g} + z = \text{const（静止流体）}$$これはベルヌーイの定理（流速ゼロの場合）の特殊形だ。
4. パスカルの原理と油圧系統
パスカルの原理：密閉された静止流体中のある点に加えた圧力変化は、流体のすべての部分に等しく伝わる。
この原理が油圧システムの基礎だ。小さい力で大きな力を発生させる「油圧増力」の原理は以下の通り。
$$\frac{F_2}{F_1} = \frac{A_2}{A_1}$$面積 $A_1 = 1\,\text{cm}^2$ のシリンダーに $F_1 = 100\,\text{N}$ の力を加えると、面積 $A_2 = 100\,\text{cm}^2$ のシリンダーでは $F_2 = 10{,}000\,\text{N} = 10\,\text{kN}$ の力が得られる。
航空機の油圧アクチュエータ
民間航空機の飛行制御（エルロン、エレベーター、ラダー）には通常、動作圧力 20〜34 MPa の油圧システムが使われる。パイロットがジョイスティックをわずかに動かすだけで、数kNに達する舵面荷重を発生させる。現代機ではこの油圧システムを3系統以上冗長化し、単一障害でも安全を確保する。
CAE流体–構造連成解析への接続
油圧系統の解析では、配管内の非定常圧力波動（ウォーターハンマー現象）をCFDで計算し、その圧力荷重を配管FEMモデルに入力する流体–構造連成（FSI）解析が必要になる。急速にバルブを閉じると $\Delta p = \rho c \Delta v$（$c$：音速）という圧力スパイクが発生し、配管を損傷させることがある。
5. アルキメデスの原理と浮力
流体中に沈んだ物体は、その物体が排除した流体の重さに等しい上向きの力（浮力）を受ける。
$$F_B = \rho_f g V_{sub}$$ここで $\rho_f$ は流体の密度、$V_{sub}$ は物体の水没体積。
浮く条件： $\rho_{obj} < \rho_f$（平均密度が流体より小さい）
船の喫水線と排水量
船は自重 = 浮力となる深さ（喫水）まで沈む。
$$W_{ship} = \rho_{seawater} g V_{displaced}$$大型タンカー（排水量 30 万トン）の場合、$V_{displaced} = 300{,}000 / 1.025 \approx 293{,}000\,\text{m}^3$ の海水を排除する。積荷を積むと喫水が増す——このシンプルな関係が船舶工学の出発点だ。
海洋プラットフォームのFEM浮力荷重設定
洋上風力発電の浮体式プラットフォームや石油掘削リグでは、浮力が構造に直接影響する。
- 水没部材には静水圧による外圧（圧縮荷重）が作用する。
- 上向きの浮力は等価な鉛直圧力荷重としてFEMに入力する。
- 波浪による動的浮力変動は、波高・波周期を考慮した非定常CFD解析または Morison 方程式（後述）で評価する。
6. ベルヌーイの定理
非粘性・非圧縮性・定常流の流線上でエネルギー保存則を適用すると、ベルヌーイの定理が得られる。
$$p + \frac{1}{2}\rho v^2 + \rho g z = \text{const}$$3項はそれぞれ、圧力エネルギー・運動エネルギー・位置エネルギー（単位体積あたり）に対応する。流速 $v$ が大きい場所では圧力 $p$ が低くなるという逆相関が直感に反するが、エネルギーの相互変換として理解できる。
forceCoeffs 関数オブジェクトで揚力係数 $C_L$ と抗力係数 $C_D$ を自動計算してくれる。乱流モデル（$k$-$\omega$ SST が一般的）の選択で精度が大きく変わるから要注意だ。
ベルヌーイの応用例
| 応用 | 原理 | 工学的意義 |
|---|---|---|
| ベンチュリ管（流量計） | 断面積縮小→流速増→圧力低下 | $\Delta p$ から流量を計算 |
| 翼の揚力（大まかな理解） | 上面高流速→低圧 | 航空機・風力発電翼の設計 |
| カルマン渦（流速測定） | 後流の圧力変動周波数 | 超音波流量計の一原理 |
| ピトー管（航空速度計） | 淀み点圧 = $p_0 + \frac{1}{2}\rho v^2$ | $v = \sqrt{2(p_0-p)/\rho}$ |
7. 表面張力・毛管現象
液体の表面では分子間力の非対称性から表面張力 $\gamma$ [N/m] が生じる。水の表面張力は 20℃ で約 0.072 N/m。
$$\Delta p = \frac{2\gamma}{r} \quad \text{（球形液滴のラプラス圧）}$$半径 $r$ が小さい気泡や液滴ほど内外の圧力差（ラプラス圧）が大きくなる。半径 1 μm の気泡では $\Delta p = 2 \times 0.072 / 10^{-6} = 144\,\text{kPa}$——大気圧をはるかに超える。
毛管現象と接触角
細管内での液面上昇高さは:
$$h = \frac{2\gamma \cos\theta}{\rho g r}$$$\theta$ は接触角（ぬれ性の指標）。$\theta < 90°$（親水面）では液面が上昇し、$\theta > 90°$（疎水面）では降下する。
マイクロ流体デバイスのCAE解析
医療診断チップ（LOC: Lab-on-a-Chip）の流路幅は 10〜100 μm。この微細スケールでは表面張力が重力や粘性力を圧倒し、流れ挙動を支配する。COMSOL Multiphysics や OpenFOAM の VOF（Volume of Fluid）法で気液界面を追跡しながら、チャンネル内の液体輸送をシミュレーションする。
8. 粘性流体と流れのCAEモデル
実在の流体は粘性を持つ。粘性係数 $\mu$ [Pa·s]（動力学的粘性係数）と動粘性係数 $\nu = \mu/\rho$ [m²/s] で流体の「粘っこさ」を表す。
Navier-Stokes方程式の圧力項
非圧縮性粘性流体の運動方程式（Navier-Stokes方程式）の連続式と運動量式：
$$\nabla \cdot \vec{v} = 0 \quad \text{（連続式、非圧縮）}$$ $$\rho \left(\frac{\partial \vec{v}}{\partial t} + (\vec{v}\cdot\nabla)\vec{v}\right) = -\nabla p + \mu \nabla^2 \vec{v} + \rho \vec{g}$$右辺第1項 $-\nabla p$ が圧力勾配力だ。流体は圧力の高い方から低い方に向かって加速される。CFDではこの圧力場と速度場を連立して解く——OpenFOAMの SIMPLE アルゴリズムがその典型的な手法。
OpenFOAMでの圧力境界条件
| 境界タイプ | 圧力境界条件 | 速度境界条件 |
|---|---|---|
| 流入口（inlet） | zeroGradient | 固定値（$v = v_{in}$） |
| 流出口（outlet） | 固定値（$p = 0$） | zeroGradient |
| 壁面（wall） | zeroGradient | no-slip（$v = 0$） |
| 対称面（symmetry） | symmetryPlane | symmetryPlane |
圧力-速度カップリングの注意点
非圧縮性流れでは圧力方程式（ポアソン方程式）が速度発散ゼロ条件を満たすように圧力を決定する。OpenFOAMの p_rgh（ゲージ圧から静水頭を引いた修正圧力）を使うと、重力の影響がある流れ（自然対流など）の設定が簡潔になる。
9. 実践：海洋風力発電モノパイル基礎の波浪荷重
洋上風力発電の基礎構造として最も普及しているのがモノパイル（単杭）方式だ。直径 6〜10 m、長さ 40〜80 m のスチールパイプを海底に打設し、風車タワーを支える。波浪荷重の計算には Morison 方程式が広く使われる。
Morison 方程式
$$F = C_M \rho \frac{\pi D^2}{4}\dot{u} + C_D \frac{1}{2}\rho D |u|u$$第1項は慣性力（$C_M$：慣性力係数、$\approx 2.0$）、第2項は抗力（$C_D$：抗力係数、$\approx 1.0$）。$u$ は水粒子水平速度、$\dot{u}$ はその時間微分（加速度）、$D$ は部材直径。
波の理論（線形波理論）から水粒子速度を計算し、Morison 式で杭全長にわたる分布荷重を求め、FEM梁モデルに入力して根本曲げモーメントや疲労荷重を評価する。
共振の回避設計
モノパイルの固有振動数は、波の卓越周期（5〜15秒）やロータブレードの回転周波数（1P: ~0.2 Hz）から離れるように設計する。固有振動数が波や回転励振に重なると、大きな動的増幅が生じ疲労寿命が激減する。
- ソフト-スティフ設計：1P と 3P（3枚翼の場合）の間に固有振動数を配置（0.25〜0.35 Hz 程度）
- FEMモード解析で確認し、パイル径・板厚・埋め込み長を最適化する