Microsoft PowerPoint - 公開用.pptx

摩擦速度は粗度層の流れを表しているか？

－粗度層を含む河床近傍の流れの非平衡性

を考慮した解析に基づく検討－

中央大学研究開発機構 内田龍彦

摩擦速度(底面せん断応力)に関連する

キーワード

流れの抵抗

土砂の掃流力

底面付近の流れ

物質交換

粒度分布

底棲生物

導入：

摩擦速度(底面せん断応力)，

粗度とは何か？

摩擦速度（底面せん断応力）とはなにか？

吉川秀夫：水理学，技法堂出版，1976．

等流状態の断面平均底面せん断応力：







gRI

明確な定義

底面せん断応力の計測)方法

1) 重力項との釣り合いから勾配と水深(径深)を用いて算出：正しいが平均値が計算され，局所的に成立しない

gRS

U

_*



2) 対数分布則：少なくともinnner layer で平衡状態

 

:

0

.

2

ln

1

_* *







h

z

layer

inner

A

z

U

U



3) (等流のせん断応力分布から)底面付近のレイノルズ応力(分布)を用いて算出：底面 で平衡状態

'

'

1

2 *

u

w

h

z

U















 

4) 粘性底層の流速分布から算出：滑面のみ(河川では用いれない)

6) 直接計測：実験室でも粗面だと難しい？

5) 乱れエネルギー法：底面近傍で乱れが平衡状態

9

.

0

,

2

.

0

,

'

2 _{1 2} 2 1 *



C

k



C

w

C



C



U

_{b b}



参考文献：

Zezu, I. & Nakagawa, H.(1993) Turbulence in open channel flow, A.A.Balkema, Rotterdam, Netherlands.

Rowinski, P., Aberle, J. and Mazurczyk, A.(2005) Shear velocity estimation in hydraulic research, Acta Geophysica Polonica 53(4), 567-583.

Bagherimiyab, F. & Lemmin, U. (2013） Shear velocity estimates in rough-bed open channel flow, Earth Surf. Process. Landforms 38, 1714-1724.

流れの解析におけるせん断応力とは？

粗面抵抗則，粗面

日野幹雄：明解水理学，丸善株式会社，1983．

流速と摩擦速度の関係式を見出す：



,

,

_



*

h

k

f

U

U

s



粗面の底面せん断応力の評価方法

(流れの境界条件)

1) 平面二次元解析：流速鉛直分布が等流

2) 三次元解析(準三次元解析)1：流速評価点より下で平衡状態

2 2

U

U

C

_i bi







)

/

ln(

)

/

1

(

/

1

1

,

6 / 1 s

k

h

Ar

h

n

g

C























2 2 b bi b bi

u

u

c







)

/

ln(

)

/

1

(

1

s b b

k

z

Ar

c







流速評価点で乱れの非平衡性を考慮1),2)するもの：

)

ln(

1

* 2 / 1 2 / 1 2 *

k

u

c

k

z

E

c

k

u

u

s b b r b b  



_

参考文献： 1) Launder, B.E. and Spalding, D.B.(1973) The numerical computation of turbulent flow, Comput. Methods Appl. Mech. Eng., 3, 269‐289, 1973. 2) Nicholas, A.P. (2001) Computational fluid dynamics modelling of boundary roughness in gravel‐bed rivers: an investigation of the effects of  random variability in bed elevation, Earth Surf. Process. Landforms 26: 345‐362. 

仮定：

鉛直方向の流れが無く，流速評価点より下では流れが平衡状態であることを仮定し，底面流速と底面せ

ん断応力の関数を幾何条件，(水理条件)などで表す．

底面せん断応力を底面の境界条件とする場合は，底面近傍で平衡状態の流れを仮定する必要がある．

粗面でそのような平衡状態で表される領域があるのか？

粗面を河床(底面の境界条件)として取り扱

わない評価方法

