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

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

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

スライド 1

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

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

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

図 宇宙論解析の流れ 次元データの CMB の例 宇宙論ゆらぎ場 F (θ) の測定 左上図 ゆらぎ場のフーリエ波数分解 右上図 右下図は パ ワースペクトル推定の結果 灰色点は各波数ビンでの測定値 エラーバーを伴う青点は 複数の波数ビンで測定値を平均した結果 エラーバーとして 有限数のフーリエモー

FEATURE Kavli IPMU CMB IPMU CMB CMB?? F(θ) CMB F(θ) [T(θ) T ]/T T(θ) θ CMB T CMB 3 F(θ) F(θ) CMB Planck 5 3 40 000 CMB 500 F(θ) 3 3 /60 40 Kavli IPMU News No 6 June 04 図 宇宙論解析の流れ 次元データの CMB の例 宇宙論ゆらぎ場

観測的宇宙論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

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] のような波動を与えたとき どのように時間変化するか調べる

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

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

, 0707

始原的ガス雲の non-biased カタログ : 始原星の初期質量関数 平野信吾 1 細川隆史 1 吉田直紀 1,2 千秋元 1 梅田秀之 1 et al 1 東京大学 2 Kavli IPMU 初代星 初代銀河研究会 2014@ 鹿児島大学 (2014/01/22-24) 始原星の質量 : 星形成過程 始原星 ( 種族 III の星 ; zero-metallicity star) 宇宙の初期進化を左右

プランクの公式と量子化

Planck の公式と量子化 埼玉大学理学部物理学科 久保宗弘 序論 一般に 量子力学 と表現すると Schrödinger の量子力学などの 後期量子力学 を指すことが多い 本当の量子概念 には どうアプローチ? 何故 エネルギーが量子化されるか という根本的な問いにどうこたえるか? どのように 量子 の扉は叩かれたのか? 序論 統計力学 熱力学 がことの始まり 総括的な動き を表現するための学問である

Autodesk Inventor Skill Builders Autodesk Inventor 2010 構造解析の精度改良 メッシュリファインメントによる収束計算 予想作業時間:15 分 対象のバージョン:Inventor 2010 もしくはそれ以降のバージョン シミュレーションを設定する際

Autodesk Inventor Skill Builders Autodesk Inventor 2010 構造解析の精度改良 メッシュリファインメントによる収束計算 予想作業時間:15 分 対象のバージョン:Inventor 2010 もしくはそれ以降のバージョン シミュレーションを設定する際に 収束判定に関するデフォルトの設定をそのまま使うか 修正をします 応力解析ソルバーでは計算の終了を判断するときにこの設定を使います

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

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

PowerPoint Presentation

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

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

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

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

物理計測法特論 No.1 第 1 章 : 信号と雑音 本講義の主題 雑音の性質を理解することで 信号と雑音の大きさが非常に近い状態での信号の測定技術 : 微小信号計測 について学ぶ 講義の Web http://www.g-munu.t.u-tokyo.ac.jp/mio/note/sig_mes/tokuron.html 物理学の基本は実験事実の積み重ねである そして それは何かを測定することから始まる

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

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

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

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

観測的宇宙論workshop.pptx

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

Microsoft Word - note02.doc

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

スライド 1

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

領域シンポ発表

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 項は無視できる この場合の振幅を

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

波動方程式と量子力学 谷村吉隆 京都大学理学研究科化学専攻 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 =

インターリーブＡＤＣでのタイミングスキュー影響のデジタル補正技術

1 インターリーブADCでのタイミングスキュー影響のデジタル補正技術 浅見幸司 黒沢烈士 立岩武徳 宮島広行 小林春夫 ( 株 ) アドバンテスト 群馬大学 2 目次 1. 研究背景 目的 2. インターリーブADCの原理 3. チャネル間ミスマッチの影響 3.1. オフセットミスマッチの影響 3.2. ゲインミスマッチの影響 3.3. タイミングスキューの影響 4. 提案手法 4.1. インターリーブタイミングミスマッチ補正フィルタ

