ISEGORIA / 数学百科事典
空気力学:翼はどう揚力を生むか
翼は通り過ぎる空気を下へ曲げ、空気は翼を押し返します。どれだけ曲げるかを説明するには三つの考えが要ります。翼型断面のまわりでは、薄い境界層の外の流れはほとんど理想流体で、揚力はただ一つの数、循環で決まり、それを鋭い後縁が選びます(クッタ条件)。ジューコフスキーは、この問題を厳密に解く等角写像を見つけました。有限の翼幅の翼では空気が翼端を回り込むので、渦が後ろへ流れ出して空気を吹き下ろし、揚力を後ろへ傾けて誘導抗力を生みます。プラントルの揚力線理論はこれを一つの積分方程式にしました。最後に、機体の抗力極曲線が、どこまで滑空でき、どれだけゆっくり沈下できるかを決めます。グライダーも旅客機も、この計算を中心に設計されています。
前提知識: 複素数、ベクトル解析、フーリエ級数
予想してから操作し、各実験の練習問題で理解を確かめてください。練習問題は三つの難易度で毎回新しく生成されます。グラフは数式の図示であり、証明ではありません。
1. ジューコフスキー翼型:等角写像による揚力
円柱を過ぎる流れは簡単です。一様流と二重わき出し、それに必要なら循環 \(\Gamma\) の渦を加えます。写像 \(z=\zeta+1/\zeta\) は \(\zeta=\pm1\) を除くどこでも等角なので、その流れを、円が写った形を過ぎる流れに移します。\(\zeta=1\) を通り、中心が原点の少し左かつ上にある円は、丸い前縁と反りをもち、\(z=2\) に鋭い後縁をもつ翼型になります。どの \(\Gamma\) も正しい理想流を与えますが、鋭い縁で速度を有限に保つのはただ一つ、円の後ろのよどみ点をちょうど \(\zeta=1\) に置くものだけです。これがクッタ条件で、それを使うとクッタ・ジューコフスキーの定理が揚力を与えます。円と迎角を動かし、圧力を読んでください。
計算例. 円の中心を \(\mu=-0.08+0.08i\) に置きます。半径は \(a=|1-\mu|=1.0830\)、後縁は中心から見て角度 \(\beta=\arcsin(0.08/1.083)=4.24^\circ\) 下にあるので、翼型の厚さは 9.7% で、迎角 \(-4.24^\circ\) で揚力がゼロになります。\(\alpha=5^\circ\) ではクッタ条件により \(\Gamma=4\pi\times1.083\times\sin9.24^\circ=2.184\)(\(U=1\))です。翼弦は 4.022 なので \(C_L=2\Gamma/c=1.086\) です。薄翼理論の \(2\pi(\alpha+\beta)=1.013\) は 7% 小さく、厚さはおよそ \(1+0.77\,t/c\) に比例して揚力を増やします。表面のまわりで圧力を積分しても同じ 1.086 が得られます。
注意点. これは非粘性・非圧縮の流れです。抗力はまったく予言されず(ダランベールのパラドックス)、失速もありません。実際の翼型は 12° から 18° のあたりで上面から流れがはく離して揚力を失いますが、ここでは揚力が増えつづけます。ジューコフスキー翼型の後縁は角度ゼロのとがった形で、実用の断面とは違い、前縁近くの圧力のピークも現代の断面より鋭くなります。
なぜ鋭い後縁が循環を決めるのですか?
理想流はどんな循環も許しますが、誤った循環だと空気は鋭い縁を回り込まなければならず、そこで速度が無限大になります。実際の空気にはそれができません。粘性の境界層がはく離し、流れが縁からなめらかに離れるようになるまで出発渦を放出します。ケルビンの定理により、翼には大きさが等しく向きが逆の循環が残り、それがちょうどクッタの値です。
2. 有限翼:後流渦と誘導抗力
実際の翼では、下面の高い圧力が翼端を回り込むので、揚力はそこでゼロにならなければなりません。ヘルムホルツの定理によれば渦線は空中で終われないので、循環 \(\Gamma(y)\) が翼幅に沿って変わるところでは必ず、強さ \(-d\Gamma/dy\) の渦が下流へ流れ出します。これらがあわせて翼に下向きの速度、つまり吹き下ろしを誘導し、各断面の有効迎角を \(\alpha_i\) だけ下げ、揚力を後ろへ傾けます。これが誘導抗力で、有限の翼幅の翼で揚力を作るための代価です。プラントルは \(\Gamma\) を正弦級数で書き、問題を連立一次方程式にしました。平面形、アスペクト比、ねじりを変え、分布を理想の楕円と比べてください。
計算例. ねじりのないアスペクト比 8 の楕円翼は純粋な楕円分布をもち、\(n>1\) で \(A_n=0\) なので \(e=1\) です。揚力傾斜は翼型の \(2\pi\) ではなく \(2\pi/(1+2/8)=5.03\)(1/rad)です。\(5^\circ\) では \(C_L=0.439\)、\(C_{D_i}=0.439^2/(8\pi)=0.0077\) です。同じアスペクト比の矩形翼は、翼端に翼弦の割に重すぎる荷重をかけるので \(e=0.937\)、\(C_L=0.422\) です。翼端の翼弦を翼根の 0.4 にテーパーをつけると \(e\) は 0.987 に戻り、ほぼ楕円になります。多くの翼にテーパーがあるのはこのためです。
注意点. 揚力線理論は、後退角がなくアスペクト比の大きい(およそ 4 以上の)翼で、流れが付着し、どの断面も揚力傾斜が \(2\pi\)(1/rad)であると仮定します。後流渦はまっすぐとしていますが、実際には翼幅の数倍のうちに巻き上がって二本の翼端渦になります。描いた渦の太さは、各帯から放出される強さの平方根に比例します。断面は対称なので揚力ゼロの迎角はゼロで、ウォッシュアウトは翼端を直線的に機首下げにねじります。
なぜグライダーの翼はあんなに細長いのですか?
誘導抗力は \(C_L^2/(\pi e\,AR)\) なので、同じ揚力ならアスペクト比に反比例して小さくなります。グライダーは大きな \(C_L\) でゆっくり飛び、そこでは誘導抗力が支配的なので、アスペクト比 20 から 40 は重さと構造の難しさに見合います。物理的には、長い翼幅はより大きな質量の空気をよりおだやかに押し下げます。抗力が支払っているのは、後流に残る運動エネルギーです。
3. 滑空の極曲線
機体の抗力は、おおよそ一定の部分(摩擦抗力と形状抗力)と、揚力係数の二乗で増える誘導抗力の和です。この放物線の抗力極曲線が性能を決めます。定常滑空では揚抗比が滑空比なので、最良の滑空は誘導抗力と有害抗力が等しくなる揚力係数で起こり、原点から接線を引けば見つかります。サーマルの中で高度を保ちたいグライダーは、代わりに最小の沈下率を求め、それはより遅い速度で得られます。重さは翼面荷重 \(W/S\) を通して、速度極曲線全体を形を変えずに速い方へずらします。
計算例. アスペクト比 18、\(e=0.9\)、\(C_{D0}=0.012\)、翼面荷重 350 N/m² のグライダーを海面高度で考えます。\(K=1/(\pi\times0.9\times18)=0.0196\) で、最良の滑空比は \(C_L=0.781\) で \(1/(2\sqrt{0.0196\times0.012})=32.6\) です。このときの速度は \(\sqrt{2\times350/(1.225\times0.781)}=27.0\) m/s、沈下率は \(27.0/32.6=0.83\) m/s です。最小沈下は \(C_L=\sqrt3\times0.781=1.354\)、速度 20.5 m/s で起こり、そこでは \(C_D=4C_{D0}\)、沈下率は 0.73 m/s です。バラストを積んで翼面荷重を二倍にしても滑空比は 32.6 のままで、どの速度も \(\sqrt2\) 倍になります。
注意点. 放物線の極曲線は理想化です。実際には \(C_{D0}\) は \(C_L\) や速度(レイノルズ数)によって変わり、本当の極曲線は測定で得ます。失速の限界 \(C_{L,\max}=1.4\) はフラップなしの典型的な値です。ここでは多くのグライダーと同じく最小沈下の点が失速の近くにあり、パイロットがサーマルの中を式より少し速く旋回するのはこのためです。空気の密度は海面高度の値です。
なぜ重さは最良の滑空比を変えないのですか?
滑空比は \(L/D=C_L/C_D\) で、揚力係数だけで決まり、重さにはよりません。重いグライダーは同じ \(C_L\) で速く飛ぶだけです。揚力 \(\tfrac12\rho V^2SC_L\) が大きな重さに等しくなければならないので \(V\propto\sqrt{W/S}\) となります。高度 1 メートルあたり同じ距離を、より速く進みます。条件のよい日にはそれがレースに勝ちます。