ものづくりに HPC を活用するためのツールとケーススタディーQ U A RTE R LY N E WSLE TTE R / SU MME R 2015

提供：日本原子力研究開発機構

電極の超大規模解析による燃料電池設計プロセスの高度化

VO L .

東京大学 鹿園直毅

熱流体解析システム

FrontFlow/violet-HC

オープンソース遠隔可視化ソフトウェア

PBVR

7

計算工学ナビ・レポート

電極の超大規模解析による燃料電池設計プロセスの高度化
ポスト京に向けた研究開発が平成26年度にスタートしました。重点課題６『革新的クリーンエネルギーシステムの実用化』（本誌 Vol.6参照）のサブ課題のひとつとして設定されたのが、『気液二相流および電極の超大規模解析による燃料電池設計プロセスの高度化』というテーマで、燃料電池システムの実用化に大きく影響する技術開発が進められています。このプロジェクトで固体酸化物形燃料電池（SOFC）に関連する領域を担当する東京大学生産技術研究所の鹿園直毅教授に話を聞きました。

はじめに鹿園さんのバックグラウンドを教どんなものでしょうか？

えてください。

東京大学生産技術研究所エネルギー工学連携研究センター

鹿園直毅 教授

SOFC の開発に大きく貢献できると考えています。

SOFC は700度から800度という高温
私のもともとの専門は「熱」です。メーで動作します。燃料極はセラミックスとこれまではそのようなシミュレーションをカーで蒸気サイクルや冷凍サイクルのコニッケルの混合物なんですが、ニッケルは実現することは難しかったのですか？

ンポーネントを作る仕事をしていました。

融点が高くないので長い時間のなかでどう

2002年に大学へ移ってから、燃料電池のしても動いてしまい、最悪の場合、凝集し、従来、我々が使うことができた計算機でなかでも熱の問題が重要な SOFC の研究、粒成長してネットワークが切れてしまいまは10μ m 強の領域に対し100nm 程度のなかでも一番困っている電極の研究にシフす。そうすると、電子が流れなくなって、メッシュが限界でした。これでは解像度不トしてきたという経緯があります。

電極が機能しなくなります。

足で精度が出ず、満足できる水準まで到達

こ の プ ロ ジ ェ ク ト で は SOFC と発電設備として使われる SOFC の運転できなかったのです。ポスト京を使うことPEFC( 固 体 高 分 子 燃 料 電 池 ) の 両 方 を期間は最低でも10年間。長期間の安定性で、20μｍの領域に対して10nm の解像

扱 い ま す が、 私 は 主 に SOFC を 担 当 し、という視点が大事です。しかし、初期の構度を実現したいと考えています。

PEFC はみずほ情報総研と京都大学を中造を変えるたびに10年かけて実験するわ心とするチームが担当します。SOFC とけにはいかない。これが加速実験をしたい

PEFC の両方に共通の課題は電極で、イオという企業のニーズが強い理由です。

ン、電子、物質の輸送、そして電気化学を現在のモデルでも温度を上げれば加速は連成させて電極の反応全体をシミュレーするのですが、温度依存性が異なる現象がションすることが目標のひとつです。

あるため、通常の温度による実験と一致しそれができて初めて、高い精度でモデルに関する議論を行えるようになります。
もう少し電極について教えてください。どんな電極が良い電極なのでしょう？

ません。真の加速実験を行うためには、要燃料極はニッケルとイットリア安定化ジ

この記事では主に SOFC の側から見た課因ごとに正しくモデル化したシミュレー

ルコニア（YSZ）からなる多孔体で、こ題について伺いたいと思います (PEFC 関ション技術を作り、それを検証する必要

こをイオンはくねくね曲がりながら流れて連の開発概要については次ページのカコミがあります。10年を短期間で予測するシいきます。まっすぐ流れている状態を基準

参照 )。SOFC の電極に特有の難しさとはミュレーションができれば、より高性能なにすると、7倍から8倍くらい遠回りして［図1］固体酸化物形燃料電池

燃料 H2,CO

排気 H2O,CO22e- 電子

燃料極
電解質

O2- 酸素イオン

空気極
O2 空気中の酸素

2

計算工学ナビ Vol.7(SOFC)の発電原理

電極は燃料極、電解質、空気極の3層で構成され、それぞれに組成の異なるセラミック系材料が使われる。