( 慣性抵抗 ) 速度の 2 乗に比例流体中を進む物体は前面にある流体を押しのけて進む. 物 aaa 体の後面には流体が付き従う ( 渦を巻いて ). 前面にある速度 0 の流体が後面に移動して速度 vとなったと考えてよい. この流体の質量は単位時間内に物体が押しのける体積に比例するので,v に比例

空気抵抗があるときの自由落下 抵抗が速度に比例する場合 1. 絵を描く, 座標と情報, 記号を記入する x F0 v

銀河風の定常解

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

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

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

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

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

計算機シミュレーション

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

ニュートン重力理論.pptx

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

DVIOUT-SS_Ma

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

極めて軽いダークマターの 新しい検出方法 In preparation

極めて軽いダークマターの新しい検出方法 In preparation Hajime Fukuda, T.T. Yanagida, S. Matsumoto Kavli IPMU, U. Tokyo August 1, 2017 Introduction DM は最も確立した BSM の一つ 質量は? Particle DM Mass Range dsph m > M Pl Vast Region!

Microsoft PowerPoint - ILEsympo_matsukiyo.pptx

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

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 アインシュタイン表記

航空機の縦系モデルに対する、非線形制御の適用例

制御システム工学研究グルプ 航空機の縦系モデルに対する非線形最適制御の適用例 菊池芳光 * * 名古屋大学 MBD 中部コンファレンス @2014 年 12 月 18 日 目次 はじめに 先行研究 提案手法 縦系航空機モデル シミュレーション結果 おわりに はじめに PIO(Pilot Induced Oscillation) Category II 速度飽和 位相遅れ PIO 事故 PIOにより墜落するGripen

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

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

FEM原理講座　（サンプルテキスト）

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

Microsoft PowerPoint - Ppt ppt[読み取り専用]

Astroparticle physics 富山大学 松本重貴 1. 暗黒物質問題 2. 暗黒物質の正体? 3. 暗黒物質の探査 Astroparticle physics って何? 素粒子 物理学 ニュートリノ暗黒物質暗黒エネルギー宇宙のバリオン数インフレーション 宇宙 物理学 宇宙の暗黒物質問題暗黒物質の存在は確立したが その正体 ( 質量 スピン 量子数や相互作用 ) については不明であるという問題!

Kerr 時空における球対称流に対するコリメーション効果 ( CQG, 26, , 2009 ) 髙見健太郎 ( 広島大学 / Albert-Einstein-Institute) 共同研究者 : 小嶌康史 ( 広島大学 ) 2009 年 10 月 01 日駒場宇宙コロキウム

Kerr 時空における球対称流に対するコリメーション効果 ( CQG, 26, 085013, 2009 ) 髙見健太郎 ( 広島大学 / Albert-Einstein-Institute) 共同研究者 : 小嶌康史 ( 広島大学 ) 2009 年 10 月 01 日駒場宇宙コロキウム 目 次 導入 Kerr 時空と測地線方程式 粒子のコリメーション条件 粒子流に対するコリメーション効果 まとめ

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

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

Formulation and constraints on decaying dark matter with finite mass daughter particles

有限質量をもつ娘粒子へ崩壊する 暗黒物質モデルとその観測的制限 名古屋大学大学院理学研究科理論天体物理学研究室 (AT 研 ) 修士 2 年 青山尚平 共同研究者 : 市來淨與 新田大輔 杉山直 今回の発表結果は下記の研究成果に基づくものです S.A., K.Ichiki, D.Nitta and N.Sugiyama, ArXiv(2011) [1106.1984] 目次 1. Introduction

Microsoft PowerPoint - シミュレーション工学-2010-第1回.ppt

シミュレーション工学 ( 後半 ) 東京大学人工物工学研究センター 鈴木克幸 CA( Compter Aded geerg ) r. Jaso Lemo (SC, 98) 設計者が解析ツールを使いこなすことにより 設計の評価 設計の質の向上を図る geerg の本質の 計算機による支援 (CA CAM などより広い名前 ) 様々な汎用ソフトの登場 工業製品の設計に不可欠のツール 構造解析 流体解析

