ものづくりにHPCを活用するためのツールとケーススタディー

Vol.

2

Quarterly News Letter Winter 2014

Supercomputing Conference 2013 Best Poster Award

強スケーリング・プラズマ乱流シミュレーション
日本 原 子 力 研 究 開 発 機 構

シミュレーションの現場からソニー 先端マテリアル研究所
ソフトウェアライブラリ

PHASE

提供：日本造船技術センター

計算工学ナビ・レポート
自然科学研究機構
核融合科学研究所 教授

渡邉智彦氏

核融合炉の閉じ込め性能を解明する
プラズマ乱流シミュレーション
去る2013年11月にデンバーで開催されたSupercomputing Conference2013（SC13）
で、日本原子力研究開発機構のプロジェクトが高い評価を受け、ベストポスター賞を獲得しました。
このプロジェクトはHPCI戦略プログラム
分野4「次世代ものづくり」で開発された共通基盤技術の一部です。

日本原子力研究開発機構
システム計算科学センター
研究主幹

井戸村泰宏氏日本原子力研究開発機構
核融合研究開発部門
博士研究員

前山伸也氏

重水素と三重水素をプラズマ化し、原子核を融合させること径が乱流の最小スケールや駆動源のスペクトル 分布を決めるで発生する膨大なエネルギーを利用する核融合。2020年頃のため、イオンスケールと電子スケールが大きく乖離したマルチス実験炉稼働に向けて国際的な開発が進められている。

ケール乱流が存在する。

核融合のためには、原子核同士の 衝突 頻 度を上げるため、これまではイオンスケールの乱流と電子スケールの乱流は磁場中に約1億度の超高温プラズマを閉じ込める必要がある。

同時に計算することができず、別々に扱われていた。両方を同いかに効率よく熱を閉じ込めるかが重要だが、プラズマの温度時に解くためには数十万コア以上の並列化が必要だったからと密度が 高くなると乱流が発 生する。乱流は中心部の高温プだ。そのため両者の影 響を含んだマルチスケール乱流が 支配ラズマを外側へと拡散させてしまうため、高温状態を維持するする電子熱輸送の特性は良く分かっていなかった。には乱流の効率的な制御・抑制が重要な課 題である。これは日本原子力研究開発機構の前山伸也氏、同井戸村泰宏氏ら

現在フランスで建設中の国際熱核融合実験炉「ITER」の炉心は、通信コストを最小にして、約70万コアの「京」をフル活用性能やその先の原型 炉、そして実 用炉の経済 性を左右する課する手法を新規に開発し、電子スケールからイオンスケールま題でもある。核融 合プラズマのマルチスケール乱流 輸 送のシでを同時に解 析することを可能にした。並列 化 率 9 9. 9 9 9 9%ミュレーション解析が重要な理由はここにある。

（シックスナイン）以上を達 成したこの成果はSC13で高評 価

核 融 合プラズマの計 算には他の 乱 流とは異なる問 題 がある。まずひとつ目は、イオンと電子のスケールの違いだ。

を獲得する。解析コードの開発者である自然科学 研究機構 核融合科学研究所の渡邉智彦氏を交え、その詳細を聞いた。

燃料 粒子はプラズマ状 態なのでイオンと電子に電離している。これらの荷電粒子は磁力線の周りを旋回運動するが、イオ並列化率 99.9999%を達成

ンと電子の質量は3 桁も異なるため、その旋回半径はイオンが5mm程 度となるのに対し、電子は0.1mmとなる。この旋回半核融合プラズマの計算が通常の流体の乱流とは異なるもうITER 実験炉

イオン系乱流電子系乱流© ITER Organization, http://www.iter.org/

I T E R（イーター、国 際 熱 核 融 合 実 験 炉 ）は日本 、E U 、ロシア、米 国 、韓 国 、中 国 、インドの七 極によって進められている超 大 型 国 際プロジェクト。
2 0 2 0 年 頃の運 転 開 始を目指して現 在 、フランスのカダラッシュに建 設 中 。
「 I T E Rが成 功すれば 核 融 合は実 現できる」と関 係 者たちは語る。

1

9 9 . 9 9 9 9％の強スケーリング

ジャイロ運動論

京に対する最適化により極めて高いスケーラビリティーを実現

実効性能（TFLOPS）

1,000

6 次元位相空間

ジャイロ運 動 論

5 次元位相空間

京(6180億格子)

100
京(340億格子)

10
BX900(43億格子)

1
1,000

