官公需情報ポータルサイト

公告内容

公告期間： ～()1．競争入札に付する事項仕様書のとおり２．入札書等の提出場所等入札説明書等の交付場所及び入札書等の提出場所並びに問い合わせ先(ダイヤルイン)入札説明書等の交付方法上記２．(1)に記載の交付場所または電子メールにより交付する。

ただし、交付は土曜,日曜,祝日及び年末年始(12月29日～1月3日)を除く平日に行う。

電子メールでの交付希望の場合は、「 公告日,入札件名,当機構担当者名,貴社名,住所,担当者所属,氏名,電話,FAX,E-Mail 」を記載し、上記２．(1)のｱﾄﾞﾚｽに送信。

交付の受付期限は の17:00までとする。

入札説明会の日時及び場所入札書の提出期限※技術審査がある場合には、別途、技術審査資料の提出期限があります。入札説明書をご確認下さい。

R3.3.23〒039-3212(1)(2)下記のとおり一般競争入札に付します。

入札公告(郵便入札)青森県上北郡六ヶ所村大字尾駮字表舘２番地１６６請負令和3年3月3日R3.3.3トカマク炉心プラズマ解析用3次元電磁流体シミュレーションコードの整備・改良作業令和3年10月22日件名内容記令和3年4月22日(木)(3)履行期限(2)(4)管理部契約課0175-71-6541履行場所五嶌 由衣六ヶ所核融合研究所(3)(火)ＴＥＬ ＦＡＸ 0175-71-650110時00分令和3年3月23日E-mail：実施しない(4)nyuusatsu_rokkasho@qst.go.jp国立研究開発法人量子科学技術研究開発機構 六ヶ所核融合研究所国立研究開発法人量子科学技術研究開発機構核融合エネルギー部門 六ヶ所核融合研究所管 理 部 長 鈴木 偉久青森県上北郡六ヶ所村大字尾駮字表舘２番地１６６R03RE-01147(1)開札の日時及び場所３．競争に参加する者に必要な資格当機構から指名停止措置を受けている期間中の者でないこと。

全省庁統一競争入札参加資格を有する者であること。

当機構が別に指定する誓約書に暴力団等に該当しない旨の誓約をできること。

４．入札保証金及び契約保証金 免除５．入札の無効入札参加に必要な資格のない者のした入札入札の条件に違反した者の入札６．契約書等作成の要否 要７．落札者の決定方法８．その他 中に当機構ホームページにおいて掲載する。

以上 公告する。

(5)技術審査に合格し、予定価格の制限の範囲内で、最低価格をもって有効な入札を行った入札者を落札者とする。 (最低価格落札方式)管理研究棟令和3年4月22日(木)開札時の立会いは不要とし、開札結果は別途通知する。開札の結果、落札者がなかった場合には再度の入札書の提出期限及び開札日時について別途通知する。

国立研究開発法人量子科学技術研究開発機構 契約事務取扱細則第10条の規定に該当しない者であること。ただし、未成年者、被保佐人又は被補助人であって、契約締結のために必要な同意を得ている者についてはこの限りでない。

六ヶ所核融合研究所上記問い合わせ先宛てに質問書を提出すること。 なお、質問に対する回答は令和3年3月17日 (水)(2)(1)(2)(3)(4)(1)(１)(2) 落札決定に当っては、入札書に記載した金額に当該金額の10パーセントに相当する額を加算した金額(当該金額に1円未満の端数があるときは、その端数を切り捨てた金額とする)をもって落札価格とするので、入札者は、消費税に係る課税事業者であるか免税事業者であるかを問わず、見積もった金額の110分の100に相当する金額を入札書に記載すること。

(4) 令和3年3月10日 (水) 11:00までに国立研究開発法人量子科学技術研究開発機構 契約事務取扱細則第11条第１項の規定に該当しない者であること。

(1) この入札に参加を希望する者は、入札書の提出時に、当機構が別に指定する暴力団等に該当しない旨の誓約書を提出しなければならない。

(2)本入札に関して質問がある場合には、前項の誓約書を提出せず、又は虚偽の誓約をし、若しくは誓約書に反することとなったときは、当該者の入札を無効とするものとする。

(3) その他、詳細については、入札説明書によるため、必ず上記２．(2)により、入札説明書の交付を受けること。