Unit 1

Unit 2. 流体としてのプラズマ Klimontovich 方程式系 ( 空間 ) 分布関数 ˆ f (r,v,t) = N [r r j (t)] [v v j (t)] ( = e, i) j1 f ˆ /t + vf ˆ /r + K ˆ f ˆ /v = 0 ˆ E = ˆ B /t ˆ B = 1/c 2 ˆ E /t + 0 ˆ J [ ˆ K (r,v,t) = q /m ( ˆ

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

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

4 19

I / 19 8 1 4 19 : : f(e, J), f(e) Phase mixing Landau Damping, violent relaxation : 2 2 : ( ) http://antwrp.gsfc.nasa.gov/apod/ap950917.html ( ) http://www-astro.physics.ox.ac.uk/~wjs/apm_grey.gif

<4D F736F F F696E74202D208D E9197BF288CF68A4A B8CDD8AB B83685D>

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

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

数値計算法 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

PowerPoint Presentation

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

1/10 平成 29 年 3 月 24 日午後 1 時 37 分第 5 章ローレンツ変換と回転 第 5 章ローレンツ変換と回転 Ⅰ. 回転 第 3 章光速度不変の原理とローレンツ変換 では 時間の遅れをローレンツ変換 ct 移動 v相対 v相対 ct - x x - ct = c, x c 2 移動

/ 平成 9 年 3 月 4 日午後 時 37 分第 5 章ローレンツ変換と回転 第 5 章ローレンツ変換と回転 Ⅰ. 回転 第 3 章光速度不変の原理とローレンツ変換 では 時間の遅れをローレンツ変換 t t - x x - t, x 静止静止静止静止 を導いた これを 図の場合に当てはめると t - x x - t t, x t + x x + t t, x (5.) (5.) (5.3) を得る

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

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

Microsoft PowerPoint - 卒業論文 pptx

時間に依存するポテンシャルによる 量子状態の変化 龍谷大学理工学部数理情報学科 T966 二正寺章指導教員飯田晋司 目次 はじめに 次元のシュレーディンガー方程式 3 井戸型ポテンシャルの固有エネルギーと固有関数 4 4 中央に障壁のある井戸型ポテンシャルの固有エネルギーと固有関数 3 5 障壁が時間によって変化する場合 7 6 まとめ 5 一次元のシュレディンガー方程式量子力学の基本方程式 ψ (

Microsoft Word - NumericalComputation.docx

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

大宇宙

大宇宙 銀河団 大規模構造 膨張宇宙 銀河群 数個 ~ 数十個の銀河の群れ 天の川銀河 250 万光年 アンドロメダ銀河 局所銀河群 http://www.astronomy.com/en/web%20extras/2005/02/ Dominating%20the%20Local%20Group.aspx 銀河団 100 個程度以上の集まり 銀河群との明確な区別はない 天の川銀河 6200 万光年

PowerPoint プレゼンテーション

原始惑星系円盤内でロスビー波不安定性によって形成される渦 小野智弘 ( 京都大 ), 武藤恭之 ( 工学院大 ), 富田賢吾 ( 大阪大 ), 野村英子 ( 東工大 ) Dec. 20th, 2016 理論懇シンポジウム 2016@ 東北大 1 様々な原始惑星系円盤構造 若い星の周りにあるガス円盤 円盤内のダストが合体成長し 惑星を形成 近年 詳細な円盤構造が明らかになってきている ALMA によるダスト連続光観測

FFT

ACTRAN for NASTRAN Product Overview Copyright Free Field Technologies ACTRAN Modules ACTRAN for NASTRAN ACTRAN DGM ACTRAN Vibro-Acoustics ACTRAN Aero-Acoustics ACTRAN TM ACTRAN Acoustics ACTRAN VI 2 Copyright

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

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

Freak wave estimation by WAVEWATCH III and HOSM

2017 年 9 月 25 日第 2 回理研データ同化ワークショップ アンサンブル 4 次元変分法を活用した巨大波浪の再現 藤本航 早稲田卓爾 東京大学新領域創成科学研究科海洋技術環境学専攻 1 導入 : フリーク波とは 海洋に突発的に生じる巨大波 波浪場の平均的な振幅 (H s : 有義波高 ) を大きく超える波 H > 2 or η crest > 1.3 H s H s (Dysthe et

平成 22 年度 革新的な三次元映像技術による超臨場感コミュニケーション技術研究開発 の開発成果について 1. 施策の目標 人体を収容できる大きさの 3 次元音響空間についてリアルタイムに音響レンダリングできるシステム ( シリコンコンサートホール ) を 2013 年までに開発する 具体的には,

平成 22 年度 革新的な三次元映像技術による超臨場感コミュニケーション技術研究開発 の開発成果について 1. 施策の目標 人体を収容できる大きさの 3 次元音響空間についてリアルタイムに音響レンダリングできるシステム ( シリコンコンサートホール ) を 2013 年までに開発する 具体的には, 直方体領域 (2m 2m 4m 程度 ) の室内音場を想定し, 音声周波数帯域 (3kHz まで )

Kumamoto University Center for Multimedia and Information Technologies Lab. 熊本大学アプリケーション実験 ~ 実環境における無線 LAN 受信電波強度を用いた位置推定手法の検討 ~ InKIAI 宮崎県美郷

熊本大学アプリケーション実験 ~ 実環境における無線 LAN 受信電波強度を用いた位置推定手法の検討 ~ InKIAI プロジェクト @ 宮崎県美郷町 熊本大学副島慶人川村諒 1 実験の目的 従来 信号の受信電波強度 (RSSI:RecevedSgnal StrengthIndcator) により 対象の位置を推定する手法として 無線 LAN の AP(AccessPont) から受信する信号の減衰量をもとに位置を推定する手法が多く検討されている

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

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

Microsoft PowerPoint - 2_FrontISTRと利用可能なソフトウェア.pptx

東京大学本郷キャンパス 工学部8号館2階222中会議室 13:30-14:00 FrontISTRと利用可能なソフトウェア 2017年4月28日 第35回FrontISTR研究会 FrontISTRの並列計算ハンズオン 精度検証から並列性能評価まで 観測された物理現象 物理モデル ( 支配方程式 ) 連続体の運動を支配する偏微分方程式 離散化手法 ( 有限要素法, 差分法など ) 代数的な数理モデル

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 流速 圧力を用いた解法

相対性理論入門 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 弨弱弩弨弲弩弨弳弩弨弴弩 が成立する. 図 : 相対速度

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 Ñ Ñ

(Microsoft Word - \216\221\227\277\201i\220\333\223\256\201jv2.doc)

宇宙工学基礎講義資料摂動 ( 松永担当分 ) ベクトル行列演算 ) 微分演算の定義 [ ] ) 微分公式 ( ベクトル記法と行列記法 ) E E ここで E は単位行列 チルダ演算は外積演算と等価の反対称行列を生成する演算 : ( ) ) 恒等演算式 : 次元列ベクトル ( ) ( ) ( ) ( ) ( ) E E ) ( ( )( ) ( )( ) ( )( ) ( ) ( ) ( ) ( )

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

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

Laplace2.rtf

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

風力発電インデックスの算出方法について 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)