10,000

旋回半径

ベストポスター 賞

旋回中心軌道SC13でベストポスター賞を受賞した前山

100,000

1,000,000

コア数氏。
「 通信のオーバーヘッドを減らし、計算配分を工夫して極限まで計算の部分を高磁化プラズマにおける不安定性や乱流を記述する理論モデルがジャイロ運動論である。第一原理計算に基づく乱流シミュレーションを行うために用いられている。
磁力線を旋回する粒子運動に対して粒子分布関数を平均化できる利点がある。

めて性 能を引き出したことに価 値がある」と井 戸 村 氏はコメントした。コア数が増 加するエクサスケールコンピュータでは「エイトナイン」を達成する必要がある。

ひとつの点は、空間3次元に加えて速度2次元の、5次元の位「いろいろなタイプの通信があるわけです。隣接するノードとの相空間上の問題を解かなければならないところだ。そのうちの通信や、縦方向でガサッと集めるタイプ、ブロック（部分領域）

2次元を10 0×10 0のメッシュで表すとすると、3次元だけのモ内の通信とか。それらをトーラスのなかで配置するにあたってはデルに比べて1万倍の計算量が必要になる。
「京」クラスのコン様々なトレードオフがある。その最適化を70万コアに対してフルピュータが登場するまで、このような約1兆格子のマルチスケーに行ったのが彼のやり方の新しいところです」（井戸村氏）

ル乱流を計算することはできなかったのだ。
ではなぜ速度次元が必要なのか。通常の流体の場合は粒子エクサスケールと実験炉間の衝突が大きく、各粒子は熱力学的平 衡状態にあると仮 定されている。そのため速度分布を解く必要はない。だが核融合

これまでの研究で、磁場構造や駆動源のスケール、プラズマプラズマのような超高温状態では粒子間の衝突が小さく熱平の圧 力などによって異なる乱れのパターンが 現れることが分衡状態からずれる。だから粒子の速度分布関数の時間発展をかっている。また乱流からは乱れだけではなく「帯状流」と呼解かなければならなくなるのだ。

ばれる層流も非線 形効果によって生まれてくること、さらに乱本来は速度も3次元なのだが、プラズマを記述するために用流輸送は帯状流によって抑えられることも分かっている。

いられる理論モデル「ジャイロ運動論」では回転は平均化でき「流れの非線 形性ダイナミクスをうまく活かして、最適な閉じるため速度次元をひとつ減らせる。この特徴によりジャイロ運込め条件をシミュレーションで見つけていく。それが私たちの動論は核融合プラズマのほか、宇宙プラズマ乱流 解 析等にも大事な役割です」（渡邉氏）

使われている。

これまでの装置では、実際の核融合反 応が起こらない状態

SC13での受賞は、スケールが異なる2つの駆動源があるマでプラズマ閉じ込め実験を行ってきた。だがITERでは核融合ルチスケール乱流の問題を解いており、物 理学的にも興味 深反 応が起こる。核融合による高エネルギー粒子は主に電子へいという点と、京の3次元トーラスネットワークで5次元の問題エネルギーを与えながらプラズマを加熱するため、電子の熱輸を解くために通信最適化を行うという数値 計算 上の成果の2送が炉心性能を決める重要な要因となる。今回の前山氏らの点が高く評価されたからではないかと、実際に発表を行った前研 究は 、I T E R の 稼 働に先 駆けて電子 熱 輸 送 の 特 性を検 証山氏は語る。では、5次元の問題をどう3次元トーラスで効率化し、炉心性能を予測することに繋がるだろう。

したのだろうか。

ITERの運転開始は2020年頃と考えられている。次世代のエ「通信コストが 高いものが最小距離で通信を行うように、並クサスケールスーパーコンピュータの実用化と時期が重なる。

列 化に用いているM P Iのプロセス配 置 の 最 適 化を行 いまし「エクサスケールでITERの解析をすることになるだろうと期待た。さらに削減された通信をマスクするように、通信と演算のしています」（井戸村氏）

同時処理を行って、通信中にCPUがさぼらないようにしています。これで京をフル活用して、数十万コアでシミュレーションがITERの実 用化にシミュレーションが果たす役割は大きく、今後も井戸村氏らの研究は続けられていく。

できるようになりました」
（前山氏）

2

ものづくりの現場からシミュレーションがもたらす新知見が柔らかいエレクトロニクスを実現する
ソニーの先端マテリアル研究所は製品開発の基礎となる新材料の研究開発を行う部門です。柔らかいディスプレイの実現に不可欠な有機半導体の研究におけるシミュレーションの役割について取材しました。