(5)11時00分

トカマク炉心プラズマ解析用3次元電磁流体シミュレーションコードの整備・改良作業引合仕様書令和3年3月国立研究開発法人量子科学技術研究開発機構核融合エネルギー部門 六ヶ所核融合研究所核融合炉システム研究開発部プラズマ理論シミュレーショングループ1目次1 一般仕様1.1 件名・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 21.2 目的・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 21.3 契約範囲・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 21.4 作業場所・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 21.5 納期・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 31.6 納入場所・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 31.7 検査条件・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 31.8 納入条件・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 31.9 打合せ・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 31.10 貸与品・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 41.11 大型計算機の使用・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 41.12 プログラム作成等・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 41.13 情報セキュリティの確保・・・・・・・・・・・・・・・・・・・・・・・・・ 41.14 産業財産権等・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 41.15 グリーン購入法の推進・・・・・・・・・・・・・・・・・・・・・・・・・・ 51.16 協議・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 52. 技術仕様2.1 業務概要・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 62.2 IFERC-CSCスーパーコンピュータシステムの主要性能・・・・・・・・・・・・・62.3 シミュレーションコードの概要・・・・・・・・・・・・・・・・・・・・・・・ 62.4 要求項目・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 72.5 提出資料の作成・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・ 8添付資料別添 1 コンピュータプログラム作成等業務特約条項別添 2 本契約において遵守すべき「情報セキュリティの確保」に関する事項別添 3 産業財産権特約条項別添4 抵抗性電磁流体方程式及び拡散方程式に対するBR2法の実装方法21 一般仕様1.1 件名トカマク炉心プラズマ解析用3次元電磁流体シミュレーションコードの整備・改良作業1.2 目的本件は、国立研究開発法人量子科学技術研究開発機構(以下「量研」という。)プラズマ理論シミュレーショングループが計算機シミュレーション活動の一環として開発を進めているトカマク型核融合炉の炉心プラズマ挙動を解析するための3次元電磁流体シミュレーションコードの整備を実施する。トカマク型核融合炉では、圧力分布や電流分布、流速分布の制御によって炉心プラズマの電磁流体力学不安性を回避・抑制し、安定な電気出力を確保する必要がある。信頼性の高いプラズマ制御を実現するためにはテアリングモードや抵抗性壁モード、内部キンクモードなどの比較的成長率の低い不安定モードの時間発展を把握し、閉じ込めへの影響を評価するとともに、モードが発生した場合の抑制手段を確立する必要がある。このような研究開発を推進するにあたり、近年、3 次元電磁流体シミュレーションが実機での計測データの解釈や実験シナリオの事前予測、物理機構の研究による科学的発見のための有効なツールと認識されている。こうした背景からプラズマ理論シミュレーショングループではJT-60SA、ITER、核融合原型炉の開発に資するため、トカマク装置の実機形状を考慮した非構造メッシュを採用した3次元電磁流体シミュレーションコードMUSESの開発を進めている。本件は、開発中のMUSESコードについて、電気抵抗、粘性、熱伝導などの散逸効果を考慮した抵抗性電磁流体モデルを取り扱うための整備及び改良作業を実施する。1.3 契約範囲本件では第2章に定める技術仕様に基づき、以下の作業を行うものとする。1) 課題分析2) 1次元抵抗性電磁流体コードの試作及び動作検証3) MUSESコードへの散逸項の導入及び動作検証4) 報告書作成1.4 作業場所本作業は受注者の環境を利用して実施する。受注者は第 1.11 項に定めるとおり、インターネットを介した利用形態によって IFERC-CSC の計算機を使用して作業を実施することができる。31.5 納期令和3年10月22日(金)1.6 納入場所青森県上北郡六ヶ所村大字尾駮字表舘2番地166量研 六ヶ所核融合研究所 核融合炉システム研究開発部プラズマ理論シミュレーショングループ1.7 検査条件第 1.8項に定める提出図書の確認及び第 1.6項に定める納入場所に第 1.8項に定める品目を納入し、第2章に定める技術仕様を満足することの確認をもって検査とする。1.8 納入条件1) 提出図書受注者は、下表に定める図書を量研担当者と合意した期日に遅延なく提出すること。