参考文献： 1) Olsen, N.R.B and Stokseth, S.(1995) Three‐dimensional numerical modelling of water flow in a river with large bed roughness, Journal of  Hydraulic Research, Vol.33(4), pp.571‐581. 2) Nikora, V., McLean, S., Coleman, S., Pokrajac, D., McEwan, I., Campbell, L., Aberle, J., Clunie, D. and Koll, K. (2007) Double‐Averaging Concept  for Rough‐Bed Open‐Channel and Overland Flows: Applications, J. Hydraul. Eng 133: 884‐895.

壁法則を用いた評価方法と透過性の抵抗

体として評価する方法1)

RANS方程式を空間平均したダブルアベレージングを用いて樹

木などと同じように透過性の抵抗体として評価する方法2)

Nikoraら(2007)による開水路粗面流れの領域区分

具体的にどのように評価，解析するのか？従来の抵抗則，底面せん断応力との関係は？

粗面抵抗則と多重粗度スケール性：

粗度と形状の区別

水深全体を考える場合は粗面

河床近傍を考える場合は形状

解析方法：

粗度の多重スケール性を考慮した解析法

（考え方）

粗度のスケール分離と解析領域

主計算領域（水深

h

）

浸透層

底面 z

_b

河床高(粗度表面)z

_t

水面 z

_s

粗度要素(SSR)

河床の起伏(LSR)

Flow

抵抗則

粗度の多重スケール性を考慮した解析方法

小さな粗度(SSR) ：

計算で直接評価できない形状⇒抵抗則(壁法則)でモデル化

大きな粗度(形状)(LSR) ：

SSRを滑らかな河床面としたときに現れる形状⇒流体の

運動方程式から直接評価

大型粗度の抵抗：相対水深が小さいと相対水深に依存する．

大きな粗度を有する浅い流れ：対数分則が適用できない．

水深を用いて粗度のスケール分離

水深

Ks <<h

小さな粗度(SSR)

大きな粗度(LSR)

Ks ～h

Ks ＞h

粗度のスケール分離のイメージ

砂

礫

巨石

砂堆

砂漣

交互砂州

粗度スケール大

河床材料

河床形態・河床形状

粗度のスケール分離

大きな粗度（LSR）の評価方法

抵抗則/壁法則

主計算領域

Flow

LSR

SSR

抗力

剥離

主計算領域を細かな計算格子を用いた三次元モデルで計算することは

計算負荷が大きい（水深より少し小さなスケールまで解ければよい）

 抵抗則領域の流れの非平衡性を考慮した非静水圧準三次元解析法（底

面流速解析法）

小さな粗度（SSR）の評価方法

壁法則

主計算領域

Flow

LSR

SSR

LSRによる流れの乱れが起こす

付加的な運動量交換

従来の抵抗評価法（粗面壁法則）

Ar

k

z

u

u

s b b





)

ln(

1

*





一般化

...)

,

(

* s b b

k

z

f

u

u

_

_

壁法則を用いる薄い領域において流

速分布が相似 (平衡状態) 

LSRによる底面近傍の複雑な流れを評価するためには，

非平衡粗面抵抗則

が必要

解析方法：

粗度の多重スケール性を考慮した解析法

（定式化）

0













j j

x

h

U

t

h





i j i j j i t ij i b b i j ij bi i s j j i i

u

u

x

U

x

U

v

v

T

x

z

dp

x

hdp

x

hT

x

z

gh

x

h

U

U

t

h

U

'

'

,

0





























































































水深積分された連続式

水深積分された運動方程式

底面せん断応力

非静水圧分布項

流速鉛直分布による

運動量交換項

水平せん断応力

2: 圧力非静水圧成分の鉛直分布









h

z

z

u

u

u

U

u

u

U

u

u

u

u

u

s bi si i i si i i i i i i i

/

)

(

,

,

,

'

3

4

1

12

12

'

3 2 3 2















































1: 流速鉛直分布:

)