ソニー
先端マテリアル研究所

2013年7月、ソニーからある論文が発表された。
『第一原理

有機物の多くは「インク」に加工して、印刷という単純な工程計算によるペンタセン、ルブレン、C 8 -BTBTのホッピング伝導で回路 化することが 可能だ。大面積の柔軟な基板にも対応しとバンド伝導の解析 』という題名から、有機材料に関するもや す い。また、印 刷 の原 料 使 用 効 率 は 、最 適 化 することでのだと分かる。その結論は、フレキシブル有機ELディスプレイ9 0%以上となる。無機半導体の場 合は数％とされるので、環の実現につながる重要な知見を示している。

境資源問題の観点でも有利な製造方法と言える。

*1

執筆者の小林一氏はソニー先端マテリアル研究所で有機材さらに、有機材料だけで構成されたディスプレイによって、料の解析を担当する研究者だ。

まったく新しい体験が生じると小林氏は語る。

「美しい動画を映すことができるフレキシブル有機 ELディス「オール有機ならば、印刷によってデバイスを製造できます。そプレイを実現するためには、動 作の速い有機半導体が必要でうすると、大きな面積のものでも非常に安くなるでしょう。さらす。別の言い方をすると、高いキャリア移動度を持つ有機材料に、今、私たちが重視しているのは、フレキシブルという性 質を作り出すことが製品開発の第1歩となります」

です。駆動回路も含めて全体を有機材料だけで作れれば、フィこれまで有機材料におけるキャリア（電子と正孔）の伝わりルム状の柔らかい素 材に印刷することができ、腕時計のよう方は「ホッピング伝導」が支配的であると考えられてきた。熱に手首に巻きつけたり、鉛 筆にくるくると巻き取るといった使励起されたキャリアが分 子から分 子へとランダムに飛び 移るい方が可能になります。これは、目に見える、新しい市場 価値現象であるホッピング伝導は遅い反 応である。つまり、この原です」

理を元に高速な半導体を作ることはできない。もうひとつの伝ソニーからは、ウェアラブルスマートフォン、インテリジェント

わり方「バンド伝導」は高速で、シリコンに代表される無機系ペーパー、ロールスクリーンディスプレイといった製品イメージの半導体はこの性質を利用している。

が 有 機トランジスタのアプリケーションとして提 案されてい遅い有機と速い無機を組み合わせてデバイスを作ることはる。巻き取り可能な有機TFT駆動・有機ELフレキシブルディス可能だ。従来の有機ELディスプレイはそうなっている。ではなプレイについては、すでに試作品が公開されている。
「柔らかぜ、ソニーは有機半導体だけからなるディスプレイの実現を目いエレクトロニクス」は将 来 のソニーを形 作るキーテクノロ指しているのだろうか？

ジーのひとつだ。

印刷で作る柔らかいディスプレイ

ファンデルワールス相互作用

無機半導体の製造には大 規模な設備が必要となる。一方、従来の「固い」ディスプレイに匹敵する表示能力を持つ「柔今後は電機メーカーでも物性シミュレーションがものづくりに活用されます

巻きとれる有機トランジスタ駆動フレキシブル有機ELディスプレイ

ソニー 先端マテリアル研究所
材料解析センター
シニアリサーチャー
小林一氏

2010年に発表された試作機。小さく巻き取れるほど柔らかでありながら、動画を滑らかに再生することができる。
http://www.sony.co.jp/SonyInfo/technology/technology/theme/or
ganictransistor_01.html

3*1 “Hopping and band mobilities of pentacene, rubrene, and 2,7-dioctyl[1]benzothieno[3,2-b][1]benzothiophene (C8-BTBT) from first principle calculations”
http://scitation.aip.org/content/aip/journal/jcp/139/1/10.1063/1.4812389

ホッピング伝導とバンド伝導

有機EL素子と有機TFT駆動素子を組み合わせたフレキシブルディスプレイの模式図

ホッピング伝導画素R

画素G

表1画素B
上面発光方式
有機EL部

有機半導体のバンド移動度

有機材料

ペンタセン*

ルブレン*

C8 -BTBT *実験値354031

計算/
vｄWなし

0.55

0.21

0.15

計算/
vｄWあり

585136

プラスチックフィルム

バンド伝導プラスチックフィルム

0.3mm

有機TFT部

