電 話：076-445-6698
ＦＡＸ：076-445-6705ホーム ＞ 受講生の皆様へ ＞ 講義内容(シラバス)＞ IT･数値解析特論Ⅱ
講義内容（シラバス）
IT･数値解析特論Ⅱ
授業科目名 |
IT･数値解析特論Ⅱ |
||
担当教員(所属) |
|||
開講日程 |
10月－12月 木曜日 6、7限 |
単位数 |
２単位 |
 |
|
|
|
連絡先(研究室、電話番号、電子メール等) |
佐藤雅弘（電話076-445-6736, |
オフィスアワー(自由質問時間) |
|
|
|
授業のねらいとカリキュラム上の位置付け(一般学習目標) |
|||
近年，コンピュータの性能向上が著しい。それに伴って，製品の開発に数値シミュレーションを利用することが常識となっている。本授業では，流体力学，振動，波動，電磁気界解析から化学計算まで，その理論と使い方を基礎から解き明かす。さらに，解析可能な物理現象の拡大や，コンピュータのパワーアップの現状を探る。 |
|||
達成目標 |
|||
１）数値解析の基礎的知識を学ぶ。
|
|||
|
|
|
|
授業計画（授業の形式、スケジュール等） |
||
第１回 |
微分方程式と数値計算と流体力学 |
|
第２回 |
コンピュータ操作と第１回の例題演習 |
|
第３回 |
電界、磁界、熱などの定常場問題の解析法 |
|
第４回 |
第３回の例題演習 |
|
第５回 |
量子化学計算の概要 |
|
第６回 |
実習 |
|
第７回 |
実習 |
|
第８回 |
実習 |
|
第９回 |
振動・波動のシミュレーション，講義１ |
|
第10回 |
振動・波動のシミュレーション，講義２ |
|
第11回 |
振動・波動のシミュレーション，実習１ |
|
第12回 |
振動・波動のシミュレーション，実習２ |
|
第13回 |
最近のシミュレーション技術 |
|
第14回 |
ＨＰＣ(high performance computing)の現在１ |
|
第15回 |
ＨＰＣ(high performance computing)の現在２ |
|
|
|
|
キーワード |
コンピュータシミュレーション，流体力学，振動・波動，量子化学，半経験的分子軌道法，非経験的分子軌道法，化学結合，化学反応熱，分子構造，分子振動解析 |
履修上の注意 |
|
教科書・参考書等 |
|
成績評価の方法 |
|
関連科目 |
|
備考 |
|
|
|
IT･数値解析特論Ⅱ：授業計画
|
|
|
回 |
主題と位置付け |
学習方法と内容 |
１ |
微分方程式と数値計算と流体力学 (松島） |
微分方程式で記述される物理現象を計算機（数値計算）で再現するとはどういうことか？を概観する。特に, 流体力学における数値計算（一般にCFDと呼ばれる）について解説する。また，数値計算を実行するに当たり注意すべきことの例を，体験を通して確認する。授業は座学であるがEXCELの使えるノートパソコンの持参が望ましい。 |
２ |
流体力学方程式のオイラー座標に基づいた離散化 (松島） |
流体力学現象を記述するNavier-Stokes（N-S）方程式の意味を、保存則の観点から考察する。N-S方程式は偏微分方程式であるが，その連続の方程式を数値計算で扱えるように離散化していく手法の紹介を行う。（座学で行う） |
３ |
流体力学方程式のラグランジュ座標に基づいた離散化（I） (松島） |
其々の速度で流れている流体粒子に乗った座標系を用いて，流体方程式を解く手法の代表的なものであるCIP法を紹介する。CIP法の基礎と考え方のエッセンスを述べ, 簡単なCIP法流体解析プログラムを解説する。（座学で行う） |
４ |
流体力学方程式のラグランジュ座標に基づいた離散化（II） (松島） |
CIP法を用いた複雑な流れの解析について議論する。 |
５ |
量子化学計算の概要 (吉村） |
イントロダクション |
６ |
[実習] (吉村） |
パソコンを用いてChemDrawやChem3D等のソフトを用いて分子を描き，計算の入力データを作成し，WinMOPAC等の半経験的分子軌道計算プログラムで計算を行い，いろいろな物質のdipole moment, ionization potential, 結合距離，結合角や生成熱を計算し，実測値との比較を行う。また電荷やフロンティア軌道を図示し，反応性指数からいろいろな試薬がどこに反応し易いかを予測する。 |
７ |
[実習] (吉村） |
ＨＰＣでGAUSSIANプログラムを実行する準備と上記の入力データをGAUSSIAN用に変換し，実際に非経験的分子軌道法(ab initio法)によりいろいろな物質の構造最適化を行い，出力結果からdipole moment, ionization potential, 結合距離，結合角や生成熱を読み取り，実測値との比較を行う。また電荷やフロンティア軌道を図示し，反応性指数からいろいろな試薬がどの原子に反応し易いかを予測する。 |
８ |
[実習] (吉村） |
同様にいろいろな化合物の構造を最適化し，それぞれ振動解析によりIRスペクトルを予想し，またTDDFT法を用いて励起エネルギーを計算して、UVスペクトルをシミュレーションする。また化学反応をシミュレートするために反応式のそれぞれの物質を構造最適化し，エネルギー及び，振動解析から零点エネルギーを求め，また必要に応じて遷移状態の構造を最適化し，反応のエネルギー図を完成させる。 |
９ |
振動・波動のシミュレーション1 (佐藤） |
以降４回の講義で，音波や固体の波動（以降，音波に対して弾性波と呼ぶ）・振動の解析を行う。特に，時間とともにどういう動きをするかを解析し，簡単なアニメーションを作ることを目的としている。この講義は，一部，電磁波の解析にも利用できる。１回目は，解析の基本となる物理式について講義する。 |
10 |
振動・波動のシミュレーション２ |
解析手法にfinite difference time domain (FDTD) 法を使う。従ってFDTD法の基礎について講義する。前回講義した物理式をどのように近似して解析するかを具体的に説明する。コンピュータに実行させる際の言語は，FORTRANを使う。大規模な計算の場合は，FORTRANが一般的に使われる。この講義では，FORTRANの基本的な規則と簡単な命令について講義する。よく知られているC言語とそう大きな違いはなく，計算に特化されているのですぐに覚えることができると思う。 |
11 |
[実習] (佐藤） |
以降２回，実際にプログラミングを行う。また，その解析結果から，アニメーションを作る実習を行う。今回は，音波に焦点を当てて実習を行う。１・２次元の計算モデルを設定し，入力信号の入れ方，境界条件の設定法，吸収境界（音波が跳ね返ってこない境界）の設定法，音波の伝播の様子などを，具体的な実習を通して学ぶ。 |
12 |
[実習] (佐藤） |
最後は，固体に衝撃を与えたときの振動の様子を解析する。前回の音波の解析と同様に，計算モデルを設定し，プログラミングを行う。音波と弾性波は，その基本物理式が異なるので，境界条件の設定法や，吸収境界の設定法に違いがある。解析結果からアニメーションを作製し，弾性波の挙動を観察する。また，FDTD法による，固体のモード解析の方法について説明する予定である。 |
13 |
スーパーコンピュータ技術１ (西） |
最近，スーパーコンピュータ並みの高速計算の出来るコンピュータが安価に手に入るようになった。この後２回にわたって，高速計算機の現状について概観する。１回目は，並列計算，GPU( graphics processing unit) による高速計算について解説する。 |
14 |
スーパーコンピュータ技術２ (西） |
インターネットなどの広域のネットワーク上にあるコンピュータ資源を結びつけ，ひとつの複合したコンピュータシステムとしてサービスを提供する仕組みをグリッド・コンピューティングと呼ぶ。最近グリッド・コンピューティングのHPCへの応用が試みられている。この講義では，グリッド・コンピューティング技術について解説する。 |
15 |
最近のシミュレーション技術 (西） |
シミュレーション技術は，日々進化している。解析できる分野が大きく広がっている。この講義では，最近のシミュレーション技術を概観する。 |
|
|
|