この節の作者: 清水康行 <yasu@i-ric.org>

Appendix V (Engelundによる主流および2次流の流速分布式)

注釈

下記の式展開はEngelund [Ref:38] による原著論文を北海道大学 岩崎理樹准教授 が式展開・整理して頂いたものを ご厚意で提供頂き,本テキスト用に再整理したものです. 2次元河床変動計算で良く用いられる「2次流強度=7」の根拠が こちら になります.

主流の鉛直方向流速分布

水深平均流の流れ方向を \(s\) ,鉛直上向き方向を \(z\) とすると, \(s\) 方向の等流流れの運動方程式は次式で表されます.

(453)\[0=-g{{\partial H}\over{\partial s}}+{\partial \over{\partial z}} \left(\nu_t{{\partial u_s}\over{\partial z}}\right)\]

ここで,\(g\) は重力加速度,\(H\) は水位,\(s\) は主流の流下方向, \(u_s\)\(s\) 方向の流速,\(z\) は鉛直方向座標軸, \(\nu_t\) は渦動粘性係数です. 鉛直方向距離を水深 \(h\) で無次元化し,河床高を \(z_b\) ,無次元鉛直方向距離を \(\zeta\) (河床で0, 水面で1)とし,次式で表します.

(454)\[\zeta=\cfrac{z-z_b}{h}\]

等流状態を仮定し,エネルギー勾配(=水面勾配)を \(I_e\) とすると, \(I_e\) は次式で 表されます.

(455)\[I_e=-{{\partial H}\over{\partial s}}\]

主流の水深平均流速を \(<u_s>\) とし, 次式で表します.

(456)\[u_s(\zeta)=<u_s>f_s(\zeta)\]

ただし, f_s は これを運動方程式 (453) に代入します.

(457)\[\cfrac{\partial^2 f_s}{\partial \zeta^2}=-\cfrac{gI_eh^2}{\nu_t <u_s>}\]

\(\zeta\) で積分すると,

(458)\[{{\partial f_s}\over{\partial \zeta}}=-\cfrac{gI_eh^2}{\nu_t<u_s>}\zeta+C_1\]

となります. ただし, \(C_1\) は積分定数です. 水面でせん断力がゼロなので,\(\zeta=1\)\(\cfrac{\partial f_s}{\partial \zeta}=0\) より \(C_1\) は次式となり, これを \(\beta\) と置くと,

(459)\[C_1=\cfrac{gI_eh^2}{\nu_t<u_s>} \equiv \beta\]

となり, (458) 式は次式となります.

(460)\[{{\partial f_s}\over{\partial \zeta}}=\beta(1-\zeta)\]

これをもう一度 \(\zeta\) で積分すると, \(C_2\) を積分定数として \(f_s\)

(461)\[f_s=\beta\left(\zeta-\cfrac{1}{2}\zeta^2\right)+C_2\]

となり, 平均流の定義より,

(462)\[\int_0^1 f_s d\zeta =1= \left[ \beta\left(\cfrac{1}{2}\zeta^2-\cfrac{1}{6}\zeta^3\right)+C_2\zeta \right]^1_0 =\cfrac{1}{3}\beta+C_2\]

なので,

(463)\[C_2=1-\cfrac{1}{3}\beta\]

となります. これを (462) 式に戻すと次式が得られます.

(464)\[f_s=\left(-{1 \over 2}\zeta^2+\zeta-{1 \over 3}\right)\beta+1\]

渦動粘性係数を \(\nu_t=\alpha u_\ast h\) で表し, 摩擦速度 \(u_\ast\) が, \(u_\ast=\sqrt {ghI_e}\) であることを考慮すると,

(465)\[\beta=\cfrac{gI_eh^2}{\nu_t<u_s>}=\cfrac{gI_eh}{\alpha u_\ast <u_s>}=\cfrac{u_\ast}{\alpha<u_s>}\]

となり, 底面での流速(Slip Velocity)を \(u_s^b\) とすると

(466)\[u_s^b=<u_s>f_s(0)=-\cfrac{u_\ast}{3\alpha}+<u_s>\]

