Vol.52 No.3

2019-7኱つᶍ⛉Ꮫィ⟬ࢩࢫࢸ࣒㛵㐃᱌ෆ㸺኱つᶍ⛉Ꮫィ⟬ࢩࢫࢸ࣒㛵㐃ᴗົࡣࠊࢧ࢖ࣂ࣮ࢧ࢖࢚ࣥࢫࢭࣥࢱ࣮ᮏ㤋ෆࡢ᝟ሗ㒊᝟ሗᇶ┙ㄢࡀᢸᙜࡋ࡚࠸ࡲࡍࠋ㸼KWWSVZZZVVFFWRKRNXDFMS

㝵

ಀ䞉ᐊྡ

฼⏝┦ㄯᐊ

㟁ヰ␒ྕ䠄ෆ⥺䠅 
HPDLO  
VRGDQ#FFWRKRNXDFMS┦ㄯဨ୙ᅾ᫬

  䝃䞊䝡䝇᫬㛫ᖹ᪥ィ⟬ᶵ฼⏝඲⯡䛻㛵䛩䜛┦ㄯ 䡚኱ุ䝥䝸䞁䝍䚸฼⏝⪅➃ᮎ➼䛾 䡚฼⏝୺䛺䝃䞊䝡䝇ෆᐜ୍㝵ྛ䝉䞁䝍䞊ᗈሗ䛾㜀ぴ฼⏝⪅ㄯヰᐊ

 

ᒎ♧ᐊ ศᩓ䡶䢙䢇䢛䡩䡬䡼༤≀㤋 ぢᏛ䜢䛤ᕼᮃ䛾᪉䛿ඹྠ฼⏝Ṕ௦䛾኱ᆺィ⟬ᶵ➼䛾ᒎ♧ᨭ᥼ಀ䜎䛷䛤㐃⤡䛟䛰䛥䛔䚹䠖䡚䠖ྍど໬ᶵჾᐊ

 ୕ḟඖྍど໬䝅䝇䝔䝮䛾฼⏝䡚

⥲ົಀ  
VRP#FFWRKRNXDFMS⥲ົ䛻㛵䛩䜛䛣䛸䡚

఍ィಀ  
NDLNHL#FFWRKRNXDFMS఍ィ䛻㛵䛩䜛䛣䛸䚸㈇ᢸ㔠䛾ㄳồ䛻㛵䛩䜛䛣䛸䡚

䠄䠅ඹྠ฼⏝ᨭ᥼ಀ
  ཷ㻌 ௜ XNHWXNH#FFWRKRNXDFMS฼⏝ᡭ⥆䛝䚸฼⏝┦ㄯ䚸ㅮ⩦఍䚸䝷䜲䝤䝷䝸䚸ぢᏛ䚸䜰䝥䝸䜿䞊䝅䝵䞁䛻㛵䛩䜛䛣䛸䡚

ඹྠ◊✲ᨭ᥼ಀ

  
UVVHF#FFWRKRNXDFMSඹྠ◊✲䚸ィ⟬ᶵ䝅䝇䝔䝮䛻㛵䛩䜛䛣䛸

䡚

䝛䝑䝖䝽䞊䜽ಀ  
QHWVHF#FFWRKRNXDFMS䝛䝑䝖䝽䞊䜽䛻㛵䛩䜛䛣䛸䡚⮬㈍ᶵ

䡚୕㝵᝟ሗ䝉䜻䝳䝸䝔䜱   ᝟ሗ䝉䜻䝳䝸䝔䜱䛻㛵䛩䜛䛣䛸ಀ
LVHFXULW\#JUSWRKRNXDFMSᅄ◊✲㛤Ⓨ㒊
㝵  ஬㝵 ➃ᮎᶵᐊ

 

3& ➃ᮎᶵ䠄; ➃ᮎ䠅䡚