負荷SOFC: Solid Oxide Fuel Cell
PEFC: Polymer Electrolyte Fuel Cell
YSZ: Yttria-stabilized Zirconia
FIB: Focused Ion Beam
SEM: Scanning Electron Microscope
TPB: Triple Phase Boundary

PEFC 関連のアプリ開発の概要
大規模気液二相流解析コード

大規模 PEFC 電極解析コード

・OpenFOAM、FrontFlow/blue 等のオープンソー・実寸電極サイズへの拡大（界面格子サイ

路界面の水滴挙動、マニホールドの二相流動を再現

・電極形成プロセスシミュレーション、発電性能予測

スをベースに表面張力・壁面濡れ性の影響を含む流・数百億グリッド( 数百枚セル積層スタック、セル1枚あたり1億グリッド規模 )の気液二相流解析の達成

ズ20nm 以下、電極厚さ20μ m）

PEFC 全体シミュレータ（工学モデルの高度化）

高精度二相流動モデル、高精度発電モデル、流動モデ
ル、伝熱モデル

コスト低減に向けたスタック設計指針いる。これを屈曲度ファクターと言います。

この技術を産業界が利用できるのはいつ頃測にも効いてきます。

でしょうか？

遠回りするほど損失は大きくなりますから、本当はまっすぐスーッと流れてほしい。つ

長期間運転を前提とする電極の劣化予測、まり、屈曲度ファクターが低いほどいい電製造時の焼結予測、そして構造の最適化と設計の現場でポスト京クラスの計算機が極と言えます。

いう3つの課題を解決するコードの開発が使えるようになるのは10年後でしょうか。

始まったところですね。

それまでには間に合わせたい。

このプロジェクトの課題のひとつは、電極構造を最適化する技術の開発で、適切なSOFC の技術競争は10年どころか、20評価関数を与えてシミュレーションを行う

まだ名前を付けていないのですが、我々年後、30年後まで続くと思います。さらことによって、たとえば屈曲度ファクターが「大規模 SOFC 電極解析コード」と呼にその先の炭化水素燃料が枯渇したあとはの低い構造を見つけることができると考えんでいるものをベースに、今後はポスト京水素社会となって PEFC が活躍するでしょています。に対応するための並列化を進めていきます。

う。その頃になると、SOFC は廃棄物やただし、シミュレーションによっていい世界的にみて、SOFC の技術開発に大バイオマス等の炭化水素燃料の発電と可逆形がわかっても、それが非現実的な構造で規模計算を適用する我々のようなアプロー発熱器としての利用が中心になるでしょう。

はあまり意味がありません。実際に作れなチをとっている人はほとんどいません。

世界的な化石燃料の利用が減少するまいといけない。実はここが SOFC の難し我々もまだ分からないことだらけです。電で、少なくともこの50年くらいの間、日いところなのです。

極の内部にはいろんな不純物があって、高本のメーカーは競争力を維持してほしいと

SOFC の電極は瀬戸物のお茶碗と同じ温でありとあらゆることが起こります。反

願っています。そのためにはポスト京のよように、セラミックスの粉をこねて焼いて応によって元素が動き、その動きによってうな計算機を使ったシミュレーションがす作ります。そうしないと安くできない。そ別の反応が起こる。その組み合わせはものごく効いてくるでしょう。もし我々がやらのかわり結果は「出来なり」です。良い構すごい数ですが、ひとつひとつ潰していっなくても、誰かがやるはずです。早くやら造の電極を作るためには、焼結のプロセスて本質に迫りたい。これが我々のアプローなければ意味がない。そういう競争のなかを予測する技術が必要で、これは劣化の予チです。にいると思っています。

［図2］実電極構造を用いた3次元分極シミュレーション

FIB-SEM は集束イオンビーム (FIB) により試料を高精度に研磨し，走査電子顕微鏡 (SEM)により断面画像を連続的に撮像する技法。これを用いて電極の3次元構造を明らかにし、その構造データに酸化物イオン ･ 電子 ･ ガス種の拡散と、三相界面(TPB)での電気化学反応を連成した3次元格子ボルツマン法数値シミュレーションを適用する技術を開発した。ポスト京に向けて、より大規模なシミュレーションを可能にし、モデルの精緻化につなげる。

［図3］耐久性能９万時間を前提とした SOFC 電極構造の劣化予測