分子から分子へキャリアが飛び移るホッピング伝導は遅い

プラスチックフィルムを基板とした有機材料ベースのフルカラー発光運動。
それに対して、波の性質を持つバンド伝導は速い運デバイス。将来的には、これを印刷工程によって製造することで、大動。高速な有機半導体の実現にはバンド移動度の高い画面のディスプレイを低環境負荷・低コストに実現できると考えられて材料の選択が重要。

いる。

vdW：ファンデルワールス相互作用

単位： cm 2/Vs

らかい」ディスプレイを実現するためには、先述のキャリア移しました。つまり、ファンデルワールス力のおかげで分子は思っ

動 度 の 問 題 を 片 付 ける 必 要 が あ る 。具 体 的 に は 移 動 度たほど振動していないこと、その結果、バンド移動度は3 6 - 58

10cm 2 / Vs以上という性能を持つ、速い有機材料の開発が必

c m 2 / Vsという実 験 値に近い高い 値となることがわかりまし

要だ。

た」表 1

「有機半導体ではホッピング伝導が主と言われていました。その移動度は5cm 2 / Vs以下で、半導体にしたときの動作が遅く高精細の動画は映せません。一方のバンド伝導なら、もっと高有機半導体の開発 者を勇気づける知見がシミュレーションによって得られた瞬間である。
今後、電機メーカーの研究開発の現場でも、シミュレーショい値が望めるのですが、理論限界がよく分かっていなかったのンが活用される場面は増えるはずと小林氏は語る。

です。有 機 材 料でバンド伝 導は起こらないという人も多かっ「第一原理計算の利用はまず化学メーカーで立ち上がり、次た。そんななか、実験によって35や40といった移動度を示す例

いで製薬の分野で普及しました。電機メーカーはどちらかといが出てきました。でも、それを理論に基づいた計算で説明するうと遅かったので す。その 理 由は、デバイスの開 発には化 学ことはできませんでした」

メーカーや製薬会社が求めるものよりも大きい系が必要だか

従 来の手法で計算すると、バンド伝導が 過小評 価され、材らでしょう。計算機の性能が上がり、PHASEのようなソフトウ料の実力を測ることができない状 況だった。その原因は分 子エアが登場したことで、これから電機業界からのアウトプット間の相互作用にある。

が増えてくると思います」「なかなかバンド伝 導の計算ができなかった原因は、電子と今後、電機メーカーにおいてシミュレーションはどんな役割フォノンの相互作用、それからファンデルワールス力という分を果たすのだろうか？

子 間 相 互作用の取り扱いが 難しいことにありました。とくに「計算によってムダな試作を省くことができれば、開 発 のスファンデルワールス力は遠くまで影響が及び、裾野のほうまでピードがあがるはずです。デバイスの試作にはたいへんなコスキレイに計算しないと正しい結果が出ないのです。複 雑な計

トがかかりますが、材料をスクリーニングする際、明らかにダメ算で すが、これを省 略すると、分 子がフワフワ動きまわってなものをはねるだけでも、ずいぶん喜ばれるでしょう。また、材キャリアがすぐに散 乱し、その結果移動度が 下がってしまう。

料だけでなく、デバイスや、それを組み立てたセットのレベルま従来の計算手法でバンド伝導の値が過小評 価された原因は、で通してシミュレーションできる、一種のマルチスケールシミュここにありました」

レーションが 実現すると、製品開発の様 相は変わると思います。10 年後くらいには、計算機の 性能もそれが 可能なレベルPHASE から得られた新たな知見
この状 況を解決したのが、物質・材料研究機構の大野隆央氏を中心に開発が 進められている第一原理バンド計算ソフト
ウェア PHASE である。PHASEがファンデルワールス相互作用に対応したことで、実 験値に近いバンド移動度が計算で求められるようになった。
「ファンデルワールス力を実 際に計算できるようになったのに達するのではないでしょうか」★用語解説

ペンタセン、ルブレン、

C 8 - BTBT

有機半導体の材料となる化学物質の名称。このなかで最も新しいC 8 - BTBTは印刷プロセスとの親和性が高く、インクジェットで塗布し半導体薄膜を形成する例が産業技術総合研究所からは、計算機の進歩と大野さんたちの開発の成果でしょう。新し発表されている。

いPHASEを使って、3 種類の有機材料についてバンド移動度

http://www.aist.go.jp/aist_j/press_release/pr2011/pr20110714を計算すると、どれも従来の結果より100倍程度大きな値を示/pr20110714.html