䡚㸦 㸧ෆࡣᮾ໭኱Ꮫෆࡢࡳࡢෆ⥺㟁ヰ␒ྕ࡛ࡍࠋ㟷ⴥᒣ࣭ᕝෆᆅ༊௨እ࠿ࡽࡣ㢌࡟  ࢆຍ࠼ࡲࡍࠋ
ᮏㄅ䛾ྡ๓䛂㻿㻱㻺㻭㻯䛃䛾⏤᮶㻌᫛࿴  ᖺ䛻ᮾ໭ᆅ༊䛾᭱ึ䛾㟁Ꮚィ⟬ᶵ䛸䛧䛶䚸ᮾ໭኱Ꮫ㟁Ẽ㏻ಙ◊✲ᡤ䛻䛚䛔䛶᏶ᡂ䛥䜜䛯䝟䝷䝯䝖䝻䞁ᘧィ⟬ᶵ䛾ྡ๓䛷 6(1$& 6(1GDL$XWRPDWLF&RPSXWHU 䛛䜙䛸䛳䛶࿨ྡ䛥䜜䛯䚹㻌SENAC Vol. 52, No. （3 2019. 7）― 1 ―［共同研究成果］

Flamelet approach に基づくガス燃料を対象とした燃焼シミュレーション
赤尾津翔大：東北大学大学院工学研究科
松下 洋介：東北大学大学院工学研究科
青木 秀之：東北大学大学院工学研究科
Weeratunge Malalasekera：Loughborough University

本研究では，ガス燃料を対象とした乱流燃焼シミュレーションのターゲット・フレームとして最も実績のある乱流拡散火炎の 1 つである Sandia Flame D を対象に Flamelet approach の 1 つである Flame/Progress-Variable approach を用いて乱流燃焼の Large Eddy Simulation を実施し，測定結果と比較することで本乱流燃焼シミュレーションの妥当性を検討した．その結果，混合分率，温度
およびラジカルを含む化学種の時間平均値とその分散について解析結果は測定結果とほぼ完全に一致した．今後，同手法を拡張することで NO やすすなどの大気汚染物質の予測精度について検討する予定である．

1.

緒言

調整パラメータを必要としない詳細化学反応機構を用いた燃焼シミュレーションは理想的である．しかしながら，最も単純な炭化水素であるメタンの燃焼ですら，例えばこのメタンの燃焼を対象に対して最も実績のある詳細化学反応機構の 1 つである GRI-Mech 3.0 [1]には 53 の化学種と353 の化学反応が含まれる．そのため，詳細化学反応機構を用いて三次元の燃焼シミュレーショ
ンを実施すると，解くべき化学種の保存式の数が多く，計算負荷が膨大となってしまう．また，一般に詳細化学反応機構を用いた反応解析の時間刻みは流体解析の時間刻みより数オーダー小さい．そのため，時間スケールの小さい化学反応の時間刻みに合わせて三次元の燃焼シミュレーシ
ョンを実施すると，やはり計算負荷が膨大となってしまう．さらに，詳細化学反応機構に限らず，有限の反応速度を用いて燃焼シミュレーションを実施する場合，化学種の保存式の生成項にアレ
ニウス型の反応速度が現れる．Reynolds-Averaged Navier-Stokes (RANS)あるいは Large Eddy
Simulation (LES)を用いて乱流燃焼の計算を行う場合，化学種の保存式に時間平均あるいは空間平均を施すが，この化学種の保存式の生成項に現れるアレニウス型の反応速度に時間平均あるいは
空間平均を施す方法論が確立されていない．これはアレニウス型の反応速度式に温度に関する指数関数が含まれ，非線形性が強いためである．LES を用いた乱流燃焼では化学反応速度の項にス
ケール相似則を適用する Scale Similarity Filtered Reaction Rate (SSFRR) [2]モデルを用いた計算が報告されているものの，その妥当性は十分に確認されているとは言えず，また，解くべき化学種の保存式の数が多くなることと時間刻みが小さくなることにより計算負荷が膨大となってしまう
問題は解決されない．
Flamelet approach [3,4]では，拡散燃焼において形成される火炎が一次元的な火炎片(Flamelet)の集合で表現できるとする．また，化学反応の時間スケールが流体の時間スケールと比較して極めて小さいことを利用し，化学反応と流体をスケール分離する．すなわち，流体の時間スケールでは化学反応が定常状態に達しているとし，あらかじめ 2，3 のパラメータに対してデータベースを作成し，三次元の燃焼シミュレーションではこのパラメータを求め，データベースを参照することで興味ある残りのパラメータすべてを決定する．そのため，Flamelet approach を用いることで，比較的低い計算負荷で間接的ではあるものの詳細化学反応機構を考慮した三次元の燃焼シミュレーションを実施することが可能となる．

― 2 ―

SENAC Vol. 52, No. （3 2019. 7）

本研究では，ガス燃料を対象とした乱流燃焼シミュレーションのターゲット・フレームとして最も実績のある乱流拡散火炎の 1 つである Sandia Flame D [5–7]を対象に Flamelet approach の 1 つである Flamelet/Progress-Variable approach [8]を用いて乱流燃焼の Large Eddy Simulation を実施し，測定結果と比較することで本乱流燃焼シミュレーションの妥当性を検討する．

2.

Flamelet approach に基づく燃焼シミュレーション

Flamelet approach では，三次元の燃焼シミュレーションに先立ち，あらかじめ一次元の対向流
拡散火炎や予混合火炎などの単純な火炎を対象に詳細化学反応機構を用いて一次元の燃焼シミュ
レーションを実施し，混合分率 Z や反応の進行を表す Progress Variable (C)と呼ばれるパラメータに対して，温度，化学種の質量分率，密度や粘度などの物性値および C の正味の生成速度を
Flamelet table と呼ばれるデータベースに保存する．三次元の燃焼シミュレーションにおいてこの混合分率 Z や C の保存式を解き，あらかじめ作成した Flamelet table を参照することで，燃焼場における温度，化学種の質量分率，密度や粘度などの物性値および正味の C の生成速度などの変数を決定する．なお，Flamelet table から参照する変数のうち，三次元の燃焼シミュレーションに影響を及ぼす変数は密度や粘度などの物性値および C の正味の生成速度である．

2.1 解析対象と解析条件
本研究では，ガス燃料を対象とした乱流燃焼シミュレーションのターゲット・フレームとして最も実績のある Sandia Flame D を対象とする．図 1 に解析対象と解析領域の概念図を示す．Sandia
Flame D は，開放空間において(a)内側のノズルから温度 294 K，体積分率でメタン 0.25 と空気 0.75の混合気を断面平均流速 49.6 m/s で，(b)その外側からメタンを当量比 0.77 で燃焼させた温度 1880
K のパイロット火炎の燃焼ガスを断面平均流速 11.4 m/s で，(c)周囲から co-flow として温度 291 Kの空気を断面平均流速 0.9 m/s で供給することで形成される部分予混合火炎である．なお，燃料ノズルを基準とすると，レイノルズ数は Re = 22,400 である．混合分率は酸化剤側の境界面で 0，燃料側の境界面で 1 となる変数であり，境界条件として，(a)では Z = 1，C = 0，(b)では Z = 0.27，C= 0.21，(c)では Z = 0, C = 0 となる．なお，計算格子には六面体格子のみを用い，総計算格子数を
128 万あるいは 300 万分割とした．

図 1 解析対象と解析領域の概念図 [5–7]Flamelet approach に基づくガス燃料を対象とした燃焼シミュレーション― 3 ―2.2 Sandia Flame D に対する Flamelet table の作成
一次元対向流拡散火炎を対象に，詳細化学反応機構にメタンの燃焼に対して多くの実績がある
GRI-Mech 3.0 [1]を採用し，FlameMaster V3.3.10 [9]を用いて一次元の燃焼シミュレーションを実施した．熱損失ゼロおよびいわゆる unity Lewis number を仮定し，混合分率空間において種々のス
カラー消散率について一次元対向流拡散火炎の定常解を求め，混合分率 Z とスカラー消散率に対して温度，化学種の質量分率，密度や粘度などの物性値および正味の Progress Variable C の生成速度などの変数で表される二次元のデータベースを保存する．その後，混合分率 Z とスカラー
消散率に対する変数で表される二次元のデータベースを混合分率と Progress Variable C に対する変数で表される二次元のデータベースに変換する．さらに，次節で述べるとおり，Large Eddy
Simulation を用いて三次元の燃焼シミュレーションを実施するため，混合分率 Z の確率密度関数に関数，Progress Variable C の確率密度関数に関数を仮定し，Favre 平均を施した混合分率，その分散と Progress Variable C に対する変数で表される三次元のデータベースを作成する．なお，
Flamelet table の作成方法の詳細は著者らの総説[10]で解説したとおりである．
2.3 三次元の燃焼シミュレーションの解析方法
Sandia Flame D を対象とした三次元燃焼シミュレーションの基礎式は低マッハ数近似に Favre
フィルタを施した連続の式 (1)，運動量保存式 (2)，混合分率の保存式 (3)と Progress Variable (PV)の保存式 (4)である．
��̅����̅ ��� � � ��� ���

�
��
�
��̅
�
���� ���� 2 �����̅ �� � ���̅ ��� ��� � ���̅ ����� ��� �� � �
�
�� �
�
� ���
��
� �� � ������ ��� 3��������� ���������� � ����� � ����� �� ���̅ ��� ���̅ ����� � ���̅ ���
�
�
��
���
���
��� �� ������� � ����� � �� ��� �� ���̅ �� � ���̅ �� ��� � ���̅ ����� � ����
���
��
���
���
��� �� ���

(1)
(2)
(3)
(4)

ここで，式 (2)–(4)の右辺第 3 項をそれぞれ式 (5)–(7)で近似する．
������ ���� 2 ����̅ ����� ��� � � �� �
�
� �
�
(5)
� �� � ���� ��� 3 �� ���� ����� � ����� � � ���̅ ���
(6)����� ���� ����� � �� ��� � � ���̅ ���
(7)����� ���式 (5)–(7)中のt は乱流粘性係数であり，式 (8)に示す Smagorinsky モデル[11]を用いて推算した．
(8)�� � ���� ∆�� �
なお，CS は Smagorinsky モデルの唯一のモデル定数であり，本研究では 0.1 で一定とした．また，ScZ,t と ScC,t,は混合分率と Progress Variable の乱流シュミット数であり，本研究では 0.7 で一定とした．混合分率の分散は式 (9)を用いて推算した．
�
���
�
���
�
� � �� ∆ � �
���

(9)

以上より，求めた混合分率，混合分率の分散と Progress Variable を用いて前述した Flamelet table
を参照し，密度と Progress Variable の正味の生成速度を求めるとともに温度や化学種の質量分率を決定する．
非構造格子の有限体積法に基づき，運動量保存式の対流項は三次風上差分法を 5%ブレンドした二次中心差分法，拡散項二次中心差分法を用いてそれぞれ離散化し，時間進行法には二次のAdams-Bashforth 法を用いた．計算時間を短縮するため時間刻みを可変とし，クーラン数が 0.4 と― 4 ―

SENAC Vol. 52, No. （3 2019. 7）なるように決定した．なお，乱流燃焼場が定常的な挙動を示す際の時間刻みは約 1 s である．圧力の解法には Simplified Marker And Cell (SMAC)法[12]を適用し，SMAC 法における圧力補正値に関するポアソン方程式の解法には Algebraic Multigid Solver [13]を用い，各タイム・ステップにおいて収束解が求められるまで反復計算を実施した．混合分率と Progress Variable の保存式の対流
項は流束制限関数に min-mod 関数[14]を適用した Total Variation Diminishing を，拡散項は二次中心差分法を用いて離散化し，時間進行法には一次の陰解法を用い，多項式前処理付き安定化双共役勾配法[15]を用い，各タイム・ステップにおいて収束解が求められるまで反復計算を実施した．
計算開始から時刻 0.2 s まで乱流燃焼場を発達させ，時刻 0.2–0.3 s のデータを平均化することで各変数の時間平均値を求めた．METIS-5.1.0 [16]を用いて解析領域を分割し，Message Passing
Interface (MPI)を用いて領域分割に基づく並列計算を実施した．また，通信の多い領域間から可能
な限り同時に通信する独自のアルゴリズム[17]を採用し，通信に要する時間を最小化し，オーバ
ヘッドを最小化することで並列化効率を向上した．なお，東北大学サイバーサイエンスセンター
所有の並列コンピュータ LX 406Re-2 を用い，4 ノード 48 コアを用いた並列計算を実施した．

3.

結果と考察

3.1 混合分率，温度および主要な化学種の質量分率分布の瞬時値
図 2 に t = 0.2 s において，z = 0 の断面における混合分率，温度および主要な化学種の質量分率の分布を示す．混合分率は酸化剤側の境界面で 0，燃料側の境界面で 1 となる変数であり，ノズルから供給された燃料の噴流が減衰するとともに周囲のガスと混合することが示されている．温度については，上流において燃料ノズルの周囲から供給された高温のパイロット火炎の燃焼ガスが高温を示し，また，ノズルから供給された燃料がこの燃焼ガスと混合することで加熱されるとともにパイロット火炎の周囲から供給された co-flow と混合することで燃焼が進行し，中流から下流に向かって高温領域が形成されている．また，燃焼の中間生成物である CO は比較的燃料過濃な領域で，OH は比較的燃料希薄な領域で生成している．さらに，燃焼の最終生成物である CO2が中流から下流に向かって生成している．これらの分布は定性的な比較ではあるものの既往の研究[18]と概ね一致している．次節において解析結果を測定結果と詳細に比較する．

図 2 速度，混合分率，温度および主要な化学種の質量分率の瞬時値の分布

Flamelet approach に基づくガス燃料を対象とした燃焼シミュレーション― 5 ―3.2 中心軸上における各変数の時間平均値の比較
図 3 に中心軸上における混合分率，その Root-Mean-Square (RMS)，温度および化学種(CH4, CO,
CO2, OH, H2O, H2)の質量分率の時間平均値を測定結果とともに示す．解析結果および測定結果はいずれも上流から下流に向かって 1 から徐々に減衰している．これはノズルから供給された燃料が拡散しながら周囲のガスと混合するためである．混合分率の時間平均値の解析結果は測定結果とほぼ完全に一致していることから，本解析はノズルから供給された燃料の噴流が乱流燃焼場において減衰し，周囲のガスと混合する過程を良好に表現していることを意味する．
混合分率の RMS については解析結果および測定結果はいずれも上流から下流に向かって上昇
し，その後減少している．解析結果は測定結果を概ね再現しているものの，上流において解析結
果は測定結果より若干低い値を示している．これは，解析では流入境界面において空間的な速度分布を考慮しているものの，その時間的な変動を考慮していないためであると考えられる．理想的には本解析の流入境界面より上流の燃料を供給するノズルを解析対象に含めることで，本解析の流入境界面における時間的な速度の変動も再現すべきである．しかしながら，LES において壁
せん断が支配となるノズル内の乱流流れの解析を実施するためにはさらに多くの計算格子を必要とし，計算負荷が大きくなってしまうため現実的ではない．ノズルの下流に位置するある断面における速度などの物理量をノズルの上流に位置する流入境界面に複製するいわゆる Recycle
boundary を用いることで解析すべきパイプの長さを短くすることはできるものの，それでもなお多くの計算格子を必要とする．一方，Pierce and Moin [19]は事前にノズルのみを対象とした解析を実施し，実際の解析の流入境界面に相当する断面における速度などの空間的な分布とその時間的な変動を含む物理量をすべて保存し，実際の解析でこの保存したデータを流入境界条件として用いることを提案している．この方法では，流入境界面において妥当な空間的な分布と時間的な変動が求められることが期待されるが，保存すべきデータ量が膨大となるだけでなく，異なる解析対象間に生じるメッシュの不整合を処理しなければならない．また，パイプ内の流動が下流の影響を受けないことも確認しなければならないため，用いるためには解決すべき課題も多い．実
用的な方法として，流入境界面において人工的に時間的な変動を生成する手法[20]が数多く提案されている．しかしながら，本解析対象と同じ Sandia Flame D を対象に流入境界面に人工的に時間的な変動を生成する手法を適用した解析結果が報告されているものの，解析対象ごとに調整すべきパラメータが存在し，汎用性に欠ける．そのため，LES において妥当な流入境界条件を与えるのは極めて困難であり，今後の研究課題であると言える．
温度の時間平均値については解析結果および測定結果はいずれも上流から下流に向かって上昇
し，最高温度を示した後に緩やかに減少している．これは，ノズルから供給された燃料がその周囲から供給されたパイロット火炎の高温の燃焼ガスと混合するとともにパイロット火炎のさらに周囲の co-flow と混合し，燃焼するためである．解析結果は測定結果とほぼ完全に一致している
ものの，最高温度を示した後に解析結果は測定結果と比較して若干過大評価している．これは，解析結果と測定結果に差が生じる領域では，燃焼ガスの主成分であり吸収性ガスである CO2 とH2O の質量分率が高く，ふく射伝熱により熱損失を考慮していないためであると考える．実際，
本解析対象と同じ Sandia Flame D を対象とした研究においても同様の傾向が報告されており[18]，
今回用いた燃焼モデルの範疇ではむしろ妥当な解析結果が求められていると考える．しかしながら，NO など温度に敏感な化学種の質量分率を正確に予測するためには温度をさらに正確に予測する必要があり，ふく射伝熱を連成し，Enthalpy defect [21]と呼ばれる熱損失を考慮可能な Flamelet
approach を実施する必要があると考える．

4.

結言

本研究では，ガス燃料を対象とした乱流燃焼シミュレーションのターゲット・フレームとして最も実績のある乱流拡散火炎の 1 つである Sandia Flame D を対象に Flamelet approach の 1 つであ― 6 ―

SENAC Vol. 52, No. （3 2019. 7）

る Flame/Progress-Variable approach を用いて乱流燃焼の Large Eddy Simulation を実施し，測定結果と比較することで本乱流燃焼シミュレーションの妥当性を検討した．その結果，混合分率，温度
およびラジカルを含む化学種の時間平均値とその分散について解析結果は測定結果とほぼ完全に一致した．今後，同手法を拡張することで NO やすすなどの大気汚染物質の予測精度について検討する予定である．

(a) Z とその RMS(b) 温度(c) CH4(d) CO(e) CO2(f) OH(g) H2(h) H2O
図 3 中心軸上における混合分率，その RMS，温度および化学種の質量分率 [5–7]謝辞
本研究は，東北大学サイバーサイエンスセンターの並列コンピュータ LX 406Re-2 を利用することで実現することができた．また，本研究の一部は特別研究員奨励費(18J11135)，JSPS 科研費(JP18K03964)，東燃ゼネラル石油研究奨励・奨学財団および東北大学若手リーダー研究者海外派
遣プログラムの助成を受けたものであり，ここに謝意を表す．
参考文献[1] Smith, G. P. et al., GRI-Mech 3.0, http://www.me.berkeley.edu/gri_mech/[2] DesJardin, P. E. and Frankel, S. H., Large eddy simulation of a nonpremixed reacting jet: application
and assessment of subgrid-scale combustion model, Phys. Fluids A, 10, 2298–2314 (1998)[3] Peters, N., Laminar diffusion flamelet models in non-premixed turbulent combustion, Prog. Energy
Combust. Sci., 10(3), 319–339 (1984)

Flamelet approach に基づくガス燃料を対象とした燃焼シミュレーション― 7 ―[4] Peters, N., Laminar flamelet concepts in turbulent combustion, Symp. (Int.) on Combust., 21, 1231–
1250 (1998)[5] Barlow, R. S. and Frank, J. H., Effects of turbulence on species mass fractions in methane/air jet
flames, Proc. Combust. Inst. 27, 1087–1095 (1998)[6] Barlow, R. S. et al., Piloted methane/air jet flames: Scalar structure and transport effects, Combust.
Flame, 143, 433–449 (2005)[7] Schneider, Ch. et al., Flow field measurements of stable and locally extinguishinghydrocarbon-fuelled jet flames, Combust. Flame, 135, 185–190 (2003)[8] Pierce, C. and Moin, P., Progress-variable approach for large-eddy simulation of non-premixed
turbulent combustion, J. Fluid Mech., 504 73–97 (2004)[9] Pitsch, H., and Bollig, M., FlameMaster; A Computer Code for Homogeneous Combustion and
One-Dimensional Laminar Flame Calculations, Institut für Technische Mechanik, RWTH, Aachen.
[10] 松下洋介ら, Flamlelet Model に基づく乱流燃焼シミュレーション, 金属, 85(11), 915–921(2015)[11] Smagorinsky, J., General circulation experiments with the primitive equations, Mon. Weather Rev.,
91(3), 99–164 (1963)[12] Harlow, F. H. and Welch, J. E., Numerical calculation of time-dependent viscous incompressible flowof fluid with a free surface, Phys. Fluids, 8 2182–2189 (1965)[13] 藤井昭宏ら, 領域分割による並列 AMG アルゴリズム, 情報処理学会論文誌コンピューティ
ングシステム（ACS）, 44, SIG06(ACS1), 9–17 (2003)[14] Roe, P. L., Characteristic-based schemes for the Euler equations, Annu. Rev. Fluid Mech., 18, 337–
365 (1986)[15] Van der Vorst, H. A., Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the
Solution of Nonsymmetric Linear Systems, SIAM J. Sci. Stat. Comput., 13(2), 631–644 (1992)[16] Karypis, G. and Kumar, V., A Fast and Highly Quality Multilevel Scheme for Partitioning Irregular
Graphs, SIAM J. Sci. Comput., 20(1), 359–392, (1999)[17] Matsushita, Y. et al., Efficient communication strategy in parallel computation based on domain
partitioning, J. Chem. Eng. Jpn., 51(1), 79–82 (2018)[18] Raman, V. and Pitsch, H, A consistent LES/filtered-density function formulation for the simulation of
turbulent flames with detailed chemistry, Proc. Combust. Inst., 31(2), 1711–1719 (2007)[19] Pierce, C. and Moin, P., Method for generating equilibrium swirling inflow conditions, AIAA J., 36(7),1325–1327 (1998)[20] Klein, M., A digital filter based generation of inflow data for spatially developing direct numerical or
large eddy simulations, J. Comput. Phys., 186(2), 652–665 (2003)[21] Hossain, M. et al., Modelling of a Bluff-Body Nonpremixed Flame using a Coupled
Radiation/Flamelet Combustion Model, Flow, Turbul. Combust., 67(3), 217–234 (2001)

― 8 ―

SENAC Vol. 52, No. （3 2019. 7）

民間航空機開発における大規模 CFD 解析の適用（その 3）― 9 ―― 10 ―

SENAC Vol. 52, No. （3 2019. 7）

民間航空機開発における大規模 CFD 解析の適用（その 3）― 11 ―

― 12 ―

SENAC Vol. 52, No. （3 2019. 7）

民間航空機開発における大規模 CFD 解析の適用（その 3）― 13 ―

― 14 ―

SENAC Vol. 52, No. （3 2019. 7）

SENAC Vol. 52, No. （3 2019. 7）― 15 ―[お 知 ら せ ]サイバーサイエンスセンター講習会(夏期・秋期開催分)のご案内
サ イバー サイ エンス セン ターで は、利 用者の みな さまに 当セ ンター の計 算資源 を効 率的に 利用し ていた だく ことを 目的 に講習 会を 開催し てお ります。初め てスー パー コンピ ュー タを利 用される 方や現 在活 用され てい る方を 対象 に、幅 広い カリキ ュラ ムを用 意し ていま す。
今 後スー パー コンピ ュー タによ るシ ミュレ ーシ ョンを 計画 してい る研 究者、学生、技術系職 員
等 広い範 囲の 皆様の 参加 をお待 ちし ており ます 。

No.

講 習 会 名開 催 日 時

募集
人数

講 師

内 容

9

Gaussian 入門

8 月 21 日(水)20

岸本・Gaussian の基本的な使い方

13:00-17:00
10

Mathematica 入門

9 月 6 日(金)(理学研究科)20

13:00-16:30
11

はじめての Linux

9 月 9 日(月)はじめてのスパコン 9 月 10 日(火)20
20

13:00-15:30
13

SX-ACE の性能分析・ 9 月 11 日(水)高速化・Mathematica の基本的な使い方(尚絅学院大学)

13:00-16:00
12

横井
佐々木

・Linux システムの基本的な使い方(情報部情報基盤課)

・エディタの使い方

大泉・スーパーコンピュータの紹介と(情報部情報基盤課)

20

13:00-17:00

江川(サイバーサイエンス

利用法入門（見学あり）
・スーパーコンピュータでの性能解析から最適化まで

センター)

14

並列プログラミン 9 月 12 日(木)20

グ入門Ⅰ（OpenMP） 13:00-17:00

小松・並列プログラミングの概要(サイバーサイエンス ・OpenMP による並列プログラミングの
センター)

基礎
・利用法

15

並列プログラミング 9 月 13 日(金)入門Ⅱ（MPI）

13:00-17:00

20

小松(サイバーサイエンスセンター)

・MPI による並列プログラミングの

基礎
・利用法

備考：・申し込みは、ウェブページ https://www.ss.cc.tohoku.ac.jp/guide/kosyu.cgi からお願いします。
・プログラムは予定のものです。若干変更になる場合がありますのでお含みおきください。
・サイバーサイエンスセンターでは、希望に応じてスーパーコンピュータに関する出張講習会の開催を検討いたします。希望される方は共同利用支援係までご連絡ください。
問合せ先：共同利用支援係（022-795-3406,uketuke@cc.tohoku.ac.jp）

― 16 ―

SENAC Vol. 52, No. （3 2019. 7）

講習会受講者の感想

はじめての Linux・初心者にとって、基本的かつ実用的なものであってとても参考になった。
・入門として非常にわかりやすかった。
・コマンドの復習になってよかった。
・資料がよかった。

はじめてのスパコン
・スーパーコンピュータのことはまったく知らなかったので興味深かった。
・実際に手を動かす講習だったので判り易く大変ためになる講習だった。
・スパコンを実際に見ることができて良かった。

並列プログラミング入門Ⅰ（OpenMP）
・全く OpenMP を知らない人にもわかりやすい
内容でした。
・すぐ応用に使えそうだ。
・演習が充実していてとてもよかった。

Gaussian 入門
・講義から実践まで新鮮でした。
・化学にほとんど詳しくない自分にとっても理解でした、とてもありがたい。
・先生の説明は判り易かったです。
・コンピュータに関する知識がうといのでかな
り難しい内容だった。

＜講習会会場＞
東北大学サイバーサイエンスセンター
仙台市青葉区荒巻字青葉 6-3
仙台市地下鉄東西線

青葉山駅出口より徒歩 3 分

並列プログラミング入門Ⅱ（MPI）
・講習会用の資料がすごくわかりやすかったです（プログラム例、図）。
・とても判り易く参考になった。充実していて面白かった。

Mathematica 入門
・基本的なことがよくわかった。丁寧に解説
いただきありがとうございました。
・自分で操作してそれをチェックしてもらえるのでよかった。

SENAC Vol. 52, No. （3 2019. 7）― 17 ―[利用相談室便り]

2019年度の利用相談について

今年度も 5 月よりサイバーサイエンスセンター本館利用相談室で利用相談を行っています。日程
等詳細は次頁をご覧ください。相談内容によってはメーカ等に問い合わせる場合や、時間を要する場合もありますが、利用者の問題解決にむけて努めております。直接面談のほかに、メールや電話での相談も受けておりますのでお気軽にご相談ください。

・プログラムを高速化するにはどうしたらいいの？
・プログラムを並列化してもっと速く計算したい！
・スパコンでプログラムを動かしても速さがPCと変わらないんだけど、どうして？
・研究室のコンピュータではメモリが足りない！
・研究室の電気代高騰で困っている。
・コンピュータの管理は面倒。研究に専念したい。
・サービスしているアプリケーションを研究室から利用するにはどうすればいいの？
このような、スーパーコンピュータ利用に関する疑問や問題をお持ちの方、これから利用してみ
たいとお考えの方、一度相談してみてはいかがでしょうか。また、サイバーサイエンスセンター本館相談室には、各種マニュアル、書籍も揃えています。相談室での閲覧、貸し出し（一部の書籍､マニュアルを除く）も可能ですのでご活用ください。

東北大学サイバーサイエンスセンター本館1階 利用相談室
所在地：
Tel
:

仙台市青葉区荒巻字青葉 6-3 （地下鉄東西線青葉山駅 北 1 出口より徒歩 3 分）
022-795-6153 学内内線 92-6153
不在の場合は 022-795-3406(3 階共同利用支援係)
e-mail : sodan@cc.tohoku.ac.jp
e-mail の場合、曜日、時間帯によっては回答に少しお時間をいただくことがあります。あらかじめご了承ください。

サイバーサイエンスセンター本館（右）、2 号館（左）利用相談室

― 18 ―

SENAC Vol. 52, No. （3 2019. 7）

2019 年度利用相談日程と主な担当分野
曜日・時間

月

14～16 時

火

14～16 時水木

金

15～17 時
要事前予約

14～16 時

14～16 時
要事前予約

要事前予約

ﾃｸﾆｶﾙｱｼｽﾀﾝﾄ（ 所 属 ）主な担当分野・スーパーコンピュータ
・並列コンピュータ
佐々木大輔（情報基盤課共同研究支援係）
・Fortran
・大判プリンタ
・可視化システム・スーパーコンピュータ
齋藤 敦子（情報基盤課共同研究支援係）
・並列コンピュータ
・大判プリンタ
山崎

森谷

小松

山下

馨

（金属材料研究所）

・アプリケーション（Gaussian）

・スーパーコンピュータ
・並列コンピュータ
友映（情報基盤課共同研究支援係）
・Fortran
・大判プリンタ・スーパーコンピュータ
・並列コンピュータ
一彦（サイバーサイエンスセンター）・高速化（ベクトル化、並列化）
・Fortran ・C/C++
・アプリケーション全般
・高速化（ベクトル化、並列化）
毅 （情報基盤課共同利用支援係）
・Fortran
・負担金*上記以外の時間帯に面談・電話での相談を希望の方は、3 階窓口（共同利用支援係）まで
相談内容をお申し出ください。センター内担当者に取り次ぎます。

SENAC Vol. 52, No. （3 2019. 7）― 19 ―［JHPCN］

JHPCN 学際大規模情報基盤共同利用・共同研究拠点
第 11 回シンポジウム報告
江川 隆輔
東北大学サイバーサイエンスセンター スーパーコンピューティング研究部【報告記事】
令和元年 7 月 11，12 日，第 11 回 学際大規模情報基盤共同利用・共同研究拠点(JHPCN)シンポ
ジウムが，東京品川 The Grand Hall において開催されました．今年度のシンポジウムでは昨年度採択された 51 課題に関する講演と今年度採択された課題 58 件のポスター発表に加えて，各構成
拠点から推薦された萌芽課題 30 件のポスター発表があり，活発な議論が展開されました．
学際大規模情報基盤共同利用・共同研究は，東北大学，北海道大学，東京大学，東京工業大学，
名古屋大学，京都大学，大阪大学，九州大学にそれぞれ附置するスーパーコンピュータを持つ 8
つの共同利用の施設を構成拠点とする「学際大規模情報基盤共同利用・共同研究拠点」が, 超大規模数値計算系応用分野，超大規模データ処理系応用分野，超大容量ネットワーク技術分野，およびこれらの技術分野を統合した大規模情報システム関連研究分野，更には分野間に亘る複合分
野の研究者らと取り組む学際的な共同利用・共同研究です．
平成 22 年度から令和元年度の 8 年間において約 350 件を超える課題が学際大規模情報基盤共同
利用・共同研究として採択されており，そのうち当センターとの共同研究課題は約 50 件にとなっております．今年度は，7 件が当センターとの共同研究課題として採択されております．本年度採択された課題のポスターに関しましては，本報告末尾にポスターを再掲しております．是非，高性能計算を用いた多岐にわたる共同研究活動をご覧頂ければと存じます．その他，これまでの採択課題に関する情報は以下の URL で公開されておりますので，併せてご参考にして頂ければ幸いです. (JHPCN URL : https://jhpcn-kyoten.itc.u-tokyo.ac.jp/dg/sympo/11th/)
来年度の学際大規模情報基盤共同利用・共同研究拠点の研究公募は 10 月頃に公開予定です．ご
興味のあるかたは奮ってご応募ください．また，応募に際しまして当センターの計算機科学を専門とする教員との共同研究の可能性を検討したい，手続き方法が分からない等，本応募に関して不明な点があります場合は，お気軽に当センターまでお問い合わせください．

― 20 ―

SENAC Vol. 52, No. （3 2019. 7）【平成 31 年度学際大規模情報基盤共同利用・共同研究拠における当センターとの共同研究課題】

熱中症リスク評価シミュレータの開発と応用
研究代表者

平田晃正 (名古屋工業大学)

日本全土の洪水氾濫被害推定の高精度化研究代表者 風間聡 (東北大学)

大規模津波浸水被害推計シミュレーションのマルチプラットフォーム向け最適化手法の研究
研究代表者 撫佐昭裕 (東北大学)

Investigation of Sound-Flow Interaction of Acoustic Liner using CFD/CAA Hybrid Approach
研究代表者

佐々木大輔 (金沢工業大学)

研究代表者

松尾亜紀子 (慶應義塾大学)

気液二相デトネーションに対する大規模数値解析

GW space-time コードの大規模な有機-金属界面への適用に向けた高効率化
研究代表者 柳澤将 (琉球大学)

HPC と高速通信技術の融合による大規模データの拠点間転送技術開発と実データを用いたシステム
実証試験
研究代表者

シンポジウム会場

村田健史 (情報通信研究機構)

パネル展示会場

JHPCN 学際大規模情報基盤共同利用・共同研究拠点
第 11 回シンポジウム報告― 21 ―

― 22 ―

SENAC Vol. 52, No. （3 2019. 7）

JHPCN 学際大規模情報基盤共同利用・共同研究拠点
第 11 回シンポジウム報告― 23 ―

― 24 ―

SENAC Vol. 52, No. （3 2019. 7）

JHPCN 学際大規模情報基盤共同利用・共同研究拠点
第 11 回シンポジウム報告― 25 ―

― 26 ―

SENAC Vol. 52, No. （3 2019. 7）

JHPCN 学際大規模情報基盤共同利用・共同研究拠点
第 11 回シンポジウム報告― 27 ―

― 28 ―

SENAC Vol. 52, No. （3 2019. 7）［報 告］

平成 31 年度科学技術分野の文部科学大臣表彰において
後藤英昭准教授が科学技術賞（開発部門）を受賞
平成 31 年度科学技術分野の文部科学大臣表彰において、本センター・後藤英昭准教授らの研究グループが科学技術賞（開発部門）を受賞しました。
後藤准教授らの研究は、全国の大学の共通の認証基盤として、ウェブ認証連携の国際標準に基づく学術認証フェデレーション「学認」、国際無線 LAN 相互利用規格「eduroam」による学術無線
LAN ローミング、オープンドメイン認証局による「UPKI 電子証明書発行サービス」を組み合わせ、統合的な認証連携アーキテクチャを設計・構築し実用化したもので、これらの技術が大学の学術コンテンツや学術ネットワーク資源等の安全な共有・共同利用、遠隔講義や単位互換などの大学間学生交流など国内の大学間だけでなく、国際連携や商用サービスなどに活用されていることが高く評価されたものです。
なお、業績名、受賞者は以下のとおりです。
〇

科学技術分野の文部科学大臣表彰 科学技術賞（開発部門）

業績名：大学間連携のための学術認証フェデレーションの開発
受賞者：岡部寿男（京都大学教授／国立情報学研究所客員教授）（筆頭者）
西村健（国立情報学研究所特任研究員）

佐藤周行（東京大学准教授）
後藤英昭（東北大学准教授／国立情報学研究所客員准教授）
曽根原登（津田塾大学教授／国立情報学研究所客員教授・名誉教授）

受賞者の方々。（左から２人目が後藤准教授）

SENAC Vol. 52, No. （3 2019. 7）― 29 ―［報 告］

一般社団法人・電気学会第 107 回通常総会において
江川隆輔准教授が第 75 回電気学術振興賞（進歩賞）を受賞
一般社団法人・電気学会第 107 回通常総会において、本センター・江川隆輔准教授らの研究グループが「第 75 回電気学術振興賞（進歩賞）」を受賞しました。
電気学術振興賞（進歩賞）は、電気に関する学術・技術において新規な概念・理論・材料・デ
バイス・システム・方法等を新たに提案あるいはこれらの提案を実証した者、および電気に関する製品・設備等を新たに完成または改良し，堅調な成果をあげた者に与えられる賞で、江川准教授らの研究グループは、名古屋工業大学・平田晃正教授らとの共同研究を通して、熱中症リスクを複合物理解析と温熱生理モデルを組み合わせた年代別等の人体モデルに基づき発汗量や体温上昇量など定量的なリスクを表現する技術の開発と、「熱中症セルフチェック」として日本気象協会が推進する「熱中症ゼロへ」プロジェクトの Web サイトにて公開を通して、個人を考慮した熱中症リスクの評価を可能とした新規性、一般の方でも活用しやすい実用性が高く評価されたものです。
同准教授らが所属するスーパーコンピューティング研究部（大学院情報科学研究科・高性能計
算論講座）では今後、計算科学を専門とする研究者の皆さまとスーパーコンピュータを活用、その成果を社会に還元できる共同研究を進めていきます。
受賞内容
・題 目 ：「大規模人体複合物理解析に基づく熱中症リスク管理システムの開発と実用化」
・対象者：平田晃正（名古屋工業大学）、江川隆輔（東北大学）、柏 達也（北見工業大学）、
Laakso, Ilkka （アールト大学）、堀江 祐圭（日本気象協会）
なお、賞の詳細については、下記 HP に掲載されております。
電気学会ページ
https://www.iee.jp/blog/award2019/
研究室ページ
https://www.sc.cc.tohoku.ac.jp/1846.html

― 30 ―

SENAC Vol. 52, No. （3 2019. 7）［報 告］

大型計算機センター法制化 50 周年記念シンポジウム報告滝沢
東北大学サイバーサイエンスセンター寛之
スーパーコンピューティング研究部

東北大学 サイバーサイエンスセンターは，以下の 6 センターと共に 2019 年 7 月 10 日(水)に品川 The Grand Hall で「大型計算機センター法制 50 周年記念シンポジウム」を開催しました．
-

北海道大学 情報基盤センター-

東京大学 情報基盤センター-

名古屋大学 情報基盤センター-

京都大学 学術情報メディアセンター-

大阪大学 サイバーメディアセンター-

九州大学 情報基盤研究開発センター

歴史を紐解くと，上述の主催 7 センターが発足したタイミングには多少の違いがあるのですが，全国の学術研究者が利用できるように大型計算機センターの制度が整備された 1969 年から数えて 50 周年となることを記念して本シンポジウムは執り行われました．当日のプログラムを以下に示します．
プログラム：
開会の挨拶

東京大学情報基盤センター 田浦 健次朗 センター長

来賓祝辞

文部科学省研究振興局 橋爪 敦 参事官(情報担当)

招待講演

米国テネシー大学 Jack Dongarra 教授
High Performance Computing and Big Data: Challenges for the Future

招待講演

高度情報科学技術研究機構 小柳 義夫 氏
大学計算センターの歩み

パネルディスカッション
テーマ: 「10 年後の情報基盤センターの在り方を考える」
司

会: 北海道大学情報基盤センター

パネリスト:
北海道大学

杉木 章義 准教授

東北大学

江川 隆輔 准教授

東京大学

星野 哲也 助教

名古屋大学

嶋田 創 准教授

京都大学

深沢 圭一郎 准教授

大阪大学

伊達 進 准教授

九州大学

嶋吉 隆夫 准教授

岩下 武史 教授

大型計算機センター法制化 50 周年記念シンポジウム報告― 31 ―

主催者を代表して東京大学情報基盤センターの田浦センター長より開会の挨拶があり，文部科
学省からは橋爪参事官より法制化 50 周年記念の祝辞をいただきました．
1 人目の招待講演者である Dongarra 教授からは，日本の大学のスパコン整備が中長期的に計画されていることや，スパコンのランキングとして有名な TOP500 List の現状分析， 科学技術計算とビッグデータ処理のエコシステムを構築するうえでの課題などについて興味深い講演がありました．2 人目の招待講演者である小柳氏からは，今となっては知っている人が少なくなってきた
大型計算機センターの歴史を，その前史から含めてご講演いただきました．また，大学の計算機センターはサービスセンターから研究センターに変わっていくべきだという方向性も示されました．
シンポジウムの後半は「10 年後の情報基盤センターの在り方を考える」というテーマでパネル
ディスカッションが行われました．北海道大学の岩下教授が司会を担当し，主催の 7 センターから若手・中堅の教員がパネリストとして登壇しました．パネリストはそれぞれの立場や専門分野からの観点で，10 年後の情報基盤センターの在り方を述べました．その後に会場から多くの発言があり，活発なディスカッションが行われました．「仕事は増える一方だが人は増えない」など，
現在のセンター教職員が抱える悩みを共有しつつ，小柳先生が招待講演で示唆したとおりに研究センターとして発展していくべきだという方向で意見が一致していたように感じました．
以上，本シンポジウムは現在，過去，未来の情報基盤センターの在り方を学び，考えるうえで大変有意義であり，当初の予想以上に盛会となりました．

― 32 ―

SENAC Vol. 52, No. （3 2019. 7）

東京大学情報基盤センター 田浦健次朗 センター長

文部科学省研究振興局 橋爪敦 参事官

米国テネシー大学 Jack Dongarra 教授 1

米国テネシー大学 Jack Dongarra 教授 2（質疑応答）

高度情報科学技術研究機構 小柳義夫氏 1

高度情報科学技術研究機構 小柳義夫氏 2

パネルディスカッション「10 年後の情報基盤センターの在り方を考える」

会場の様子

SENAC Vol. 52, No. （3 2019. 7）― 33 ―[報 告]

EMC Sapporo & APEMC 2019 報 告東北大学情報部情報基盤課
小野敏2019 年 6 月 3 日〜7 日に北海道札幌市の札幌コンベンションセンターにおいて、EMC
Sapporo & APEMC 2019(2019 Joint International Symposium Electromagnetic
Compatibility and Asia-Pacific International Symposium on Electromagnetic
Compatibility, Sapporo)が開催されました。EMC は環境電磁工学に関する国際会議であり、開催された札幌コンベンションセンターでは、企業や大学、研究所などが展示ブースを構え、開発した製品、研究内容等の情報発信をしていました。
サイバーサイエンスセンターでは、本学の高性能計算に関する研究成果の発信を目的に、展示ブースを設置しました。展示ブースでは、本センターのシステム紹介、各研究部の研究内容の紹介、スーパーコンピュータで利用・開発されているアプリケーションの紹介、将来の高性能システムのための研究開発に関する展示・成果発表を行いました。
また、将来のユーザ獲得を目指し無料利用アカウントの配布などを行い、展示ブースに多くの訪問者を迎え、広く当センターの活動をアピールできただけでなく、大変有意義な議論と情報交換を行う事ができました。

― 34 ―

SENAC Vol. 52, No. （3 2019. 7）［紹 介］

サイバーサイエンスセンター展示室の紹介

サイバーサイエンスセンターの 1 階にある展示室をご紹介します。この展示室は 2010 年 3 月に情報処理学会分散コンピュータ博物館として認定されました。
当センターではコンピュータ技術の発展を広く知っていただくため、センターで実際に使用した計算機等を平成 10 年頃から収集を始め、平成 14 年からセンター1 階の展示室に展示しています。
展示品の主なものとしては、コンピュータの黎明期に東北大学と日本電気が共同開発した科学技術
用パラメトロン計算機の資料と、歴代のコンピュータを中心としたその部品や装置、歴代のネット
ワーク関連機器と資料などです。
展示室はサイバーサイエンスセンター本館 1 階にあり、2 室に分かれています。展示室 1 にはコンピュータの黎明期の資料、当センターの歴代のコンピュータや関連資料を展示しています。展示室 2 には、ネットワーク関連の展示と大学関係者から寄贈いただいた Mac 関連の機器を展示しています。
展示室は平日 9:00～17:00 の公開となっております。見学を希望される方は、センター1F カウン
ターの電話で内線 801 へご連絡ください。どなたでも見学していただけます。係員による説明をご希望の方、見学の人数が 10 名以上の場合は事前のご連絡をお願いいたします。

パラメトロン演算ユニット

NEAC シリーズ 2200

計算機室（写真展示）

分散コンピュータ博物館（センター展示室）の紹介― 35 ―

SX-9 の正面と側面

ルータ(Cisco 社製 Catalyst 6500)＜過去の紹介記事一覧＞
分散コンピュータ博物館認定式（SENAC vol.43,No.2）
サイバーサイエンスセンター展示室の紹介（SENAC vol.43,No.3）
展示室便り①（SENAC vol.44,No.3）

SENAC-1

展示室便り②（SENAC vol.44,No.4）

SENAC-1（続）展示室便り③（SENAC vol.45,No.1）

大泉充郎先生展示室便り④（SENAC vol.45,No.2）

NEAC シリーズ 2200

展示室便り⑤（SENAC vol.45,No.3）

ACOS シリーズ

展示室便り⑥（SENAC vol.45,No.4）

スーパーコンピュータ SX シリーズ

展示室便り⑦（SENAC vol.46,No.1）

センター刊行物

展示室便り⑧（SENAC vol.46,No.2）

スカラ並列コンピュータ

展示室便り⑨（SENAC vol.46,No.3）

磁気ディスク装置と磁気記録媒体

展示室便り⑩（SENAC vol.46,No.4）

1960 年代の計算機に関する資料

展示室便り⑪（SENAC vol.47,No.1）

地球シミュレータ

展示室便り⑫（SENAC vol.47,No.2）

古 Mac、ワークステーション、PDA

展示室便り⑬（SENAC vol.47,No.4）

ネットワーク

展示室便り⑭（SENAC vol.48,No.3）

スーパーコンピュータ SX-9

＜サイバーサイエンスセンター展示室 URL＞
https://www.cc.tohoku.ac.jp/museum/index.html＜お問合せ先＞
東北大学サイバーサイエンスセンター内共同利用支援係
022-795-3406
uketuke@cc.tohoku.ac.jp

― 36 ―

SENAC Vol. 52, No. （3 2019. 7）[Web 版大規模科学計算システムニュースより]
大規模科学計算システムニュースに掲載された記事の一部を転載しています。 https://www.ss.cc.tohoku.ac.jp/tayori/

利用負担金額の表示コマンドについて (No.278)

本センター大規模科学計算システムでは、利用者の利用額とプロジェクトごとに集計した負担額、請求情報を表示するためのコマンドとして ukakin, pkakin があります。また、利用者のジャーナル
情報とプロジェクトごとに集計したジャーナル情報を CSV 形式で出力するコマンド ulist, plist があります。これらのコマンドは、並列コンピュータ(front.cc.tohoku.ac.jp)にログインして使用し
ます。
コマンド名

機

能

ukakin

利用者ごとの利用額を各システム、月ごとに表示

pkakin

プロジェクトごとに集計した負担額、請求情報を表示

ulist

利用者ごとのジャーナルを CSV 形式で出力

plist

プロジェクトごとに集計したジャーナルを CSV 形式で出力

いずれも、前日までご利用いただいた金額を表示します。コマンド使用例は大規模科学計算システ
ムウェブページをご覧ください。
負担金の確認https://www.ss.cc.tohoku.ac.jp/utilize/academic.html#負担金の確認(共同利用支援係)

利用者番号（アカウント）管理の徹底について （No.278）
当センターの大規模科学計算システムでは、利用に関する内規に記載された利用資格に該当し、利用を承認された本人のみが、利用者番号（アカウント）を利用できます。利用者番号を他者へ貸与することは厳しく禁止されておりますので、改めてご注意をお願いいたします。
利用者番号を登録しても、利用しなければ負担金は発生しませんので、実際に計算機システムを利用する方全員の利用申請を行ってください。また、不正アクセスの危険性が高まりますので、以下のことを行わないよう秘密鍵の管理の徹底もお願いいたします。
・パスフレーズなしの秘密鍵の使用
・秘密鍵、パスフレーズの使い回し
・秘密鍵のメール添付、USB メモリやホームディレクトリに保存
情報セキュリティ対策の強化のため、皆様のご理解とご協力をお願いいたします。
（共同利用支援係）

SENAC Vol. 52, No. （3 2019. 7）― 37 ―― SENAC 執筆要項 ―１．お寄せいただきたい投稿内容
サイバーサイエンスセンターでは、研究者・技術者・学生等の方々からの原稿を募集しております。
以下の内容で募集しておりますので、皆さまのご投稿をお待ちしております。なお、一般投稿いただ
いた方には、謝礼として負担金の一部を免除いたします。
・一般利用者の方々が関心をもたれる事項に関する論説・センターの計算機を利用して行った研究論文の概要
・プログラミングの実例と解説・センターに対する意見、要望
・利用者相互の情報交換

２．執筆にあたってご注意いただく事項(1)原稿は横書きです。
(2)術語以外は、「常用漢字」を用い、かなは「現代かなづかい」を用いるものとします。
(3)学術あるいは技術に関する原稿の場合、200 字～400 字程度のアブストラクトをつけてください。
(4)参考文献は通し番号を付し末尾に一括記載し、本文中の該当箇所に引用番号を記入ください。
・雑誌：著者,タイトル,雑誌名,巻,号,ページ,発行年
・書籍：著者,書名,ページ,発行所,発行年
３．原稿の提出方法
原稿のファイル形式は Word を標準としますが、PDF での提出も可能です。サイズ*は以下を参照してください。ファイルは電子メールで提出してください。
－Word の場合－
・用紙サイズ：A4
・余白：上＝30mm 下＝25mm 左右＝25mm 綴じ代＝0
・標準の文字数（45 文字 47 行）＜文字サイズ等の目安＞
・表題＝ゴシック体 14pt 中央 ・副題＝明朝体 12pt 中央
・氏名＝明朝体 10.5pt 中央
・所属＝明朝体 10.5pt 中央
・本文＝明朝体 10.5pt
・章・見出し番号＝ゴシック体 11pt～12pt*余白サイズ、文字数、文字サイズは目安とお考えください。
４．その他(1)執筆者には、希望により本誌（10 部以内の希望部数）と本誌 PDF 版を進呈します。
(2)一般投稿を頂いた方には謝礼として、負担金の一部を免除いたします。免除額は概ね 1 ページ
1 万円を目安とします。詳細は共同利用支援係までお問い合わせください。
(3)投稿予定の原稿が 15 ページを超す場合は共同利用支援係まで前もってご連絡ください。
(4)初回の校正は、執筆者が行って、誤植の防止をはかるものとします。
(5)原稿の提出先は次のとおりです。
東北大学サイバーサイエンスセンター内 情報部情報基盤課共同利用支援係
e-mail uketuke@cc.tohoku.ac.jp
TEL022-795-3406

― 38 ―

SENAC Vol. 52, No. （3 2019. 7）

スタッフ便り
令和になって最初の SENAC をお送りします。
昭和から平成になったときには、当然ですが国全体が暗い自粛ムードだったのを今でもよく覚えています。しっくりくる表現が思い浮かびませんが、今回の改元はまるでお正月のような、とにかく前回とはまるで違う雰囲気でした。ゴールデンウィークと合わせて異例の長期休暇だったことも、そんなある種のお祭りムードを盛り上げる要因になったのかもしれません。平成生まれの方々にとっては、改元というのは明るく楽しい雰囲気の出来事として記憶に残ることでしょう。
開始当初から明るい幸先のよい令和、これから先も明るくよい時代にしていきましょう！昭和、平成に引き続き、令和も東北大学サイバーサイエンスセンターをよろしくお願いいたします。
（H.T）

去年の夏は猛暑でしたが、今年の夏は平年並みの暑さらしいです。7 月にはオープンキャンパスがあります。いつもは寒いと感じるくらいのスパコン部屋も去年はちょっと涼しいくらいにしか感じられなかったのは、やはり猛暑のせいだったからでしょうか。果たして今年は本当に平年並みなのか、
熱中症対策もしっかりしながら、オープンキャンパスを大いに盛り上げていければと思います。
暑い夏を乗り切るためにも、日頃から適度な運動は大事だと思い、一日 8,000 歩を目標に歩くことにしてみました。ですが、車通勤、デスクワークとなるとなかなか達成できず、休日のみなんとか達成できている状況です。
とりあえず、天気の良い日には、昼休みに散歩をしながら体を動かしつつ、健康に留意して、今年の夏も乗り切りたいものです。(Y.I)

情報部情報基盤課スタッフ異動のお知らせ
2019.7.1 付け
[転入]
高橋 雄一
総務係主任（国際文化研究科総務係主任から）
佐藤 充
情報セキュリティ係主任（情報推進課事務情報係主任から）
[転出]
加藤 美久
小野﨑伸久

総務係主任（加齢医学研究所総務係主任へ）
情報セキュリティ係主任(情報推進課事務情報係主任へ)

SENAC 編集部会
滝沢寛之 水木敬明 後藤英昭 江川隆輔
伊藤昭彦 吉田貴子 大泉健治 小野 敏
斉藤くみ子

編集・発行

サイバーサイエンスセンター
本館（右）、2 号館(左)印

刷

令和元年 7 月発行
東北大学
サイバーサイエンスセンター
仙台市青葉区荒巻字青葉 6-3
郵便番号 980-8578
株式会社 東誠社

ࢩࢫࢸ୍࣒ぴィ⟬ᶵࢩࢫࢸ࣒ᶵ ✀ࢫ࣮ࣃ࣮ࢥࣥࣆ࣮ࣗࢱ

6;$&(୪ิࢥࣥࣆ࣮ࣗࢱ

/;5H

ࢧ࣮ࣂ࡜࣍ࢫࢺྡ
ࣇࣟࣥࢺ࢚ࣥࢻࢧ࣮ࣂIURQWFFWRKRNXDFMS

66+ ࢔ࢡࢭࢫㄆド㘽⏕ᡂࢧ࣮ࣂ

NH\FFWRKRNXDFMSࢧ࣮ࣅࢫ᫬㛫฼⏝ࢩࢫࢸ࣒ྡ➼฼⏝᫬㛫ᖏ

ࢫ࣮ࣃ࣮ࢥࣥࣆ࣮ࣗࢱ

㐃

⥆

㐠

㌿

୪ิࢥࣥࣆ࣮ࣗࢱ

㐃

⥆

㐠

㌿

ࢧ ࣮ ࣂ

㐃

⥆

㐠

㌿

ྍど໬ᶵჾᐊ

ᖹ᪥ 㹼

㤋ෆ฼⏝ᖹ᪥ 㹼

ࢫ࣮ࣃ࣮ࢥࣥࣆ࣮ࣗࢱ㸦6;$&(㸧ࡢ฼⏝ᙧែ࡜ไ㝈್฼⏝ᙧែ฼⏝ࣀ࣮ࢻᩘ ᐇ⾜᫬㛫ไ㝈࣓ࣔࣜࢧ࢖ࢬไ㝈 T ࢜ࣉࢩࣙࣥ E ࢜ࣉࢩࣙࣥ
ͤ㸦⤒㐣᫬㛫㸧つᐃ್㸸 㐌㛫㏻ᖖࠥ↓ᩱ ᫬㛫

ࠥ ᫬㛫

ࠥ ᫬㛫

ࢹࣂࢵࢢ᭱኱್㸸 ࣨ᭶*%ࣀ࣮ࢻᩘ

฼⏝ࣀ࣮ࢻᩘ

V[*%
*%ࣀ࣮ࢻᩘI
GHEXJ

฼⏝ࣀ࣮ࢻᩘͤ ࣀ࣮ࢻ௨ୖࢆ฼⏝ࡋࡓ୪ิᐇ⾜࡟ࡣ03,ࡢ฼⏝ࡀᚲ⏝୪ิࢥࣥࣆ࣮ࣗࢱ㸦/;5H㸧ࡢ฼⏝ᙧែ࡜ไ㝈್฼⏝ᙧែ฼⏝ࣀ࣮ࢻᩘ
ͤ

㏻ᖖ

ࠥ

࢔ࣉࣜ
ࢣ࣮ࢩࣙࣥᐇ⾜᫬㛫ไ㝈㸦⤒㐣᫬㛫㸧つᐃ್㸸 㐌㛫᭱኱್㸸 ࣨ᭶࡞ࡋ࣓ࣔࣜࢧ࢖ࢬไ㝈
*%ࣀ࣮ࢻᩘ

㸦 ࢥ࢔ࡲ࡛㸧฼⏝ࣀ࣮ࢻᩘ
O[*% ᫬㛫
*%
㸦&38 ᫬㛫ྜィ㸧ͤ ࣀ࣮ࢻ௨ୖࢆ฼⏝ࡋࡓ୪ิᐇ⾜࡟ࡣ03,ࡢ฼⏝ࡀᚲ⏝఍ヰᆺT ࢜ࣉࢩࣙࣥ E ࢜ࣉࢩࣙࣥ

D




Vol.52 No.3 2019-7［共同研究成果］
Flamelet approach に基づくガス燃料を対象とした燃焼シミュレーション…………………………………………………………………………赤尾津翔大
松下 洋介
青木 秀之
Weeratunge Malalasekera
民間航空機開発における大規模 CFD 解析の適用（その３）…………………………………………………………………………前田

1

一郎

8

…………………

15［利用相談室便り］
2019 年度の利用相談について ………………………………………………………………

17［JHPCN シンポジウム］
JHPCN 学際大規模情報基盤共同利用・共同研究拠点第 11 回シンポジウム報告………………………………………………………………………………………

19

［報 告］
平成 31 年度科学技術分野の文部科学大臣表彰において
後藤英昭准教授が科学技術賞（開発部門）を受賞 ………………………………………

28

一般社団法人・電気学会第 107 回通常総会において
江川隆輔准教授が第 75 回電気学術振興賞（進歩賞）を受賞 ……………………………

29

大型計算機センター法制化 50 周年記念シンポジウム報告 ………………………………

30

EMC Sapporo & APEMC 2019 報告 ………………………………………………………

33

［紹 介］
サイバーサイエンスセンター展示室の紹介 ………………………………………………

34［Web 版大規模科学計算システムニュースより］
利用負担金額の表示コマンドについて（No.278）…………………………………………

36

利用者番号（アカウント）管理の徹底について（No.278）………………………………

36

執筆要項 …………………………………………………………………………………………

37

スタッフ便り ……………………………………………………………………………………

38［お知らせ］
サイバーサイエンスセンター講習会（夏期・秋期開催分）のご案内