クッタ・ジュコーフスキーの定理

Heuristic argument

コード長 ${\displaystyle c}$ スパン長無限の薄い翼体が、密度 ${\displaystyle \rho }$ の空気中を移動する。このとき翼体を傾けて2つの翼面の一方の速度が ${\displaystyle V}$ もう一方の面の速度が ${\displaystyle V+v}$ となったとき、

循環はつぎのように表される。

${\displaystyle \Gamma =Vc-(V+v)c=-vc.\,}$

翼体上下の圧力差 ${\displaystyle \Delta P}$ は ベルヌーイ式によりつぎのように導かれる。

${\displaystyle {\frac {\rho }{2}}(V)^{2}+(P+\Delta P)={\frac {\rho }{2}}(V+v)^{2}+P,\,}$

${\displaystyle {\frac {\rho }{2}}(V)^{2}+\Delta P={\frac {\rho }{2}}(V^{2}+2Vv+v^{2}),\,}$

${\displaystyle \Delta P=\rho Vv\qquad {\text{(ignoring }}{\frac {\rho }{2}}v^{2}),\,}$

単位スパン長あたりの揚力は

${\displaystyle L'=c\Delta P=\rho Vvc=-\rho V\Gamma \,}$

となる。

この式の 微分形 を平板の要素に当て嵌めたものが thin-airfoil theory の基礎となる。

厳密な導出

先に要点をまとめると、

1. 揚力を物体表面上の圧力場の積算に書き表す。
2. ベルヌーイ則にしたがい圧力を速度平方に直す。ブラジウス式が得られる。
3. 速度場を級数で一般化。遠方速度Uと物体表面速度uの和となる
4. 速度平方 UU+Uu+uu の周回積分によりUUとuuが消え、Uuの積算が残る。これはUΓに等しい。

以上は、特定の条件下について複素解析により証明可能である。

さらに見る z軸単位長当たりの力 (以降は単に力と呼ぶ)を ...

複素解析による証明

対象となる回転円筒の回転軸がZ軸に一致するものとし、任意のXY断面を考えて、ここに働く[2] z軸単位長当たりの力 (以降は単に力と呼ぶ)を ${\displaystyle \mathbf {F} \,}$とする。粘性無視でき圧力のみ作用する場合に次のように表される。 　${\displaystyle \mathbf {F} =-\oint _{C}p\mathbf {n} \,ds}$ ここで C はXY面における円筒の境界線、${\displaystyle p}$ は 静圧、${\displaystyle \mathbf {n} \,}$ は 円筒面上の外向き単位法線、ds は断面上境界線の微小円弧要素長とする。ここで ${\displaystyle \phi }$ をx軸と微少要素のなす角度とすると、上述の力は成分ごとにつぎのようにあらわされる。 ${\displaystyle F_{x}=-\oint _{C}p\sin \phi \,ds\quad ,\qquad F_{y}=\oint _{C}p\cos \phi \,ds.}$ ここで対象のXを実数軸、Yを虚数軸として複素平面に持ち込む。 前述の力は次のように表される。${\displaystyle \phi }$ はX軸すなわち実数軸に対する角度であるから複素平面上の偏角として読み替えられる。 ${\displaystyle F=F_{x}+iF_{y}=-\oint _{C}p(\sin \phi -i\cos \phi )\,ds.}$ 複素ポテンシャルの定義に合わせて ${\displaystyle F_{y}}$ を正負反転する。ここでは ${\displaystyle F}$ の複素共役で表す。 ${\displaystyle {\bar {F}}=-\oint _{C}p(\sin \phi +i\cos \phi )\,ds=-i\oint _{C}p(\cos \phi -i\sin \phi )\,ds=-i\oint _{C}pe^{-i\phi }\,ds.}$ 面要素長 ds に対応する複素要素は ${\displaystyle dz=dx+idy=ds(\cos \phi +i\sin \phi )=ds\,e^{i\phi }\qquad \Rightarrow \qquad d{\bar {z}}=e^{-i\phi }ds.}$ まとめるとつぎの式が得られる。 ${\displaystyle {\bar {F}}=-i\oint _{C}p\,d{\bar {z}}.}$ ここでベルヌーイ則を適用し、圧力を速度式に置き換える。ここでは外力なしが前提である。空気の密度 ${\displaystyle \rho }$ および圧力 ${\displaystyle p}$ と速度 ${\displaystyle v=v_{x}+iv_{y}}$ は次の関係にある。 ${\displaystyle p=p_{0}-{\frac {\rho |v|^{2}}{2}}.}$ これにより力 ${\displaystyle F}$ は、 ${\displaystyle {\bar {F}}=-ip_{0}\oint _{C}d{\bar {z}}+i{\frac {\rho }{2}}\oint _{C}|v|^{2}\,d{\bar {z}}={\frac {i\rho }{2}}\oint _{C}|v|^{2}\,d{\bar {z}}.}$ となる。 複素ポテンシャル ${\displaystyle w=f(z)}$ を導入する。 速度成分との関係は ${\displaystyle w'=v_{x}-iv_{y}={\bar {v}}}$ ここでアポストロフィは複素変数 z での微分演算を示す。 ここで速度は境界 C に対し接線方向であるため、${\displaystyle v=\pm |v|e^{i\phi }.}$ したがって ${\displaystyle v^{2}d{\bar {z}}=|v|^{2}dz,\,}$ これにより次のブラジウス式と呼ばれる式が導かれる。 ${\displaystyle {\bar {F}}={\frac {i\rho }{2}}\oint _{C}w'^{2}\,dz,}$ ジューコフスキー式に至るために、この積分量を評価する。wはC上およびCより外側で正則でありwの導関数ｗ’も同じく正則でまたローラン級数で表されてよい。複素ポテンシャル ${\displaystyle w}$ の導関数は次のとおり、z=0周りのローラン級数で表される: ${\displaystyle w'(z)=a_{0}+{\frac {a_{1}}{z}}+{\frac {a_{2}}{z^{2}}}+\dots .}$ w'の速度場がC上で接線方向のみであることから、a0とa1の項のみ残して他は消える。（係数を求めるためのC上の積分操作でゼロとなる。[要確認]） （　${\displaystyle a_{0}\,}$は無限遠方速度と対応することから：　${\displaystyle a_{0}=v_{x\infty }-iv_{y\infty }\,}$.） ${\displaystyle a_{1}\,}$についてみるとローラン級数の定義より ${\displaystyle a_{1}={\frac {1}{2\pi i}}\oint _{C}w'\,dz}$ であるが、w'の積分は定義どおり循環値Γであるから　${\displaystyle a_{1}={\frac {\Gamma }{2\pi i}}}$ となる。 ${\displaystyle a_{1}\,}$を元の級数に戻すと ${\displaystyle w'(z)=a_{0}+{\frac {a_{1}}{z}}=a_{0}+{\frac {\Gamma }{2\pi iz}}}$ 平方をとると ${\displaystyle w'^{2}(z)=a_{0}^{2}+{\frac {a_{0}\Gamma }{\pi iz}}-{\frac {\Gamma ^{2}}{4\pi ^{2}z^{2}}}.}$ これを前述のブラジウス式に戻す。C上の周回積分により第2項だけ残して他は消える。 ${\displaystyle {\bar {F}}={\frac {i\rho }{2}}\left[2\pi i{\frac {a_{0}\Gamma }{\pi i}}\right]=i\rho a_{0}\Gamma =i\rho \Gamma (v_{x\infty }-iv_{y\infty })=\rho \Gamma v_{y\infty }+i\rho \Gamma v_{x\infty }=F_{x}-iF_{y}.}$ これにてクッタ・ジューコフスキー式: ${\displaystyle F_{x}=\rho \Gamma v_{y\infty },\quad \qquad F_{y}=-\rho \Gamma v_{x\infty }.}$ が得られる。

閉じる