潮流によって形成される海底境界層の不安定とその混合効果

Slides:



Advertisements
Similar presentations
海洋流体力学 2014 海洋流体力学とは、海洋に関する流体力学。本講義では、 海洋のみならず、大気も含めた地球流体力学について学ぶ。 Fluid Dynamics( 流体力学 ) Geophysical Fluid Dynamics (地球流体力学) 目標 海洋・大気大循環のイメージを描けるようにする。
Advertisements

熊野灘海流予測システム開発 進捗状況報告 (株)三菱総合研究所. 熊野灘海流予測システム 内容 – 熊野灘で作業中である、地球深部探査船「ち きゅう」のために海流予測を行う 黒潮の変動を数キロメートルのオーダーで予測 –JCOPE をネスティング » 日本近海 1/36 度モデル(同化あり)
宇宙ジェット形成シミュレー ションの 可視化 宇宙物理学研究室 木村佳史 03S2015Z. 発表の流れ 1. 本研究の概要・目的・動機 2. モデルの仮定・設定と基礎方程式 3. シンクロトロン放射 1. 放射係数 2. 吸収係数 4. 輻射輸送方程式 5. 結果 6. まとめと今後の発展.
Determining Optical Flow. はじめに オプティカルフローとは画像内の明る さのパターンの動きの見かけの速さの 分布 オプティカルフローは物体の動きの よって変化するため、オプティカルフ ローより速度に関する情報を得ること ができる.
YohkohからSolar-Bに向けての粒子加速
自己重力多体系の 1次元シミュレーション 物理学科4年 宇宙物理学研究室  丸山典宏.
数値気象モデルCReSSの計算結果と 観測結果の比較および検討
コリオリ力の復習資料 見延 庄士郎(海洋気候物理学研究室)
プロセス制御工学 3.伝達関数と過渡応答 京都大学  加納 学.
スペクトル法による数値計算の原理 -一次元線形・非線形移流問題の場合-
潮流によって形成される海底境界層の不安定とその混合効果
加藤真理子1、藤本正樹2、井田茂1 1) 東京工業大学 2) JAXA/ISAS
成層圏突然昇温の 再現実験に向けて 佐伯 拓郎 神戸大学 理学部 地球惑星科学科 4 回生 地球および惑星大気科学研究室.
資源の空間的不均一性がプランクトン群集の共存に与える影響: 格子モデルシミュレーションによる予測
大阪工業大学 情報科学部 情報システム学科 宇宙物理研究室 B 木村悠哉
渦位(Potential Vorticity)と角運動量の保存
北海道大学 理学部 地球科学科 惑星物理学研究室 4年 馬場 健聡
相対論的輻射流体力学における 速度依存変動エディントン因子 Velocity-Dependent Eddington Factor in Relativistic Photohydrodynamics 福江 純@大阪教育大学.
風成海洋大循環 (準地衡流渦位方程式+エクマン層の力学)
北海道大学大学院理学研究科地球惑星科学専攻 地球流体力学研究室 M1 山田 由貴子
輻射優勢円盤のMHD数値実験 千葉大学宇宙物理学研究室 M2 松尾 圭 Thu.
反応性流体力学特論  -燃焼流れの力学- 燃焼の流体力学 4/22,13 燃焼の熱力学 5/13.
流体のラグランジアンカオスとカオス混合 1.ラグランジアンカオス 定常流や時間周期流のような層流の下での流体の微小部分のカオス的運動
輻射圧駆動風の臨界点ついて 非相対論的領域 Radiatively Driven Spherical Wind Critical Points and Curves Nonrelativistic Regime 福江 純@大阪教育大学.
Lorenz modelにおける 挙動とそのカオス性
原子核物理学 第4講 原子核の液滴模型.
バングラデシュにおける対流活動と局地風に関する研究
圧力発展格子ボルツマン法による大規模気液二相流GPUコードの開発 ならびに多孔体浸潤液滴シミュレーション
磁気リコネクション研究 における MHDシミュレーションの限界と解析的研究の展望
半無限領域のスペクトル法による竜巻を模した渦の数値実験に向けた研究開発
正規分布確率密度関数.
Fermi Bubble と銀河中心の巨大構造
磁気リコネクション (Craig-Henton解の安定性) ~シミュレーションサマースクール@千葉大より~
磯部洋明 京都大学花山天文台 波動加熱勉強会 2004年2月23日
YT2003 論文紹介 荻原弘尭.
流体の粘性項を 気体分子運動論の助けを借りて、 直感的に理解する方法
川崎浩司:沿岸域工学,コロナ社 第2章(pp.12-22)
論文紹介 坂本圭 Finite-amplitude crossflow vortices, secondary instability and transition in the rotating-disk boundary layer Pier, B. J. Fluid Mech.
燃焼の流体力学 4/22 燃焼の熱力学 5/13 燃焼流れの数値解析 5/22
メンバー 梶川知宏 加藤直人 ロッケンバッハ怜 指導教員 藤田俊明
回転系において潮流が形成する 海底境界層の不安定
講義ノート(ppt)は上記web siteで取得可 #但し、前日に準備すると思われるのであまり早々と印刷しない方が身の為
回転系における潮流海底境界層の 乱流に関する数値的研究
(昨年度のオープンコースウェア) 10/17 組み合わせと確率 10/24 確率変数と確率分布 10/31 代表的な確率分布
開放端磁場における低温プラズマジェットに関する研究
渦位(Potential Vorticity)と角運動量の保存
2009年秋の北極海ラジオゾンデ観測によって観測された 大気の順圧不安定とメソ渦列
4章 開水路における不等流(2) 漸変流 4-1漸変流とは ① 断面形状や底面形状が緩やかに変わる流れ。
連続体とは 連続体(continuum) 密度*が連続関数として定義できる場合
音声分析 フーリエ解析の定性的理解のために.
Large-eddy simulations of the wind-induced turbulent Ekman layer
随伴解析を用いた物体表面形状最適化による抵抗低減
Iida, O., N. Kasagi and Y. Nagano
速度ポテンシャルと 流線関数を ベクトルで理解する方法
Chapter 26 Steady-State Molecular Diffusion
バリオン音響振動で探る ダークエネルギー ~非線形成長と赤方偏移歪みの影響~
九州大学 猿渡元彬 共同研究者 橋本正章 (九州大学)、江里口良治(東京大学)、固武慶 (国立 天文台)、山田章一(早稲田理工)
竜巻状渦を伴う準定常的なスーパーセルの再現に成功
潮流によって形成される海底境界層の不安定とその混合効果
潮流によって形成される海底境界層の不安定とその混合効果
落下水膜の振動特性に関する実験的研究 3m 理工学研究科   中村 亮.
対象:せん断補強筋があるRCはり(約75万要素)
第 5 章 :周波数応答 5.1 周波数応答と伝達関数 周波数伝達関数,ゲイン,位相 キーワード : 5.2 ベクトル軌跡 ベクトル軌跡
卒論中間発表 2001/12/21 赤道の波動力学の基礎 北海道大学理学部 地球科学科 4年 山田 由貴子.
磁気リコネクションによる Alfven波の発生
北大MMCセミナー 第23回 Date:2014年3月6日(木) 16:30~18:00 ※通常と曜日が異なります
従来研究 本研究 結果 南極大型大気レーダーPANSYで観測された大気重力波の数値モデル再現実験による力学特性の解明
北大MMCセミナー 第94回 附属社会創造数学センター主催 Date: 2019年1月25日(金) 16:30~18:00
Presentation transcript:

潮流によって形成される海底境界層の不安定とその混合効果 ○坂本圭、秋友和典  (京都大学大学院・理学研究科)

1 はじめに(1) 潮流海底境界層 潮流によって形成される海底境界層(潮流海底境界層)の流速構造 定常エクマン層と相似 鉛直スケール 1 はじめに(1) 潮流海底境界層 Fang and Ichiye 1983 潮流によって形成される海底境界層(潮流海底境界層)の流速構造   定常エクマン層と相似   鉛直スケール σ:潮流振動数(反時計回りの潮流楕円を正とする) f:コリオリ・パラメータ、ν:粘性係数 流速の鉛直プロファイル(1時間毎)  ν=50cm2/s U V 海底 潮流周期12時間 慣性周期8時間 潮流周期12時間 慣性周期12.5時間 定常エクマン層 慣性周期12.5時間

1 はじめに(2) 定常エクマン層不安定 変曲点不安定(タイプI) Uの変曲点から直接u’へ擾乱エネルギー コリオリ型不安定(タイプII) 1 はじめに(2) 定常エクマン層不安定 Kaylor and Faller 1972 変曲点不安定(タイプI)  Uの変曲点から直接u’へ擾乱エネルギー コリオリ型不安定(タイプII)  V(内部流方向)のシアーからv’へエネルギーが供給され、その後コリオリ力によってu’へ エクマン層 流速鉛直構造 U V 内部定常流  (y方向) 海底 エクマン螺旋擾乱(u’,v’,w’) 安定性はレイノルズ数に依存  ~100   : 安定  100~150 : コリオリ型不安定  150~   : 変曲点不安定が卓越

