クッタ・ジューコフスキーの定理は、空気力学における基本定理であり、一様な流体中を一定の速度で移動する翼(および円筒を含む任意の2次元物体)の揚力を計算するために使用される。この定理は、翼が流体中を移動する速度、流体の密度、および翼の周りの循環に、翼によって生成される揚力を関連付ける。循環は、ループに接する流体の速度成分の翼を囲む閉ループの周りの線積分として定義される。 [1]この定理は、20世紀初頭にその主要なアイデアを最初に開発したマーティン・クッタとニコライ・ジューコフスキー(またはジューコフスキー)にちなんで名付けられました。クッタ・ジューコフスキーの定理は非粘性理論ですが、一般的な空気力学アプリケーションにおける実際の粘性流の近似として優れています。[2]
クッタ・ジューコフスキーの定理は、マグヌス効果が横方向の力(マグヌス力と呼ばれる)を回転に関連付けるのとよく似て、揚力を循環に関連付けます。 [3]しかし、ここでの循環は翼の回転によって引き起こされるものではありません。翼の存在下での流体の流れは、並進流と回転流の重ね合わせと考えることができます。この回転流は、キャンバー、迎え角、および翼の鋭い後縁の影響によって引き起こされます。翼を囲む竜巻のような渦と混同しないでください。翼からかなり離れている場合、回転流は線渦(回転線が2次元平面に垂直な線)によって引き起こされるものと見なすことができます。クッタ・ジューコフスキーの定理の導出では、通常、翼は円柱上にマッピングされます。多くの教科書では、定理は円柱とジューコフスキー翼について証明されていますが、一般的な翼にも当てはまります。
揚力の計算式
この定理は、固定翼(または無限翼幅の任意の形状)の周りの2次元流れに適用されます。翼の単位翼幅あたりの揚力は[4]で与えられます。
ここで、およびは翼の上流側の流体の密度と流速であり、およびは線積分として定義される循環である。
翼を囲む閉じた輪郭の周りを負の方向(時計回り)にたどる。以下で説明するように、この経路はポテンシャルフロー領域内にある必要があり、シリンダーの境界層内ではない。積分関数は、曲線に接する方向の局所流体速度の成分であり、曲線上の無限小長さである。式(1)は、クッタ・ジューコフスキーの定理の一形式である。
Kuethe と Schitzer は、Kutta-Joukowski の定理を次のように述べています: [5]
- いかなる断面の直円筒に作用する単位長さあたりの力は、次の式に等しく、その方向に垂直である。
循環とクッタ状態
揚力を生み出す翼型は、キャンバーを持つか、正の迎え角(翼弦線と翼型よりかなり上流の流体の流れとの間の角度)で動作します。さらに、翼型は鋭い後縁を持たなければなりません。
実際の流体は粘性があり、これは流体の速度が翼上でゼロになることを意味します。プラントルは、 と定義される 大きなレイノルズ数と小さな迎え角に対して、薄い翼の周りの流れは、物体に近い境界層と呼ばれる狭い粘性領域と、外側の非粘性流れ領域で構成されることを示しました。クッタ・ジューコフスキーの定理を適用するには、この境界層の外側にループを選択する必要があります。(たとえば、翼の表面に対応するループを使用して計算された循環は、粘性流体の場合はゼロになります。)
鋭い後縁要件は、物理的には、翼の下部と上部の表面に沿って移動する流体が滑らかに合流し、翼の後縁の周囲に流体が移動しない流れに対応します。これは、クッタ条件として知られています。
クッタとジュコウスキーは、レイノルズ数が大きく迎え角が小さい流れにおける薄い翼の圧力と揚力を計算する場合、クッタ条件が課せられていれば、翼の外側の領域全体で流れが非粘性であると仮定できることを示しました。これはポテンシャル流れ理論として知られ、実際に非常にうまく機能します。
導出
以下に 2 つの導出を示します。1 つ目は物理的な洞察に基づいたヒューリスティックな議論です。2 つ目は基本的なベクトル解析と複素解析を必要とする形式的かつ技術的な議論です。
ヒューリスティックな議論
経験的な議論として、弦長 が薄く翼幅が無限の翼が密度の空気中を移動することを考えてみましょう。翼を向かい合う流れに対して傾けると、翼の片側に空気速度が生じ、反対側に空気速度が生じます。循環は次のようになります。
翼の両側の圧力差はベルヌーイの式を適用することで求めることができます。
空気に働く下向きの力は、単位スパンあたり、
そして翼の上向きの力(揚力)は
この定理の微分バージョンはプレートの各要素に適用され、薄翼理論の基礎となります。
正式な導出
まず、任意の断面を持つ円筒の各単位長さに働く力を計算します。[6]この単位長さあたりの力(以下、単に力と呼ぶ)を とします。すると、全体の力は次のようになります。
ここで、C はシリンダーの境界線、は流体の静圧、はシリンダーに垂直な単位ベクトル、 ds は断面の境界線の円弧要素です。ここで、法線ベクトルと垂直線の間の角度をとします。上記の力の成分は次のようになります。
ここで重要なステップが来ます。使用する 2 次元空間を複素平面として考えます。したがって、すべてのベクトルは複素数として表すことができ、その最初の要素は複素数の実数部に等しく、2 番目の要素は虚数部に等しくなります。すると、力は次のように表すことができます。
次のステップは、力の複素共役を取って何らかの操作を行うことです。
表面セグメントds は、それに沿った変化dzと次のように関連しています。
これを積分に戻すと、結果は次のようになります。
ここで、積分から圧力を取り除くためにベルヌーイ方程式が使用されます。解析全体を通して、外部力場は存在しないものと仮定します。流れの質量密度は、圧力と速度の関係は次のようになります。
これにより、力は次のようになります。
残されたステップは 1 つだけです。流れの複素ポテンシャルを導入することです。これは速度成分と次のように関連しています。ここで、アポストロフィは複素変数zに関する微分を表します。速度は境界線Cに接しているので、次のようになります。したがって、力の目的の式が得られます。
これはブラシウスの定理と呼ばれます。
ジュコウスキーの公式を得るには、この積分を評価する必要があります。複素解析から、正則関数はローラン級数として表すことができることがわかっています。問題の物理的性質から、複素ポテンシャルの導関数は次のようになると推測されます。
速度は無限大では有限のままなので、関数には高次の項は含まれません。したがって、 は無限大における複素ポテンシャルの導関数を表します。次のタスクは の意味を見つけることです。上記の級数に 留数定理を使用すると、次のようになります。
次に、上記の統合を実行します。
最初の積分は、次のように表される循環として認識されます。2 番目の積分は、いくつかの操作を行うことで評価できます。
ここで、流れ関数は、円筒のC境界自体が流線であるため、流れ関数は変化せず、 となります。したがって、上記の積分はゼロになります。結果として、
次の級数の平方をとります。
これを Blasius–Chaplygin の公式に代入し、留数定理を使用して積分を実行します。
したがって、クッタ・ジューコフスキーの公式は次のようになります。
より複雑な状況における揚力
非粘性ポテンシャル流理論の枠組み内でクッタ・ジューコフスキー定理によって予測される揚力は、流れが定常かつ分離していない限り、実際の粘性流に対しても非常に正確です。[7] クッタ・ジューコフスキー定理を導く際には、非回転流という仮定が使用されました。多くの非定常流の場合のように、物体の外側に自由渦がある場合、流れは回転します。流れが回転している場合、揚力を導くにはより複雑な理論を使用する必要があります。以下にいくつかの重要な例を示します。
- 小さな迎え角で衝動的に流れが始まった
- 翼を急加速したり迎え角を設定したりすることで得られる衝撃的に始まる流れの場合、後縁で渦シートが連続的に放出され、揚力は不安定または時間に依存します。小さな迎え角の開始流れの場合、渦シートは平面経路をたどり、時間の関数としての揚力係数の曲線はワグナー関数によって与えられます。[8]この場合、初期揚力はクッタ-ジュコウスキーの式で与えられる最終揚力の半分です。[9]翼が約7弦長の距離を移動したとき、揚力は定常値の90%に達します。
- 大きな迎え角で衝動的に流れが始まった
- 迎え角が十分に高い場合、後縁渦シートは最初は螺旋状になり、揚力は初期時点で特異(無限大)になります。[10] 通常想定される単調に増加する揚力曲線に達する前に、揚力は非常に短い期間低下します。
- 鋭い前縁を持つ翼の大きな迎え角での流れの開始
- 平板の場合と同様に、前縁も鋭い場合は、前縁でも渦が発生し、前縁渦の役割は2つあります。1) 前縁に近いときに揚力を増加させ、ワグナー揚力曲線を上昇させます。2) 前縁に対流すると揚力に悪影響を及ぼし、揚力が減少する方向に移動する新しい後縁渦螺旋を引き起こします。このタイプの流れでは、渦力線(VFL)マップ[11]を 使用して、さまざまな状況(開始流れよりも多くの状況を含む)でのさまざまな渦の影響を理解することができ、揚力を増強または減少させるための渦制御を改善するために使用できます。渦力線マップは、渦力線が表示される2次元マップです。流れのどの地点の渦でも、その揚力の寄与は、その速度、循環、および流線と渦力線の間の角度の余弦に比例します。したがって、渦力線マップは、特定の渦が揚力を生み出すか、揚力を阻害するかを明確に示します。
- ラガリー定理
- (質量)源が物体の外側に固定されている場合、この源による力の補正は、外部源の強度と、この源以外のすべての原因によってこの源に誘導される速度の積として表すことができます。これはラガリー定理として知られています。[12] 2次元の非粘性流れの場合、古典的なクッタ・ジュコウスキー定理は抗力がゼロであると予測します。ただし、物体の外側に渦がある場合、誘導揚力に似た形で渦誘導抗力が発生します。
- 一般化されたラガリー定理
- 自由渦および束縛渦度がなく渦生成のない物体外部の物体については、一般化されたラガリー定理が成り立ち、[13]これによれば、力は内部特異点(各物体内部の像渦、源、二重点)の強度と、この物体内部以外のすべての原因によってこれらの特異点に誘導される速度の積として表される。各内部特異点による寄与を合計すると、全体の力となる。外部特異点の運動も力に寄与し、この寄与による力の成分は特異点の速度に比例する。
- 多体回転流における各体の個々の力
- 複数の自由渦と複数の物体に加えて、物体表面に束縛渦と渦生成がある場合、一般化されたラガリー定理は依然として成立しますが、渦生成による力が存在します。この渦生成力は、渦生成率と生成中の渦対間の距離に比例します。このアプローチでは、すべての原因(内部特異点、外部の渦と物体、すべての特異点と物体の運動、および渦生成)を考慮した明示的かつ代数的な力の式が、追加の特異点によって表される他の物体の役割とともに、各物体に対して個別に成立します[14]。したがって、物体に応じた力の分解が可能です。
- 一般的な3次元粘性流れ
- 一般的な3次元の粘性非定常流れの場合、力の式は積分形式で表されます。渦度モーメントなどの特定の流量の体積積分は力に関連しています。現在、さまざまな形式の積分アプローチが、境界のない領域[9] [15] [16]と人工的に切り取られた領域[17]に対して利用可能です。Kutta Joukowskiの定理は、2次元翼に適用され、流れが定常で分離されていない場合、これらのアプローチから復元できます。
- 翼の揚力線理論、翼端渦、誘導抗力
- 翼には有限の幅があり、翼のどの部分でも循環は翼幅方向によって変化します。この変化は、渦度保存則またはケルビンの循環保存定理により、後流渦と呼ばれる縦渦の放出によって補正されます。これらの縦渦は、翼幅に近い距離で分離された 2 つの反対方向に回転する強い螺旋に合流し、相対湿度が高い場合はその中心が見えることがあります。後流渦を半無限の直線渦の列として扱うと、よく知られている揚力線理論が導かれます。この理論によると、翼の揚力は、クッタ-ジュコフスキーの定理を使用した純粋に 2 次元の理論で予測されるものよりも小さくなります。これは、後流渦によって追加されたダウンウォッシュが翼の迎え角に上流から影響するためです。これにより、翼の有効迎え角が減少し、特定の迎え角で生成される揚力の量が減少し、この失われた揚力を回復するためにより高い迎え角が必要になります。この新しいより高い迎え角では、抗力も増加しています。誘導抗力は、2D 翼の揚力曲線の傾斜を効果的に減少させ、迎え角を増加させます(同時に の値も減少します)。
参照
参考文献
- ^ アンダーソン、JD Jr. (1989)。「圧力、温度、密度高度」。フライト入門(第3版)。ニューヨーク:マグロウヒル。pp. 100–103。ISBN 0-07-001641-0。
- ^ Liu, LQ; Zhu, JY; Wu, JZ (2015). 「 2次元定常粘性流および圧縮性流における揚力と抗力」。流体力学ジャーナル。784 :304–341。Bibcode : 2015JFM ...784..304L。doi :10.1017/jfm.2015.584。S2CID 125643946 。
- ^ 「回転シリンダーの揚力」NASAグレン研究センター。2010年11月9日。2014年1月11日時点のオリジナルよりアーカイブ。2013年11月7日閲覧。
- ^ Clancy, LJ (1975).空気力学. ロンドン: ピットマン. セクション 4.5. ISBN 0-273-01120-0。
- ^ キューテ、AM; JD シェッツァー (1959)。空気力学の基礎。ニューヨーク:ジョン・ワイリー&サンズ。セクション4.9。ISBN 0-471-50952-3。
- ^ Batchelor, GK (1967).流体力学入門. ケンブリッジ大学出版局. p. 406.
- ^ Anderson, J. (2010).空気力学の基礎. 航空宇宙工学におけるマグロウヒルシリーズ. ニューヨーク: マグロウヒル教育.
- ^ ワーグナー、H. (1925)。 「Über die Entstehung des dynamischen Auftriebes von Tragflügeln」。Z.アンジュー。数学。メカ。 5 (1): 17-35。Bibcode :1925ZaMM....5...17W。土井:10.1002/zamm.19250050103。
- ^ ab Saffman, PG (1992). Vortex Dynamics . ニューヨーク: Cambridge University Press. ISBN 0-521-42058-X。
- ^ Graham, JMR (1983). 「始動流における翼の揚力」. Journal of Fluid Mechanics . 133 : 413–425. Bibcode :1983JFM...133..413G. doi :10.1017/S0022112083001986. S2CID 123501457.
- ^ Li , J .; Wu, ZN (2015). 「追加の前縁後縁渦が存在する場合のワグナー問題に対する不安定揚力」。流体力学ジャーナル。769 :182–217。Bibcode : 2015JFM ...769..182L。doi :10.1017/jfm.2015.118。S2CID 121892071 。
- ^ Milne-Thomson, LM (1968).理論流体力学. 香港: Macmillan Education. p. 226.
- ^ Wu, CT; Yang, FL; Young, DL (2012). 「自由渦を伴う一般化された2次元ラガリー定理と流体-物体相互作用問題へのその応用」(PDF) . Journal of Fluid Mechanics . 698 : 73–92. Bibcode :2012JFM...698...73W. doi :10.1017/jfm.2012.45. S2CID 120656935.
- ^ Bai, CY; Li, J.; Wu, ZN (2014). 「渦発生を伴う多重渦および多重翼流に対する一般化クッタ・ジュコウスキー定理 - 一般モデル」.中国航空学誌. 27 (5): 1037–1050. doi : 10.1016/j.cja.2014.03.014 .
- ^ Wu, JC (1981). 「粘性流における空気力とモーメントの理論」. AIAA Journal . 19 (4): 432–441. Bibcode :1981AIAAJ..19..432W. doi :10.2514/3.50966.
- ^ Howe, MS (1995). 「非圧縮性流体中の物体にかかる力とモーメントについて、高レイノルズ数での剛体と気泡への応用」Quarterly Journal of Mechanics and Applied Mathematics . 48 (3): 401–425. doi :10.1093/qjmam/48.3.401.
- ^ Wu, JZ; Lu, XY; Zhuang, LX (2007). 「局所的な流れ構造により物体に作用する積分力」. Journal of Fluid Mechanics . 576 : 265–286. Bibcode :2007JFM...576..265W. doi :10.1017/S0022112006004551. S2CID 122293574.
文献
- ミルン・トムソン、LM (1973)理論空気力学、ドーバー出版、ニューヨーク
