学際大規模情報基盤共同利用・共同研究拠点 平成 26 年度共同研究 最終報告書 2015 年 5 月

14-NA07

階層分割型数値計算フレームワークを用いた波源から地上構造物までの実地形津波解析
室谷浩平（東京大学）
概要 2011 年 3 月 11 日の東日本大震災によって引き起こされた津波により，東日本は甚大な被害にみまわれた．構造物への被害は，水圧によるだけでなく，浮遊する瓦礫が構造物にぶつかり構造物に被害を及ぼすケースが多くあった．我々は，津波による構造物への被害を検討する際には，津波の遡上計算を行うだけでは不十分と考えており，浮遊物の衝突を含んだ構造解析が必要であると考えている．平成 25 年度は，津波が地上構造物へ及ぼす影響を解析するために 3 段階にわけて解析を行うことに成功した．平成 26 年度
は，より現実に近い浮遊物の漂流解析を行うために，粒子法による流体解析と有限要素法による構造解析の連成計算を行えるフレームワークの開発を目指す．有限要素法による弾性体解析を行うことで，定量的な応力評価が可能となり，津波動的荷重と浮遊物衝突の影響を考慮した沿岸構造物の地震時健全性評価が可能となる．さらに，分散メモリ環境で高いパフォーマンスで実行できるチューニングも行う．

1.

F) 和田義孝（近畿大学 理工学部 機械工学科）：

共同研究に関する情報

VSCG ライブラリの開発(1) 共同研究を実施した拠点名G) 越塚誠一（東京大学大学院工学系研究科シス名古屋大学情報基盤センター

テム創成学専攻）：津波の沿岸から遡上の解析(2) 共同研究分野
のための物理モデルの構築

超大規模数値計算系応用分野

2.

研究の目的と意義

2011 年 3 月 11 日の東日本大震災によって引き(3) 参加研究者の役割分担

A) 室谷浩平（東京大学大学院工学系研究科シス起こされた津波により，東日本は甚大な被害にみ

テム創成学専攻）：津波の沿岸から遡上の解析まわれた．構造物への被害は，水圧によるだけでのための物理モデルの構築，構造流体連成機なく，浮遊する瓦礫が構造物にぶつかり構造物に能の開発被害を及ぼすケースが多くあった．我々は，津波B) 荻野正雄（名古屋大学情報基盤センター大型による構造物への被害を検討する際には，津波の計算支援環境研究部門）：階層領域分割型ソル遡上計算を行うだけでは不十分と考えており，浮バーライブラリの開発，VSCG ライブラリの遊物の衝突を含んだ構造解析が必要であると考え開発ている．

C) 塩谷隆二（東洋大学総合情報学部総合情報学

本研究では，津波が地上構造物へ及ぼす影響を科）：階層領域分割型ソルバーライブラリの開解析するためのマルチスケールマルチフィジック

発，連続体力学向け DSL の開発，地上構造物ス解析を 3 つに分けて行う．第 1 の解析は，震源の解析のための物理モデルの構築で発生する波源から沿岸部までの津波伝播計算でD) 金山寛（名古屋大学情報基盤センター大型計ある．この解析は，数十～数百キロ四方の広範囲

算支援環境研究部門）：津波の波源から沿岸まな解析が必要となるため，粘性浅水波方程式を有での解析のための物理モデルの構築限要素法で解くことにする．第 2 の解析は，沿岸E) 河合浩志（諏訪東京理科大学システム工学部に押し寄せた津波が地上へ遡上し浮遊物を運搬部）：連続体力学向け DSL の開発，構造流体する計算である．この解析は，遡上計算と浮遊物連成機能の開発の取り扱いが容易な粒子法を用いて計算を行う．
1

学際大規模情報基盤共同利用・共同研究拠点 平成 26 年度共同研究 最終報告書 2015 年 5 月

粒子法は，解析領域に比例する粒子数を必要とす

