4. ロスビー波と 1.5 層モデル見延庄士郎 ( 海洋気候物理学研究室 ) 予習ビデオ : NASA JPL の海面高度偏差でエルニーニョを見るページ. 赤道ケルビン波が見えるか? の動画 ht

Size: px
Start display at page:

Download "4. ロスビー波と 1.5 層モデル見延庄士郎 ( 海洋気候物理学研究室 ) 予習ビデオ : NASA JPL の海面高度偏差でエルニーニョを見るページ. 赤道ケルビン波が見えるか? の動画 ht"

Transcription

1 4. ロスビー波と 1.5 層モデル見延庄士郎 ( 海洋気候物理学研究室 ) 予習ビデオ : NASA JPL の海面高度偏差でエルニーニョを見るページ. 赤道ケルビン波が見えるか? hps://sealevel.jpl.nasa.gov/elnino015/inde.hml の動画 hps://sealevel.jpl.nasa.gov/elnino015/1997vs015-animaed-800.gif 大気のロスビー波 :hps:// 室内実験でロスビー波 :hps:// 式番号について : 異なるスライドで式を引用するときには,(3.1) などと ( 章番号. 式番号 ) とする. 一つのスライドの中, あるいはすぐ次のスライドだけで使うなら 1 などとする. 第 3 回のまとめ おおざっぱにいうと, 時間スケールが 1 日よりも短い現象では, コリオリ項を無視することができる. その場合に重要となる波は, ( 非回転 ) 重力波である. 層厚がHの1 層モデルでは重力波の伝播速度は である. 上の波が1/fの時間で伝播する距離を, ロスビーの変形半径と言い, これよりも大きいスケールでは一般に地衡流バランスが成り立つ ( 摩擦無しの場合 ) 岸があることで生じる波が, 沿岸ケルビン波で,1 層モデルの伝播速度は で, 岸を右に見る方向に伝播する. それと近い関係にあるのが, 赤道ケルビン波で, 1 層モデルの伝播速度は で赤道に沿って東向きに伝播する. 1 本日の目的 今日は, ロスビー波, そして 1.5 層モデルを導入する. 地球流体の特徴は, コリオリ力を受けることで, しかもコリオリパラメータ (f=ω sinθ) が緯度の関数であることから明白なように, コリオリ力の効果が緯度によって違う. この緯度変化の効果から生じる つの波がある. その一つが赤道を東に進む赤道ケルビン波で, もう一つが西に進むロスビー波である. またエルニーニョを含めて多くの風が海を駆動する現象で, 海洋の運動は上層に限られるので, 海が上から下まで動くという 1 層モデルでは都合が悪い. この問題を解決する, 上層だけ動く 1.5 層モデルを導入しよう. 海に見えるロスビー波 海面高度 ( 海の表面のでこぼこ ) 3 4

2 赤道で東へ伝播 赤道ケルビン波 赤道外では西へ伝播 ロスビー波 5N と 1N で, 後者の波の速度が速いのは 1N である. ロスビー波 5 6 地衡流の収束 発散?(f 一定の場合 ) 1) 初期水位 ) 流れは? ロスビー波のメカニズム :f が北で大きい場合 1) 初期水位 ) 流れは? 北 東 赤道 北 東 赤道 3) 収束発散? 4) 移動方向? 3) 収束発散? 4) 移動方向? 収束 発散 収束 発散 収束発散はゼロ, 伝播はしない 伝播方向は, 西向き 7 8

3 定量的に理解するために式を使おう. 前回学んだ, 摩擦 風応力を考えない浅水方程式系は, 東向き運動方程式 北向き運動方程式 連続の式 u fv = gh v + fu = ghy h = H( u + v ) 簡素化のために, 運動方程式の時間変化項を無視 ( 消去 ) したい 前回のスケーリングの議論を踏まえると, 時間変化項を無視できるのは, 以下の二つの等価な場合である 時間スケールが 1 日よりも十分長い 空間スケールがロスビーの変形半径よりも十分大きい 今回導出するロスビー波は, 長波ロスビー波と呼ばれる 1 3 長波ロスビー波の伝播速度 東向き運動方程式 北向き運動方程式連続の式 d f /d y β 0 に注意 fv = gh 1 + fu = gh y h = H( u + v ) 3 クロス微分の定跡 : / - 1/ yで, 必ず圧力項を 消すことができる. f ( u + vy ) + βv = 0 β これを3に代入して h = Hv f これに1を代入してvを消すと gh h β h f これは一階線形偏微分方程式 ( 伝播方程式 ). 伝播方向は西. 伝播速度は gh β = β R O f = 0 (4.1) 前回登場のロスビーの変形半径 R O を使うと 9 10 Nino3.4 と 00 hpa 等圧面高度の相関 第一回スライド再掲 エルニーニョの影響は中高緯度に及ぶ 11 大気のロスビー波 中緯度の大気では, 東西方向に地球を一周する西風 ( 東向き ) の成分が強い. 陸があるので, 海流は地球を東西方向に一周できない. その強い西風と, 静止流体では西に伝播しようとするロスビー波の速度が逆方向で大きさが等しいと, 波の伝播が風に押しとどめられ, 定在ロスビー波が生じる. 実際には, ロスビー波は様々な速度 ( とそれに対応した波長 ) を持つので, 定在する波が選択的に励起される. 定在ロスビー波の存在には西風が不可欠なので, 貿易風などの東風領域では, ロスビー波は存在できない. 定在ロスビー波はエネルギーを必ず東に伝えるので, エルニーニョに伴う大気のロスビー波はアメリカ大陸を直撃するけれど, 日本には来ない. 日本へのエルニーニョの影響は間接的なメカニズムによる. 1

4 テレ コネクション 温暖 寒冷 テレ コネクション ( 遠隔結合 ) ある場所の気候状態が離れた地点の気候状態に影響 代表格が, エルニーニョに伴って生じる太平洋 / 北米パターン ( 左図 ) この図から, エルニーニョの影響が最も強いのは, 北米 南米 日本の中では, 北米である ロスビー波の速度を計算しよう Q. 水深 4000 m, 重力加速度 10 m/s, として 5N でのロスビー波の伝搬速度を求めよ. なお,f=Ω sinθ = s -1, β=(ω/a) cosθ = m -1 s -1 を用いよ (a は地球半径 ). また, スライド 5 の 5N の伝搬速度を, 図から求めよ. この際, 経度 1 度当たりの距離が,5N で 100 km であることを用いる. 両者の結果は一致するといえるか? A ロスビー波の伝搬速度は, β gh / f =0 m/s. 一方, 図の伝搬は経度 10 度を 6 年で伝播するので, c=10*100*1000/(6*365*4*60*60)=0.063 m/s となって, 全然合わない 海洋成層と運動の鉛直構造 混合層 0 m 0 上層 温度躍層 層モデル m 0 C 10 C 0 C WOA94,SODA7 海洋の運動の多くは, 上から下まで一様 (1 層 ) ではなく, 上 ( 上層 ) だけ速い. 層のうち上層だけ動く,1.5 層モデルが有効 N 45N C 16

5 1.5 層モデル 1/3 密度が少し異なる, 層流体を考える. 層目の厚さが, 無限だと仮定する. 層目は動かない. この場合, 上層の流速をu, v とすると, 以下の浅水方程式と同じ形の u fv = g v + fu = g h + H ( u + v ) = 0 流体表面 境界面 方程式系が得られることを示そう. ここで,g は低減重力と呼ばれる量である. H y η η 1 1 g ' g ρ = ρ (4.4) ρ 0 (4.a) (4. b) (4.3) h p, ρ z=0 u,v, p 1, ρ 1 上層 深層 層モデル /3 上層の ( 線形 ) 運動方程式は 1 層モデルと一緒で H η u fv = gη 1a 1 1 ab v + fu = gη 1b 境界面 1y 一方連続の式は明らかに次式が成り立つ h + H ( u + v ) = 0, h = H + η η したがって, η 1 を h のみで表せればよい. この関係を出すのに, 層目が動かないという条件を使う. まず p は ( 水平変化成分のみに着目して ) p = p1 + ρ0g ' η = ρ0gη1 + ρ0g ' η なので, これがゼロとなるには η = η 1 g / g ', gはg よりもずっと大きいので したがって, 流体表面 1 η >> η 1 1 η 1 h p, ρ z=0 u,v, p 1, ρ 1 上層 深層 下層が動かないので,p はゼロ. 四角に η 1 か η を入れよう h η, η = hg '/ g 層モデル 3/3 一つ前スライドの運動方程式 1abに,と3を代入すると u fv = g v + fu = g となって求める運動方程式 (4.ab) が得られた.( 導出完了 ) 1.5 層モデル方程式系は低減重力モデルとも呼ばれる. 1 層浅水方程式の, 重力加速度が低減重力に, 全水深が上層の厚さに, 置き換わっただけ. y 19 1 層モデル vs1.5 層モデル 1 層モデルの層厚を H=4000 m, 重力加速度を g=10 m/s とし,1.5 層モデルの層厚さを H=300 m, 低減重力を g =0.03 m/s とする. 共通のパラメータとして, f=ω sinθ = s -1, b=(ω/a) cosθ =10-11 m -1 s -1 とする.(a は地球半径 ) Q1 1.5 層モデルで重力波の速度を表す式を書き, 速度を求めよ. A1 1.5 層モデルの重力波の速度は, =3 m/s Q 1.5 層モデルで長波ロスビー波の速度を表す式を書き, 速度を求めよ. この速度は,18 枚目のスライドで求めたロスビー波速度の観測地と一致すると言えるか? A1 1.5 層モデルの長波ロスビー波の速度は, = /( ) = 9/ =0.05 m/s 0