線積分.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+

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

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

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

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

宇宙機工学　演習問題

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

On the X-ray and Mass Distribution in the Merging Galaxy Cluster 1E

衝突銀河団に関する話題 : 質量分布 質量評価 磁場進化 滝沢元和 ( 山形大学理学部物理学科 ) 研究会 : マクロでミクロな銀河団 (2007.10.24 26@ 26@ 山形蔵王 たかみや瑠璃倶楽リゾート ) 目次 Introduction ダークマター分布 vs ガス分布 質量評価の不定性について 銀河団磁場の進化 まとめ Introduction(1): 銀河団衝突の痕跡 (X 線 weak

微分方程式 モデリングとシミュレーション

1 微分方程式モデリングとシミュレーション 2018 年度 2 質点の運動のモデル化 粒子と粒子に働く力 粒子の運動 粒子の位置の時間変化 粒子の位置の変化の割合 速度 速度の変化の割合 加速度 力と加速度の結び付け Newtonの運動方程式 : 微分方程式 解は 時間の関数としての位置 3 Newton の運動方程式 質点の運動は Newton の運動方程式で記述される 加速度は力に比例する 2

/ 18 12 4 ( ) ( ) : : f(e, J), f(e) Phase mixing Landau Damping, violent relaxation : 2 2 : ? m i d 2 x i dt 2 = j i f ij (1) x i m i i f ij j i f ij G f ij = Gm i m j x j x i x j x i 3, (2) ( ) 2 (

PowerPoint プレゼンテーション

多倍長計算手法 平成 年度第 四半期 今回はパラメータ の設定と精度に関してまとめて記述しました ループ積分と呼ばれる数値積分計算では 質量 の光子や質量が非常に小さい事はわかっているが その値は不明なニュートリノに対して赤外発散を防ぐため微小量を与えて計算しています この設定する微少量の値により 結果の精度及び反復に要する時間が大きく作用したり 誤った値を得る事があります ここでは典型的な つのケースで説明します

シミュレーション物理4

シミュレーション物理 4 運動方程式の方法 運動方程式 物理で最もよく出てくる そもそも物理はものの運動を議論する学問から出発 ( つり合いは運動を行わないという意味で含まれる ) 代表例 ニュートンの運動方程式 波動方程式 シュレーディンガー方程式 運動方程式 ( 微分方程式の解法 ) 高次の微分方程式を 1 階微分方程式に変形 N 変数の 階微分方程式 N 変数の 1 階微分方程式 dy/dt=f(t,y)

Chap2.key

. f( ) V (V V ) V e + V e V V V V ( ) V V ( ) E. - () V (0 ) () V (0 ) () V (0 ) (4) V ( ) E. - () V (0 ) () V (0 ) O r θ ( ) ( ) : (r θ) : { r cos θ r sn θ { r + () V (0 ) (4) V ( ) θ θ arg( ) : π π

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

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

解析力学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

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 この事から 線形写像の性質を用いると 次の格子上の点全ての写像先が求まる

ଗȨɍɫȮĘർǻ 図 : a)3 次元自由粒子の波数空間におけるエネルギー固有値の分布の様子 b) マクロなサイズの系 L ) における W E) と ΩE) の対応 として与えられる 周期境界条件を満たす波数 kn は kn = πn, L n = 0, ±, ±, 7) となる 長さ L の有限