提出場所は、第 1.6 項に定める納入場所とし、紙媒体の不要な図書は電子メールによる提出を可とする。受注者は1週間程度の確認期間を考慮して事前に提出すること。

量研担当者は受理しない場合には修正を指示し、受領後、1 週間以内に修正等を指示しないときは、承認したものとする。図書名 内容 提出時期 様式実施計画書工程表、体制表を含む作業実施計画を記述する。変更が生じた際は速やかに変更後の実施計画書を提出し、量研の承認を得ること。契約後速やかに電子ファイル打合せ議事録 議事録打合せ後1週間以内電子ファイル中間報告書工程表、作業実績報告書、本仕様の要求事項に対する作業内容のレビューを記述する。契約後3ヶ月電子ファイル作業報告書 第2.5項に定める提出資料 納品時 紙媒体1部2) 納入品第 1.3 項に定める契約範囲において開発されたコード一式と提出図書の電子ファイルを格納したCD-R等メディア媒体41.9 打合せ本契約に関する打合せは量研六ヶ所核融合研究所で実施するものとし、担当者が認めた場合、Web会議に代えることができるものとする。1.10 貸与品本作業の実施に当たり、量研から以下のものを無償貸与する。1) ソフトウェアMUSESコード 1式(コード本体、メッシュジェネレータ、入力データを含む)不連続ガレルキン法による1次元理想電磁流体コード 1式不連続ガレルキン法による拡散方程式ソルバー 1式2) 図書関連報告書類 1式貸与品は、作業完了時に全て返却すること。1.11 大型計算機の使用受注者は、IFERC-CSCの並列計算機システム JFRS-1を無償で使用できるものとする。

