PowerPoint プレゼンテーション

Size: px
Start display at page:

Download "PowerPoint プレゼンテーション"

Transcription

1 原始惑星系円盤内でロスビー波不安定性によって形成される渦 小野智弘 ( 京都大 ), 武藤恭之 ( 工学院大 ), 富田賢吾 ( 大阪大 ), 野村英子 ( 東工大 ) Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 1

2 様々な原始惑星系円盤構造 若い星の周りにあるガス円盤 円盤内のダストが合体成長し 惑星を形成 近年 詳細な円盤構造が明らかになってきている ALMA によるダスト連続光観測 非軸対称 リング スパイラル HD HL Tau Elias 2-27 Fukagawa+13 ALMA Partnership+15 Pérez+16 Dec. 20th, 2016 理論懇シンポジウム 東北大 2

3 Fukagawa+13 Fukagawa+13 ガスの非軸対称構造 HD ダスト連続光観測 Muto+15 実線 : ガス分布破線 : ダスト分布赤線 : 高面密度側青線 : 低面密度側 ガス輝線観測 理論モデル ダスト分布は強く非軸対称 ( 約 70 倍 ) ガス分布は弱く非軸対称 (3-10 倍 ) Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 3

4 どうやって非軸対称構造を作るのか ガス渦 を用いるモデルが有力 1. 流体不安定性によるガス渦形成 ロスビー波不安定性 (e.g., Lovelace+99, Li+00, 01) Baroclinic 不安定性 (e.g., Klahr & Bodenheimer 03) Vertical shear 不安定性 (e.g., Goldreich & Schubert 67) Zombie vortex 不安定性 (Marcus+13, 15) 2. ガス渦によるダスト捕獲 r ガス ダスト Fu+14 Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 4

5 面密度 ロスビー波不安定性 [RWI] 差動回転円盤中で起こるシアー不安定性 動径方向の密度分布に強い変化がある時 不安定 ex.) 密度バンプ 軸対称密度バンプ 20 半径 不安定時 ガス渦を形成 複数の渦ができた後 合体により一つの渦に 最終的に渦移動が起こり中心星に落下 Dec. 20th, 2016 理論懇シンポジウム 東北大 5

6 面密度 ロスビー波不安定性 [RWI] 差動回転円盤中で起こるシアー不安定性 動径方向の密度分布に強い変化がある時 不安定 ex.) 密度バンプ 軸対称密度バンプ 20 半径 不安定時 ガス渦を形成 複数の渦ができた後 合体により一つの渦に 最終的に渦移動が起こり中心星に落下 Dec. 20th, 2016 理論懇シンポジウム 東北大 6

7 面密度 ロスビー波不安定性 [RWI] 差動回転円盤中で起こるシアー不安定性 動径方向の密度分布に強い変化がある時 不安定 ex.) 密度バンプ 軸対称密度バンプ 800 半径 不安定時 ガス渦を形成 複数の渦ができた後 合体により一つの渦に 最終的に渦移動が起こり中心星に落下 Dec. 20th, 2016 理論懇シンポジウム 東北大 7

8 面密度 ロスビー波不安定性 [RWI] 差動回転円盤中で起こるシアー不安定性 動径方向の密度分布に強い変化がある時 不安定 ex.) 密度バンプ 軸対称密度バンプ 800 半径 不安定時 ガス渦を形成 複数の渦ができた後 合体により一つの渦に 最終的に渦移動が起こり中心星に落下 Dec. 20th, 2016 理論懇シンポジウム 東北大 8

9 不安定成長率 面密度 線形安定性解析 (Ono+16) 背景流 2 次元 平衡 軸対称 平行流 ガウシアン面密度分布高さ A, 幅 w がパラメータ A 1 w r 仮定 一様エントロピー ( 順圧 ) 高い 不安定成長率の分布レイリー不安定 結果 ( 円盤スケールハイト =0.1) A 非軸対称モードに対するバンプの不安定性を調べた 中立安定線を求めた 低い #Poster 9 Muto: 非圧縮の場合 小さい w 大きい Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 9

10 不安定成長率 レイリー不安定 レイリー条件 比角運動量 (j = r 2 Ω) の動径勾配が負になる場所が存在する レイリー条件を満たすと不安定となって動径勾配が負となる領域を消す 面密度の動径変化が強い時 圧力勾配力の効果でレイリー条件を満たし得る 高い A 不安定成長率の分布 レイリー不安定 低い小さい w 大きい Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 10

11 本研究の目的 目標 観測される円盤非軸対称構造の成因がロスビー波不安定性であるかを明らかにする 本研究の目的 1. 線形安定性解析の結果を検証する 2. ロスビー波不安定性によってどのような渦が形成されるのかを調べる 使う道具 2 次元流体シミュレーション (athena++) Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 11

12 面密度 数値計算の設定 Code : Athena++ Coordinates : cylindrical (2D: r φ) Reimann Solver : HLLC Spatial Reconstruction : PLM(2 次 ) Hydro time-integration : VL2(2 次 ) 仮定 : 一様エントロピー流 ( 順圧 ) 計算領域 (r) : [0.3:3.1] dr: 渦領域 0.002, 領域外 計算領域 (φ) : [-π:π] grid: 渦領域 1024, 領域外 512 計算時間 : r=1 の場所で200 回転 (t=400π) 減衰領域 : [0.3,0.375], [2.6:3.1] で速度を初期値へ減衰境界条件 (r) : 初期条件で固定 境界条件 (φ): 周期境界計算機 : 基礎物理研究所 CRAY XC40, CfCA Cray XC30 初期条件 2 次元 平衡 軸対称 平行流ガウシアン面密度分布 # パラメータ ( 高さ A, 幅 w) A 1 ersity.github.io/athe na/index.html w r Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 12

13 解析方法面密度分布 方位角平均面密度をガウシアンでフィッティング Σ=1+A exp r r c Δw 高い 2 phase2 1 つの渦が移動 方位角平均面密度分布 A 低い phase1 数 ~ 数十回転で渦が合体 小さい w 大きい Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 13

14 結果 1: 進化線 円盤スケールハイト =0.1 レイリー中立安定の状況を初期条件に 10run A 平均面密度分布の進化線 phase1 phase2 高い レイリー不安定 phase 転換点 全 run で渦形成 w A 初期バンプが大きいほど phase2 が見えなくなる 中立安定線で進化は止まるわけではない 低い 中立安定線 小さい w : 初期条件 : phase 転換点 大きい Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 14

15 結果 2: 規格化した進化線 case A の (A, w) を初期値で規格化した場合の進化線 Phase 転換点は A/A 0 = ( w/ w 0 ) 1.25 の線上にのる 高い 平均面密度分布の規格化進化線 初期バンプが小さいほど 大きい割合で進化する A/A 0 低い 緑線 : A ( w) 1 青線 : A ( w) 1.25 小さい w/ w 0 大きい Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 15