(

,

dp

P

g

z

z

dp

dp





_b







_s



3: 底面の壁法則 （平衡状態）























s b b b b

k

z

z

Ar

c

u

u

u

0 * 2 *

ln

1

,











底面流速解析法

z

x

h



z

_b bx

u

sx

u

y b

P

_ y



x

U

Flow

y s

P

_

z

_s

z

_b

k

x

u



W

P

_b

流れの解析法の枠組み

主計 算領 域（ 水深 h ） 浸透層 底面 zb 河床高(粗度表面)zt 水面 zs 粗度要素(SSR) 河床の起伏(LSR) Flow 抵抗則

水深積分連続式: h 水深積分乱れエネルギーの輸送方程式: k 鉛直方向流速の方程式（水深二重積分連続式）: Wh 底面流速の方程式（水深積分渦度）.： u_bi 0        b j j w x h U t h   _                k i k i j j P x k vh x h x k U t k 1 i b b i s s i j ij si bi x z w x z w x Wh h u u              ₃ 流速 と圧 力の 鉛直分 布 水深積分運動方程式.: U_i                                   i b b i j ij bi i s bi b j j i i x z dp x hdp x h x z gh u w x h U U t h U _{ }  ₀  _ 水面の運動方程式： u_si 水深積分渦度方程式：_i bi b j ij i i i _w x hD P R t h      _         si i s z z j si sj si _P x z z dp g x u u t u s                            







j j



j j m j j j u k u k h x w x z U x h U h Wh _ 2 ₂  ₁                    水深積分鉛直方向運動方程式: dp_b j zj j b bj b b j j b x h x z w w x W hU t hW dp                  z x h zb bx u sx u y  x U Flow y s P z_s pb k

w

_

_b

W

主計算領域の解析法：

大きな粗度による流速と圧力の鉛直分布を

考慮できる非静水圧準三次元解析法

渦層，粗度層で平衡状態の流速分布を仮定し，これらを通過するフラックスを無視























s b b b b

k

z

z

Ar

c

u

u

u

0 * 2 *

ln

1

,











底面せん断応力:

渦度の生産項:



























0 0 *

ln

2

,

3

,

z

z

h

z

z

h

u

C

h

v

C

P

b b be p bi bei tb p i b





















  

粗度層

u

_b

渦層



z

₀



z

_b



_b

, P

_b

_



_b

主計算領域（BVC法）

z

u

z

_b



_be

z

_t

平衡流れ（等流）

非平衡流れ

18

従来の抵抗則とその問題

仮定

:

非平衡粗面抵抗則

b ti b bi i b b j vi vj vi

z

z

x

gz

p

x

u

u

t

u





































1

(

)

渦層の運動方程式

粗度層の運動方程式

r i r ti i b b j ri rj ri

z

D

z

x

gz

p

x

u

u

t

u



































1

(

)

粗度層



_z

_r

_=a

_r

_k

_s

u

_b

渦層



z

_b



_b



_b

u

_v

u

_r



t

w

_

_t

w

_

_b

主計算領域（BVC法）

D

浸透層（無視）

渦層と粗度層の連続式

,

,

i ri t t i vi b t b

x

u

z

w

x

u

z

w

w















_



_



 19

: 粗度層の空隙率

せん断応力項と流体力項の評価方法





h

u

u

A

v

z

u

v

_{tb b} bi vi b i t bi





























渦層に作用するせん断応力

:

平衡状態で対数分布を仮定: 

v

_tb

: 水深平均値に換算された底面の渦動

粘性係数

)

(

1

v b b

c

c

A



































































h

z

z

h

C

u

u

c

Ar

k

z

z

u

u

c

Ar

k

z

z

u

u

c

r b D r r s b v v s b b b



















8

,

2

/

ln

1

,

ln

1

* 0 * 0 *





h

z

h

h

u

u

A

v

z

u

v

_{tt t} ri ri b b i t ti







_

_























粗度層上面に作用するせん断応力

:

v

_tt

: 水深平均値に換算された粗度層上面

の渦動粘性係数

)

