ベンチュリ管 を空気が流れている。管の太さが小さくなると速度が増加するが、それには圧力の減少を伴う。圧力の変化は水柱の高さの差に現れる。
ベルヌーイの定理 (ベルヌーイのていり、英語 : Bernoulli's principle )またはベルヌーイの法則 とは、完全流体 のいくつかの特別な場合において、ベルヌーイの式 と呼ばれる運動方程式 の第一積分 が存在することを述べた定理 である。
概要
ベルヌーイの式は流体 の速さ と圧力 と外力 のポテンシャル の関係を記述する式で、力学 的エネルギー保存則 に相当する。この定理により流体の挙動を平易に表すことができる。
ダニエル・ベルヌーイ (Daniel Bernoulli 、1700年 - 1782年)によって1738年 に発表された。なお、運動方程式からのベルヌーイの定理の完全な誘導はその後の1752年 にレオンハルト・オイラー により行われた[ 1] 。ベルヌーイの定理は粘性のない流体である完全流体のとき成り立つ。同一流線上で運動エネルギーと圧力の二つの力の和が一定であることを言う定理であり、流体の流れが速い場所では圧力が低いことを言う定理ではない。[ 2] 身近なベルヌーイの定理の使用例として、鳥や飛行機、ビル風 の一部、霧吹き と内燃機関のキャブレター 、自動車部品のウイング やディフューザー 、野球ボール やゴルフボール が曲がる現象、速い列車 が側を通過するときに吸い寄せられる現象などがある。
分類
ベルヌーイの定理は適用する非粘性流体 の分類に応じて様々なタイプに分かれるが、大きく二つのタイプに分類できる。
外力が保存力 であること 、バロトロピック性 (密度が圧力のみの関数となる)という条件に加えて、
定常 流という条件で成り立つ法則(I)
渦なしの流れ という条件で成り立つ法則(II)
である。
(I)の法則は流線 上(正確にはベルヌーイ面上)でのみベルヌーイの式が成り立つという制限があるが、(II)の法則は全空間 で式が成立する。
最も典型的な例である
外力のない非粘性 ・非圧縮性流体 の定常な流れに対して
1
2
v
2
+
p
ρ
=
c
o
n
s
t
a
n
t
{\displaystyle {\frac {1}{2}}v^{2}+{\frac {p}{\rho }}=\mathrm {constant} }
が流線上で成り立つ。ただし、
v
{\displaystyle v}
は流体の速さ、
p
{\displaystyle p}
は圧力、
ρ
{\displaystyle \rho }
は密度 を表す。
や
一様重力のもとでの非粘性・非圧縮流体の定常な流れに対して
1
2
v
2
+
p
ρ
+
g
z
=
c
o
n
s
t
a
n
t
{\displaystyle {\frac {1}{2}}v^{2}+{p \over \rho }+gz=\mathrm {constant} }
が流線上で成り立つ。ただし、
v
{\displaystyle v}
は速さ、
p
{\displaystyle p}
は圧力、
ρ
{\displaystyle \rho }
は密度、
g
{\displaystyle g}
は重力加速度 の大きさ、
z
{\displaystyle z}
は鉛直方向の座標を表す。
は(I)のタイプに属する。
(II)を「一般化されたベルヌーイの定理」と呼ぶこともある。
基本形
完全流体の運動方程式からベルヌーイの定理を導出する[ 3] 。
オイラー方程式
バロトロピック性 ρ =ρ (p ) と外力が保存力であることを仮定すると、非粘性流体の運動を記述するオイラー方程式
D
v
D
t
=
−
1
ρ
∇
p
+
f
{\displaystyle {D{\boldsymbol {v}} \over Dt}=-{1 \over \rho }\nabla p+{\boldsymbol {f}}}
は
∂
v
∂
t
−
v
×
(
∇
×
v
)
+
∇
{
v
2
2
+
∫
d
p
ρ
+
Ω
}
=
0
{\displaystyle {\partial {\boldsymbol {v}} \over \partial t}-{\boldsymbol {v}}\times \left(\nabla \times {\boldsymbol {v}}\right)+\nabla \left\{{v^{2} \over 2}+\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}\right\}=0}
と変形できる。ただし、
v
{\displaystyle {\boldsymbol {v}}}
は速度 ベクトル、
p
{\displaystyle p}
は圧力、
ρ
{\displaystyle \rho }
は密度、
Ω
{\displaystyle \Omega }
は外力のポテンシャル
f
=
−
∇
Ω
{\displaystyle {\boldsymbol {f}}=-\nabla \Omega }
である。
なお、
B
=
v
2
2
+
∫
d
p
ρ
+
Ω
{\displaystyle B={v^{2} \over 2}+\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}}
をベルヌーイ関数と呼ぶ。更に、右辺第2項を圧力関数 と呼ぶ。
オイラー方程式の変形の導出
非粘性流体の運動はオイラー方程式で記述される。
D
v
D
t
=
−
1
ρ
∇
p
+
f
{\displaystyle {D{\boldsymbol {v}} \over Dt}=-{1 \over \rho }\nabla p+{\boldsymbol {f}}}
ただし、
v
{\displaystyle {\boldsymbol {v}}}
は速度、
ρ
{\displaystyle \rho }
は密度、
p
{\displaystyle p}
は圧力、
f
{\displaystyle {\boldsymbol {f}}}
は外力である。
バロトロピック性
ρ
=
ρ
(
p
)
{\displaystyle \rho =\rho (p)}
と外力が保存力であることを仮定すると、
D
v
D
t
=
−
∇
{
∫
d
p
ρ
+
Ω
}
{\displaystyle {D{\boldsymbol {v}} \over Dt}=-\nabla \left\{\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}\right\}}
と書き換えられる。ただし、
Ω
{\displaystyle \Omega }
は外力のポテンシャルである。
左辺は速度の物質微分 、すなわち、加速度 であるが、加速度の回転形表示 を使うと、
D
v
D
t
=
∂
v
∂
t
+
v
⋅
∇
v
=
∂
v
∂
t
+
∇
(
v
2
2
)
−
v
×
(
∇
×
v
)
{\displaystyle {\begin{aligned}{D{\boldsymbol {v}} \over Dt}&={\partial {\boldsymbol {v}} \over \partial t}+{\boldsymbol {v}}\cdot \nabla {\boldsymbol {v}}\\&={\partial {\boldsymbol {v}} \over \partial t}+\nabla \left({v^{2} \over 2}\right)-{\boldsymbol {v}}\times \left(\nabla \times {\boldsymbol {v}}\right)\end{aligned}}}
と変形できるので、オイラー方程式は
∂
v
∂
t
−
v
×
(
∇
×
v
)
+
∇
{
v
2
2
+
∫
d
p
ρ
+
Ω
}
=
0
{\displaystyle {\partial {\boldsymbol {v}} \over \partial t}-{\boldsymbol {v}}\times \left(\nabla \times {\boldsymbol {v}}\right)+\nabla \left\{{v^{2} \over 2}+\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}\right\}=0}
となる。
これより、以下の二つの定理が導出できる。
(I) 定常流におけるベルヌーイの定理
外力が保存力 である非粘性バロトロピック流体 の定常 な流れでは、流線 と渦線 から作られるベルヌーイ面 上で
(
B
=
)
v
2
2
+
∫
d
p
ρ
+
Ω
=
c
o
n
s
t
a
n
t
{\displaystyle (B=)\,{v^{2} \over 2}+\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}=\mathrm {constant} }
が、成り立つ。
なお、簡単のため、「ベルヌーイ面上」でなく「流線 上」とすることが多い。
定常流におけるベルヌーイの定理の導出
定常流の場合、オイラー方程式の左辺第1項は消え、両辺に
v
{\displaystyle {\boldsymbol {v}}}
を内積 でかけると左辺第2項も消え、
v
⋅
∇
{
v
2
2
+
∫
d
p
ρ
+
Ω
}
=
0
{\displaystyle {\boldsymbol {v}}\cdot \nabla \left\{{v^{2} \over 2}+\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}\right\}=0}
となる。流線上の道のりを s で表すと、速度ベクトルが流線に接していることと方向微分 の考え方により、
v
⋅
∇
=
v
∂
∂
s
{\displaystyle {\boldsymbol {v}}\cdot \nabla =v{\partial \over \partial s}}
となるので、
∂
∂
s
{
v
2
2
+
∫
d
p
ρ
+
Ω
}
=
0
{\displaystyle {\partial \over \partial s}\left\{{v^{2} \over 2}+\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}\right\}=0}
すなわち流線上で
v
2
2
+
∫
d
p
ρ
+
Ω
{\displaystyle {v^{2} \over 2}+\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}}
は一定値をとる。
なお、渦度ベクトル
∇
×
v
{\displaystyle \nabla \times {\boldsymbol {v}}}
を結んで得られる渦線上でも一定値をとることも同様に示される。すなわち、流線と渦線から作られる面(ベルヌーイ面)上で
v
2
2
+
∫
d
p
ρ
+
Ω
=
c
o
n
s
t
a
n
t
{\displaystyle {v^{2} \over 2}+\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}=\mathrm {constant} }
が成り立つ。
(II) 非定常・渦なし流れにおけるベルヌーイの定理
∇
×
v
=
0
{\displaystyle \nabla \times {\boldsymbol {v}}=0}
となる流れを渦なしの流れ と呼ぶが、このとき、速度ポテンシャル と呼ばれる関数
ϕ
{\displaystyle \phi }
が存在して、
v
=
∇
ϕ
{\displaystyle {\boldsymbol {v}}=\nabla \phi }
と表せる。
渦なしの流れにおいては以下の定理(一般化されたベルヌーイの定理 )
外力が保存力 である非粘性バロトロピック流体 の渦なし の流れでは、全空間 において
∂
ϕ
∂
t
+
|
∇
ϕ
|
2
2
+
∫
d
p
ρ
+
Ω
=
f
(
t
)
{\displaystyle {\partial \phi \over \partial t}+{\left|\nabla \phi \right|^{2} \over 2}+\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}=f(t)}
(圧力方程式)が成り立つ。ただし、
f
(
t
)
{\displaystyle f(t)}
は任意の関数である。
が導ける。
(I)のタイプと違って、全空間で成り立つのが大きな特徴である。
ϕ
′
=
ϕ
−
∫
f
(
t
)
d
t
{\displaystyle \phi '=\phi -\int f(t)\mathrm {d} t}
と変更しても速度場は変わらないので、圧力方程式より
f
(
t
)
{\displaystyle f(t)}
を消去することは可能である。
渦なし流れにおけるベルヌーイの定理の導出
渦なしの流れでは
v
=
∇
ϕ
{\displaystyle {\boldsymbol {v}}=\nabla \phi }
と表せるので、オイラー方程式は
∇
{
∂
ϕ
∂
t
+
|
∇
ϕ
|
2
2
+
∫
d
p
ρ
+
Ω
}
=
0
{\displaystyle \nabla \left\{{\partial \phi \over \partial t}+{\left|\nabla \phi \right|^{2} \over 2}+\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}\right\}=0}
となり、これを積分すると
∂
ϕ
∂
t
+
|
∇
ϕ
|
2
2
+
∫
d
p
ρ
+
Ω
=
f
(
t
)
{\displaystyle {\partial \phi \over \partial t}+{\left|\nabla \phi \right|^{2} \over 2}+\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}=f(t)}
となる。ただし、
f
(
t
)
{\displaystyle f(t)}
は任意の関数である。
ベルヌーイの定理の適用条件
渦なしの流れであれば(II)のタイプ「一般化されたベルヌーイの定理」により、異なる流線間でも圧力や速さの比較ができる。非粘性流体においては上流が一様流である流れや静止状態から出発した流れは渦なし流れであるので、例えば、一様流の中の翼 の問題では異なる流線でも比較ができる。(ただし、不連続流や噴流領域を跨いではいけない。また、実在流体の場合、後述のように境界層や伴流領域は除かれる。)
上述のように「渦なしの流れ」の性質を使わなくても、一様流中の翼の問題においてベルヌーイの式を全空間で使えることを証明できる。以下簡単のため重力は無視する。
翼の近傍を通る任意の異なる2つの流線 A, B を考える。流線 A, B はともに上流の一様流まで伸びること、さらに、一様流中では速度だけでなく圧力、密度も一定、つまり、ベルヌーイ関数も一定値をとることを考慮すると、流線A上のベルヌーイ関数の値と流線B上のベルヌーイ関数の値とは等しいことが導かれる。これより、全空間でベルヌーイの式(ベルヌーイ関数の値=一定)が成立することが導かれた。
一般には、(I)のタイプの定理では異なる流線間の比較はできないが、流線曲率の定理 を使えば異なる流線間での比較ができる。流線上で成り立つベルヌーイの定理と流線曲率の定理は運動方程式の流線に関する接線 成分と主法線 成分にそれぞれ対応する。
粘性流体であっても、境界層 外部や伴流 外部の層流 領域のように、非圧縮・渦なし流れであれば粘性項の寄与を無視できるので、その領域ではベルヌーイの定理を適用可能である。
(0) 静水圧平衡
一般にベルヌーイの定理に含まれることはないが、静止流体における圧力と保存力の関係(静水圧平衡 )も運動方程式の第一積分である。
v
=
0
{\displaystyle {\boldsymbol {v}}=0}
をオイラー方程式に代入すると
∇
{
∫
d
p
ρ
+
Ω
}
=
0
∴
∫
d
p
ρ
+
Ω
=
c
o
n
s
t
a
n
t
{\displaystyle {\begin{aligned}\nabla \left\{\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}\right\}&=0\\\therefore ~\int {\mathrm {d} p \over \rho }+{\mathit {\Omega }}&=\mathrm {constant} \end{aligned}}}
が全空間で成り立つ。これより外力の等ポテンシャル面(
Ω
=
c
o
n
s
t
a
n
t
{\displaystyle \Omega =\mathrm {constant} }
)の上では
p
=
c
o
n
s
t
a
n
t
{\displaystyle p=\mathrm {constant} }
,
ρ
=
c
o
n
s
t
a
n
t
{\displaystyle \rho =\mathrm {constant} }
であることが導かれる。一様重力の等ポテンシャル面である水平面に水面(
p
{\displaystyle p}
= 大気圧 = 一定)が一致するのはこのためである[ 3] 。
一様重力のもとでの非圧縮非粘性定常流の場合
非圧縮性バロトロピック流体では密度一定だから
∫
d
p
/
ρ
=
p
/
ρ
+
c
o
n
s
t
a
n
t
{\textstyle \int dp/\rho =p/\rho +\mathrm {constant} }
とでき、一様重力のポテンシャルは
Ω
=
g
z
{\displaystyle \Omega =gz}
となるので、(I)の基本形から以下の定理が導かれる。
一様重力 のもとでの非粘性・非圧縮流体 の定常 な流れでは、流線 上で
v
2
2
+
p
ρ
+
g
z
=
c
o
n
s
t
a
n
t
{\displaystyle {v^{2} \over 2}+{p \over \rho }+gz=\mathrm {constant} }
が、成り立つ(
v
{\displaystyle v}
は速さ、
p
{\displaystyle p}
は圧力、
ρ
{\displaystyle \rho }
は密度、
g
{\displaystyle g}
は重力加速度 の大きさ、
z
{\displaystyle z}
は鉛直方向の座標である)。
エネルギー保存則に基づいた導出
右方向に流れる流体の作る流管 。
非粘性・非圧縮性・定常流におけるベルヌーイの定理は、体積の保存則(質量保存則 )、および、仕事 とエネルギーの関係(エネルギー保存則)から導くことができる。
初期に断面A 1 とA 2 に挟まれた流体に着目する。断面A 1 から流入した流体粒子 は時間Δt の間に、s 1 =v 1 Δt だけ移動し、
断面A 2 から流出した流体粒子はs 2 =v 2 Δt だけ移動する。
(1)体積の保存。断面 A 1 から流入した体積と断面 A 2 から流出した体積はそれぞれ A 1 s 1 と A 2 s 2 となり、定常な非圧縮性流体を考えているので、
A
1
s
1
=
A
2
s
2
(
=
Δ
V
)
{\displaystyle A_{1}s_{1}=A_{2}s_{2}(=\Delta V)}
が成り立つ。
(2) 系の力学的エネルギーの増分は系になされた仕事に等しい。
(2-1) 接触力(圧力由来)は、断面 A 1 では正の向きに、断面 A 2 では負の向きに、挟まれた流体に対して仕事をするので、
W
=
(
p
1
A
1
)
s
1
−
(
p
2
A
2
)
s
2
=
Δ
V
(
p
1
−
p
2
)
{\displaystyle W=(p_{1}A_{1})s_{1}-(p_{2}A_{2})s_{2}=\Delta V(p_{1}-p_{2})}
(2-2) 重力の位置エネルギー U の変化は、高さ z 1 にある質量 ρ ΔV の流体が、高さ z 2 に移動したと考えれば、
Δ
U
=
ρ
Δ
V
(
g
z
2
−
g
z
1
)
{\displaystyle \Delta U=\rho \Delta V\,(gz_{2}-gz_{1})}
(2-3) そして、運動エネルギー K の変化は、速度 v 1 である質量 ρ ΔV の流体が、速度 v 2 になると考えれば、
Δ
K
=
ρ
Δ
V
(
v
2
2
2
−
v
1
2
2
)
{\displaystyle \Delta K=\rho \Delta V\left({v_{2}^{2} \over 2}-{v_{1}^{2} \over 2}\right)}
(2-4) エネルギー保存則
Δ
K
+
Δ
U
=
W
{\displaystyle \Delta K+\Delta U=W}
より、
ρ
Δ
V
(
v
2
2
2
−
v
1
2
2
)
+
ρ
Δ
V
(
g
z
2
−
g
z
1
)
=
ρ
Δ
V
(
p
1
ρ
−
p
2
ρ
)
∴
v
2
2
2
+
p
2
ρ
+
g
z
2
=
v
1
2
2
+
p
1
ρ
+
g
z
1
{\displaystyle {\begin{aligned}&\rho \Delta V\left({v_{2}^{2} \over 2}-{v_{1}^{2} \over 2}\right)+\rho \Delta V\,(gz_{2}-gz_{1})=\rho \Delta V\left({p_{1} \over \rho }-{p_{2} \over \rho }\right)\\\therefore ~&{v_{2}^{2} \over 2}+{p_{2} \over \rho }+gz_{2}={v_{1}^{2} \over 2}+{p_{1} \over \rho }+gz_{1}\end{aligned}}}
が得られる。
(3) これは流管内の任意の断面で成り立つものであり、断面積を小さくとると流線上の任意の点で成り立つと考えてよい。
これより、流線上で
v
2
2
+
p
ρ
+
g
z
=
c
o
n
s
t
a
t
n
t
{\displaystyle {\frac {v^{2}}{2}}+{\frac {p}{\rho }}+gz=\mathrm {constatnt} }
が成り立つことが導けた。
表現の違い
水頭 による表現
v
2
2
g
+
p
ρ
g
+
z
=
H
0
=
c
o
n
s
t
a
n
t
{\displaystyle {\frac {v^{2}}{2g}}+{\frac {p}{\rho g}}+z=H_{0}=\mathrm {constant} }
[m ]
第一項
v
2
/
(
2
g
)
{\displaystyle v^{2}/(2g)}
を速度ヘッド (英 : velocity head )という。第二項
p
/
ρ
g
{\displaystyle p/\rho g}
を 圧力ヘッド (英 : pressure head )という。第三項
z
{\displaystyle z}
を位置ヘッド (英 : potential head )という。これら全てを足しあわせた値
H
0
{\displaystyle H_{0}}
を全ヘッド (英 : total head )という。
エネルギーによる表現
1
2
m
0
v
2
+
p
m
0
ρ
+
m
0
g
z
=
c
o
n
s
t
a
n
t
{\displaystyle {\frac {1}{2}}m_{0}v^{2}+p{\frac {m_{0}}{\rho }}+m_{0}gz=\mathrm {constant} }
[J ]
圧力による表現
1
2
ρ
v
2
+
p
+
ρ
g
z
=
c
o
n
s
t
a
n
t
{\displaystyle {\frac {1}{2}}\rho v^{2}+p+\rho gz=\mathrm {constant} }
[Pa ]
単純形
位置エネルギー の変化が無視できる場合、
非粘性・非圧縮流 の定常 な流れでは、流線 上で
1
2
ρ
v
2
+
p
=
p
0
(
=
c
o
n
s
t
a
n
t
)
{\displaystyle {\frac {1}{2}}\rho v^{2}+p=p_{0}(=\mathrm {constant} )}
が、成り立つ(
v
{\displaystyle v}
は速さ、
p
{\displaystyle p}
は圧力、
ρ
{\displaystyle \rho }
は密度)。
となる。なお、非圧縮流とは非圧縮性流体(液体)のことではなく低マッハ数 の流れを指す。
左辺第一項を動圧 、第二項を静圧 、右辺の値を総圧 という。
静圧 (static pressure):
p
{\displaystyle p}
流体が実際に外界に及ぼす圧力。
動圧 (dynamic pressure):
q
=
1
2
ρ
v
2
{\displaystyle q={\frac {1}{2}}\rho v^{2}}
動圧 は流体要素の運動エネルギーに相当する量であり、次元 が圧力に一致するものの、流体要素が速度を保つ限りは周囲の流体要素を押すような効果はない。仮想的には流体要素を静止させられればその瞬間に生じる圧力であるが実際測定 はできない。よどみ点圧(=総圧)と静圧の差や、密度と流速から算出される。
総圧(total pressure):
p
0
=
p
+
q
{\displaystyle p_{0}=p+q}
総圧 は動圧と静圧の和。よどみ点以外では総圧を直接測定することはできない。全圧 ともよぶが、「全圧」は分圧 に対しても使われる。
流速が増すと動圧は増すが、上記条件の総圧が一定の系 では、そのぶん静圧が減る。
なお、「総圧」も「動圧」もベルヌーイ式の保存性を説明するために使われる言葉で圧力としてはそれ以上の意味はない。これらと区別するために付けられた「静圧」も「圧力」以上の意味は無い。
よどみ点圧
相対的な流れの中の物体表面で流速が0になる点(よどみ点 )での圧を、よどみ点圧 と呼ぶ。よどみ点では動圧が0なので、よどみ点圧は静圧であり総圧でもある。
直感的解釈
ベルヌーイの定理は非粘性流体の支配方程式であるオイラー方程式から直接導出できるが、ベルヌーイの定理(I)の物理的解釈は流体粒子 に対する力と加速度 の関係(ニュートンの運動の第2法則)で以下のように解釈が可能である。
流体粒子が圧力の高い領域から低い領域へと水平に流れていくとき、流体粒子が後方から受ける圧力は前方から受ける圧力より大きい。よって流体粒子全体には流線に沿って前方へと加速する力が働く。つまり、粒子の速さは移動につれて大きくなる[ 4] 。
よって流線上で、相対的に圧力が低い所では相対的に運動エネルギー が大きく、相対的に圧力が高い所では相対的に運動エネルギーが小さい。これは粒子の位置エネルギーと運動エネルギーの関係に相当する。
圧縮性流体
ベルヌーイは液体の実験より法則を導き出したので、彼のオリジナルの理論は非圧縮性流体あるいはマッハ数の小さな圧縮性流体にしか適用できない。しかし、バロトロピック性を仮定すれば、一般の圧縮性流体に対しても適用可能である。
特に気体の場合、断熱過程として適用することが多い。以下、タイプ(I)のベルヌーイの定理の応用例について解説する。タイプ(II)のベルヌーイの定理の応用例については文献を参照[ 3] 。
気体の定常流の場合
気体の運動では、重力が無視でき、また、運動の時間スケールが熱伝導 の時間スケールに比べて十分小さく断熱過程 と見なせる場合が多い。このときポアソンの法則 により
p
∝
ρ
γ
{\displaystyle p\propto \rho ^{\gamma }}
(
γ
=
C
p
/
C
v
{\displaystyle \gamma =C_{p}/C_{v}}
は比熱比 )とできるので、ベルヌーイの定理は
断熱過程 に従う非粘性気体 の定常 な流れでは、流線 上で
v
2
2
+
(
γ
γ
−
1
)
p
ρ
=
(
γ
γ
−
1
)
p
s
ρ
s
{\displaystyle {\frac {v^{2}}{2}}+\left({\frac {\gamma }{\gamma -1}}\right){\frac {p}{\rho }}=\left({\frac {\gamma }{\gamma -1}}\right){\frac {p_{s}}{\rho _{s}}}}
が、成り立つ。ただし、
v
{\displaystyle {\boldsymbol {v}}}
は速度ベクトル、
p
{\displaystyle p}
は圧力、
ρ
{\displaystyle \rho }
は密度、
γ
=
C
p
/
C
v
{\displaystyle \gamma =C_{p}/C_{v}}
は比熱比、
p
s
,
ρ
s
{\displaystyle p_{s},\rho _{s}}
はよどみ点 における圧力と密度である。
と書き換えられる。
音速 の定義
a
=
γ
(
p
/
ρ
)
{\displaystyle a={\sqrt {\gamma \left({p/\rho }\right)}}}
を用いれば、ベルヌーイの定理は
v
=
2
γ
−
1
(
a
s
2
−
a
2
)
{\displaystyle v={\sqrt {{\frac {2}{\gamma -1}}(a_{s}^{2}-a^{2})}}}
が流線上で成り立つ。(as はよどみ点における音速)
となる。
真空 では a =0 となるので、そのとき流速は最大値
v
m
a
x
=
2
γ
−
1
a
s
{\displaystyle v_{max}={\sqrt {\frac {2}{\gamma -1}}}a_{s}}
に到達する。例えば、大きな容器に封入された気体が器壁の小さな孔から真空中に噴出する場合の流速がそれにあたる。容器の中が1気圧、15℃の空気の場合、
γ
=
1.4
,
a
s
=
340
m
/
s
{\displaystyle \gamma =1.4,a_{s}=340m/s}
であるから、
v
m
a
x
=
760
m
/
s
{\displaystyle v_{max}=760m/s}
となる[ 3] 。
ベルヌーイの定理と流線曲率の定理
ベルヌーイの定理(I)と流線曲率の定理 とは運動方程式の流線に関する接線 成分と主法線 成分に対応する。
外力が無視できる非粘性バロトロピック流体の定常な流れの運動方程式
v
⋅
∇
v
=
−
∇
∫
d
p
ρ
{\displaystyle {\boldsymbol {v}}\cdot \nabla {\boldsymbol {v}}=-\nabla \int {\mathrm {d} p \over \rho }}
の接線成分と主法線成分は、定常流における加速度の分解 により、
∂
∂
s
(
v
2
2
+
∫
d
p
ρ
)
=
0
∂
p
∂
r
=
ρ
v
2
r
(
o
r
∂
∂
r
∫
d
p
ρ
=
v
2
r
)
{\displaystyle {\begin{aligned}&{\partial \over \partial s}\left({v^{2} \over 2}+\int {\mathrm {d} p \over \rho }\right)=0\\&{\partial p \over \partial r}=\rho {v^{2} \over r}\quad \left(\mathrm {or~} {\partial \over \partial r}\int {\mathrm {d} p \over \rho }={v^{2} \over r}\right)\end{aligned}}}
となる。
ただし、s は流線上の道のり(接線方向)、r は流線を円弧と近似したときの中心(曲率中心 )からの距離(主法線方向)を表す。
第1式がベルヌーイの定理、第2式が流線曲率の定理に対応する。
一般には、(I)のタイプのベルヌーイの定理では異なる流線間の比較はできないが、流線曲率の定理を使えば異なる流線間での比較ができる。
揚力とベルヌーイの定理
ベルヌーイの定理は十分に検証された理論である[ 5] [ 6] [ 7] 。翼の周りの流体の速度分布 が正しくわかれば、翼に発生する揚力 の大きさをベルヌーイの定理を使って十分に良い精度で計算できる 。しかし、ベルヌーイの定理では翼の形から流体の速度分布を求めることはできないので、翼の周りの流体の速度分布を説明する理論は別途必要である。その理論について誤解がある。
同着の原理にまつわる誤解
揚力についての一般向けの解説には、
「同着の原理」のため翼の上の流れが下の流れより速くなり、ベルヌーイの定理により翼の上の圧力が下の圧力より小さくなり、よって上向きの揚力が発生する
と説明しているものがある[ 8] 。
「同着の原理」とは、「翼の前縁で上下に別れた流体は翼の後縁に同着する。」という原理である。この原理により、翼の上の経路長が下の経路長より長い場合、「翼の上を流れる速さが下の速さより大きくなる」という翼の周りの流体の速度分布が「導かれる」。しかし、実際には、上面の流れの方が後縁により早く到着し、同着の原理は成り立たない [ 9] 。
現在、「同着の原理」が間違いであることは広く知られるようになった。しかし、「ベルヌーイの定理を揚力の説明に使うのは誤りで、流線曲率の定理やニュートンの運動方程式を使うべきだ」という誤解も見られるようになった。
一般向けの説明で誤っているのは「同着の原理」のみであり、「同着の原理」はベルヌーイの定理とは無関係である。むしろ、同着原理の不成立に導いた、上面の流れの方が後端により早く到着するという実験事実は、ベルヌーイの定理による揚力の発生を「補強こそすれ、否定的な意見とはならない」[ 10] 。
また、ベルヌーイの定理が間違いで流線曲率の定理やニュートンの運動方程式が正しいというのは矛盾を含む。翼の周りの流体の速度分布が正しくわかれば 、ベルヌーイの定理でも、流線曲率定理でも、運動量変化と力積(あるい反作用)の関係でも、正しく適用する限り、同じ結果が得られる 。なぜなら、これらはいずれもニュートン力学に起源を持つ理論だからである [ 11] [ 12] 。
コアンダ効果とクッタの条件
翼の周りの速度分布を説明する理論としては、「同着の原理」のほかに、コアンダ効果 [ 13] [ 14] やクッタの条件 などがある。
コアンダ効果は「粘性の効果によって翼の形に沿うように流れる」というもので、これと作用反作用則を使った揚力の原理の説明はベルヌーイの定理を使わない説明として知られている。ただし、コアンダ効果は本来噴流 (ジェット)が物体に沿う性質であり、通常の翼には噴流は発生しないので、コアンダ効果を通常の翼の速度分布の説明に使うのは不適切であるとの意見もある[ 15] ので、コアンダ効果を使った揚力の説明には疑問がある。
クッタの条件は「粘性の効果によって翼の後端のエッジにおいて気流が翼から離れる」というものである。適切な形の翼に対して、クッタの条件に基づき循環 量を決定(= 速度分布を決定)し、クッタ・ジュコーフスキーの定理 を使って循環量と速度から計算した揚力が、実験ともよく合うことが知られている[ 16] 。なお、「クッタ・ジュコーフスキーの定理」の導出にはベルヌーイの定理 が使われている。
ベルヌーイの定理の拡張
系外からのエネルギーのやりとりを考えた拡張されたベルヌーイの定理も存在する[要出典 ] 。
脚注
出典
^ 日野幹雄『流体力学』朝倉書店、1992年。ISBN 4254200668 。
^ “ベルヌーイの定理:楽しい流れの実験教室 ”. 日本機械学会流体工学部門:楽しい流れの実験教室 . 2025年5月14日閲覧。
^ a b c d
巽友正『流体力学』培風館、1982年。ISBN 456302421X 。
^ Babinsky, Holger (November 2003). “How do wings work?” (PDF). Physics Education 38 (6): 497. doi :10.1088/0031-9120/38/6/001 . http://www3.eng.cam.ac.uk/outreach/Project-resources/Wind-turbine/howwingswork.pdf .
^
Batchelor, G.K. (1967). An Introduction to Fluid Dynamics . Cambridge University Press. ISBN 0-521-66396-2
Sections 3.5 and 5.1
^
Lamb, H. (1993). Hydrodynamics (6th ed.). Cambridge University Press. ISBN 978-0-521-45868-9
§17–§29
^
ランダウ&リフシッツ『流体力学』東京図書、1970年。ISBN 4489011660 。
^ 飛行機はなぜ飛ぶかのかまだ分からない?? - NPO法人 知的人材ネットワーク・あいんしゅたいん - 松田卓也 による解説。
^
Glenn Research Center (2006年3月15日). “Incorrect Lift Theory ”. NASA. 2012年4月20日閲覧。
^
早川尚男. “飛行機の飛ぶ訳 (流体力学の話in物理学概論) ”. 京都大学OCW. 2013年4月8日閲覧。
^
“Newton vs Bernoulli ”. NASA. 2012年4月20日閲覧。
^
Ison, David. Bernoulli Or Newton: Who's Right About Lift? Retrieved on 2009-11-26
^
David Anderson; Scott Eberhardt,. "Understanding Flight, Second Edition" (2 edition (August 12, 2009) ed.). ,McGraw-Hill Professional. ISBN 0071626964
^
日本機械学会『流れの不思議』(2004年8月20日第一刷発行)講談社ブルーバックス。ISBN 4062574527 。
^ Report on the Coandă Effect and lift , オリジナル の2011年7月14日時点におけるアーカイブ。, https://web.archive.org/web/20110714172646/http://newfluidtechnology.com/THE_COANDA_EFFECT_AND_LIFT.pdf
^
Kundu, P.K. (2011). Fluid Mechanics Fifth Edition . Academic Press. ISBN 0123821002
関連項目