より,

(467)\[\cfrac{<u_s>}{u_\ast} = \cfrac{u_s^b}{u_\ast}+\cfrac{1}{3\alpha}\]
(468)\[\cfrac{u_s^b}{u_\ast}=2+{{1}\over{\kappa}} \ln {{h}\over{k_s}} =r_\ast\]

とおくと,

(469)\[u_\ast^2h(1-\xi)=\alpha u_* h{{\partial u}\over{\partial \xi}}\]

したがって,

(470)\[\cfrac{<u_s>}{u_\ast} = r_\ast+\cfrac{1}{3\alpha}\]
(471)\[\cfrac{u_\ast}{<u_s>}= \cfrac{1}{r_\ast + \cfrac{1}{3\alpha}} =\cfrac{3\alpha}{3\alpha r_\ast +1}\]

これを (464) 式に代入すると

(472)\[\beta=\left(\cfrac{3\alpha}{3\alpha r_\ast +1}\right)\cfrac{1}{\alpha} =\cfrac{1}{\alpha r_\ast+\cfrac{1}{3}}\]
(473)\[\cfrac{1}{\beta} = \alpha r_\ast + \cfrac{1}{3}\]

ゆえに,

(474)\[f_s=\left( -\cfrac{1}{2}\zeta^2+\zeta-\cfrac{1}{3}\right)\beta+1 =\left( -\cfrac{1}{2}\zeta^2+\zeta-\cfrac{1}{3}+\cfrac{1}{\beta} \right)\beta =\cfrac{\alpha r_\ast + \zeta - \cfrac{1}{2} \zeta^2}{\alpha r_\ast + \cfrac{1}{3}}\]

ここで, \(r_\ast \alpha=\chi\) , \(\chi_1=\alpha r_\ast +\cfrac{1}{3}\) とすると,

(475)\[f_s=\cfrac{\chi+\zeta-\cfrac{\zeta^2}{2}}{\chi_1}\]

したがって, \(u_s(\zeta)\) は次式となります.

(476)\[u_s(\zeta)=<u_s>\cfrac{\chi+\zeta-\cfrac{\zeta^2}{2}}{\chi_1}\]

(476) 式が 放物線分布 と呼ばれる流線方向の流速鉛直分布式です.

2次流の流速分布

_images/s-n-coord.png

Figure 156 : 座標系の定義

Figure 156 に示すように, 主流方向の流れが曲率半径 \(r_s\) でカーブしている場合(主流方向を \(s\) 軸とし, \(s\) 軸に向かって右にカーブしている場合を \(r_s>0\) とする, Figure 156 では \(r_s<0\) となる), \(n\) 軸を \(s\) 軸に直交する方向で \(s\) が増加する方向に向かって左側が正になる方向に取り,\(n\) 方向の流速を \(u_n\) とすると,\(n\) 軸方向の等流流れの運動方程式は次式で表されます.

(477)\[{{u_s^2}\over{r_s}} = -g{{\partial H}\over{\partial n}} +{{\partial}\over{\partial z}} \left( \nu_t {{\partial u_n}\over{\partial z}} \right)\]

ただし,

(478)\[\cfrac{1}{r_s}= \cfrac{\partial \theta_s}{\partial s}\]

ここで,\(\theta\)\(x\) 軸と主流方向 \(s\) とのなす角度です. \(n\) 方向の流速に以下の分布を仮定します.

(479)\[u_n(\zeta)=A_n f_n(\zeta)\]

ただし, \(A_n\) は2次流強度を表す係数,\(f_n\) は無次元流速分布関数です. これを (477) に代入して整理すると,

(480)\[\cfrac{\partial^2 f_n}{\partial \zeta^2}=\cfrac{g h^2}{\nu_t A_n} \cfrac{\partial H}{\partial n} +\cfrac{<u_s>^2 h^2}{\nu_t A_n r_s} f_s^2 =\cfrac{<u_s>^2h^2}{\nu_t A_n r_s} \left( \cfrac{g r_s}{<u_s>^2}\cfrac{\partial H}{\partial n} +f_s^2 \right)\]