(

1

r t t

c

c

A







流体力項

2

1

r r ri i r r i

c

u

u

A

F

z

z

D









2 2 2

,

4

,

2

A

d

d

A

u

A

C

F

D r D D











平衡状態で対数分布を仮定: 

基礎方程式から見た主計算領域と

現象スケールの定義

：対象スケール

（水平スケール）

：代表水深

0

0

L

h

s





： 代表流速

浅水流パラメータ

0

0

U

W



： 代表鉛直流速

0 0

L

h

s





対象スケール

対象スケール

大

小

対象スケール

小スケールの現象（狭域）

大スケールの現象（広域）

1



s



1

.

0

01

.

0





_s



001

.

0



s



水深平均流の運動方程式

最も大きい項を1としたときの

各項の大きさ

0.0 

0.2 

0.4 

0.6 

0.8 

1.0 

10

‐3

₁₀

‐2

₁₀

‐1

₁

0 0

L

h

s





2D

3D

準三次元解析法

1

.

0

0 *



U

u

j ij i b b i i i i b j j i i

x

h

x

z

dp

x

hdp

x

gh

x

z

gh

x

h

U

U

t

h

U





















































































_

_







₀ 2

/

2

水深積分された運動方程式

移流項

重力項，底面

せん断応力項

非静水圧項

非静水圧

底面流速の方程式

z

x

y

構造物



_x

_

y

u

_sx

u

_bx

U

_x

U

_y

u

_sy

u

_by

x

y

W

渦線

渦線

流線

流線

h

x

z

w

x

z

w

x

Wh

h

u

u

_{bx sx y s} s _b b

























底面流速の方程式

si i s j i s sj i s

P

x

z

g

x

u

u

t

u





















'

j ij i i i

x

hD

P

ER

t

h

















  

水深積分渦度方程式

水表面の運動方程式

(底面流速場)=(水深平均流速場)＋(渦度分布の影響：水面流速，渦度の移流項)＋…

＋(圧力分布の影響：鉛直方向流速の場所的変化)

底面流速の決定要因

底面流速の方程式

10

‐3

₁₀

‐2

₁₀

‐1

₁

₀ 0

L

h

s





2D

3D

0.0

0.2

0.4

0.6

0.8

1.0

準三次元解析法

鉛直流速

(底面流速場)=(水深平均流速場)＋(流速鉛直分布の変形：水面流速，渦度の移流項)＋…

＋(鉛直方向流速の場所的変化)

2D

3D

準三次元解析法

非静水圧

[水深平均流速場の場合]

底面流速を用いて従来の抵抗則から求めた

底面せん断応力は河床に作用する力か？

粗度層

底面流速

u

_b

渦層

底面せん断応力



_b



_b

u

_v

u

_r



t

w

_

_t

w

_

_b

主計算領域

流体力

D

河床極近傍

を拡大

不等流

等流



_b

u

_b

z

u

_u

z

渦層の運動方程式の各項の大きさ

0.0 0.5 1.0 1.5 2.0 2.5 3.0 重力項と底面せん断力項 0.0 0.5 1.0 1.5 2.0 2.5 3.0 重力項と底面せん断力項

重力項，せん断力項に対する

各項の大きさ

k

_s

/h=0.01

_k

s

/h=0.1

3D

2D for u

_b

2D for U

10

‐2

₁₀

‐1

₁

0 0

L

h

s





10

‐2

₁₀

‐1

₁

0 0

L

h

s





3D

2D for u

_b

2D for U

b ti b bi i b i b i j vi vj vi

z

z

x

z

g

x

dp

x

h

g

x

u

u

t

u













































1

重力項・せん断応力項

（平衡粗面抵抗則で考慮される項）

非静水圧項

移流項

水深変化項

e

_s

>0.1：渦層内の流れの非平衡性は無視できない．

相対粗度が小さいほど渦層内の非平衡性は大きくなる．

⇒渦層内の非平衡流れは砂河床の局所洗掘の解析において特に重要となる．

（洗掘部では水深が大きく，流れスケールは逆に小さいため）

0.0 0.5 1.0 1.5 2.0 2.5 3.0 重力項と底面せん断力項 0.0 0.5 1.0 1.5 2.0 2.5 3.0 重力項と底面せん断力項

3D

2D for u

_b

10

‐2

₁₀

‐1

₁

0 0

L

h

s





10

‐2

₁₀

‐1

₁

0 0

L

h

s





3D

2D for u

_b

重力項，せん断力項に対する

各項の大きさ

k

_s

/h=0.01

_k

s

/h=0.1

粗度層の運動方程式の各項の大きさ

r i r ti i b i b i j vi vj vi

z

D

z

x

z

g

x

dp

x

h

g

x

u

u

t

u











































1

重力項・せん断応力・流体力項

（平衡粗面抵抗則で考慮される項）

非静水圧項

移流項

水深変化項

2D for U

e

_s

>0.1：渦層内の流れの非平衡性は無視できない．

渦層とは逆に相対粗度が大きいほど粗度内の非平衡性は大きくなる．

⇒粗度層は通常河床と定義される面より下であるため，

礫床河川の局所流解析では非平衡粗面抵抗則が必要であることが分かる．

水深 h

₀

計測範囲の平均水深(0.193 m)

SSR：河床材料 (d

₀

=0.019 m) 

相当粗度 k

_s

=0.9d

₀

，原点位置



z

₀

=0.3d

(巨石が無い条件で算出)

LSR：巨石 (d

_c

=0.055 m)

河床高の変化

-0.06 -0.04 -0.02 0 0.02 0.04 0.06

z （ｍ）：

B=0.91m

流量 Q (0.127m

3

_/s)

実験条件

アスペクト比: h/B

0.21

水路勾配: S

0.003

Fr 数

0.52

相対粗度高さ: d

₀

/h  

0.099

相対巨石高さ: d

_c

/h  

0.28 

粗度設置間隔 : l/d

_c

6.0

h: 水深, B: 水路幅, Fr2_=U2_{/gh, U: 平均流速, d} 0:  河床材料直 径, d_c: 巨石直径, l: 巨石対角線間隔

孤立化した水没巨石が点在する流れの実験

Papanicolaou, A.N., C.M. Kramer, A.G. Tsakiris, T. Stoesser, S. Bomminnayuni, and Z. Chen. 2012. Effects of a fully submerged boulder within a  boulder array on the mean and turbulent flow fields: Implications to bedload transport. Acta Geophysica. 60(6):1502‐1546.

0.00

0.02

0.04

0.06

0.08

0.10

0.12

0.14

0.16

0.18

0.20

0

0.2

0.4

0.6

0.8

z 

(m

)

u (m/s)

LDV

Log‐law

BVC

2D

巨石なしの場合の底面粗度：

k

_s

=0.9d，z

₀

=0.3d

巨石が無い場合の相当粗度の算出





















s b

k

z

z

z

Ar

u

u

₀ *

ln

1





Fig.1 Definitions of heights of bottom and sediment

particles z z0 Bottom level zb Water level zs h zb Bottom velocity ub

Water surface velocity us

Roughness surface level zt

z0

Velocity u

粗度を通過する縦断面の流速鉛直分布の比較 (y/D=0)

0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2

x/d

_c

=‐5

x/d

_c 

=‐4

x/d

_c 

=‐3

x/d

_c 

=‐2

x/d

_c 

=‐1

x/d

_c 

=0

x/d

_c 

=1

x/d

_c 

=2

x/d

_c 

=3

x/d

_c 

=4

x/d

_c 

=5

0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2

EXP.

BVC‐NWL

BVC‐EWL

2D

32

U m/s

z (m)

粗度間の縦断面の流速鉛直分布の比較 (y/D=1)

33 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2

x/D=‐5

x/D=‐4

x/D=‐3

x/D=‐2

x/D=‐1

x/D=0

x/D=1

x/D=2

x/D=3

x/D=4

x/D=5

0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2 0.00 0.05 0.10 0.15 0.20 ‐0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 exp Δ=D/4 Δ=D/3 Δ=D/2

EXP.

BVC‐NWL

BVC‐EWL

2D

U m/s

z (m)

二次元解析法(k

_s

/d=0.9, z

₀

/d

₀

=0.3) 

BVC‐EWL (底面流速解析法＋平衡抵抗則)

(k

_s

/d=0.9, z

₀

/d

₀

=0.3) 

-0.3 -0.2 -0.1 0 0.1 0.2 0.3 -0.1 0 0.1 0.2 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 -0.1 0 0.1 0.2 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 -0.1 0 0.1 0.2 Us0: 0 0.2 0.4 0.6 0.8 1 1.2 1

実測結果 (=0.1, z

_s

z)/h)

