スライド 1

Size: px
Start display at page:

Download "スライド 1"

Transcription

1 相対論的プラズマにおける PIC シミュレーションに伴う数値チェレンコフ不安定の特性ついて 宇宙物理学研究室 4 年池谷直樹

2 研究背景と目的 0 年 Ie Cube 国際共同実験において超高エネルギーニュートリノを検出 780Tev-5.6PeV 890TeV-8.5PeV 相互作用が殆んど起こらないため銀河磁場による軌道の湾曲が無く 正確な到来方向の情報 を得られる可能性がある

3 ニュートリノから高エネルギー宇宙線の起源を追う ニュートリノ生成過程 超高エネルギー宇宙線 宇宙背景輻射 超高エネルギーニュートリノ 効率よく粒子を超高エネルギーまで加速させる機構 が存在するはず 無衝突衝撃波における粒子加速に注目する

4 無衝突衝撃波 平均自由行程 >> 電磁場による散逸のサイズ このようなプラズマを無衝突プラズマといい 無衝突プラズマ中を伝搬する衝撃波を無衝突衝撃波という 相対論的な速度を持つ無衝突衝撃波 AGN 活動銀河核 ジェット 線バースト 超高エネルギー粒子の起源と考えられている

5 考えられる加速機構 衝撃波統計加速 フェルミ 次加速 上流 下流 電磁場による粒子の散乱 衝撃波面を往復 粒子 往復の度にエネルギーを得る 相対論的無衝突衝撃波における物理過程電磁場とプラズマを構成する粒子が複雑に関係する非線形現象 数値シミュレーションによる解析が有効

6 シミュレーション法 フェルミ 次加速のような一部の粒子を高エネルギーまで加速させる非断熱的なプロセスの理解には PIC Parie-In-Ce 法が有効 磁気流体近似 MHD 法 粒子シミュレーション PIC 法 粒子をまとまった系として みなし 流体として扱う 熱的な系に有効 個々の粒子の運動方程式を解く

7 PIC シミュレーション法の特徴粒子 相対論的な運動方程式電磁場 グリット上での Mawe 方程式 B u E m q d du E B J B E 4 B 0 v γ CGS 系 E 4

8 Numeria Cherenov Radiaion PIC 法において電磁場をグリッド上で解く際 電磁波の位相速度が光速より落ちる 数値分散 B 相対論的な系を扱う際 数値チェレンコフ放射と呼ばれる 数値的な不安定性が発生する / 例 方向に速度 v のプラズマを流す / 数値チェレンコフ放射発生の理論と その対処法を追う

9 なぜ数値チェレンコフ? 自然界 粒子速度 > 光速 v 干渉 = チェレンコフ光発生 数値上 プラズマ速度 > 電磁波の位相速度 v 数値不安定として発生 = 数値チェレンコフ放射

10 Mawe 方程式の差分解法に伴う数値分散 PIC 法における標準的な電磁場の差分解法 FDTD 法陽的 epii な解法 B B E E B B E E E E E E B B m n m n m n m n m n n m n m n m m n m n n m n m m n m n E B E m n n m

11 FDTD 法における電磁場の位相速度の分散関係離散化させた平面波の式を代入することで得られる FDTD 法における分散関係分散公式と異なる形で現れる ep ep 0 0 n m i E E n m i B B in in arin

12 pcans では FDTD 法と異なる陰的 Impii な解法を採用している +/ における物理量を + と における情報より定義している 陰的解法による 分散関係は F F F シミュレーションパッケージ pcans における Mawe 方程式の差分解法 in in aran

13 分散関係の比較パラメータ FDTD 法 epii arin in in pcans 陰的解法 Impii aran in in

14 分散関係の比較 パラメータ epii impii

15 Numeria Cherenov Radiaion 発生の理論 aran in in v 数値チェレンコフ放射が発生する理論線は arin an v in

16 数値チェレンコフ発生の理論を確認 6 in an arin v 実空間 B フーリエ空間 /

17 無衝突衝撃波における物理過程 PIC 法による解析 電磁場の位相速度が数値分散により落ちる 相対論的な流れを扱う際には Numeria Cherenov Radiaion 発生 PSATD 法フーリエ空間で Mawe 方程式を解き 解析解を得る方法 特定のクーラン数を選択し 不安定性の成長率を抑える方法 クーラン数 エイリアスや 計算コストの問題

18 近年における数値チェレンコフに対する発表 Va e a. J. Comp. Ph 0 Numeria mehod for inabii in he modeing of aer waefied aeeraor in a Loren booed frame 特定のクーラン数において 不安定性の成長率が落ちる Godfre and Va J. Comp. Ph. 03 Xu e a. Comp. Ph Commun. 03 クーラン数が 0.5 において 成長率が落ちると発表された 以上の論文とは異なったスキームである pcans を用いて同様の結果が得られるか確認

19 growh rae.δ/ pcans による不安定性の成長率の変化 磁場エネルギーの時間変化 = 成長率とする CFL Δ/Δ

20 クーラン数.0 において Numeria Cherenov Radiaion を抑えられているか 実空間 / /

21 クーラン数.0 において Numeria Cherenov Radiaion を抑えられているか 0.5 フーリエ空間.0

22 相対論的衝撃波への応用シミュレーション設定 v h 背景磁場なし 上流速度のローレンツ因子 熱速度 n = 0 セルあたりの平均粒子数 電子 陽電子 粒子を反射させる壁

23 相対論的衝撃波への応用 磁場 0.5 B 4nm pe.0 pe B 4nm 数値チェレンコフ放射 pe

24 相対論的衝撃波への応用 温度 T / m 0.5 pe.0 pe pe

25 まとめ 考察 陰的解法による分散関係を求めることにより数値チェレンコフ放射発生の理論を理解 pcansスキームにおいて クーラン数を.0 付近にとることで 数値的な不安定性の成長率が抑えられる 紹介した論文におけるクーラン数 0.5 と異なるのは 数値分散関係の違いによるものと考えられる.0 をとることで 0. 5 と比較して を大きくとれる 計算コスト削減 空間 時間 次精度 クーラン数を 倍にすると誤差の大きさに 4 倍の違い

26 今後の課題 今回得た数値チェレンコフが解消される クーラン数.0 を用いて 多次元の相対論的衝撃波における粒子加速 の解析にアプローチしていく ご清聴ありがとうございました

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

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