で 26％まで出るようになった[a10,b1-6]．

るため，数 km 四方程度の中規模な範囲の解析し④ SC14 名大ブースでの成果公開サポートか行うことができない．第 3 の解析は，構造物がSC14 の名大ブースで本 JHPCN 課題の成果を受ける水圧と浮遊物の衝突による構造解析である． ポスターとフル HD モニターで展示し，成果公開
この解析は，複数材料の弾塑性非線形解析であるに協力してもらえた[b1]．

ため，有限要素法を用いて計算を行う．

4. 前年度までに得られた研究成果の概要

我々は，津波が地上構造物へ及ぼす影響を解析するために，それぞれの解析ステージで適した解平成 25 年度は，津波が地上構造物へ及ぼす影響析手法を用いる．これらの異なる 3 つの解析手法を解析するために 3 段階にわけて解析を行うことを効率的に組み合わせるために，我々は「階層分に成功し，津波に運搬された多数の浮遊物がどの割型数値計算フレームワーク」を開発し，そのフ様に市街地に押し寄せるかを実現できるようになレームワークを用いて，それぞれの解析ステージ

った[a3,4]．

において高性能で高精度な解析を実現する．

5. 今年度の研究成果の詳細

今年度は，これまで解析を行ってきた石巻市に加え，新たに気仙沼湾の解析を行う．構造物は，現在までに，次の項目が完了している．

昨年度までは剛体として扱ってきたが，今年度は① 分 散 メ モ リ 並 列 EMPS ラ イ ブ ラ リ
LexADV_EMPS v0.1b[e1]の公開【中間報告書

線形弾性体として扱う．

記載済】

3. 当拠点公募型共同研究として実施した意義

塩谷 CREST プロジェクトで開発が進められて共同研究先である名古屋大学情報基盤センター

いる LexADV シリーズの一つとして，2014 年 10からは次のサポートを受けた．

月に LexADV_EMPS v0.1b を① 大規模高精細可視化サポート【http://adventure.sys.t.u-tokyo.ac.jp/lexadv】

これまでは，1 億粒子を超える解析結果は，メに公開した．LexADV_EMPS v0.1b は，高並列なモリ不足のため可視化結果を間引いて可視化して分散メモリ環境での粒子法ソルバー開発支援を目いた．名大の大規模共有メモリ計算機 UV2000 に的としており，「ParMETIS を用いた動的負荷分

リモートディスクトップ NVC を利用することで，散機能」，「隣接間通信のための Halo 領域交換機解析結果を間引くことなく可視化することが可能能」，「高次精度空間離散化モデルのマクロ関数」，となった．また，16 面 8K モニターを利用しての「流体剛体連成機能」の機能を有している．今後

大規模高精細可視化を行うことができるようにな

定期的に本課題の成果を取り込み，機能を拡充しった．ていく予定である．第 9 回 ADVENTURE 定期セ② JHPCN-DF の導入

ミナーや学会発表[c1-4,d1-7]で，LexADV_EMPSJHPCN 採択課題（jh140004-NA02）である萩の紹介を行った．

田（防衛大学）・荻野(名古屋大学情報基盤センター)らの研究成果 JHPCN-DF を利用して，解析結② 「MPS 半陰解法を開発し，MPS 陽解法とのト果出力ファイルの大幅な圧縮に成功した[c5]．

レードオフを検討する．」【中間報告書から進捗③ MPS 陽解法コードのチューニングサポートあり】

昨年度までは，カーネル関数のピーク性能比が，陽解法と半陰解法の精度を検証しつつ，計算時CX400 の Xeon E5-2697 v2 で 7％程度だったのが， 間やメモリ使用量を比較した．
最 大 で 18 ％ ま で 出 る よ う に な り ， FX10 のMPS 陽解法と MPS 半陰解法の圧力分布を比較