インターネットを介した計算機へのアクセス及び計算機利用に必要となる機材は受注者が用意し、計算機の使用に当たっては IFERC-CSC 計算機の利用規則及び量研の情報セキュリティ管理規定を遵守するものとする。1.12 プログラム作成等プログラム作成等については、別添１「コンピュータプログラム作成等業務特約条項」に定められたとおりとする。1.13 情報セキュリティの確保情報セキュリティの確保については、別添２『本契約において遵守すべき「情報セキュリティの確保」に関する事項』のとおりとする。1.14 産業財産権等1) 産業財産権の取扱い本契約に関して発生する知的財産権の取扱いについては、別添３「産業財産権特約条項」に定められたとおりとする。2) 技術情報の開示制限受注者は、本契約を実施することにより得た技術情報を第三者に対して開示しよ5うとするときは、あらかじめ書面により量研の承認を得なければならないものとする。量研が本契約に関し、その目的を達成するため受注者の保有する技術情報を了知する必要が生じた場合は、量研と受注者協議の上、決定するものとする。3) 成果の公開受注者は、本契約に基づく業務の内容及び成果について、発表若しくは公開し、又は特定の第三者に提示しようとするときは、あらかじめ書面により量研の承認を得なければならないものとする。1.15 グリーン購入法の推進1) 本契約において、グリーン購入法(国等による環境物品等の調達等に関する法律)に適用する環境物品(事務用品、OA 機器等)が発生する場合は、これを採用するものとする。2) 本仕様に定める提出図書(納入印刷物)については、グリーン購入法の基本方針に定める「紙類」の基準を満たしたものであること。1.16 協議本仕様書に記載されている事項及び本仕様書に記載のない事項について疑義が生じた場合は、量研と協議の上、その決定に従うものとする。62. 技術仕様2.1 業務概要本件ではトカマク炉心プラズマ解析用 3 次元電磁流体シミュレーションコードの整備・改良作業に関し、以下の業務を実施する。1) 仕様分析(打合せを含む)2) プログラム設計3) プログラミング4) デバッグ5) テスト・ラン6) 処理性能評価7) 作業報告書等の納入資料の作成2.2 IFERC-CSCスーパーコンピュータシステムの主要性能1) 大規模並列演算部(1,370ノード)3.072 TFLOPS/ノード1ノード2CPU: Intel Xeon Gold 6148 ,2.4GHz/20コアメモリ192GB2) フロントシステム部(4ノード)3.072TFLOPS/ノード1ノード 2CPU: Intel Xeon Gold 6148 (2.4GHz/20コア)メモリ768GB3) データ処理サーバ部(2ノード)3.072 TFLOPS/ノード1ノード 2CPU: Intel Xeon Gold 6148 (2.4GHz/20コア)メモリ 1.5TB2.3 シミュレーションコードの概要1) 名称MUSES2) プログラミング言語Fortran903) 概要現行の MUSES コードは 3 次元保存型理想電磁流体方程式を Runge-KuttaDiscontinuous Galerkin 法により離散化しており、ソレノイダル条件に関する拘束条件は課していない。要素形状はくさび形 (wedge) 形状であり、1次の基底関数を7用いた空間 2次精度、体積分及び面積分には 2次精度のガウス求積、時間積分には2次精度TVD Runge-Kutta法を使用する。また、MUSESコードは領域分割法によるMessage Passing Interface (MPI) 並列計算に対応しており、逐次計算と並列計算の切り替えはプリプロセッサにより達成される。コード行数は前処理 5,785 行(内メッシュジェネレータ2,916行)、本体1,663行、後処理710行の合計8,158行である。2.4 要求項目2.4.1 一般事項プログラミングには本仕様において特に指定する以外は Fortran 言語(Fortran90 以降の仕様)を用いること。プログラム開発に当たって必要が生じた場合には LAPACKなどの数値計算ライブラリを用いてもよい。ただし、使用するライブラリはIFERC-CSCの並列計算機システムで利用できるものに限る。受注者は作業に使用したライブラリについて、IFERC-CSC及び一般的なUNIX環境でのライブラリの導入方法や使用方法に関するリファレンスマニュアルを作業報告書に記載すること。また、MUSESコード本体の改良作業はプリプロセッサによる逐次処理・並列処理の切り替えに対応すること。2.4.2 1次元抵抗性電磁流体コードの試作及び動作検証2.4.2.1 1次元抵抗性電磁流体コードの試作受注者は、別添4に示す抵抗性電磁流体方程式に対する1次元コードを量研の指定する方法に従い試作する。量研が貸与する不連続ガレルキン法に基づく1次元理想電磁流体コードに散逸物理として、粘性・熱伝導・電気抵抗を導入する。1次元理想電磁流体コードは0次及び1次の基底関数を採用することで空間精度を可変としており、空間積分は2次のガウス求積、時間積分は2次精度TVD Runge-Kutta法、要素境界における数値流束は移流項にてHLLD法を採用しており、同コードにBR2法 [F.

Bassi et al., "Discontinuous Galerkin solution of Reynolds-averaged Navier–Stokes and k-εturbulence model equation," Comput. Fluids 34, 507-540 (2005)] に基づく散逸項を実装する。受注者は方程式の離散化を含めたプログラム設計を実施すること。BR2法の実装に際し、別添4に概要を示す量研が貸与する拡散方程式のBR2法ソルバーを参考として良いが、受注者がより有利な散逸項の離散化手法を提示できる場合、量研と協議の上離散化手法を変更することができる。2.4.2.2 1次元抵抗性MHDコードの動作検証受注者は、本項で整備した機能の動作確認として、アルフベン波の線形移流問題に8よる数値実験を行う。理想電磁流体との定性的な差異及び格子収束性について論ずる。また、衝撃波管問題による数値実験を行い、数値安定性について論ずる。ただし、2つの問題においてBR2法の安定化パラメータは同一とし、衝撃波管問題については空間1次精度モードのみ検証を行うものとする。得られた結果は作業報告書に記載すること。量研が貸与する1次元電磁流体コードにはアルフベン波の線形移流問題と衝撃波管問題がテスト問題として実装済みであり、受注者は必要に応じて該当箇所を変更し、動作検証に利用することができる。2.4.3 MUSESコードへの散逸項の導入及び動作検証2.4.3.1 MUSESコードへの散逸項の導入前項で動作検証を行った離散化手法に基づき、現行のMUSESコードへ粘性・熱伝導・電気抵抗の導入を行う。MUSESコードはMPIによる並列化及びプリプロセッサによる逐次・並列切り替え機能が実装されている。そのため、BR2法の1段目で計算する物理量の勾配ベクトル若しくは散逸流束ベクトルを通信できるように、既に実装されているMPI通信モジュールを改修する必要がある。2.4.3.2 改修したMUSESコードの動作検証量研が提示する理想電磁流体ベンチマーク問題に対して散逸係数を導入して数値実験を実施する。数値実験結果は可視化の上、作業報告書に記載すること。また、逐次計算と並列計算で計算結果が変わらないことを検証すること。量研が貸与するメッシュジェネレータにはベンチマーク問題は既に実装済みであり、受注者は抵抗性を考慮する改修を行った上で動作検証に利用することができる。2.5 提出資料の作成1) 提出資料の作成作業報告書には、本仕様に基づき、受注者が実施した作業内容及び動作試験、性能測定、可視化結果に加え、プログラムのインストールマニュアル及び使用手引きを含めること。2) 提出資料に関わる電子ファイル作業報告書は紙媒体及びCD-R等のメディア媒体として提出する。開発したコード一式と合わせて、可視化や性能試験、動作試験に使用したコンパイラオプションの設定ファイルやジョブ投入スクリプト類も合わせて納入すること。別添1コンピュータプログラム作成等業務特約条項(目的物)第１条 この契約の目的物は、次の各号の一又は二以上の組み合せに該当するコンピュータプログラムの著作物(データ、データベース、マニュアル及びドキュメンテーションを含む。以下同じ。)及び当該コンピュータプログラムによる計算結果であって、仕様書に定める範囲のものとする。(1) コンピュータプログラム(コンピュータプログラムの設計を含む。)著作物(2) 量研機構が提供するコンピュータプログラムの著作物により得られた計算結果(3) 受注者が所有するコンピュータプログラムの著作物及びこれにより得られた計算結果(権利義務の譲渡等)第２条 受注者は、この契約に基づく権利又は義務を、第三者に譲渡し、若しくは承継させ、又は担保その他の目的に供してはならない。ただし、あらかじめ書面により量研機構の承認を得た場合は、この限りではない。(権利の帰属等)第３条 この業務により作成された目的物(第1条各号に掲げるものをいう。以下同じ。)に係る著作権その他この目的物の使用、収益及び処分(複製、翻訳、翻案、変更、譲渡・貸与及び二次的著作物の利用を含む。)に関する一切の権利は量研機構に帰属するものとする。