: Email: mizushima@mp.es.osaka-u.ac.jp, D38 0 08 5 S = k B ln W ) W n [] [] 5 N. 6 d h m dx ϕ nx) = E n ϕ n x) ) L 5 ϕ n x = 0) = ϕ n x = L) = 0, N k n ϕ n = N sink n x), E n = h k n m 3) k n = nπ, n =,,

Microsoft PowerPoint - siryo7

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

Microsoft PowerPoint - 東大講義09-13.ppt [互換モード]

物性物理学 IA 平成 21 年度前期東京大学大学院講義 東京大学物性研究所高田康民 2009 年 4 月 10 日 -7 月 17 日 (15 回 ) 金曜日 2 時限 (10:15-11:45) 15 11 理学部 1 号館 207 号室 講義は自己充足的 量子力学 ( 第 2 量子化を含む ) 統計力学 場の量子論のごく初歩を仮定 最後の約 10 分間は関連する最先端の研究テーマを雑談風に紹介する

<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

Microsoft PowerPoint _量子力学短大.pptx

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

講演番号

平成 25 年度先端的計算科学研究プロジェクト流体型宇宙プラズマシミュレーションコードの性能チューニング 梅田隆行 ( 名古屋大学太陽地球環境研究所 ) 深沢圭一郎 ( 九州大学情報基盤研究開発センター ) 1. 研究の目的と意義太陽から地球に至るジオスペース環境の変動を理解することは 人類の活動が宇宙へと拡大しつつある今日 極めて重要な課題である 人類の活動に影響を与えるジオスペースの変動現象としては

Microsoft Word - thesis.doc

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

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

