波の回折係数K_Dとは何ですか？

回折係数K_Dは、防波堤等の構造物背後における波高を入射波高で無次元化した値です。K_D=1は入射波と同じ波高、K_D=0は完全に遮蔽された状態を表します。ソマーフェルト理論（Penney-Price解）では、半無限防波堤背後のK_Dは回折角・波長・距離の関数として解析的に求められます。

港湾設計における回折計算の重要性は？

港湾内の静穏度は船舶の安全な係留・荷役作業に直結します。防波堤開口部からの回折波は港湾奥部に伝播し、K_D=0.3以下の条件で良好な静穏度が得られるとされています。数値波動計算（BEM、FEM）の前段階として、ソマーフェルト理論による解析解は初期設計段階での防波堤配置・開口幅の概算に広く利用されています。

波エネルギーフラックスの計算式はどのようなものですか？

波エネルギーフラックスP（単位クレスト幅あたり）はP = ρg²H²T/(32π)で表されます。ここでρは海水密度（1025 kg/m³）、gは重力加速度、Hは波高、Tは周期です。波速C=gT/(2π)（深水）を用いると、P = (1/8)ρgH²×C_g（群速度）とも表せます。港湾内ではK_Dを用いてP_harbor = K_D² × P_incident となります。

このツールはCAE・数値波動解析とどう連携できますか？

本ツールで得られるK_D分布は、OpenFOAMやMIKE 21等の数値波動計算の境界条件設定や検証データとして利用できます。特に初期設計段階での防波堤開口幅Bの最適化、入射角感度解析に有効です。BEMコード（WAMIT等）との比較においても、遠方場の回折係数は良好な一致を示します。

波の回折・港湾内波浪計算機 — 無料オンライン計算機

パラメータ設定

モード

波周期 T

水深 h

開口幅 B

入射角 θ

0°=正面入射，±60°=斜め入射

入射波高 H₀

計算結果

—

K_D（港湾中心）

—

H_harbor [m]

—

遮蔽率 [%]

—

波エネルギー [kW/m]

—

波形勾配 H/L

可視化

K_D vs 角度（固定距離）

中心線波高 vs 距離

理論・主要公式

半無限防波堤による回折（Penney-Price簡易解）：

$$K_D(\xi) = \frac{1}{2}\left|F(\xi_+) + F(\xi_-)\right|$$

フレネル積分：$F(\xi) = \frac{1+i}{2}\int_\xi^\infty e^{-i\pi t^2/2}dt$

$\xi_\pm = \sqrt{\frac{2kr}{\pi}}\sin\frac{\theta\pm\alpha}{2}$，　$k=2\pi/L$（波数），$r$=距離，$\alpha$=防波堤方向角

波エネルギーフラックス：$P = \dfrac{\rho g^2 H^2 T}{32\pi}$ [W/m]

波の回折・港湾内波浪計算とは

🙋

防波堤の後ろって、波が入ってくるのに、どうやって波の高さを計算するんですか？シミュレーターの「入射波高 H₀」を変えると、中の波が全部大きくなるだけに見えます。

🎓

大まかに言うと、防波堤の「影」の部分に回り込む波の高さを計算するんだ。シミュレーターで「開口幅 B」を小さくしてみて。防波堤の隙間が狭まると、中の波高マップが赤（波高い）から青（波高低い）に変わるのがわかるよ。これが「回折」の効果で、隙間から回り込んだ波が広がっていく様子を計算しているんだ。

🙋

え、そうなんですか！「回折係数K_D」って、この青や赤の値そのものなんですか？

🎓

その通り！K_Dは「中での波高」÷「外の波高(H₀)」を表す無次元の値だ。例えばK_D=0.5なら、外の波の半分の高さしかないということ。実務では、K_Dが0.3以下のエリアを「静穏域」として、安全に船を停められる場所と考えることが多いね。上の「波周期 T」や「水深 h」を変えると、波長が変わるから、回折の広がり方も変わってくる。確認してみて。

🙋

なるほど！でも、こんな複雑な波の広がりを、どうやって数式一つで計算できるんですか？

🎓