左図の灰色の塊はニッケルが連結している初期の状態を示している。時間が経過するにつれ、ニッケルが凝集・粒成長して連結が切れた状態（右図の赤い部分）が発生する。この変化過程を高速かつより大規模にシミュレーションすることで、10年間とNi-YSZ 燃料極内酸化物イオン電気化学ポテンシャル分布

Ni-YSZ 燃料極内イオン電流および電子電流分布（赤＝イオン電流、青＝電子電流）

いう長い運転期間における電極劣化を短期間で予測可能にするのが、本プロジェクトの課題のひとつとなっている。

3

計算工学ナビ・レポート

オープンソース遠隔可視化ソフトウェア

PBVR

数値シミュレーションの規模が大きくなるにつれ、可視化システムに対する要求も厳しくなっています。とくに京コンピュータのような大規模並列環境に適したクライアント・サーバ型の可視化システムの開発は重要なテーマと言えるでしょう。この課題に取り組んでいる日本原子力研究開発機構の河村拓馬さんに、可視化ソフトウエア "PBVR"について解説をお願いしました。

日本原子力研究開発機構
システム計算科学センター 高度計算機技術開発室 ( 任期付 )研究員
博士
（工学）

河村拓馬

従来の遠隔可視化システムは、可視化処

マ乱流シミュレーションの結果データに対そこで、豊富な計算資源を持つサーバ上理速度やメモリ使用量、データ転送量がボして、右上の伝達関数エディタにより磁場で可視化処理を行い、クライアントに可視トルネックとなって、可視化作業に膨大なポテンシャル・静電ポテンシャル・密度揺化要素を転送するクライアント・サーバ型時間が必要とされていました。

動の三変量に関する描画を制御し、左上の可視化が重要視されるようになりました。

遠 隔 可 視 化 ソ フ ト ウ ェ ア PBVR

スクリーンに可視化結果を表示しています。

しかし、既存の商用ソフトによるクライ(Particle-Based Volume Rendering) はク

PBVR は以下の Web ページからダウン

アント・サーバ型可視化では、ボリュームライアント・サーバシステムを採用し、大ロードが可能になっています。

データの大規模化に比例して可視化データ

規模な結果データをサイズの小さい可視化

http://ccse.jaea.go.jp/ja/download/software.htmlとして使用されるポリゴンの数が爆発し、用粒子データに変換して転送することで、これらの問題を解決し、ボリュームレンダリングによるインタラクティブな可視化を実現しています。

データ転送の増加、および、描画処理速度従来のクライアント・サーバ型可視化システムの問題点の低下という問題が発生します。特に、高度な奥行き表現手法であるアルファブレン
ディングを用いる場合には、可視化データ

PBVR はスーパーコンピュータ「京」を京に代表されるペタスケールスーパーコの探索や並べ替え処理に起因するノード間はじめとする様々な超並列計算機に対応し、ンピュータの登場によってシミュレーショコミュニケーションが多発します。そのたクロスプラットフォーム環境で利用できまン規模が増大し、得られる結果データ（ボめ描画処理速度のストロングスケーリングす。また高度な伝達関数設計機能を実装しリュームデータ）も飛躍的に大規模化してが達成できず、大規模並列環境を利用してており、ボリュームレンダリングだけでないます。そのため、昔ながらの「ユーザもインタラクティブな可視化は困難というく、等値面、断面、そして多変量データ解PC にボリュームデータを転送して、商用問題がありました。

析にも対応しています。

ソフトを利用してポスト処理」という方法

これらの課題を解決するため、京都大

図1に PBVR の使用例を示します。

は、データ転送量の増大やユーザ PC のメ

学 小 山 田 研 究 室 で 開 発 さ れ た Particle-GT5D(Gyrokinetic Toroidal 5

モリ・ストレージ容量の不足という問題にBased Rendering（PBR）の技術を利用より、非現実的なものとなりました。

し、大規模並列分散環境に適したクライア

Dimensional Eulerian code) によるプラズ［図1］遠隔可視化システム PBVR の GUI

PBVR は GNU Lesser General Public License(LPGL)で公開されている。クライアントプログラ
ムは Linux、Mac、Windows の各 OS で動作する。

［ 図 3］GT5D デ ー

タの静電ポテンシャ
ルの可視化結果

4

計算工学ナビ Vol.7

［図2］PBVR の処理の概要

クライアント（PC）
粒子データ

ボリュームデータ

粒子データ

サーバ
（スーパーコンピュータ）

要素(セル)