ただし、本契約遂行のために使用するプログラム等のうち、本契約締結以前から、受注者が所有するものについては、その著作権は受注者に帰属するものとする。(氏名の表示の制限)第４条 受注者は、第1条に規定する著作物に著作者氏名を表示しないものとする。(第三者の権利の保護)第５条 受注者は、この業務の実施に関し第三者(著作者を含む。)の著作権その他の権利を侵害することのないよう必要な措置を自らの責任において講じなければならない。(技術情報)第６条 量研機構が、この業務の実施に関し、受注者の保有する技術情報を知る必要が生じた場合には、受注者は、この契約の業務に必要な範囲内において当該技術情報を量研機構に無償で 提供しなければならない。２ 量研機構は、受注者からの書面による事前の同意を得た場合を除き前項により知り得た技術情報を第三者に提供しないものとする。(プログラム開発に必要な技術情報)第７条 量研機構は、仕様書に定めるところにより、受注者がこの業務の実施に必要な計算コードその他必要な技術情報を受注者に使用させることがある。(危険負担)第８条 この業務の完了前に目的物について生じた損害その他この業務の実施に関して生じた損害は、受注者の負担とする。ただし、その損害のうち受注者の責めに帰し難い事由により生じたものについては、この限りでない。別添２本契約において遵守すべき「情報セキュリティの確保」に関する事項１．受注者は、契約の履行に関し、情報システム(情報処理及び通信に関わるシステムであって、ハードウェア、ソフトウェア及びネットワーク並びに記録媒体で構成されるものをいう。)を利用する場合には、量研機構の情報及び情報システムを保護するために、情報システムからの情報漏洩、コンピュータウィルスの侵入等の防止その他必要な措置を講じなければならない。２．受注者は、次の各号に掲げる事項を遵守するほか、量研機構の情報セキュリティ確保のために、量研機構が必要な指示を行ったときは、その指示に従わなければならない。1.(1) 受注者は、契約の業務に携わる者(以下「業務担当者」という。)を特定し、それ以外の者に作業をさせてはならない。1.(2) 受注者は、契約に関して知り得た情報(量研機構に引き渡すべきコンピュータプログラム著作物及び計算結果を含む。以下同じ。)を取り扱う情報システムについて、業務担当者以外が当該情報にアクセス可能とならないよう適切にアクセス制限を行うこと。1.(3) 受注者は、契約に関して知り得た情報を取り扱う情報システムについて、ウィルス対策ツール及びファイアウォール機能の導入、セキュリティパッチの適用等適切な情報セキュリティ対策を実施すること。1.(4) 受注者は、P2P ファイル交換ソフトウェア(Winny、WinMX、KaZaa、Share 等)及び SoftEtherを導入した情報システムにおいて、契約に関して知り得た情報を取り扱ってはならない。1.(5) 受注者は、量研機構の承諾のない限り、契約に関して知り得た情報を量研機構又は受注者の情報システム以外の情報システム(業務担当者が所有するパソコン等)において取り扱ってはならない。1.(6) 受注者は、委任をし又は下請負をさせた場合は、当該委任又は下請負を受けた者の契約に関する行為について、量研機構に対し全ての責任を負うとともに、当該委任又は下請負を受けた者に対して、情報セキュリティの確保について必要な措置を講ずるように努めなければならない。1.(7) 受注者は、量研機構が求めた場合には、情報セキュリティ対策の実施状況についての監査を受け入れ、これに協力すること。1.(8) 受注者は、量研機構の提供した情報並びに受注者及び委任又は下請負を受けた者が契約業務のために収集した情報について、災害、紛失、破壊、改ざん、き損、漏えい、コンピュータウィルスによる被害、不正な利用、不正アクセスその他の事故が発生、又は生ずるおそれのあることを知った場合は、ただちに量研機構に報告し、量研機構の指示に従うものとする。契約の終了後においても、同様とする。なお、量研機構の入札に参加する場合、または量研機構からの見積依頼を受ける場合にも、上記事項を遵守していただきます。以 上別添３産業財産権特約条項(受注者が単独で行った発明等の産業財産権の帰属)第１条 受注者は、本契約に関して、受注者が単独でなした発明又は考案(以下「発明等」という。)に対する特許権、実用新案権又は意匠権(以下「特許権等」という。)を取得する場合は、単独で出願できるものとする。ただし、出願するときはあらかじめ出願に際して提出すべき書類の写しを添えて量研機構に通知するものとする。(受注者が単独で行った発明等の特許権等の譲渡等)第２条 受注者は、受注者が前条の特許権等を量研機構以外の第三者に譲渡又は実施許諾する場合には、本取扱いの各条項の規定の適用に支障を与えないよう当該第三者と約定しなければならない。(受注者が単独で行った発明等の特許権等の実施許諾)第３条 量研機構は、第１条の発明等に対する特許権等を無償で自ら試験又は研究のために実施することができる。量研機構が量研機構のために受注者以外の第三者に製作させ、又は業務を代行する第三者に再実施権を許諾する場合は、受注者の承諾を得た上で許諾するものとし、その実施条件等は量研機構、受注者協議の上、決定する。(量研機構及び受注者が共同で行った発明等の特許権等の帰属及び管理)第４条 量研機構及び受注者は、本契約に関して共同でなした発明等に対する特許権等を取得する場合は、共同出願契約を締結し、共同で出願するものとし、出願のための費用は、量研機構、 受注者の持分に比例して負担するものとする。(量研機構及び受注者が共同で行った発明等の特許権等の実施)第５条 量研機構は、共同で行った発明等を試験又は研究以外の目的に実施しないものとする。ただし、量研機構は量研機構のために受注者以外の第三者に製作させ、又は業務を代行する第三者に実施許諾する場合は、無償にて当該第三者に実施許諾することができるものとする。２ 受注者が前項の発明等について自ら商業的実施をするときは、量研機構が自ら商業的実施をしないことに鑑み、受注者の商業的実施の計画を勘案し、事前に実施料等について量研機構、受注者協議の上、別途実施契約を締結するものとする。(秘密の保持)第６条 量研機構及び受注者は、第１条及び第４条の発明等の内容を出願により内容が公開される日まで他に漏えいしてはならない。

