なぜ変化球が曲がるのか

> 國破れて　山河在り
城春にして　草木深し
時に感じて　花にも涙を濺ぎ
別れを恨んで　鳥にも心を驚かす
峰火　三月に連なり
家書　萬金に抵る
白頭掻いて　更に短かし
渾べて簪に　勝えざらんと欲す

こんにちわ. ブミブミ. 私わ[**東京大学理学部物理学科のJohn. T**](https://x.com/johnt_299792458/status/1724040226711216385?s=20)でわありません. Los Angels Dodgers#17 大谷翔平でもないです. **学部満期退学**をしたエセ科学者をご存知ですか？ 20231219は変化球がなぜ曲がるのか, 審議していきますよ. 
# 基礎方程式とベルヌーイの定理
ここでわ, 簡単のため空気は非圧縮で粘性のない理想流体とみなす. 空気の密度$\rho$, 速度ベクトル$\vec{u}=(u, v, w)$, 圧力$p$対して理想流体の基礎方程式
```math
\frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\vec{u})=0
```
```math
\frac{\partial \vec{u}}{\partial t}+(\vec{u}\cdot\nabla)\vec{u}=-\frac{\nabla p}{\rho}
```
が成り立つ. 一つ目の式は粒子数の保存
, 2つ目の式は運動量の保存(運動方程式)である. ここでは, 開いた流体の一部しか考えないのでエネルギーの保存則は使わない. ここでベクトル解析の公式より**渦度ベクトル**$\nabla\times\vec{u}=\vec{w}$として
```math
(\vec{u}\cdot\nabla)\vec{u}=\nabla(\frac{u^2}{2})-\vec{u}\times\vec{w}
```
が成り立つので, 定常な流れ$\frac{\partial \vec{u}}{\partial t}=0$の場合, 流線$\vec{u}$に沿って($\vec{u}\cdot\vec{u}\times\vec{w}=0$より)
```math
\frac{p}{\rho}+\frac{u^2}{2}=const. 
```
と表せる. これわ**ベルヌーイの定理**と呼ばれる.

# 流れ関数と速度ポテンシャル
簡単のため2次元で考えることにする. 今$\rho=const.$としているので, 基礎方程式から, $\nabla\cdot\vec{u}=0$がわかる. 
![IMG_0755.jpeg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/3626265/bfaa0e0c-da69-9b35-09d7-546fa1ed63b9.jpeg)
この時, 任意の2点AとPを結ぶ曲線について, 接線$\vec{dr}$に垂直な方向の単位ベクトルを, $\vec{n}$として次の量を考える. ただし, $\vec{n}$は$\vec{dr}$から左に90度回転したものとする.
```math
\Psi=\int_{A\rightarrow P}\vec{u}\cdot\vec{n}dr
```
この量は$A\rightarrow P$を横切る流れの総量を表し, その経路によらず, 流れ関数と呼ばれる.
流れ関数は次を満たす. 
```math
\frac{\partial \Psi}{\partial y}=u, -\frac{\partial \Psi}{\partial x}=v
```
流線であることと, $\Psi=const.$は等価である.
渦がない場合, (会っても良いが多価関数になるだけ)速度場$\vec{u}$に対してスカラーのポテンシャルを定義できる. それは**速度ポテンシャル**$\Phi$と呼ばれ, 次を満たす.
```math
u=\frac{\partial \Phi}{\partial x}, v=\frac{\partial \Phi}{\partial y}
```
この二つは*Caushy-Riemann*の関係式を満たしている.
```math
\begin{align}
u&=\frac{\partial \Phi}{\partial x}=\frac{\partial \Psi}{\partial y} \\
v&=\frac{\partial \Phi}{\partial y}=-\frac{\partial \Psi}{\partial x}
\end{align}
```
従って$\Phi$を実部, $\Psi$を虚部に持つ正則関数$f$を定義できる.　この正則関数は**複素速度ポテンシャル**と呼ばれる.

```math
f=\Phi+i\Psi
```
として定義できる. このポテンシャルで表せる流れを**ポテンシャル流れ**という.

# 複素速度ポテンシャルの例
![IMG_0757.jpeg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/3626265/1d4e2ad0-0c25-68c3-b072-9e304ba960ee.jpeg)

・一様流　$f=Uz$
・原点からの湧き出し $f=m\ln z$
・原点からの二重湧き出し(ダブレット)$f=-\frac{\mu}{z}$
・渦糸$f=ik\ln z$ $k$は渦の強さに対応している $k>0$が右回り, $k<0$が左回りである.

極座標表示すると, 
```math
f=ik\ln(re^{i\theta})=ik\ln r-k\theta
```
半径方向速度$u_r$と接線方向速度$u_\theta$はそれぞれ
```math
u_r=\frac{1}{r}\frac{\partial k\ln r}{\partial \theta}=0, 
u_\theta=-\frac{\partial k\ln r}{\partial r}=-\frac{k}{r}
```
となるので, 半径$r$の円周に沿った速度の積分は
```math
\oint u_\theta\cdot  2\pi rd\theta=-2\pi k\equiv\Gamma
```
と表せるので, 渦の強さを$k$が表しているとも言える. $\Gamma>0$は, こちらの方が便利なので$k$の代わりにこちらを使うことにする. $\Gamma>0$が左回り, $\Gamma<0$が右回りである.


# 動く円盤周りの流れ
半径$a$の十分薄い円盤が$x$方向に速さ$U$で動いている時のポテンシャルを求めよう.
極座標を用いると, 境界条件として, 円盤表面上では
```math
u_r(r=a)=\frac{\partial \Phi}{\partial r}=U\cos\theta
```
とあるので, 適当にローラン展開して考えれば, 複素速度ポテンシャル$f$が$r>a$では
```math
f=-\frac{Ua^2}{z}  
```
となることがわかる. これは実際原点周りの２重湧き出しで$\mu=Ua^2$と置いたものである.

![IMG_0758.jpeg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/3626265/df477bb1-70f2-f060-f7a7-5ae103021fa3.jpeg)


# 一様流れ中の回転しながら動く円盤
$x$方向の速度$-U$の一様流に, 速度$U$で動く円盤を入れる時の複素速度ポテンシャルは, 一様流のポテンシャルと, 動く円盤のポテンシャルを重ね合わせることで,
```math
f=-Uz-\frac{Ua^2}{z}
```
のように表せる.
この流れに, 円盤を左向きに回した循環を加えると, そのポテンシャルは
```math
f=-Uz-\frac{Ua^2}{z}-i\frac{\Gamma}{2\pi}\ln z
```
と表される.
極座標表示すると, $f=\Phi+i\Psi$は実部と虚部は
```math
\Phi=-U\left\{ r+\frac{a^2}{r} \right\}\cos\theta+\frac{\Gamma\theta}{2\pi}
```
```math
\Psi=-U\left\{ r-\frac{a^2}{r} \right\}\sin\theta-\frac{\Gamma\ln r}{2\pi}
```
となる.
円盤表面$(r=a)$での接線方向速度$u_\theta$は
```math
u_\theta=\frac{1}{a}\frac{\partial \Phi}{\partial \theta}=2U\sin\theta+\frac{\Gamma}{2\pi a}=2U\left\{\sin\theta+\frac{\Gamma}{4\pi Ua} \right\}
```
となり, $\theta$に依存している.
また, 流線は$r=a$で$\Psi=const.$となっていることから, 円盤そのものである.
ベルヌーイの定理を流線$r=a$に沿って用いると
```math
p(\theta)+\frac{\rho\left\{ (u_\theta)^2+(u_r)^2\right\}}{2}=const.
```
$r=a$では$u_r=const.$なので結局
```math
p(\theta)+2U^2\left\{ \sin\theta+\frac{\Gamma}{4\pi Ua}\right\}^2=p_0=const.
```
ともかける. 従って圧力$p$は$\theta$の関数として
```math
p(\theta)=p_0-\left\{ \sin\theta+\frac{\Gamma}{4\pi Ua} \right\}^2 
```
円盤表面で法線ベクトル$\vec{n}=(\cos\theta, \sin\theta)$をかけて積分すると円盤の受ける単位厚みあたりの力になる.
```math
\vec{F}=-\oint p(\theta)\vec{n}ad\theta
```
簡単な計算によって自明に次のように示される(**check by yourself**)
```math
F_x=0
```
```math
F_y=+\rho U\Gamma
```
これは[**クッタ・ジューコフスキーの定理**](https://ja.wikipedia.org/wiki/%E3%82%AF%E3%83%83%E3%82%BF%E3%83%BB%E3%82%B8%E3%83%A5%E3%82%B3%E3%83%BC%E3%83%95%E3%82%B9%E3%82%AD%E3%83%BC%E3%81%AE%E5%AE%9A%E7%90%86), あるいは**マグヌス効果**と呼ばれる. 動きながら回転する半径$a$の円盤は, 単位厚みあたりに進行方向$(+x)$と回転軸方向$(+z)$の外積の逆方向に, 大きさ$ \rho U\Gamma$の力を受ける. $\Gamma$は回転数に比例した量である. 回転軸方向と回転数, 速度が変化球の曲がる方向を決める, というわけだ. もし半径$a$の球の場合は, これを積分することで, 単位時間あたりの回転数を$n$とすると, 正味の力$\frac{4\pi a^3}{3}Un\rho$を受けると求められる.


# ストレート・フォーク・ナックル
![IMG_0766.jpeg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/3626265/4dcc7bec-7d73-1b67-aace-25a7724c9388.jpeg)

![IMG_0762.jpeg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/3626265/178761aa-ca1c-f892-9ed7-cd6813cc0a2a.jpeg)

マウンドからホームベースに向かう方向を$+x$方向, 鉛直上向き$+y$方向, マウンドから３塁ベースに向かう方向に$+z$方向にとる. 
ストレートは進行方向$+x$と回転軸方向$+z$を持つので, 重力だけでなく, 回転数に比例した上向きの揚力が働くことがわかる. これが俗に言う**ホップする**ストレートである.　フォークやナックルはストレートとほぼ同じだが, 回転数を抑えることで揚力を少なくし, ストレートに比べ, **落とす**ことができる.

# スライダー・カットボール・カーブ
![IMG_0765.jpeg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/3626265/4a489549-6654-4f8e-b2b7-23ce48640ed0.jpeg)

右投手ならスライダー・カットボール系のは回転軸が$+y$方向をむく. つまり, この揚力は$-x$方向, つまり左側に向く. (左投手だと逆になる) つまり**利き手と逆側**に向かって曲がっていることがわかる. 右投手のカーブの場合は, 回転軸が$+y$方向よりやや$+x$方向の成分も持つため, カット系の球よりも鉛直上むき成分の揚力を持つ. そのため, カーブは一度浮き上がるような軌道に見えるわけだ.

# シュート・シンカー・チェンジアップ
![IMG_0767.jpeg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/3626265/5c2a1592-ce74-2293-e540-39fb31763cfe.jpeg)

カット系と逆向きの回転軸を持つのでカット系の球と逆向きの揚力を受ける. つまり利き手側に変化していく. チェンジアップはシュートやストレートに比べて回転数を落とすので, 受ける揚力の大きさも小さくなるため, シュートに比べて落ちる.

# ジャイロボール
回転軸が進行方向に並行な回転はジャイロ回転と言われる. 回転軸が完全に進行方向と一致している場合はこのマグヌス力が働かない. そのため変化量は小さくなる.

# 球全体の受ける正味の揚力
![IMG_0768.jpeg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/3626265/9c6fb2b6-8be6-0d61-c309-5ef81499b8f8.jpeg)

単位時間あたりの回転数を$n$とする. 回転軸に垂直な面で切っているものとして考えよう. 図のように角度$\phi$をとると, 回転数と接線方向速度の関係から次が成り立つ.
```math
u_\theta(\phi)=2\pi n a\cos\phi=\frac{\Gamma}{a\cos\phi}
```
従って$\Gamma$と回転数$n$の間の関係が
```math
\Gamma=2\pi na^2\cos^2 \phi
```
とわかる. よって単位厚みあたりの力$\rho U\Gamma$を全ての$\phi$で積分すれば, 球全体にかかる力の大きさ$F_{total}$は
```math
F_{total}=\int_{0}^{\frac{\pi}{2}}\rho U 2\pi na^2\cos^2\phi\cdot a\sin\phi d\phi=\frac{4\pi a^3}{3}Un\rho
```


# 終わりに
ここまで長々変化球を流体力学で扱ってきたが, 筆者[John T(UT)](https://twitter.com/johnt_299792458/status/1727333273565802592?s=61&t=nUno8HDqjipVRSRosXDOlQ)のできることは, 所詮ボールが投げられたあとの話である. どう投げるかではない. しかし[Justin Verlander(HOU)](https://www.mlb.com/player/justin-verlander-434378)や[Shohe Ohtani(LAD)](https://www.mlb.com/player/shohei-ohtani-660271), [Spencer Strider(ATL)](https://www.mlb.com/player/spencer-strider-675911)のようなどんなに優れた選手・WBC決勝のような緊張した場面でも, 物理学の下にある.　そのような場面は彼らのような優れた選手と, 途方もない努力が必要であろうが, 物理学はそこまでではない.(少なくとも身体的な才能や野球の技術は必要ない) みなさんもこの記事をきっかけに物理学に興味を持って, 理学部物理学科に進学すると良いではないか. もちろん, MLBを目指しても良いね.

# Reference
Referenceという単語, eが多すぎないか. そんなんだから英語で最も使われるアルファベットがいつまでも**e**のままなんだ(審議あり)
杉山弘, 遠藤剛, 新井隆景. 流体力学 第２版. 森北出版株式会社, 2014
巽友正, 新物理学シリーズ21 流体力学. 培風館, 1982