本文へ移動
ISEGORIABenjamin Haire
English

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° のあたりで上面から流れがはく離して揚力を失いますが、ここでは揚力が増えつづけます。ジューコフスキー翼型の後縁は角度ゼロのとがった形で、実用の断面とは違い、前縁近くの圧力のピークも現代の断面より鋭くなります。

なぜ鋭い後縁が循環を決めるのですか?

理想流はどんな循環も許しますが、誤った循環だと空気は鋭い縁を回り込まなければならず、そこで速度が無限大になります。実際の空気にはそれができません。粘性の境界層がはく離し、流れが縁からなめらかに離れるようになるまで出発渦を放出します。ケルビンの定理により、翼には大きさが等しく向きが逆の循環が残り、それがちょうどクッタの値です。

参考文献: MIT OpenCourseWare · 16.100 Aerodynamics (Fall 2005)

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 は重さと構造の難しさに見合います。物理的には、長い翼幅はより大きな質量の空気をよりおだやかに押し下げます。抗力が支払っているのは、後流に残る運動エネルギーです。

参考文献: MIT OpenCourseWare · 16.100 Aerodynamics (Fall 2005)

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 メートルあたり同じ距離を、より速く進みます。条件のよい日にはそれがレースに勝ちます。

参考文献: MIT OpenCourseWare · 16.100 Aerodynamics (Fall 2005)

次に読む

線形代数:変換の幾何フーリエ解析:波で関数を作る力学系:安定性とカオス最適化:最良の選択の幾何確率:不確実性から学ぶ群論:対称性を代数にする代数的位相幾何:穴を見つける数値解析:計算が誤解を招くとき数論:整数のパターングラフ理論:経路・木・ネットワーク偏微分方程式:動く場情報理論:不確実性と符号変分法:経路と原理古典力学:運動と力電磁気学:場と誘導光学:光線・波・色熱力学:エネルギー・仕事・エントロピー量子力学:振幅とスピン複素解析:写像・留数・調和場流体力学:流れ・圧力・渦度統計と推測:データの中の信号特殊相対性理論:空間・時間・光微分幾何:曲率と形統計力学:ミクロ状態と温度測度論:大きさ・近似・収束マルコフ連鎖:遷移・定常性・吸収微分形式:循環・カール・引き戻し一般相対性理論:曲率・時計・光関数解析:ノルム・射影・作用素プラズマ物理:遮蔽・軌道・波リー群とリー代数:連続的な対称性ハミルトン力学:相空間とその幾何確率過程:ブラウン運動とノイズ固体物理:結晶の中の波制御理論:フィードバック・極・安定性論理と計算可能性:何が計算できるか双曲幾何:平行線が増える世界剛体力学:回転・転倒・歳差楕円曲線:足し算ができる幾何原子物理:軌道とスペクトル四元数:掛け算としての回転押し出しと引き戻し:写像を通した積分結び目理論:絡まりを見分けるフラクタル幾何:整数の間の次元量子情報:もつれとその限界宇宙論:膨張する宇宙セル・オートマトン:局所規則から生まれる計算複雑系:単純な多数の部品から生まれる秩序ガロア理論:方程式の対称性磁性とイジング模型:整列から生まれる秩序振動子と脱進機:時計はどうやって時を刻むか歯車と機構:運動を正確に伝える水晶振動子:時を刻む結晶スピーカー:ムービングコイル型ドライバーフィルターとクロスオーバー:音をドライバーに振り分ける室内音響:部屋はスピーカーの一部ウェーブレット:信号を拡大して見るレーザー物理:自らを複製する光表現論:行列として作用する群半導体物理:バンド・ドーピング・接合組合せ論:組み立てて数える音響学:弦、板、動く音源スイッチング電源:刻んで変換するゲーム理論:相手も選ぶときの戦略原子核物理:結合、崩壊、トンネル効果ピストンエンジン:クランクからサイクルまでエルゴード理論:時間平均が空間平均に等しくなるとき連続体力学:固体はどう伸び、曲がり、座屈するか電動機:電流をトルクに変える重力波:空間の幾何のさざ波ロボットアーム:運動学とヤコビ行列リーマン面:多価関数のすみか射影幾何:平行線が交わる場所恒星天体物理:星の輝きと死構造力学:橋はどう荷重を支えるか符号理論:誤りを訂正する超伝導:ゼロ抵抗と磁束の排除アダプティブサイクルエンジン:一基で二つの性格能動流れ制御:空気の噴流で操る

数学百科事典に戻る