画像面

グラフィックカードによる粒子投影と画像生成処理

並列環境におけるデータの生成

グラフィック
カード(例) 4.2GBのボリュームデータ(例) 120MBの粒子データ

CPUメモリストレージ

粒子データの転送

ント・サーバ型可視化システムを開発しま送され、インタラクティブな可視化が可能

X5570・48コア)、クライアント (Xeonました。それが遠隔可視化システム PBVRとなります。また、PBVR は前処理としてE5 + Quadro K5000)を用いました。

なのです。

フィルタリングによるボリュームデータの画像が表示されるまでの時間（合計）に分割を行い、並列読み込みとマスター・スついて、PBVR は EnSight の約30倍高速レーブ型の動的負荷分散により効率的な並で、60frame/sec という、いわゆるイン列処理を実現します。

タラクティブフレームレートを達成してい効率的な並列処理を実現する
遠隔可視化システム PBVRます。またクライアント使用メモリについPBVR の重要な特徴として、ポリゴンベースの可視化手法と異なり、可視化用の粒子データのサイズが画像解像度で決定さ京コンピュータに対応しストロングスケーリングを達成れるため、元の大規模データのサイズと比粒子データの処理に京 ( 約1000コア )べて十分に小さくなるということが挙げらを利用し、GT5D から得られた大規模デーれます。

タ ( 約3億格子、100タイムステップ ) の

図2に PBVR の 処 理 概 要 を 示 し ま す。

可視化を行いました [ 図3]。ても PBVR が小さく、効率的な可視化システムだということがわかります。

より高画質な可視化へ
遠隔可視化システム PBVR は、可視化インターフェースや大規模可視化手法の両サーバは大規模データを読み込み、伝達関その結果、従来の可視化システムでは観点から現在も開発を続けています。可視数を参照してボリュームデータに対して粒ノード間通信がボトルネックとなって困難化インターフェースとしては、キーフレー子生成を行います。PBVR ではサーフェスだったストロングスケーリングを達成する

ムアニメーションやクロッピングのようなも粒子で表現するため、大規模データに対という成果が得られています [ 図4]。

従来的な基本機能の追加とバージョンアッしても大量のポリゴンを生成することはありません。
粒 子 デ ー タ の サ イ ズ は 通 常 10 ～プを行っています。さらに昨今ますます複商用可視化ソフトとの比較

雑化してゆくシミュレーション結果から知見を得るためのチャレンジングな機能開発

100MB のオーダーであり、ソケット通信

表１に商用可視化ソフト EnSight のクとして、多次元データの可視化機能やマルによりクライアント PC に転送されます。

ライアント・サーバモードによるボリューチスケール可視化手法の開発を進めていま

クライアント PC 上で粒子データはメイン

ムレンダリングとの比較結果を示します。

す。エクサスケールを見据えた更なる要素メモリに格納され、レンダリング時にグラ実 験 に は GT5D デ ー タ を 使 用 し、 実技術の開発として、超高並列環境下での最フィックスコントローラ（GPU）に格納

験 環 境 と し て 帯 域 幅 〜3.4MB/s の イ ー

適化はもちろんのこと、より高画質な可視されます。粒子データは時系列順に逐次転