6 第 4 回のまとめ 1/ 赤道に沿って東に伝播する波が赤道ケルビン波, 1 層モデルの伝播速度は. 岸に沿うのは沿岸ケルビン波. 長波ロスビー波は, 西に伝播し, その速度はロスビーの変形半径をR O とすると, β R O である. ロスビーの変形半径は, 重力波が1 日のスケールで伝播できる距離を表し,1 層モデルではである. 流体が動く上層と, 動かない無限の厚さを持つ 層からなるモデルを,1.5 層モデルという. このモデルでの重力波の伝播速度は低減重力をで表すと, である. 海面水位の凹凸が赤道外で西へ伝播する様子は,1.5 層モデルのロスビー波でよく説明できる. 大気のテレコネクションを担う波は定在ロスビー波である. 第 4 回のまとめ / 1.5 層モデルは1 層浅水方程式の, 重力加速度が低減重力に, 全水深が上層の厚さに, 置き換わったものであり, 外力 摩擦がない場合の式は, 以下のとおりである. u fv = g 低減重力は, v + fu = g h + H ( u + v ) = 0 ρ 1 g ' = g ρ ρ (4.4) 0 y (4.a) (4. b) (4.3) 1 おまけ : 主要英単語 ( テストには出しません ) 波 重力波 :graviy wave 赤道ケルビン波 : equaorial Kelvin wave 沿岸ケルビン波 :coasal Kelvin wave ロスビー波 :Rossby wave 方程式 運動方程式 :momenum equaion 連続の式 :coninuiy equaion 浅水方程式 :shallow waer equaions 1.5 層モデル :1.5-layer model または reduced graviy model 定数 コリオリパラメータ :Coriolis parameer 重力加速度 :graviy acceleraion 低減重力 :reduced graviy ロスビーの変形半径 :Rossby s deformaion radius 現象 地衡風 ( 地衡流 ):geosrophic wind (geosrophic flow) 伝播 :propagaion 3

4. ロスビー波と 1.5 層モデル見延庄士郎 ( 海洋気候物理学研究室 ) 予習ビデオ : NASA JPL の海面高度偏差でエルニーニョを見るページ. 赤道ケルビン波が見えるか? の動画 ht

4. ロスビー波と 1.5 層モデル見延庄士郎 ( 海洋気候物理学研究室 ) 予習ビデオ : NASA JPL の海面高度偏差でエルニーニョを見るページ. 赤道ケルビン波が見えるか?   の動画 ht 4. ロスビー波と 1.5 層モデル見延庄士郎 ( 海洋気候物理学研究室 ) 予習ビデオ : NASA JPL の海面高度偏差でエルニーニョを見るページ. 赤道ケルビン波が見えるか? hps://sealevel.jpl.nasa.gov/elnino015/inde.hml の動画 hps://sealevel.jpl.nasa.gov/elnino015/1997vs015-animaed-800.gif

More information

運動方程式の基本 座標系と変数を導入 (u,v) ニュートンの第一法則 力 = 質量 加速度 大気や海洋に加わる力を, 思いつくだけ挙げてみよう 重力, 圧力傾度力, コリオリ力, 摩擦力 水平方向に働く力に下線をつけよう. したがって水平方向の運動方程式は 質量 水平加速度 = コリオリ力 + 圧

運動方程式の基本 座標系と変数を導入 (u,v) ニュートンの第一法則 力 = 質量 加速度 大気や海洋に加わる力を, 思いつくだけ挙げてみよう 重力, 圧力傾度力, コリオリ力, 摩擦力 水平方向に働く力に下線をつけよう. したがって水平方向の運動方程式は 質量 水平加速度 = コリオリ力 + 圧 2. 潜水方程式系の導出 見延庄士郎 ( 海洋気候物理学研究室 ) [email protected] 第 1 回まとめ 1/2 二つの変数の関係の強さを表す統計量は相関であり, 最小値は -1, 最大値は +1, 無相関は である. 過去数十年間の ( 気象庁は 3 年 ) 月ごとの平均値を, 月平均データの平年値または気候値という. 観測値から平年値を引いたものが, 偏差である.

More information

運動方程式の基本 ニュートンの第一法則 力 = 質量 加速度 大気や海洋に加わる力を, 思いつくだけ挙げてみよう 重力, 圧力傾度力, コリオリ力, 摩擦力 水平方向に働く力に下線をつけよう. したがって水平方向の運動方程式は 質量 水平加速度 = コリオリ力 + 圧力傾度力 + 摩擦力 流体の運動

運動方程式の基本 ニュートンの第一法則 力 = 質量 加速度 大気や海洋に加わる力を, 思いつくだけ挙げてみよう 重力, 圧力傾度力, コリオリ力, 摩擦力 水平方向に働く力に下線をつけよう. したがって水平方向の運動方程式は 質量 水平加速度 = コリオリ力 + 圧力傾度力 + 摩擦力 流体の運動 2. 浅水方程式系の導出 見延庄士郎 ( 海洋気候物理学研究室 ) 第 1 回まとめ 1/2 二つの変数の関係の強さを表す統計量は相関であり, 最小値は -1, 最大値は +1, 無相関は である. 過去数十年間の ( 気象庁は 3 年 ) 月ごとの平均値を, 月平均データの平年値または気候値という. 観測値から平年値を引いたものが, 偏差である. 連続する n 個のデータを平均して, 中央のデータの値に置き換える平滑化が,

More information

第五回まとめ 2/2 東岸域では, 北風がエクマン流を通じて湧昇をもたらす. これを沿岸湧昇と呼ぶ. 沿岸湧昇域での, 局所的な大気海洋結合変動現象であるカリフォルニアニーニョ, ニンガルーニーニョなどが発見されている. エクマン流と地衡流の関係の仮説 前回学んだエクマン流が, どう地衡流と関係する

第五回まとめ 2/2 東岸域では, 北風がエクマン流を通じて湧昇をもたらす. これを沿岸湧昇と呼ぶ. 沿岸湧昇域での, 局所的な大気海洋結合変動現象であるカリフォルニアニーニョ, ニンガルーニーニョなどが発見されている. エクマン流と地衡流の関係の仮説 前回学んだエクマン流が, どう地衡流と関係する 6. 海洋風成循環と海面高度 見延庄士郎 ( 海洋気候物理学研究室 ) 復習課題 http://www.sci.hokudai.ac.jp/~minobe/class/po_clm/05_review_ekman.pdf を印刷解答の上, 授業冒頭で提出してください. 冒頭のみ提出を受け付けます. 東岸域でのエクマン流 Q. 大洋の東岸域では, 岸に沿う風が吹くと, エクマン流が岸に直交するので,

More information

浅水方程式 順圧であるためには, 静水圧近似が必要 Dw Dt + コリオリ力 = 1 p + 粘性 g ρ z w が u, v に比べて小さい 運動の水平距離に対して水深が浅い 浅水 海は深いが, 水平はさらに広い 最大 1 万 km 浅水方程式 : u, v, の式 水平 2 次元の解 D D

浅水方程式 順圧であるためには, 静水圧近似が必要 Dw Dt + コリオリ力 = 1 p + 粘性 g ρ z w が u, v に比べて小さい 運動の水平距離に対して水深が浅い 浅水 海は深いが, 水平はさらに広い 最大 1 万 km 浅水方程式 : u, v, の式 水平 2 次元の解 D D 流体地球科学第 11 回 東京大学大気海洋研究所准教授藤尾伸三 ttp://ovd.aori.u-tokyo.ac.jp/ujio/2015ciba/ [email protected] 2016/1/8 順圧流の運動方程式 流体の密度が一様ならば, 圧力 静水圧 の水平勾配は鉛直一様 海面の高さによる水平圧力勾配のみ ηx,y px, y, z = ρ g dz = ρgη z p x

More information

パソコンシミュレータの現状