16 結果 3: 進化線上を初期条件におく A ( w) 1.25 線上を初期条件として計算 レイリー中立安定 ~ ロスビー中立安定 : 4 run ロスビー中立安定の時以外で渦形成 全 run で phase 転換点は A/A 0 = ( w/ w 0 ) 1.25 の線上にのる 初期バンプが大きい時ほど phase2 が見えなくなる 高い A 平均面密度分布の進化線 レイリー不安定 低い 小さい w 大きい Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 16

17 結果まとめ 1. 線形安定性解析の検証 不安定領域では渦形成が起こる 中立安定 安定領域では渦形成は起こらない 2. ロスビー波不安定性でどんな渦が形成されるのか 最終的には渦は1つ 方位角平均密度分布の形が緑領域のパラメータを持つ 高い A 平均面密度分布の進化線 レイリー不安定 低い 小さい w 大きい Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 17

18 結果まとめ 1. 線形安定性解析の検証 不安定領域では渦形成が起こる 中立安定 安定領域では渦形成は起こらない 2. ロスビー波不安定性でどんな渦が形成されるのか 最終的には渦は1つ 方位角平均密度分布の形が緑領域のパラメータを持つ 高い A 平均面密度分布の進化線 レイリー不安定 低い 小さい w 大きい Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 18

19 渦中心の面密度 渦の強さ 面積 実際の観測で重要なのは渦中心での面密度 渦の面積 数値シミュレーションの結果できる渦のアスペクト比 (χ b/a) は常に χ~6-8 であった 渦のアスペクト比が一定の時 渦中心での面密度 A/ w 渦の面積 w 2 φ 高い b a 形成される渦の分布 r 小さいが高い渦 大きいが低い渦は作ることができない 低い 小さい 渦の面積 大きい Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 19

20 渦中心の面密度 観測への示唆 ロスビー波不安定性でできる渦は サイズが大きいものが選択的に観測されると期待される ロスビー波不安定性によってできた渦なら緑領域の中に入る 緑領域に入らない渦はロスビー波不安定性で作れない 他のメカニズムか他の物理効果が必要 高い 形成される渦の分布 低い 小さい 渦の面積 大きい Dec. 20th, 2016 理論懇シンポジウム 東北大 20

21 まとめ & 課題 まとめ ロスビー波不安定性について数値シミュレーションを用いて調べた Ono+16 の線形安定性解析が正しいことを確認した ロスビー波不安定性によってできた渦が満たすべき 方位角平均密度分布に制限をつけた 課題 異なる背景 熱的性質を持つ円盤への拡張 実際の観測モデルへの適用 3D, 粘性, 自己重力, 磁場などの影響を考慮 より深い物理的解釈 Dec. 20th, 2016 理論懇シンポジウム 東北大 21

2011 年度第 41 回天文 天体物理若手夏の学校 2011/8/1( 月 )-4( 木 ) 星間現象 18b 初代星形成における水素分子冷却モデルの影響 平野信吾 ( 東京大学 M2) 1. Introduction 初代星と水素分子冷却ファーストスター ( 初代星, PopIII) は重元素を

2011 年度第 41 回天文 天体物理若手夏の学校 2011/8/1( 月 )-4( 木 ) 星間現象 18b 初代星形成における水素分子冷却モデルの影響 平野信吾 ( 東京大学 M2) 1. Introduction 初代星と水素分子冷却ファーストスター ( 初代星, PopIII) は重元素を 2011 年度第 41 回天文 天体物理若手夏の学校 2011/8/1( 月 )-4( 木 ) 星間現象 18b 初代星形成における水素分子冷却モデルの影響 平野信吾 ( 東京大学 M2) 1. Introduction 初代星と水素分子冷却ファーストスター ( 初代星, PopIII) は重元素を含まない原始ガスから形成される 宇宙で最初に誕生する星である 初代星はその後の星形成や再電離など宇宙初期の天文現象に強く関係し

More information

有限密度での非一様なカイラル凝縮と クォーク質量による影響

有限密度での非一様なカイラル凝縮と  クォーク質量による影響 空間的に非一様なカイラル凝縮に対する current quark mass の影響 東京高専 前段眞治 東京理科大学セミナー 2010.9.6 1 1.Introduction 低温 高密度における QCD の振る舞い 中性子星 compact star クォーク物質の理解に重要 T 0 での QCD の基底状態 カイラル対称性の破れた相 カラー超伝導相 μ 2 有限密度において fermionic

More information

領域シンポ発表

