この記事は著者が韓国語で書いた原文の翻訳です。 原文(한국어)を読む
回転するボールはなぜ曲がるのか — マグヌス効果と循環
スピンを消さず飛行軸にそろえると外積が0になり、横力が消えます。曲がる量はスピン比Sが決めます。
1742年、ベンジャミン・ロビンスはマスケット弾がまっすぐ飛ばないのを見ていました。同じ銃、同じ火薬、同じ狙いです。ある弾は左へ、ある弾は右へ曲がりました。銃口を出た球形の弾が不規則に回っていました。回転と曲がり。二つは揃って出ます。
100年後、グスタフ・マグヌスがその横力を実験で固定しました。回転する物体は流れの中で横に押されます。その力の大きさを書く量が循環(circulation)です。
スピンのスライダーを0からバックスピンへ上げると、上側の流線が密になり、力の矢印が上を向きます。
よどみ点が下へずれるかを見ます。スピンを逆にすると矢印も反転します。
曲がる砲弾という古くからの悩み
ロビンス以前から砲手たちは知っていました。砲弾は的を外します。風のせいだけではありません。鋳造の公差のせいで砲身と砲弾のあいだに隙間(当時の用語で windage)があり、発射の瞬間に弾は砲身の壁を跳ねて予測できない回転を得ました。その回転が飛行中の軌道を曲げたのです。
問題は深刻でした。弾道学は国家の存亡がかかる計算だったからです。ニュートンは『プリンキピア』で放物線軌道を導きましたが、実際の弾は放物線を外れました。ロビンスはそのずれを実験で切り分けました。一部は空気抵抗、一部は回転でした。
核心の問いはこれです。なぜ回転が横向きの力を生むのか。
横に押す力は循環から生まれる
まず力の大きさから見ましょう。回転する円柱が速度 の一様流の中に置かれると、単位長さあたりの横向きの力は驚くほど単純です。
ここで は流体密度、 は自由流速度、 は循環です。この関係を 揚力を循環に結ぶクッタ–ジューコフスキーの定理 と呼びます。揚力は循環だけに比例します。物体の形は、どれだけ循環を作るかという点でのみ関わります。
循環は、閉じた経路に沿って速度を積分した量です。
は流速、 は経路の接線方向の微小ベクトルです。循環がゼロでなければ、物体のまわりを回る正味の旋回流があるという意味です。半径 、角速度 で回る円柱は、表面が流体を引きずってこの正味の旋回を作ります。
向きは外積で決まります。
は回転角速度ベクトルです。この式が、後でライフリングの秘密を解く鍵になります。
なので、 が2倍になると循環も2倍です。定常の横力 も2倍です。形はこの式に出てきません。半径 は二乗で入ります。同じ角速度でも球が大きいと循環は早く大きくなります。向きは です。二つのベクトルが平行なら積は0で、横力は消えます。ライフリングが後で使う整列が、この一行です。
回転が循環を作る本当の理由
循環はなぜ生じるのでしょうか。答えは粘性と滑りなし条件(no-slip condition、壁で流体速度が壁の速度に等しくなる)にあります。
回転するボールの片面は流れと同じ向きに回り、反対面は流れに逆らって回ります。流れと同じに回る側では、表面が境界層を押してやるので流れが遅れて剥離します(剥離の遅れ)。逆に回る側では流れが早く剥離します。結果として後流が片側に傾きます。
後流が下に傾けば、流体は下向きの運動量を得ます。その反作用でボールは上に押されます。これがマグヌス力の本質です。ボールが流れを片側へ曲げ、その代償に反対向きの力を受けるのです。
どれだけ曲がるかはスピン比が決める
レイノルズ数はいったん脇に置きましょう。マグヌス効果の強さを支配する無次元数はスピン比(spin ratio)です。
分子は表面速度、分母は自由流速度です。 は「ボールの表面が風の何倍速く回るか」を表します。 が大きくなるほど循環が大きくなり、揚力係数も大きくなります。野球のカーブはおよそ 、サッカーの強い曲げ蹴りはそれより高くなります。
簡単な2次元の飛行を積分して軌道を描いてみましょう。抵抗とマグヌス力だけを入れます。
import numpy as np
def magnus_trajectory(u0, spin, cm=0.16, cd=0.006, dt=0.012, steps=900):
# 上から見た平面。x は前進、y は横方向。
pos = np.array([0.0, 0.0])
vel = np.array([u0, 0.0]) # 最初はまっすぐ発射
path = [pos.copy()]
for _ in range(steps):
sp = np.hypot(*vel) + 1e-9
# マグヌス: 速度に垂直 (-vy, vx)、符号はスピンが決める
a_magnus = cm * spin * np.array([-vel[1], vel[0]])
a_drag = -cd * sp * vel # 抵抗: 速度と逆向き
vel = vel + (a_magnus + a_drag) * dt
pos = pos + vel * dt
path.append(pos.copy())
if pos[0] > 30:
break
return np.array(path)
for s in (0.0, 1.0, 2.0):
p = magnus_trajectory(14.0, s)
print(f"spin={s}: 横方向のずれ = {p[-1, 1]:6.2f} m")
# spin=0.0: 横方向のずれ = 0.00 m
# spin=1.0: 横方向のずれ = 3.71 m
# spin=2.0: 横方向のずれ = 6.90 mスピンが2倍になっても、ずれはちょうど2倍にはなりません。軌道が曲がると速度の向きが変わり、力の向きも回り続けるからです。非線形です。
コードは 、、、、最大900ステップ、 で止まります。出力は spin 0 で横ずれ 0.00 m、spin 1 で 3.71 m、spin 2 で 6.90 m です。3.71 の2倍は 7.42 です。6.90 はそれより 0.52 m 短いです。軌道が曲がった瞬間、マグヌス加速度 も一緒に回ります。力が最初の横方向に固定されないので、ずれは線形に増えません。スピン比 では分母が減ると が大きくなります。発射速度を下げると、同じ でもより曲がります。野球のカーブの はこの比の一区間です。レイノルズを脇に置いたのは、揚力の式に Re が無いからです。
スピンのスライダーを左端から右端へ動かし、発射速度も下げると、破線を境に青い軌道が反転し、同じスピンのずれが大きくなります。
スピンを左端から右端へ動かすと、青い軌道が破線(無回転の基準線)を境に反対側へ反転します。発射速度を下げると、同じスピンでもより大きく曲がります。表面速度が流れに対して相対的に大きくなり、 が上がるためです。
ライフリングという逆さまの解法
ここで逆説が現れます。ロビンスが見つけた解は、回転をなくすことではなく、回転を制御することでした。銃身の内側にらせん状の溝(ライフリング、rifling)を刻み、弾を飛行方向の軸まわりに回転させます。
なぜこれが効くのでしょうか。マグヌス力は に比例します。ライフリングが与える回転軸は飛行方向、つまり が と平行です。二つのベクトルが平行なら外積はゼロです。横向きの力は消えます。
同時にこの軸回転はコマのようなジャイロ安定性を与えます。弾がでたらめに転がって得ていたランダムな回転を、一つの定まった軸に束ねてしまうのです。ランダムなマグヌスが消え、弾はまっすぐ飛びます。回転をなくすからではなく、横向きの力がゼロになる角度に回転をそろえて勝ったのです。
ライフリングの軸を飛行方向に合わせ損ねると横力が残る
マグヌス力は です。ライフリングが回転軸を飛行方向と平行に置くと、二つのベクトルは平行で外積は0です。横力は消えます。同時にこの軸回転はコマのようなジャイロ安定で、砲身の壁で跳ねて生まれていた無作為なスピンを一つの軸に束ねます。回転を消したのではありません。横力が0になる角度へ回したのです。
循環が0でなければ、物体のまわりに正味の旋回があります。円柱は で、その旋回を表面から引きずります。後流が下へ傾けば流体は下向き運動量を得て、ボールは上へ押されます。レイノルズをこの力の大きさを書く場所に置くと、 とすれ違います。Reynolds数が流れのスケールを分ける側は剥離には入っても、 の行には入りません。
関連記事
役に立ったらシェアしてください。