1 / 21

回転系における潮流海底境界層の 乱流に関する数値的研究

日本地球惑星科学連合 大会 2007.05.24. 回転系における潮流海底境界層の 乱流に関する数値的研究. ○ 坂本圭、秋友和典　 ( 京都大学大学院・理学研究科 地球惑星科学専攻 ). 1 はじめに (1). 世界海洋の底層の多くを占める「 南極底層水 」 (Schmitz 1995) 南極陸棚上の重い海水 が、沈降していく際に 周囲の海水との混合 を通して変質し形成される (Fahrbach et al. 2001). 海面冷却. 混合. 混合. 南極陸棚. 南極底層水. 陸棚斜面. 外洋. 1 はじめに (2).

mira-buckley
Download Presentation

回転系における潮流海底境界層の 乱流に関する数値的研究

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 日本地球惑星科学連合 大会 2007.05.24 回転系における潮流海底境界層の乱流に関する数値的研究 ○坂本圭、秋友和典　 (京都大学大学院・理学研究科地球惑星科学専攻)

  2. 1 はじめに(1) 世界海洋の底層の多くを占める「南極底層水」(Schmitz 1995) 南極陸棚上の重い海水が、沈降していく際に周囲の海水との混合を通して変質し形成される(Fahrbach et al. 2001) 海面冷却 混合 混合 南極陸棚 南極底層水 陸棚斜面 外洋

  3. 1 はじめに(2) 他の海域の潮流海底境界層はせいぜい厚さ数十mなのに対し、コリオリ・パラメータfがM2潮の振動数σに一致する臨界緯度(74.5度)が近いため、数百mの厚い境界層が形成される(Furevik and Foldvik 1996) →不安定による混合は海底からはるか上方にまで及ぶ この混合過程において、潮流によって形成される海底境界層(潮流海底境界層)のシアー不安定が重要 (Foster et al. 1987) 南極陸棚上でのM2潮(半日周潮) 海面 ◇：観測 さらに、潮流構造の解析から、乱流粘性係数は300~1000cm2/sと見積もり(Nost 1994) 海底(水深580m) (Pereira et al. 2002)

  4. 1 はじめに(3)　境界層に関する先行研究 潮流の振動と地球の回転の効果が同程度となる臨界緯度付近の場合も含めて、両者が存在する下での、境界層の不安定力学については既に調べられている ○Aelbrecht et al. (1999): 水槽実験 ○Sakamoto and Akitomo (2006):数値実験 振動と回転の効果の比を示す「時間ロスビー数」Rot = σ/fに応じて、非回転系での振動流による境界層(ストークス層)の不安定や回転系での定常流による境界層(エクマン層)の不安定が現れる 一方、不安定の発達に伴って引き起こされる乱流の特性とその混合効果については明らかでない そこで本研究では、まず密度一様の下での、回転系における潮流海底境界層の乱流に関する数値実験を行う

  5. 1 はじめに(4)　本研究の内容 これまでの研究から、乱流エクマン層では、次のouter scaleでスケーリングを行うことで乱流が相似性を持つ(Coleman 1999)　時間：1/|f|　速さ：摩擦速度u* =(海底応力/密度)1/2　長さ：u*/|f| これを参考に　時間：T=1/|f+σ|(回転と振動を考慮)　速さ：摩擦速度u*　長さ：δ=u*/|f+σ|という新たなouter scaleを導入すると、潮流海底境界層における乱流の特性と混合効果が相似性を持つという結果が得られた 本報告： 2節:数値モデル3節:乱流特性の相似性4節:混合効果5節:まとめと課題

  6. 2 数値モデル：領域、支配方程式系 モデル領域 Lx×Ly×Hの矩形海領域 支配方程式系 回転系、密度一様、非圧縮、非静水圧、リジッド・リッド条件 変数を基本潮流場(、後述)と擾乱場( )に分ける ナヴィエ・ストークス方程式 渦粘性係数　ν =1cm2/s (等方)　、標準密度　ρ0=1.027g/cm3南半球を想定し f < 0 変数のチルダは有次元量であることを示す

  7. 2 境界条件、初期条件、差分 境界条件 海面：リジッド･リッド、非粘着 海底：粘着条件 水平：周期条件 初期場：微小擾乱　　積分期間：12潮流周期(いくつかのケースでは低分解能モデルで長期積分した後、本実験) 実験領域とグリッド間隔: 　 基本潮流場 の鉛直スケールHtide=(2ν/ |f+σ|)1/2で無次元化した値でLx=Ly=64, H=256 ⊿x=⊿y=0.125 ⊿z=0.02-10 (160グリッド) Htideと潮流振幅を用いて方程式を無次元化し実験を行うが、実験結果にはouter scale (T=1/|f+σ|,u*,δ=u*/|f+σ|)でスケーリングした無次元量(チルダなし)を示す。 e.g.

  8. 2 実験ケース、基本潮流場 時間ロスビー数Rotに対する乱流の依存性に注目潮流振幅は全て一定(8.53cm/s) M2潮に対する緯度67°53° 29° Ek ABCDEF St エクマン層 ストークス層 乱流特性の解明に主眼を置くため、現実より低いレイノルズ数 潮流ベクトル (utide, vtide) ケースEk ADSt x方向に、潮汐に伴う圧力勾配が働くとした場合の解析解を与える ｙ方向 x方向 ケースEk,Stを除いて、ケースA~Fの全てで、コリオリ力(f<0)により反時計回り

  9. 3 結果　渦運動エネルギーEKEの時間発展 統計的に定常解析に用いる (領域平均) ケース：Ek, A, B, CD, E, F, St(破線) Time (tidal cycle)

  10. 摩擦速度u*は大きく変化しない(潮流振幅の3.5～5.4％)摩擦速度u*は大きく変化しない(潮流振幅の3.5～5.4％) 長さスケールδは、fとσの絶対値が近づくRot～1で大きい(~450m) 3 乱流場, outer scale D Case A 実験終了時wの鉛直断面分布 Outer scale Ek ABCDEF St 摩擦速度時間スケール長さスケール

  11. 3 有次元における乱流特性 ケース：Ek, A, B, C, D, E, F, St(破線) 平均流 (水平・時間平均) 全応力(cm2/s2) (レイノルズ応力+モデル粘性) (m) 直交方向 潮流方向 潮流方向 直交方向 (Rot>1では直交方向の符号を反転させて表示) 有次元では、各ケースの鉛直分布は大きく異なる

  12. 3 乱流特性の相似性 (outer scale) ケース：Ek, A, B, C, D, E, F, St(破線) 平均流 全応力 z 直交方向 潮流方向 潮流方向 直交方向 (Rot>1では直交方向の符号を反転させて表示) 平均流、応力： ストークス層(St、破線)を除いて、outer scaleで無次元化すれば鉛直分布はエクマン層の結果にほぼ重なる →潮流海底境界層の乱流特性は、エクマン層と同じ“相似性”を持つ

  13. 3 相似性の理由 • Q.なぜ、エクマン層と同じ相似性が、新たなouter scale(fのかわりにf+σ)の導入によって現れる？ • 一般に、潮流ベクトル(utide, vtide)は2つの回転成分に分解できる(Makinson et al. 2006)南半球(f<0)では潮流ベクトルは反時計回り　→　反時計回り成分が卓越よって潮流は、一定の流速を保ち、反時計回りに向きを時々刻々変える流れと近似できる。 • この流れに従って回転する新たな座標系を考えると　潮流は定常流とみなせ、その座標系のコリオリ・パラメータはf+σとなる。 • よって、南半球における潮流海底境界層の乱流は、f+σで回転する座標系において定常流が形成する境界層、つまりエクマン層と相似となる。 反時計回り　＋　時計回り ※北半球では、時計回り成分が卓越するため、fのかわりにf-σが適切である。

  14. 4 混合効果：トレーサーによる見積もり St(破線)以外のケースでκapもほぼ相似形 　→緯度(f)、潮流周期(σ)、潮流振幅(u*)から乱流拡散係数を決定できる Ek, A, B, C, D, E, F, St(破線) 鉛直に線形な初期値を持つトレーサーの時間発展を計算する →「見かけの鉛直拡散係数」κapを評価する z C:トレーサー濃度 κap (outer scale)

  15. 4 乱流混合：有次元 緯度29°(ケースF)：　最大60cm2/s、　最大の半分となる高さ40m緯度53°(E)：　最大170cm2/s、高さ160m緯度67°(D):　最大600cm2/s、高さ350m Ek, A(Rot=0.5), B(0.8), C(0.95), D(1.05), E(1.2), F(2.0), St(破線) 振幅8.53cm/sのM2潮について　有次元でκapを示す Rotが1に近いほど、長さスケールδの増大に伴い、広い範囲で強い混合が起こる (m) (cm2/s) ケースDは、極海域における潮流観測からの乱流混合の見積もり(厚さ:～300m, 強さ:300～1000cm2/s)をよく説明する結果

  16. 5 まとめと課題 • 回転系(コリオリ・パラメータf)における潮流海底境界層(振動数σ)の乱流に関して、　　時間:T=1/|f+σ|速さ：摩擦速度u*　長さ：δ=u*/|f+σ|という「outer scale」でスケーリングすることで、ストークス層を除いて、平均流や応力といった乱流特性が相似性を持つことが分かった • 　　北半球では、時間：T=1/|f-σ|　長さ：δ=u*/|f-σ| • 混合効果についても、outer scaleによって相似性が現れた　○緯度(f)、潮流周期(σ)、潮流振幅(u*)の値から乱流拡散係数を決定できる　○現実の海洋へ適用すると、fとσが近づく臨界緯度で顕著な乱流混合 • 今後の課題： • 現実的な高レイノルズ数実験(数十万から数百万) • 安定成層の影響

  17. 4 乱流混合：安定成層 緯度67°(ケースD):　最大110cm2/s、高さ35m緯度53°(E)：　最大85cm2/s、高さ28m緯度29°(F)：　最大53cm2/s、高さ18m 現実的な安定成層の下でも、Rotが1に近いほど広い範囲(～2倍)で強い(～2倍)混合 Ek, A, B, C, D, E, F, St(破線) 高緯度海域：弱いながらも安定成層　→初期場に浮力振動数N=1.0×10-3s-1の密度成層を与えた追加実験 見かけの鉛直拡散係数κapを評価※ κapは時間変化する (m) (cm2/s) 実験開始から12~24潮流周期の期間

  18. 粘性係数を一定とした場合の潮流海底境界層の解析解粘性係数を一定とした場合の潮流海底境界層の解析解 振動数ωの潮流楕円 V(z,t)を反時計回り成分(振幅R+、初期位相φ+)と時計回り成分 (R-、φ-)に分解する。 それぞれの回転成分に対する境界層の厚さHtide+,Htide- ν,fは鉛直渦粘性係数、コリオリ・パラメータを示す。 Htide+ Htide- f > 0　→　潮流楕円は時計回り　　→　R-が支配的f < 0　→　潮流楕円は反時計回り　→　R+が支配的 Prandle (1982)

  19. 議論：Rot~1における慣性波 群速度の見積もり 粘性の時間スケール 波の到達距離 境界層では波長δの擾乱が発生 Rot~1　→　群速度大　→　波長δの波がエネルギーを維持したまま上層へその他　→　群速度小　→　δより波長が大きい波、エネルギーは小さい l:波長

  20. 3 慣性波　ケースC 慣性波群速度での、エネルギーの上方への輸送と反射・Rot~1では波長δの慣性波が盛んに上方へと伝播・他のケース(St以外)では、慣性波の波長は大きくエネルギーは小さい　→上層のqとlに差 渦運動エネルギー水平平均の、t-zダイアグラム 海面 EKE z 海底 t

  21. 4 乱流混合：l×q 乱流理論：　乱流拡散の強さ ∝ l×q Ek, A, B, C, D, St(点線) l×q：オーダーとしては同程度同じような分布だが、 波の影響によって、相似性があるとは言い難い z しかし、砕波が起こらず波の影響が無視できれば、引き起こされる乱流混合も相似性を持つはず 　→トレーサーによる混合効果の見積もり l×q

More Related