パソコンシミュレータの現状 第 2 章微分 偏微分, 写像 豊橋技術科学大学森謙一郎 2. 連続関数と微分 工学において物理現象を支配する方程式は微分方程式で表されていることが多く, 有限要素法も微分方程式を解く数値解析法であり, 定式化においては微分 積分が一般的に用いられており. 数学の基礎知識が必要になる. 図 2. に示すように, 微分は連続な関数 f() の傾きを求めることであり, 微小な に対して傾きを表し, を無限に

More information

2 図微小要素の流体の流入出 方向の断面の流体の流入出の収支断面 Ⅰ から微小要素に流入出する流体の流量 Q 断面 Ⅰ は 以下のように定式化できる Q 断面 Ⅰ 流量 密度 流速 断面 Ⅰ の面積 微小要素の断面 Ⅰ から だけ移動した断面 Ⅱ を流入出する流体の流量 Q 断面 Ⅱ は以下のように

2 図微小要素の流体の流入出 方向の断面の流体の流入出の収支断面 Ⅰ から微小要素に流入出する流体の流量 Q 断面 Ⅰ は 以下のように定式化できる Q 断面 Ⅰ 流量 密度 流速 断面 Ⅰ の面積 微小要素の断面 Ⅰ から だけ移動した断面 Ⅱ を流入出する流体の流量 Q 断面 Ⅱ は以下のように 3 章 Web に Link 解説 連続式 微分表示 の誘導.64 *4. 連続式連続式は ある領域の内部にある流体の質量の収支が その表面からの流入出の合計と等しくなることを定式化したものであり 流体における質量保存則を示したものである 2. 連続式 微分表示 の誘導図のような微小要素 コントロールボリューム の領域内の流体の増減と外部からの流体の流入出を考えることで定式化できる 微小要素 流入

More information

実践海洋情報論 57 ここで mは流体の質量 (kg) u は流体の速度 (m sec -1 ) である しかし 緯度 φにおける角速度を直観的に理解することが困難である そこで 地球儀を北極上空から観察しよう 北極点に十字のマークを置くと 地球儀を反時計回りに回転させると 極点上の十字のマークも反

実践海洋情報論 57 ここで mは流体の質量 (kg) u は流体の速度 (m sec -1 ) である しかし 緯度 φにおける角速度を直観的に理解することが困難である そこで 地球儀を北極上空から観察しよう 北極点に十字のマークを置くと 地球儀を反時計回りに回転させると 極点上の十字のマークも反 実践海洋情報論 56 3-3 大気 海洋へ及ぶ力と流れ 3-3-1 コリオリ力 コリオリ力は 自転する惑星上の流体に働く見かけの回転力であり 惑星渦度として定義される 惑星渦度は 惑星の回転を示す角速度に依存し 極域で最大となり 赤道上では0となる 北半球において 低気圧が反時計回りに回るのはコリオリ力によるものと説明されている しかし なぜ反時計回りになるか 十分に説明のなされた図書が少ない ここでは

More information

風力発電インデックスの算出方法について 1. 風力発電インデックスについて風力発電インデックスは 気象庁 GPV(RSM) 1 局地気象モデル 2 (ANEMOS:LAWEPS-1 次領域モデル ) マスコンモデル 3 により 1km メッシュの地上高 70m における 24 時間の毎時風速を予測し

風力発電インデックスの算出方法について 1. 風力発電インデックスについて風力発電インデックスは 気象庁 GPV(RSM) 1 局地気象モデル 2 (ANEMOS:LAWEPS-1 次領域モデル ) マスコンモデル 3 により 1km メッシュの地上高 70m における 24 時間の毎時風速を予測し 風力発電インデックスの算出方法について 1. 風力発電インデックスについて風力発電インデックスは 気象庁 GPV(RSM) 1 局地気象モデル 2 (ANEMOS:LAWEPS-1 次領域モデル ) マスコンモデル 3 により 1km メッシュの地上高 70m における 24 時間の毎時風速を予測し 2000kW 定格風車の設備利用率として表示させたものです 数値は風車の定格出力 (2000kW)

More information

数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュ

数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュ 数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュレーションによって計算してみる 4.1 放物運動一様な重力場における放物運動を考える 一般に質量の物体に作用する力をとすると運動方程式は

More information

3 数値解の特性 3.1 CFL 条件 を 前の章では 波動方程式 f x= x0 = f x= x0 t f c x f =0 [1] c f 0 x= x 0 x 0 f x= x0 x 2 x 2 t [2] のように差分化して数値解を求めた ここでは このようにして得られた数値解の性質を 考

3 数値解の特性 3.1 CFL 条件 を 前の章では 波動方程式 f x= x0 = f x= x0 t f c x f =0 [1] c f 0 x= x 0 x 0 f x= x0 x 2 x 2 t [2] のように差分化して数値解を求めた ここでは このようにして得られた数値解の性質を 考 3 数値解の特性 3.1 CFL 条件 を 前の章では 波動方程式 f x= x = f x= x t f c x f = [1] c f x= x f x= x 2 2 t [2] のように差分化して数値解を求めた ここでは このようにして得られた数値解の性質を 考える まず 初期時刻 t=t に f =R f exp [ik x ] [3] のような波動を与えたとき どのように時間変化するか調べる

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 第 2 回 : 海洋数値モデルの概要 海洋大循環モデル海洋プリミティブ方程式圧力項 粘性 拡散項の取り扱い境界条件初期値 境界値問題海洋モデルの種類地球流体力学的現象地衡流平衡と温度風平衡海洋アイソスタシーとリモートセンシングどこまでダウンスケールできるか 海洋大循環モデル 地球流体力学的な現象を数値的に表現する数値モデル地球流体力学的現象 : 地球回転と成層の効果が支配的な現象 海洋大循環モデルの適用範囲

More information

ÿþŸb8bn0irt

ÿþŸb8bn0irt 折戸の物理 スペシャル補習 http://orito-buturi.com/ NO.3 今日の目的 : 1 微分方程式をもう一度 三角関数の近似について学ぶ 3 微分の意味を考える 5. 起電力 の電池, 抵抗値 の抵抗, 自己インダクタンス のコイルとスイッチを用いて右図のような回路をつくった 始めスイッチは 開かれている 時刻 t = でスイッチを閉じた 以下の問に答えよ ただし, 電流はコイルに

More information

以下 変数の上のドットは時間に関する微分を表わしている (ex. 2 dx d x x, x 2 dt dt ) 付録 E 非線形微分方程式の平衡点の安定性解析 E-1) 非線形方程式の線形近似特に言及してこなかったが これまでは線形微分方程式 ( x や x, x などがすべて 1 次で なおかつ

以下 変数の上のドットは時間に関する微分を表わしている (ex. 2 dx d x x, x 2 dt dt ) 付録 E 非線形微分方程式の平衡点の安定性解析 E-1) 非線形方程式の線形近似特に言及してこなかったが これまでは線形微分方程式 ( x や x, x などがすべて 1 次で なおかつ 以下 変数の上のドットは時間に関する微分を表わしている (e. d d, dt dt ) 付録 E 非線形微分方程式の平衡点の安定性解析 E-) 非線形方程式の線形近似特に言及してこなかったが これまでは線形微分方程式 ( や, などがすべて 次で なおかつそれらの係数が定数であるような微分方程式 ) に対して安定性の解析を行ってきた しかしながら 実際には非線形の微分方程式で記述される現象も多く存在する

More information

線積分.indd

線積分.indd 線積分 線積分 ( n, n, n ) (ξ n, η n, ζ n ) ( n-, n-, n- ) (ξ k, η k, ζ k ) ( k, k, k ) ( k-, k-, k- ) 物体に力 を作用させて位置ベクトル A の点 A から位置ベクトル の点 まで曲線 に沿って物体を移動させたときの仕事 W は 次式で計算された A, A, W : d 6 d+ d+ d@,,, d+ d+

More information

Microsoft PowerPoint - 1章 [互換モード]