-0.3 -0.2 -0.1 0 0.1 0.2 0.3 -0.1 0 0.1 0.2

BVC‐NWL (底面流速解析法＋非平衡抵抗則)

(C

_D

=0.4, z

_r

/k

_s

=1.0) (k

_s

/d=0.9, z

₀

/d

₀

=0.3)  

₃₄

相対巨石高さ: d

_c

/h=0.28  

実験と計算結果の比較：水面流速

-0.3 -0.2 -0.1 0 0.1 0.2 0.3 -0.1 0 0.1 0.2 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 -0.1 0 0.1 0.2 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 -0.1 0 0.1 0.2 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 -0.1 0 0.1 0.2 Uabs: -0.1 0 0.1 0.2 0.3 0.4 0.5 1 35

二次元解析法(k

_s

/d=0.9, z

₀

/d

₀

=0.3) 

BVC‐EWL (底面流速解析法＋平衡抵抗則)

(k

_s

/d=0.9, z

₀

/d

₀

=0.3) 

実測結果 (=0.1, z

_s

z)/h)

BVC‐NWL (底面流速解析法＋非平衡抵抗則)

(C

_D

=0.4, z

_r

/k

_s

=1.0) (k

_s

/d=0.9, z

₀

/d

₀

=0.3)  

実験と計算結果の比較：底面流速