ただし、あらかじめ書面により出願を行った者の了解を得た場合はこの限りではない。(委任・下請負)第７条 受注者は、本契約の全部又は一部を第三者に委任し、又は請け負わせた場合においては、その第三者に対して、本取扱いの各条項の規定を準用するものとし、受注者はこのために必要な措置を講じなければならない。２ 受注者は、前項の当該第三者が本取扱いに定める事項に違反した場合には、量研機構に対し全ての責任を負うものとする。(協議)第８条 第１条及び第４条の場合において、単独若しくは共同の区別又は共同の範囲等について疑義が生じたときは、量研機構、受注者協議して定めるものとする。(有効期間)第９条 本取扱いの有効期限は、契約締結の日から当該特許権等の消滅する日までとする。以上別添41 抵抗性電磁流体方程式抵抗性電磁流体方程式は次に示す保存形式の方程式により表現される．∂Q∂t+∂Ea∂x+∂Fa∂y+∂Ga∂z=∂Ed∂x+∂Fd∂y+∂Gd∂z, Q=ρρuρvρwBxByBzε, (1)Ea =ρuρu2 +pT −B2xρuv−BxByρuw−BxBz0uBy −vBxuBz −wBx(ε+pT)u−(u · B)Bx, Fa =ρvρuv−BxByρv2 +pT −B2yρvw−ByBzvBx −uBy0vBz −wBy(ε+pT)v −(u · B)By, Ga =ρwρuw−BxBzρvw−ByBzρw2 +pT −B2zwBx −uBzwBy −vBz0(ε+pT)w−(u · B)Bz, (2)Ed =0τxxτxyτxz0ηJz−ηJyαx +η(ByJz −BzJy), Fd =0τyxτyyτyz−ηJz0ηJxαy +η(BzJx −BxJz), Gd =0τzxτzyτzzηJy−ηJx0αz +η(BxJy −ByJx), (3)αx = uτxx +vτxy +wτxz −qx, (4)αy = uτyx +vτyy +wτyz −qy, (5)αz = uτzx +vτzy +wτzz −qz. (6)ここで，ρは質量密度，u = [u v w]T は流速，B= [Bx By Bz]T は磁束密度，εは全エネルギー密度，pT は熱的圧力と磁気圧の和，η は抵抗率である．また，粘性応力τ，熱流束q，電流Jは，τij = µv (∂ui∂xj+∂uj∂xi )+(βv −23µv)(∇ · u), (7)q =−κ∇T, (8)J =∇×B, (9)である．ここで，µv は剪断粘性係数，βv は体積粘性率，κは熱伝導係数，T は温度である．また，pT 及びT は理想気体に対する状態方程式，pT = p+12B · B, (10)p =ρRT = (γ −1)(ε−12ρu · u−12B · B), (11)により定義される．ただし，γ は比熱比，Rは気体定数である．現行のコードには左辺の移流項のみが実装されており，本案件では右辺の散逸項の実装を行う．12 拡散方程式に対するBR2法の実装方法例として，次に示す放物型方程式を考える．∂Q∂t=∂2Q∂x2. (12)ここで，Qは保存変数である．この方程式を放物型及び楕円型方程式に対する標準解法の1つであるBR2法により離散化するために，次のような保存形式に分解する．∂Q∂t+∂F∂x= 0, (13)F =−∂Q∂x. (14)ここで，F は流束である．BR2法はまず流束F を求め，次にQの時間発展を計算するという2段階の手続きにより構成される．流束F を得るためには式(14)を次のように離散化する．∫ xi+1xiψFdx =−[ψQ∗] +∫ xi+1xiψ′Qdx, (15)Q∗ =Q++Q−2. (16)ここで，ψは試験関数(不連続ガレルキン法においては基底関数と同義)，Q∗ は保存変数Qの数値流束，上付き添字の+と−はそれぞれ隣接セルと自セルを意味する．したがって，数値流束Q∗ は中心差分的に定義される．また，時間発展の式(13)は次のように離散化される．ddt ∫ xi+1xiψQdx =−[ψF∗]xi+1xi+∫ xi+1xiψ′Fdx, (17)F∗ =−(∂Q∂x)++(∂Q∂x)−2+ηδ++δ−2. (18)ここで，数値流束F∗ に含まれるδ はリフティング演算子であり，次式で定義される．δ =F +∂Q∂x. (19)ただし，式(14)よりリフティング演算子は本来δ = 0となるべき物理量であるが，離散化誤差によって有限の値を持ちうることに注意を要する．また，η は安定化パラメータであり，一般にセル境界面の数以上であれば安定に計算することが可能である．2