ここで,

(481)\[\cfrac{<u_s>^2h^2}{\nu_t A_n r_s} \equiv A\]
(482)\[\cfrac{g r_s}{<u_s>^2}\cfrac{\partial H}{\partial n} \equiv B\]

とおくと, (480) 式は,

(483)\[\cfrac{\partial^2 f_n}{\partial \zeta^2}=A(B+f_s^2)\]

と書けます,これを, \(\zeta\) で積分しますと,

(484)\[\begin{split}\cfrac{\partial f_n}{\partial \zeta} =AB\zeta+ A\int \left[\left\{{ {1}\over{\chi_1}} \left( \chi + \zeta -{{1}\over{2}} \zeta^2 \right)\right\}^2 \right] d\zeta + C_3 \\ =AB \zeta + \cfrac{A}{{\chi_1}^2} \left[ \chi^2 \zeta + \chi \zeta^2 +{{1}\over{3}}(1-\chi)\zeta^3-{{1}\over{4}}\zeta^4+ {{1}\over{20}}\zeta^5 \right] +C_3\end{split}\]

となる.ただし \(C_3\) は積分定数です. 水面では \(\cfrac{\partial f_n}{\partial \zeta}=0\) (Slip条件)なので,

(485)\[C_3=-AB - \cfrac{A}{{\chi_1}^2} \left(\chi^2+{{2}\over{3}}\chi+{{2}\over{15}} \right)\]

よって,

(486)\[\begin{split}\cfrac{\partial f_n}{\partial \zeta} =AB \zeta + \cfrac{A}{{\chi_1}^2} \left[ \chi^2 \zeta + \chi \zeta^2 +{{1}\over{3}}(1-\chi)\zeta^3-{{1}\over{4}}\zeta^4+ {{1}\over{20}}\zeta^5 \right] \\ -\left[ AB + \cfrac{A}{{\chi_1}^2} \left(\chi^2+{{2}\over{3}}\chi+{{2}\over{15}} \right) \right]\end{split}\]

これをもう1度 \(\zeta\) に関して積分します.

(487)\[\begin{split}f_n={{1}\over{2}} AB\zeta^2 +\cfrac{A}{{\chi_1}^2} \left[ {{1}\over{2}}\chi^2 \zeta^2 + {{1}\over{3}}\chi \zeta^3 +{{1}\over{12}}(1-\chi)\zeta^4-{{1}\over{20}}\zeta^5+ {{1}\over{120}}\zeta^6 \right] \\ -\left[ AB + \cfrac{A}{{\chi_1}^2} \left(\chi^2+{{2}\over{3}}\chi+{{2}\over{15}} \right) \right]\zeta + C_4\end{split}\]

ただし, \(C_4\) は積分定数です. 2次流はその定義から水深方向積分値はゼロになります. すなわち \(\displaystyle{\int_0^1 f_n d\zeta =0}\) なので,

(488)\[\begin{split}\int_0^1 f_n d \zeta = {{1}\over{6}} AB\zeta^3 +\cfrac{A}{{\chi_1}^2} \left[ {{1}\over{6}}\chi^2 \zeta^3 + {{1}\over{12}}\chi \zeta^4 +{{1}\over{60}}(1-\chi)\zeta^5-{{1}\over{120}}\zeta^6+ {{1}\over{840}}\zeta^7 \right] \\ -\left[ AB + \cfrac{A}{{\chi_1}^2} \left(\chi^2+{{2}\over{3}}\chi+{{2}\over{15}} \right) \right]\cfrac{\zeta^2}{2} + C_4\zeta =0\end{split}\]

これより,

(489)\[C_4={{1}\over{3}}AB + \cfrac{A}{{\chi_1}^2} \left[ {{1}\over{3}} \chi^2 +{{4}\over{15}}\chi +{{2}\over{35}} \right]\]

これを (487) 式に代入しますと,

