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

Size: px
Start display at page:

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

Transcription

1 乱流とは? 不規則運動であり, 速度の時空間的な変化が複雑であり, 個々の測定結果にはまったく再現性がなく, 偶然の値である. 渦運動 3 次元流れ 非定常流 乱流は確率過程 (Stochastic Process) である. 乱流工学 1 レイノルズの実験 UD = = ν 慣性力粘性力 乱流工学

2 F レイノルズ数 U L / U 3 = mα = ρl = ρ 慣性力 L U u U A = µ A = µ L = µ L U y L τ.. せん断応力 面積 レイノルズ数 = 慣性力粘性力 = ρlu µ L U = UL UL = µ / ρ ν 乱流工学 3 r Du Dt ナビエ ストークス方程式 = gradp + µ ρ r u 慣性力 = 圧力 + 粘性力 乱流工学 4

3 ナビエ ストークス方程式 * u * p * t * u =, p =, t =, x = U ρu L / U x L r Du Dt * * = gradp * + 大 : 慣性力 = 圧力小 : 粘性力 = 圧力 1 乱流工学 5 r u * 代表速度 乱流の場合, 代表速度は主流の速度ではなく, どれだけ流れが変動しているかを考える必要がある 速度変動の 乗平均値は乱流のもつエネルギーの大きさを表す尺度である. 速度の単位にするにはその平方根をとればよい 代表速度 = u' 一般に実効値を平均速度で無次元化し, 平均値に対する割合で表すことが多い. 論文等では乱流強度, 乱れ強さ (Turbulent intensity) と表現される. 乱流工学 6

4 代表寸法 代表寸法はどのようにとるべきか? 乱流工学 7 どちらのレイノルズ数も 同じである 乱流工学 8

5 エネルギースペクトル 乱流発生装置を用いて 乱流場の統計的性質を 解明する 乱流工学 9 積分特性距離 Rˆ ij ( r) = u i( x) u j( x) [ u i ( x) u j ( x ) ] 1 / L 1 ij R( r) i j = dr u ( x) u ( x) 1. Rij 1 L f(r) Integeral scale r 乱流工学 1

6 乱流工学 11 乱流レイノルズ数 1 = x u u λ ν λ = λ ' u ν = L u L ' 積分特性距離を代表寸法 t U x c = 1 テイラーのマイクロスケールテイラーの凍結性仮説微分特性距離を代表寸法乱流工学 1 乱流発生装置 7 5

7 乱流レイノルズ数を用いた応用例 トラックなどの大型車両は走行抵抗に占める空気抵抗の割合が高い 空気抵抗の軽減.7 風洞実験および数値解析による車両空力特性の評価 実車 8 km/h でのレイノルズ数 燃費の向上 CO 排出量の削減 1 5 模型実験 小型風洞中で大型風洞と同等の遷移流れを模擬する 乱流工学 13 風洞設備 大型風洞財団法人鉄道総合技術研究所米原風洞技術センター ノズル寸法 :5 m 3 m 最大流速 :83.3 m/s 非一様性 : % 以内 乱れ強度 :. % 以下 小型風洞工学院大学 =5 1 6 ノズル寸法 :.3 m.3 m 最大流速 :3 m/s 非一様性 :1 % 以内 乱れ強度 :1 % 以下 =3 1 5 乱流工学 14

8 トラックキャビンモデル H キャビンルーフ前端部 W L M 曲率半径 R 縮尺比 W L M H / / / 模型サイズ [m] キャビンルーフ前端部の R 寸法 [mm] 縮尺比 R R1 R 1 1 1/ / / R 形状による空気抵抗依存性を調べる 乱流工学 15 乱流発生装置 ステッピングモータ 矩形翼 3 矩形翼を振動させることにより気流中に渦を励起して乱れを生成 Integlal Scale [mm] 1 乱流格子 乱流発生装置 Turbulent Intensity [%] 乱流発生装置と乱流格子を用いて広範囲の乱れ強度, 渦スケールの流れ場で実験を行った 乱流工学 16

9 大型風洞実験結果.9.7 R /H =.7 R /H = 代表長さ : 模型全長 L M 代表速度 : 主流速度 U としたレイノルズ数を本研究では幾何学的レイノルズ数 と定義 ルーフ前端部の R を大きくすることで 値が低下 ルーフ前端部の R を大きくすることで crit が低下 (R/H=.7 : crit =. 1 6,R/H=.14 : crit = ) 乱流工学 17 小型風洞実験結果 小型風洞 : 乱流格子および乱流発生装置非装着時 1. R /H = Swt R /H = Lwt R /H =.7 Lwt R /H =.7 Swt Swt Lwt 小型風洞の実験結果は臨界レイノルズ数以下 乱流工学 18

10 抗力係数に対する主流乱れの影響.9.7 = R /H = R /H =.7 R /H =3.9.7 L /L M=.184 L /L M=.65 L /L M=.153 = u' /U [%] R /L 主流の乱れ強度を増加することで 値が低下 ルーフ前端部の R が大きくなるにつれ 値が低下 R/L が 1 に近くなるにつれ 値が低下する傾向がある 乱流工学 19 キャビンモデル表面の流れ場 (PIV) 小 乱れ小 大 乱れ大 トラックキャビン周りの流れ (PIV) 乱流工学

11 キャビンモデル表面の流れ場 Z 計測範囲 Ū /U : u /U =% X R/H=.14 Ū /U : 主流乱れを増加 u /U =3.63% 模型表面がはく離流れから付着流れへ遷移する傾向がある 値が低下 Ū /U : u /U =1% 乱流工学 1 A 主流乱れとレイノルズ数効果 R/H =.14 R/H =.14 Ū /U : A Ū /U : = ,u /U =% = ,u /U =1% A A = ,u /U =% = ,u /U =.8% 乱流工学

12 レイノルズ数推定モデル式 幾何学的レイノルズ数 に対して実効レイノルズ数 e を以下のように定義 e = {1+ f (u /U ) g(r/l)} 鉄道総合技術研究所の大型風洞は主流の乱れ強度, 渦スケールが十分小さい e = L = S {1+ f (u /U ) g(r/l)} L : 大型風洞の幾何学的 S : 小型風洞の幾何学的 乱流工学 関数 f および関数 g の算出 u' /U [%] R /H = R /H =.7 R /H =.7 R /H = S= S= R /H = L R /H = 主流に乱れを与えた場合の 値と大型風洞の 値が一致する点を求める R/L 乱流工学 4

13 L L 関数 f および関数 g の算出 u' /U [%] R /H =.7 R /H =.7 R /H = S = R /H = g (R /L ) = 11.(R /L ) -.9(R /L ) R /L S = ( L- S)/ S ( L- S)/ S f (u'/u ) = -.35(u'/U ) + 4(u'/U ) u' /U [%] R /H =.7 R /H = R /H =.7 R /H = S = S = R /L 乱流工学 5 実効レイノルズ数と抗力係数の関係.9.7 Swt Exp. L wt Exp. R/H=.9.7 Swt Exp. R/H=.7 L wt Exp Swt Exp. L wt Exp. 1 6 e e R/H= e プロット塗り潰し : 小型風洞の実験結果プロット白抜き : 大型風洞の実験結果 小型風洞の実験から大型風洞における高レイノルズ数の流れ場を模擬できる 乱流工学 6

