カオス水車のシミュレーションと その現象解析

Slides:



Advertisements
Similar presentations
Absolute Orientation. Absolute Orientation の問題 二つの座標系の間における剛体 (rigid body) 変換を復元す る問題である。 例えば: 2 台のステレオカメラから得られた3次元情報の間の関 係を推定する問題。 2 台のステレオカメラから得られた3次元情報の間の関.
Advertisements

Division of Process Control & Process Systems Engineering Department of Chemical Engineering, Kyoto University
22 ・ 3 積分形速度式 ◎ 速度式: 微分方程式 ⇒ 濃度を時間の関数として得るためには積分が必要 # 複雑な速度式 数値積分 (コンピューターシミュ レーション) # 単純な場合 解析的な解(積分形速度式) (a)1 次反応 1次の速度式 の積分形 [A] 0 は A の初濃度 (t = 0 の濃度.
Determining Optical Flow. はじめに オプティカルフローとは画像内の明る さのパターンの動きの見かけの速さの 分布 オプティカルフローは物体の動きの よって変化するため、オプティカルフ ローより速度に関する情報を得ること ができる.
『モデリング・シミュレーション入門』 井庭 崇 第5回 オートマトン(状態機械)
CGアニメーションの原理 基本技術 対象物体の動きや変形の設定方法 レンダリング技術
電磁気学C Electromagnetics C 7/27講義分 点電荷による電磁波の放射 山田 博仁.
・力のモーメント ・角運動量 ・力のモーメントと角運動量の関係
Akio Arimoto March 7,2011 Seminar at Tokyo City University
(Fri) Astrophysics Laboratory MATSUO Kei
プロセス制御工学 3.伝達関数と過渡応答 京都大学  加納 学.
エージェントモデル シミュレーション.
スペクトル法による数値計算の原理 -一次元線形・非線形移流問題の場合-
数楽(微分方程式を使おう!) ~第5章 ラプラス変換と総仕上げ~
常微分方程式と偏微分方程式 1.常微分方程式 独立変数が一個のもの 振動の運動方程式 2.偏微分方程式 独立変数が二個以上のもの
水中で落下する球体の運動.
Korteweg-de Vries 方程式のソリトン解 に関する考察
シミュレーション論Ⅰ 第4回 基礎的なシミュレーション手法.
大阪工業大学 情報科学部 情報システム学科 宇宙物理研究室 B 木村悠哉
次に 円筒座標系で、 速度ベクトルと加速度ベクトルを 求める.
工業力学 補足・復習スライド 第13回:偏心衝突,仕事 Industrial Mechanics.
大阪市立大学数学研究所 孝森洋介 共同研究者: 大川、諏訪(京大基研)、 高本(京大理)
1.Atwoodの器械による重力加速度測定 2.速度の2乗に比例する抵抗がある場合の終端速度 3.減衰振動、強制振動の電気回路モデル
流体のラグランジアンカオスとカオス混合 1.ラグランジアンカオス 定常流や時間周期流のような層流の下での流体の微小部分のカオス的運動
需要の価格弾力性 価格の変化率と需要の変化率の比.
数楽(微分方程式を使おう!) ~第4章 他分野への応用(上級編)~
動力学(Dynamics) 運動方程式のまとめ 2008.6.17
2.伝送線路の基礎 2.1 分布定数線路 2.1.1 伝送線路と分布定数線路 集中定数回路:fが低い場合に適用
シミュレーション演習 G. 総合演習 (Mathematica演習) システム創成情報工学科
Lorenz modelにおける 挙動とそのカオス性
準光速ロケットでのブラックホール旅行における時間の遅れ
大阪工業大学 情報科学部 情報科学科 学生番号 A03-017 犬束 高士
北大MMCセミナー 第38回 Date: 2015年2月13日(金)16:30~18:00 Speaker: 宮路 智行(明治大学)
二輪型倒立振子ロボットの製作 箱木研究室 T20R020 三留隼.
東京大学大学院情報理工学系研究科 Y.Sawa
SPH法を用いた重力の潮汐力効果のシミュレーション
奈良女子大集中講義 バイオインフォマティクス (9) 相互作用推定
独立成分分析 5 アルゴリズムの安定性と効率 2007/10/24   名雪 勲.
メンバー 梶川知宏 加藤直人 ロッケンバッハ怜 指導教員 藤田俊明
電磁気学C Electromagnetics C 7/17講義分 点電荷による電磁波の放射 山田 博仁.
システム制御基礎論 システム工学科2年後期.
深海底は海嶺軸から遠ざかるにつれ、規則正しく深くなる。 なぜか?
(昨年度のオープンコースウェア) 10/17 組み合わせと確率 10/24 確率変数と確率分布 10/31 代表的な確率分布
動力学(Dynamics) 力と運動方程式 2008.6.10
連続体とは 連続体(continuum) 密度*が連続関数として定義できる場合
全長 (Loa) 6.45m 水線長 (Lwl) 5.3m 最大幅 (Bmax) 2.48m 排水量 (Disp) 1045kgf
電磁気学Ⅱ Electromagnetics Ⅱ 8/11講義分 点電荷による電磁波の放射 山田 博仁.
ウィルスって どの位感染しているのかな? 菊池研究室  小堀智弘.
シミュレーションパラメータの設定 一次系の時間応答 二次系の時間応答
「データ学習アルゴリズム」 第3章 複雑な学習モデル 報告者 佐々木 稔 2003年6月25日 3.1 関数近似モデル
北大MMCセミナー 第7回 Date : 2013年3月6日(水) 16:30~18:00
シミュレーション論 Ⅱ 第1回.
データ解析 静岡大学工学部 安藤和敏
22・3 積分形速度式 ◎ 速度式: 微分方程式 ⇒ 濃度を時間の関数として得るためには積分が必要
22・3 積分形速度式 ◎ 速度式: 微分方程式 ⇒ 濃度を時間の関数として得るためには積分が必要
Massive Gravityの基礎と宇宙論
ガウス分布における ベーテ近似の理論解析 東京工業大学総合理工学研究科 知能システム科学専攻 渡辺研究室    西山 悠, 渡辺澄夫.
大阪工業大学 情報科学部 情報システム学科 学生番号 B02-014 伊藤 誠
【第六講義】非線形微分方程式.
目で見る一次変換 河合塾 数学科 生越茂樹 オゴセ シゲキ.
電気回路学Ⅱ 通信工学コース 5セメ 山田 博仁.
卒論中間発表 2001/12/21 赤道の波動力学の基礎 北海道大学理学部 地球科学科 4年 山田 由貴子.
シミュレーション演習 G. 総合演習 (Mathematica演習) システム創成情報工学科
宿題を提出してください. 配布物:ノート 3枚 (p.49~60), 中間アンケート, 解答用紙 3枚 (1枚は小テスト,2枚は宿題用)
How shall we do “Numerical Simulation”?
Massive Gravityの基礎と宇宙論
北大MMCセミナー 第23回 Date:2014年3月6日(木) 16:30~18:00 ※通常と曜日が異なります
信号データの変数代入と変数参照 フィードバック制御系の定常特性 フィードバック制御系の感度特性
(昨年度のオープンコースウェア) 10/17 組み合わせと確率 10/24 確率変数と確率分布 10/31 代表的な確率分布
FUT 原 道寛 学籍番号__ 氏名_______
Presentation transcript:

カオス水車のシミュレーションと その現象解析 カオス水車のシミュレーションと その現象解析 数理情報科学コース 中山研究室 北 いづみ

はじめに カオスとは、簡単な規則にしたがって起こるのに、不規則に揺れ動いて先の予測ができない動きのことである。 複雑な動きをするもの 水車モデル

水車モデルの解析を詳しく行い、カオス性を検証する。 参考文献 科学シミュレーション研究会:「パソコンで見る 複雑系・カオス・量子」 より導出 問題点 導出されている式に疑問がある カオス性の証明がない 研究目的 水車モデルの解析を詳しく行い、カオス性を検証する。

Java による現実水車のシミュレーション 回転速度、回転の減衰率、かごの穴の直径、注水量を初期値   として与える。 回転速度のグラフを表示     左回りを正、右回りを負の回転とし、時間推移による     速度の動きを表す。 回転の減衰率・・・水車の回転が遅くなる割合。 かごの穴の直径・・・かごの底には穴が開いていて                水はそこから抜け落ちている。

現実水車数学モデル 各かごの慣性モーメント 力のモーメント

現実水車モデルの問題点 かごに注水が行われるかどうかによって 水の変化量の微分方程式が切り替わって しまい、解析が難しい。  水の変化量の微分方程式が切り替わって  しまい、解析が難しい。 かごの中の水の総和を一定にできない。 理想的現実水車モデルを立てる

理想的現実水車モデルのしくみ かごの数はn個 注水は高さに比例してすべてのかご  に行われる。

Java による理想的現実水車のシミュレーション 回転速度、回転の減衰率、かごの穴の直径、注水量を初期値として与える。 回転速度のグラフを表示     左回りを正、右回りを負の回転とし、時間推移による     速度の動きを表す。 シミュレーションにおいてはかごの数は8個とする。 回転の減衰率・・・水車の回転が遅くなる割合。 かごの穴の直径・・・かごの底には穴が開いていて                水はそこから抜け落ちている。

理想的現実水車数学モデル 水量の総和

理想的現実水車の問題点 変数の数が多くなってしまうことから  解析が難しい。 かごが無限個あるとする理想水車を考える

理想水車数学モデル(W.Malkusモデル) 新しい関数を定義する

Java による理想水車のシミュレーション 初期値として、回転速度、注水量の高さに関する  定数を入力する。 回転速度のグラフを表示     左回りを正、右回りを負とし、時間推移による速さの動きを     表す。

ローレンツ方程式の導出 より導出

ローレンツ方程式の導出 線形変換 ローレンツ方程式

ローレンツアトラクタ

まとめ 現実水車モデルから導出した微分方程式系 よりローレンツ方程式が導かれ、そのローレンツ アトラクタがバタフライのように描かれた。  まとめ  現実水車モデルから導出した微分方程式系 よりローレンツ方程式が導かれ、そのローレンツ   アトラクタがバタフライのように描かれた。 水車はカオス的な振る舞いをすると考えられる。

終わり

角速度の変化量(角加速度)

理想水車数学モデル  水の変化量 水車全体の水の変化量 特異解より、水車全体の質量は一定とする。

角速度の変化量 水車全体の角運動量 角運動量を時間tで微分したものは 力のモーメントに等しい