1 はじめに(3) ストークス層不安定 コリオリ力が0 → 潮流海底境界層はストークス層 ストークス層の流速鉛直プロファイル(1時間毎) 1 はじめに(3) ストークス層不安定 Akhavan et al. 1990  コリオリ力が0 → 潮流海底境界層はストークス層 ストークス層の流速鉛直プロファイル(1時間毎) 振動周期12時間 振幅40cm/s 流速プロファイルに様々な変曲点 →振動流振幅がある値を超えれば、ストークス層も変曲点不安定 海底

1 はじめに(4) 目的 潮流海底境界層も、他の粘性境界層と同様に、潮流流速が大きくなれば不安定となると考えられる。 1 はじめに(4) 目的  潮流海底境界層も、他の粘性境界層と同様に、潮流流速が大きくなれば不安定となると考えられる。  しかし、その不安定性に関する研究はほとんどない。  そこで本研究では、鉛直2次元数値実験によって、潮流海底境界層の不安定を再現しその力学と混合効果について調べる。  潮流の向きは時々刻々変化するため、不安定は本来3次元的な構造を持つと考えられるが、今回はある鉛直平面で発達するロールに注目する。 1.密度一様実験  相似な流速構造を持つ定常エクマン層の不安定との違いを明らかにする(3章)。 2.成層実験  海洋陸棚上においてどの程度の混合効果を持つかを見積もる(4章)。

2 モデル領域 鉛直2次元、 14km×500mの矩形海。

2 支配方程式系 運動方程式 連続の式 移流拡散方程式 鉛直2次元、非圧縮、非静水圧、ブシネスク近似、リジッド・リッド条件。 数値計算には渦度と流線関数を用いる。 渦粘性・拡散係数 ν =50cm2/s,κ=5cm2/s 重力加速度 g=980cm/s2 標準密度 ρ0=1.027g/cm3

2 境界条件、初期条件 海面:リジッド・リッド、no-flux 海底:粘着条件、no-flux 潮流として左右境界で時間振動流 U0(t)=-Utidecos( 2π×t / Ttide ) 振動周期Ttide =12時間 初期条件:密度一様(3章)または線形成層(4章)、静止状態 積分期間:40日間 実験の制御パラメータは、  1:慣性周期Ti (=2π/f)  2:振動流振幅Utide コリオリ力によって、内部の流れは潮流楕円を描く。   (軸:Utide、 Utide×(Ttide / Ti) )

3 実験結果 基本流Uの鉛直プロファイル(1時間毎) 慣性周期Ti =8時間 Utide =30cm/s Ti =12.5時間 海底 慣性周期Ti =8時間 Utide =30cm/s Ti =12.5時間 Utide =20cm/s Ti =20時間 Utide =50cm/s 擾乱に伴う流線関数 (等値線間隔:200cm2/s) 海底

3 鉛直2次元擾乱運動エネルギー方程式 基本流Uのシアーからu’へ →変曲点不安定 →コリオリ型不安定 コリオリ力によってv’からu’へ 粘性による消失 再分配 左右境界での流出入 変曲点不安定 コリオリ型不安定 (ストークス層型) 擾乱成長期 12時間平均値 海底 Ti =8時間 Utide =30cm/s Ti =12.5時間 Utide =20cm/s Ti =20時間 Utide =50cm/s