(490)\[\begin{split}f_n={{A}\over{2}} \left(B+\cfrac{\chi^2}{{\chi_1}^2}\right)\zeta^2 +\cfrac{A}{{\chi_1}^2} \left[ {{1}\over{3}}\chi \zeta^3 +{{1}\over{12}}(1-\chi)\zeta^4-{{1}\over{20}}\zeta^5+ {{1}\over{120}}\zeta^6 \right] \\ -\left[ AB + \cfrac{A}{{\chi_1}^2} \left(\chi^2+{{2}\over{3}}\chi+{{2}\over{15}} \right) \right]\zeta +{{1}\over{3}}AB+\cfrac{A}{{\chi_1}^2} \left({{1}\over{3}}\chi^2+{{4}\over{15}}\chi+{{2}\over{35}}\right)\end{split}\]

一方, 底面流速ベクトルの方向とせん断力ベクトルの方向は同ですので,

(491)\[\cfrac{u_n^b}{u_s^b} = \cfrac{\tau_n^b}{\tau_s^b}\]

ただし, \(u_s^b, u_n^b, \tau_s^b, \tau_n^b\) はそれぞれ, \(s\) 方向よよび, \(n\) 方向の底面流速および河床せん断力です. (491) 式中の各値は下記 (492) \(\sim\) (495) です.

(492)\[u_s^b=<u_s>f_s(0)=<u_s> \cfrac{\chi}{\chi_1}\]
(493)\[u_n^b=A_n f_n(0) = A A_n \left[{{1}\over{3}}B + {{1}\over{{\chi_1}^2}} \left( {{1}\over{3}}\chi^2+{{4}\over{15}}\chi+{{2}\over{35}} \right) \right]\]
(494)\[\cfrac{\tau_s^b}{\rho} = u_\ast^2\]
(495)\[\cfrac{\tau_n^b}{\rho} = \nu_t \left. \cfrac{\partial u_n}{\partial z} \right|_{z=0} =\nu_t \cfrac{A_n}{h} \left. \cfrac{\partial f_n}{\partial \zeta} \right|_{\zeta=0} =-\alpha u_\ast A_n A \left[ B+{{1}\over{{\chi_1}^2}} \left(\chi^2+{{2}\over{3}}\chi +{{2}\over{15}} \right) \right]\]

これらを (491) 式に代入します.

(496)\[\begin{split}\cfrac{ A A_n \left[\cfrac{1}{3}B + \cfrac{1}{{\chi_1}^2} \left( \cfrac{1}{3}\chi^2+\cfrac{4}{15}\chi+\cfrac{2}{35} \right) \right]} {<u_s> \cfrac{\chi}{\chi_1}} \\ =-\cfrac{\alpha u_\ast A_n A \left[ B+\cfrac{1}{{\chi_1}^2} \left(\chi^2+\cfrac{2}{3}\chi +\cfrac{2}{15} \right) \right]} {u_\ast^2} \\ \cfrac{1}{3}B+\cfrac{1}{{\chi_1}^2} \left( \cfrac{1}{3}\chi^2+\cfrac{4}{15}\chi+\cfrac{2}{35} \right) \\ =-\alpha \cfrac{<u_s>}{u_\ast} \cfrac{\chi}{\chi_1} \left[ B+\cfrac{1}{{\chi_1}^2} \left(\chi^2+\cfrac{2}{3}\chi +\cfrac{2}{15} \right) \right]\end{split}\]

ここで,

(497)\[\cfrac{<u_s>}{u_\ast} = r_\ast +\cfrac{1}{3\alpha}=\cfrac{\chi_1}{\alpha}, \; \; \chi_1=\chi+\cfrac{1}{3}\]

の関係を用いると,

(498)\[\cfrac{1}{3}B+\cfrac{1}{{\chi_1}^2} \left( \cfrac{1}{3}\chi^2+\cfrac{4}{15}\chi+\cfrac{2}{35} \right) =-\chi \left[ B+\cfrac{1}{{\chi_1}^2} \left(\chi^2+\cfrac{2}{3}\chi +\cfrac{2}{15} \right) \right]\]
(499)\[\left(\chi+\cfrac{1}{3}\right)B=-\cfrac{1}{{\chi_1}^2} \left(\chi^3+\chi^2+\cfrac{2}{5}\chi+\cfrac{2}{35}\right)\]
(500)\[B=-\cfrac{1}{{\chi_1}^3}\left(\chi^3+\chi^2+\cfrac{2}{5}\chi+\cfrac{2}{35}\right)\]