<4D F736F F D20332E322E332E819C97AC91CC89F090CD82A982E78CA982E9466F E393082CC8D5C91A291CC90AB945C955D89BF5F8D8296D85F F8D F5F E646F63>

3.2.3. 流体解析から見る Fortran90 の構造体性能評価 宇宙航空研究開発機構 高木亮治 1. はじめに Fortran90 では 構造体 動的配列 ポインターなど様々な便利な機能が追加され ユーザーがプログラムを作成する際に選択の幅が広がりより便利になった 一方で 実際のアプリケーションプログラムを開発する際には 解析対象となる物理現象を記述する数学モデルやそれらを解析するための計算手法が内包する階層構造を反映したプログラムを作成できるかどうかは一つの重要な観点であると考えられる

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

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

第 5 章 構造振動学 棒の振動を縦振動, 捩り振動, 曲げ振動に分けて考える. 5.1 棒の縦振動と捩り振動 まっすぐな棒の縦振動の固有振動数 f[ Hz] f = l 2pL である. ただし, L [ 単位 m] は棒の長さ, [ 2 N / m ] 3 r[ 単位 Kg / m ] E r

第 5 章 構造振動学 棒の振動を縦振動, 捩り振動, 曲げ振動に分けて考える 5 棒の縦振動と捩り振動 まっすぐな棒の縦振動の固有振動数 f[ Hz] f l pl である ただし, L [ 単位 m] は棒の長さ, [ N / m ] [ 単位 Kg / m ] E は (5) E 単位は棒の材料の縦弾性係数 ( ヤング率 ) は棒の材料の単位体積当りの質量である l は境界条件と振動モードによって決まる無

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

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

衝突銀河団のＮ体＋ 流体シミュレーション

衝突銀河団の N 体 + 流体シミュレーション 滝沢元和 Introduction 銀河団 銀河団 : 模式図 可視光 ( 銀河 ) X 線 ( 高温ガス ) 暗黒物質の重力ポテンシャル中に束縛された高温ガス (T~10 7-8 K) と銀河のかたまり 宇宙で最大のビリアライズした天体 (R ~Mpc,, M ~10 15 太陽質量 ) 宇宙の構造形成の ( 観測可能な ) 現場 プラズマ物理の実験場

PowerPoint プレゼンテーション

斎藤貴之東京工業大学地球生命研究所 1 Contents 研究背景 ASURA: 並列 N 体 /SPHコード Density Independent SPH まとめ 2 3 内海洋輔(広島大学) NASA/JPL-Caltech/ESO/R. Hurt 4 Credit: S. Beckwith & the HUDF Working Group (STScI), HST, ESA, NASA 我々の宇宙

(Microsoft PowerPoint - \221\34613\211\361)

計算力学 ~ 第 回弾性問題の有限要素解析 (Ⅱ)~ 修士 年後期 ( 選択科目 ) 担当 : 岩佐貴史 講義の概要 全 5 講義. 計算力学概論, ガイダンス. 自然現象の数理モデル化. 行列 場とその演算. 数値計算法 (Ⅰ) 5. 数値計算法 (Ⅱ) 6. 初期値 境界値問題 (Ⅰ) 7. 初期値 境界値問題 (Ⅱ) 8. マトリックス変位法による構造解析 9. トラス構造の有限要素解析. 重み付き残差法と古典的近似解法.

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

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

<4D F736F F F696E74202D20906C8D488AC28BAB90DD8C7689F090CD8D488A D91E F1>

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

ハートレー近似（Hartree aproximation)

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

Microsoft Word - koubo-H26.doc

平成 26 年度学際共同利用プログラム 計算基礎科学プロジェクト 公募要項 - 計算基礎科学連携拠点 ( 筑波大学 高エネルギー加速器研究機構 国立天文台 ) では スーパーコンピュータの学際共同利用プログラム 計算基礎科学プロジェクト を平成 22 年度から実施しております 平成 23 年度からは HPCI 戦略プログラム 分野 5 物質と宇宙の起源と構造 の協力機関である京都大学基礎物理学研究所

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