4シミュレーションソフトウェア・ライブラリ

第 一 原 理 電 子 状 態 計 算ソフトウェア

PHASE のシステム構成
統合環境 GUI
擬ポテンシャル作成

第一原理電子状態計算

CIAO

PHASE

擬ポテンシャルデータベース

PHASE - Viewer

量子伝導解析

ASCOT

第一原理計算とは、分子の構造を経験的データを用いずにすでにPHASEが知見を提供した事例を2つ挙げる。

基本的な方程式だけを使って求めることを意味する。PHASE

ひとつは、立体的な構造によってリーク電流を低減するトラは密度汎関数理論という計算手法を基に、金 属・絶 縁 体・半ンジスタ「FinFE T」、もうひとつは原子1個 分の厚さしかない導体といったあらゆる材料に対して第一原理電子状態計算を

炭 素シート＝グラフェンからなる低 消費 電 力 デバイス「グラ適 用し、高精度かつ大 規模にその性質を分析することができフェン・トランジスタ」の例である。

るプログラムだ。
同じ分野の他のプログラムと比較してPHASEがとくに優れ

FinFETの事例では、走査型トンネル顕微鏡を用いても正確には見ることができないシリコンの表面構造をPHASEによっている点は、その強力なスケーラビリティーである。て可視化した。表面構造モデルを探索的に計算して（1,0 0 0原

「世界 的に見ると同 様の 機 能を持つソフトはいくつかありま子の系に対して約1,0 0 0のモデルを適用）、最も安定なモデルす。しかし、パソコンから京クラスのスパコンにまで対応していを見つけることで、実際の原子配置を解明した。

る無償ソフトはPHASEだけでしょう。スケーラビリティーに関グラフェン・トランジスタの事例では、静的解 析では妥当性しては十分時間をかけてやってきたので、かなり先を行っていに疑問が生じる薄膜成長過程の解析を第一原理計算によってると自負しています」と語るのは物質・材料研究機構（NIMS）行い、大面積で厚さが均一な高品質グラフェンの製造につながの大野隆央氏。

る新たな知見を提供した。

数万原子規模の系に対応可能な高並列性能と強スケーリン
グ性能を持つPHASEは、革新的な特性を持つ次世代ナノデバ

こうしたシミュレーションによって得られた知見がデバイス開発を助けると大野氏は語る。

イスの実現につながる研究開発に活用されている。

次世代半導体ナノデバイス
2025年にIT機器の消費電力は全消費電力の約25%に達すると言われ、半導体の消費電力低減が社会的に求められている。また、プロセスの微細化が限界に近づくなか、新しい半導体素子の構造と材料の開発は急務だ。こうした課題に対して、5物質・材料研究機構（NIMS）
大野隆央氏

シリコン表面構造の大規模構造探索

高精度電子状態解析

5員環（白線）の陰に隠れるよ

ファンデルワールス相互作用によるたわみ

うに、7員環（黄色の円内）が存 在 することが P H A S Eを使った解析によって判明した。

トンネル走 査 顕 微 鏡で観測したシリコン表面には、5員（5角形の化合物）が並んでいる。
シリコンカーバイド・グラフィン系（446原子）における、ファンデルワールス相互作用の影響を計
フィン電界効果トランジスタ素子（FinFET）で重要な110面と呼ばれる表面の構造を、探索的に算した例。従来の計算手法（Dion法）に対して2100倍の高速化を実現している。これにより、解析した例。走査型トンネル顕微鏡による観測では分からなかった詳細な構造が明らかになった。

有機デバイスへの適用が可能となった。

「デバイスのなかで 何が 起こっているかを把 握したいときに理でできるとは思われていなかったファンデルワールス力のよ

PHASEは有効です。デバイスの生成過程だけでなく、劣化す

うな相 互作用も、最 新のPH AS Eでは記 述することができまる原因を計算で確かめるためにも使えます。そうした解析の結す。それによって2桁から3桁増えた計算量には、新しいアルゴ果をもとに技術者はアイデアを生み出し、プロセスや機能を良リズムを導入することで対応しました」

くすることができるでしょう」

PHASEで 解決可能な問題は着実に増えているという大野氏。続いて今後の開発の方向性について聞いた。

次世代エネルギーデバイス
ピコ秒からナノ秒へ
次世代半導体の研究開発と並んでPHASEの使用事例が増えているのが、次世代エネルギーデバイスの分野である。

