水（20°C）
石油（鉱油）
水銀
— 判定中 —
圧力波形 P(t)
弁閉鎖時間 vs ΔP
管路長 vs 波速・水撃圧
弁位置での圧力-時間波形（簡易矩形波近似）。赤点線 = 許容圧力（PN16=1.6MPa）目安。
弁閉鎖時間 tc に対する最大水撃圧 ΔP の変化。縦点線 = 現在の tc、横点線 = 臨界閉鎖時間 Tr。
管路長 L に対する臨界閉鎖時間 Tr と最大水撃圧（同一 V₀）の変化。
急閉鎖時（tc ≤ 2L/a）の最大水撃圧上昇（Joukowsky 式）：
\(\Delta P = \rho \cdot a \cdot \Delta V\)
弾性管内の圧力波速度（Halliwell 式）：
\(a = \sqrt{\frac{K/\rho}{1 + (K \cdot D)/(E \cdot e)}}\)
ここで K = 液体の体積弾性係数、E = 管材のヤング率、D = 内径、e = 管壁肉厚。
臨界（往復）時間と緩閉鎖時の水撃圧（Allievi 近似）：
\(T_r = \frac{2L}{a}, \quad \Delta P_{slow} \approx \Delta P_{Joukowsky} \cdot \frac{T_r}{t_c} \quad (t_c > T_r)\)
水撃（ウォーターハンマー）の仕組み — 会話で理解する
🙋 「水撃」って水道の蛇口を急に閉めると「ドン！」っていう音がするやつですよね。でも工場配管で問題になるのはなぜですか？
🎓 そう、あれが水撃だ。家庭の蛇口なら音がする程度だけど、工場や発電所の大口径配管では流速も大きくて、弁を急閉鎖すると何MPaもの圧力衝撃波が発生することがある。これで配管の破断・継手の損傷・ポンプの逆転が起きる。過去に水力発電所や化学プラントで大事故が起きているんだ。
🙋 Joukowsky の式 ΔP = ρaΔV、この「a」は何ですか？音速と同じですか？
🎓 水中の音速に近いけど、正確には「管内の圧力波速度」で、液体の体積弾性係数 K と管壁の弾性 E の両方が影響する。純水の a は約 1480 m/s、でも鋼管内では管壁の弾性で少し遅くなって 1200〜1400 m/s 程度になる。ゴム管や HDPE 管だとさらに遅くて 200〜500 m/s 程度になるため水撃圧が下がる。
🙋 「臨界閉鎖時間」って何ですか？これより遅く閉めれば大丈夫なんですか？
🎓 Tr = 2L/a は圧力波が管路を往復する時間だ。弁閉鎖時間 tc > Tr なら「緩閉鎖」で、反射波が戻ってくる前に部分的に弁が開いたままなので水撃圧が下がる。tc < Tr だと「急閉鎖」で Joukowsky の最大値 ΔP = ρaΔV が発生する。例えば L = 100m、a = 1400m/s なら Tr = 0.14s。これより遅く閉めれば完全水撃は避けられる。
🙋 水撃を防ぐには弁をゆっくり閉めるだけでいいんですか？他に対策はありますか？
🎓 緩閉鎖が最も基本的な対策だ。他にはサージタンク（管路に設けた開放タンクで圧力を吸収）、エアチャンバー（密閉空気タンクが圧力衝撃を緩和）、安全弁・リリーフバルブ（過圧時に逃がす）、フライホイール付きポンプ（急停止防止）なども使われる。複合的に対策することが多い。
🙋 水撃の「負圧」って何ですか？水が「引っ張られる」んですか？
🎓 正圧波（ΔP > 0）の後に負圧波が来て、管内圧力が大気圧以下に下がることがある。これが「カラム分離（column separation）」を引き起こすと、水が気化・空洞化して、その後再び衝突するときに「二次水撃」として非常に大きな衝撃が発生する。長い管路や起伏のある配管で特に問題になる危険な現象だ。
🙋 MOC（方法of Characteristics）解析ってよく聞きますが、このシミュレーターとどう違うんですか？
🎓 このツールは Joukowsky 式を使った「単純な1次元解析」だ。実際の設計ではMOC（特性曲線法）という数値解法を使って、管路の各区間での圧力・流速の時間変化を逐次計算する。分岐管・ポンプ・弁の特性曲線を組み込んで複雑な管網全体を精密にシミュレーションできる。AFT Impulse や HAMMER（Bentley）などの専用ソフトウェアがある。
よくある質問
水撃圧の「急閉鎖」と「緩閉鎖」の切り替え境界は？
臨界閉鎖時間 Tr = 2L/a です。tc ≤ Tr で急閉鎖（最大 ΔP = ρaΔV 発生）、tc > Tr で緩閉鎖（ΔP が 1/Tr×tc 倍に低減）。L=200m、a=1400m/s なら Tr=0.29s。これより遅く弁を閉めれば最大水撃圧を抑制できます。実設計では Tr の 2〜3 倍以上の閉鎖時間を目安とすることが多いです。
管材別の圧力波速度の目安は？
鋼管（E=210GPa）: 1200〜1400 m/s、鋳鉄管: 1000〜1200 m/s、コンクリート管: 1000〜1200 m/s、HDPE 管（E=1GPa）: 200〜400 m/s、ゴム管: 100〜300 m/s。弾性の低い材料ほど a が小さく、Joukowsky 式から水撃圧も下がります。HDPE や可撓管を使う水撃低減設計もあります。
カラム分離（Column Separation）とは何ですか？
負圧波が大気圧以下（絶対圧で液体の蒸気圧以下）になったとき、液体が気化して空洞（vapor cavity）が形成される現象です。その後空洞が崩壊するときに「二次水撃」として非常に大きな衝撃圧が発生します。長い管路、起伏のある管路、高速ポンプ停止時に特に危険です。MOC 解析での空洞モデルによるシミュレーションが必要です。
ポンプ停止時の水撃はどう評価する？
ポンプが急停止すると逆流防止弁が閉じることで弁閉鎖水撃と同様の現象が起きます。ΔV はポンプ定格流速と同程度になります。対策としてフライホイール（慣性力で停止を緩やかにする）、緩閉逆止弁、バイパス配管などが使われます。ポンプ保護のためにチェックバルブの動特性も MOC 解析に組み込む必要があります。
エアチャンバー（空気室）による水撃緩和の仕組みは？
管路に設置した密閉容器に空気を封入したものです。圧力波が来ると空気が圧縮されてエネルギーを吸収し、圧力ピークを低減します。また負圧時には空気が膨張して液体を補給し、カラム分離を防止します。設計では空気容積・初期圧力・位置の最適化が必要で、圧縮性ガスの状態方程式（多変過程）を考慮します。
日本の配管設計基準での水撃の扱いは？
JISB 8265（圧力容器の構造）や JEAG 4601（原子力）では過渡荷重として水撃を考慮します。上水道ではWSP（日本水道協会規格）で弁操作時間に制限が設けられています。石油プラントでは API 651 等の国際規格、発電所では ASME B31.1 が参照されます。設計では MOC 解析結果の最大圧力に安全率をかけてフランジ定格圧力（ANSI PN）を選定します。
水撃（ウォーターハンマー）解析シミュレーターとは
本シミュレーターの物理モデルは、非定常流れにおける水撃現象をJoukowskyの基礎理論に基づき解析します。管内の流速が急変する際、圧力波が発生し、その最大水撃圧 \(\Delta p\) は \(\Delta p = \rho c \Delta v\) で与えられます。ここで \(\rho\) は流体密度、\(c\) は圧力波速度、\(\Delta v\) は流速変化量です。圧力波速度 \(c\) は管の弾性を考慮し、\(c = \sqrt{\frac{K}{\rho}} \left(1 + \frac{K D}{E e}\right)^{-1/2}\) と表されます。\(K\) は流体の体積弾性率、\(D\) は管内径、\(E\) は管材のヤング率、\(e\) は管厚です。弁閉鎖時間 \(t_c\) が圧力波の往復時間 \(2L/c\)（\(L\) は管路長）より短い場合、急閉鎖と判定され最大水撃圧が直接作用します。本モデルはこれらの関係式を基に、入力パラメータからリアルタイムで最大水撃圧・圧力波速度・臨界閉鎖時間を算出し、安全設計の指標を提供します。
実世界での応用
産業での実際の使用例
研究・教育での活用
CAE解析との連携や実務での位置付け
よくある誤解と注意点
「弁閉鎖時間が臨界閉鎖時間より短ければ必ず危険」と思いがちですが、実際にはJoukowsky式は急閉鎖（瞬時閉鎖）を前提としており、弁閉鎖時間が臨界時間より短くても、閉鎖パターン（直線的か指数関数的か）や管路の弾性、流体の圧縮性によって最大水撃圧は理論値より低くなる場合があります。また、「流速が速いほど水撃圧が高くなる」という理解は正しいものの、圧力波速度（水撃波の伝播速度）は管内径や管材のヤング率、流体の体積弾性率に依存するため、同じ流速でも配管材料や支持条件によって水撃圧が大きく変わる点に注意が必要です。さらに、「シミュレーション結果が安全判定を満たせば実機も安全」と過信しないでください。本ツールは単一の直管と単純な弁閉鎖をモデル化しており、実際のプラント配管では分岐やエルボ、バルブの非線形特性、キャビテーションの影響が無視できず、安全率を適切に考慮した設計が不可欠です。