14 まとめ 乱流発生装置を備えた小型風洞を用いて, トラックなどの大型車両の空気抵抗を評価する手法について検討し以下の知見を得た 小型風洞で乱流発生装置を用いることで, 臨界レイノルズ数近傍におけるトラックキャビンの空気抵抗を推定することができる 主流乱れは実効レイノルズ数を増加させる因子である可能性が高い 主流の渦スケール L とルーフ前端部の R の比が 1 に近い場合 値が小さくなる傾向があることがわかった. 主流乱れを増加することで, 模型表面がはく離流れから付着流れに遷移し 値が低下する 乱流工学 7 境界層と模型設置位置の関係 1 小型風洞 : 乱流格子, 乱流発生装置非装着時 1. R /H = Swt R /H =.7 Swt Swt R /H = Lwt Lwt R /H =.7 Lwt Z /δ Z /δ = u' /U =% U =3m/s u' /U =.% 1/1.6 縮尺模型 Ū /U 1/3.6 縮尺模型 1/7. 縮尺模型 Ū /U 乱流工学 8

15 境界層と模型設置位置の関係 小型風洞 : 乱流格子, 乱流発生装置非装着時 1. Lwt R /H = Swt R /H =.7 Swt Swt R /H = Lwt R /H =.7 Lwt Z /δ Z /δ = u' /U =% Ū /U U =3m/s u' /U =.% 1/1.6 縮尺模型 1/3.6 縮尺模型 1/7. 縮尺模型 Ū /U 乱流工学 9 模型設置位置による流体力の変化.9 地上高 mm R /H =.9 地上高 mm R /H = 地上高 mm = u' /U [%] 地上高 mm = u' /U [%] 地上高 : 模型下面と地面板との距離.7 地上高 mm 地上高 mm = 模型下面と地面板との距離が離れることで 値が低下 u' /U [%] 乱流工学 3

16 乱れ強度と大型風洞の幾何学的レイノルズ数の関係.9 R /H = R /H = R /H =.7.7 L S = u' /U [%] R /H =.7 R /H = 1 5 S = u' /U [%].7 f (u /U ) = L - S S L 乱流工学 31 渦スケールと大型風洞の幾何学的レイノルズ数の関係 R /H = R /H =.7 R /H =.7.7 L S = R /L R /H = S = R /L R /H =.7 L - S g (R /L) = S L 乱流工学 3

