植生帯を有する複雑な地形上の風速場及び海塩粒子拡散の数値実験の試み: University of the Ryukyus Repository

Title

植生帯を有する複雑な地形上の風速場及び海塩粒子拡散

_{の数値実験の試み}

Author(s)

仲座, 栄三; 津嘉山, 正光; 山路, 功祐

Citation

琉球大学工学部紀要(48): 57-66

Issue Date

1994-09-01

URL

http://hdl.handle.net/20.500.12000/15869

Rights

57

The Simulating of Sea Winds of Coastal Areas with Complex Topography and

Windbreak Forests, and the Analysis of Salt Particles Contained in Such Winds

Eizoh N

AKAZA·

Seikoh

TSUKA YAMA*

Kohsuke Y

AMAJI"

Abstract

The Okinawa Islands are frequently hit by typhoons. There are many more sea salt particles in summer typhoon winds than in the winds that blow during other seasons. Disasters caused by these salty winds occur much more often in the Okinawa prefec-tureth~Ulany other prefectu re of Japan.

Civil engineers plan and design artificial windbreaks made as countermeasures against the disasters caused by the salty sea winds. Such windbreaks arc usually areas of plant life, such as trees, planted by the coast. The location and size of these areas of plant life depend completely on the judgement of the engineers. If the judgement of these engineers is inaccurate. the windbreaks would simply be a waste of time, effort. and money; therefore. an understanding of the salty sea winds and an analysis of the sea salt particles they contain must be done.

This paper describes and explains a method of simulating the salty sea winds. Us-ing the numerical simulation, the sea winds can be studied and the sea salt particles can be analyzed. Finally, this method can help engineers understand the way the sal· ty sea winds flow over coastal areas that have windbreak forests and complex topography.

Key words

Navier-Stoks equation. Vv'indbreak forest. Direct numerical simulation, Grid generation

ifij:i~!:i.p

G

lI~*-ef ~ ihjJ.~f~:-hlftlrrq)tJi:ltXl: rul-~- ~ fi)f'~'E

li,

!tl~~(I~t!J·< ,j'G~~il-(I.-\~ (f~Uitli, X~;J;J'

III

r.J:lG,

1966.1967a.b,~'?i·~'f.1970). L,j'L~i.li"G. JIl:

~q)1e'I:.q) j.iJ.:.:.

7;

j..,

"1

~

b

Iffsh-;ib.h-, ~h ~M

c) i.I'

q)-Jf~-~.:c-f

JHt

L

-C

tl~Wi ~ ;flt.:if!ttU.1':-T·iJ:V!tQ)fJtt~·

r·

ifilJJ-\;;

li

ffi

C'

7&:'

<,

;£

t.:

~I~~(7).ltl!Jf~ ~

it

~

c

L

t.:

iiHl.:c--T

}vli.:

*1."1

~tJ£~ ~

tt

-C

I.-\~I.-

\, .:

(7)J: oj~.:.

Ci.I'

G,

JJL;(£q)

C

'::'1::>IW}£t:t+q)~~U~~tli. illtnl-#q)iitU~

~~ ~Hfrt

Lt.:

I), MU:.I:i;)lt~:g·jffiQ)f1tt~c

L-Cq)

f&1.-

\~~rf;f)h~!J)jfti.li"~I.-\.

')j, 1&j5:Q)~1it!iifL{~)Jq:q)~1Mlicl6~""1

L

<,

IhJ

IIU

c

!'L~ ~

it

i.I'l ttlli, .:..:

-~IUHm

c

L

-C

1,,\~jIVMq)tu:

ItXIllJ

~Ii~I~~J-."9fiC'rJ1~

un

~

<

M.f1j-

L

1~~ ~

l'stllW

I.:.it

L-(

It\7.J

t

r:J

it

~(Nishikawa,1985,Shima, 1988).

L

lP

L

~

-/}i"

c

y , ,]fHi1~ ~ Wf.11~ ~;it&')

h.li*

&')~

11

c\

~!dit!aH':

Jtl:i.I{1eIlXEt·H.:m* L,

j..:c-I)-:#Ht:~i;I·n:p;j"Uq)iIiii.l'

c~ ~).£.m~:~ffl Lt!f.~

<

~ ~

CI.- \

-jPrHm:ii.l{~

t.::

o.

"1

t:..

JJltEcf;~?iLmt.:c--T}v'i, ,t~!a6-C~f1!~ttA~<Ji> t),

liJf?l:

'~J1(:1994iF 5

fJ

178

•

I~$JJ;ttjU1!IDtIq:*;j- Dept. of Civil Engg. and architecture. Fac. of Engg.