SPARC64TM IXfx で 9％程度であったのが，最大するために，流体圧力分布(図 1)と構造解析でメッ2

学際大規模情報基盤共同利用・共同研究拠点 平成 26 年度共同研究 最終報告書 2015 年 5 月

シュが受け取る圧力分布(図 2)を比較した．本解析

た．半陰解法は，並列効率が悪い上，行列の保持では，粒子サイズ：メッシュサイズを 1:4 に設定にメモリを多く使うため，1 ノード当たりの割りしている．図 1(a)の陽解法では圧力振動が激しく当てることができる粒子数が 1/10 程度に少なく起こっているが，図 1(b)の半陰解法では圧力振動なる．

は起こっていない．本来なら半陰解法を用いたい

図 4 は，1 ノード当り 851,402 粒子の 1 ステッところであるが，図 2 のメッシュ面が受け取る圧プ当たりの計算時間である．1 ノードから 240 ノ

力分布を比較すると図 2(a)の陽解法と図 2(b)の半ードまで比較を行い，240 ノードでは 204,336,469

陰解法の差はかなり小さくなる．これは，メッシ粒子の解析となった．ポアソン方程式の解法は，面に所属する壁粒子の平均圧力を，メッシュ面が対角スケーリング前処理 BiCGSTAB 法を用い，受け取るため，壁粒子の圧力が平均化されるため収束判定値は 1e-6 とした．1 ノードでは，陰解法である[a2,6-9]．の計算時間は陽解法の 35 倍のであった．ノード数が増えるにつれて計算時間の差は開き，240 ノードでは 147 倍にまで増えた．これはポアソン方程式の次元数が増え，反復解法の反復回数が増えるためである．半陰解法は，約 10 倍メモリの使用量と，100 倍以上の計算時間が必要であるので大規模な解析を行う際には，計算コストの面から陽解(a)陽解法法の方が推奨される．

(b)半陰解法
図 1．静水圧解析の流体圧力分布

図 3．MPS 陽解法と MPS 半陰解法のストロングスケーリング(a)陽解法(b)半陰解法

図 2．静水圧解析でメッシュ面が受け取る圧力分布

一方で，東大 FX10 上で陽解法と半陰解法のストロングケーリングと計算コストを測定した．図3 は，12 ノードから 4800 ノードのストロングスケーリングの結果である．陽解法は省メモリであ

図 4．MPS 陽解法と MPS 半陰解法の計算時間るため，314,572,800 粒子のストロングスケール

以上結果，構造物が受ける流体圧力を考えた場を行うことができ，半陰解法は行列を保存するた合，陽解法の圧力振動はある程度抑えられること

め，40,867,293 粒子のストロングスケールを行っが確認でき，さらに，問題が大規模になるにつれ3

学際大規模情報基盤共同利用・共同研究拠点 平成 26 年度共同研究 最終報告書 2015 年 5 月て陽解法と半陰解法の並列性能と計算時間の差が④ 大規模高精細可視化と可視化用データの圧縮大きく開くことが確認できため，本研究課題では，従来通り MPS 陽解法を用いることにした[c2,d4]．

【中間報告書記載済】

一方で，今回は，対角スケーリング前処理

名古屋大学情報基盤センターのサポートを受けBiCGSTAB を用いたが，DDM の実装が完了し，て，図 7 の様に，UV2000 を利用した 8K モニタ

今後パフォーマンス測定を行う予定である[d8]．

ーでの可視化に成功した[d6-7]．JHPCN-DF を利更に，実装中の DDM-BDD や AMG を試し，半陰用することで，ほぼ同程度の品質の可視化を実現

解法の高速化に努める予定である．

しつつデータを 60％前後圧縮することに成功した[c5]．