領域シンポ発表 1 次元の減衰運動の中の強制振動 ) ( f d d d d d e f e ce ) ( si ) ( 1 ) ( cos ω =ω -γ とおくと 一般解は 外力 f()=f siω の場合 f d d d d si f ce f ce si ) cos( cos si ) cos( この一般解は 1 φ は外力と変位との間の位相差で a 時間が経つと 第 1 項は無視できる この場合の振幅を

More information

超新星残骸Cassiopeia a と 非球対称爆発

超新星残骸Cassiopeia  a と 非球対称爆発 物理学専攻 松尾康秀 宇宙物理理論 指導教員 : 橋本正章 < 超新星残骸 > 星の外層が超新星爆発により吹き飛ばされ 爆発の際の衝撃波によって周囲の物質 ( 星周物質 ) を加熱し 輝いている天体 かに星雲 Kepler Cas A http://www.spacetelescope.o rg/images/large/heic0515a.j pg http://apod.nasa.gov/apod/i

More information

観測的宇宙論WS2013.pptx

観測的宇宙論WS2013.pptx ì コンテンツ イントロダクション 球対称崩壊モデル ビリアル平衡 結果 まとめ イントロダクション 宇宙磁場 銀河や銀河団など様々なスケールで磁場が存在 起源や進化について未だに謎が多い 宇宙の構造形成に影響 P(k)[h -3 Mpc 3 ] 10 6 10 5 10 4 10 3 10 10 1 10 0 10-1 10-10 -3 10-4 10-4 10-3 10-10 -1 10 0 10

More information

銀河風の定常解

銀河風の定常解 2011年 国立天文台プラズマセミナー 2011/12/02 球対称定常銀河風の遷音速解 銀河の質量密度分布との関係 筑波大学 教育研究科 教科教育専攻 つちや まさみ 理科教育コース 2年 土屋 聖海 共同研究者 森正夫 筑波大学 新田伸也 筑波技術大学 発表の流れ はじめに 銀河風とは 流出過程 エネルギー源 周囲に及ぼす影響 研究内容 問題の所在 研究の目的 方法 理論 銀河の質量密度分布 研究成果

More information

スライド 1

スライド 1 非線形数理秋の学校 パターン形成の数理とその周辺 - 反応拡散方程式理論による時 空間パターンの解析を中心に - 2007 年 9 月 25 日 -27 日 モデル方程式を通してみるパターン解析ー進行波からヘリカル波の分岐を例としてー 池田勉 ( 龍谷大学理工学部 ) 講義概要, 講義資料, 講義中に使用する C 言語プログラムと初期値データ, ヘリカル波のアニメーションをウェブで公開しています :

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

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

様々なミクロ計量モデル†

様々なミクロ計量モデル† 担当 : 長倉大輔 ( ながくらだいすけ ) この資料は私の講義において使用するために作成した資料です WEB ページ上で公開しており 自由に参照して頂いて構いません ただし 内容について 一応検証してありますが もし間違いがあった場合でもそれによって生じるいかなる損害 不利益について責任を負いかねますのでご了承ください 間違いは発見次第 継続的に直していますが まだ存在する可能性があります 1 カウントデータモデル

More information

ポリトロープ、対流と輻射、時間尺度

ポリトロープ、対流と輻射、時間尺度 宇宙物理学 ( 概論 ) 6/6/ 大阪大学大学院理学研究科林田清 ポリトロープ関係式 1+(1/) 圧力と密度の間にP=Kρ という関係が成り立っていると仮定する K とは定数でをポリトロープ指数と呼ぶ 5 = : 非相対論的ガス dlnp 3 断熱変化の場合 断熱指数 γ, と dlnρ 4 = : 相対論的ガス 3 1 = の関係にある γ 1 等温変化の場合は= に相当 一様密度の球は=に相当

More information

Microsoft PowerPoint - siryo7

Microsoft PowerPoint - siryo7 . 化学反応と溶液 - 遷移状態理論と溶液論 -.. 遷移状態理論 と溶液論 7 年 5 月 5 日 衝突論と遷移状態理論の比較 + 生成物 原子どうしの反応 活性錯体 ( 遷移状態 ) は 3つの並進 つの回転の自由度をもつ (1つの振動モードは分解に相当 ) 3/ [ ( m m) T] 8 IT q q π + π tansqot 3 h h との並進分配関数 [ πmt] 3/ [ ] 3/

More information

大気環境シミュレーション

大気環境シミュレーション 第 3 回 (Q) 各自 eelを用いて 次の漸化式 + = の解の初期値依存性を調べよ.は50まで () 0 =.0 () 0 =.5 (3) 0 =.0 締切 04 年 月 6 日 ( 月 ) 夕方まで 提出先 347 室 オーバーフロー失敗ゴメンなさい (Q) 各自 eelを用いて 次の漸化式 + = の解の初期値依存性を調べよ.は50まで () 0 =.330 () 0 =.33 (3) 0

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

プラズマ バブルの到達高度に関する研究 西岡未知 齊藤昭則 ( 京都大学理学研究科 ) 概要 TIMED 衛星搭載の GUVI によって観測された赤道異常のピーク位置と 地上 GPS 受信機網によって観測されたプラズマ バブルの出現率や到達率の関係を調べた 高太陽活動時と低太陽活動時について アジア

プラズマ バブルの到達高度に関する研究 西岡未知 齊藤昭則 ( 京都大学理学研究科 ) 概要 TIMED 衛星搭載の GUVI によって観測された赤道異常のピーク位置と 地上 GPS 受信機網によって観測されたプラズマ バブルの出現率や到達率の関係を調べた 高太陽活動時と低太陽活動時について アジア プラズマ バブルの到達高度に関する研究 西岡未知 齊藤昭則 ( 京都大学理学研究科 ) 概要 TIMED 衛星搭載の GUVI によって観測された赤道異常のピーク位置と 地上 GPS 受信機網によって観測されたプラズマ バブルの出現率や到達率の関係を調べた 高太陽活動時と低太陽活動時について アジア地域とアメリカ地域においてそれらの関係を調べたところ 赤道異常高度とプラズマ バブルの出現頻度に強い相関が見られたのは

More information

A Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical

A Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical A Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical calculation method of the gradient as a differential

More information

DVIOUT

DVIOUT 第 章 離散フーリエ変換 離散フーリエ変換 これまで 私たちは連続関数に対するフーリエ変換およびフーリエ積分 ( 逆フーリエ変換 ) について学んできました この節では フーリエ変換を離散化した離散フーリエ変換について学びましょう 自然現象 ( 音声 ) などを観測して得られる波 ( 信号値 ; 観測値 ) は 通常 電気信号による連続的な波として観測機器から出力されます しかしながら コンピュータはこの様な連続的な波を直接扱うことができないため

More information

ÿþŸb8bn0irt

ÿþŸb8bn0irt 折戸の物理 演習プリント N.15 43. 目的 : 電磁誘導は, 基本を理解すれば問題はそれほど難しくない! ということを学ぶ 問 1 の [ ] に適切な数値または数式を入れ, 問 に答えよ 図 1 のように, 紙面に垂直で一様な磁界が 0 の領域だけにある場合について考える 磁束密度は Wb/m で, 磁界は紙面の表から裏へ向かっている 図のように,1 辺の長さが m の正方形のコイル を,

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/ ujio@aori.u-tokyo.ac.jp 2016/1/8 順圧流の運動方程式 流体の密度が一様ならば, 圧力 静水圧 の水平勾配は鉛直一様 海面の高さによる水平圧力勾配のみ ηx,y px, y, z = ρ g dz = ρgη z p x

More information

2 Hermite-Gaussian モード 2-1 Hermite-Gaussian モード 自由空間を伝搬するレーザ光は次のような Hermite-gaussian Modes を持つ光波として扱う ことができる ここで U lm (x, y, z) U l (x, z)u m (y, z) e

2 Hermite-Gaussian モード 2-1 Hermite-Gaussian モード 自由空間を伝搬するレーザ光は次のような Hermite-gaussian Modes を持つ光波として扱う ことができる ここで U lm (x, y, z) U l (x, z)u m (y, z) e Wavefront Sensor 法による三角共振器のミスアラインメント検出 齊藤高大 新潟大学大学院自然科学研究科電気情報工学専攻博士後期課程 2 年 214 年 8 月 6 日 1 はじめに Input Mode Cleaner(IMC) は Fig.1 に示すような三角共振器である 懸架鏡の共振などにより IMC を構成する各ミラーが角度変化を起こすと 入射光軸と共振器軸との間にずれが生じる

More information

ニュートン重力理論.pptx

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

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

untitled

untitled 熱対流現象 山中透 2005 年 3 月 概要 流体を熱源に接触させ, 流体に温度傾度を与えたときを考える. 流体の温度傾度が小さいときは, 熱拡散のみが起こるが, 流体の温度傾度が閾値を越えると, 熱拡散だけでは温度傾度を解消できなくなって不安定となり, 対流が生じる. これをベナール対流とよぶ. ここでは, ベナール対流を記述する非線型方程式の線型安定性の解析によって, 流体が不安定化する条件を求め,

More information

オープン CAE シンポジウム @ 名古屋 C17 遷 速における OpenFOAM の圧縮性ソルバーの 較 2017 年 12 9 松原 輔 ( オープンCAE 勉強会 @ 関 ) 1 お詫びと訂正 講演概要集で誤記がありました 記載されている計算結果は 粘性 の速度の発散項はupwind で った結果となっております 境界条件にも誤記があります ( 後に します ) 申し訳ありません 2 目次

More information

D論研究 :「表面張力対流の基礎的研究」

D論研究 :「表面張力対流の基礎的研究」 D 論研究 : 表面張力対流の基礎的研究 定常 Marangoni 対流 及び非定常 Marangoni 対流に関する実験及び数値解析による検討 Si 単結晶の育成装置 Cz 法による Si 単結晶育成 FZ 法による Si 単結晶育成 気液表面 るつぼ加熱 気液表面 大きな温度差を有す気液表面では表面張力対流 (Marangoni 対流 ) が顕著 プロセス終了後のウエハ Cz 法により育成した

More information

フィードバック ~ 様々な電子回路の性質 ~ 実験 (1) 目的実験 (1) では 非反転増幅器の増幅率や位相差が 回路を構成する抵抗値や入力信号の周波数によってどのように変わるのかを調べる 実験方法 図 1 のような自由振動回路を組み オペアンプの + 入力端子を接地したときの出力電圧 が 0 と

フィードバック ~ 様々な電子回路の性質 ~ 実験 (1) 目的実験 (1) では 非反転増幅器の増幅率や位相差が 回路を構成する抵抗値や入力信号の周波数によってどのように変わるのかを調べる 実験方法 図 1 のような自由振動回路を組み オペアンプの + 入力端子を接地したときの出力電圧 が 0 と フィードバック ~ 様々な電子回路の性質 ~ 実験 (1) 目的実験 (1) では 非反転増幅器の増幅率や位相差が 回路を構成する抵抗値や入力信号の周波数によってどのように変わるのかを調べる 実験方法 図 1 のような自由振動回路を組み オペアンプの + 入力端子を接地したときの出力電圧 が 0 となるように半固定抵抗器を調整する ( ゼロ点調整のため ) 図 1 非反転増幅器 2010 年度版物理工学実験法

More information

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

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

More information

コロイド化学と界面化学

コロイド化学と界面化学 環境表面科学講義 http://res.tagen.tohoku.ac.jp/~liquid/mura/kogi/kaimen/ E-mail: mura@tagen.tohoku.ac.jp 村松淳司 分散と凝集 ( 平衡論的考察! 凝集! van der Waals 力による相互作用! 分散! 静電的反発力 凝集 分散! 粒子表面の電位による反発 分散と凝集 考え方! van der Waals

More information

太陽系における地球型惑星の水の起源ーー惑星形成の大域シミュレーションーー 小南淳子 ( 東京工業大学地球生命研究所 ) 台坂博 ( 一橋 ) 似鳥啓吾 ( 理研 ) 牧野淳一郎 ( 東工大 )

太陽系における地球型惑星の水の起源ーー惑星形成の大域シミュレーションーー 小南淳子 ( 東京工業大学地球生命研究所 ) 台坂博 ( 一橋 ) 似鳥啓吾 ( 理研 ) 牧野淳一郎 ( 東工大 ) 太陽系における地球型惑星の水の起源ーー惑星形成の大域シミュレーションーー 小南淳子 ( 東京工業大学地球生命研究所 ) 台坂博 ( 一橋 ) 似鳥啓吾 ( 理研 ) 牧野淳一郎 ( 東工大 ) 太陽系 http://www.astroarts.co.jp/news/2006/08/28planet_5/ 惑星の形成シナリオ 微惑星 (~km サイズ ) 暴走成長 寡占成長 原始惑星 ( 火星サイズ

More information

Microsoft PowerPoint - 2_6_shibata.ppt [互換モード]

Microsoft PowerPoint - 2_6_shibata.ppt [互換モード] 圧密問題への逆問題の適用 一次元圧密と神戸空港の沈下予測 1. 一次元圧密の解析 2. 二次元圧密問題への適用 3. 神戸空港の沈下予測 1. 一次元圧密の解析 一次元圧密の実験 試験システムの概要 分割型圧密試験 逆解析の条件 未知量 ( 同定パラメータ ) 圧縮指数 :, 透水係数 :k 初期体積ひずみ速度 : 二次圧密係数 : 観測量沈下量 ( 計 4 点 ) 逆解析手法 粒子フィルタ (SIS)

More information

計算機シミュレーション

計算機シミュレーション . 運動方程式の数値解法.. ニュートン方程式の近似速度は, 位置座標 の時間微分で, d と定義されます. これを成分で書くと, d d li li とかけます. 本来は が の極限をとらなければいけませんが, 有限の小さな値とすると 秒後の位置座標は速度を用いて, と近似できます. 同様にして, 加速度は, 速度 の時間微分で, d と定義されます. これを成分で書くと, d d li li とかけます.

More information

19年度一次基礎科目計算問題略解

19年度一次基礎科目計算問題略解 9 年度機械科目 ( 計算問題主体 ) 略解 基礎科目の解析の延長としてわかる範囲でトライしてみたものです Coprigh (c) 7 宮田明則技術士事務所 Coprigh (c) 7 宮田明則技術士事務所 Ⅳ- よってから は許容荷重として は直径をロ - プの断面積 Ⅳ- cr E E E I, から Ⅳ- Ⅳ- : q q q q q q q q q で絶対値が最大 で絶対値が最大モーメントはいずれも中央で最大となる

More information

三重大学工学部

三重大学工学部 反応理論化学 ( その 軌道相互作用 複数の原子が相互作用して分子が形成される複数の原子軌道 ( または混成軌道 が混合して分子軌道が形成される原子軌道 ( または混成軌道 が混合して分子軌道に変化すると軌道エネルギーも変化する. 原子軌道 原子軌道は3つの量子数 ( nlm,, の組合せにより指定される量子数の取り得る値の範囲 n の値が定まる l の範囲は n の値に依存して定まる m の範囲は

More information

PowerPoint Presentation

PowerPoint Presentation 付録 2 2 次元アフィン変換 直交変換 たたみ込み 1.2 次元のアフィン変換 座標 (x,y ) を (x,y) に移すことを 2 次元での変換. 特に, 変換が と書けるとき, アフィン変換, アフィン変換は, その 1 次の項による変換 と 0 次の項による変換 アフィン変換 0 次の項は平行移動 1 次の項は座標 (x, y ) をベクトルと考えて とすれば このようなもの 2 次元ベクトルの線形写像

More information

電磁気学 A 練習問題 ( 改 ) 計 5 ページ ( 以下の問題およびその類題から 3 題程度を定期試験の問題として出題します ) 以下の設問で特に断らない限り真空中であることが仮定されているものとする 1. 以下の量を 3 次元極座標 r,, ベクトル e, e, e r 用いて表せ (1) g

電磁気学 A 練習問題 ( 改 ) 計 5 ページ ( 以下の問題およびその類題から 3 題程度を定期試験の問題として出題します ) 以下の設問で特に断らない限り真空中であることが仮定されているものとする 1. 以下の量を 3 次元極座標 r,, ベクトル e, e, e r 用いて表せ (1) g 電磁気学 A 練習問題 ( 改 ) 計 5 ページ ( 以下の問題およびその類題から 題程度を定期試験の問題として出題します ) 以下の設問で特に断らない限り真空中であることが仮定されているものとする. 以下の量を 次元極座標,, ベクトル e, e, e 用いて表せ () gad () ot A (). 以下の量を 次元円柱座標,, z 位ベクトル e e, e, z 用いて表せ () gad ()

More information

第6章 実験モード解析

第6章 実験モード解析 第 6 章実験モード解析 6. 実験モード解析とは 6. 有限自由度系の実験モード解析 6.3 連続体の実験モード解析 6. 実験モード解析とは 実験モード解析とは加振実験によって測定された外力と応答を用いてモードパラメータ ( 固有振動数, モード減衰比, 正規固有モードなど ) を求める ( 同定する ) 方法である. 力計 試験体 変位計 / 加速度計 実験モード解析の概念 時間領域データを利用する方法

More information

OpenFOAM 掲示版のまとめ 2012/12/01 富山県立大学中川慎二

OpenFOAM 掲示版のまとめ 2012/12/01 富山県立大学中川慎二 OpenFOAM 掲示版のまとめ 2012/12/01 富山県立大学中川慎二 Q1. 管内流の周期境界条件 パイプ内部の流れを解析するとき, 上流の流入面と下流の流出面を周期境界条件として, 発達した流れを計算したい 単純に cyclic 境界を使うと, 流入面と流出面とが同一圧力になり, 流れがなくなってしまう どうすれば良いか? A1-1. 管内流の周期境界条件 cyclicjump から派生した

More information

PowerPoint Presentation

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

More information

観測的宇宙論workshop.pptx

観測的宇宙論workshop.pptx 名古屋 大学宇宙論論研究室 嵯峨承平 ( 共同研究者 : 市來來淨與, 杉 山直 ) 2013/12/4 観測的宇宙論論 workshop 1/20 目次 1. イントロ 2. 2 次摂動論論 3. 重 力力波 ( 線形摂動 ) 4. 重 力力波 (2 次摂動 ) 5. まとめ 2/20 1. イントロ 非ガウス性 重 力力レンズ効果 2 次ドップラー効果 2 次重 力力波 磁場 Mode coupling

More information

スーパー地球の熱進化と 磁場の寿命 立浪千尋 千秋博紀 井田茂 衛星系形成小研究会 2012 小樽

スーパー地球の熱進化と 磁場の寿命 立浪千尋 千秋博紀 井田茂 衛星系形成小研究会 2012 小樽 スーパー地球の熱進化と 磁場の寿命 立浪千尋 千秋博紀 井田茂 衛星系形成小研究会 2012 夏 @ 小樽 地球型惑星 岩石マントル 金属コア 岩石マントル 金属コア (e.g. Ida and Lin, 2008) HARPS CoRoT Kepler 観測された系外惑星と スーパー地球候補 赤 : トランジット法緑 : 視線速度法 惑星質量 ( 地球質量 ) 平均密度 (g/cm 3 ) 軌道長半径

More information

高軌道傾斜角を持つメインベルト 小惑星の可視光分光観測

高軌道傾斜角を持つメインベルト 小惑星の可視光分光観測 高軌道傾斜角を持つメインベルト小惑星の可視光分光観測 天文 天体物理夏の学校 @ 福井神戸大学 M2 岩井彩 背景 小惑星岩石質の太陽系小天体であり 彗星活動を行わない 分類軌道長半径による空間分布可視光波長域のスペクトル形状 ( 大きく 5 種類 ) 空間分布による分類 メインベルト ( 小惑星帯 ) 太陽から 2.1-3.3AU 離れた環状の領域軌道が確定した小惑星の約 9 割が存在 トロヤ群木星のラグランジュ点

More information

第 4 週コンボリューションその 2, 正弦波による分解 教科書 p. 16~ 目標コンボリューションの演習. 正弦波による信号の分解の考え方の理解. 正弦波の複素表現を学ぶ. 演習問題 問 1. 以下の図にならって,1 と 2 の δ 関数を図示せよ δ (t) 2

第 4 週コンボリューションその 2, 正弦波による分解 教科書 p. 16~ 目標コンボリューションの演習. 正弦波による信号の分解の考え方の理解. 正弦波の複素表現を学ぶ. 演習問題 問 1. 以下の図にならって,1 と 2 の δ 関数を図示せよ δ (t) 2 第 4 週コンボリューションその, 正弦波による分解 教科書 p. 6~ 目標コンボリューションの演習. 正弦波による信号の分解の考え方の理解. 正弦波の複素表現を学ぶ. 演習問題 問. 以下の図にならって, と の δ 関数を図示せよ. - - - δ () δ ( ) - - - 図 δ 関数の図示の例 δ ( ) δ ( ) δ ( ) δ ( ) δ ( ) - - - - - - - -

More information

アンデン株式会社第 1 技術部 DE 開発藤井成樹 < 業務内容 > アンデンとして CAE 解析を強化するために 10/1 月に DE(Degital Engineering) 開発が 5 名で発足 CAE 開発 活用が目的 解析内容は 構造解析 ( 動解析 非線形含む ) 電場 磁場 音場 熱流

アンデン株式会社第 1 技術部 DE 開発藤井成樹 < 業務内容 > アンデンとして CAE 解析を強化するために 10/1 月に DE(Degital Engineering) 開発が 5 名で発足 CAE 開発 活用が目的 解析内容は 構造解析 ( 動解析 非線形含む ) 電場 磁場 音場 熱流 アンデン株式会社第 1 技術部 DE 開発藤井成樹 < 業務内容 > アンデンとして CAE 解析を強化するために 10/1 月に DE(Degital Engineering) 開発が 5 名で発足 CAE 開発 活用が目的 解析内容は 構造解析 ( 動解析 非線形含む ) 電場 磁場 音場 熱流 流体解析など様々 項目 04 05 06 07 08 09 10 11

More information

3. 重力波と沿岸 赤道ケルビン波 2014 年 9 月 30 日 16:35 見延庄士郎 ( 海洋気候物理学研究室 ) 予習課題 : 以下の you tube のビデオを見ておくこと. 個々のビデオは全部は見ずに, 雰囲気がつかめる程度見ればいい.

3. 重力波と沿岸 赤道ケルビン波 2014 年 9 月 30 日 16:35 見延庄士郎 ( 海洋気候物理学研究室 ) 予習課題 : 以下の you tube のビデオを見ておくこと. 個々のビデオは全部は見ずに, 雰囲気がつかめる程度見ればいい. 3. 重力波と沿岸 赤道ケルビン波 2014 年 9 月 30 日 16:35 見延庄士郎 ( 海洋気候物理学研究室 ) minobe@sci.hokudai.ac.jp 予習課題 : 以下の you ube のビデオを見ておくこと. 個々のビデオは全部は見ずに, 雰囲気がつかめる程度見ればいい. 大気の重力波 : hp://www.youube.com/wach?v=yxnkzecu3be 津波シミュレーション

More information

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

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

More information

NLMIXED プロシジャを用いた生存時間解析 伊藤要二アストラゼネカ株式会社臨床統計 プログラミング グループグルプ Survival analysis using PROC NLMIXED Yohji Itoh Clinical Statistics & Programming Group, A

NLMIXED プロシジャを用いた生存時間解析 伊藤要二アストラゼネカ株式会社臨床統計 プログラミング グループグルプ Survival analysis using PROC NLMIXED Yohji Itoh Clinical Statistics & Programming Group, A NLMIXED プロシジャを用いた生存時間解析 伊藤要二アストラゼネカ株式会社臨床統計 プログラミング グループグルプ Survival analysis using PROC NLMIXED Yohji Itoh Clinical Statistics & Programming Group, AstraZeneca KK 要旨 : NLMIXEDプロシジャの最尤推定の機能を用いて 指数分布 Weibull

More information

Microsoft Word _鹿児島用(正式版).docx

Microsoft Word _鹿児島用(正式版).docx 超巨大ブラックホールを取り巻くドーナツ構造の正体を暴く 1/7 概要 国立天文台の泉拓磨氏 鹿児島大学の和田桂一氏を中心とする研究チームは アルマ望遠鏡を使ってコンパス座銀河の中心に位置する超巨大ブラックホールを観測し その周囲のガスの分布と動きをこれまでになく詳細に明らかにすることに成功しました 活動的な超巨大ブラックホールの周囲にはガスや塵のドーナツ状構造が存在すると考えられてきましたが その成因は長年の謎でした

More information

Probit , Mixed logit

Probit , Mixed logit Probit, Mixed logit 2016/5/16 スタートアップゼミ #5 B4 後藤祥孝 1 0. 目次 Probit モデルについて 1. モデル概要 2. 定式化と理解 3. 推定 Mixed logit モデルについて 4. モデル概要 5. 定式化と理解 6. 推定 2 1.Probit 概要 プロビットモデルとは. 効用関数の誤差項に多変量正規分布を仮定したもの. 誤差項には様々な要因が存在するため,

More information

スライド 1

スライド 1 相対論的プラズマにおける PIC シミュレーションに伴う数値チェレンコフ不安定の特性ついて 宇宙物理学研究室 4 年池谷直樹 研究背景と目的 0 年 Ie Cube 国際共同実験において超高エネルギーニュートリノを検出 780Tev-5.6PeV 890TeV-8.5PeV 相互作用が殆んど起こらないため銀河磁場による軌道の湾曲が無く 正確な到来方向の情報 を得られる可能性がある ニュートリノから高エネルギー宇宙線の起源を追う

More information

i i = 1, 2,..., i r, t t r i 1 r, tt i r, t r Σ φ i r, t =1. i =1 φ i t = 2 n r, t Σj i s ij M ij δ F δ F δφ i δφ j, i =1 n F M ij δ F δφ i = Σj i 2 ε

i i = 1, 2,..., i r, t t r i 1 r, tt i r, t r Σ φ i r, t =1. i =1 φ i t = 2 n r, t Σj i s ij M ij δ F δ F δφ i δφ j, i =1 n F M ij δ F δφ i = Σj i 2 ε Phase-field Simulation of Grain Growth Abstract Grain growth occurring in the steel manufacturing process is one of the most important phenomena in order to control the polycrystalline microstructure.

More information

構造力学Ⅰ第12回

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

More information

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

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

More information

³ÎΨÏÀ

³ÎΨÏÀ 2017 12 12 Makoto Nakashima 2017 12 12 1 / 22 2.1. C, D π- C, D. A 1, A 2 C A 1 A 2 C A 3, A 4 D A 1 A 2 D Makoto Nakashima 2017 12 12 2 / 22 . (,, L p - ). Makoto Nakashima 2017 12 12 3 / 22 . (,, L p

More information

横浜市環境科学研究所

横浜市環境科学研究所 周期時系列の統計解析 単回帰分析 io 8 年 3 日 周期時系列に季節調整を行わないで単回帰分析を適用すると, 回帰係数には周期成分の影響が加わる. ここでは, 周期時系列をコサイン関数モデルで近似し単回帰分析によりモデルの回帰係数を求め, 周期成分の影響を検討した. また, その結果を気温時系列に当てはめ, 課題等について考察した. 気温時系列とコサイン関数モデル第 報の結果を利用するので, その一部を再掲する.

More information

Microsoft PowerPoint - H22制御工学I-10回.ppt

Microsoft PowerPoint - H22制御工学I-10回.ppt 制御工学 I 第 回 安定性 ラウス, フルビッツの安定判別 平成 年 6 月 日 /6/ 授業の予定 制御工学概論 ( 回 ) 制御技術は現在様々な工学分野において重要な基本技術となっている 工学における制御工学の位置づけと歴史について説明する さらに 制御システムの基本構成と種類を紹介する ラプラス変換 ( 回 ) 制御工学 特に古典制御ではラプラス変換が重要な役割を果たしている ラプラス変換と逆ラプラス変換の定義を紹介し

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

2

2 1 2 3 4 5 6 7 8 9 10 11 12 13 14 10 6 15 7 13 7 14 13 31 13 13 9 30 13 10 14 13 10 25 13 11 29 13 11 30 14 111 14 223 14 2 23 3 17 14 3 17 14 3 21 14 3 26 13 11 29 5 2 2 14 1 11 1 3 28 1 3 14 2 28 1 28

More information

えられる球体について考えよ 慣性モーメント C と体積 M が以下の式で与えられることを示せ (5.8) (5.81) 地球のマントルと核の密度の平均値を求めよ C= kg m 2, M= kg, a=6378km, rc=3486km 次に (5.82) で与えら

えられる球体について考えよ 慣性モーメント C と体積 M が以下の式で与えられることを示せ (5.8) (5.81) 地球のマントルと核の密度の平均値を求めよ C= kg m 2, M= kg, a=6378km, rc=3486km 次に (5.82) で与えら 5.5 慣性モーメント (5-42) 式で与えられたマッカラーの公式は 扁球惑星体の重力加速度とその主な慣性モーメントを関連づけている その公式を使うことで 探査飛行や軌道上を周回する宇宙船によって 例えば慣性モーメントを束縛している惑星の重力場を計測することができる 慣性モーメントは惑星全体の形や内部の密度分布を反映するため 惑星の内部構造を調べるために慣性モーメントの数値を利用することができる

More information

ブラックホール近傍の相対論的光軌道

ブラックホール近傍の相対論的光軌道 ラックホールに落下するガスの blob 2014 年 2 月 2 日京都大学宇宙物理学教室修士 1 年森山小太郎 本研究 遠方 S 降着円盤 B ガスの塊 ( 以降 spot) 最内縁安定円軌道からずれて BH に落ち込むガスの塊について考える 遠方からどう観測されるか理論的に研究する a の決定に用いる 円運動する Spot からの 光のエネルギーフラックス スタート地点 B S a=0.9981m

More information

サブ課題Cの目標 大規模な宇宙論的構造形成シミュレーションの共分散解析による広域銀 河サーベイの統計解析 (吉田 石山) ブラックホール降着円盤の一般相対論的輻射磁気流体シミュレーション及 びグローバルシミュレーション 松元 大須賀 大規模なプラズマ粒子シミュレーションによる磁気再結合と高エネルギー

サブ課題Cの目標 大規模な宇宙論的構造形成シミュレーションの共分散解析による広域銀 河サーベイの統計解析 (吉田 石山) ブラックホール降着円盤の一般相対論的輻射磁気流体シミュレーション及 びグローバルシミュレーション 松元 大須賀 大規模なプラズマ粒子シミュレーションによる磁気再結合と高エネルギー 多次元高精度ブラソフソルバーの開発 素粒子 原子核 宇宙 京からポスト京に向けて シンポジウム 2017年2月17日 筑波大学 東京キャンパス 筑波大学 計算科学研究センター 吉川 耕司 サブ課題Cの目標 大規模な宇宙論的構造形成シミュレーションの共分散解析による広域銀 河サーベイの統計解析 (吉田 石山) ブラックホール降着円盤の一般相対論的輻射磁気流体シミュレーション及 びグローバルシミュレーション

More information

×××××××××× ×××××××××××××××

×××××××××× ××××××××××××××× Hoizon-penetating Tansonic Accetion Disks aound Rotating Black Holes with Causal Viscosity 高橋労太 ( 東大総合文化 ) ホライズンの内側まで解かれた ADAF の遷音速流のサンプル解 (4 元速度の 成分 ) 要旨 ブラックホール周りの定常降着流の遷音速解を外側の領域からホライズンの中まで計算できるようになった

More information

カイ二乗フィット検定、パラメータの誤差

カイ二乗フィット検定、パラメータの誤差 統計的データ解析 008 008.. 林田清 ( 大阪大学大学院理学研究科 ) 問題 C (, ) ( x xˆ) ( y yˆ) σ x πσ σ y y Pabx (, ;,,, ) ˆ y σx σ y = dx exp exp πσx ただし xy ˆ ˆ はyˆ = axˆ+ bであらわされる直線モデル上の点 ( ˆ) ( ˆ ) ( ) x x y ax b y ax b Pabx (,

More information

T2K 実験 南野彰宏 ( 京都大学 ) 他 T2Kコラボレーション平成 25 年度宇宙線研究所共同利用成果発表会 2013 年 12 月 20 日 1

T2K 実験 南野彰宏 ( 京都大学 ) 他 T2Kコラボレーション平成 25 年度宇宙線研究所共同利用成果発表会 2013 年 12 月 20 日 1 T2K 実験 南野彰宏 ( 京都大学 ) 他 T2Kコラボレーション平成 25 年度宇宙線研究所共同利用成果発表会 2013 年 12 月 20 日 1 T2K 実験 J- PARC でほぼ純粋な ν µμ ビームを生成 生成点直後の前置検出器と 295km 離れたスーパーカミオカンデでニュートリノを観測 ニュートリノ振動の精密測定 T2K 実験における振動モード 1. ν µμ ν e (ν e

More information

物性物理学I_2.pptx

物性物理学I_2.pptx The University of Tokyo, Komaba Graduate School of Arts and Sciences I 凝縮系 固体 をデザインする 銅()面上の鉄原子の 量子珊瑚礁 IBM Almaden 許可を得て掲載 www.almaden.ibm.com/vis/stm/imagesstm5.jpg&imgrefurl=http://www.almaden.ibm.com/vis/

More information

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

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

More information

<4D F736F F D D082B882DD90AC89CA95F18D908F F312D385F895E896388CF2E646F63>

<4D F736F F D D082B882DD90AC89CA95F18D908F F312D385F895E896388CF2E646F63> 3.1.8 伸張場におけるひずみ集中メカニズムに関する研究 業務の目的日本海拡大による伸張とその後の圧縮によって形成されたひずみ集中帯の機構を解明するためには 圧縮場のみならず伸張場における応力の境界条件とこれに対する地殻の応答を比較することが重要であるため 現在日本列島の中でも顕著な伸張場を形成している別府島原地溝帯を含む九州地域において 地殻の挙動を正確に把握する検知システム 及び局所的な構造を調べるための観測システムを整備し

More information

ÿþŸb8bn0irt

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

More information

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

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

More information

WTENK5-6_26265.pdf

WTENK5-6_26265.pdf 466 2014年秋季 極域 寒冷域研究連絡会 の報告 海 カラ海 北大西洋 北米大陸の北部 東アジアで が多重に見られることが多い 南極昭和基地 69.0 S, 寒気質量の減少傾向が 中央シベリアの内陸部とベー 39.6 E における PANSY レーダー Sato et al.2014 リング海で寒気質量の増加傾向が5つの再解析データ のデータは このような小さな に共通して見られた 中央シベリアの内陸部の寒気質

More information

Microsoft PowerPoint - 第3回2.ppt

Microsoft PowerPoint - 第3回2.ppt 講義内容 講義内容 次元ベクトル 関数の直交性フーリエ級数 次元代表的な対の諸性質コンボリューション たたみこみ積分 サンプリング定理 次元離散 次元空間周波数の概念 次元代表的な 次元対 次元離散 次元ベクトル 関数の直交性フーリエ級数 次元代表的な対の諸性質コンボリューション たたみこみ積分 サンプリング定理 次元離散 次元空間周波数の概念 次元代表的な 次元対 次元離散 ベクトルの直交性 3

More information

Microsoft PowerPoint - ce07-13b.ppt

Microsoft PowerPoint - ce07-13b.ppt 制御工学 3 第 8 章 : フィードバック制御系の設計法 8. 設計手順と性能評価キーワード : 設計手順, 性能評価 8. ID 補償による制御系設計キーワード : ( 比例 ),I( 積分 ),D( 微分 ) 8.3 進み 遅れ補償による制御系設計キーワード : 遅れ補償, 進み補償 学習目標 : 一般的な制御系設計における手順と制御系の性能評価について学ぶ. ループ整形の考え方を用いて, 遅れ補償,

More information

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

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

More information

Microsoft PowerPoint - elast.ppt [互換モード]

Microsoft PowerPoint - elast.ppt [互換モード] 弾性力学入門 年夏学期 中島研吾 科学技術計算 Ⅰ(48-7) コンピュータ科学特別講義 Ⅰ(48-4) elast 弾性力学 弾性力学の対象 応力 弾性力学の支配方程式 elast 3 弾性力学 連続体力学 (Continuum Mechanics) 固体力学 (Solid Mechanics) の一部 弾性体 (lastic Material) を対象 弾性論 (Theor of lasticit)

More information

Microsoft PowerPoint - chap8.ppt

Microsoft PowerPoint - chap8.ppt 第 8 章 : フィードバック制御系の設計法 第 8 章 : フィードバック制御系の設計法 8. 設計手順と性能評価 キーワード : 設計手順, 性能評価 8. 補償による制御系設計 キーワード : ( 比例 ),( 積分 ),( 微分 ) 学習目標 : 一般的な制御系設計における手順と制御系の性能評価について学ぶ. 学習目標 : 補償の有効性について理解し, その設計手順を習得する. 第 8 章

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 - 4_構造特性係数の設定方法に関する検討.doc

Microsoft Word - 4_構造特性係数の設定方法に関する検討.doc 第 4 章 構造特性係数の設定方法に関する検討 4. はじめに 平成 年度 年度の時刻歴応答解析を実施した結果 課題として以下の点が指摘 された * ) 脆性壁の評価法の問題 時刻歴応答解析により 初期剛性が高く脆性的な壁については現在の構造特性係数 Ds 評価が危険であることが判明した 脆性壁では.5 倍程度必要保有耐力が大きくなる * ) 併用構造の Ds の設定の問題 異なる荷重変形関係を持つ壁の

More information

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

国土技術政策総合研究所 研究資料 1. 概要本資料は, 重力式コンクリートダムの地震時における挙動の再現性を 多数の地震計が配置され 地震記録を豊富に有する札内川ダムをモデルダムとして三次元地震応答解析を実施したものである 実際に観測された加速度時刻歴波形から 地震時における構造物 - 貯水池 - 基礎岩盤の相互作用を考慮して実施した三次元応答解析の結果と実際の地震時の観測結果を比較することで当該ダムの地震時の物性値を同定した上で

More information

syuryoku

syuryoku 248 24622 24 P.5 EX P.212 2 P271 5. P.534 P.690 P.690 P.690 P.690 P.691 P.691 P.691 P.702 P.702 P.702 P.702 1S 30% 3 1S 3% 1S 30% 3 1S 3% P.702 P.702 P.702 P.702 45 60 P.702 P.702 P.704 H17.12.22 H22.4.1

More information

土壌環境行政の最新動向(環境省 水・大気環境局土壌環境課)

土壌環境行政の最新動向(環境省 水・大気環境局土壌環境課) 201022 1 18801970 19101970 19201960 1970-2 1975 1980 1986 1991 1994 3 1999 20022009 4 5 () () () () ( ( ) () 6 7 Ex Ex Ex 8 25 9 10 11 16619 123 12 13 14 5 18() 15 187 1811 16 17 3,000 2241 18 19 ( 50

More information

Microsoft Word - thesis.doc

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

More information

Microsoft PowerPoint - H24 aragane.pptx

Microsoft PowerPoint - H24 aragane.pptx 海上人工島の経年品質変化 研究背景 目的 解析条件 ( 境界条件 構成モデル 施工履歴 材料パラメータ ) 実測値と解析値の比較 ( 沈下量 ) 将来の不等沈下予測 ケーススタディー ( 埋土施工前に地盤改良を行う : 一面に海上 SD を打設 ) 研究背景 目的 解析条件 ( 境界条件 構成モデル 施工履歴 材料パラメータ ) 実測値と解析値の比較 ( 沈下量 ) 将来の不等沈下予測 ケーススタディー

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

宇宙機工学 演習問題

宇宙機工学 演習問題 宇宙システム工学演習 重力傾度トルク関連. 図に示すように地球回りの円軌道上を周回する宇宙機の運動 を考察する 地球中心座標系を 系 { } 軌道面基準回転系を 系 { } 機体固定系を 系 { } とする 特に次の右手直交系 : 地心方向単位ベクトル 軌道面内 : 進行方向単位ベクトル 軌道面内 : 面外方向単位ベクトル 軌道面外 を取る 特に この { } Lol Horiotl frme と呼ぶ

More information

Microsoft PowerPoint - 20111029第8回ライトカーブ研究会_並木則行.ppt

Microsoft PowerPoint - 20111029第8回ライトカーブ研究会_並木則行.ppt 小 惑 星 探 査 における ダストその 場 観 察 並 木 則 行, 小 林 正 規, 千 秋 博 紀, 和 田 浩 二 2011 年 10 月 29 日 第 8 回 ライトカーブ 研 究 会 国 立 天 文 台 三 鷹 キャンパス 本 日 の 予 定 惑 星 探 査 研 究 センターの 紹 介 小 惑 星 ダスト 観 測 の 科 学 はやぶさの 成 果 を 発 展 させるために, 次 に 日 本

More information

で通常 0.1mm 程度であるのに対し, 軸受内部の表面の大きさは通常 10mm 程度であり, 大きさのスケールが100 倍程度異なる. 例えば, 本研究で解析対象とした玉軸受について, すべての格子をEHLに用いる等間隔構造格子で作成したとすると, 総格子点数は10,000,000のオーダーとなる

で通常 0.1mm 程度であるのに対し, 軸受内部の表面の大きさは通常 10mm 程度であり, 大きさのスケールが100 倍程度異なる. 例えば, 本研究で解析対象とした玉軸受について, すべての格子をEHLに用いる等間隔構造格子で作成したとすると, 総格子点数は10,000,000のオーダーとなる 論文の内容の要旨 論文題目 転がり軸受における枯渇弾性流体潤滑とマクロ流れのマルチスケール連成解析手法の開発 氏名柴﨑健一 転がり軸受は, 転動体が, 外輪および内輪上の溝を転がることにより, 軸を回転自在に支持する機械要素であり, 長寿命化, 低摩擦化が強く求められている. 軸受の摩耗や焼付を防ぎ, 寿命を延ばすため, 通常は潤滑油またはグリースなどの潤滑剤が用いられる. 潤滑油は, 転がり接触する二表面間に表面粗さよりも厚い膜を形成し,

More information

ÿþŸb8bn0irt

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

More information

Microsoft PowerPoint _量子力学短大.pptx

Microsoft PowerPoint _量子力学短大.pptx . エネルギーギャップとrllouゾーン ブリルアン領域,t_8.. 周期ポテンシャル中の電子とエネルギーギャップ 簡単のため 次元に間隔 で原子が並んでいる結晶を考える 右方向に進行している電子の波は 間隔 で規則正しく並んでいる原子が作る格子によって散乱され 左向きに進行する波となる 波長 λ が の時 r の反射条件 式を満たし 両者の波が互いに強め合い 定在波を作る つまり 式 式を満たす波は

More information

Microsoft PowerPoint - 9.pptx

Microsoft PowerPoint - 9.pptx 9. 線形写像 ここでは 行列の積によって 写像を定義できることをみていく また 行列の積によって定義される写像の性質を調べていく 行列演算と写像 ( 次変換 3 拡大とスカラー倍 p ' = ( ', ' = ( k, kk p = (, k 倍 k 倍 拡大後 k 倍拡大の関係は スカラー倍を用いて次のように表現できる ' = k ' 拡大前 拡大 4 拡大と行列の積 p ' = ( ', '

More information

Microsoft PowerPoint - 9.pptx

Microsoft PowerPoint - 9.pptx 9/7/8( 水 9. 線形写像 ここでは 行列の積によって 写像を定義できることをみていく また 行列の積によって定義される写像の性質を調べていく 拡大とスカラー倍 行列演算と写像 ( 次変換 拡大後 k 倍 k 倍 k 倍拡大の関係は スカラー倍を用いて次のように表現できる p = (, ' = k ' 拡大前 p ' = ( ', ' = ( k, k 拡大 4 拡大と行列の積 拡大後 k 倍

More information

Microsoft Word - 力学12.doc

Microsoft Word - 力学12.doc 慣性モーメント. 復習 角運動量と角速度 L p υ, L 質点の角運動量 : ( ) ( ) 剛体の角運動量 L ( ) ρ ( ) ( ) d 注 ) この積分は普通の三重積分 d d d ( ) ( ) A B C A C B A B より ベクトル三重積の公式 ( ) ( ) ( )C ( ) L ( ) ( ) R 但し 慣性モーメント (oent of net): I R( ) ρ ;

More information