Microsoft PowerPoint - 1章 [互換モード] 1. 直線運動 キーワード 速さ ( 等速直線運動, 変位 ) 加速度 ( 等加速度直線運動 ) 重力加速度 ( 自由落下 ) 力学 I 内容 1. 直線運動 2. ベクトル 3. 平面運動 4. 運動の法則 5. 摩擦力と抵抗 6. 振動 7. 仕事とエネルギー 8. 運動量と力積, 衝突 9. 角運動量 3 章以降は, 運動の向きを考えなければならない 1. 直線運動 キーワード 速さ ( 等速直線運動,

More information

流体地球科学第 7 回 力のバランス永遠に回れるバランス ( 以下, 北半球 =コリオリ力は進行方向の右向き ) 慣性振動 : 遠心力 =コリオリ力 地衡風 : コリオリ力 = 圧力傾度力 東京大学大気海洋研究所准教授藤尾伸三

流体地球科学第 7 回 力のバランス永遠に回れるバランス ( 以下, 北半球 =コリオリ力は進行方向の右向き ) 慣性振動 : 遠心力 =コリオリ力 地衡風 : コリオリ力 = 圧力傾度力 東京大学大気海洋研究所准教授藤尾伸三 流体地球科学第 7 回 力のバランス永遠に回れるバランス ( 以下, 北半球 =コリオリ力は進行方向の右向き ) 慣性振動 : 遠心力 =コリオリ力 地衡風 : コリオリ力 = 圧力傾度力 東京大学大気海洋研究所准教授藤尾伸三 http://ovd.aori.u-tokyo.ac.jp/fujio/205chiba/ [email protected] F C F A 旋衡風 : 遠心力

More information

図 7-: コリオリ力の原理 以下では 回転台の上で物体が運動したとき 物体にはたらくみかけの力を定量的に求めてみる 回転台は角速度 で回転していて 回転台に乗っている観測者から見た物体の速度ベクトルの動径方向の成分を u 接線方向の成分を v とする 図 7-3: 回転台の上での物体の運動 はじめ

図 7-: コリオリ力の原理 以下では 回転台の上で物体が運動したとき 物体にはたらくみかけの力を定量的に求めてみる 回転台は角速度 で回転していて 回転台に乗っている観測者から見た物体の速度ベクトルの動径方向の成分を u 接線方向の成分を v とする 図 7-3: 回転台の上での物体の運動 はじめ 7 大気の力学 () 7. コリオリ力 水平面内に気圧の差があると風が吹く原因となる 気圧の差によって空気塊 高にはたらく力を気圧傾度力 (pessue gaient foce) という 気圧傾度力は等 圧線と直角に 高圧側から低圧側に向かってはたらく しかし 天気図で見ら れる風向と 等圧線とのなす角は直角ではないことが多い これは 地球の自 高転の影響によって 地球上を運動する空気塊にコリオリ力

More information

物理演習問題

物理演習問題 < 物理 > =0 問 ビルの高さを, ある速さ ( 初速 をとおく,において等加速度運動の公式より (- : -= t - t : -=- t - t (-, 式よりを消去すると t - t =- t - t ( + - ( + ( - =0 0 t t t t t t ( t + t - ( t - =0 t=t t=t t - 地面 ( t - t t +t 0 より, = 3 図 問 が最高点では速度が

More information

OCW-iダランベールの原理

OCW-iダランベールの原理 講義名連続体力学配布資料 OCW- 第 2 回ダランベールの原理 無機材料工学科准教授安田公一 1 はじめに今回の講義では, まず, 前半でダランベールの原理について説明する これを用いると, 動力学の問題を静力学の問題として解くことができ, さらに, 前回の仮想仕事の原理を適用すると動力学問題も簡単に解くことができるようになる また, 後半では, ダランベールの原理の応用として ラグランジュ方程式の導出を示す

More information

DVIOUT-SS_Ma

DVIOUT-SS_Ma 第 章 微分方程式 ニュートンはリンゴが落ちるのを見て万有引力を発見した という有名な逸話があります 無重力の宇宙船の中ではリンゴは落ちないで静止していることを考えると 重力が働くと始め静止しているものが動き出して そのスピードはどんどん大きくなる つまり速度の変化が現れることがわかります 速度は一般に時間と共に変化します 速度の瞬間的変化の割合を加速度といい で定義しましょう 速度が変化する, つまり加速度がでなくなるためにはその原因があり

More information

PowerPoint Presentation

PowerPoint Presentation Non-linea factue mechanics き裂先端付近の塑性変形 塑性域 R 破壊進行領域応カ特異場 Ω R R Hutchinson, Rice and Rosengen 全ひずみ塑性理論に基づいた解析 現段階のひずみは 除荷がないとすると現段階の応力で一義的に決まる 単純引張り時の応カーひずみ関係 ( 構成方程式 ): ( ) ( ) n () y y y ここで α,n 定数, /

More information

2009 年 11 月 16 日版 ( 久家 ) 遠地 P 波の変位波形の作成 遠地 P 波の変位波形 ( 変位の時間関数 ) は 波線理論をもとに P U () t = S()* t E()* t P() t で近似的に計算できる * は畳み込み積分 (convolution) を表す ( 付録

2009 年 11 月 16 日版 ( 久家 ) 遠地 P 波の変位波形の作成 遠地 P 波の変位波形 ( 変位の時間関数 ) は 波線理論をもとに P U () t = S()* t E()* t P() t で近似的に計算できる * は畳み込み積分 (convolution) を表す ( 付録 遠地 波の変位波形の作成 遠地 波の変位波形 ( 変位の時間関数 ) は 波線理論をもとに U () t S() t E() t () t で近似的に計算できる は畳み込み積分 (convolution) を表す ( 付録 参照 ) ここで St () は地震の断層運動によって決まる時間関数 1 E() t は地下構造によって生じる種々の波の到着を与える時間関数 ( ここでは 直達 波とともに 震源そばの地表での反射波や変換波を与える時間関数

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 電磁波工学 第 5 回平面波の媒質への垂直および射入射と透過 柴田幸司 Bounda Plan Rgon ε μ Rgon Mdum ( ガラスなど ε μ z 平面波の反射と透過 垂直入射の場合 左図に示す様に 平面波が境界面に対して垂直に入射する場合を考える この時の入射波を とすると 入射波は境界において 透過波 と とに分解される この時の透過量を 反射量を Γ とおくと 領域 における媒質の誘電率に対して透過量

More information

<4D F736F F D20824F B CC92E8979D814696CA90CF95AA82C691CC90CF95AA2E646F63>

<4D F736F F D20824F B CC92E8979D814696CA90CF95AA82C691CC90CF95AA2E646F63> 1/1 平成 23 年 3 月 24 日午後 6 時 52 分 6 ガウスの定理 : 面積分と体積分 6 ガウスの定理 : 面積分と体積分 Ⅰ. 直交座標系 ガウスの定理は 微分して すぐに積分すると元に戻るというルールを 3 次元積分に適用した定理になります よく知っているのは 簡単化のため 変数が1つの場合は dj ( d ( ににします全微分 = 偏微分 d = d = J ( + C d です

More information

Microsoft Word - 8章(CI).doc

Microsoft Word - 8章(CI).doc 8 章配置間相互作用法 : Configuration Interaction () etho [] 化学的精度化学反応の精密な解析をするためには エネルギー誤差は数 ~ kcal/mol 程度に抑えたいものである この程度の誤差内に治まる精度を 化学的精度 と呼ぶことがある He 原子のエネルギーをシュレーディンガー方程式と分子軌道法で計算した結果を示そう He 原子のエネルギー Hartree-Fock

More information

ニュートン重力理論.pptx

ニュートン重力理論.pptx 3 ニュートン重力理論 1. ニュートン重力理論の基本 : 慣性系とガリレイ変換不変性 2. ニュートン重力理論の定式化 3. 等価原理 4. 流体力学方程式とその基礎 3.1 ニュートン重力理論の基本 u ニュートンの第一法則 = 力がかからなければ 等速直線運動を続ける u 等速直線運動に見える系を 慣性系 と呼ぶ ² 直線とはどんな空間の直線か? ニュートン理論では 3 次元ユークリッド空間

More information

構造力学Ⅰ第12回

構造力学Ⅰ第12回 第 回材の座屈 (0 章 ) p.5~ ( 復習 ) モールの定理 ( 手順 ) 座屈とは 荷重により梁に生じた曲げモーメントをで除して仮想荷重と考える 座屈荷重 偏心荷重 ( 曲げと軸力 ) 断面の核 この仮想荷重に対するある点でのせん断力 たわみ角に相当する曲げモーメント たわみに相当する ( 例 ) 単純梁の支点のたわみ角 : は 図 を仮想荷重と考えたときの 点の支点反力 B は 図 を仮想荷重と考えたときのB

More information

ÿþŸb8bn0irt

ÿþŸb8bn0irt 折戸の物理 スペシャル補習 http://oritobuturi.co/ NO.5(009..16) 今日の目的 : 1 物理と微分 積分について 微分方程式について学ぶ 3 近似を学ぶ 10. 以下の文を読み,[ ア ]~[ ク ] の空欄に適当な式をいれよ 物体物体に一定の大きさの力を加えたときの, 物体の運動について考え よう 右図のように, なめらかな水平面上で質量 の物体に水平に一定の大きさ

More information

自由落下と非慣性系における運動方程式 目次無重力... 2 加速度計は重力加速度を測れない... 3 重量は質量と同じ数値で kg が使える... 3 慣性系における運動方程式... 4 非慣性系における運動方程式... 6 見かけの力... 7 慣性系には実在する慣

自由落下と非慣性系における運動方程式 目次無重力... 2 加速度計は重力加速度を測れない... 3 重量は質量と同じ数値で kg が使える... 3 慣性系における運動方程式... 4 非慣性系における運動方程式... 6 見かけの力... 7 慣性系には実在する慣 自由落下と非慣性系における運動方程式 1 1 2 3 4 5 6 7 目次無重力... 2 加速度計は重力加速度を測れない... 3 重量は質量と同じ数値で kg が使える... 3 慣性系における運動方程式... 4 非慣性系における運動方程式... 6 見かけの力... 7 慣性系には実在する慣性力があるか... 7 1 2 無重力 (1) 非慣性系の住人は無重力を体感できる (a) 併進的な加速度運動をしている非慣性系の住人

More information

平面波

平面波 平面波 図.に示すように, 波源 ( 送信アンテナあるいは散乱点 ) から遠い位置で, 観測点 Pにおける波の状態を考えてみる. 遠いとは, 波長 λ に比べて距離 が十分大きいことを意味しており, 観測点 Pの近くでは, 等位相面が平面とみなせる状態にある. 平面波とは波の等位相面が平面になっている波のことである. 通信や計測を行うとき, 遠方における波の振舞いは平面波で近似できる. したがって平面波の性質を理解することが最も重要である.

More information

微分方程式による現象記述と解きかた

微分方程式による現象記述と解きかた 微分方程式による現象記述と解きかた 土木工学 : 公共諸施設 構造物の有用目的にむけた合理的な実現をはかる方法 ( 技術 ) に関する学 橋梁 トンネル ダム 道路 港湾 治水利水施設 安全化 利便化 快適化 合法則的 経済的 自然および人口素材によって作られた 質量保存則 構造物の自然的な性質 作用 ( 外力による応答 ) エネルギー則 の解明 社会的諸現象のうち マスとしての移動 流通 運動量則

More information

第1章 単 位

第1章  単  位 H. Hmno 問題解答 問題解答. 力の釣合い [ 問題.] V : sin. H :.cos. 7 V : sin sin H : cos cos cos 上第 式より これと第 式より.. cos V : sin sin H : coscos cos 上第 式より これと第 式より.98. cos [ 問題.] :. V :. : 9 9. V :. : sin V : sin 8.78 H

More information

<4D F736F F D2097CD8A7793FC96E582BD82ED82DD8A E6318FCD2E646F63>

<4D F736F F D2097CD8A7793FC96E582BD82ED82DD8A E6318FCD2E646F63> - 第 章たわみ角法の基本式 ポイント : たわみ角法の基本式を理解する たわみ角法の基本式を梁の微分方程式より求める 本章では たわみ角法の基本式を導くことにする 基本式の誘導法は各種あるが ここでは 梁の微分方程式を解いて基本式を求める方法を採用する この本で使用する座標系は 右手 右ネジの法則に従った座標を用いる また ひとつの部材では 図 - に示すように部材の左端の 点を原点とし 軸線を

More information

伝熱学課題

伝熱学課題 練習問題解答例 < 第 章強制対流熱伝達 >. 式 (.9) を導出せよ (.6) を変換する 最初に の微分値を整理しておく (.A) (.A) これを用いて の微分値を求める (.A) (.A) (.A) (.A6) (.A7) これらの微分値を式 (.6) に代入する (.A8) (.A9) (.A) (.A) (.A) (.9). 薄い平板が温度 で常圧の水の一様な流れの中に平行に置かれている

More information

剛体過去問解答例 2 1.1) 長さの棒の慣性モーメントは 公式より l I G = Ml /12 A 点のまわりは平行軸の定理より 2 2 I A = Ml /12 + M ( l / 2) = Ml 2 / 3 B y 2) 壁からの垂直抗力を R, 床からの垂直抗力と摩擦力を N,f とすると

剛体過去問解答例 2 1.1) 長さの棒の慣性モーメントは 公式より l I G = Ml /12 A 点のまわりは平行軸の定理より 2 2 I A = Ml /12 + M ( l / 2) = Ml 2 / 3 B y 2) 壁からの垂直抗力を R, 床からの垂直抗力と摩擦力を N,f とすると 剛体過去問解答例. 長さの棒の慣性モーメントは 公式より l G l A 点のまわりは平行軸の定理より A l l l B y 壁からの垂直抗力を R, 床からの垂直抗力と摩擦力を N,f とすると 運動方程式は 方向 : R f, y 方向 : y N l 回転 : G { N f R cos } A 静止しているとき 方向の力と 力のモーメントがつり合うので y ~ より R ' また 摩擦力が最大静止摩擦力より大きいとはしごは動き出すので

More information

Laplace2.rtf

Laplace2.rtf =0 ラプラスの方程式は 階の微分方程式で, 一般的に3つの座標変数をもつ. ここでは, 直角座標系, 円筒座標系, 球座標系におけるラプラスの方程式の解き方を説明しよう. 座標変数ごとに方程式を分離し, それを解いていく方法は変数分離法と呼ばれる. 変数分離解と固有関数展開法. 直角座標系における 3 次元の偏微分方程式 = x + y + z =0 (.) を解くために,x, y, z について互いに独立な関数の積で成り立っていると考え,

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 3 か月予報 (11 月 ~1 月の天候 の見通し ) とその解説 気象庁地球環境 海洋部 気候情報課 1 季節予報が対象とする大気の変動 空間スケ ル km 10 4 10 3 10 2 10 1 積乱雲 熱帯季節内変動テレコネクション定常ロスビー波ブロッキング総観規模高 低気圧 メソスケール低気圧 エルニーニョ現象 アジアモンスーンの変動 海洋の影響を強く受けた変動 十年規模変動 温暖化 10

More information

<4D F736F F D2094F795AA95FB92F68EAE82CC89F082AB95FB E646F63>

<4D F736F F D2094F795AA95FB92F68EAE82CC89F082AB95FB E646F63> 力学 A 金曜 限 : 松田 微分方程式の解き方 微分方程式の解き方のところが分からなかったという声が多いので プリントにまとめます 数学的に厳密な話はしていないので 詳しくは数学の常微分方程式を扱っているテキストを参照してください また os s は既知とします. 微分方程式の分類 常微分方程式とは 独立変数 と その関数 その有限次の導関数 がみたす方程式 F,,, = のことです 次までの導関数を含む方程式を

More information

相対性理論入門 1 Lorentz 変換 光がどのような座標系に対しても同一の速さ c で進むことから導かれる座標の一次変換である. (x, y, z, t ) の座標系が (x, y, z, t) の座標系に対して x 軸方向に w の速度で進んでいる場合, 座標系が一次変換で関係づけられるとする

相対性理論入門 1 Lorentz 変換 光がどのような座標系に対しても同一の速さ c で進むことから導かれる座標の一次変換である. (x, y, z, t ) の座標系が (x, y, z, t) の座標系に対して x 軸方向に w の速度で進んでいる場合, 座標系が一次変換で関係づけられるとする 相対性理論入門 Lorentz 変換 光がどのような座標系に対しても同一の速さ で進むことから導かれる座標の一次変換である. x, y, z, t ) の座標系が x, y, z, t) の座標系に対して x 軸方向に w の速度で進んでいる場合, 座標系が一次変換で関係づけられるとすると, x A x wt) y y z z t Bx + Dt 弨弱弩弨弲弩弨弳弩弨弴弩 が成立する. 図 : 相対速度

More information

Microsoft PowerPoint - 流体力学の基礎02(OpenFOAM 勉強会 for geginner).pptx

Microsoft PowerPoint - 流体力学の基礎02(OpenFOAM 勉強会 for geginner).pptx ~ 流体力学の基礎 ~ 第 2 回 流体静力学 2011 年 10 月 22 日 ( 土 ) 講習会のスケジュール概要 ( あくまでも現時点での予定です ) 流体力学の基礎 第 1 回目 2011.09 流体について 第 2 回目 2011.10 流体静力学 第 3 回目 2011.11/12 流体運動の基礎理論 1 第 4 回目 2012.01 流体運動の基礎理論 2 第 5 回目 2012.02

More information

はじめに 動機 IGModel プロジェクトとは 目次 IGModel プロジェクトの構成 構造 現在の開発状況 IGMBaseLib, IGMTool IGModel-SW( 正二十面体格子全球浅水モデル ) IGModel-SW のテスト計算の結果 ( 一部 ) まとめ, 今後の展望

はじめに 動機 IGModel プロジェクトとは 目次 IGModel プロジェクトの構成 構造 現在の開発状況 IGMBaseLib, IGMTool IGModel-SW( 正二十面体格子全球浅水モデル ) IGModel-SW のテスト計算の結果 ( 一部 ) まとめ, 今後の展望 正二十面体格子大気モデル IGModel プロジェクトの紹介 2011/08/21 GFD セミナー 2011 河合佑太神戸大学地球および惑星大気科学研究室 はじめに 動機 IGModel プロジェクトとは 目次 IGModel プロジェクトの構成 構造 現在の開発状況 IGMBaseLib, IGMTool IGModel-SW( 正二十面体格子全球浅水モデル ) IGModel-SW のテスト計算の結果

More information

ギリシャ文字の読み方を教えてください

ギリシャ文字の読み方を教えてください 埼玉工業大学機械工学学習支援セミナー ( 小西克享 ) 単振り子の振動の近似解と厳密解 -/ テーマ H: 単振り子の振動の近似解と厳密解. 運動方程式図 のように, 質量 m のおもりが糸で吊り下げられている時, おもりには重力 W と糸の張力 が作用しています. おもりは静止した状態なので,W と F は釣り合った状態注 ) になっています. すなわち, W です.W は質量 m と重力加速度

More information

7 章問題解答 7-1 予習 1. 長方形断面であるため, 断面積 A と潤辺 S は, 水深 h, 水路幅 B を用い以下で表される A = Bh, S = B + 2h 径深 R の算定式に代入すると以下のようになる A Bh h R = = = S B + 2 h 1+ 2( h B) 分母の

7 章問題解答 7-1 予習 1. 長方形断面であるため, 断面積 A と潤辺 S は, 水深 h, 水路幅 B を用い以下で表される A = Bh, S = B + 2h 径深 R の算定式に代入すると以下のようになる A Bh h R = = = S B + 2 h 1+ 2( h B) 分母の 7 章問題解答 7- 予習. 長方形断面であるため, 断面積 と潤辺 S は, 水深, 水路幅 B を用い以下で表される B, S B + 径深 R の算定式に代入すると以下のようになる B R S B + ( B) 分母の /B は河幅が水深に対して十分に広ければ, 非常に小さな値となるため, 上式は R ( B) となり, 径深 R は水深 で近似できる. マニングの式の水深 を等流水深 0 と置き換えると,

More information

Microsoft PowerPoint - H21生物計算化学2.ppt

Microsoft PowerPoint - H21生物計算化学2.ppt 演算子の行列表現 > L いま 次元ベクトル空間の基底をケットと書くことにする この基底は完全系を成すとすると 空間内の任意のケットベクトルは > > > これより 一度基底を与えてしまえば 任意のベクトルはその基底についての成分で完全に記述することができる これらの成分を列行列の形に書くと M これをベクトル の基底 { >} による行列表現という ところで 行列 A の共役 dont 行列は A

More information

FEM原理講座 (サンプルテキスト)

FEM原理講座 (サンプルテキスト) サンプルテキスト FEM 原理講座 サイバネットシステム株式会社 8 年 月 9 日作成 サンプルテキストについて 各講師が 講義の内容が伝わりやすいページ を選びました テキストのページは必ずしも連続していません 一部を抜粋しています 幾何光学講座については 実物のテキストではなくガイダンスを掲載いたします 対象とする構造系 物理モデル 連続体 固体 弾性体 / 弾塑性体 / 粘弾性体 / 固体

More information

Microsoft Word - thesis.doc

Microsoft Word - thesis.doc 剛体の基礎理論 -. 剛体の基礎理論初めに本論文で大域的に使用する記号を定義する. 使用する記号トルク撃力力角運動量角速度姿勢対角化された慣性テンソル慣性テンソル運動量速度位置質量時間 J W f F P p .. 質点の並進運動 質点は位置 と速度 P を用いる. ニュートンの運動方程式 という状態を持つ. 但し ここでは速度ではなく運動量 F P F.... より質点の運動は既に明らかであり 質点の状態ベクトル

More information

例 e 指数関数的に減衰する信号を h( a < + a a すると, それらのラプラス変換は, H ( ) { e } e インパルス応答が h( a < ( ただし a >, U( ) { } となるシステムにステップ信号 ( y( のラプラス変換 Y () は, Y ( ) H ( ) X (

例 e 指数関数的に減衰する信号を h( a < + a a すると, それらのラプラス変換は, H ( ) { e } e インパルス応答が h( a < ( ただし a >, U( ) { } となるシステムにステップ信号 ( y( のラプラス変換 Y () は, Y ( ) H ( ) X ( 第 週ラプラス変換 教科書 p.34~ 目標ラプラス変換の定義と意味を理解する フーリエ変換や Z 変換と並ぶ 信号解析やシステム設計における重要なツール ラプラス変換は波動現象や電気回路など様々な分野で 微分方程式を解くために利用されてきた ラプラス変換を用いることで微分方程式は代数方程式に変換される また 工学上使われる主要な関数のラプラス変換は簡単な形の関数で表されるので これを ラプラス変換表

More information

大阪大学物理 8 を解いてみた Ⅱ. 問 ( g cosq a sin q ) m - 台 B 上の観測者から見ると, 小物体は, 斜面からの垂直抗力 N, 小物体の重力 mg, 水平左向きの慣性力 ma を受け, 台 B の斜面と平行な向きに運動する したがって, 小物体は台 B の斜面に垂直な方

大阪大学物理 8 を解いてみた Ⅱ. 問 ( g cosq a sin q ) m - 台 B 上の観測者から見ると, 小物体は, 斜面からの垂直抗力 N, 小物体の重力 mg, 水平左向きの慣性力 ma を受け, 台 B の斜面と平行な向きに運動する したがって, 小物体は台 B の斜面に垂直な方 大阪大学物理 8 を解いてみた Ⅰ. 問 g 最高点の座標を y max とすると, 力学的エネルギー保存則より \ y m mgy 補足 max g max 小物体の運動方向に対する仕事は重力 ( 保存力 ) の斜面に沿った成分のみであり, 垂直抗力 ( 非保存力 ) の仕事は である よって, 力学的エネルギー保存則が成り立つ これを確かめてみよう 小物体は重力の斜面に沿った外力を受けながらその運動エネルギーを失っていく

More information

() 実験 Ⅱ. 太陽の寿命を計算する 秒あたりに太陽が放出している全エネルギー量を計測データをもとに求める 太陽の放出エネルギーの起源は, 水素の原子核 4 個が核融合しヘリウムになるときのエネルギーと仮定し, 質量とエネルギーの等価性から 回の核融合で放出される全放射エネルギーを求める 3.から

() 実験 Ⅱ. 太陽の寿命を計算する 秒あたりに太陽が放出している全エネルギー量を計測データをもとに求める 太陽の放出エネルギーの起源は, 水素の原子核 4 個が核融合しヘリウムになるときのエネルギーと仮定し, 質量とエネルギーの等価性から 回の核融合で放出される全放射エネルギーを求める 3.から 55 要旨 水温上昇から太陽の寿命を算出する 53 町野友哉 636 山口裕也 私たちは, 地球環境に大きな影響を与えている太陽がいつまで今のままであり続けるのかと疑問をもちました そこで私たちは太陽の寿命を求めました 太陽がどのように燃えているのかを調べたら水素原子がヘリウム原子に変化する核融合反応によってエネルギーが発生していることが分かった そこで, この反応が終わるのを寿命と考えて算出した

More information

Taro-解答例NO3放物運動H16

Taro-解答例NO3放物運動H16 放物運動 解答のポイント 初速度, 水平との角度 θ で 高さ の所から投げあげるとき 秒後の速度 =θ =θ - 秒後の位置 =θ 3 ( 水平飛行距離 ) =θ - + 4 ( 高さ ) ~4 の導出は 基本問題 参照 ( 地上から投げた場合の図 : 教科書参照 ) 最高点の 高さ 最高点では において = 水平到達距離 より 最高点に到達する時刻 を求め 4に代入すると最高点の高さH 地上では

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 平成 28 年 1 月 26 日 エルニーニョ現象と 世界 日本の天候 安田珠幾 エルニーニョ情報管理官気象庁地球環境 海洋部気候情報課 はじめに 1 はじめに 現在 1997-98 年のエルニーニョ現象以来の強いエルニーニョ現象が発生中 エルニーニョ現象は世界の異常気象を引き起こし 日本には 冷夏 暖冬 をもたらすと言われる エルニーニョ現象はなぜ世界の広い範囲の天候に影響を及ぼすのか? そもそもエルニーニョ現象とは?

More information

Microsoft PowerPoint - 夏の学校(CFD).pptx

Microsoft PowerPoint - 夏の学校(CFD).pptx /9/5 FD( 計算流体力学 ) の基礎理論 性能 運動分野 夏の学校 神戸大学大学院海事科学研究科勝井辰博 流体の質量保存 流体要素内の質量の増加率 [ 単位時間当たりの増加量 ] 単位時間に流体要素に流入する質量 流体要素 Fl lm (orol olm) v ( ) ガウスの定理 v( ) /9/5 = =( ) b=b =(b b b ) b= b = b + b + b アインシュタイン表記

More information

応用数学Ⅱ 偏微分方程式(2) 波動方程式(12/13)

応用数学Ⅱ 偏微分方程式(2) 波動方程式(12/13) 偏微分方程式. 偏微分方程式の形 偏微分 偏導関数 つの独立変数 をもつ関数 があるとき 変数 が一定値をとって だけが変化したとす ると は だけの関数となる このとき を について微分して得られる関数を 関数 の に関する 偏微分係数 略して偏微分 あるいは偏導関数 pil deiie といい 次のように表される についても同様な偏微分を定義できる あるいは あるいは - あるいは あるいは -

More information

木村の物理小ネタ 単振動と単振動の力学的エネルギー 1. 弾性力と単振動 弾性力も単振動も力は F = -Kx の形で表されるが, x = 0 の位置は, 弾性力の場合, 弾性体の自然状態の位置 単振動の場合, 振動する物体に働く力のつり合

木村の物理小ネタ   単振動と単振動の力学的エネルギー 1. 弾性力と単振動 弾性力も単振動も力は F = -Kx の形で表されるが, x = 0 の位置は, 弾性力の場合, 弾性体の自然状態の位置 単振動の場合, 振動する物体に働く力のつり合 単振動と単振動の力学的エネルギー. 弾性力と単振動 弾性力も単振動も力は F = -x の形で表されるが, x = の位置は, 弾性力の場合, 弾性体の自然状態の位置 単振動の場合, 振動する物体に働く力のつり合いの位置 である たとえば, おもりをつるしたばねについて, ばねの弾性力を考えるときは, ばねの自然長を x = とし, おもりの単振動で考える場合は, おもりに働く力がつり合った位置を

More information

Microsoft PowerPoint - 第7章(自然対流熱伝達 )_H27.ppt [互換モード]

Microsoft PowerPoint - 第7章(自然対流熱伝達 )_H27.ppt [互換モード] 第 7 章自然対流熱伝達 伝熱工学の基礎 : 伝熱の基本要素 フーリエの法則 ニュートンの冷却則 次元定常熱伝導 : 熱伝導率 熱通過率 熱伝導方程式 次元定常熱伝導 : ラプラスの方程式 数値解析の基礎 非定常熱伝導 : 非定常熱伝導方程式 ラプラス変換 フーリエ数とビオ数 対流熱伝達の基礎 : 熱伝達率 速度境界層と温度境界層 層流境界層と乱流境界層 境界層厚さ 混合平均温度 強制対流熱伝達 :

More information

位相最適化?

位相最適化? 均質化設計法 藤井大地 ( 東京大学 ) 位相最適化? 従来の考え方 境界形状を変化させて最適な形状 位相を求める Γ t Ω b Γ D 境界形状を変化させる問題点 解析が進むにつれて, 有限要素メッシュが異形になり, 再メッシュが必要になる 位相が変化する問題への適応が難しい Γ Γ t t Ω b Ω b Γ D Γ D 領域の拡張と特性関数の導入 χ Ω ( x) = f 0 f x Ω x

More information

破壊の予測

破壊の予測 本日の講義内容 前提 : 微分積分 線形代数が何をしているかはうろ覚え 材料力学は勉強したけど ちょっと 弾性および塑性学は勉強したことが無い ー > ですので 解らないときは質問してください モールの応力円を理解するとともに 応力を 3 次元的に考える FM( 有限要素法 の概略 内部では何を計算しているのか? 3 物が壊れる条件を考える 特に 変形 ( 塑性変形 が発生する条件としてのミーゼス応力とはどのような応力か?

More information

今週の内容 後半全体のおさらい ラグランジュの運動方程式の導出 リンク機構のラグランジュの運動方程式 慣性行列 リンク機構のエネルギー保存則 エネルギー パワー 速度 力の関係 外力が作用する場合の運動方程式 粘性 粘性によるエネルギーの消散 慣性 粘性 剛性と微分方程式 拘束条件 ラグランジュの未

今週の内容 後半全体のおさらい ラグランジュの運動方程式の導出 リンク機構のラグランジュの運動方程式 慣性行列 リンク機構のエネルギー保存則 エネルギー パワー 速度 力の関係 外力が作用する場合の運動方程式 粘性 粘性によるエネルギーの消散 慣性 粘性 剛性と微分方程式 拘束条件 ラグランジュの未 力学 III GA 工業力学演習 X5 解析力学 5X 5 週目 立命館大学機械システム系 8 年度後期 今週の内容 後半全体のおさらい ラグランジュの運動方程式の導出 リンク機構のラグランジュの運動方程式 慣性行列 リンク機構のエネルギー保存則 エネルギー パワー 速度 力の関係 外力が作用する場合の運動方程式 粘性 粘性によるエネルギーの消散 慣性 粘性 剛性と微分方程式 拘束条件 ラグランジュの未定乗数法

More information

第1章 単 位

第1章  単  位 H. Hamano,. 長柱の座屈 - 長柱の座屈 長い柱は圧縮荷重によって折れてしまう場合がある. この現象を座屈といい, 座屈するときの荷重を座屈荷重という.. 換算長 長さ の柱に荷重が作用する場合, その支持方法によって, 柱の理論上の長さ L が異なる. 長柱の計算は, この L を用いて行うと都合がよい. この L を換算長 ( あるいは有効長さという ) という. 座屈荷重は一般に,

More information

国土技術政策総合研究所 研究資料

国土技術政策総合研究所 研究資料 参考資料 崩壊の恐れのある土層厚の空間分布を考慮したがけ崩れ対策に関する検討 参考資料 崩壊の恐れのある土層厚の空間分布を考慮したがけ崩れ対策に関する検討 ここでは 5 章で示した方法により急傾斜地における崩壊する恐れがある層厚の面的分布が明らかとなった場合のがけ崩れ対策手法について検討する 崩壊する恐れがある層厚の面的な分布は 1 土砂災害警戒区域等における土砂災害防止対策の推進に関する法律( 以下

More information

Microsoft PowerPoint - 海の波.ppt [互換モード]

Microsoft PowerPoint - 海の波.ppt [互換モード] 008/5/6 風波の発生 ( 深海域 波浪について ( 初級編 波の浅水変形と砕波 波による質量輸送と海浜流系の発達 離岸流の発生 ( 浅海域 3.6m の高波が襲来 水粒子の軌道 ( スト - クスドリフト 離岸流 沿岸流 (1 波 浅海域の波浪変形の研究 1. 海岸保全に関して (A 季節風時 台風時の波と流れによる地形変化 ( 侵食 (B 漁港 港湾 海岸保全構造物による波 流れ 地形変化

More information

人間科学部研究年報平成 24 年 (1) (2) (3) (4) 式 (1) は, クーロン (Coulomb) の法則とも呼ばれる.ρは電荷密度を表し,ε 0 は真空の誘電率と呼ばれる定数である. 式 (2) は, 磁荷が存在しないことを表す式である. 式 (3) はファラデー (Faraday)

人間科学部研究年報平成 24 年 (1) (2) (3) (4) 式 (1) は, クーロン (Coulomb) の法則とも呼ばれる.ρは電荷密度を表し,ε 0 は真空の誘電率と呼ばれる定数である. 式 (2) は, 磁荷が存在しないことを表す式である. 式 (3) はファラデー (Faraday) 複素振幅をもつ球面波の人間科学部研究年報 Maxwell 平成 24 方程式年 複素振幅をもつ球面波の Maxwell 方程式 Maxwell Equation of Spherical Wave with Complex Amplitude 戸上良弘 Yoshihiro TOGAMI Abstract 複素振幅をもつ球面波に関して, マクスウェル (Maxwell) 方程式との関係を考察した. 電気的な球面波としてのスカラーポテンシャルが与えられたとき,

More information

τ-→K-π-π+ν τ崩壊における CP対称性の破れの探索

τ-→K-π-π+ν τ崩壊における CP対称性の破れの探索 τ - K - π - π + ν τ 崩壊における CP 対称性の破れの探索 奈良女子大学大学院人間文化研究科 物理科学専攻高エネルギー物理学研究室 近藤麻由 1 目次 はじめに - τ 粒子の概要 - τ - K - π - π + ν τ 崩壊における CP 対称性の破れ 実験装置 事象選別 τ - K - π - π + ν τ 崩壊の不変質量分布 CP 非対称度の解析 - モンテカルロシミュレーションによるテスト

More information

木村の物理小ネタ ケプラーの第 2 法則と角運動量保存則 A. 面積速度面積速度とは平面内に定点 O と動点 P があるとき, 定点 O と動点 P を結ぶ線分 OP( 動径 OP という) が単位時間に描く面積を 動点 P の定点 O に

木村の物理小ネタ   ケプラーの第 2 法則と角運動量保存則 A. 面積速度面積速度とは平面内に定点 O と動点 P があるとき, 定点 O と動点 P を結ぶ線分 OP( 動径 OP という) が単位時間に描く面積を 動点 P の定点 O に ケプラーの第 法則と角運動量保存則 A. 面積速度面積速度とは平面内に定点 O と動点 P があるとき, 定点 O と動点 P を結ぶ線分 OP( 動径 OP という が単位時間に描く面積を 動点 P の定点 O に関する面積速度の大きさ という 定点 O まわりを回る面積速度の導き方導き方 A ( x( + D, y( + D v ( q r ( A ( x (, y( 動点 P が xy 座標平面上を時刻

More information

2011年度 大阪大・理系数学

2011年度 大阪大・理系数学 0 大阪大学 ( 理系 ) 前期日程問題 解答解説のページへ a a を自然数とする O を原点とする座標平面上で行列 A= a の表す 次変換 を f とする cosθ siθ () >0 および0θ

More information

高校電磁気学 ~ 電磁誘導編 ~ 問題演習

高校電磁気学 ~ 電磁誘導編 ~ 問題演習 高校電磁気学 ~ 電磁誘導編 ~ 問題演習 問 1 磁場中を動く導体棒に関する問題 滑車 導体棒の間隔 L m a θ (1) おもりの落下速度が のとき 導体棒 a に生じる誘導起電力の 大きさを求めよ 滑車 導体棒の間隔 L m a θ 導体棒の速度 水平方向の速度 cosθ Δt の時間に回路を貫く磁束の変化 ΔΦ は ΔΦ = ΔS = LcosθΔt ΔΦ ファラデーの法則 V = N より

More information

0 21 カラー反射率 slope aspect 図 2.9: 復元結果例 2.4 画像生成技術としての計算フォトグラフィ 3 次元情報を復元することにより, 画像生成 ( レンダリング ) に応用することが可能である. 近年, コンピュータにより, カメラで直接得られない画像を生成する技術分野が生

0 21 カラー反射率 slope aspect 図 2.9: 復元結果例 2.4 画像生成技術としての計算フォトグラフィ 3 次元情報を復元することにより, 画像生成 ( レンダリング ) に応用することが可能である. 近年, コンピュータにより, カメラで直接得られない画像を生成する技術分野が生 0 21 カラー反射率 slope aspect 図 2.9: 復元結果例 2.4 画像生成技術としての計算フォトグラフィ 3 次元情報を復元することにより, 画像生成 ( レンダリング ) に応用することが可能である. 近年, コンピュータにより, カメラで直接得られない画像を生成する技術分野が生まれ, コンピューテーショナルフォトグラフィ ( 計算フォトグラフィ ) と呼ばれている.3 次元画像認識技術の計算フォトグラフィへの応用として,

More information

<4D F736F F F696E74202D C CC89C88A B8CDD8AB B83685D>

<4D F736F F F696E74202D C CC89C88A B8CDD8AB B83685D> 断面積 (A) を使わずに, 間隙率を使う透水係数の算定 図に示したような 本の孔を掘って, 上流側から食塩を投入した 食塩を投入してから,7 時間後に下流側に食塩が到達したことが分かった この地盤の透水係数を求めよ 地盤の間隙比は e=0.77, 水位差は 0 cmであった なお, この方法はトレーサ法の中の食塩法と呼ばれている Nacl 計測器 0 cm 0.0 m 断面積 (A) を使わずに,

More information

<4D F736F F D20824F F6490CF95AA82C696CA90CF95AA2E646F63>

<4D F736F F D20824F F6490CF95AA82C696CA90CF95AA2E646F63> 1/15 平成 3 年 3 月 4 日午後 6 時 49 分 5 ベクトルの 重積分と面積分 5 重積分と面積分 Ⅰ. 重積分 と で 回積分することを 重積分 といいます この 重積分は何を意味しているのでしょう? 通常の積分 (1 重積分 ) では C d 図 1a 1 f d (5.1) 1 f d f ( ) は 図形的には図 1a のように面積を表しています つまり 1 f ( ) を高さとしてプロットすると図

More information

<4D F736F F F696E74202D D488A778AEE B4F93B982CC8AEE A2E707074>

<4D F736F F F696E74202D D488A778AEE B4F93B982CC8AEE A2E707074> 宇宙工学基礎 ( 軌道の基礎 松永三郎 機械宇宙学科 機械宇宙システム専攻 ニュートンの法則 第 法則 力が作用作用しないしない限り 質点質点は静止静止ないしはないしは一定速度一定速度で運動するする ( 慣性の法則 慣性空間 慣性座標系慣性座標系の定義第 法則 慣性座標系におけるにおける質点質点の運動 p F ( pɺ t ( F: 全作用力, pmv: 並進運動量 ( 質量と速度速度の積 慣性系を規準規準としてとして時間微分時間微分を行うことにことに注意第

More information

s ss s ss = ε = = s ss s (3) と表される s の要素における s s = κ = κ, =,, (4) jωε jω s は複素比誘電率に相当する物理量であり ここで PML 媒質定数を次のように定義する すなわち κξ をPML 媒質の等価比誘電率 ξ をPML 媒質の

s ss s ss = ε = = s ss s (3) と表される s の要素における s s = κ = κ, =,, (4) jωε jω s は複素比誘電率に相当する物理量であり ここで PML 媒質定数を次のように定義する すなわち κξ をPML 媒質の等価比誘電率 ξ をPML 媒質の FDTD 解析法 (Matlab 版 2 次元 PML) プログラム解説 v2.11 1. 概要 FDTD 解析における吸収境界である完全整合層 (Perfectl Matched Laer, PML) の定式化とプログラミングを2 次元 TE 波について解説する PMLは異方性の損失をもつ仮想的な物質であり 侵入して来る電磁波を逃さず吸収する 通常の物質と接する界面でインピーダンスが整合しており

More information

反射係数

反射係数 平面波の反射と透過 電磁波の性質として, 反射と透過は最も基礎的な現象である. 我々の生活している空間は, 各種の形状を持った媒質で構成されている. 人間から見れば, 空気, 水, 木, 土, 火, 金属, プラスチックなど, 全く異なるものに見えるが, 電磁波からすると誘電率, 透磁率, 導電率が異なるだけである. 磁性体を除く媒質は比透磁率がで, ほとんど媒質に当てはまるので, 実質的に我々の身の回りの媒質で,

More information

CERT化学2013前期_問題

CERT化学2013前期_問題 [1] から [6] のうち 5 問を選んで解答用紙に解答せよ. いずれも 20 点の配点である.5 問を超えて解答した場合, 正答していれば成績評価に加算する. 有効数字を適切に処理せよ. 断りのない限り大気圧は 1013 hpa とする. 0 C = 273 K,1 cal = 4.184 J,1 atm = 1013 hpa = 760 mmhg, 重力加速度は 9.806 m s 2, 気体

More information

多次元レーザー分光で探る凝縮分子系の超高速動力学

多次元レーザー分光で探る凝縮分子系の超高速動力学 波動方程式と量子力学 谷村吉隆 京都大学理学研究科化学専攻 http:theochem.kuchem.kyoto-u.ac.jp TA: 岩元佑樹 [email protected] ベクトルと行列の作法 A 列ベクトル c = c c 行ベクトル A = [ c c c ] 転置ベクトル T A = [ c c c ] AA 内積 c AA = [ c c c ] c =

More information

Microsoft PowerPoint - 第5回電磁気学I 

Microsoft PowerPoint - 第5回電磁気学I  1 年 11 月 8 日 ( 月 ) 1:-1: Y 平成 年度工 系 ( 社会環境工学科 ) 第 5 回電磁気学 Ⅰ 天野浩 項目 電界と電束密度 ガウスの発散定理とガウスの法則の積分形と微分形 * ファラデーの電気力線の使い方をマスターします * 電界と電束密度を定義します * ガウスの発散定理を用いて ガウスの法則の積分形から微分形をガウスの法則の積分形から微分形を導出します * ガウスの法則を用いて

More information

ギリシャ文字の読み方を教えてください

ギリシャ文字の読み方を教えてください 埼玉工業大学機械工学学習支援セミナー ( 小西克享 ) 慣性モーメント -1/6 テーマ 01: 慣性モーメント (Momet of ietia) コマ回しをすると, 長い時間回転させるには重くて大きなコマを選ぶことや, ひもを早く引くことが重要であることが経験的にわかります. 遊びを通して, 回転の運動エネルギーを増やせば, 回転の勢いが増すことを学習できるので, 機械系の学生にとってコマ回しも大切な体験学習のひとつと言えます.

More information