③ 「HDDM 構造解析：1000 万要素のモデルで物理時間 15 分間の解析を行う．
」【中間報告書から進捗あり】
石巻市街地モデルを地上構造物と地盤を合わせたメッシュを作成し，流体の動的荷重を地上構造物へ渡す有限要素解析を行った．図 5 はメッシュ
サイズ 2m の四面体メッシュである．本メッシュ
は 1000 万要素となった．図 6 は流体圧力を受けて ADVENTURE_Solid で構造解析[a1]をした結図 7．16 面 8K モニターによる可視化

果 を相当 応力 で表示 した結 果で ある[a2-7,c3-4,d1-4]．

⑤ LexADV_EMPS のチューニング【中間報告書

更に，構造解析と流体解析の計算ノードの割当から進捗あり】てを，FX10 の Tofu インターコネクトの幾何的なLexADV_EMPS のチューニングを行い，大幅に情報と通信実測値を用いて最適配置を求めた[d2]． 性能が改善された．図 8 は，それぞれのプロセッサ/アクセラレータでのピーク性能比を示した結果である．昨年度までは，カーネル関数のピーク性能比が，CX400 の Xeon E5-2697 v2 で 7％程度だったのが，最大で 18％まで出るようになり，
FX10 の SPARC64TM IXfx で 9％程度であったのが，最大で 26％まで出るようになった[a10,b1-6]．

図 5．メッシュモデル

図 8．各プロセッサ/アクセラレータでのピーク性能比⑥ 気仙沼湾と市街地モデルの作成と津波遡上解
析【中間報告書から大幅進捗あり】
空測量データを国土地理院からレンタルし，航空写真を国際航業株式会社から購入し，海上保安

図 6．地上構造物の相当応力4学際大規模情報基盤共同利用・共同研究拠点 平成 26 年度共同研究 最終報告書 2015 年 5 月