となります. これより \(f_n\)

(501)\[\begin{split}\cfrac{f_n}{A}=\cfrac{1}{2}\left(B+\cfrac{\chi^2}{{\chi_1}^2}\right)\zeta^2 +\cfrac{1}{{\chi_1}^2} \left[ \cfrac{1}{3}\chi\zeta^3+\cfrac{1}{12}(1-\chi)\zeta^4-\cfrac{1}{20}\zeta^5+\cfrac{1}{120}\zeta^6 \right] \\ -\left[B+\cfrac{1}{{\chi_1}^2}\left(\chi^2+\cfrac{2}{3}\chi+\cfrac{2}{15}\right)\right]\zeta +\left[\cfrac{1}{3}B+\cfrac{1}{{\chi_1}^2}\left(\cfrac{1}{3}\chi^2+\cfrac{4}{15}\chi+\cfrac{2}{35}\right)\right]\end{split}\]

となります. ここで, 右辺最終項は (498) 式の関係を用いると,

(502)\[\begin{split}\cfrac{f_n}{A}=\cfrac{1}{2}\left(B+\cfrac{\chi^2}{{\chi_1}^2}\right)\zeta^2 +\cfrac{1}{{\chi_1}^2} \left[ \cfrac{1}{3}\chi\zeta^3+\cfrac{1}{12}(1-\chi)\zeta^4-\cfrac{1}{20}\zeta^5+\cfrac{1}{120}\zeta^6 \right] \\ -\left[B+\cfrac{1}{{\chi_1}^2}\left(\chi^2+\cfrac{2}{3}\chi+\cfrac{2}{15}\right)\right]\zeta -\chi \left[ B+\cfrac{1}{{\chi_1}^2} \left(\chi^2+\cfrac{2}{3}\chi +\cfrac{2}{15} \right) \right] \\ =\cfrac{1}{{\chi_1}^2} \left[ -\left(\chi^2+\cfrac{2}{3}\chi+\cfrac{2}{15}\right) \left(\zeta+\chi\right) +\cfrac{1}{2}\chi^2 \zeta^2+\cfrac{1}{3}\chi\zeta^3 \right. \\ +\left. \cfrac{1}{12}\left(1-\chi \right)\zeta^4-\cfrac{1}{20}\zeta^5+\cfrac{1}{120}\zeta^6 \right]+B\left(\cfrac{1}{2}\zeta^2-\zeta-\chi\right)\end{split}\]

となります. ここで,

(503)\[\begin{split}A=\cfrac{<u_s>^2 h^2}{\nu_t A_n r_s}= \cfrac{1}{A_n} \cfrac{<u_s>^2 h^2}{\alpha u_\ast h r_s} \\ =\cfrac{1}{A_n} \cfrac{1}{\alpha} \cfrac{<u_s>}{u_\ast} <u_s> \cfrac{h}{r_s} =\cfrac{1}{A_n} \cfrac{1}{C_f \chi_1} <u_s> \cfrac{h}{r_s}\end{split}\]

なので,2次流強度 \(A_n\)

(504)\[A_n = <u_s>\cfrac{h}{r_s}\]

と定義すると,最終的に2次流分布は下記のようになります.

(505)\[u_n=A_n f_n, \; \; f_n=\cfrac{G_0(\zeta)}{C_f \chi_1}\]

ただし,