More information

Microsoft PowerPoint - hoshino_part1.pptx

Microsoft PowerPoint - hoshino_part1.pptx 天体プラズマにおける粒子加速機構 無衝突系粒子シミュレーション 星野真弘理学系研究科 地惑惑星科学専攻 内容 ** PIC (Particle In Cell) シミュレーション ** プラズマ輸送係数の研究 (MHD では現象論的に扱う粘性 電気抵抗 熱伝導に関わる物理 ) 非熱的プラズマ ( 局所的熱平衡ではない物理 ) 電子とイオンのエネルギー分配 宇宙での高エネルギー粒子の観測 数値チェレンコフ問題

More information

Microsoft PowerPoint - hiei_MasterThesis

Microsoft PowerPoint - hiei_MasterThesis LHC 加速器での鉛鉛衝突における中性 πおよびω 中間子測定の最適化 日栄綾子 M081043 クォーク物理学研究室 目的 概要 目的 LHC 加速器における TeV 領域の鉛鉛衝突実験における中性 π および ω 中間子の測定の実現可能性の検証 および実際の測定へ向けた最適化 何故鉛鉛衝突を利用して 何を知りたいのか中性 πおよびω 中間子測定の魅力 ALICE 実験検出器群 概要予想される統計量およびバックグランドに対するシグナルの有意性を見積もった

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

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

Microsoft PowerPoint - ILEsympo_matsukiyo.pptx

Microsoft PowerPoint - ILEsympo_matsukiyo.pptx プラズマ衝撃波のマルチスケール構造の精密測定 九州大学松清修一 --- 共同研究者 --- 坂和洋一 佐野孝好 (LE) 富田健太郎 森田太智 ( 九大 ) 蔵満康浩 ( 国立中央大 ) 山崎了 ( 青学大 ) 無衝突衝撃波研究 : その意義と現状 レーザー実験のポテンシャル カギを握る局所量計測 TS 計測の有用性 V-CTS システムの構築 今後の展望 まとめ 1 無衝突衝撃波研究 : その意義と現状

More information

PowerPoint Presentation

PowerPoint Presentation 06 年 8 月 日 ( 月 )-6 日 ( 金 ) 千葉大学総合校舎 号館 4 階情報演習室 宇宙磁気流体 プラズマシミュレーションサマースクール 差分法の基礎 三好隆博 広島大学大学院理学研究科 時限目の目標 線形移流方程式 コンピュータ を計算機で解く! 内容 はじめに 差分法 移流方程式の差分法 高次精度風上差分法 はじめに はじめに 微分方程式 未知関数とその導関数を含む方程式 自然現象などを記述する基礎方程式

More information

Microsoft PowerPoint - qchem3-9

Microsoft PowerPoint - qchem3-9 008 年度冬学期 量子化学 Ⅲ 章量子化学の応用 4.4. 相対論的効果 009 年 月 8 日 担当 : 常田貴夫准教授 相対性理論 A. Einstein 特殊相対論 (905 年 ) 相対性原理: ローレンツ変換に対して物理法則の形は不変 光速度不変 : 互いに等速運動する座標系で光速度は常に一定 ミンコフスキーの4 次元空間座標系 ( 等速系のみ ) 一般相対論 (96 年 ) 等価原理

More information

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

計算機シミュレーション

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

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

<4D F736F F F696E74202D208D E9197BF288CF68A4A B8CDD8AB B83685D>

<4D F736F F F696E74202D208D E9197BF288CF68A4A B8CDD8AB B83685D> 離散化手法とスキームの基礎 と選択法 007//6 宇宙航空研究開発機構情報 計算工学センター嶋英志 本講習の目的 基礎的な計算法の性質を述べ 各手法の持つ長所短所を理解することによって 手法の背景を理解した正しい選択に近づくこと クーラン数 風上差分 等の広い範囲の CFD 技術に共通の概念について その意味とイメージを把握すること 本講習の方針 様々な流体方程式の基礎となる移流方程式を用いて色々な計算法の特徴を計算例を示しながら解説する

More information

PRESS RELEASE 2019 年 4 月 3 日理化学研究所金沢大学国立天文台 ガンマ線バーストのスペクトルと明るさの相関関係の起源 - 宇宙最大の爆発現象の理論的解明へ前進 - 理化学研究所 ( 理研 ) 開拓研究本部長瀧天体ビッグバン研究室の伊藤裕貴研究員 長瀧重博主任研究員 数理創造プ

PRESS RELEASE 2019 年 4 月 3 日理化学研究所金沢大学国立天文台 ガンマ線バーストのスペクトルと明るさの相関関係の起源 - 宇宙最大の爆発現象の理論的解明へ前進 - 理化学研究所 ( 理研 ) 開拓研究本部長瀧天体ビッグバン研究室の伊藤裕貴研究員 長瀧重博主任研究員 数理創造プ PRESS RELEASE 2019 年 4 月 3 日理化学研究所金沢大学国立天文台 ガンマ線バーストのスペクトルと明るさの相関関係の起源 - 宇宙最大の爆発現象の理論的解明へ前進 - 理化学研究所 ( 理研 ) 開拓研究本部長瀧天体ビッグバン研究室の伊藤裕貴研究員 長瀧重博主任研究員 数理創造プログラムのドナルド ウォレン研究員 金沢大学理工研究域数物科学系の米德大輔教授らの国際共同研究グループ

More information

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

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

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

Unit 1

Unit 1 Unt 3. プラズマ中の衝突過程 衝突 nutral 原子により遮られる割合 n ndx + d = (1 n ndx) d/dx = n n = xp( n nx) = xp( x / mfp) mfp = 1/(n n) man fr path = mfp / v collson tm = 1/ = n nv collson frquncy ( 電子の速度分布について平均 ) 電離 再結合水素原子を考える

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

コロイド化学と界面化学

コロイド化学と界面化学 環境表面科学講義 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

素粒子物理学2 素粒子物理学序論B 2010年度講義第4回

素粒子物理学2 素粒子物理学序論B 2010年度講義第4回 素粒子物理学 素粒子物理学序論B 010年度講義第4回 レプトン数の保存 崩壊モード 寿命(sec) n e ν 890 崩壊比 100% Λ π.6 x 10-10 64% π + µ+ νµ.6 x 10-8 100% π + e+ νe 同上 1. x 10-4 Le +1 for νe, elμ +1 for νμ, μlτ +1 for ντ, τレプトン数はそれぞれの香りで独立に保存