良い質問だ！鍵は「フレネル積分」という数学的な関数にある。防波堤の端を波が回り込む現象は、光の回折と物理的に同じなんだ。だから、光学で使われるソマーフェルトの理論を波に応用したPenney-Priceの簡易解が使える。シミュレーターは、君がパラメータを変えるたびに、この積分をリアルタイムで解いて、中の波高分布を描いているんだよ。

よくある質問

理論上、回折係数K_Dは1を超えることはありません。K_Dは防波堤背後の波高を入射波高で除した比であり、エネルギー保存則により最大値は1です。ただし、複数の波が干渉して局所的に波高が高くなる場合、見かけ上1を超えるように見えることがありますが、これは計算上の誤差や表示スケールの影響です。

開口幅は防波堤の設計図面から実寸をメートル単位で入力してください。波長Lは、水深と周期から分散関係式L = (gT²/(2π))·tanh(2πd/L)で計算されますが、本ツールでは周期と水深から自動算出します。浅海では波長が短くなるため、適切な水深を設定することが重要です。

色の急変は、波の回折パターンにおける干渉縞（回折ローブ）に対応します。防波堤開口部から放射状に広がる波面が互いに強め合う（明るい領域）と弱め合う（暗い領域）を繰り返すことで生じます。特に開口部近傍や防波堤先端付近で顕著に現れ、物理的に正しい現象です。

波エネルギーフラックスの分布から、係留船舶の動揺が大きいエリアや、護岸への波圧集中箇所を特定できます。例えば、エネルギーが集中する領域では航路設計の見直しや消波工の設置を検討する材料となります。また、静穏域の確保が必要な岸壁の配置計画にも活用できます。

実世界での応用

港湾・漁港の静穏度設計：船舶が安全に停泊し、荷役作業を行うためには、港内の波高を一定以下（目安としてK_D<0.3）に抑える必要があります。本理論を用いて防波堤の配置や開口幅を初期設計し、港のどのエリアが十分に静穏かを事前に評価します。

数値波動解析（CAE）の検証・境界条件設定：OpenFOAMやMIKE 21 SWなどで港湾全体の詳細な波動解析を行う際、計算領域の境界や防波堤近傍の結果が物理的に正しいかどうかを、この解析解と比較して検証します。また、大規模計算の境界条件として回折波の情報を与えるのにも利用されます。

海洋構造物の波浪力評価：防波堤や桟橋の背後に設置される小型船舶用係留施設や作業台などは、回折によって低減された波高を設計波高として設定します。これにより、過大な設計を防ぎ、建設コストを最適化できます。

海岸線の侵食対策計画：防波堤の建設により、背後に波の「影」ができ、海岸線への波浪エネルギー供給が不均一になります。この回折計算は、防波堤建設が下流側の海岸線に与える影響（影領域での堆積、その外側での侵食）を予測するための第一歩として用いられます。

よくある誤解と注意点

このツールを使い始めるとき、いくつか勘違いしやすいポイントがあるよ。まず、「回折係数K_Dは波高の比であって、エネルギーそのものではない」ということ。K_D=0.5なら波高は半分だけど、波のエネルギー（波力や越波量の目安）は波高の2乗に比例するから、実は0.25倍にまで減っているんだ。設計では「波高が小さいからOK」と油断せず、施設に働く力を適切に見積もろう。

次に、水深hと波周期Tは独立したパラメータではないということ。ツールでは別々に変えられるけど、現実の波は分散関係 $L = (gT^2/2\pi) \tanh(2\pi h/L)$ で水深と波長Lが結びついている。例えば、周期10秒の波が水深5mの浅い海域に来ると、波長は約50m（深海なら約150m）と短くなる。この「浅水変形」が起きると回折の広がり方も変わるから、実データを使う時は水深と周期の組み合わせに注意してね。

最後に、この計算は「定常な単一波」が前提だという根本的な限界を理解しておこう。実際の海は不規則波で、方向スペクトルを持っている。だから、ツールで計算したK_D=0.3の静穏域でも、別の方向から来る波成分が入り込む可能性はある。実務では、この結果を基にしつつ、不規則波や方向分布を考慮したより高度な数値シミュレーションで確認するのが流れだよ。

波の回折・港湾内波浪計算

波の回折・港湾内波浪計算とは

よくある質問

実世界での応用

よくある誤解と注意点

関連ツール