(506)\[\begin{split}G_0(\zeta)=\cfrac{1}{{\chi_1}^2}\left[ -\left(\chi^2+\cfrac{2}{3}\chi+\cfrac{2}{15}\right)\left(\zeta+\chi\right)+\cfrac{1}{2}\chi^2\zeta^2+\cfrac{1}{3}\chi\zeta^3 \right. \\ \left. +\cfrac{1}{12}\left(1-\chi\right)\zeta^4 -\cfrac{1}{20}\zeta^5+\cfrac{1}{120}\zeta^6 \right]+\chi_{20}\left(\cfrac{1}{2}\zeta^2-\zeta-\chi\right)\end{split}\]
(507)\[\chi_{20} = B = -\cfrac{1}{{\chi_1}^3} \left(\chi^3+\chi^2+\cfrac{2}{5}\chi+\cfrac{2}{35}\right), \: \: \cfrac{<u_s>}{u_\ast} = \cfrac{1}{\sqrt{C_f}}, \; \chi=\chi_1-\cfrac{1}{3}\]

これがEngelundによる2次流の6次式です.

底面流速

河床変動計算においては2次流の影響を考慮した底面流速の方向(流砂の方向)が重要になります. ここまで得られた流速分布から底面流速を求めます.

(508)\[\begin{split}\left. u_n\right|_{z=0}=A_n f_n(0) \\ =\cfrac{A_n}{C_f \chi_1} \left[ -\cfrac{\chi}{{\chi_1}^2}\left(\chi^2+\cfrac{2}{3}\chi+\cfrac{2}{15}\right) +\cfrac{\chi}{{\chi_1}^3} \left(\chi^3+\chi^2+\cfrac{2}{5}\chi+\cfrac{2}{35}\right) \right] \\ =\cfrac{A_n \chi}{C_f {\chi_1}^4} \left[ \left(\chi^3+\chi^2+\cfrac{2}{5}\chi+\cfrac{2}{35}\right) -\chi_1 \left(\chi^2+\cfrac{2}{3}\chi+\cfrac{2}{15}\right) \right] \\ =\cfrac{A_n \chi}{C_f {\chi_1}^4} \left[ \left(\chi^3+\chi^2+\cfrac{2}{5}\chi+\cfrac{2}{35}\right) -\left( \chi + {1 \over 3} \right)\left(\chi^2+\cfrac{2}{3}\chi+\cfrac{2}{15}\right) \right] \\ =\cfrac{A_n \chi}{C_f {\chi_1}^4} \left(\cfrac{2}{45}\chi+\cfrac{4}{315}\right)\end{split}\]

なお, 一般的に良く2次元モデルで使用される横断方向( \(n\) 方向)の底面流速式は

(509)\[\left. u_n \right|_{z=0} = \left. u_s \right|_{z=0} N_\ast \cfrac{h}{r_s}\]

のような形で主流方向( \(s\) 方向)の底面との関係で示される場合が多い. 主流の底面流速は,

ですので,

(510)\[\left. u_n \right|_{z=0} = \cfrac{\chi}{\chi_1} N_\ast <u_s> \cfrac{h}{r_s}\]

となります. 一方, (508) 式で平衡状態の \(A_n\) を与えれば,

(511)\[\left. u_n \right|_{z=0} =\cfrac{A_n \chi}{C_f {\chi_1}^4} \left(\cfrac{2}{45}\chi+\cfrac{4}{315}\right) =\cfrac{\chi}{C_f {\chi_1}^4} \left(\cfrac{2}{45}\chi+\cfrac{4}{315}\right) <u_s>\cfrac{h}{r_s}\]

ここで, (510) 式と (511) 式を比較すると,

(512)\[N_\ast =\cfrac{1}{C_f {\chi_1}^3} \left(\cfrac{2}{45}\chi + \cfrac{4}{315}\right)\]

と表すことが出来ます.

2次流強度について

例えば, \(\alpha=\cfrac{\kappa}{6}=0.077\) , \(C_f=0.01\) とすれば, \(N_\ast=7.03\) となりこれが, 2次元河床変動計算で良く使用される 2次流強度 \(N_\ast=7\) の根拠 ということに なります. 逆に \(N_\ast\) を条件として与えた場合,

(513)\[C_f =\cfrac{1}{N_\ast {\chi_1}^3} \left(\cfrac{2}{45}\chi + \cfrac{4}{315}\right)\]

となります.