BVC法と非平衡粗面抵抗則を導入すること

による底面流速場の解析精度向上の効果

0.00

0.10

0.20

0.30

0.40

_{between boulders}

across boulders

2D

BVC‐EWL

BVC‐NWL

0 exp bcal ) ( RMS U u u  _b

二次抵抗（粗度の

増加）

一次抵抗（流体力）

BVC method 効果

非平衡粗面抵抗則効果

-0.3 -0.2 -0.1 0 0.1 0.2 0.3 -0.1 0 0.1 0.2 WVst: -0.1 -0.06 -0.02 0.02 0.06 0.1 1 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 -0.1 0 0.1 0.2 WVsb: -0.2 -0.1 0 0.1 0.2 1 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 -0.1 0 0.1 0.2 Ub0: -0.1 0 0.1 0.2 0.3 0.4 0.5 1 Roughness layer zr=arks

u

_b

Vortex layer z_b



_b



_b

u

_v

u

_r



t

w

_

_t

w

_

_b

BVC calculation layer

D

底面付近の流れが解析できるために，より実際的な流砂運動の解析法を検討できる．

Ur=1 m/s Uv=1 m/s Wb(m/s) Wt(m/s) Ub=1 m/s Ub (m/s)

渦層と粗度層における流れの解析結果

まとめ

流速鉛直分布を決定する渦度の変形と圧力効果を考慮できる

非静水圧準三次

元モデル(BVC法)

に底面極近傍と河床表層の流れの非平衡性が解析できる

非

平衡粗面抵抗則(DWL)

を組み込んだ解析法を開発し，河床に作用する底面せ

ん断応力の評価方法について検討した．

運動方程式を無次元化して各項のオーダーを比較することにより，流れの三

次元性，非静水圧分布及び非平衡粗面抵抗則がどの程度の流れのスケールか

ら重要となるかを明らかにした．

開発した数値解析法は水没する孤立粗度近傍などの局所流場における底面近

傍絵流速を再現できることを示した．解析結果から，流れの三次元性が強い

局所流においては底面近傍においても流れの非平衡性が強く，抵抗則を適用

する底面流速場と粗度層内の流速場は大きく異なることを示した．