「今後はプロセスのシミュレーションが求められると考えてい「第一原理計算が半導体産業とともに発展したこともあって、ます。いま我々が京を使ってやっているのは、1,0 0 0原子程度

従来はシリコンを対象とする計算が主流でした。しかし、近年で数ピコ秒のシミュレーションです。これでもイニシャルなプ

は電池や触媒 材料といったエネルギー関連の取り組みが増えロセスは分かりますが、トータルで 何ができてくるかまでは分ています。また、半導体だとしても、もっとソフトなもの、たとからない。時間をもっと長くして、ナノ秒まで行きたいんです。

えば有機材料やバイオ的なものに適用する例が多くなってきま

エクサスケールコンピュータでは1万原子に対してナノ秒単位した。日本のお家芸である材料の分野にPHASEを導入する効の計算をすることが目標です。ただし、そのためにはさらに並果は大きいと見ています」

列 性 能をあげ て、1ステップあたりの 処 理 時 間を1/10 0 からエネルギー分野の研究は、新材料の研究と表裏一体で進め

1/10 0 0に短 縮する必要があります。それができれば、次世代られていることが多い。特性がまだよく分かっていない材料を半導体デバイスや高効率二次電池のようなエネルギーデバイ扱うとき、シミュレーションは欠かせない手段だ。様々な材料

スの開発において、材料スクリーニング、デバイス界面設 計、に対応するため、PHASEの開発は、スケーラビリティーだけプロセス制御 設 計にPHAS Eを活用することが 可能となるででなく機能面の拡張も継続的に進められている。

しょう」「シリコン系というのは非常にきれいな系で、欠陥といっても1

現在、PHASEはCISSライセンス*で公開されており、誰で万個から数万個に1個をコントロールする話です。一方、有機やもダウンロードして無償で試用することができる。

エネルギー材料というのはシリコンよりもずっとごちゃごちゃ「Windows版とLinux版を提供していますので、興味を持ったしているんですね。計算規模が大きいだけでなく、電子の相互方はまず使い慣れている手元のPCで試してください。ダウン作用が複 雑です。ものづくりに役立てるためには、そうした複ロードしたファイルには、比較的小規模なサンプルデータがい雑さをハンドリングしなくてはいけない。数年前までは第一原くつも入っていて、すぐに動かしてみることができます」
今後の開発方針としては、並列性能の強化だけでなく、ユーザビリティーの向上も重視しているという。

★用語解説

CISS ライセンス
非営利目的の使用は自由です。営利目的の利用については商用「これまでもワークショップやユーザーの集まりを通じて、使い勝手や計算の安定性に関するフィードバックを集め、プログラムの改善につなげてきました。PHASEは誰でも安心して使ライセンスの申請が必要です（許諾に3ヶ月程度の時間がかかりまえるようにしたいと考えています。ユーザビリティーを良くするす）。問い合わせ先：東京大学生産技術研究所 革新的シミュレー

ことは口で言うほど簡単ではありませんが、今後も力を入れてション研究センター s o f t w a r e @ c i s s . i i s . u - t o k y o . a c . j p

やっていきます」

6

今号の表紙京を用いた船体まわりの流れ解析
この画は船が進むときに生じる非常に小さな渦を忠実に再現したものです。この計算では計算領域全体で320億点の計算格子が必要でした。従来のシミュレーション技術を船舶設計に用いるには精 度 、信 頼 性が低すぎましたが、F r o n t F l o w / b l u eと京コンピュータを使う事により、飛躍的に精度を向上させることができ、巨大な水槽に大型模型船を浮かべて行う曳航水槽試験と同等の精度が確保できるため、設計コストの削減が期待されています。
（一般財団法人 日本造船技術センター 西川達雄）

編集後記
SC13でベストポスターアワードの受賞式に参加した前山さんは28歳の新鋭研究者。そのときの写真を拝見すると、しっかりガッツポーズ。それを見てこちらまで嬉しくなりました。日本のHPCソフトウェア技術が世界的な発表の場で評価されたと考えると、重要な１シーンです。
（F）

計算工学ナビ オフィシャルサイト
本 誌のダウンロードや、ソフトウェアライブラリ、その他の読み物を提 供しているオフィシャルW e bサイト

http://www.cenav.org/

計算工学ナビ Vol.2
発行日： 2014年1月31日
発 行： 東京大学生産技術研究所
革新的シミュレーション研究センター〒153-8505 東京都目黒区駒場4-6-1
oﬃce@ciss.iis.u-tokyo.ac.jp