More information

( 全体 ) 年 1 月 8 日,2017/1/8 戸田昭彦 ( 参考 1G) 温度計の種類 1 次温度計 : 熱力学温度そのものの測定が可能な温度計 どれも熱エネルギー k B T を

( 全体 ) 年 1 月 8 日,2017/1/8 戸田昭彦 ( 参考 1G) 温度計の種類 1 次温度計 : 熱力学温度そのものの測定が可能な温度計 どれも熱エネルギー k B T を ( 全体 htt://home.hiroshima-u.ac.j/atoda/thermodnamics/ 9 年 月 8 日,7//8 戸田昭彦 ( 参考 G 温度計の種類 次温度計 : 熱力学温度そのものの測定が可能な温度計 どれも熱エネルギー k T を単位として決められている 9 年 月 日 ( 世界計量記念日 から, 熱力学温度 T/K の定義も熱エネルギー k T/J に基づく. 定積気体温度計

More information

差分スキーム 物理 化学 生物現象には微分方程式でモデル化される例が多い モデルを使って現実の現象をコンピュータ上で再現することをシミュレーション ( 数値シミュレーション コンピュータシミュレーション ) と呼ぶ そのためには 微分方程式をコンピュータ上で計算できる数値スキームで近似することが必要

差分スキーム 物理 化学 生物現象には微分方程式でモデル化される例が多い モデルを使って現実の現象をコンピュータ上で再現することをシミュレーション ( 数値シミュレーション コンピュータシミュレーション ) と呼ぶ そのためには 微分方程式をコンピュータ上で計算できる数値スキームで近似することが必要 差分スキーム 物理 化学 生物現象には微分方程式でモデル化される例が多い モデルを使って現実の現象をコンピュータ上で再現することをシミュレーション ( 数値シミュレーション コンピュータシミュレーション ) と呼ぶ そのためには 微分方程式をコンピュータ上で計算できる数値スキームで近似することが必要になる その一つの方法が微分方程式を差分方程式におき直すことである 微分方程式の差分化 次の 1 次元境界値問題を考える

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

Microsoft PowerPoint - em01.pptx

Microsoft PowerPoint - em01.pptx No. 基礎 ~ マクスウェルの方程式 ~ t t D H B E d t d d d t d D l H B l E 微分形積分形 電磁気学の知識からマクスウェルの方程式を導く No. ファラデーの法則 V d dt E dl t B d ストークスの定理を使って E d E ファラデー : 近接作用 界の概念を提唱 B t t B d アンペアの法則 I H rh I H dl d r dl V

More information

Microsoft Word - NumericalComputation.docx

Microsoft Word - NumericalComputation.docx 数値計算入門 武尾英哉. 離散数学と数値計算 数学的解法の中には理論計算では求められないものもある. 例えば, 定積分は, まずは積分 ( 被積分関数の原始関数をみつけること できなければ値を得ることはできない. また, ある関数の所定の値における微分値を得るには, まずその関数の微分ができなければならない. さらに代数方程式の解を得るためには, 解析的に代数方程式を解く必要がある. ところが, これらは必ずしも解析的に導けるとは限らない.

More information

解析力学B - 第11回: 正準変換

解析力学B - 第11回: 正準変換 解析力学 B 第 11 回 : 正準変換 神戸大 : 陰山聡 ホームページ ( 第 6 回から今回までの講義ノート ) http://tinyurl.com/kage2010 2011.01.27 正準変換 バネ問題 ( あえて下手に座標をとった ) ハミルトニアンを考える q 正準方程式は H = p2 2m + k 2 (q l 0) 2 q = H p = p m ṗ = H q = k(q

More information

スライド 1

スライド 1 暫定版修正 加筆の可能性あり ( 付録 ) マクスウェルの方程式 : 真空中 () 1. 電磁波 ( 光波 ) の姿 : 真空中. エネルギー密度 3. ポインティング ベクトル 4. 絵解き : ポインティング ベクトル 5. ポインティング ベクトル : 再確認 6. 両者の関係 7. 付録 : ベクトル解析 注意 1. 本付録 : マクスウェルの方程式: 微分型 を使用. マクスウェルの方程式を数学的に取扱います

More information

スライド 1

スライド 1 暫定版修正 加筆の可能性あり ( 付録 ) 屈折率と誘電率 : 金属. 復習. 電気伝導度 3. アンペールの法則の修正 4. 表皮効果 表皮深さ 5. 鏡の反射 6. 整理 : 電子振動子模型 注意 : 整理しましょう! 前回 : 付録 (4) のアプローチ. 屈折率と損失について記述するために分極振動 ( 電気双極子の集団運動 ) による電気双極子放射を考慮. 誘電率は 真空中の値 を採用 オリジナル光

More information

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

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

More information

Microsoft Word - 素粒子物理学I.doc

Microsoft Word - 素粒子物理学I.doc 6. 自発的対称性の破れとヒッグス機構 : 素粒子の標準模型 Dc 方程式.5 を導くラグランジアンは ϕ ϕ mϕϕ 6. である [H] Eu-nn 方程式 を使って 6. のラグランジア ンから Dc 方程式が導かれることを示せ 6. ゲージ対称性 6.. U 対称性 :QED ディラック粒子の複素場 ψに対する位相変換 ϕ ϕ 6. に対して ラグランジアンが不変であることを要請する これは簡単に示せる

More information

Microsoft Word - MHD-wave.doc

Microsoft Word - MHD-wave.doc 電磁流体力学波 基礎方程式と線形化 谷一郎著 流れ学 を参考にして以下説明する. 流体の運動と電磁場の作用が相互に干渉する 結果, 磁場が存在する場合の導電性気体中の小さい撹乱は, 普通の音波や電磁波とは異なる波の 形で伝播する. 特に伝播速度は, 伝播の方向が磁場の方向となす角に関係する. このような伝播 の法則を明らかにすることは, 導電性気体中の物体の運動を正しく理解するためにも重要なもの である.

More information

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

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

More information

ÿþŸb8bn0irt

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

More information

Microsoft Word - Chap17

Microsoft Word - Chap17 第 7 章化学反応に対する磁場効果における三重項機構 その 7.. 節の訂正 年 7 月 日. 節 章の9ページ の赤枠に記載した説明は間違いであった事に気付いた 以下に訂正する しかし.. 式は 結果的には正しいので安心して下さい 磁場 の存在下でのT 状態のハミルトニアン は ゼーマン項 と時間に依存するスピン-スピン相互作用の項 との和となる..=7.. g S = g S z = S z g

More information

Microsoft PowerPoint - 発表II-3原稿r02.ppt [互換モード]

Microsoft PowerPoint - 発表II-3原稿r02.ppt [互換モード] 地震時の原子力発電所燃料プールからの溢水量解析プログラム 地球工学研究所田中伸和豊田幸宏 Central Research Institute of Electric Power Industry 1 1. はじめに ( その 1) 2003 年十勝沖地震では 震源から離れた苫小牧地区の石油タンクに スロッシング ( 液面揺動 ) による火災被害が生じた 2007 年中越沖地震では 原子力発電所内の燃料プールからの溢水があり

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

Microsoft PowerPoint _量子力学短大.pptx

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

More information

天体物理特論

天体物理特論 高エネルギー宇宙ニュートリノ : 突発天体起源の可能性について 浅野勝晃 ( 東工大 ) IceCube による PeV ニュートリノ検出 2 イベント 7.8x10 5-5.6x10 6 GeV 8.9x10 5-8.5x10 6 GeV 当初は最高エネルギー宇宙線起源の 10 18 ev(eev) 程度のニュートリノ検出が期待されていた 予想を裏切って 10 15 ev(pev) のニュートリノが最初に検出された!

More information

untitled

untitled に, 月次モデルの場合でも四半期モデルの場合でも, シミュレーション期間とは無関係に一様に RMSPE を最小にするバンドの設定法は存在しないということである 第 2 は, 表で与えた 2 つの期間及びすべての内生変数を見渡して, 全般的にパフォーマンスのよいバンドの設定法は, 最適固定バンドと最適可変バンドのうちの M 2, Q2 である いずれにしても, 以上述べた 3 つのバンド設定法は若干便宜的なものと言わざるを得ない

More information

横浜市環境科学研究所

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

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

Microsoft PowerPoint - 第8章

Microsoft PowerPoint - 第8章 講義予定 案. 9/ 数値シミュレーションの手続き テキスト第 章. 9/ 9 偏微分方程式と解析解 テキスト第 章 3. 9/6 休講 4. 9/30 差分方程式とそのスキーム テキスト第 3 章 変換 テキスト第 4 章 5. 0/ 7 計算 テキスト第 5 章 連立一次方程式の解法 テキスト第 6 章 6. 0/ 流れ関数 ポテンシャルによる解法 テキスト第 7 章 7. 0/8 流速 圧力を用いた解法

More information

Problem P5

Problem P5 問題 P5 メンシュトキン反応 三級アミンとハロゲン化アルキルの間の求核置換反応はメンシュトキン反応として知られている この実験では DABCO(1,4 ジアザビシクロ [2.2.2] オクタン というアミンと臭化ベンジルの間の反応速度式を調べる N N Ph Br N N Br DABCO Ph DABCO 分子に含まれるもう片方の窒素も さらに他の臭化ベンジルと反応する可能性がある しかし この実験では

More information

社会保険料の賃金への影響について

社会保険料の賃金への影響について 社会保険料の賃金への影響について Borja,G. Labor economic, 3r e McGraw-Hill, Chapter, -3: Policy Application: payroll taxe an ubiie N グレゴリー マンキュー マンキュー経済学 Ⅰミクロ編 足立他訳 東洋経済新報社 2000 年 68-78 ページただし 保険料 ( 税金 ) のかかり方は 教科書のものと以下で扱うものとでは異なっていることに注意.

More information

Microsoft Word - プレス原稿_0528【最終版】

Microsoft Word - プレス原稿_0528【最終版】 報道関係各位 2014 年 5 月 28 日 二酸化チタン表面における陽電子消滅誘起イオン脱離の観測に成功 ~ 陽電子を用いた固体最表面の改質に道 ~ 東京理科大学研究戦略 産学連携センター立教大学リサーチ イニシアティブセンター 本研究成果のポイント 二酸化チタン表面での陽電子の対消滅に伴って脱離する酸素正イオンの観測に成功 陽電子を用いた固体最表面の改質に道を拓いた 本研究は 東京理科大学理学部第二部物理学科長嶋泰之教授

More information

銀河風の定常解

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

More information

スライド 1

スライド 1 暫定版修正 加筆の可能性あり ( 付録 ) 電磁波 ( 光 ) の角運動量. 復習 : 電磁波 ( 光 ) のエネルギー. 運動量 角運動量 ( 実空間 ) 3. 軌道 スピン角運動量 4. 円偏光状態 5. 螺旋状態 付録 8 のアプローチ. 本付録では電磁波 ( 光 ) の軌道 スピン角運動量ついて古典的に扱う. スピン角運動量は直線偏光状態では零 円偏光状態では非零 右 左回りで大きさは同じ

More information

今後の予定 6/29 パターン形成第 11 回 7/6 データ解析第 12 回 7/13 群れ行動 ( 久保先生 ) 第 13 回 7/17 ( 金 ) 休講 7/20 まとめ第 14 回 7/27 休講?

今後の予定 6/29 パターン形成第 11 回 7/6 データ解析第 12 回 7/13 群れ行動 ( 久保先生 ) 第 13 回 7/17 ( 金 ) 休講 7/20 まとめ第 14 回 7/27 休講? 今後の予定 6/29 パターン形成第 11 回 7/6 データ解析第 12 回 7/13 群れ行動 ( 久保先生 ) 第 13 回 7/17 ( 金 ) 休講 7/20 まとめ第 14 回 7/27 休講? 数理生物学演習 第 11 回パターン形成 本日の目標 2 次元配列 分子の拡散 反応拡散モデル チューリングパタン 拡散方程式 拡散方程式 u t = D 2 u 拡散が生じる分子などの挙動を記述する.

More information

(3) E-I 特性の傾きが出力コンダクタンス である 添え字 は utput( 出力 ) を意味する (4) E-BE 特性の傾きが電圧帰還率 r である 添え字 r は rrs( 逆 ) を表す 定数の値は, トランジスタの種類によって異なるばかりでなく, 同一のトランジスタでも,I, E, 周

(3) E-I 特性の傾きが出力コンダクタンス である 添え字 は utput( 出力 ) を意味する (4) E-BE 特性の傾きが電圧帰還率 r である 添え字 r は rrs( 逆 ) を表す 定数の値は, トランジスタの種類によって異なるばかりでなく, 同一のトランジスタでも,I, E, 周 トランジスタ増幅回路設計入門 pyrgt y Km Ksaka 005..06. 等価回路についてトランジスタの動作は図 のように非線形なので, その動作を簡単な数式で表すことができない しかし, アナログ信号を扱う回路では, 特性グラフのの直線部分に動作点を置くので線形のパラメータにより, その動作を簡単な数式 ( 一次式 ) で表すことができる 図. パラメータトランジスタの各静特性の直線部分の傾きを数値として特性を表したものが

More information

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

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

More information

スライド 1

スライド 1 グループ発表天体核研究室 低光度ガンマ線バーストの起源 D2 当真賢二 宇宙ひもを重力レンズで探る D3 須山輝明 2006 年度物理学第二教室教室発表会 @ 第四講義室 天体核研究室の大雑把な研究グループ 天体物理学中村 犬塚 井岡 山田 PD: 町田 石津 三浦 D3: 道越 宇宙論中村 田中 早田 D3: 須山 D2: 横山 D1: 泉 M2: 棚橋 村田 D2: 井上 ( 剛 ) 当真 D1:

More information

JPS2012spring

JPS2012spring BelleII 実験用 TOP カウンターの性能評価 2012.7.7( 土 ) 名古屋大学高エネルギー物理学研究室 (N 研究室 ) 有田義宣 BelleII に搭載する粒子識別装置 TOP カウンター 2 BelleII 実験 もっとも識別の難しい π/k 識別 BelleⅡ 実験は Belle 実験をさらに高輝度化 (40 倍 ) し 大量の B 中間子からの稀崩壊現象を探る電子陽電子コライダー

More information

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

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

More information

<4D F736F F F696E74202D208CB48E7197CD8A7789EF F4882CC91E589EF8AE989E A2E B8CDD8AB B83685D>

<4D F736F F F696E74202D208CB48E7197CD8A7789EF F4882CC91E589EF8AE989E A2E B8CDD8AB B83685D> 数値解析技術と標準 (3) 数値解析の信頼性に関する標準 平成 24 年 9 月 21 日原子力学会 2012 秋の大会標準委員会セッション5( 基盤 応用専門部会 ) 独立行政法人原子力安全基盤機構原子力システム安全部堀田亮年 AESJ MTG 2012 Autumn @Hiroshima 1 シミュレーションの信頼性 WG 報告書の構成 本文 (118 頁 ):V&Vの構造案解説 A) V&V

More information

1/17 平成 29 年 3 月 25 日 ( 土 ) 午前 11 時 1 分量子力学とクライン ゴルドン方程式 ( 学部 3 年次秋学期向 ) 量子力学とクライン ゴルドン方程式 素粒子の満たす場 y ( x,t) の運動方程式 : クライン ゴルドン方程式 : æ 3 ö ç å è m= 0

1/17 平成 29 年 3 月 25 日 ( 土 ) 午前 11 時 1 分量子力学とクライン ゴルドン方程式 ( 学部 3 年次秋学期向 ) 量子力学とクライン ゴルドン方程式 素粒子の満たす場 y ( x,t) の運動方程式 : クライン ゴルドン方程式 : æ 3 ö ç å è m= 0 /7 平成 9 年 月 5 日 ( 土 午前 時 分量子力学とクライン ゴルドン方程式 ( 学部 年次秋学期向 量子力学とクライン ゴルドン方程式 素粒子の満たす場 (,t の運動方程式 : クライン ゴルドン方程式 : æ ö ç å è = 0 c + ( t =, 0 (. = 0 ì æ = = = ö æ ö æ ö ç ì =,,,,,,, ç 0 = ç Ñ 0 = ç Ñ 0 Ñ Ñ

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 - 複素数.pptx

Microsoft PowerPoint - 複素数.pptx 00 年 月 9 日 ( 金 第 時限 平成 年度物質科学解析第 7 回 複素数 冨田知志 0. なぜ複素数か?. 虚数単位. 複素数の計算. オイラーの公式. 複素平面 5. 級数での複素数 ( オイラーの公式 の活用 6. 量子力学で出てくる複素数の例 0. なぜ複素数か? 量子論 ( 量子力学 で不可欠だから参照 : 光ナノサイエンスコアI 古典論や電気回路でも複素数は使うただしそれはあくまでも数学的道具

More information

偏微分方程式、連立1次方程式、乱数

偏微分方程式、連立1次方程式、乱数 数値計算法 011/6/8 林田清 大阪大学大学院理学研究科 常微分方程式の応用例 1 Rutherford 散乱 ( 原子核同士の散乱 ; 金の薄膜に α 粒子をあてる ) 1 クーロン力 f= 4 0 r r r Ze y からf cos, si f f f y f f 粒子の 方向 y方向の速度と座標について dv Ze dvy Ze y, 3 3 dt 40m r dt 40m r d dy

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 復習 ) 時系列のモデリング ~a. 離散時間モデル ~ y k + a 1 z 1 y k + + a na z n ay k = b 0 u k + b 1 z 1 u k + + b nb z n bu k y k = G z 1 u k = B(z 1 ) A(z 1 u k ) ARMA モデル A z 1 B z 1 = 1 + a 1 z 1 + + a na z n a = b 0

More information

ハートレー近似(Hartree aproximation)

ハートレー近似(Hartree aproximation) ハートリー近似 ( 量子多体系の平均場近似 1) 0. ハミルトニアンの期待値の変分がシュレディンガー方程式と等価であること 1. 独立粒子近似という考え方. 電子系におけるハートリー近似 3.3 電子系におけるハートリー近似 Mde by R. Okmoto (Kyushu Institute of Technology) filenme=rtree080609.ppt (0) ハミルトニアンの期待値の変分と

More information

EP7000取扱説明書

EP7000取扱説明書 EP7000 S0109-3012 3 47 811 1213 1419 2021 53 54 5560 61 6263 66 2223 2427 2830 3133 3436 3740 4142 4344 45 46 4750 5152 2 4 5 6 7 1 3 4 5 6 7 8 9 15 16 17 18 13 EP7000 2 10 11 12 13 14 19 20 21 22 23 24

More information

Microsoft PowerPoint SIGAL.ppt

Microsoft PowerPoint SIGAL.ppt アメリカン アジアンオプションの 価格の近似に対する 計算幾何的アプローチ 渋谷彰信, 塩浦昭義, 徳山豪 ( 東北大学大学院情報科学研究科 ) 発表の概要 アメリカン アジアンオプション金融派生商品の一つ価格付け ( 価格の計算 ) は重要な問題 二項モデルにおける価格付けは計算困難な問題 目的 : 近似精度保証をもつ近似アルゴリズムの提案 アイディア : 区分線形関数を計算幾何手法により近似 問題の説明

More information

2014 年 10 月 2 日 本日の講義及び演習 数値シミュレーション 2014 年度第 2 回 偏微分方程式の偏微分項をコンピュータで扱えるようにする 離散化 ( 差分化 ) テイラー展開の利用 1 階微分項に対する差分式 2 階微分項に対する差分式 1 次元熱伝導方程式に適用して差分式を導出

2014 年 10 月 2 日 本日の講義及び演習 数値シミュレーション 2014 年度第 2 回 偏微分方程式の偏微分項をコンピュータで扱えるようにする 離散化 ( 差分化 ) テイラー展開の利用 1 階微分項に対する差分式 2 階微分項に対する差分式 1 次元熱伝導方程式に適用して差分式を導出 04 年 0 月 日 本日の講義及び演習 数値シミュレーション 04 年度第 回 偏微分方程式の偏微分項をコンピュータで扱えるようにする 離散化 ( 差分化 テイラー展開の利用 階微分項に対する差分式 階微分項に対する差分式 次元熱伝導方程式に適用して差分式を導出 Ecel を利用した温度変化シミュレーション 永野 ( 熱流体システム研究室 hagao@tc.ac.p 重要! 熱の伝わり方 ( 伝熱モード

More information

ニュートン重力理論.pptx

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

More information

画像類似度測定の初歩的な手法の検証

画像類似度測定の初歩的な手法の検証 画像類似度測定の初歩的な手法の検証 島根大学総合理工学部数理 情報システム学科 計算機科学講座田中研究室 S539 森瀧昌志 1 目次 第 1 章序論第 章画像間類似度測定の初歩的な手法について.1 A. 画素値の平均を用いる手法.. 画素値のヒストグラムを用いる手法.3 C. 相関係数を用いる手法.4 D. 解像度を合わせる手法.5 E. 振れ幅のヒストグラムを用いる手法.6 F. 周波数ごとの振れ幅を比較する手法第

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

OpenFOAM(R) ソースコード入門 pt1 熱伝導方程式の解法から有限体積法の実装について考える 前編 : 有限体積法の基礎確認 2013/11/17 オープンCAE 富山富山県立大学中川慎二

OpenFOAM(R) ソースコード入門 pt1 熱伝導方程式の解法から有限体積法の実装について考える 前編 : 有限体積法の基礎確認 2013/11/17 オープンCAE 富山富山県立大学中川慎二 OpenFOAM(R) ソースコード入門 pt1 熱伝導方程式の解法から有限体積法の実装について考える 前編 : 有限体積法の基礎確認 2013/11/17 オープンCAE 勉強会 @ 富山富山県立大学中川慎二 * OpenFOAM のソースコードでは, 基礎式を偏微分方程式の形で記述する.OpenFOAM 内部では, 有限体積法を使ってこの微分方程式を解いている. どのようにして, 有限体積法に基づく離散化が実現されているのか,

More information

Microsoft PowerPoint pptx

Microsoft PowerPoint pptx 4.2 小信号パラメータ 1 電圧利得をどのように求めるか 電圧ー電流変換 入力信号の変化 dv BE I I e 1 v be の振幅から i b を求めるのは難しい? 電流増幅 電流ー電圧変換 di B di C h FE 電流と電圧の関係が指数関数になっているのが問題 (-RC), ただし RL がない場合 dv CE 出力信号の変化 2 pn 接合の非線形性への対処 I B 直流バイアスに対する抵抗

More information

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

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

More information

宇宙線の起源と加速と伝搬 大平豊青山学院大学 内容 銀河宇宙線の標準モデル 標準モデルと観測の矛盾 まとめ

宇宙線の起源と加速と伝搬 大平豊青山学院大学 内容 銀河宇宙線の標準モデル 標準モデルと観測の矛盾 まとめ 宇宙線の起源と加速と伝搬 大平豊青山学院大学 内容 銀河宇宙線の標準モデル 標準モデルと観測の矛盾 まとめ 宇宙線のエネルギースペクトル dn/de E -2.7 (E

More information

スライド 1

スライド 1 暫定版修正 加筆の可能性あり ( 付録 ) 準備 : 非線形光学効果 (). 絵解き : 第二高調波発生. 基本波の波動方程式 3. 第二高調波の波動方程式 4. 二倍分極振動 : ブランコ 5. 結合波動方程式へ 6. 補足 : 非線形電気感受率 ( 複素数 ) 付録 43 のアプローチ. 分極振動とは振動電場に誘われて伸縮する電気双極子の集団運動. 電気感受率と波動方程式の関係を明らかにする 3.

More information

目次 Boltzmann 方程式 自己重力系 Vlasov-Poisson 方程式系 6次元位相空間上でのVlasov-Poissonシミュレーション 自己重力系でのVlasovシミュレーションの応用例 まとめ

目次 Boltzmann 方程式 自己重力系 Vlasov-Poisson 方程式系 6次元位相空間上でのVlasov-Poissonシミュレーション 自己重力系でのVlasovシミュレーションの応用例 まとめ 無衝突自己重力系の 無衝突ボルツマンシミュレーション 宇宙磁気流体 プラズマシミュレーションワークショップ 2013年2月18日 千葉大学 吉川 耕司 (筑波大学 計算科学研究センター) 吉田 直紀 (東京大学 物理 Kavli IPMU) 梅村 雅之 (筑波大学 計算科学研究センター) Direct Integration of the Collisionless Boltzmann Equation

More information

<4D F736F F F696E74202D20906C8D488AC28BAB90DD8C7689F090CD8D488A D91E F1>

<4D F736F F F696E74202D20906C8D488AC28BAB90DD8C7689F090CD8D488A D91E F1> 人工環境設計解析工学構造力学と有限要素法 ( 第 回 ) 東京大学新領域創成科学研究科 鈴木克幸 固体力学の基礎方程式 変位 - ひずみの関係 適合条件式 ひずみ - 応力の関係 構成方程式 応力 - 外力の関係 平衡方程式 境界条件 変位規定境界 反力規定境界 境界条件 荷重応力ひずみ変形 場の方程式 Γ t Γ t 平衡方程式構成方程式適合条件式 構造力学の基礎式 ひずみ 一軸 荷重応力ひずみ変形

More information

ディジタル信号処理

ディジタル信号処理 ディジタルフィルタの設計法. 逆フィルター. 直線位相 FIR フィルタの設計. 窓関数法による FIR フィルタの設計.5 時間領域での FIR フィルタの設計 3. アナログフィルタを基にしたディジタル IIR フィルタの設計法 I 4. アナログフィルタを基にしたディジタル IIR フィルタの設計法 II 5. 双 次フィルタ LI 離散時間システムの基礎式の証明 [ ] 4. ] [ ]*

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション モニタリング観測からわかった電波銀河 3C111 の γ 線活動期と電波ノットの噴出時期との関係 VLBI 懇談会シンポジウム 12 月 27 日 ( 火 ) 山口大学 B4 塩谷康允共同研究者 : 藤澤健太 新沼浩太郎 導入 AGN 統一モデル AGN 電波で明るい (10 %) 超大質量 BH+ 降着円盤 電波で暗い (90 %) 莫大なエネルギー放射 (10 6-14 L ) 0 いくつかの種類に大別される

More information

画像解析論(2) 講義内容

画像解析論(2) 講義内容 画像解析論 画像解析論 東京工業大学長橋宏 主な講義内容 信号処理と画像処理 二次元システムとその表現 二次元システムの特性解析 各種の画像フィルタ 信号処理と画像処理 画像解析論 処理の応答 記憶域 入出力の流れ 信号処理系 実時間性が求められる メモリ容量に対する制限が厳しい オンラインでの対応が厳しく求められる 画像処理系 ある程度の処理時間が許容される 大容量のメモリ使用が容認され易い オフラインでの対応が容認され易い

More information

目次 2 1. イントロダクション 2. 実験原理 3. データ取得 4. データ解析 5. 結果 考察 まとめ

目次 2 1. イントロダクション 2. 実験原理 3. データ取得 4. データ解析 5. 結果 考察 まとめ オルソポジトロニウムの寿命測定による QED の実験的検証 課題演習 A2 2016 年後期 大田力也鯉渕駿龍澤誠之 羽田野真友喜松尾一輝三野裕哉 目次 2 1. イントロダクション 2. 実験原理 3. データ取得 4. データ解析 5. 結果 考察 まとめ 第 1 章イントロダクション 実験の目的 4 ポジトロニウム ( 後述 ) の崩壊を観測 オルソポジトロニウム ( スピン 1 状態 ) の寿命を測定

More information

Chapter 版 Maxima を用いた LC のインピーダンス測定について [ 目的 ] 電気通信大学 先進理工学科の2 年次後期に実施される電気 電子回路実験において L,C のインピーダンス測定を実施している この実験項目について 無料ソフトの Maxima を用い

Chapter 版 Maxima を用いた LC のインピーダンス測定について [ 目的 ] 電気通信大学 先進理工学科の2 年次後期に実施される電気 電子回路実験において L,C のインピーダンス測定を実施している この実験項目について 無料ソフトの Maxima を用い Chapter 2 2016.10.14 版 Maxima を用いた LC のインピーダンス測定について [ 目的 ] 電気通信大学 先進理工学科の2 年次後期に実施される電気 電子回路実験において L,C のインピーダンス測定を実施している この実験項目について 無料ソフトの Maxima を用いることで 理論解析と実験値の比較が可能である また 近年のパソコンの性能の向上により Maxima の実行処理速度が大幅に改善された

More information

三重大学工学部

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

More information

B. モル濃度 速度定数と化学反応の速さ 1.1 段階反応 ( 単純反応 ): + I HI を例に H ヨウ化水素 HI が生成する速さ は,H と I のモル濃度をそれぞれ [ ], [ I ] [ H ] [ I ] に比例することが, 実験により, わかっている したがって, 比例定数を k

B. モル濃度 速度定数と化学反応の速さ 1.1 段階反応 ( 単純反応 ): + I HI を例に H ヨウ化水素 HI が生成する速さ は,H と I のモル濃度をそれぞれ [ ], [ I ] [ H ] [ I ] に比例することが, 実験により, わかっている したがって, 比例定数を k 反応速度 触媒 速度定数 反応次数について. 化学反応の速さの表し方 速さとは単位時間あたりの変化の大きさである 大きさの値は 0 以上ですから, 速さは 0 以上の値をとる 化学反応の速さは単位時間あたりの物質のモル濃度変化の大きさで表すのが一般的 たとえば, a + bb c (, B, は物質, a, b, c は係数 ) という反応において,, B, それぞれの反応の速さを, B, とし,

More information

第6章 実験モード解析

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

More information

3. 重力波と沿岸 赤道ケルビン波見延庄士郎 ( 海洋気候物理学研究室 ) 予習課題 を印刷解答の上, 授業冒頭で提出してください. 予習のための課題なので

3. 重力波と沿岸 赤道ケルビン波見延庄士郎 ( 海洋気候物理学研究室 ) 予習課題   を印刷解答の上, 授業冒頭で提出してください. 予習のための課題なので 3. 重力波と沿岸 赤道ケルビン波見延庄士郎 ( 海洋気候物理学研究室 ) 予習課題 hp://www.sci.hokudai.ac.jp/~minobe/class/po_clm/03_e_propagaion.pd を印刷解答の上, 授業冒頭で提出してください. 予習のための課題なので, 冒頭のみ提出を受け付けます. また以下の you ube のビデオを見ておくこと. 個々のビデオは全部は見ずに,

More information

2018/6/12 表面の電子状態 表面に局在する電子状態 表面電子状態表面準位 1. ショックレー状態 ( 準位 ) 2. タム状態 ( 準位 ) 3. 鏡像状態 ( 準位 ) 4. 表面バンドのナローイング 5. 吸着子の状態密度 鏡像力によるポテンシャル 表面からzの位置の電子に働く力とポテン

2018/6/12 表面の電子状態 表面に局在する電子状態 表面電子状態表面準位 1. ショックレー状態 ( 準位 ) 2. タム状態 ( 準位 ) 3. 鏡像状態 ( 準位 ) 4. 表面バンドのナローイング 5. 吸着子の状態密度 鏡像力によるポテンシャル 表面からzの位置の電子に働く力とポテン 表面の電子状態 表面に局在する電子状態 表面電子状態表面準位. ショックレー状態 ( 準位. タム状態 ( 準位 3. 鏡像状態 ( 準位 4. 表面バンドのナローイング 5. 吸着子の状態密度 鏡像力によるポテンシャル 表面からzの位置の電子に働く力とポテンシャル e F z ( z z e V ( z ( Fz dz 4z e V ( z 4z ( z > ( z < のときの電子の運動を考える

More information

Microsoft Word - planck定数.doc

Microsoft Word - planck定数.doc . 目的 Plck 定数 光電効果についての理解を深める. また光電管を使い実際に光電効果を観察し,Plck 定数および仕事関数を求める.. 課題 Hg- スペクトルランプから出ている何本かの強いスペクトル線のなかから, フィルターを使い, 特定の波長域のスペクトル線を選択し, それぞれの場合について光電効果により飛び出してくる電子の最高エネルギーを測定する. この測定結果から,Plck 定数 h

More information

...Y..FEM.pm5

...Y..FEM.pm5 . 剛塑性有限要素法 名古屋大学大学院工学研究科. はじめに. 剛塑性体の構成式.. 降伏条件.. 構成方程式 ([D] マトリックス ). 節点速度 ひずみ速度関係..[B] マトリックス.. 四角形一次要素の [B] マトリックス.4 4 仮想仕事の原理 ( 剛性マトリックス ([K] マトリックス )).5 非線形方程式の解法.5. 直接代入法.5.wto-Raphso 法.6 非圧縮性の拘束と数値積分.7

More information

Microsoft Word - note02.doc

Microsoft Word - note02.doc 年度 物理化学 Ⅱ 講義ノート. 二原子分子の振動. 調和振動子近似 モデル 分子 = 理想的なバネでつながった原子 r : 核間距離, r e : 平衡核間距離, : 変位 ( = r r e ), k f : 力の定数ポテンシャルエネルギー ( ) k V = f (.) 古典運動方程式 [ 振動数 ] 3.3 d kf (.) dt μ : 換算質量 (m, m : 原子, の質量 ) mm

More information

Microsoft PowerPoint - koidenao10.ppt

Microsoft PowerPoint - koidenao10.ppt 抵抗性相対論的 MHD の因果律について < 一般化相対論的 MHD 事始め > 熊本大学理学科 小出眞路 ブラックホール周辺で磁気リコネクションを伴うプラズマの研究で用いられる標準的な抵抗性相対論的 MHD(RMHD) の因果律の問題を提起し, 一般化 RMHDの因果律について議論する 時間があれば, 一般化 GRMHD 方程式の導出とその示唆する現象について紹介する プラズマセミナー 010.06.10(

More information

<4D F736F F F696E74202D F F8F7482CC944E89EF8AE989E6835A E6F325F8CF68A4A94C55231>

<4D F736F F F696E74202D F F8F7482CC944E89EF8AE989E6835A E6F325F8CF68A4A94C55231> 日本原子力学会 2010 年春の年会茨城大学計算科学技術部会企画セッション シミュレーションの信頼性確保の あり方とは? (2) 海外における熱流動解析の信頼性評価の取り組み 平成 22 年 3 月 28 日東芝中田耕太郎 JNES 笠原文雄 調査対象 OECD/NEA CFD ガイドライン NEA/CSNI/R(2007)5 単相 CFD の使用に関する体系的なベストプラクティスガイドライン 原子炉安全解析に対する単相

More information

Microsoft PowerPoint - 10.pptx

Microsoft PowerPoint - 10.pptx m u. 固有値とその応用 8/7/( 水 ). 固有値とその応用 固有値と固有ベクトル 行列による写像から固有ベクトルへ m m 行列 によって線形写像 f : R R が表せることを見てきた ここでは 次元平面の行列による写像を調べる とし 写像 f : を考える R R まず 単位ベクトルの像 u y y f : R R u u, u この事から 線形写像の性質を用いると 次の格子上の点全ての写像先が求まる

More information

21世紀型パラメータ設計―標準SN比の活用―

21世紀型パラメータ設計―標準SN比の活用― 世紀のパラメータ設計ースイッチ機構のモデル化ー 接点 ゴム 変位 スイッチ動作前 スイッチ動作後 反転ばねでスイッチの クリック感 を実現した構造 世紀型パラメータ設計 標準 SN 比の活用 0 世紀の品質工学においては,SN 比の中に, 信号因子の乱れである 次誤差 (S res ) もノイズの効果の中に加えて評価してきた.のパラメータ設計の例では, 比例関係が理想であるから, 次誤差も誤差の仲間と考えてもよかったが,

More information

本講演の内容 1. はじめに 2. 天体 MHD 現象のおもしろさ 3. 天体 MHD シミュレーションは超困難 4. 天体 MHD シミュレーションの魅力 5. 天体 MHD シミュレーションの魔力 はまると危険 ( 落とし穴の数々 ) 6. むすび : ノーベル賞課題 ( 超難問 ) に挑戦せよ

本講演の内容 1. はじめに 2. 天体 MHD 現象のおもしろさ 3. 天体 MHD シミュレーションは超困難 4. 天体 MHD シミュレーションの魅力 5. 天体 MHD シミュレーションの魔力 はまると危険 ( 落とし穴の数々 ) 6. むすび : ノーベル賞課題 ( 超難問 ) に挑戦せよ 2000 年天文天体物理夏の学校における講演を元に少し改訂 天体 MHD 数値シミュレーション の魅力と魔力 柴田一成京大理花山天文台 本講演の内容 1. はじめに 2. 天体 MHD 現象のおもしろさ 3. 天体 MHD シミュレーションは超困難 4. 天体 MHD シミュレーションの魅力 5. 天体 MHD シミュレーションの魔力 はまると危険 ( 落とし穴の数々 ) 6. むすび : ノーベル賞課題

More information

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

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

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