サ ネ ッ ト で 接 続 さ れ た サ ー バ (Xeon化を可能にする技術を開発する予定です。

PBVREnSight-[sec]3000

フィルタ処理 [sec/step]

0.7

画像生成処理 [sec/step]

51.4

データ転送 [sec/step]75.7

合計 [sec/step]

127.8

3873

1000

描画速度 [frame/sec]

60.0

2.7500クライアントメモリ使用量 [MB]

257

900

0

［表1］全体性能の比較

3873

2500
2000
1500

3264128256[ 並列数 ]

5121024［図4］京による GT5D データの可視化処理時間※ EnSight は米国 Computational Engineering International, Inc. の商標です5ソフトウェアライブラリ

熱流体解析システム

FrontFlow/violet-HC

東京大学生産技術研究所
革新的シミュレーション研究センター特任研究員(Hierarchical Cartesian)

大西順也

FrontFlow/violet Hierarchicalて重要な性質です。なぜなら、工業製品のう課題があります。この課題を解決するこCartesian(FFV-HC) は FrontFlow/設計に流体解析を利用する際、現状ではことを目的として、FFV-HC の開発が始まviolet(FFV)の派生バージョンのひとれらの処理がボトルネックになることが多りました。

つであり、その名の通り、階層型の直交格いからです。FFV では、計算格子作成機子に最適化された熱流体解析システムです。

能をソルバー初期化過程の一部として取り一様等間隔格子の利点を継承しつつ、格子込むことで、大規模な格子を迅速、かつ、配置に柔軟性を持たせることを目的として、自動的に作成することを可能としています。

FFV-HC では、階層型直交格子の導入開発が進められています。

また、将来的には、可視化機能をソルバーにおいて、block-based adaptive mesh内部に取り込むことも検討しています。

refinement(AMR)[2] を利用しています。

[1]

ちなみに、FFV-HC の開発が始まった

階層型直交格子の導入のは2012年9月で、「京」の共用開始と同

図1には、一様等間隔格子を用いた計算

図2に示した通り、block-based AMR で時期です。いわば双子のような関係です。

例として、自動車車体周りの空力解析の結は解析領域全体を八分木構造で分割し、各果を示します。この例では解析領域全体リーフノードにブロックと呼ばれる直交格(24.5m ×12.2m ×6.14m) に4mm の格子を割り当てます。このとき、ブロック内子を配置しました。特筆すべき点は、車体の直交格子には一様等間隔のものが用いら

一様等間隔格子の利点・課題
直交格子の特徴は、データ構造、および、の外部だけではなく、内部 ( エンジンルーれます。また、それぞれの直交格子の大き計算アルゴリズムが単純であることです。

ム、車室など ) にも格子が配置されている

さは物理的には異なるものの、論理的にはこのことは、多方面で様々な利点をもたらところです。なぜなら、従来の格子作成手同じになるようにします ( 計算セル数を同します。

法では、エンジンルームのような複雑な形じにするという意味です )。

状を扱うことは不可能ではないにしても、▶ 計算プログラムの生産性、保守性、可搬

これらの性質は、現代のスパコン事情にとって、非常に都合の良いものと言えます。

非常に難しいからです。
一方で、この解析で用いた計算セルの総例えば、一様等間隔格子を用いることはメ▶ 計算プログラムの単体性能、並列性能が数は、約290億個に達します。非構造格

モリ使用量の削減につながりますし、メモ良好子を用いた解析では、同程度の格子規模で、リアクセスが連続になることから、キャッ▶ 物理現象のモデリングが容易格子間隔を1mm 以下にまで小さくできてシュの有効活用という点でも有利です。ま▶ プリ・ポスト処理が簡便。ソルバーと一いることを考えると、直交格子がいかに非た、各ブロックの計算セル数が同じである体的な実行が可能。

効率的であるかが分かります。つまり、直という点は、プロセスに割り当てるブロッ交格子には、格子配置の柔軟性が低く、必ク数を揃えるだけで計算負荷の均衡を実現要な計算セルが膨大になりがちであるといできることから、分散メモリ環境において性が良好

特に最後の項目は、工学的な応用におい

Cartesian grid

［図2］block-based AMR［図1］290億セルによる計算例
協力：スズキ株式会社

項目

時間 [ 秒 ]１．ポリゴンの読込

85.14

２．
八分木作成・通信

85.61

３．ポリゴン通信

14.03

４．格子—ポリゴン交差判定

14.04

５．計算セル分類

14.16［表1］格子作成所要時間の内訳
6

計算工学ナビ Vol.7

高い並列性能を実現するために重要となり力解析に、階層型格子を適用した結果です。

である FFV-C のように、様々な機能が実ます。

このとき、最小格子間隔が一様直交格子と装されていたり、豊富な応用事例が蓄積さ図3は FFV-HC の計算速度性能を評価

同じ (4mm) になるように階層型格子を作れていたりするわけではありません。加えした結果です。横軸は計算プロセス数で、成したところ、総計算セル数は約4.0億にて、マニュアル等の整備についても、不十何台の計算機を並列動作させたかを示しなりました。つまり、階層型格子を導入す分な点が残っています。また、自動車の例ています。縦軸は FFV-HC の計算速度で、ることで、計算セルを約72分の1にまでで示した通り、階層型格子の導入で計算セ1秒当りに処理した演算数のプロセス間平削減することができたのです。また、最小ル数の削減には成功しましたが、それだけ均値を示しています。なお、各プロットは、格子間隔を半分 (2mm) にしたところ、総ですぐに実用に結びつくわけではありませブロックに割り当てる計算セル数による違計算セル数は約32億となり、それでもまん。より実用的な流体解析を実施するためいを示しています (G5から G1の順に多くだ9分の1程度に抑えられることも分かりには、例えば以下のような、様々な要素技なります )。また、このテストはいわゆるました。

術の開発が必要と考えています。

ウィークスケーリングを調査したものです。

表1には、約32億セルの格子を作成すしたがって、計算速度が計算プロセス数にるのに要した時間をまとめました。ポリゴ▶ 直交格子上で物体形状を正確に表現する依らず、高い水準で一定になるのが理想的

ン ( 車体形状を表す STL ファイル。ファための形状近似機能な結果です。

イルサイズは約3.3GB、ポリゴン数は約▶ 流体 - 構造連成や気液二相流解析のため

FFV-HC の性能はブロックに割り当て1240万 ) の読み込みと八分木の作成・通の移動境界捕獲機能

る計算セル数によって大きく異なります

信がホットスポットになっていることが分▶ 回 転 体 を 考 慮 す る た め の Multipleが、概して言えば、ブロック当りの計算セかります。しかしながら、いずれの処理時

Reference Frame 機能

ル数が多いほど性能が良いと言えます。特間もたかだか数分のオーダーであること、▶ 物体表面の解像度を向上させるためのレに、ブロック当りの計算セル数が最大の場

また、格子作成全体を通しても10分以内

イヤー格子機能合 (G1) は、単体性能がピーク性能比 ( 京に処理が完了することも見て取れます。しでは128GFLOPS) で約10.2%、並列化たがって、従来の手作業による格子作成に

現在は、主に最初の項目に関して研究・効率が32,768ノード使用時で約72.8%比べれば、格段に速くなっていると言えま開発を進めています。そして、HPCI 戦略となり、非常に良い性能を示すことが分かす。

プログラム分野４・課題３の活動の一環とります。逆に、ブロック当りの計算セル数が少なくなると、性能は急激に低下することも見て取れます。この原因は、ブロック間のデータ同期にかかるコストの増加です。

して実施している車体空力とエンジン冷却今後の展望のトレードオフ解析を通して、FFV-HCの実用性を実証することに取り組んでいまFFV-HCは、オープンソース・ソフトウェ
アとして公開されており、いつでも試用が階層型直交格子の導入効果
図4は、先に示したのと同様の自動車空可能な状態となっています

す。次のステップでは、気液二相流解析への適用を進める予定です。

。しかしな

[3]

がら、FFV-HC は開発が始まってから日が浅いため、もうひとつの派生バージョン［図4］階層型直交格子
協力：スズキ株式会社［図3］FFV-HC の単体・並列性能[1] 小野, 三次元非定常熱流体解析プログラム FrontFlow/violet, 計算工学ナビ Vol.5.
[2] 国内では直交格子積み上げ法(Building-Cube Method, BCM)と呼ばれることが多いです[3] FFV-HC ソースコード配布サイト, https://github.com/avr-aics-riken/FFV-HC

7

ｖっｖ v

今号の表紙

GT5D データの静電ポテンシャルの可視化「京」で実行した核融合プラズマ乱流シミュレーションから得られた静電ポテンシャルの値を遠隔可視化システム PBVR で可視化したもの。従来手法
ではデータ転送と可視化作業に長時間必要ですが、PBVR によりインタラクティブに京のストレージ上のデータを可視化できます。図では従来的なスカラー場のボリュームレンダリングに加え、PBVR の多変量データ解析機能によって座標値とスカラー場で不透明度と色を与えて抽出した断面を合成表示しています。
日本原子力研究開発機構 河村拓馬

編集後記

この号の表紙のテーマであるプラズマ乱流シミュレーションについては、本紙 vol.2にも記事があります。京で99.9999％の強スケーリング性能を実現した開発者たちへのインタビューです。現在も下記サイトで公開中ですので、併せてご覧ください。（F）

計算工学ナビ オフィシャルサイト
本誌の PDF 版やソフトウェアライブラリ、ニュースなどのコンテンツを提供している Web サイトは下記の URL からアクセスできます

http://www.cenav.org/

計算工学ナビ Vol.7
発行日：2015年6月25日
発

行：東京大学生産技術研究所
革新的シミュレーション研究センター〒153-8505 東京都目黒区駒場4-6-1
office@ciss.iis.u-tokyo.ac.jp