3 全ケーススタディ 様々な慣性周期Ti、潮流振幅Utideにおける不安定のタイプ 慣性周期が潮流周期の 0.9倍以下: 変曲点不安定  変曲点不安定 1.0-1.1倍:  コリオリ型不安定 1.2倍以上:  ストークス層不安定 (Ttide=12時間) ×安定  □エクマン層変曲点不安定 △コリオリ型不安定  ◇変曲点不安定とコリオリ型が同程度  ※ストークス層不安定

3 Ti < Ttide エクマン層変曲点不安定  慣性周期 < 潮流周期 (点線:レイノルズ数) ×安定 □変曲点不安定 △コリオリ型不安定  ◇変曲点不安定とコリオリ型が同程度 定常 エクマン層 ただし、潮流が一方向を向いている間に不安定は成長しなければならない  →成長率の小さいコリオリ型不安定は現れない (Ttide=12時間) 慣性周期が短く、潮流が反転する前にエクマン層が発達  →定常エクマン層の不安定に帰着

3 Ti ~ Ttide エクマン層コリオリ型不安定  慣性周期:潮流周期の1.0-1.1倍 (点線:レイノルズ数) 定常 エクマン層 ×安定 □変曲点不安定 △コリオリ型不安定  ◇変曲点不安定とコリオリ型が同程度 ※ストークス層不安定 広い範囲のUtideでコリオリ型不安定が支配的 (Ttide=12時間) Ti~Ttide → Hの増大  1.粘性の寄与の低下  2.コリオリ型の寄与が変曲点不安定より相対的に上昇 擾乱エネルギー各項のスケール     鉛直スケール:H= 変曲点不安定: U/H コリオリ型不安定: f 粘性による消失:ν/H2 (U:流速スケール, f:コリオリ・パラメータ, ν:粘性係数)

3 Ti > Ttide ストークス層不安定  慣性周期 > 潮流周期 ×安定 △コリオリ型不安定  ※ストークス層不安定 不安定となる潮流振幅:完全なストークス層と同じ  擾乱の流れ場などもほぼ同じ →ストークス層の不安定に帰着できる Ti→∞ (f=0) ストークス層 (Ttide=12時間)

4 極と中緯度海域における混合効果 線形成層実験 (潮流周期:12時間、Utide=50cm/s) ケースA:慣性周期Ti=12.5時間(74°S)、成層なし ケースB:Ti=12.5h(74°S)、浮力振動数N=0.001s-1 ケースC:Ti=18.8h(41°S)、N=0.001s-1 ケースD:Ti=18.8h(41°S)、N=0.01s-1 A,B C,D 実験結果 ケースA,B:変曲点とコリオリ型不安定が同程度        ケースC,D:ストークス層不安定

4 極と中緯度海域における混合効果 トレーサー濃度鉛直分布(52日目水平平均) 初期濃度 海底~25m:1.0、25m以上:0.0  (講演要旨より鉛直分解能を上げた実験結果) 海面 A:Ti=12.5h(74°S)、N=0 B:Ti=12.5h(74°S)、N=0.001s-1 海底から約130mまでよく混合  見かけの拡散係数:168cm2/s C:Ti=18.8h(41°S)、N=0.001s-1 D:Ti=18.8h(41°S)、N=0.01s-1 海底 ケースCとD:擾乱流速は10-4cm/s程度  トレーサーの変化はモデルに与えた拡散(拡散係数5cm2/s)によるもの トレーサー濃度

5 まとめと課題 潮流海底境界層の不安定力学 慣性周期が潮流周期の0.9倍以下: エクマン層変曲点不安定  慣性周期が潮流周期の0.9倍以下: エクマン層変曲点不安定          1.0から1.1倍:  コリオリ型不安定          1.2倍以上  : ストークス層不安定 半日周潮の混合効果  慣性周期12.5時間(74°S):   海底から約130mまで一様化、混合効果は168cm2/s    → 極域陸棚上で形成される高密度水に影響  18.8時間(41°S):不安定は非常に弱く混合効果はほとんどない 課題:  1.成層によって不安定が抑制される効果  2.3次元実験   擾乱はそもそも3次元構造を持つはずである。   潮流楕円に関する制限を外す。