17 ( L- S)/ S ( s- L)/ L 実効レイノルズ数と抗力係数の関係 f (u'/u )= -.8(u'/U ) +.65(u'/U ) R /H =.7 S = u' /U [%] g (R /L ) = -3.4(R /L ) + 1.1(R /L ) R /H =.7 S = 各模型の値を算術平均し関数 f および関数 g を算出 e = {1+ f (u /U ) g(r/l)} S = e R/L 乱流工学 R /H = R /H =.7 大型風洞試験結果との比較.9 L wt Exp. R/H=.9 R/H=.7 L wt Exp..7 Swt Exp..7 Swt Exp e e.9 R/H=.14.7 Swt Exp. 1 5 L wt Exp 乱流工学 e 34

18 流体力計測結果 R /H = u' /U [%] R /H = R /H =.7 R /H =.7 S = R /H =.7 R /H = L s = R /L 乱流工学 35 レイノルズ数推定モデル式の要素 ( ) L R /H =.7 s= u' /U [%] f u / U ) ( ( L- S)/ S f (u' /U ) =.15(u' /U ) - 93(u' /U ) + 93 R /H =.7 s = u' /U [%] L 1 6 R /H = R /L s = g( R / L) ( L- S)/ S 1 5 g (R /L ) = -45.3(R /L ) + 3(R /L ) + R /H =.7 s = R /L 乱流工学 36

19 実効レイノルズ数と抗力係数の関係.9 Swt Exp. R/H=.9 R/H=.7.7 L wt Exp..7 Swt Exp. L wt Exp e e R/H=.14 Swt Exp L wt Exp 乱流工学 e 37 表面圧力計測結果 -1 R /H =.7 C P min s = u' /U [%] C P min R /H =.7 R /H = -5 C P min R /H = L -5 s = R/L 乱流工学 38

20 レイノルズ数推定モデル式の要素 (C P ) L 1 6 R /H = u' /U [%] s= f u / U ) ( ( L- S)/ S f (u'/u ) = -.8(u'/U ) +.6(u'/U ) R /H =.7 s = u' /U [%] L R /H = R /L s= g( R / L) ( L- S)/ S 1 5 g (R /L ) =.6(R /L ) (R /L ) R /H =.7 s = R /L 乱流工学 39 実効レイノルズ数と圧力係数の関係 -1 R/H=.7-1 R/H=.14 - Swt Exp. - C p min -3-4 L wt Exp. C p min -3-4 L wt Exp. Swt Exp e e 乱流工学 4

21 15 測定部における境界層付近の流速分布 u' /U =1% Z [mm] 1 5 u' /U =%. 1. Ū /U Y Z X 1. u' /U =% Z /δ. u' /U =1%.. 1. Ū /U 乱流工学 41 地面板の圧力場.1.5 フラップ調整後 C p. -.5 フラップ調整前 X [mm] 乱流工学 4

22 大型風洞の測定部における流速分布 Z [mm] Uo=3m/s Uo=55m/s Uo=8m/s Z /δ 1. Uo=3m/s Uo=55m/s Uo=8m/s Z [mm] BLS On Ū /U Uo=3m/s Uo=55m/s Uo=8m/s Z /δ.. 1. BLS On Ū /U Uo=3m/s Uo=55m/s Uo=8m/s 5. BLS Off BLS Off Ū /U Ū /U 乱流工学 43 流体力および表面圧力に対する境界層の影響.9 BLS On C P min -3-4 BLS Off 1 5 BLS Off BLS On 乱流工学 44

23 キャビンモデル周りの流れ場の比較 Ū /U : Ū /U : u /U =% Ū /U : Ū /U : u /U =3.63% Ū /U : Ū /U : u /U =1% R/H=.7 R/H=.14 乱流工学 45 抗力係数に対する主流乱れの影響.9 R /H = R /H =.7.7 = u' /U [%] 乱流工学 46

24 PIV 計測結果との比較 A Ū /U A : Ū /U : = ,u /U =% = ,u /U =1% A A = ,u /U =% = ,u /U =.8% 乱流工学 47 キャビンモデル表面の境界層厚さの比較 u'/u =3.63% δ /H δ /H u'/u =% u'/u =1% u'/u =3.63% R /H =.7 = X /H u'/u =% δ /H δ max/h.. u'/u =.8% u'/u =.8% u'/u =1%.1.1 = X /H 乱流工学 u'/u =% u'/u =1% = X /H u'/u =% u'/u =3.63% u'/u =%

25 曲率半径 R に対する渦スケールの影響 C d.9.7 R/H= R/H=.7 R/H=.9 R/H=.11 R/H=.13 R/H=.14 R/H=.16 R/H=.18 R/H=.1 R/H=.9 R/H=.36 R/H=3 = R /L C p min R/H=.7 R/H=.9 R/H=.11 R/H=.13 R/H=.14 R/H=.16 R/H=.18 R/H=.1 R/H=.9 R/H=.36 R/H=3 = R /L 乱流工学 49 乱流発生装置の速度分布 Z X Y 風軸 u'/u [%] U =3 m/s Laminar 1.1 U =3 m/s GT ATG_1deg. Ū /U.9 Laminar GT ATG_1deg X /D Z 乱流工学 /D 5

技術資料 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

車体まわり非定常流れの制御による空気抵抗低減技術の開発 プロジェクト責任者 加藤千幸 国立大学法人東京大学生産技術研究所 著者加藤千幸 * 1 鈴木康方 * 2 前田和宏 * 3 槇原孝文 * 3 北村任宏 * 3 高山務 * 4 廣川雄一 * 5 西川憲明 * 5 * 1 国立大学法人東京大学生産

車体まわり非定常流れの制御による空気抵抗低減技術の開発 プロジェクト責任者 加藤千幸 国立大学法人東京大学生産技術研究所 著者加藤千幸 * 1 鈴木康方 * 2 前田和宏 * 3 槇原孝文 * 3 北村任宏 * 3 高山務 * 4 廣川雄一 * 5 西川憲明 * 5 * 1 国立大学法人東京大学生産 車体まわり非定常流れの制御による空気抵抗低減技術の開発 プロジェクト責任者 加藤千幸 国立大学法人東京大学生産技術研究所 著者加藤千幸 * 1 鈴木康方 * 2 前田和宏 * 3 槇原孝文 * 3 北村任宏 * 3 高山務 * 4 廣川雄一 * 5 西川憲明 * 5 * 1 国立大学法人東京大学生産技術研究所 * 2 日本大学理工学部機械工学科 * 3 トヨタ自動車株式会社 * 4 みずほ情報総研株式会社

More information

第 3 章二相流の圧力損失

第 3 章二相流の圧力損失 第 3 章二相流の圧力損失 単相流の圧力損失 圧力損失 (/) 壁面せん断応力 τ W 力のバランス P+ u m πd 4 τ w 4 τ D u τ w m w πd : 摩擦係数 λ : 円管の摩擦係数 λ D u m D P τ W 摩擦係数 層流 16/Re 乱流 0.079 Re -1/4 0.046 Re -0.0 (Blasius) (Colburn) 大まかには 0.005 二相流の圧力損失液相のみが流れた場合の単相流の圧力損失

More information

1 抗力 揚力の計測 Ⅰ 18 年度用 はじめに 機械応用実験であることから, 意図的に親切なテキストとはしていない. 説明を良く聞き, 自分で考え, 実験を進めること. また, レポートには 1. 目的,. 実験方法,3. 結果,4. 考察,5. 検討 等を記すこと. このため, 実験を進めながらメモを残してゆき, このメモを基にしてまとめることが必要となる. なお, この実験の HP(http://www.cce.kanagawa-it.ac.jp/~t514/experiment/index.html)

More information

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

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

More information

研究の背景これまで, アルペンスキー競技の競技者にかかる空気抵抗 ( 抗力 ) に関する研究では, 実際のレーサーを対象に実験風洞 (Wind tunnel) を用いて, 滑走フォームと空気抵抗の関係や, スーツを含むスキー用具のデザインが検討されてきました. しかし, 風洞を用いた実験では, レー

研究の背景これまで, アルペンスキー競技の競技者にかかる空気抵抗 ( 抗力 ) に関する研究では, 実際のレーサーを対象に実験風洞 (Wind tunnel) を用いて, 滑走フォームと空気抵抗の関係や, スーツを含むスキー用具のデザインが検討されてきました. しかし, 風洞を用いた実験では, レー 報道関係者各位 平成 29 年 1 月 6 日 国立大学法人筑波大学 アルペンスキー競技ダウンヒルにおいてレーサーが受ける空気抵抗は下腿部が最大 ~ 身体部位ごとの空力特性を初めて解明 ~ 研究成果のポイント 1. アルペンスキー競技ダウンヒルにおける レーサーの身体全体と, 各身体部分の空気抵抗 ( 抗力 ) を, 世界に先駆けて明らかにしました. 2. 風洞実験と数値流体解析の結果, クラウチング姿勢におけるレーサー身体各部位の抵抗の大きさは,

More information

Microsoft PowerPoint - 第3回MSBS研究会.pptx

Microsoft PowerPoint - 第3回MSBS研究会.pptx 2013 年 3 月 1 日第 3 回 MSBS 研究会 アーチェリー矢の空力特性 MSBS 風洞実験と飛翔実験 電気通信大学大学院宮嵜武 JAXA 杉浦裕樹 円柱境界層 理論解 ( 境界層近似 ): 円柱側面の境界層 ( べき級数解 ) Seban & Bond (1951) J. Aero. Sci. 18 先端部べき級数解 Kelly (1954) J. Aero. Sci. 21 修正版べき級数解

More information

報道関係者各位 平成 26 年 5 月 29 日 国立大学法人筑波大学 サッカーワールドカップブラジル大会公式球 ブラズーカ の秘密を科学的に解明 ~ ボールのパネル構成が空力特性や飛翔軌道を左右する ~ 研究成果のポイント 1. 現代サッカーボールのパネルの枚数 形状 向きと空力特性や飛翔軌道との

報道関係者各位 平成 26 年 5 月 29 日 国立大学法人筑波大学 サッカーワールドカップブラジル大会公式球 ブラズーカ の秘密を科学的に解明 ~ ボールのパネル構成が空力特性や飛翔軌道を左右する ~ 研究成果のポイント 1. 現代サッカーボールのパネルの枚数 形状 向きと空力特性や飛翔軌道との 報道関係者各位 平成 26 年 5 月 29 日 国立大学法人筑波大学 サッカーワールドカップブラジル大会公式球 ブラズーカ の秘密を科学的に解明 ~ ボールのパネル構成が空力特性や飛翔軌道を左右する ~ 研究成果のポイント 1. 現代サッカーボールのパネルの枚数 形状 向きと空力特性や飛翔軌道との関係を明らかにしました 2. 風洞実験の結果 ブラズーカ ( ワールドカップ 2014 公式球 ) は

More information

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

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

More information

伝熱学課題

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

More information

粒子画像流速測定法を用いた室内流速測定法に関する研究

粒子画像流速測定法を用いた室内流速測定法に関する研究 可視化手法を用いた室内気流分布の測定法に関する研究 -PIV を用いた通風時及び空調吹出気流の測定 - T08K729D 大久保肇 指導教員 赤林伸一教授 流れの可視化は古来より流れの特性を直感的に把握する手法として様々な測定法が試みられている 近年の画像処理技術の発展及び PC の性能向上により粒子画像流速測定法 (PIV ) が実用化されている Particle Image Velocimetry

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 反応工学 Raction Enginring 講義時間 ( 場所 : 火曜 限 (8-A 木曜 限 (S-A 担当 : 山村 火 限 8-A 期末試験中間試験以降 /7( 木 まで持ち込みなし要電卓 /4( 木 質問受付日講義なし 授業アンケート (li campus の入力をお願いします 晶析 (crystallization ( 教科書 p. 濃度 溶解度曲線 C C s A 安定 液 ( 気

More information

<4D F736F F D2091E6358FCD31328B438FDB A5182F08ADC82DE816A2E646F6378>

<4D F736F F D2091E6358FCD31328B438FDB A5182F08ADC82DE816A2E646F6378> 注 : 大阪管区気象台における平成 15 年から平成 24 年の観測データを元に作成 図 5-12-3 日最大平均風速の風向出現頻度 346 5.12. 2 施設の存在に伴う影響の予測 評価 (1) 予測内容施設の存在に伴う影響として 建築物の出現が事業計画地周辺の風環境に及ぼす影響について 風洞実験により予測した 予測内容は表 5-12-2 に示すとおりである 表 5-12-2 予測内容 予測項目予測範囲

More information

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

オープン CAE 関東 数値流体力学 輪講 第 4 回 第 3 章 : 乱流とそのモデリング (3) [3.5~3.7.1 p.64~75] 日時 :2013 年 11 月 10 日 14:00~ 場所 : 日本 新宿 2013/11/10 数値流体力学 輪講第 4 回 1 オープン CAE 勉強会 @ 関東 数値流体力学 輪講 第 4 回 第 3 章 : 乱流とそのモデリング (3 [3.5~3.7.1 p.64~75] 日時 :2013 年 11 月 10 日 14:00~ 場所 : 日本 ESI@ 新宿 1 数値流体力学 輪講に関して 目的 数値流体力学の知識 ( 特に理論ベース を深め OpenFOAM の利用に役立てること 本輪講で学ぶもの 数値流体力学の理論や計算手法の概要

More information

64 3 g=9.85 m/s 2 g=9.791 m/s 2 36, km ( ) 1 () 2 () m/s : : a) b) kg/m kg/m k

64 3 g=9.85 m/s 2 g=9.791 m/s 2 36, km ( ) 1 () 2 () m/s : : a) b) kg/m kg/m k 63 3 Section 3.1 g 3.1 3.1: : 64 3 g=9.85 m/s 2 g=9.791 m/s 2 36, km ( ) 1 () 2 () 3 9.8 m/s 2 3.2 3.2: : a) b) 5 15 4 1 1. 1 3 14. 1 3 kg/m 3 2 3.3 1 3 5.8 1 3 kg/m 3 3 2.65 1 3 kg/m 3 4 6 m 3.1. 65 5

More information

Microsoft PowerPoint - siryo7

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

More information

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

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

More information

Sample function Re random process Flutter, Galloping, etc. ensemble (mean value) N 1 µ = lim xk( t1) N k = 1 N autocorrelation function N 1 R( t1, t1

Sample function Re random process Flutter, Galloping, etc. ensemble (mean value) N 1 µ = lim xk( t1) N k = 1 N autocorrelation function N 1 R( t1, t1 Sample function Re random process Flutter, Galloping, etc. ensemble (mean value) µ = lim xk( k = autocorrelation function R( t, t + τ) = lim ( ) ( + τ) xk t xk t k = V p o o R p o, o V S M R realization

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

untitled

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

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

講義「○○○○」

講義「○○○○」 講義 信頼度の推定と立証 内容. 点推定と区間推定. 指数分布の点推定 区間推定 3. 指数分布 正規分布の信頼度推定 担当 : 倉敷哲生 ( ビジネスエンジニアリング専攻 ) 統計的推測 標本から得られる情報を基に 母集団に関する結論の導出が目的 測定値 x x x 3 : x 母集団 (populaio) 母集団の特性値 統計的推測 標本 (sample) 標本の特性値 分布のパラメータ ( 母数

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

オープン 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

II ( ) (7/31) II ( [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Re

II ( ) (7/31) II (  [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Re II 29 7 29-7-27 ( ) (7/31) II (http://www.damp.tottori-u.ac.jp/~ooshida/edu/fluid/) [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Reynolds [ (4.6), (45.8)] [ p.186] Navier Stokes I Euler Navier

More information

1

1 新幹線車両の空力騒音シミュレーション プロジェクト責任者 栗田健 東日本旅客鉄道株式会社 JR 東日本研究開発センター 著者水島文夫 *1 栗田健 *1 山出吉伸 *2 加藤千幸 *2 上原均 *3 *3 廣川雄一 *1 東日本旅客鉄道株式会社 JR 東日本研究開発センター *2 国立大学法人東京大学生産技術研究所 *3 独立行政法人海洋研究開発機構計算システム計画 運用部 新幹線の車間部から発生する空力騒音の発生メカニズムを解明するため

More information

20年度一次基礎略解

20年度一次基礎略解 年度一次機械問題略解 計算問題中心 orih c 0 宮田明則技術士事務所 正解番号 Ⅳ- Ⅳ- Ⅳ- Ⅳ- Ⅳ- Ⅳ-6 Ⅳ-7 Ⅳ-8 Ⅳ-9 Ⅳ-0 Ⅳ- Ⅳ- Ⅳ- Ⅳ- Ⅳ- Ⅳ-6 Ⅳ-7 Ⅳ-8 orih c 0 宮田明則技術士事務所 Ⅳ-9 Ⅳ-0 Ⅳ- Ⅳ- Ⅳ- Ⅳ- Ⅳ- Ⅳ-6 Ⅳ-7 Ⅳ-8 Ⅳ-9 Ⅳ-0 Ⅳ- Ⅳ- Ⅳ- Ⅳ- Ⅳ- 特定入力関数と応答の対応の組み合わせフィードバック制御に関する記述の正誤正弦波入力に対する定常出力の計算フィードバック系の特性根を求める計算比熱等に関する

More information

Untitled

Untitled II 14 14-7-8 8/4 II (http://www.damp.tottori-u.ac.jp/~ooshida/edu/fluid/) [ (3.4)] Navier Stokes [ 6/ ] Navier Stokes 3 [ ] Reynolds [ (4.6), (45.8)] [ p.186] Navier Stokes I 1 balance law t (ρv i )+ j

More information

大型輸送車両の空気抵抗低減の研究

大型輸送車両の空気抵抗低減の研究 平成 25 年度 修士論文 大型輸送車両の空気抵抗低減装置の実験及び開発 指導教員 筒井康賢教授 高知工科大学 知能機械システム工学科 筒井研究室 学生番号 1165045 矢野佑樹 目次 第 1 章序論 1.1 研究背景 1.1.1 走行抵抗 1.1.2 境界層とレイノルズ数 Re 1.2 研究目的 第 2 章従来の研究 第 3 章実験装置 3.1 風洞 3.2 計測機器 3.2.1 差圧計 DSA3217/16Px

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

構造力学Ⅰ第12回

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

More information

スライド 1

スライド 1 移動体観測を活用した交通 NW の リアルタイムマネジメントに向けて : プローブカーデータを用いた動的 OD 交通量のリアルタイム推定 名古屋大学山本俊行 背景 : マルチモード経路案内システム PRONAVI 2 プローブカーデータの概要 プローブカー : タクシー 157 台 蓄積用データ収集期間 : 22 年 1 月 ~3 月,1 月 ~23 年 3 月 データ送信はイベントベース : 車両発進

More information

Probit , Mixed logit

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

More information

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

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

More information

資料

資料 平成 28 年度地球シミュレータ利用報告会 2017/4/20 コクヨホール 研究分野 : 計算科学分野 高精度流体シミュレーションによる小型ファンまわりの流体現象理解と空力騒音発生メカニズムの解析 : 騒音 という環境課題の改善にむけて 青野 高橋 浅田 関本 立川 藤井 東京理科大学阿部 東京工業大学村上 ミネベア株式会社板倉 廣川 JAMSTEC 研究概要 小型ファンから発生する空力騒音を近年

More information

Microsoft Word - 補論3.2

Microsoft Word - 補論3.2 補論 3. 多変量 GARC モデル 07//6 新谷元嗣 藪友良 対数尤度関数 3 章 7 節では 変量の対数尤度を求めた ここでは多変量の場合 とくに 変量について対数尤度を求める 誤差項 は平均 0 で 次元の正規分布に従うとする 単純化のため 分散と共分散は時間を通じて一定としよう ( この仮定は後で変更される ) したがって ij から添え字 を除くことができる このとき と の尤度関数は

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

Techniques for Nuclear and Particle Physics Experiments Energy Loss by Radiation : Bremsstrahlung 制動放射によるエネルギー損失は σ r 2 e = (e 2 mc 2 ) 2 で表される為

Techniques for Nuclear and Particle Physics Experiments Energy Loss by Radiation : Bremsstrahlung 制動放射によるエネルギー損失は σ r 2 e = (e 2 mc 2 ) 2 で表される為 Techniques for Nuclear and Particle Physics Experiments.. Energy Loss by Radiation : Bremsstrahlung 制動放射によるエネルギー損失は σ r e = (e mc ) で表される為 質量に大きく依存する Ex) 電子の次に質量の小さいミューオンの制動放射によるエネルギー損失 m e 0.5 MeV, m

More information

Microsoft PowerPoint - 01_内田 先生.pptx

Microsoft PowerPoint - 01_内田 先生.pptx 平成 24 年度 SCOPE 研究開発助成成果報告会 ( 平成 22 年度採択 ) 塩害劣化した RC スラブの一例 非破壊評価を援用した港湾コンクリート構造物の塩害劣化予測手法の開発 かぶりコンクリートのはく落 大阪大学大学院鎌田敏郎佐賀大学大学院 内田慎哉 の腐食によりコンクリート表面に発生したひび割れ ( 腐食ひび割れ ) コンクリート構造物の合理的な維持管理 ( 理想 ) 開発した手法 点検

More information

Microsoft Word - planck定数.doc

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

More information

流束の大きさは濃度勾配に比例すると見なせ ( フィックの法則 ) その比例係数 D を拡散係 数と呼ぶ J = D 拡散定数は [ 面積 ]/[ 時間 ] の次元を持つ ある地点の濃度の変化に注目すると 化学反応など が起きなければ 濃度変化は流束の変化に等しく 次の偏微分方程式が成立する ( 拡散

流束の大きさは濃度勾配に比例すると見なせ ( フィックの法則 ) その比例係数 D を拡散係 数と呼ぶ J = D 拡散定数は [ 面積 ]/[ 時間 ] の次元を持つ ある地点の濃度の変化に注目すると 化学反応など が起きなければ 濃度変化は流束の変化に等しく 次の偏微分方程式が成立する ( 拡散 化学実験法 II 2014.5.15 吉村洋介 内容 拡散と混合の話...1 物質輸送の構成 : 対流と拡散...1 フィックの法則と拡散方程式...1 定常的な拡散...2 非定常な拡散...2 拡散の一般的挙動と拡散定数...3 対流と分散による拡散 混合...4 管中の流れの中の拡散 (Taylor 分散 )...4 乱流中の混合 拡散...4 化学反応 相互作用をともなう場合の拡散...5

More information

A

A A04-164 2008 2 13 1 4 1.1.......................................... 4 1.2..................................... 4 1.3..................................... 4 1.4..................................... 5 2

More information

<4D F736F F F696E74202D208BAB8A458FF08C8F82CC8AEE916282C68C8892E896402E707074>

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

More information

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

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

More information

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

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

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

PRESS RELEASE (2016/11/29) 九州大学広報室 福岡市西区元岡 744 TEL: FAX: URL:

PRESS RELEASE (2016/11/29) 九州大学広報室 福岡市西区元岡 744 TEL: FAX: URL: PRESS RELEASE (2016/11/29) 九州大学広報室 819-0395 福岡市西区元岡 744 TEL:092-802-2130 FAX:092-802-2139 MAIL:koho@jimu.kyushu-u.ac.jp URL:http://www.kyushu-u.ac.jp 地形起因の大気乱流が大型風車の構造強度に与える影響の評価に成功 - 大規模風力発電の普及 拡大に期待

More information

Microsoft Word 卒業論文2.doc

Microsoft Word 卒業論文2.doc 平成 6 年度 卒業論文 狭窄部を有する血管内の血流の有限要素解析 高知工科大学工学部知能機械システム工学科知能流体力学研究室 清水昌彦 目次 第 章緒言 - 本研究を行う背景と目的 - 血液の性質 -3 数値計算 - 有限要素法の概要 第 章基礎方程式 - 支配方程式 -- 連続の式 5 -- コーシーの運動方程式 6 --3 血液の構成方程式 6 - 無次元化 7 第 3 章解析手法 3- 有限要素解析

More information

Microsoft PowerPoint - 12_2019裖置工�榇諌

Microsoft PowerPoint - 12_2019裖置工å�¦æ¦‡è«Œ 1 装置工学概論 第 12 回 蒸留装置の設計 (3) 流動装置の設計 (1) 東京工業大学物質理工学院応用化学系 下山裕介 2019.7.15 装置工学概論 2 第 1 回 4 /15 ガイダンス : 化学プロセスと装置設計 第 2 回 4 /22 物質 エネルギー収支 第 3 回 5 /6( 祝 ) 化学プロセスと操作変数 5 /13 休講 第 4 回 5 /20 無次元数と次元解析 第 5 回

More information

研究成果報告書

研究成果報告書 様式 C-19 科学研究費補助金研究成果報告書 研究種目 : 若手研究 ( スタートアップ ) 研究期間 :7~8 課題番号 :19861 研究課題名 ( 和文 ) 大規模火災旋風の性状予測および被害評価 平成 1 年 5 月 9 日現在 研究課題名 ( 英文 ) Understandingthenatureof,andevaluatingthedamagecaused bylarge-scalefirewhirls

More information

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

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

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

平成 23 年度 JAXA 航空プログラム公募型研究報告会資料集 (23 年度採用分 ) 21 計測ひずみによる CFRP 翼構造の荷重 応力同定と損傷モニタリング 東北大学福永久雄 ひずみ応答の計測データ 静的分布荷重同定動的分布荷重同定 ひずみゲージ応力 ひずみ分布の予測 or PZT センサ損

平成 23 年度 JAXA 航空プログラム公募型研究報告会資料集 (23 年度採用分 ) 21 計測ひずみによる CFRP 翼構造の荷重 応力同定と損傷モニタリング 東北大学福永久雄 ひずみ応答の計測データ 静的分布荷重同定動的分布荷重同定 ひずみゲージ応力 ひずみ分布の予測 or PZT センサ損 平成 3 年度 JAXA 航空プログラム公募型研究報告会資料集 (3 年度採用分 1 計測ひずみによる CFRP 翼構造の荷重 応力同定と損傷モニタリング 東北大学福永久雄 ひずみ応答の計測データ 静的分布荷重同定動的分布荷重同定 ひずみゲージ応力 ひずみ分布の予測 or PZT センサ損傷発生位置の推定発表内容 (1 荷重同定 1:11 点衝撃荷重同定 ( 荷重同定 : 分布荷重同定 (3 今後の予定

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

gr09.dvi

gr09.dvi .1, θ, ϕ d = A, t dt + B, t dtd + C, t d + D, t dθ +in θdϕ.1.1 t { = f1,t t = f,t { D, t = B, t =.1. t A, tdt e φ,t dt, C, td e λ,t d.1.3,t, t d = e φ,t dt + e λ,t d + dθ +in θdϕ.1.4 { = f1,t t = f,t {

More information

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

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

More information

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

第 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 は境界条件と振動モードによって決まる無

More information

水の粘度 (mpa s) y R φd p τ 図 3. 円管内流速分布の解析モデル 応力を τ とすると 円筒の全側面に作用するせん 断応力による力は となる 定常状態 では この力と圧力による力 () 式が釣り合うので これから τ を求めると () (3) 壁から円筒側面までの距離を y とす

水の粘度 (mpa s) y R φd p τ 図 3. 円管内流速分布の解析モデル 応力を τ とすると 円筒の全側面に作用するせん 断応力による力は となる 定常状態 では この力と圧力による力 () 式が釣り合うので これから τ を求めると () (3) 壁から円筒側面までの距離を y とす 圧損型粘度計による日用品流体の粘度測定 望月聡 小椋勝仁 山田重良 桜木俊一 概要 : 私達の日常生活で利用される 醤油やソース ドレッシングやマヨネーズなどの食品流体 さらに シャンプーやリンス 液体洗剤などの日用品流体は 充填機を利用して容器に自動充填され製品として出荷されている しかし この充填プロセスにおいて 日用品流体の流体物性 ( 粘度または粘性係数 ) に起因する様々な技術課題が存在している

More information

数値シミュレーションによる次世代高速新幹線用低騒音 パンタグラフの開発 プロジェクト責任者栗田健東日本旅客鉄道株式会社 JR 東日本研究開発センター 著者 * 水島文夫 1 * 栗田健 1 * 山出吉伸 2 * 加藤千幸 2 * 上原均 3 * 3 廣川雄一 * 1 東日本旅客鉄道株式会社 JR 東

数値シミュレーションによる次世代高速新幹線用低騒音 パンタグラフの開発 プロジェクト責任者栗田健東日本旅客鉄道株式会社 JR 東日本研究開発センター 著者 * 水島文夫 1 * 栗田健 1 * 山出吉伸 2 * 加藤千幸 2 * 上原均 3 * 3 廣川雄一 * 1 東日本旅客鉄道株式会社 JR 東 数値シミュレーションによる次世代高速新幹線用低騒音 パンタグラフの開発 プロジェクト責任者栗田健東日本旅客鉄道株式会社 JR 東日本研究開発センター 著者 * 水島文夫 1 * 栗田健 1 * 山出吉伸 2 * 加藤千幸 2 * 上原均 3 * 3 廣川雄一 * 1 東日本旅客鉄道株式会社 JR 東日本研究開発センター * 2 国立大学法人東京大学生産技術研究所 * 3 独立行政法人海洋研究開発機構計算システム計画

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 不飽和土の力学を用いた 締固めメカニズムの解明 締固めとは 土に力を加え 間隙中の空気を追い出すことで土の密度を高めること 不飽和土 圧縮性の減少透水性の減少せん断 変形抵抗の増大 などに効果あり 締固め土は土構造物の材料として用いられている 研究背景 現場締固め管理 締固め必須基準 D 値 施工含水比 施工層厚 水平まきだし ( ρdf ) 盛土の乾燥密度 D値 = 室内締固め試験による最大乾燥密度

More information

第6章 実験モード解析

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

More information

. .. る試みがなされており, 橋本ら [6] は, 浮揚距離は数 μm~ 数百 μ W 崎町川 ~ 々 ~ 約匂デ対 rp たときで, 移動速度は O.7~ l. Omm/m 凶で あった. 実験において, 苅 ~ 毛焚書己情酸合会 z -...--~ :...",. 今 ~eo 細管を垂直に固定して, 印加周波数を 40kHz とし,

More information

Microsoft Word - 演習問題9,10

Microsoft Word - 演習問題9,10 流体力学 - 演習問題 9- 境界層の運動量方程式 境界層厚さ (1) 流速 の一様流れに平行に置かれた平板の層流境界層の流速分布を u= f (y/δ)= f (η) とおけば 境界層の厚さδ 壁面摩擦応力 平均摩擦抗力係数 C f はそれぞれ下式にて表されることを証明せよ 1 x x x 3 C 0 f 但し α= f 0 1 f d,β= d 0 df である () 平板に沿う一様流の境界層内の速度分布を

More information

数値計算:有限要素法

数値計算:有限要素法 ( ) 1 / 61 1 2 3 4 ( ) 2 / 61 ( ) 3 / 61 P(0) P(x) u(x) P(L) f P(0) P(x) P(L) ( ) 4 / 61 L P(x) E(x) A(x) x P(x) P(x) u(x) P(x) u(x) (0 x L) ( ) 5 / 61 u(x) 0 L x ( ) 6 / 61 P(0) P(L) f d dx ( EA du dx

More information

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

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

More information

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

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

More information

第 42 回流体力学講演会 航空宇宙数値シミュレーション技術シンポジウム 2010 論文集 161 低レイノルズ数における矩形翼とデルタ翼の空力特性比較野々村拓 1, 小嶋亮次 2, 福本浩章 2, 大山聖 1, 藤井孝蔵 1 1. 宇宙航空研究開発機構宇宙科学研究所,2. 東京大学大学院 Comp

第 42 回流体力学講演会 航空宇宙数値シミュレーション技術シンポジウム 2010 論文集 161 低レイノルズ数における矩形翼とデルタ翼の空力特性比較野々村拓 1, 小嶋亮次 2, 福本浩章 2, 大山聖 1, 藤井孝蔵 1 1. 宇宙航空研究開発機構宇宙科学研究所,2. 東京大学大学院 Comp 第 42 回流体力学講演会 航空宇宙数値シミュレーション技術シンポジウム 2 論文集 6 低レイノルズ数における矩形翼とデルタ翼の空力特性比較野々村拓, 小嶋亮次 2, 福本浩章 2, 大山聖, 藤井孝蔵. 宇宙航空研究開発機構宇宙科学研究所,2. 東京大学大学院 Comparative Study of Aerodynamic Characteristics of Rectangular and

More information

スライド 1

スライド 1 Tokyo Univesity of Agicultue and Technology パラフォイル型飛翔体の飛行安定に関する研究 15.1.11 (Fi) 金丸拓樹 ( 農工大 B4) 前川啓 ( 東海大 M) 森吉貴大 ( 農工大 M1) 山田和彦 (ISAS/JAXA) 東野伸一郎 ( 九大 ) 長崎秀史 ( 九大 ) 西田浩之 ( 農工大 ) 1/3 OUTIE 1. 背景, 目的. 安定評価の方法

More information

コンクリート工学年次論文集 Vol.27

コンクリート工学年次論文集 Vol.27 論文アクティブ赤外線法における照射光源の影響に関する基礎的研究 田中寿志 *1 仁平達也 * 鳥取誠一 *3 *4 栗田耕一 要旨 : アクティブ赤外線法に用いる照射設備は, はく離検知の程度に大きな影響を及ぼす そこで, 本研究では, 遠赤外線, キセノンランプ, およびハロゲンランプを用いた場合のコンクリート平板の照射試験および非定常熱伝導解析を行い, 熱伝導の挙動を確認した また, 照射条件を検討するために,

More information

PowerPoint Presentation

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

More information

密集市街地における換気・通風性能簡易評価ツールの開発 (その2 流体計算部分の開発)」

密集市街地における換気・通風性能簡易評価ツールの開発 (その2 流体計算部分の開発)」 OpenCAE ワークショップ 2013 2013.6.21 密集市街地における換気 通風性能簡易評価ツールの開発その 2 : 流体計算部分の開発 福本雅彦 ( 株式会社森村設計 ) 小縣信也 ( 株式会社森村設計 ) 勝又済 ( 国土交通省国土技術政策総合研究所 ) 西澤繁毅 ( 国土交通省国土技術政策総合研究所 ) 岩見達也 ( 国土交通省国土技術政策総合研究所 ) 概要 換気 通風性能簡易評価ツール

More information

していることが報告されている 1)2) しかし このような共鳴現象の解明のために 圧縮機内部空間の流れや音の伝播を実験的に調べるのは容易でない 一方 これまでの解析的な騒音現象解明では 羽根車の回転や高速の流れを無視しており 圧力変動の定量的な検討や予測ができるまでには至っていない また 解析上圧縮

していることが報告されている 1)2) しかし このような共鳴現象の解明のために 圧縮機内部空間の流れや音の伝播を実験的に調べるのは容易でない 一方 これまでの解析的な騒音現象解明では 羽根車の回転や高速の流れを無視しており 圧力変動の定量的な検討や予測ができるまでには至っていない また 解析上圧縮 遠心圧縮機の空力騒音低減の研究 プロジェクト責任者三浦治雄株式会社日立プラントテクノロジー 著者 * 三浦治雄 1 * 小林博美 1 * 三橋利玄 2 * 西村香純 2 * 大友洋 2 * 西川憲明 3 * 3 廣川雄一 * 1 株式会社日立プラントテクノロジー * 2 アドバンスソフト株式会社 * 3 独立行政法人海洋研究開発機構 利用施設 : 独立行政法人海洋研究開発機構地球シミュレータ利用期間

More information

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

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

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

層流撹拌槽解析 Laminar 循環流の概略図 流れ特性 撹拌レイノルズ数撹拌羽根のレイノルズ数旋回流の平均流速循環流量 動力 撹拌動力トルク動力数撹拌翼に与えられる軸方向の力 エネルギー散逸 平均エネルギー散逸率壁付近でのエネルギー散逸率高せん断ゾーンでのエネルギー散逸率全動力に対する翼付近でのエネルギー散逸の割合 翼付近のせん断 撹拌羽根付近のシアレート翼付近の有効粘性撹拌羽根付近の乱流せん断応力翼付近の高せん断領域の相対体積

More information

本日話す内容

本日話す内容 6CAE 材料モデルの VV 山梨大学工学部土木環境工学科吉田純司 本日話す内容 1. ゴム材料の免震構造への応用 積層ゴム支承とは ゴムと鋼板を積層状に剛結 ゴム層の体積変形を制限 水平方向 鉛直方向 柔 剛 加速度の低減 構造物の支持 土木における免震 2. 高減衰積層ゴム支承の 力学特性の概要 高減衰ゴムを用いた支承の復元力特性 荷重 [kn] 15 1 5-5 -1-15 -3-2 -1 1

More information

データ解析

データ解析 データ解析 ( 前期 ) 最小二乗法 向井厚志 005 年度テキスト 0 データ解析 - 最小二乗法 - 目次 第 回 Σ の計算 第 回ヒストグラム 第 3 回平均と標準偏差 6 第 回誤差の伝播 8 第 5 回正規分布 0 第 6 回最尤性原理 第 7 回正規分布の 分布の幅 第 8 回最小二乗法 6 第 9 回最小二乗法の練習 8 第 0 回最小二乗法の推定誤差 0 第 回推定誤差の計算 第

More information

<4D F736F F F696E74202D A957A A8EC0895E8D7182C982A882AF82E EF89FC915082CC82BD82DF82CC A83808DC5934B89BB A2E >

<4D F736F F F696E74202D A957A A8EC0895E8D7182C982A882AF82E EF89FC915082CC82BD82DF82CC A83808DC5934B89BB A2E > Techno Forum 2012 実運航における燃費改善のためのトリム最適化 株式会社 MTI 技術戦略グループ上級研究員堀正寿 1 目次 1. はじめに 2. 最適トリムの評価手法 2-1. オペレーションプロファイル調査 2-2. 水槽試験とトリム影響解析 2-3. 実船検証 3. トリムチャートと運用 4. まとめ 2 1-1 トリムの定義 1. はじめに 船尾喫水 (da) と船首喫水 (df)

More information

1

1 問題を解こう. 熱力学の基礎 問題. 容積 [m ] の密閉容器内に 温度 0[ ] 質量 0[kg] の酸素が含まれている この容器内の圧力を求めよ ただし 酸素の気体定数を R= 59.8[J/kg K] とする 解答 酸素の体積 V=m 質量 m=0kg なので 酸素の比容積 v=/0 m /kg である 式 (.) において ガス定数 R=59.8 温度 T=(0+7)K であるので 圧力

More information

9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L) du (L) = f (9.3) dx (9.) P

9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L) du (L) = f (9.3) dx (9.) P 9 (Finite Element Method; FEM) 9. 9. P(0) P(x) u(x) (a) P(L) f P(0) P(x) (b) 9. P(L) 9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L)

More information

問題 2-1 ボルト締結体の設計 (1-1) 摩擦係数の推定図 1-1 に示すボルト締結体にて, 六角穴付きボルト (M12) の締付けトルクとボルト軸力を測定した ボルトを含め材質はすべて SUS304 かそれをベースとしたオーステナイト系ステンレス鋼である 測定時, ナットと下締結体は固着させた

問題 2-1 ボルト締結体の設計 (1-1) 摩擦係数の推定図 1-1 に示すボルト締結体にて, 六角穴付きボルト (M12) の締付けトルクとボルト軸力を測定した ボルトを含め材質はすべて SUS304 かそれをベースとしたオーステナイト系ステンレス鋼である 測定時, ナットと下締結体は固着させた 問題 2-1 ボルト締結体の設計 (1-1) 摩擦係数の推定図 1-1 に示すボルト締結体にて, 六角穴付きボルト (M12) の締付けトルクとボルト軸力を測定した ボルトを含め材質はすべて SUS304 かそれをベースとしたオーステナイト系ステンレス鋼である 測定時, ナットと下締結体は固着させた 測定データを図 1-2 に示す データから, オーステナイト系ステンレス鋼どうしの摩擦係数を推定せよ

More information

輸送現象まとめファイル2017_01

輸送現象まとめファイル2017_01 第 5 章流れの解析 5.1 平行平板内流れ 流体の連続の式と運動量方程式を連立して解き 流体内部の速 度分布を求める 第一の例は 平行平板内流れの定常状態速度分 布である 非圧縮性流体 ( 密度 ρ が一定 ) 定常状態 ( t の項 が0) 重力以外の外力が働かないとき デカルト座標系の -! 二次元速度分布 v, ( ) を決める方程式は次の様に書ける 連続の式 : u + v = 0 (5-1)

More information

Microsoft PowerPoint - zairiki_3

Microsoft PowerPoint - zairiki_3 材料力学講義 (3) 応力と変形 Ⅲ ( 曲げモーメント, 垂直応力度, 曲率 ) 今回は, 曲げモーメントに関する, 断面力 - 応力度 - 変形 - 変位の関係について学びます 1 曲げモーメント 曲げモーメント M 静定力学で求めた曲げモーメントも, 仮想的に断面を切ることによって現れる内力です 軸方向力は断面に働く力 曲げモーメント M は断面力 曲げモーメントも, 一つのモーメントとして表しますが,

More information

Microsoft PowerPoint - 卒業論文 pptx

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

More information

Microsoft PowerPoint - 資料04 重回帰分析.ppt

Microsoft PowerPoint - 資料04 重回帰分析.ppt 04. 重回帰分析 京都大学 加納学 Division of Process Control & Process Sstems Engineering Department of Chemical Engineering, Koto Universit manabu@cheme.koto-u.ac.jp http://www-pse.cheme.koto-u.ac.jp/~kano/ Outline

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 反応工学 Reacio Egieerig 講義時間 場所 : 火曜 限 8- 木曜 限 S- 担当 : 山村 補講 /3 木 限 S- ジメチルエーテルの気相熱分解 CH 3 O CH 4 H CO 設計仕様 処理量 v =4.8 m 3 /h 原料は DME のみ 777K 反応率 =.95 まで熱分解 管型反応器の体積 V[m 3 ] を決定せよ ただし反応速度式反応速度定数 ラボ実験は自由に行ってよい

More information

Microsoft PowerPoint - qchem3-11

Microsoft PowerPoint - qchem3-11 8 年度冬学期 量子化学 Ⅲ 章量子化学の応用.6. 溶液反応 9 年 1 月 6 日 担当 : 常田貴夫准教授 溶液中の反応 溶液反応の特徴は 反応する分子の周囲に常に溶媒分子が存在していること 反応過程が遅い 反応自体の化学的効果が重要 遷移状態理論の熱力学表示が適用できる反応過程が速い 反応物が相互に接近したり 生成物が離れていく拡散過程が律速 溶媒効果は拡散現象 溶液中の反応では 分子は周囲の溶媒分子のケージ内で衝突を繰り返す可能性が高い

More information

NS NS Scalar turbulence 5 6 FEM NS Mesh (A )

NS NS Scalar turbulence 5 6 FEM NS Mesh (A ) 22 3 2 1 2 2 2 3 3 4 NS 4 4.1 NS............ 5 5 Scalar turbulence 5 6 FEM 5 6.1 NS.................................... 6 6.2 Mes A )................................... 6 6.3.....................................

More information

線形粘弾性 a.応力緩和とマクスウェル模型

線形粘弾性 a.応力緩和とマクスウェル模型 地殻 上部マントルのレオロジー 高温 長時間では, 岩石は流動する. 線形粘弾性体 ( 歪速度はせん断応力に比例 ) 流体力学 ( ナビエ - ストークスの式 ) べき乗型流動則 ( 歪速度はせん断応力のべき乗に比例 ) 比例係数の温度依存性が大きい. 参考書岩波講座地球科学 2 地球の物質科学 Ⅰ レオロジーと地球科学唐戸俊一郎グローバルテクトニクス杉村新 色々な深度の断層岩 ( 産業技術総合研究所,

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション デジタルエンジニアリング演習 流体 CAE 演習 (1) 2014 年 10 月 10 日 ( 金 ) 担当 : 杵淵郁也, 波田野明日可石川明克, 諸山稔員 流れの数値シミュレーション CFD = Computational Fluid Dynamics ( 数値流体力学 ) v t ( v ) v v 0 1 p 2 v F http://jda.jaxa.jp/ http://ansys.jp/applications/industry/auto/index_a.html

More information

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

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

More information

Microsoft PowerPoint SIGAL.ppt

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

More information

<4D F736F F F696E74202D C CC89C88A B8CDD8AB B83685D>

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

More information

Microsoft PowerPoint - 6.PID制御.pptx

Microsoft PowerPoint - 6.PID制御.pptx プロセス制御工学 6.PID 制御 京都大学 加納学 Division of Process Control & Process Systems Engineering Department of Chemical Engineering, Kyoto University manabu@cheme.kyoto-u.ac.jp http://www-pse.cheme.kyoto-u.ac.jp/~kano/

More information

物体の自由落下の跳ね返りの高さ 要約 物体の自由落下に対する物体の跳ね返りの高さを測定した 自由落下させる始点を高くするにつれ 跳ね返りの高さはただ単に始点の高さに比例するわけではなく 跳ね返る直前の速度に比例することがわかった

物体の自由落下の跳ね返りの高さ 要約 物体の自由落下に対する物体の跳ね返りの高さを測定した 自由落下させる始点を高くするにつれ 跳ね返りの高さはただ単に始点の高さに比例するわけではなく 跳ね返る直前の速度に比例することがわかった 物体の自由落下の跳ね返りの高さ 要約 物体の自由落下に対する物体の跳ね返りの高さを測定した 自由落下させる始点を高くするにつれ 跳ね返りの高さはただ単に始点の高さに比例するわけではなく 跳ね返る直前の速度に比例することがわかった (1) 目的球技において必ず発生する球の跳ね返りとはどのような規則性に基づいて発生しているのかを調べるために 4 種類の物体を用い様々な床の上で実験をして跳ね返りの規則性を測定した

More information

解析センターを知っていただく キャンペーン

解析センターを知っていただく キャンペーン 005..5 SAS 問題設定 目的 PKパラメータ (AUC,Cmax,Tmaxなど) の推定 PKパラメータの群間比較 PKパラメータのバラツキの評価! データの特徴 非反復測定値 個体につき 個の測定値しか得られない plasma concentration 非反復測定値のイメージ図 測定時点間で個体の対応がない 着目する状況 plasma concentration 経時反復測定値のイメージ図

More information

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

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

More information

7 π L int = gψ(x)ψ(x)φ(x) + (7.4) [ ] p ψ N = n (7.5) π (π +,π 0,π ) ψ (σ, σ, σ )ψ ( A) σ τ ( L int = gψψφ g N τ ) N π * ) (7.6) π π = (π, π, π ) π ±

7 π L int = gψ(x)ψ(x)φ(x) + (7.4) [ ] p ψ N = n (7.5) π (π +,π 0,π ) ψ (σ, σ, σ )ψ ( A) σ τ ( L int = gψψφ g N τ ) N π * ) (7.6) π π = (π, π, π ) π ± 7 7. ( ) SU() SU() 9 ( MeV) p 98.8 π + π 0 n 99.57 9.57 97.4 497.70 δm m 0.4%.% 0.% 0.8% π 9.57 4.96 Σ + Σ 0 Σ 89.6 9.46 K + K 0 49.67 (7.) p p = αp + βn, n n = γp + δn (7.a) [ ] p ψ ψ = Uψ, U = n [ α

More information