5８仲座・津嘉山・[11路：植生帯を有する複雑な地形上の風速場及び海塩粒子拡散の数値実験の試み

帥一町

十

助一欧

ｌｌ Ｌ 者の数ほど乱流モデルが存在するといっても過言でな い．苔らには，それらのモデルからシミュレートきれ る結果の適合性は，計算ケースに依存して変化し，例 えばuが合えばvで合わなくなるといったものが多い いずれにしてもこれまで提案されている乱流モデルの 多くは，膨大な計算時間を要する上に，完壁なものと は言えない．本研究はⅡ海陸風や飛塩拡散に対する完 結モデル(closermodel)を提案するものではない．む しろ，現在あるモデルの中から最も簡単な計算モデルを 選択し，それで現実に起きている現象を解析する努力 を行う．それを基に，現地や実験では測定し難い現象の 理解と予測を行う．さらに，複雑な地形条件下での風 の場や飛塩拡散現象を明らかにすると共に，植生の配 置に関する基礎的なデータを得ることを目的とする． である． ここで 本数値計算ではこれらの基礎方程式を一般座標系に 対応できるように座標変換を行ったものを用いてい る．基礎方程式の差分化は空間微分項には中心差分を 用い，時間微分項には前進差分を用いた．非線型,慣性 項には計算安定を図る目的でKawamura&Kuwahara (1985)らの高精度風上差分を適用した． 2.2海塩粒子渡度の拡散 海塩粒子濃度の拡散を規定する式としては，以下に 示すFickの拡散方程式を用いる． Ｍ （０託 ＋一．

（『西「「》調》，（一釧弓一引制》

ＭＭ

ｆｌＷ〃Ⅵい－、ＪＩＩｌｌｌＩ

８｜帥へ。｜へ⑪

＋

｜ｘ｜町

０ Ｗ 一一

一忙一砂

。ｖ● ＋

一死｜弘

一必 十

一誕一趾

2．基礎理論及び計算手法 2.1風速場 数値計算で基礎となる理論は，以下に示す連続の式 とNavier-StokeSの方程式である． ｄｉｖＵ＝0 ここで，Ｃは時間平均濃度，Ｋはそれぞれ乱流拡散 係数，、ｍは分子拡散，ｗＯは海塩粒子の降下速度， Ｑｕｌは植生の飛来塩分に対するfilter効果を意味する． 飛塩拡散に関しては、砕波iIlf内に海塩粒子の大きな 供給源が集中するものとした（村上ら，1992,仲座， 1993） 飛来塩分に対する植生の影響は－種のフィルターと して作用ものと考えられるここで問題となるのが植 生によるフィルター効果をどうモデル化するかである が，本数値計算では次式のようにモデル化を行った．

qww薑‐;v砕呈lKl卦昊ﾄ卦鼎

ここでＵ二Ｍとする．

Ｐｂの項は植生の効果を考慮した項である．植生の 影響は流れの場に対しては流体抵抗として作用すると 考え次式で与えた．

Ｅ=-Q1sl叩

Ｃｌ卿＝ｓＹｃｌｄＵ(」

ここで，ｃ,"は抵抗係数であり，ｓは葉而穂密度であ る．葉、積齋度とは葉の繁})の程度を表すもので，よ く繁っている程この値は大きいということになる通 常の植生帝は5.0程度であると言われてお}〕，ｌＩｌ１繩県 のようなIIIi熱梢地域における樅生補は１０．０程度の値を 持つものと考えられる．この葉面横密度に関する研究は Kondo・Akashi（1967)，神H1ら（1992）のものがある． 通常，実際の数値計算では，連続の方程式が直接)Ⅱ いられるのではなく，次に示すPoisson方程式が解か れる．この方程式は運動方程式の発散をとることに よって導かれる． ここで，ＣＭは植生による単位時|Ｈ１単位体祇当たり の飛塩除去1ii，γごは定数，Ｃ化)及び〔/(x)はそれぞ゛れ 海塩粒子濃度及び風速の関数であるここで先ず， Ｃﾋﾉ及び【ﾉ(ｽﾉの関数形を決定しなければならないが， これらの関数形を決定し得る程の詳細な現地観測等の データはこれまで得られていない．そこで，本研究で は，肢も簡単な関数形を用いることとし，一次の関数 形とした． 2.3基礎方程式の￣般座標系による表示 本研究では，複雑な地形上の風速場や飛塩拡散を計 算対象としているため，地形の境界条件が比較的簡単 に与えられるように，座標系に一般座標系を用いた解 析手法を用いた．上述の基礎方程式は，座標系のみに

V2P＝divq＋divq-2qSUgradU

-ｄＭＵｗ芸-筈

琉球大学J:学部紀要第48号，1994年 5９ 図－１に示すような変換を施すと以下のとお')にな る． （ｘ方1,1のみ）

叫十'`に';+叩壯に』`:十,’１`J＝

-;(剛Ⅲ派にい,K胸,ｗ÷

に火;十Ｍm)に｣`:+,Luq1十ｋ▽}`+ﾉｱ；

ここで:今いぅ凡=ｻﾞ､諄

▽経姜鶚弄

であり，Ｊはヤコビアンで,ノーｘ５ｙｉｍ－為ｙ§である。

Choi・ＭｏｉｎａｎｄＫｉｍ,1993)．本研究では，植生による 乱れのプロダクションのモデル化や気温変化をも孝噸 する予定であるので，現在のところLESが本研究でめ ざすl汁飾手法としては有力である．しかしながら，計 算によって得られる結果を現地観測値や室内実験値と 比較すると，かならずしもいずれが勝るとも言えない 鰯合が多い．要は数値計算の対象とするものが何なの か？あるいはどこまで詳細に見たいかによる．ここに 示す研究に限れば，数値解析手法の梢度向上が目的で なく，複雑な地形上の大気の流れや飛来塩分の拡散の 理解が当miBlhrすところであり，また植生の配Iiilの違 いによる飛来塩分砒の空間的な分布の仕方に興味があ る．また。実務被ii↑での利用を考えているのでワーク ステーションによる解析で２日以内に計算が終了し得 る程度の計算モデルを採用することとした．そうする と，やはり現実的な選択として，乱流場を渦動粘性係 数を川いることによって，運動flt交換係数を空llI1的に 一定とし，実limりなReynolds数を下げ，ｊＩｉｉ流的な解析 を１jう手法が選'１(きれる．ここで示す研究に限っては このモデルをｌｌ１いた解析を示すが，現在LESを川い， またワークステーション等でも実務に耐え得るような モデルの開発'１１であり，その結来については別の機会 に紹介する予定である．

§－〃qyFだ、

工

§

３．計算結果及び考察 ｜ｘｌ－２(a)，(b)，(c)に計算対象とした海岸地形の概念 図を示す．図－２(a)は，汀線背後に幅約30.0ｍの一様 斜iliiの前浜部分と砂丘からなる後浜があり，ざらにそ のｆｆ後に械生l11r7Wr在する場合を示している．lil・算で は，ｈｈｌＭＭｒをlx1-3に示すように，砂丘の典｣2やiiil後 に移鋤させ,それぞれについて数値iil算を行った.図一 ２（a)は，土手を越える海風の変形を対象としたiil算 に川いた．なお，これらの計算は１１K１－２に示すスケー ルをju1定して行われているが，これは代表災さを変え て兄ることで，磁さ数百メートルの丘やベル型の''1を 越える海風及び飛塩濃度の変形|M1題として読み換える ことも可能である． 図－１CoordinaLeSystcm ここに示すとおり，地表面など境界の正確な表現のた めに･股ＩＪＩ：標系を擁川した反動として，鵬礎ｿjｲi1式に は爆発的に項数が噸えた形とな})．デカルト114標系に よる鵬本形の場合と比較して膨大なii付):IlIFIlllが必婆と される． 2.4揮励最拡散係数の取り扱い 本研究では，乱流場の解析が対象となるため，数値 計算上lBl題となるのはNavier-Stokcs中の慣性〕Ｉ､の取 り扱い方と，乱れが作り１１}す応力をどのようにモデル化 するかである．乱流場の代表的な数恢iif算手法（IIiF11I］ 平均をとらない）としては，ＤＮＳとLESとが挙げられ る（SumagriluSky,1963．Deadorffl970，MoiIu＆Ｋｉｍ・ '982,KHlwamurn＆Kuwahara，1984)．，1冊Reyn(’１(１s数 流れのDNSは，一般にｌＩＩｉ性項にlili次の風｣冒溌分を)|}い て実現ざれるがけ数'１{（的な粘性Ⅲ(の評価が難しい（し かし，IllGRcyn(’1〔l呂数ijlLｵしに関してはこのlUlI)でない，

尋

6０仲座・津噺山・山路：植生梢をｲｿﾞする複雑な地形上の風速場及び海塩粒子拡散の数値実験の試み

ｌＸｌ－３ＩｉＭＬ１ｉｌｒの配悩例 (b)ベルlllIl1

封

図－４に，iil，):格子の形成例を示す．計算格子は， 地iii条件など境界条件を与えやすくするために直行さ せてある(Thompson,1974)－格子数はいずれの場合も 85×45である．また，iil算領域のiniざを代表長苔とし たときの乱流ReynoIds数は`lx10000である． (c)平担地 |叉１－２対象とした海岸地形Model 、■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■l■■■■■■ ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■Ⅱ■■■ロ■ロ■ロ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■Ⅱ■■■ロ■ロ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■Ⅱ■■■■■■■■■■ロ■■■■■■Ⅱ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■ロ■■ロ■■Ⅱ■ロ■■ロ■■■ロ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■ロロロ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■ロロロ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■ロロ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■I■■■■■■ ロロ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ｌ■■■■■■ ロ－口￣￣■￣￣￣￣￣＝--■■■■＝￣■■■■＝--■■■■■■■■－－■■■■￣■■■■－－－－－－－－ 口 ■■￣￣ロー￣￣￣￣■■＝■■￣■■■■■■■■■■■■■■￣■■■■■■■■--■■■■■■■■－－－－－－－－■■ ￣ロー■■ロー￣￣￣－－＝＝--＝￣＝＝■■－－■■－－－－－－－－－－－－－－－－■■＝ ￣￣＝＝＝＝＝－－－＝￣＝－－＝＝＝－－－－－－－－－－－－－－－－－－－－－－－ ---■－－■---■－■－■---－－－－－－－－－－－－￣---￣--- --■－－－－■￣＝＝■￣＝■￣■－－■￣■￣■￣■￣■￣■￣￣■￣￣■---－－－－－－－￣-－---￣￣￣￣￣-－－￣ ＝--■￣■--■--■￣■￣■---＝■￣＝■￣■￣■￣■￣■￣－－￣■－－－－－－－￣￣￣￣￣---￣￣----－－－－ －－－－－■-－－－■￣■￣■￣■￣■￣－－￣■￣－－￣■￣-－￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣－－ロー￣￣ ■－■－■－■￣－－－－-－－－－－－－--￣￣----－－－－￣ 凶－４（a）格子形成例（，vl11地） BIUUU■ﾛ■■■■■■■■■■■■■■■■■■■■■■■ﾛ■ﾛ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■I■■■■＝■■I■■I■■I■■U■■U■■■■■ﾛ■■ ﾛ,ﾛﾛﾛ百口g｢~--~－５己ｮ自白：皀已ﾆﾆ目ﾆﾆﾆﾆﾆｰﾆ二号ﾆﾆﾆﾆﾆﾆﾆﾆﾆﾆ手======二二＝=－－－－~･~･~･~･~･~･･･ ■、、

！。;:ﾛ･占呂ﾖﾆ目壺且亘旦E二王弓三壺舌＝=ﾆｰｰｰｰｰｪ---.---・・!.‐：’二二二二一二コ＝P￣~~~－－￣－－－－－

：：‐さ=ﾆｰﾆ２ｺ圭ﾆﾆ旨二二F-ｺZ古=_=ﾆﾆﾆｰﾆﾆﾆ~~---~=~￣~~~~~~~~~ョ

，百Ｊ５－．ﾛｰｰﾆﾛ■ﾆｰｰｮﾆ’’百ｉ==三三三二＝－－－－－－－－－－－－－=－－－－－－－－－－－－－－－ ・・￣ＩL二三=二言～＝_二，=÷－－－．－－=－－－－－－－－－－－－－－－－-－－－．－－－－－－－－

_{..-:皀宰ﾆｰﾆﾆﾆﾆｰ二言デー~ﾏｰﾆｰｰｰｰ=====================￣=.~~~~~~~￣}

－生ﾆｰﾆﾆｰ乞一_=ﾛｰﾆｮｰｰｰｰﾏ=---===============ニーニー======･---‐平Ｈ￣￣－～Hﾆﾆﾆ号=ﾆﾆF司嵜￣'二皀ﾆｰﾆｰｰｰﾆニニーーーーーーーーーーーーー＝＝＝＝＝＝＝＝￣￣￣￣￣

;ｉｌ－ﾆｦﾆﾋｰEﾆｰﾆｰﾆ三_全==ﾆｰﾆﾆ=二丁ﾆｰｰｰｰｰﾆ__---_---1i＿ミニニヒヨニニニニーーニ白＝＝=---可一一---－

１Ｌ二二三Ｅｺﾆｰﾆﾆﾆﾆ=三三二二三三三ﾆﾆｰﾆｰｰﾆ=＝==＝三三=---=---ニニニ

ｉｉ三三二三二-毛三二三三二三三ﾆﾆﾆﾆﾆﾆｰﾆﾆｰｰﾆｰｰﾆﾆｰﾆﾆﾆﾆﾆﾆｰｰｰｰﾆﾆﾆﾆｰ~~~ﾆｰｰ

：:_二三三=二二二.ﾆﾆｰﾆﾆｰｰ=ﾆﾆｰｰﾆｰﾆｰｰﾆｰﾆｰｰｰﾆｰｰ三一三二三＝－－=＝‐－=＝

＿:.一三二三二=菫二三三ﾆﾆｰｰｰｰｰｰｰｰﾆｰﾆｰｰﾆｰﾆｰｰｰﾆｰｰｰｰﾆﾆﾆ＝－－ニーーーーー

ー￣ ロ■ ■＝－￣＝－ ＿ニニーニーーー_{ニーーニ---}＝ニーニー＝＝=－－－元－－＝＝＝ニニーーーーー･-＝-.--三一一一一一一一一一一一＝＝ ■■－ － 図－４（b）格子形成例（ステップ!(ｿ地形）

琉球大学工学部紀要第48号，1994年 6１ ⅡⅡⅡⅡⅡⅡⅡ１１ⅡⅡⅡⅡⅡⅡⅡnmⅡⅡⅡⅢllIllⅡⅡIⅡﾛﾛﾛ■■■■■■■■■■■■■■■■ ⅡUHUIIⅡIⅡ11ⅡⅡⅡⅡⅢⅡⅡⅡⅡnmⅡIⅡⅢⅢ11ⅡllImHU■■■■■■■■■■■■■■■■■■■■■■■■■■■ ﾛﾛﾛﾛﾛⅡⅡⅡIⅡⅡ11ⅡⅢⅢⅡmUUIuⅡⅡⅡⅢlⅡⅡⅡⅡIⅡﾛﾛﾛ■■■■■■■■■■■■■■■■■■■■■■■■■■ BHUⅡⅡⅡⅡIⅡⅡⅡIⅢⅡⅡⅡⅡUmUⅡⅡⅡIⅡllIllⅢⅡⅡImUU■■■■■■■■■■■■■■■■■■■■■■■■■■■ UmUHlIlIⅡⅡⅡlⅡⅢⅡⅡⅡⅡⅡⅡⅢⅡIⅡⅢⅢⅡⅡllIﾛﾛﾛﾛﾛﾛ■■■■■■■■■■■■■■■■■■■■■■■■■ ppmⅡⅡlⅡlⅡIⅡⅡⅡⅡⅡⅡⅡⅡⅡmUIⅡⅡIⅢⅡlⅡlⅡⅡlﾛﾛﾛﾛﾛﾛﾛ■■■■■■■■■■■■■■■■■■■■■■■■ ⅢﾛⅢ1111ⅡIⅡⅡIUmmⅡUmUⅡUlⅡⅡIⅡⅢlIllⅡⅡInUU■■■■■■■■■■■■■■■■■■■■■■■■■■■■ mpmⅡⅡIUIIIUmⅡUNUⅡUUUUⅡﾛⅡⅡⅡ11ⅡIⅡIⅡ11ⅡⅢﾛ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■ 図－４（C）格子形成例（ベル型111） 図－５風速ベクトル（平担地） では流れの縮流効果で境界層がかなりの勢いで先走っ ている．しかしそれは，下流に行くにつれ徐々に拡散 し，下流側では消滅している． 図－５に，水平な地面条件での計算結果を示す．風 の上流側では一様流を与えているため，境界の次の計 算点から内部境界層が発達し始める.図－６に,ステッ プ地形を越える海風の変形を示す．ステップ先端付近

6２（''１座・ilt臘山・'１１路：植生桁をｲLi｢する複雑な地形上の風速場及び海塩粒子拡散の数値実験の試み

図－７風速ベクトル（ベル剛llLt＝5.0s）

薑篝篝薑薑iii篝

譲讓=－－＝

図－８風速ベクトル（ベル型111,t＝15.qＪ 現象は，乱れのⅡ等'111スケールを過かに越えた平均流の 振舞いである．また，これ程の非定常組織渦が造り出 す乱流応力については。やはり何らかの形でモデル化 し計算へ反映させる必要がある．しかし，図－６に示す ステップ地形を越える場合にはこのような組織渦の剥 離は２ｋじず.この場合流れは定常的なものであった． ただしこの結来はReynoldS数が10000瀧皮の場合であ り，それを越えるReyI1ol〔1s数に対してはその限りでない． 図－７，８に領域内にベル型lIlを有する場合の海風 の変化を示す。それぞれ実時間で5.0.15.0秒後の様 ]亀である．図示のとおり，111頂付近から組織化した渦 の剥離が生じている．さらに時ＩｉＩｊが経過すると，111頂 からは新たな渦がストローハル数でﾉﾌﾞ･えられる周1UIで掃 きＨ１される．この所たに掃きlllされた渦は，１１１の下手側 で移流中の渦とマージングし，－(１Ｍの渦としての渦庇を ．時闘めた後，拡散しながら下流に流きれる．これらの

琉球大学工学部紀要第48号，1994年 6３ 図－９に，汀線付近で飛塩肚の鉛直分布を与えた場 合の飛塩武の分布を示す．ただし，瓶塩撒は汀線での 鹸大飛塩iiiを用いた無次元値で示してある（実際の計 算での境界条件や飛塩趾の与え力については,仲座ら， 1991,｛【|'座・津嚇１１ルlll路，1993,を参照きれたい)． 図－１０に，ステップ型の地形における飛塩j1itの拡散 の計算結果を示す.後にiiiLい繊論を行うことにして， 図-11にステップ上にＩｉｌｌ小１１Fをi没けた｣illi合のii↑算結果 を示す．また，｜XI-l2にN賑Ilfの配悩パターンと計算 結果の比較位置を示す．ざらに，図一]３(a)－(e)に。 ステップ１２の幾つかの比較位枇における飛塩並の鉛直 分布を示す．図11,には，IiMijlMrがないj〕附合（S=O)， 植生密度がS=5及びS=10の場合について示してある LineAにおいては，いずれの場合にも大差は認められ ない．しかし，ＬｉｎｅＢ以降については，下流に行くに つれて，比較的iIfIi庇が低いmiで述いが見られる． 特に，植生;''1Fから約109ｍ離れたところ（LilleE)では， 植生密度がS二10の場合，植生511fを設けない場合に比 較して約1/6となっている．図－１４に，植生梢を設け た場合の風速分布を示す．ＬｉｎｅＡ及びＢについては他 のLineと比較して，地iiiの近くで風速がかなり大きく なっている．飛塩濃度と風速との機で与えられる飛塩 の輸送鼠（fIux)で考えると，植生背後の地表面付近 の風速の落ち込みは，ＩｉＭ里による飛塩除去効果力稿い ことを意味する． 図－１５(a)－(c)に，１１Ｍﾆｲｉｆの彼|流位償の通いによる 飛塩並の鉛DII分布の変化を刀くす．械生１|ｻﾞの設置位置に よる飛塩101分布の述いは．ＬｉｎｅＥでのみ大きな違いを 示している．また，その述いも他ﾉﾋ'H1rの高きとほぼ等 しい範'111となっている．これらの結果から，植生帯の 効果が顕蒋に生じるのは，ＩｉＭ２の満さの約７倍から下 流１１１１にあると判断きれる．

H(､）

4０ 2０ 0

L(､）

ｌｘｌ－９塩分濃度分布図（平柧地）

H(､）

4０ 2０ 0.01-－０．０３ ） 00 １ 130 5０

L(､）

０ ｌｘｌ－１０塩分濃度分布図（ステップ型地形）

6４仲座．津嘉'１１．１１１路：植生ilIfを有する複雑な地形上の風速場及び海INK粒子拡散の数仙実験の試み

H(､）

4０ 2０ ₀₀₁ ０．０３

E､￣～。

５０ 100 130

L(､）

図－１１塩分機度分布ＩｘＩ（IiMijllfをｲ｢する場合）

LjncE

LｉｎｃＡＬｍｃＢＬｉｎｃＣＬｉｎｃＤ

植生帯

Ｐ、 ＣａｓｃＣ ＣａｓｃＡ 109ｍ 73ｍ 28ｍ３６ｍ５２ｍ 図一ｌ２ｉｉＩｉ１:Iihの比較位ｉｌｆ 勿苑”猶勿 宕ミミ讐也電 ‘０ ３５ ０５０ ３２２ ｛Ｅ｝』急一週国 IS ID １５ １０ ０．００ＯＤ０ＳＯ０ＪＯ ノVhcl （a)ＬｉｎｅＡ ０．１s０６２０ (Ｃ／ＣＯ） ０．，００．ｍ0.0‘、､０６０．０８０.ＪＯＯｏＪ２ ハMzcl（Ｃ／ＣＯ） （[)ＬｉｎｅＣ 勿錨”巧” 盲こ』督懲酉 ‘０ ３５

ミョ，

::：

｣Ｓ ＩＯ １５ １０ 。』 0．〃０．”０．１２０６２６０２０初 Ｎｂｃｌ（Ｏ/ＣＯ）jVbDノ （b)LineB（｡)Line、 図－１３塩分濃度の鉛直分布（次ページへつづく） ０．００ “８肱１０ (C／ＣＯ） ｈ ■ 0 ..･･･････････;････｡｡ ＄ ｉ ｉ ● ･－゜･・・・･････Ｃｌ一･････

F二ｍ！

_{。￣･￣．ご＝Ｓ･Ｕ} ･－－ｓ＝１０．０ ｡●● ●●● ：！； ０ OB・・ ：：：ｏ 、｡：

、_

月霞雰雪i電で旨篭鶯

Ｉ､、

=更三更

、 － ，￣G～

｜川口凶

１コ１ｍ

●ＧＧ Ｃ ■ ■ ■ 。 ● ■ ●醇 ●。 ●函 ●■ ■● の● －ｓ＝０．０ －．－．－s＝５．０ －－－s＝１０．０ エ勾一梺華。_{－－－－口} s← ｂ Ｄ ｂ ８ ： ： ．....－...｡･･･､-...8･･･b ●

｛

!.… ご●■●■ ､●■中 ■ －－－s＝０．０ ．．－．－．ｓ＝５．０ ・・－－s＝１０．０ ～ ５ほロ｡引 ： ｡ ｡~●●の｡●●●● ● ●■-■●▲●●｡▲の●-■の-■■●0 ………i………÷..…・・■ ■ ＢＯ Ｃ ● 》》 ■、０■●■●●■Ｇ０Ｂ０ｄ■□●●●■０□■ＵＧｑ０ｄｑ０ｑＯ‐● ● ● ● ● 》》 ● ■ ● ● ● ● －－－－ｓ＝０．０ ･－－．ｓ＝5.0 ･－－ｓ＝１０．０

、一

－－ －コ_鐸莞

言壱

=宝定 ～ Ｕ８0口8０

琉球大学工学部紀要第48号，1994年 6５ ０５０５０５０ イ３３２２１１ ０５０５０ “３３２２ 宕昌芭急一遍国

百ミミ響四国

１５ ．１０

抑jr１０．３

０１０２０ ＭｚｃＪ Ｎｂｃｌ（C／CO） (eLine-E (b)ＬｉｎｅＤ－ 図－１３塩分濃度の鉛直分布 イ０

初苑勿麺、遁皿

3s

（Ｅ『ミ碑凹国

３０ 亙彗… ２Ｓ ２０ １s １０ ３ ０ “ ８１２ェ１０． ０．００．２００イ０．６ 0.8Ｊ｡０１．２１.‘１．６

U７ＵＯ

_Ｍｚｃｌ （c)Line-E 図－１５植生帯の配置の違いによる塩分濃度の変化 図－１４風速分布（植生帯を有する場合） 4．おわりに 従来，埋め立て地や塩害の激しい地域での殖栽は， 技術者の経験や景観の観点から決定きれている．任意 地形上の海風や海塩粒子濃度の拡散の予測が定量的に 行えれば，塩害などへの有効な対策が行えるようにな る．本研究では，そうした実務設計への手がかりとし て，実務に耐えるような風速場・海塩粒子拡散に対す る計算コードの開発を行った．また，それを幾つかの 海岸モデルに適用して，複雑な地形上における海風の 変化及び飛塩の拡散特性を具体的に検討すると共に， 海岸部に植生帯を設けた場合の風速場及び拡散の変化 を検討した．計算はワークステーション(EWS4800-350)を用いて行った．計算に要したCPU-timeは，１ ０５０５０５０ イ３３２２Ｊ１

百ミ這碑色国

３０６Ｄｒｌｍ lVUzcJ (a)Line-C ０５１０１５２０３０rｌが ０２０イ０６０JｒＩＯ ３ B●．．●・・・・｡●・Ｓｃロ゛・・０ ･･･－･ロ････････－．．９．．－．－．．･･･・・・・･・・・１．－．．．.・・・｡｡・・・･・・・･Ｉ －－－ＣａｓｅＡ －ＣａｓｅＢ ……ＣａｓｅＣ l＿ .－‐ＣａｓｅＡ －ＣａｓｅＢ -.---‐ＣａｓｅＣ

匙

⑤》》 ● ● ● ■ ■ ■ Ｐ ● 一》句■

彙鞠

』■ ■ ● ● ■ ● ■ ｇ ● ■ ● ぜ ● ● ■ Ｐ

》一

Ｐ》 ●函 ■酩 一』 ご如 已

』》

一戸 一一 一四 〕～￣ Ｉ ■ ■ ■ －Ｌ１ｎｅＡ ………･ＬｌｎｅＢ －－－ＬＩｎｅＤ -.-.-Ｌ１ｎｅＥ

…j=』鐘;:1二；

=iPlijf

●･■￣･■■●｡Ｕ ● ●酌釦 ●■の●■● ● 。●●のＣＧ◆●●●ＳＦ●●●●●●■？●●●●■■●●●●●ｏ９ＣＣ●甲の●●●● ／’一

＞

垣

～ ザ ｡■ .~｡｡￣･~｡..｡￣ｳﾛｰｰ｡｡~一･･.・・Ｔ~･-...-..F･･･.｡￣｡｡｡･･･● １１ －ｓ=０．０ ･－．－．s=５．０ －－ｓ＝１０．０ ■ ● ■ ■ ロ ■ ● ■ ■ ●

穐

昌昌ミニ望 －－－ＣａｓｅＡ －ＣａｓｅＢ ……ＣａｓｅＣ 塞置白=_＝書豊 ＝＝＝

弓司雲，

？◆ ０

6６仲座・津溌１１｣・山路：植生補を有する複雑な地形上の風速場及び海塩粒子拡散の数値実験の試み ケースにつき約45時間であった．乱流モデルをなんら 用いない計算モデルに対してもこの程度の計算時間を 要する．これは，一つには一般座標系を用いた計算手 法に原因があるが，これも複雑な地形条件を想定して の現象理解の犠牲の一つといえる.先に述べたとおり， 現在LESによるよ})現実的なモデルのUiI発中であり， それをもとにした現象解明については別の機会に紹介 したい． 幸喜蕃橘（1970）：海岸付近の空中塩素錘に関する 研究，琉球大学農学部学術報告第17号． 神111学（1993）：土壌および水面の相互作用を考慮 した水文気象に関する研究，東京工業大学学位論文． Nishikawa@ｓ（1985）：Threedimentionalboundarv-Iayer-separatio【loI1apar＆boid・Tralbs.”n.ＳＯＣ､ Aeron・ＳｐＳｃｉ．，２８，ｐｐ､27-43． Shima，Ｅ・（1988)：NumericaIanalysisofmuItipleele、 menthighliftdevicesbvNavier-Stokesequationusing implicitTVDfiｎｉｔｅｖｏ１ｕｍｅｍｅｔｈｏｄ・AIAApaper 88-2574-cp Th(〕、I)sonＪ.Ｆ，（1974):Automaticnunlericalgridgen-erationoIbody-fittedcurviIinearcoordinatesvstem 【ｏｒ【ieldcontainillgandynumberofarbitrarvtwo- dimensionaIbodies．］、Ｃｏｍｐ、Phys．，24,pp274-302 Kawamura，Ｔ、andKuwahalra，Ｋ・（1985）：Direct simulaIionofaturbulentinnerⅡowbvfinite-diflcrencemethod0AIAA･ｐｐ,１－１０． Smagrinsky．』．（1963）：Mon・WeatherRev、９１， ｐｐ９９－１６４、 ＫａｗａｍｕｒａＴ、ａｎｄＫｕｗａｈａｒａＫ．（1984)：Computation ofhigbRCynoldsnumberfIowaroundacircularcvlin、 derwithsurlaceroughness，ＡＩＡＡ､84-0340． DeadorfI」．Ｗ・（1970）：Anumericalstudvofthree-dimensi()llalUlrbulentcbanmeIfIowatIargeRGynoIdS nunlbers.』・FIuidMech．‘Ｖ０１．４Ｌｐｌ).453-480． MoiI】Ｐ、ａｍｄＫｉｍｊ．（1982)：Numericalinvestigationof turbuIentchannelfIowJ・FIuidMech．,ＶＯＬ１１８.ｐｐ、 ３４１－３７７． ChoilL，ＭｏｉｎＰ、ａｎｄＫｉｍｊ．（1993)：Directnumeric alsimulationoIturbulentflowoverribletsJ、FIuid Mech．Ｖ０１．２５５，ｐｐ､503-539． ＫｏｎｄｌＩｏＪａｎｄＡ.，Akashi(1967)：NumericaIstudyon thetwodinlensionalflowinhorizontallyhomogeneous canopylayers，Boundary-LayerMeteor．、ＶＯＬ１０， pp255-272 謝辞 本研究で用いたN-S方程式の数値計算コードは’ 1987年当時，東京工業大学工学部日野研究室でロ野． 神田・岸が開発した基本コードをベースに発展させた ものである．ここに記して深甚なる感謝の意を表した い． 参考文献 鳥羽良明（1966）：海塩粒子一大気と海洋との相互 作用の－要素として－，海と空，４１，３．４併合，ｐ､７１ ～１１８． 鳥羽良明・'1]中正昭（1967a）：塩答に|則する基礎 研究（第１報）海塩粒子の生成と陸上への輸送モデル， 京大防災研究所年報第10号Ｂ・ 鳥羽良明・H1中正昭（1967b）：塩害に閲する基礎 研究（第３報）海塩粒子の三次元分布の解析,京大防 災研究所年報第10号Ｂ・ 仲座栄二・津溌lll正光・照屋雅彦(1991）：大気環境 アメニテイの－要素としての飛塩最特性，第38回海岸 工学講演会論文集，ｐ､896-900． 仲座栄三・津粥山正光・照屋雅彦(1990）：沿岸微気 候の一要紫としての飛塩逓特性，環境システム研究 Ｖ０１．１８，ｐ､８８～９３． 村上和男・加藤一正・清水勝義・福Ｈ１興人・宮崎啓 示(1992）：砕波梢内における海塩粒子の発生に関する 現地調査，雛39回海岸工学講淡会論文染，ｐ､1046-１０５０．