オープン CAE 関東 数値流体力学 輪講 第 4 回 第 3 章 : 乱流とそのモデリング (3) [3.5~3.7.1 p.64~75] 日時 :2013 年 11 月 10 日 14:00~ 場所 : 日本 新宿 2013/11/10 数値流体力学 輪講第 4 回 1

Size: px
Start display at page:

Download "オープン CAE 関東 数値流体力学 輪講 第 4 回 第 3 章 : 乱流とそのモデリング (3) [3.5~3.7.1 p.64~75] 日時 :2013 年 11 月 10 日 14:00~ 場所 : 日本 新宿 2013/11/10 数値流体力学 輪講第 4 回 1"

Transcription

1 オープン CAE 関東 数値流体力学 輪講 第 4 回 第 3 章 : 乱流とそのモデリング (3 [3.5~3.7.1 p.64~75] 日時 :2013 年 11 月 10 日 14:00~ 場所 : 日本 ESI@ 新宿 1

2 数値流体力学 輪講に関して 目的 数値流体力学の知識 ( 特に理論ベース を深め OpenFOAM の利用に役立てること 本輪講で学ぶもの 数値流体力学の理論や計算手法の概要 2

3 書籍 数値流体力学 第 2 版 原著 : H. K. Versteeg & W. Malalasekera 共訳 : 松下洋介 斎藤泰洋青木秀之 三浦隆利 出版社 : 森北出版株式会社 出版年月 : 2011 年 7 月価格 : 9975 円 高い ページ数 : 544ページ 量が多い 有限体積法を説明した書籍 ( 和書 の中では 最も丁寧に記述されている 3

4 本日 日程パート部分ページ 第 3 章 : 乱流とそのモデリング担当セクション :3.5~3.7.1 p.64~75 3.5~3.6 ごめんなさい! 全部まとめられませんでした p.64~69 今回は 3.5~3.6 節 (p に絞って説明したいと思います 3.7 節以降は次回ということで 4

5 内容 乱流流れの計算 3.6 (p.68~69 レイノルズ平均ナビエ - ストークス方程式 (RANS 3.5 (p.64~68 5

6 内容 乱流流れの計算 3.6 (p.68~69 レイノルズ平均ナビエ - ストークス方程式 (RANS 3.5 (p.64~68 6

7 乱流流れの特徴 流速計プローブ 流速の実測値 流速の時間平均値 ローパスフィルターをかけた流速 乱流噴流 画像 :wikipedia より 乱流の速度変動の要因 スムージングの一種 ( 空間平均化 t 乱流中には 大小様々な渦構造が存在するため 7

8 エネルギーカスケードと コルモゴロフスケール エネルギーカスケード エネルギー供給 乱流の運動エネルギーは 大きな渦から小さな渦へ ( ほとんど散逸なく 輸送され そのエネルギーは最小渦に至って熱に代わる コルモゴロフスケール 上記の最小渦の大きさを次のように見積もることができる 分子動粘性 l d 散逸率 3 1/ 4 最小渦 ( 大きさ l d エネルギーカスケード l d をコルモゴロフスケールと呼び 乱流の最小渦スケールとして認識されている 散逸 ( 熱 8

9 乱流流れの計算 (1 直接数値シミュレーション (Direct Nmerical Simlation: DNS コルモゴロフスケールまで計算格子を細かく切り 乱流構造のすべてを直接計算する方法 ラージエディシミュレーション (Large Eddy Simlation: LES 設定した計算格子よりも細かい変動スケールはモデル化し 計算格子よりも大きなスケールを直接計算する方法 レイノルズ平均ナビエ - ストークス方程式 (Reynolds-Averaged Navier-Stokes eqation: RANS 乱流のすべての変動成分をモデル化し 乱流の平均流を計算する方法 9

10 乱流流れの計算 (2 DNS のターゲット 流速の実測値 RANS のターゲット 流速の時間平均値 LES のターゲット ローパスフィルターをかけた流速 スムージングの一種 ( 空間平均化 t 乱流計算手法モデル化の有無モデル化のターゲット DNS なし LES 格子より小さなスケールの乱流変動成分 RANS すべての乱流変動成分 10

11 内容 乱流流れの計算 3.6 (p.68~69 レイノルズ平均ナビエ - ストークス方程式 (RANS 3.5 (p.64~68 11

12 レイノルズ分解 流体力学の基礎 ( 第 8 回 の資料より 乱流中の瞬時流体速度 :, v また v 時間平均を取るという意味 v 平均流速 v 流速の変動成分 1 T 1 dt v T vdt ( 0 T 0 T 0 1 ' T 1 ' dt 0 v' T v' dt 0 T 0 T 0 上記のように と物理量の平均値と変動成分を 分解する考え方を レイノルズ分解と呼ぶ 12

13 乱流の速度変動の性質 流体力学の基礎 ( 第 8 回 の資料より ランダムな変動 v 0 ランダムな変動 v 実際の流体の変動 v 0 実際の流体の変動 v 負の相関 流体速度の変動成分の積 ( 二次モーメント は 0 ではない 13

14 レイノルズ応力 (1 流体力学の基礎 ( 第 8 回 の資料より y 方向の流体の運動量輸送を考える ( 右図 流体塊が有する単位時間 単位体積当たりの x 方向の運動量 ( ' 単位時間当たりに y 方向に輸送される流体塊の体積 V vs ( v v' S S=1 v=0 なので V v' y 2 y y 1 0 y v = 0 2 v 2 v S: 単位面積 14

15 レイノルズ応力 (2 流体力学の基礎 ( 第 8 回 の資料より y 方向の流体の運動量輸送を考える ( 右図 単位時間当たりに y 方向に輸送される流体塊の (x 方向の 運動量 V ( ' v' 1 T Vdt 0 T T 時間平均 T 0 T ( T 0 ' v' dt v' dt T T 0 ' v' dt y 2 y y 1 0 y ' v' v = 0 2 v 2 v 応力発生 2 S: 単位面積 15

16 レイノルズ応力 (3 流体力学の基礎 ( 第 8 回 の資料より 今 y 2 の面には y の面から y 方向に v (>0 の速度を持つ流体塊が到達する この流体塊は x 方向速度が である面から 更に大きい速度 2 を持つ面に移動している この結果 y 2 面での x 方向速度が減速し 2 (<0 の変動速度が発生する 16 y 2 y y 1 0 y v = 0 2 v 2 v S: 単位面積 つまり 互いに関連のある と v の符号は逆であるから v の時間平均値 v は負である したがって 乱流の速度変動がもたらす応力は次のように示す xy 'v' レイノルズ応力

17 レイノルズ応力 (4 前ページでは xy の導出を行ったが レイノルズ応力は応力テンソルと同様に 2 階テンソルを有しているので 具体的には次のように示される xx yx zx xy yy zy xz yz zz v w v vv wv w vw w w v なので v w w xy yx xz zx vw xz zx つまり レイノルズ応力は対象テンソルを取る wv 17

18 交換則 (1 ある変数 ( スカラー変数 j, yのレイノルズ分解 j j y y j, yの時間平均では 次のような関係 ( 交換則 が成立する 基本 : j y 0 j j 0 加減即 : 微積分 : 積 : 平均値変動平均値変動 j y j y j j s s s j 0 18 jds jds ds j ( j jy ( j( y yj jy jy jy 0

19 交換則 (2 ベクトル量についても前ページと同じ交換則が成立する あるベクトル a のレイノルズ分解 a A a とする ベクトル解析演算子 : div a div a div A div( ja div div(grad j div(grad j a div( A div j a 19

20 レイノルズ平均ナビエ - ストークス方程式 (RANS の導出 (1 直交座標系 ( カルテシアン座標系 における 非圧縮性流体の基礎方程式 連続の式 : ナビエ - ストークス方程式 div 0 1 p div( div(grad t x v 1 p div( v div(grad v t y w 1 p div( w div(grad w t z =(, v, w: (x, y, z 方向の流速 t: 時間 : 密度 ( 定数 p: 圧力 : 動粘性係数 ( 定数 20

21 レイノルズ平均ナビエ - ストークス方程式 (RANS の導出 (2 2013/11/10 21 数値流体力学 輪講第 4 回 =(, v, w p についてレイノルズ分解をし 時間平均をとる w v W V w v 各変数のレイノルズ分解 p P p, div 0 div(grad 1 div( div( x P t div(grad 1 div( div( V y P v V t V div(grad 1 div( div(w W z P w t W 新たな項が発生レイノルズ平均ナビエ - ストークス方程式 (RANS 連続の式 : 平均値変動平均値変動平均値変動

22 レイノルズ平均ナビエ - ストークス方程式 (RANS の導出 (3 2013/11/10 22 数値流体力学 輪講第 4 回レイノルズ応力を用いて表現すると div(grad 1 div( x P t div(grad 1 div( V y P V t V div(grad 1 div(w W z P t W レイノルズ平均ナビエ - ストークス方程式 (RANS 連続の式 : 0 div xz xy xx z y x 1 yz yy yx z y x 1 zz zy zx z y x 1

23 スカラー変数の輸送方程式 前ページの式に発生した項を レイノルズ応力を用いて表現すると RANS 方程式と類似の形式をとる ある変数 ( スカラー変数 j を考える レイノルズ分解 : j j 平均値 変動 t div( 1 div( 拡散係数 grad div( j S ソース項 23

24 圧縮性 ( 密度変動 の考慮 (1 これまでは 非圧縮性流体として 密度一定として RANS 方程式を取り扱ってきたが ここからは密度変動がある場合を想定して RANS 方程式を考える 密度のレイノルズ分解は厳しい 密度のレイノルズ分解を考慮した連続の式 t 密度のレイノルズ分解を考慮したナビエ - ストーク方程式 t ( div( div( ( ( div( div( div( t 取扱難取扱難 P div(grad x 取り扱いの難しい項が発生し いろいろ大変 24 0 湧き出しの発生

25 圧縮性 ( 密度変動 の考慮 (2 Morkovin の仮説 境界層流れでは主流のマッハ数が 5 以下 噴流ではマッハ 1.5 以下であれば 圧縮性による乱流構造の変化は無視でき その範囲内であれば ファーブル平均が使える ファーブル平均 ( 密度加重平均 変数 x のファーブル平均 ~ x x ファーブル平均を用いることで 密度と変数の平均値の処理を分離することができる 25

26 圧縮性流れにおけるレイノルズ平均 ナビエ - ストークス方程式 (1 レイノルズ平均ナビエ - ストークス方程式 (RANS ~ div t ~ ~ ~ ~ P ~ div( div( grad t x ( ( v ( w x y z ~ ~ V ~~ P ~ div( V div( gradv t y ( v ( vv ( vw x y z ~ ~ W ~ ~ P ~ div( W div( gradw t z ( w ( wv ( ww x y z 連続の式 : 0 S Mx S My S Mz 26

27 圧縮性流れにおけるレイノルズ平均 ナビエ - ストークス方程式 (2 スカラー変数の輸送方程式 ~ t ~ ~ div( x div( ~ gradw ( j ( jv ( jw y z S 以上の方程式は 低マッハ数での圧縮性乱流解析でよく用いられる 27

28 次回 日程パート部分ページ 第 3 章 : 乱流とそのモデリング担当セクション :3.7~3.7.1 p.69~75 次回は 今回できなかったところ ( 乱流モデルの概要と混合長理論 をやりたいと思います 引き続き 私が輪講を担当します 28

オープン CAE 関東 数値流体力学 輪講 第 6 回 第 3 章 : 乱流とそのモデリング (5) [3.7.2 p.76~84] 日時 :2014 年 2 月 22 日 14:00~ 場所 : 日本 新宿 2013/02/22 数値流体力学 輪講第 6 回 1

オープン CAE 関東 数値流体力学 輪講 第 6 回 第 3 章 : 乱流とそのモデリング (5) [3.7.2 p.76~84] 日時 :2014 年 2 月 22 日 14:00~ 場所 : 日本 新宿 2013/02/22 数値流体力学 輪講第 6 回 1 オープン CAE 勉強会 @ 関東 数値流体力学 輪講 第 6 回 第 章 : 乱流とそのモデリング (5) [.7. p.76~84] 日時 :04 年 月 日 4:00~ 場所 : 日本 ESI@ 新宿 本日 日程パート部分ページ 04.0 第 章 : 乱流とそのモデリング担当セクション :.7. p.76~84 今回は北風が担当しました ご質問 記述ミス等に関するご指摘がありましたら 以下までご連絡下さい

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

ナビエ・ストークス方程式

ナビエ・ストークス方程式 Vol. 4 ナビエ ストークス方程式 序 Dt = F 1 ρ gradp 1 3 νgradθνδv 会員堀城之 これが, ナビエ ストークス方程式である 主に流体力学において用いられる 土木工学, 建築工学, 機械工学, 航空工学, 船舶工学, 物理学を履修した方は学んだかもしれない 流体に係る発明も少なくないであろう しかし, 大学院の流体研究室或いは, 当該数学研究室に入らない限り圧縮項

More information

<4D F736F F F696E74202D208BAB8A458FF08C8F82CC8AEE916282C68C8892E896402E707074>

<4D F736F F F696E74202D208BAB8A458FF08C8F82CC8AEE916282C68C8892E896402E707074> No.07-131 講習会 ( 流体工学部門企画 ) 境界条件の基礎と決定法 千葉科学大学 戸田和之 講演の流れ 数値解析とは何か 境界条件の役割と目的 境界の分類 計算法による 設定の違い 非圧縮流れ解析における境界条件の設定法 乱流解析における境界条件の設定法 圧縮性流れ解析における境界条件の設定法 1 流れの数値解析とは 偏微分型で書かれた基礎方程式を解く作業 連続の式 υ = 0 υ: 速度ベクトル

More information

応力とひずみ.ppt

応力とひずみ.ppt in yukawa@numse.nagoya-u.ac.jp 2 3 4 5 x 2 6 Continuum) 7 8 9 F F 10 F L L F L 1 L F L F L F 11 F L F F L F L L L 1 L 2 12 F L F! A A! S! = F S 13 F L L F F n = F " cos# F t = F " sin# S $ = S cos# S S

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 中部 CAE 懇話会 流体伝熱基礎講座 第 3 回午後 名古屋工業大学大学院 創成シミュレーション工学専攻 後藤俊幸 粘性流体 H y U A F u(y,t) -F x 単位面積当たりのせん断応力 Newton 流体 t 線形関係 応力テンソル t ij 力 力の方向 面 ( 法線 ) z n=(0,0,1) t zz t yz t xz n=(0,1,0) y t yy t zy t xy t

More information

変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy,

変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy, 変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy, z + dz) Q! (x + d x + u + du, y + dy + v + dv, z +

More information

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

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

More information

テンソル ( その ) テンソル ( その ) スカラー ( 階のテンソル ) スカラー ( 階のテンソル ) 階数 ベクトル ( 階のテンソル ) ベクトル ( 階のテンソル ) 行列表現 シンボリック表現 [ ]

テンソル ( その ) テンソル ( その ) スカラー ( 階のテンソル ) スカラー ( 階のテンソル ) 階数 ベクトル ( 階のテンソル ) ベクトル ( 階のテンソル ) 行列表現 シンボリック表現 [ ] Tsor th-ordr tsor by dcl xprsso m m Lm m k m k L mk kk quott rul by symbolc xprsso Lk X thrd-ordr tsor cotrcto j j Copyrght s rsrvd. No prt of ths documt my b rproducd for proft. テンソル ( その ) テンソル ( その

More information

物性物理学 I( 平山 ) 補足資料 No.6 ( 量子ポイントコンタクト ) 右図のように 2つ物質が非常に小さな接点を介して接触している状況を考えましょう 物質中の電子の平均自由行程に比べて 接点のサイズが非常に小さな場合 この接点を量子ポイントコンタクトと呼ぶことがあります この系で左右の2つ

物性物理学 I( 平山 ) 補足資料 No.6 ( 量子ポイントコンタクト ) 右図のように 2つ物質が非常に小さな接点を介して接触している状況を考えましょう 物質中の電子の平均自由行程に比べて 接点のサイズが非常に小さな場合 この接点を量子ポイントコンタクトと呼ぶことがあります この系で左右の2つ 物性物理学 I( 平山 ) 補足資料 No.6 ( 量子ポイントコンタクト ) 右図のように つ物質が非常に小さな接点を介して接触している状況を考えましょう 物質中の電子の平均自由行程に比べて 接点のサイズが非常に小さな場合 この接点を量子ポイントコンタクトと呼ぶことがあります この系で左右のつの物質の間に電位差を設けて左から右に向かって電流を流すことを行った場合に接点を通って流れる電流を求めるためには

More information

Microsoft PowerPoint - 乱流の数値解析2010_02.ppt

Microsoft PowerPoint - 乱流の数値解析2010_02.ppt 00 数値流体力学 目次 00 数値流体力学 乱流の数値解析 ~ 乱流の物理モデル ~ 平成 年 月 日筑波大学大学院システム情報工学研究科金子暁子 aneo@zsbaa 乱流とは何か 乱流の特徴 乱流の統計的表現 乱流の表現 流れのエネルギーとエネルギー方程式 乱流の特性と数値シミュレーション シミュレーション法の分類 渦粘性の概念 0- 方程式モデル - 方程式モデル - 方程式モデル LES

More information

Microsoft Word - 9章(分子物性).doc

Microsoft Word - 9章(分子物性).doc 1/1/6 9 章分子物性 1 節電気双極子モーメント (Electric Dipole Moment) 電子双極子モーメント とは 微小な距離 a だけ離れて点電荷 q が存在する状態 絶対値は aq で 負電荷 q から正電荷 q へ向かうベクトルである 例えば 水分子は下右図のような向きの電気双極子モーメントをもち その大きさは約 1.85D である このように元々から持っている双極子モーメントを

More information

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

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

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 Word - thesis.doc

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

More information

技術資料 JARI Research Journal OpenFOAM を用いた沿道大気質モデルの開発 Development of a Roadside Air Quality Model with OpenFOAM 木村真 *1 Shin KIMURA 伊藤晃佳 *2 Akiy

技術資料 JARI Research Journal OpenFOAM を用いた沿道大気質モデルの開発 Development of a Roadside Air Quality Model with OpenFOAM 木村真 *1 Shin KIMURA 伊藤晃佳 *2 Akiy 技術資料 176 OpenFOAM を用いた沿道大気質モデルの開発 Development of a Roadside Air Quality Model with OpenFOAM 木村真 *1 Shin KIMURA 伊藤晃佳 *2 Akiyoshi ITO 1. はじめに自動車排出ガスの環境影響は, 道路沿道で大きく, 建物など構造物が複雑な気流を形成するため, 沿道大気中の自動車排出ガス濃度分布も複雑になる.

More information

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

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

More information

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

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

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

応用数学A

応用数学A 応用数学 A 米田 戸倉川月 7 限 1930~2100 西 5-109 V を :x 2 + y 2 + z 2 = 4 で囲まれる内部とする F = ye x xe y + ze z FdV = V e x e y e z F = = 2e z 2e z dv = 2e z 3 23 = 64π 3 e z y x z 4π V n Fd = 1 F nd 2 F nd 法線ベクトル n g x,

More information

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

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

More information

vecrot

vecrot 1. ベクトル ベクトル : 方向を持つ量 ベクトルには 1 方向 2 大きさ ( 長さ ) という 2 つの属性がある ベクトルの例 : 物体の移動速度 移動量電場 磁場の強さ風速力トルクなど 2. ベクトルの表現 2.1 矢印で表現される 矢印の長さ : ベクトルの大きさ 矢印の向き : ベクトルの方向 2.2 2 個の点を用いて表現する 始点 () と終点 () を結ぶ半直線の向き : ベクトルの方向

More information

.....Z...^.[.......\..

.....Z...^.[.......\.. 15 10 16 42 55 55 56 60 62 199310 1995 134 10 8 15 1 13 1311 a s d f 141412 2 g h j 376104 3 104102 232 4 5 51 30 53 27 36 6 Y 7 8 9 10 8686 86 11 1310 15 12 Z 13 14 15 16 102193 23 1712 60 27 17 18 Z

More information

座標変換におけるテンソル成分の変換行列

座標変換におけるテンソル成分の変換行列 座標変換におけるテンソル成分の変換行列 座標変換におけるテンソル成分の変換関係は 次元数によらず階数によって定義される変換行列で整理することができる 位置ベクトルの変換行列を D としてそれを示そう D の行列式を ( = D ) とするとき 鏡映や回映といった pseudo rotation に対しては = -1 である が問題になる基底は 対称操作に含まれる pseudo rotation に依存する

More information

<4D F736F F F696E74202D E94D58B9393AE82F AC82B782E982BD82DF82CC8AEE E707074>

<4D F736F F F696E74202D E94D58B9393AE82F AC82B782E982BD82DF82CC8AEE E707074> 地盤数値解析学特論 防災環境地盤工学研究室村上哲 Mrakam, Satoh. 地盤挙動を把握するための基礎. 変位とひずみ. 力と応力. 地盤の変形と応力. 変位とひずみ 変形勾配テンソルひずみテンソル ひずみテンソル : 材料線素の長さの 乗の変化量の尺度 Green-Lagrange のひずみテンソルと Alman のひずみテンソル 微小変形状態でのひずみテンソル ひずみテンソルの物理的な意味

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

<4D F736F F D E B82CC89DF8B8E81458CBB8DDD814596A297882E646F63>

<4D F736F F D E B82CC89DF8B8E81458CBB8DDD814596A297882E646F63> 5 特集 RANS RANS モデルの過去 現在 未来 * 東京理科大学工学部山本誠 Pas Presen and Fuure of RANS Model Maoo YAMAMOTO Faculy of Engneerng Toyo Unversy of Scence はじめに乱流を計算するために様々な計算手法が開発 利用されているが レイノルズ平均 ( あるいは時間平均 ) に基づくものを Reynolds-Averaged

More information

1/30 平成 29 年 3 月 24 日 ( 金 ) 午前 11 時 25 分第三章フェルミ量子場 : スピノール場 ( 次元あり ) 第三章フェルミ量子場 : スピノール場 フェルミ型 ボーズ量子場のエネルギーは 第二章ボーズ量子場 : スカラー場 の (2.18) より ˆ dp 1 1 =

1/30 平成 29 年 3 月 24 日 ( 金 ) 午前 11 時 25 分第三章フェルミ量子場 : スピノール場 ( 次元あり ) 第三章フェルミ量子場 : スピノール場 フェルミ型 ボーズ量子場のエネルギーは 第二章ボーズ量子場 : スカラー場 の (2.18) より ˆ dp 1 1 = / 平成 9 年 月 日 ( 金 午前 時 5 分第三章フェルミ量子場 : スピノール場 ( 次元あり 第三章フェルミ量子場 : スピノール場 フェルミ型 ボーズ量子場のエネルギーは 第二章ボーズ量子場 : スカラー場 の (.8 より ˆ ( ( ( q -, ( ( c ( H c c ë é ù û - Ü + c ( ( - に限る (. である 一方 フェルミ型は 成分をもち その成分を,,,,

More information

96 7 1m =2 10 7 N 1A 7.1 7.2 a C (1) I (2) A C I A A a A a A A a C C C 7.2: C A C A = = µ 0 2π (1) A C 7.2 AC C A 3 3 µ0 I 2 = 2πa. (2) A C C 7.2 A A

96 7 1m =2 10 7 N 1A 7.1 7.2 a C (1) I (2) A C I A A a A a A A a C C C 7.2: C A C A = = µ 0 2π (1) A C 7.2 AC C A 3 3 µ0 I 2 = 2πa. (2) A C C 7.2 A A 7 Lorentz 7.1 Ampère I 1 I 2 I 2 I 1 L I 1 I 2 21 12 L r 21 = 12 = µ 0 2π I 1 I 2 r L. (7.1) 7.1 µ 0 =4π 10 7 N A 2 (7.2) magnetic permiability I 1 I 2 I 1 I 2 12 21 12 21 7.1: 1m 95 96 7 1m =2 10 7 N

More information

線形弾性体 線形弾性体 応力テンソル とひずみテンソルソル の各成分が線形関係を有する固体. kl 応力テンソル O kl ひずみテンソル

線形弾性体 線形弾性体 応力テンソル とひずみテンソルソル の各成分が線形関係を有する固体. kl 応力テンソル O kl ひずみテンソル Constitutive equation of elasti solid Hooke s law λδ μ kk Lame s onstant λ μ ( )( ) ( ) linear elasti solid kl kl Copyright is reserved. No part of this doument may be reprodued for profit. 線形弾性体 線形弾性体

More information

多変量解析 ~ 重回帰分析 ~ 2006 年 4 月 21 日 ( 金 ) 南慶典

多変量解析 ~ 重回帰分析 ~ 2006 年 4 月 21 日 ( 金 ) 南慶典 多変量解析 ~ 重回帰分析 ~ 2006 年 4 月 21 日 ( 金 ) 南慶典 重回帰分析とは? 重回帰分析とは複数の説明変数から目的変数との関係性を予測 評価説明変数 ( 数量データ ) は目的変数を説明するのに有効であるか得られた関係性より未知のデータの妥当性を判断する これを重回帰分析という つまり どんなことをするのか? 1 最小 2 乗法により重回帰モデルを想定 2 自由度調整済寄与率を求め

More information

Microsoft PowerPoint - 第3回OpenCAE初歩情報交換会@北東北_若嶋2.pptx

Microsoft PowerPoint - 第3回OpenCAE初歩情報交換会@北東北_若嶋2.pptx 調査報告 乱流モデルの選択および設定について 一関高専 若嶋 OpenFOAM 2.3.x についてのみ調査 2014/12/5 第 3 回 OpenCAE 初歩情報交換会 @ 北東北 1 OpenFOAM で設定できる乱流モデル http://www.openfoam.org/features/turbulence.php Incompressible Compressible RAS(RANS)

More information

<4D F736F F D E4F8E9F82C982A882AF82E98D7397F1>

<4D F736F F D E4F8E9F82C982A882AF82E98D7397F1> 3 三次における行列 要旨高校では ほとんど 2 2 の正方行列しか扱ってなく 三次の正方行列について考えてみたかったため 数 C で学んだ定理を三次の正方行列に応用して 自分たちで仮説を立てて求めていったら 空間における回転移動を表す行列 三次のケーリー ハミルトンの定理 三次における逆行列を求めたり 仮説をたてることができた. 目的 数 C で学んだ定理を三次の正方行列に応用する 2. 概要目的の到達点として

More information

回転.rtf

回転.rtf ベクトルの回転の定義は A A rot A ΔS 0 n ΔS (5.) : ounter dl 図 5. ベクトルの回転 であり, 回転量を調べる演算子である. ローテーションA, カールA,Aの回転とも読む. 図 5.のように, 閉曲線 に沿ってベクトル Aの線積分を行うものとする. 線積分はベクトル Aと線素 dl の内積だから, ある大きさ ( スカラー量 ) が得られる. その大きさをもち,

More information

物理学 II( 熱力学 ) 期末試験問題 (2) 問 (2) : 以下のカルノーサイクルの p V 線図に関して以下の問題に答えなさい. (a) "! (a) p V 線図の各過程 ( ) の名称とそのと (& きの仕事 W の面積を図示せよ. # " %&! (' $! #! " $ %'!!!

物理学 II( 熱力学 ) 期末試験問題 (2) 問 (2) : 以下のカルノーサイクルの p V 線図に関して以下の問題に答えなさい. (a) ! (a) p V 線図の各過程 ( ) の名称とそのと (& きの仕事 W の面積を図示せよ. #  %&! (' $! #!  $ %'!!! 物理学 II( 熱力学 ) 期末試験問題 & 解答 (1) 問 (1): 以下の文章の空欄に相応しい用語あるいは文字式を記入しなさい. 温度とは物体の熱さ冷たさを表す概念である. 物体は外部の影響を受けなければ, 十分な時間が経過すると全体が一様な温度の定常的な熱平衡状態となる. 物体 と物体 が熱平衡にあり, 物体 と物体 が熱平衡にあるならば, 物体 と物体 も熱平衡にある. これを熱力学第 0

More information

スライド 1

スライド 1 暫定版修正 加筆の可能性あり ( 付録 ) マクスウェルの応力テンソル (). ある領域に作用する力 2. 応力テンソル 3. 力の総和と応力テンソル 4. ローレンツ力 5. マクスウェルの方程式 6. 孤立系 注意. 本付録 : マクスウェルの応力テンソル(stress tesor) 2. 簡単のため 個々の電荷が真空中をバラバラに運動する孤立系を考えます 3. 背景は真空とします 真空中の誘電率と透磁率を使用します

More information

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

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

More information

Microsoft PowerPoint - Š’Š¬“H−w†i…„…C…m…‰…Y’fl†j.ppt

Microsoft PowerPoint - Š’Š¬“H−w†i…„…C…m…‰…Y’fl†j.ppt 乱流とは? 不規則運動であり, 速度の時空間的な変化が複雑であり, 個々の測定結果にはまったく再現性がなく, 偶然の値である. 渦運動 3 次元流れ 非定常流 乱流は確率過程 (Stochastic Process) である. 乱流工学 1 レイノルズの実験 UD = = ν 慣性力粘性力 乱流工学 F レイノルズ数 U L / U 3 = mα = ρl = ρ 慣性力 L U u U A = µ

More information

<4D F736F F D B89B98D758F4B89EF5F B89B991CE8DF42E646F63>

<4D F736F F D B89B98D758F4B89EF5F B89B991CE8DF42E646F63> 空力騒音の予測と低減技術 豊橋技術科学大学機械工学系教授飯田明由 (iida@mech.tut.ac.jp) はじめにファンなどの工業製品の空力音の解析では非圧縮性の流れ場の計算結果から空力音を予測する分離解法が用いられる場合が多い. この要因の一つとして, 流れ場と音場の代表スケールが大きく異なることがあげられる. たとえば, 時速 100km 程度で走行する自動車周りの流れと空力騒音について考える場合,

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

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

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

More information

伝熱学課題

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

More information

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

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

More information

No δs δs = r + δr r = δr (3) δs δs = r r = δr + u(r + δr, t) u(r, t) (4) δr = (δx, δy, δz) u i (r + δr, t) u i (r, t) = u i x j δx j (5) δs 2

No δs δs = r + δr r = δr (3) δs δs = r r = δr + u(r + δr, t) u(r, t) (4) δr = (δx, δy, δz) u i (r + δr, t) u i (r, t) = u i x j δx j (5) δs 2 No.2 1 2 2 δs δs = r + δr r = δr (3) δs δs = r r = δr + u(r + δr, t) u(r, t) (4) δr = (δx, δy, δz) u i (r + δr, t) u i (r, t) = u i δx j (5) δs 2 = δx i δx i + 2 u i δx i δx j = δs 2 + 2s ij δx i δx j

More information

Microsoft Word - EM_EHD_2010.doc

Microsoft Word - EM_EHD_2010.doc H のための電磁気学 機能材料工学科阿部洋 . 電磁気学電磁気学電磁気学電磁気学の基礎基礎基礎基礎 - マクスウェルマクスウェルマクスウェルマクスウェルの応力応力応力応力静電場の条件は e div ρ ( ) ot ( ) である 体積 V で電荷密度 ρ e に働く力はクーロン力から ρ dv F e ( 3) と表せる ( 3) 式に ( ) を代入すると ( ) dv div F ( 4) となる

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

Laplace2.rtf

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

More information

7-12.dvi

7-12.dvi 26 12 1 23. xyz ϕ f(x, y, z) Φ F (x, y, z) = F (x, y, z) G(x, y, z) rot(grad ϕ) rot(grad f) H(x, y, z) div(rot Φ) div(rot F ) (x, y, z) rot(grad f) = rot f x f y f z = (f z ) y (f y ) z (f x ) z (f z )

More information

技術者のための構造力学 2014/06/11 1. はじめに 資料 2 節点座標系による傾斜支持節点節点の処理 三好崇夫加藤久人 従来, マトリックス変位法に基づく骨組解析を紹介する教科書においては, 全体座標系に対して傾斜 した斜面上の支持条件を考慮する処理方法として, 一旦, 傾斜支持を無視した

技術者のための構造力学 2014/06/11 1. はじめに 資料 2 節点座標系による傾斜支持節点節点の処理 三好崇夫加藤久人 従来, マトリックス変位法に基づく骨組解析を紹介する教科書においては, 全体座標系に対して傾斜 した斜面上の支持条件を考慮する処理方法として, 一旦, 傾斜支持を無視した . はじめに 資料 節点座標系による傾斜支持節点節点の処理 三好崇夫加藤久人 従来, マトリックス変位法に基づく骨組解析を紹介する教科書においては, 全体座標系に対して傾斜 した斜面上の支持条件を考慮する処理方法として, 一旦, 傾斜支持を無視した全体座標系に関する構造 全体の剛性マトリックスを組み立てた後に, 傾斜支持する節点に関して対応する剛性成分を座標変換に よって傾斜方向に回転処理し, その後は通常の全体座標系に対して傾斜していない支持点に対するのと

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

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

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

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

Microsoft Word - 中村工大連携教材(最終 ).doc

Microsoft Word - 中村工大連携教材(最終 ).doc 音速について考えてみよう! 金沢工業大学 中村晃 ねらい 私たちの身の回りにはいろいろな種類の波が存在する. 体感できる波もあれば, できない波もある. その中で音は体感できる最も身近な波である. 遠くで雷が光ってから雷鳴が届くまで数秒間時間がかかることにより, 音の方が光より伝わるのに時間がかかることも経験していると思う. 高校の物理の授業で音の伝わる速さ ( 音速 ) は約 m/s で, 詳しく述べると

More information

Microsoft Word - Chap17

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

More information

PowerPoint Presentation

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

More information

スライド 1

スライド 1 流体の力学の基礎 三重大学 大学院生物資源学研究科共生環境学専攻地球環境気候学研究室教授立花義裕 イントロ ゆく河の流れは絶えずして しかももとの水にあらず 淀みに浮ぶうたかたは かつ消え かつ結びて 久しくとどまりたる例なし 世の中にある人と 栖とまたかくのごとし ( 方丈記 鴨長明 ) 美しい地球大気の流れ 海洋の流れ 海氷の流れ 川の流れ 水道管やガス管の水やガスの流れ 血管の中の血液の流れ

More information

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

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

More information

<4D F736F F D20837D834E B95FB92F68EAE>

<4D F736F F D20837D834E B95FB92F68EAE> マクスウエルの方程式 Akio Arimoto, Monday, November, 7. イントロ長野 []p.4 に証明抜きで以下のような解説がある 次節以下これを証明していきたいと思う grad f «df d dx =,, rot «( i i), [ ] div «d ( dx dx + dx dx + dx dx ) æ f f f æ f f f rot grad f = rot( df

More information

<4D F736F F D2097CD8A7793FC96E582BD82ED82DD8A E6318FCD2E646F63>

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

More information

平面波

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

More information

ニュートン重力理論.pptx

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

More information

スライド 1

スライド 1 ディスカッション ( データ同化相談会 ) 縮約モデルと粒子フィルタを用いた 低層風擾乱の流れ場予測 菊地亮太 ( 東北大学流体科学研究所 ) 1 研究背景低層風擾乱の運航影響 空港周辺での低層風擾乱による運航影響が問題となっている 着陸復行 代替地着陸 ハードランディング 地形や ビルなどの影響を強く受ける 小規模な気象現象 地形性乱気流 2 航空機運航において重要な課題となっている 研究背景低層風擾乱の解析

More information

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

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

More information

RIST ニュース No.64(2018) ビット演算による CFD( 数値流体力学 ) と等価な高精度流体解析手法 AFluidAnalysisMethodbybitwiseoperations forachievinghighaccuracyofcfd 高度情報科学技術研究機構松岡浩 流体解析の

RIST ニュース No.64(2018) ビット演算による CFD( 数値流体力学 ) と等価な高精度流体解析手法 AFluidAnalysisMethodbybitwiseoperations forachievinghighaccuracyofcfd 高度情報科学技術研究機構松岡浩 流体解析の ビット演算による CFD( 数値流体力学 ) と等価な高精度流体解析手法 AFluidAnalysisMethodbybitwiseoperations forachievinghighaccuracyofcfd 高度情報科学技術研究機構松岡浩 流体解析の時間発展計算をビット演算で超高速に実行する方法として 格子ガス法 がある まず 流体が存在する空間中に格子を張り 多数の仮想粒子を格子点上に配置する

More information

Microsoft Word - Freefem減ページ原稿.doc

Microsoft Word - Freefem減ページ原稿.doc 月刊下水道 2015 11 月号 VOL.38 No.13 有限要素法 (FreeFem++) による三次元流体解析 - 手軽に流れを観察するその2 - 中日本建設コンサルタント ( 株 ) 中根進 1. まえがき筆者は 本誌 Vol.36.No.10 (2013 年増刊号 ) で 格子ボルツマン法による下水流れの可視化 - 手軽に流れを観察する- と題して フリーソフト(Blender:Fluid)

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

Microsoft PowerPoint - 第2回半導体工学

Microsoft PowerPoint - 第2回半導体工学 17 年 1 月 16 日 月 1 限 8:5~1:15 IB15 第 回半導体工学 * バンド構造と遷移確率 天野浩 項目 1 章量子論入門 何故 Si は光らず GN は良く光るのか? *MOSFET ゲート SiO / チャネル Si 界面の量子輸送過程 MOSFET には どのようなゲート材料が必要なのか? http://www.iue.tuwien.c.t/ph/vsicek/noe3.html

More information

1.1 1 A

1.1 1 A . A..2 2 2. () (xyz) ( xyz) ( xy z) = (x x)yz ( xy z) = yz ( xy z) = y(z ( x z)) = y((z x)(z z)) = y( x z) (2) (3) M aj (x, y, M aj ( x, ȳ, z)) = xy ȳm aj ( x, ȳ, z) M aj ( x, ȳ, z)x M aj (x, y, z) x =

More information

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

Microsoft PowerPoint - LectureB1handout.ppt [互換モード] 本講義のスコープ 都市防災工学 後半第 回 : イントロダクション 千葉大学大学院工学研究科建築 都市科学専攻都市環境システムコース岡野創 耐震工学の専門家として知っていた方が良いが 敷居が高く 入り口で挫折しがちな分野をいくつか取り上げて説明 ランダム振動論 地震波形に対する構造物応答の理論的把握 減衰と地震応答 エネルギーバランス 地震動の各種スペクトルの相互関係 震源モデル 近年では震源モデルによる地震動予測が良く行われている

More information

PowerPoint Presentation

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

More information

ベクトル公式.rtf

ベクトル公式.rtf 6 章ラプラシアン, ベクトル公式, 定理 6.1 ラプラシアン Laplacian φ はベクトル量である. そこでさらに発散をとると, φ はどういう形になるであろうか? φ = a + a + a φ a + a φ + a φ = φ + φ + φ = 2 φ + 2 φ 2 + 2 φ 2 2 φ = 2 φ 2 + 2 φ 2 + 2 φ 2 = 2 φ したがって,2 階の偏微分演算となる.

More information

untitled

untitled CFD JAXA CFD CFD Navier-Stokes -1- -2-1 CFD CFD Navier-Stokes 1. 2. 2,500,000 CFD Copyright Boeing CFD Boeing B777 HP -3- -4-2 CFD European Transonic Wind tunnel (-163 ) JAXA 2mx2m CFD (Computational Fluid

More information

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

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

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 米田 戸倉川月 7 限 193~21 西 5-19 応用数学 A 積分定理 Gaussの定理 divbd = B nds Stokesの定理 E bds = E dr Green の定理 g x f y dxdy = fdx + gdy = f e i + ge j dr Gauss の発散定理 S n FdS = Fd 1777-1855 ドイツ Johann arl Friedrich Gauss

More information

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

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

More information

116 宇宙航空研究開発機構特別資料 JAXA-SP cc ww1 = cc bb1 κκ 2 + (1 + cc bb2), σσ cc ww2 = 0.3, cc ww3 = 2, cc vv1 = 7.1, cc vv2 = 5 (4.7) 図 1 翼型モデルの外観 4. プリプ

116 宇宙航空研究開発機構特別資料 JAXA-SP cc ww1 = cc bb1 κκ 2 + (1 + cc bb2), σσ cc ww2 = 0.3, cc ww3 = 2, cc vv1 = 7.1, cc vv2 = 5 (4.7) 図 1 翼型モデルの外観 4. プリプ 第 45 回流体力学講演会 / 航空宇宙数値シミュレーション技術シンポジウム 2013 論文集 115 OpenFOAM を用いた NACA0012 翼型まわりの準二次元解析 中谷淳, 村澤杏樹岐阜工業高等専門学校 Quasi-2D Flow Analysis around NACA0012 Airfoil using OpenFOAM by Jun NAKAYA and Anju MURASAWA

More information

スライド 1

スライド 1 データ解析特論第 5 回 ( 全 15 回 ) 2012 年 10 月 30 日 ( 火 ) 情報エレクトロニクス専攻横田孝義 1 をもっとやります 2 第 2 回 3 データマイニングの分野ではマクロ ( 巨視的 ) な視点で全体を捉える能力が求められる 1. コンピュータは数値の集合として全体を把握していますので 意味ある情報として全体を見ることが不得意 2. 逆に人間には もともと空間的に全体像を捉える能力が得意

More information

2018年度 岡山大・理系数学

2018年度 岡山大・理系数学 08 岡山大学 ( 理系 ) 前期日程問題 解答解説のページへ 関数 f ( x) = ( + x) x について, 以下の問いに答えよ () f ( x ) = 0 を満たす x の値を求めよ () 曲線 y = f ( x ) について, 原点を通るすべての接線の方程式を求めよ (3) 曲線 y = f ( x ) について, 原点を通る接線のうち, 接点の x 座標が最大のものを L とする

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

Microsoft PowerPoint - GPUシンポジウム _d公開版.ppt [互換モード]

Microsoft PowerPoint - GPUシンポジウム _d公開版.ppt [互換モード] 200/0/9 数値流体解析の並列効率とその GPU による高速化の試み 清水建設 ( 株 ) 技術研究所 PHAM VAN PHUC ( ファムバンフック ) 流体計算時間短縮と GPU の活用の試み 現 CPUとの比較によりGPU 活用の可能性 現 CPU の最大利用 ノード内の最大計算資源の利用 すべてCPUコアの利用 適切なアルゴリズムの利用 CPU コア性能の何倍? GPU の利用の試み

More information

Microsoft Word - 11 進化ゲーム

Microsoft Word - 11 進化ゲーム . 進化ゲーム 0. ゲームの理論の分類 これまで授業で取り扱ってきたゲームは 協 ゲームと呼ばれるものである これはプレイヤー同士が独立して意思決定する状況を表すゲームであり ふつう ゲーム理論 といえば 非協力ゲームを表す これに対して プレイヤー同士が協力するという前提のもとに提携形成のパタンや利得配分の在り方を分析するゲームを協 ゲームという もっとも 社会現象への応用可能性も大きいはずなのに

More information

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

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

More information

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

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

More information

Matrix and summation convention Kronecker delta δ ij 1 = 0 ( i = j) ( i j) permutation symbol e ijk = (even permutation) (odd permutation) (othe

Matrix and summation convention Kronecker delta δ ij 1 = 0 ( i = j) ( i j) permutation symbol e ijk = (even permutation) (odd permutation) (othe Matr ad summato covto Krockr dlta δ ( ) ( ) prmutato symbol k (v prmutato) (odd prmutato) (othrs) gvalu dtrmat dt 6 k rst r s kt opyrght s rsrvd. No part of ths documt may b rproducd for proft. 行列 行 正方行列

More information

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

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

More information

スライド 1

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

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

数学 ⅡB < 公理 > 公理を論拠に定義を用いて定理を証明する 1 大小関係の公理 順序 (a > b, a = b, a > b 1 つ成立 a > b, b > c a > c 成立 ) 順序と演算 (a > b a + c > b + c (a > b, c > 0 ac > bc) 2 図

数学 ⅡB < 公理 > 公理を論拠に定義を用いて定理を証明する 1 大小関係の公理 順序 (a > b, a = b, a > b 1 つ成立 a > b, b > c a > c 成立 ) 順序と演算 (a > b a + c > b + c (a > b, c > 0 ac > bc) 2 図 数学 Ⅱ < 公理 > 公理を論拠に定義を用いて定理を証明する 大小関係の公理 順序 >, =, > つ成立 >, > > 成立 順序と演算 > + > + >, > > 図形の公理 平行線の性質 錯角 同位角 三角形の合同条件 三角形の合同相似 量の公理 角の大きさ 線分の長さ < 空間における座漂とベクトル > ベクトルの演算 和 差 実数倍については 文字の計算と同様 ベクトルの成分表示 平面ベクトル

More information

<4D F736F F F696E74202D CD C E688EA95945F CC8AEE91622E707074>

<4D F736F F F696E74202D CD C E688EA95945F CC8AEE91622E707074> CFD の基礎 九州大学応用力学研究所 内田孝紀 非定常 非線形風況 拡散シミュレータ RIAM-COMPACT Research Institute for Applied Mechanics, Kyushu University, COMputational Prediction of Airflow over Complex Terrain 数 ( 十 )km 以下の局所域局所域スケールに的を絞り,

More information

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

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

More information

<4D F736F F D20824F B834E835882CC92E8979D814690FC90CF95AA82C696CA90CF95AA2E646F63>

<4D F736F F D20824F B834E835882CC92E8979D814690FC90CF95AA82C696CA90CF95AA2E646F63> 1/10 平成 23 年 6 月 1 日午後 4 時 33 分 07 ストークスの定理 : 線積分と面積分 07 ストークスの定理 : 線積分と面積分 ストークスの定理はガウスの定理とともに 非常に重要な定理であり 線積分と面積分の関係を表します つまり ガウスの定理 : 面積分と体積分 ( 体積を囲む閉じた面 = 表面 ) の関係 ストークスの定理 : 線積分と面積分 ( 面積を囲む外周の線 )

More information

PowerPoint Presentation

PowerPoint Presentation Embedded CFD 1D-3D 連成によるエンジンコンパートメント熱収支解析手法の提案 June 9, 2017 . アジェンダ Embedded CFD 概要 エンコパ内風流れデモモデル 他用途への適用可能性, まとめ V サイクルにおける,1D-3D シミュレーションの使い分け ( 現状 ) 1D 機能的表現 企画 & 初期設計 詳細 3D 形状情報の無い段階 1D 1D 空気流れ計算精度に限度

More information

物性基礎

物性基礎 水素様原子 水素原子 水素様原子 エネルギー固有値 波動関数 主量子数 角運動量 方位量子数 磁気量子数 原子核 + 電子 個 F p F = V = 水素様原子 古典力学 水素様原子 量子力学 角運動量 L p F p L 運動方程式 d dt p = d d d p p = p + dt dt dt = p p = d dt L = 角運動量の保存則 ポテンシャルエネルギー V = 4πε =

More information

平成 年 月 7 日 ( 土 第 75 回数学教育実践研究会アスティ 45 ビル F セミナールーム A 札幌医科大学 年 P ab, を正の定数とする 平面上において ( a, を中心とする円 Q 4 C と (, b を中心とする円 C が 原点 O で外接している また P を円 C 上の点と

平成 年 月 7 日 ( 土 第 75 回数学教育実践研究会アスティ 45 ビル F セミナールーム A 札幌医科大学 年 P ab, を正の定数とする 平面上において ( a, を中心とする円 Q 4 C と (, b を中心とする円 C が 原点 O で外接している また P を円 C 上の点と 平成 年 月 7 日 ( 土 第 75 回数学教育実践研究会アスティ 45 ビル F セミナールーム 微分積分の拡張 変数関数問題へのアプローチ 予選決勝優勝法からラグランジュ未定乗数法 松本睦郎 ( 札幌北高等学校 変数関数の最大値 最小値に関する問題には多様なアプローチ法がある 文字を固定した 予選決勝優勝法, 計算のみで解法する 文字消去法, 微分積分を利用した ラグランジュ未定乗数法 がある

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 2009年11月7日 第2回オープンソースCAEワークショップ 1 オープンCAE DEXCS-OpenFOAMの紹介 株式会社デンソー 開発部 野村悦治 今川洋造 DEXCS http://dexcs.gifu-nct.ac.jp/ 2 オープンCAE DEXCS 2007/10より公開中 ADVENTUREによるお手軽なCAE環境です DEXCSの構成 3 Linux(Ubuntu)上に 構造解析にADVENTUREを活用し

More information

Microsoft Word - kogi10ex_main.docx

Microsoft Word - kogi10ex_main.docx 機能創造理工学 Ⅱ 期末試験 追試験問題 ( 病欠等による ) 途中の計算を必ず書こう 答えのみでは採点できない 問. 二次元面内を運動する調和振動子のラグランジアン L ( ) ( ) を 極座標, に変換し 極座標でのオイラーラグランジュ方程式を書こう ( 解く必要はない ) 但し, は定数であり また 極座標の定義は cos, sin である 問. 前問において極座標, に共役な一般化運動量,

More information

Xamテスト作成用テンプレート

Xamテスト作成用テンプレート 気体の性質 1 1990 年度本試験化学第 2 問 問 1 次の問い (a b) に答えよ a 一定質量の理想気体の温度を T 1 [K] または T 2 [K] に保ったまま, 圧力 P を変える このときの気体の体積 V[L] と圧力 P[atm] との関係を表すグラフとして, 最も適当なものを, 次の1~6のうちから一つ選べ ただし,T 1 >T 2 とする b 理想気体 1mol がある 圧力を

More information

Chap2.key

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( ) : π π

More information

2014年度 名古屋大・理系数学

2014年度 名古屋大・理系数学 04 名古屋大学 ( 理系 ) 前期日程問題 解答解説のページへ空間内にある半径 の球 ( 内部を含む ) を B とする 直線 と B が交わっており, その交わりは長さ の線分である () B の中心と との距離を求めよ () のまわりに B を 回転してできる立体の体積を求めよ 04 名古屋大学 ( 理系 ) 前期日程問題 解答解説のページへ 実数 t に対して 点 P( t, t ), Q(

More information