庁から海底測定データ(海図)を購入して，気仙沼湾と市街地モデル(図 9）を作成した．全長 60 メートル船舶(図 10)を 3ds MAX を用いて作成した．

(a) 震源付近
図 9．気仙沼湾と市街地モデル(b) 気仙沼湾と市街地

図 13．TSUNAMI-K による解析結果

図 10．全長 60 メートル船舶モデル

津波波源から沿岸部までを構造計画研究所のソフトウェア TSUNAMI-K を用いて浅水長波解析を行っ

図 14．宮城北部沖の TSUNAMI-K による解析結果

た．断層パラメータに藤井佐竹の 55 小断層モデル

図 15 と図 16 は気仙沼湾と市街地解析のVer. 8.0(図 11)を用い観測点を図 12 の様に設定しLexADV_EMPS よる解析結果である．本解析は，た．
図 13 は TSUNAMI-K による解析結果である．

解析領域 4km×2km を直径 1m の粒子でモデル化

図 13 (b)の気仙沼湾周辺では，TSUNAMI-K 標準し，30 分の解析で最大 2 億粒子になった．本解析のデータベースは 500m 間隔の水深データを用いを行うために，本研究課題により提供を受けた名ているため，気仙沼湾周辺の水深データを，海図

古屋大学情報基盤センターFX10 を 96 ノード用いを元に手入力で修正を行った．図 14 は宮城北部沖て 1 ヶ月かかった．

における GPS 観測結果と TSUNAMI-K による解

6.

析結果の波高を比較した結果である．
TSUNAMI-K による解析結果が GPS 観測結果に今年度の進捗状況と今後の展望

今年度の進捗状況は 80％である．まず，JHPCNの開発項目を取り入れた LexADV_EMPS を 2014

比べて波高が若干低めの結果になった．

年度 10 月に公開できたのが最大の成果であった．
更に，LexADV_EMPS を用いた MPS 半陰解法の実装と精度検証とパフォーマンス測定を行い，現在公開準備中である．今後，DDM-BDD や AMGを導入することで，半陰解法の高速化に努める予定である．
名古屋大学情報基盤センターとの共同研究では，
UV2000 と 8K モニターによる大規模高精細可視化サポート，JHPCN 採択課題（jh140004-NA02）
図 11．断層パラメータの研究成果 JHPCN-DF の導入によるデータ圧縮

図 12．測定点5学際大規模情報基盤共同利用・共同研究拠点 平成 26 年度共同研究 最終報告書 2015 年 5 月

Publishing, Vol. 33, pp. 171-186, 2014 (DOI:10.1007/978-3-319-06136-8_8).
[a6] N. Mitsume, S. Yoshimura, K. Murotani and T.
Yamada : Improved MPS-FE Fluid-StructureInteraction Coupled Method with MPS Polygon
Wall Boundary Model, Computer Modeling inEngineering & Sciences (CMES), Vol.101, No. 4,
pp. 229-247, 2014
(DOI:10.3970/cmes.2014.101.229).
[a7] Naoto Mitsume, Shinobu Yoshimura, Kohei
Murotani, Tomonori Yamada: ExplicitlyRepresented Polygon Wall Boundary Model for
Explicit-MPS Method, Computational Particle
Mechanics(CPM),
13May
2015.
(DOI:10.1007/s40571-015-0037-8)[a8] Kazuya Shibata,Issei Masaie,Masahiro
Kondo, Kohei Murotani, Seiichi Koshizuka:
Improved pressure calculation for the moving
particle semi-implicit method, Computational
Particle Mechanics (CPM), 08 May 2015.
(DOI:10.1007/s40571-015-0039-6)[a9] Kazuya Shibata, Seiichi Koshizuka, Kohei
Murotani, Mikio Sakai, Issei Masaie: Boundary
Conditions for Simulating Karman Vortices
Using the MPS Method, Journal of Advanced
Simulation in Science and Engineering(JASSE), accepted.
[a10] Kohei MUROTANI, Issei MASAIE, Takuya
MATSUNAGA, Seiichi KOSHIZUKA, Ryuji
SHIOYA,Masao
OGINO,Toshimitsu
FUJISAWA: Performance Improvements ofDifferential Operators Code for MPS method
on GPU, Computational Particle Mechanics(CPM), accepted.

手法の導入，カーネル関数のピーク性能比の大幅改善に大きな協力があった．
今年度の未達成項目は，気仙沼解析の流体構造連成解析である．解析上の問題点が幾つか明らかになったため，次年度に研究開発が引き続き行う予定である．

7.

研究成果リスト(1) 学術論文[a1] A.M.M. Mukaddes, Masao Ogino, and Ryuji
Shioya, “Performance Evaluation of Domain
Decomposition Method with Sparse Storage
Schemes
inModern
Supercomputer”,
International Journal of ComputationalMethods, Vol. 11, Suppl. 1, 1344007, pp.1-14,
2014 (DOI: 10.1142/S0219876213440076).
[a2] N. MITSUME, S. YOSHIMURA, K.
MUROTANI and T. YAMADA: MPS-FEM
PARTITIONED COUPLING APPROACH FOR
FLUID-STRUCTURE INTERACTION WITH
FREE SURFACE FLOW, International Journalof Computational Methods (IJCM), Vol. 11,
issue 4, pp.1350101, 16 pages, 2014 (DOI:
10.1142/S0219876213501016).
[a3] Kohei Murotani, Seiichi Koshizuka, Tasuku
Tamai, Kazuya Shibata, Naoto Mitsume,
Shinobu Yoshimura, Satoshi Tanaka, Kyoko
Hasegawa, Eiichi Nagai and Toshimitsu
Fujisawa:
Developmentof
HierarchicalDomain Decomposition Explicit MPS Methodand Application to Large-scale Tsunami
Analysis with Floating Objects, Journal ofAdvanced
SimulationinScienceandEngineering (JASSE), Vol. 1, No. 1, pp.16—35,
2014 (DOI:10.15748/jasse.1.16).
[a4] Kentaro Tanaka, Satoshi Tanaka, Kyoko
Hasegawa, Kohei Murotani, Seiichi Koshizuka:
Translucent Visual Analysis of Large Scale 3D
Point Data Generated by Particle Fluid,
Journal of Advanced Simulation in Science andEngineering (JASSE), Vol. 1, No. 1, pp.5—15,
2014 (DOI: 10.15748/jasse.1.5).
[a5] Yokoyama, M., Murotani, K., Yagawa, G., andMochizuki, O., Some Considerations on Surface
ConditionofSolidin
Computational
Fluid-Structure Interaction. In Numerical
SimulationsofCoupledProblemsinEngineering (Computational Methods inApplied Sciences), Springer International(2) 国際会議プロシーディングス[b1] Kohei Murotani, Seiichi Koshizuka, Masao
Ogino, Ryuji Shioya, Yasushi Nakabayashi:
Development of distributed parallel explicit
Moving Particle Simulation (MPS) method andzoom up tsunami analysis on urban areas,
SC14, poster session, New Orleans, Nov. 16-21,
2014.
[b2] Kohei Murotani, Seiichi Koshizuka, Naoto
Mitsume, Shinobu Yoshimura, Kazuya Shibata,Tasuku Tamai, Satoshi Tanaka, Kyoko
Hasegawa, Toshimitsu Fujisawa: Large-Scale
Tsunami Analysis on Urban Area using
MPS-FE Fluid-Structure Interaction Coupled
Method, The 14th Asia Simulation Conference& The 33rd JSST Annual Conference:6学際大規模情報基盤共同利用・共同研究拠点 平成 26 年度共同研究 最終報告書 2015 年 5 月

International Conference on Simulation
Technology, Kitakyushu, Japan, October 26-30,
2014.

層型領域分割 MPS 陽解法を用いた多数の浮遊物が漂流する市街地津波遡上解析，第 19 回計算工学講演会，広島国際会議場所, 2014 年 6 月 11-13日.
[d2] 室谷浩平，越塚誠一，柴田和也，三目直登，吉村忍，藤澤智光：階層型領域分割 MPS-FEM 流体構造連成解析のための最適なロードバランスの検討，日本機械学会 第 27 回計算力学講演会（CMD2014），岩手大学，2013 年 11 月 22-24日．
[d3] 室谷浩平: 大規模分散メモリ並列 MPS 法の開発と市街地浸水津波解析, 先駆的科学計算に関するフォーラム 2014, 福岡, 九州大学情報基盤研究開発センター, 2014 年 8 月 5-6 日.
[d4] 室谷浩平: 分散メモリ並列 MPS 陽解法・半陰解法の圧力振動抑制, 計算精度と計算コストの比
較, ワークショップ 自由表面や気液界面を含む流れの数値解析 ―粒子法と格子法, 数理と実践―, 福岡, リファレンス駅東ビル, 2014 年 9 月 10日.
[d5] 室谷浩平: 粒子法シミュレーションコードの実
装 を通 した Intel x86 CPU, Fujitsu Sparc64
CPU, Intel Xeon Phi, GPU コンピューティングの比較, Prometech Simulation Conference 2014,東京コンファレンスセンター・品川, 2014 年 9 月26 日.
[d6] 室谷浩平: スーパーコンピュータを用いた大規模粒子法解析とその可視化, 第 20 回ビジュアリ
ゼ ー シ ョン カン フ ァレ ンス , タ イム 24,東 京 ,2014 年 11 月 7 日.
[d7] 室谷浩平，越塚誠一，塩谷隆二, 荻野正雄 : 分
散メモリ並列 EMPS ライブラリ LexADV_EMPSの使い方と大規模可視化システム UV2000 上での ParaView, AVS, Ensight による大規模粒子系シミュレーション結果の可視化, 第 30 回粒子法
コードユーザグループ会合, 海上技術安全研究所
本館 A 会議室,2014/11/26.
[d8] 本橋佑一, 室谷浩平, 柴田和也, 越塚誠一，荻野正雄，田上大助：圧力ポアソン方程式への部分構
造反復法の実装とその性能検証，日本機械学会
第 27 回計算力学講演会（CMD2014），岩手大学，2013 年 11 月 22 日-24 日．

(3) 国際会議発表[c1] Kohei MUROTANI, Seiichi KOSHIZUKA,Hiroshi KANAYAMA, Kazuya SHIBATA,Tasuku TAMAI, Naoto MITSUME, Shinobu
YOSHIMURA, Satoshi TANAKA, Kyoko
HASEGAWA, and Toshimitsu FUJISAWA:
Large-scale run-up tsunami analysis withmany floating objects on urban area by threeanalyses stages using hierarchical domain
decomposition in explicit MPS method, 1stInternational Conference on Computational
Engineering and Science for Safety andEnvironmental Problems (COMPSAFE 2014),Sendai, Japan , 13-16 April 2014.
[c2] (Invited Lecture) Kohei MUROTANI, Seiichi
KOSHIZUKA: Comparisons of calculation costand accuracy between the explicit andsemi-implicit distributed memory parallel MPS
method (Invited), 5th International Conference
on Computational Methods (ICCM2014),
Fitzwilliam College Cambridge, England,
28-30th July 2014.
[3] (Invited Lecture) Kohei MUROTANI andSeiichi KOSHIZUKA: Zoom up tsunami
analysis with many floating objects on urban
areas by three analyses stages using
distributed parallel explicit MPS method, 5thAICS International Symposium, Kobe, Japan,
December 8--9, 2014.
[c4] Kohei Murotani, Ryuji Shioya, Seiichi
Koshizuka, Masao Ogino, Hiroshi Kawai:
Development of a Numerical Library based on
Hierarchical Domain Decomposition for Post
Petascale Simulation, SC14, ATIP Workshop
onJapanese
ResearchToward
Next-Generation Extreme Computing, NewOrleans, Nov. 16-21, 2014.
[c5] Katsumi Hagita, Manabu Omiya, Takashi
Honda, Kohei Murotani, Masao Ogino, Study ofEfficient Data Compression by JHPCN-DF,
Annual Meeting on Advanced Computing
System and Infrastructure (ACSI) 2015, つくば国際会議場, 2015 年 1 月 26 日－28 日(5) その他（特許，プレス発表，著書等）[e1]【受賞】日本機械学会計算力学部門 2013 年度優秀講演表彰, 和田義孝, 2014 年 7 月 1 日.
[e2]【著書】越塚誠一，柴田和也，室谷浩平：粒子法
入門 「MPS 法による流体シミュレーションの基礎から並列計算と可視化まで」 （サンプルプ
ログラムのソースコード付き），丸善出版，256ページ，2014 年 6 月 25 日初版出版．

(4) 国内会議発表[d1] 室谷浩平，越塚誠一，玉井佑，柴田和也，三目直登，吉村忍，田中覚，長谷川恭子，藤澤智光: 階7

学際大規模情報基盤共同利用・共同研究拠点 平成 26 年度共同研究 最終報告書 2015 年 5 月[e3]【プログラム登録】著作者：室谷浩平，越塚誠一，
塩谷隆二，題号：分散メモリ並列 EMPS ライブ
ラリ LexADV_EMPS v0.1-b，登録番号：P 第10453 号-1, 平成 26 年 10 月 14 日，内容：分散メモリ並列 MPS 陽解法解析プログラム．

(a) 0 s(b) 600 s(c) 1200 s(b) 1792 s

図 15．気仙沼湾と市街地の LexADV_EMPS による解析結果（真上から）(a) 0 s(b) 600 s(c) 1200 s(b) 1792 s

図 16．気仙沼湾と市街地の LexADV_EMPS による解析結果（高台から）

8