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

公告内容

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

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

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

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

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

R4.11.10履行期限(2)(4)(3)下記のとおり一般競争入札に付します。

入札公告(郵便入札)請負 R4.10.21管理部経理・契約課量子エネルギー部門 六ヶ所研究所〒039-3212(1)(2)令和4年10月21日量子エネルギー部門 六ヶ所研究所青森県上北郡六ヶ所村大字尾駮字表舘２番地１６６近似ニュートン法による陰解法を実装した１次元抵抗性磁気流体コードの性能分析及び改良作業令和5年2月22日神田 健志件名内容記(3)(木)ＴＥＬ ＦＡＸ 0175-71-650112時00分令和4年11月10日E-mail：令和4年12月5日(月)実施しない国立研究開発法人量子科学技術研究開発機構(4)nyuusatsu_rokkasho@qst.go.jp国立研究開発法人量子科学技術研究開発機構 六ヶ所研究所管 理 部 長 前田 勝青森県上北郡六ヶ所村大字尾駮字表舘２番地１６６R04RE-15726(1)0175-71-6538履行場所開札の日時及び場所３．競争に参加する者に必要な資格当機構から指名停止措置を受けている期間中の者でないこと。

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

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

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

以上 公告する。

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

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

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

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

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

前項の誓約書を提出せず、又は虚偽の誓約をし、若しくは誓約書に反することとなったときは、当該者の入札を無効とするものとする。

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

(2) 本件以外にも、当機構ホームページの調達情報において、今後の「調達予定情報」を掲載しておりますのでご確認下さい。 (URL : https://www.qst.go.jp/site/procurement/(5)13時30分

 ㏆ఝࢽ࣮ࣗࢺࣥἲ࡟ࡼࡿ㝜ゎἲࢆᐇ⿦ࡋࡓ1ḟඖ᢬ᢠᛶ☢Ẽὶయࢥ࣮ࢻࡢᛶ⬟ศᯒཬࡧᨵⰋసᴗᘬྜ௙ᵝ᭩௧࿴4ᖺ10᭶ᅜ❧◊✲㛤Ⓨἲே㔞Ꮚ⛉Ꮫᢏ⾡◊✲㛤Ⓨᶵᵓ㔞Ꮚ࢚ࢿࣝࢠ࣮㒊㛛 භࣨᡤ◊✲ᡤ᰾⼥ྜ⅔ࢩࢫࢸ࣒◊✲㛤Ⓨ㒊ࣉࣛࢬ࣐⌮ㄽࢩ࣑࣮ࣗࣞࢩࣙࣥࢢ࣮ࣝࣉ1 ୍⯡௙ᵝ1.1 ௳ྡ㏆ఝࢽ࣮ࣗࢺࣥἲ࡟ࡼࡿ㝜ゎἲࢆᐇ⿦ࡋࡓ 1ḟඖ᢬ᢠᛶ☢Ẽὶయࢥ࣮ࢻࡢᛶ⬟ศᯒཬࡧᨵⰋసᴗ1.2 ┠ⓗᅜ❧◊✲㛤Ⓨἲே㔞Ꮚ⛉Ꮫᢏ⾡◊✲㛤Ⓨᶵᵓ㸦௨ୗࠕ㔞◊ࠖ࡜࠸࠺ࠋ㸧ࣉࣛࢬ࣐⌮ㄽࢩ࣑࣮ࣗࣞࢩࣙࣥࢢ࣮ࣝࣉ࡛ࡣࠊ୪ิィ⟬ࢆࡣࡌࡵ࡜ࡋࡓ㏆ᖺࡢィ⟬⛉Ꮫࡢⴭࡋ࠸Ⓨᒎࢆ⫼ᬒ࡟ࠊࣉࣛࢬ࣐ࡢ᣺ࡿ⯙࠸ࢆつᐃࡍࡿᇶ♏᪉⛬ᘧ࡟࡛ࡁࡿࡔࡅᛅᐇ࡟ᚑ࠺ࡇ࡜࡟ࡼࡾࠊࣉࣛࢬ࣐୰࡛ᒎ㛤ࡋ࡚࠸ࡿ⌧㇟ࢆィ⟬ᶵୖ࡟෌⌧ࡋࠊ᰾⼥ྜࣉࣛࢬ࣐ࡢ」㞧⌧㇟ࡢ≀⌮ⓗゎ᫂ࢆᅗࡿᩘ್ࢺ࣐࢝ࢡィ⏬㸦௨ୗࠕNEXTィ⏬ࠖ࡜࠸࠺ࠋ㸧ࢆ᥎㐍ࡋ࡚࠸ࡿࠋᮏ௳ࡣࠊNEXTィ⏬࡟࠾࠸࡚◊✲ࢆ㐍ࡵ࡚࠸ࡿ☢Ẽὶయࢩ࣑࣮ࣗࣞࢩࣙࣥᡭἲ㛤Ⓨࡢ୍⎔࡜ࡋ࡚ࠊ㏆ఝࢽ࣮ࣗࢺࣥἲࢆ⏝࠸࡚㝜ゎἲࢆᐇ⿦ࡋࡓ1ḟඖ᢬ᢠᛶ☢Ẽὶయࢥ࣮ࢻ࡟ᑐࡋࠊ࢔ࣝࢦࣜࢬ࣒ࡢᨵⰋసᴗࢆ⾜࠺ࠋᨵⰋ๓ࡢࢥ࣮ࢻ࡜ィ⟬㔞ࡸᩘ್Ᏻᐃᛶࡢẚ㍑ࢆ⾜࠺ࡇ࡜࡛ࠊࣉࣟࢢ࣒ࣛࡢၥ㢟ᢳฟ࡜ᚲせ࡞ಟṇసᴗࢆᐇ᪋ࡍࡿࠋ1.3 ዎ⣙⠊ᅖᮏ௳࡛ࡣ➨2❶࡟ᐃࡵࡿᢏ⾡௙ᵝ࡟ᇶ࡙ࡁࠊ௨ୗࡢసᴗࢆ⾜࠺ࡶࡢ࡜ࡍࡿࠋ1) 㠀⥺ᙧ࣑ࣜࢵࢱ࣮ࢆྵࡵ࡚ᒁᡤ⥺ᙧ໬ࡋࡓ㞳ᩓ᪉⛬ᘧࡢᑟฟ2) ᕪศ㏆ఝࡋࡓࣖࢥࣅ⾜ิ࡟ᇶ࡙ࡃࢽ࣮ࣗࢺࣥἲࡢᐇ⿦సᴗ3) ᩘ್ᐇ㦂4) ሗ࿌᭩సᡂ1.4 సᴗሙᡤᮏసᴗࡣཷὀ⪅ࡢ⎔ቃࢆ฼⏝ࡋ࡚ᐇ᪋ࡍࡿࠋཷὀ⪅ࡣ➨ 1.11 㡯࡟ᐃࡵࡿ࡜࠾ࡾࠊ࢖ࣥࢱ࣮ࢿࢵࢺࢆ௓ࡋࡓ฼⏝ᙧែ࡟ࡼࡗ࡚ IFERC-CSC ࡢィ⟬ᶵࢆ౑⏝ࡋ࡚సᴗࢆᐇ᪋ࡍࡿࡇ࡜ࡀ࡛ࡁࡿࠋ1.5 ⣡ᮇ௧࿴5ᖺ2᭶22᪥㸦Ỉ㸧1.6 ⣡ධሙᡤ㟷᳃┴ୖ໭㒆භࣨᡤᮧ኱Ꮠᑿ㥨Ꮠ⾲⯓2␒ᆅ166㔞◊ භࣨᡤ◊✲ᡤ ᰾⼥ྜ⅔ࢩࢫࢸ࣒◊✲㛤Ⓨ㒊ࣉࣛࢬ࣐⌮ㄽࢩ࣑࣮ࣗࣞࢩࣙࣥࢢ࣮ࣝࣉ1.7 ᳨ᰝ᮲௳➨ 1.8㡯࡟ᐃࡵࡿᥦฟᅗ᭩ࡢ☜ㄆཬࡧ➨ 1.6㡯࡟ᐃࡵࡿ⣡ධሙᡤ࡟➨ 1.8㡯࡟ᐃࡵࡿရ┠ࢆ⣡ධࡋࠊ➨2❶࡟ᐃࡵࡿᢏ⾡௙ᵝࢆ‶㊊ࡍࡿࡇ࡜ࡢ☜ㄆࢆࡶࡗ᳨࡚ᰝ࡜ࡍࡿࠋ1.8 ⣡ධ᮲௳1) ᥦฟᅗ᭩ཷὀ⪅ࡣࠊୗ⾲࡟ᐃࡵࡿᅗ᭩ࢆ㔞◊ᢸᙜ⪅࡜ྜពࡋࡓᮇ᪥࡟㐜ᘏ࡞ࡃᥦฟࡍࡿࡇ࡜ࠋᥦฟሙᡤࡣࠊ➨ 1.6 㡯࡟ᐃࡵࡿ⣡ධሙᡤ࡜ࡋࠊ⣬፹యࡢ୙せ࡞ᅗ᭩ࡣ㟁Ꮚ࣓࣮ࣝ࡟ࡼࡿᥦฟࢆྍ࡜ࡍࡿࠋ㆟஦㘓ࡣཷ⌮ࡋ࡞࠸ሙྜ࡟ࡣಟṇࢆᣦ♧ࡋࠊཷ㡿ᚋࠊ㔞◊ᢸᙜ⪅ࡀ1㐌㛫௨ෆ࡟ಟṇ➼ࢆᣦ♧ࡋ࡞࠸࡜ࡁࡣࠊᢎㄆࡋࡓࡶࡢ࡜ࡍࡿࠋసᴗሗ࿌᭩ࡣ1㐌㛫⛬ᗘࡢ☜ㄆᮇ㛫ࢆ⪃៖ࡋ࡚஦๓࡟ᥦฟࡍࡿࡇ࡜ࠋᅗ᭩ྡ ෆᐜ ᥦฟ᫬ᮇ ᵝᘧᐇ᪋ィ⏬᭩ ᕤ⛬⾲ࠊయไ⾲ࢆྵࡴసᴗᐇ᪋ィ⏬ࢆグ㏙ࡍࡿࠋኚ᭦ࡀ⏕ࡌࡓ㝿ࡣ㏿ࡸ࠿࡟ኚ᭦ᚋࡢᐇ᪋ィ⏬᭩ࢆᥦฟࡋࠊ㔞◊ࡢᢎㄆࢆᚓࡿࡇ࡜ࠋዎ⣙ᚋ㏿ࡸ࠿࡟㟁Ꮚࣇ࢓࢖ࣝᡴྜࡏ㆟஦㘓 ㆟஦㘓 ᡴྜࡏᚋ1㐌㛫௨ෆ㟁Ꮚࣇ࢓࢖ࣝసᴗሗ࿌᭩ ➨2.5㡯࡟ᐃࡵࡿᥦฟ㈨ᩱ ⣡ရ᫬ ⣬፹య1㒊2) ⣡ධရ➨ 1.3 㡯࡟ᐃࡵࡿዎ⣙⠊ᅖ࡟࠾࠸࡚㛤Ⓨࡉࢀࡓࢥ࣮ࢻ୍ᘧ࡜ᥦฟᅗ᭩ࡢ㟁Ꮚࣇ࢓࢖ࣝࢆ᱁⣡ࡋࡓCD-R➼࣓ࢹ࢕࢔፹య1.9 ᡴྜࡏᮏዎ⣙࡟㛵ࡍࡿᡴྜࡏࡣ㔞◊භࣨᡤ◊✲ᡤ࡛ᐇ᪋ࡍࡿࡶࡢ࡜ࡋࠊᢸᙜ⪅ࡀㄆࡵࡓሙྜࠊWeb఍㆟࡟௦࠼ࡿࡇ࡜ࡀ࡛ࡁࡿࡶࡢ࡜ࡍࡿࠋ1.10 ㈚୚ရᮏసᴗࡢᐇ᪋࡟ᙜࡓࡾࠊ㔞◊࠿ࡽ௨ୗࡢࡶࡢࢆ↓ൾ㈚୚ࡍࡿࠋ1) ࢯࣇࢺ࢙࢘࢔1ḟඖ᢬ᢠᛶ☢Ẽὶయࢥ࣮ࢻ 1ᘧFortran࡟ࡼࡿࣈࣟࢵࢡ୕㔜ᑐゅ⾜ิࢯࣝࣂ࣮1ᘧ2) ᅗ᭩㛵㐃ሗ࿌᭩㢮 1ᘧ㈚୚ရࡣࠊసᴗ᏶஢᫬࡟඲࡚㏉༷ࡍࡿࡇ࡜ࠋ1.11 ኱ᆺィ⟬ᶵࡢ౑⏝ཷὀ⪅ࡣࠊIFERC-CSCࡢ୪ิィ⟬ᶵࢩࢫࢸ࣒ JFRS-1ࢆ↓ൾ࡛౑⏝࡛ࡁࡿࡶࡢ࡜ࡍࡿࠋ࢖ࣥࢱ࣮ࢿࢵࢺࢆ௓ࡋࡓィ⟬ᶵ࡬ࡢ࢔ࢡࢭࢫཬࡧィ⟬ᶵ฼⏝࡟ᚲせ࡜࡞ࡿᶵᮦࡣཷὀ⪅ࡀ⏝ពࡋࠊィ⟬ᶵࡢ౑⏝࡟ᙜࡓࡗ࡚ࡣ IFERC-CSC ィ⟬ᶵࡢ฼⏝つ๎ཬࡧ㔞◊ࡢ᝟ሗࢭ࢟ࣗࣜࢸ࢕⟶⌮つᐃࢆ㑂Ᏺࡍࡿࡶࡢ࡜ࡍࡿࠋ1.12 ࣉࣟࢢ࣒ࣛసᡂ➼ࣉࣟࢢ࣒ࣛసᡂ➼࡟ࡘ࠸࡚ࡣࠊูῧ㸯ࠕࢥࣥࣆ࣮ࣗࢱࣉࣟࢢ࣒ࣛసᡂ➼ᴗົ≉⣙᮲㡯ࠖ࡟ᐃࡵࡽࢀࡓ࡜࠾ࡾ࡜ࡍࡿࠋ1.13 ᝟ሗࢭ࢟ࣗࣜࢸ࢕ࡢ☜ಖ᝟ሗࢭ࢟ࣗࣜࢸ࢕ࡢ☜ಖ࡟ࡘ࠸࡚ࡣࠊูῧ㸰ࠗᮏዎ⣙࡟࠾࠸࡚㑂Ᏺࡍ࡭ࡁࠕ᝟ሗࢭ࢟ࣗࣜࢸ࢕ࡢ☜ಖࠖ࡟㛵ࡍࡿ஦㡯࠘ࡢ࡜࠾ࡾ࡜ࡍࡿࠋ1.14 ⏘ᴗ㈈⏘ᶒ➼1) ⏘ᴗ㈈⏘ᶒࡢྲྀᢅ࠸ᮏዎ⣙࡟㛵ࡋ࡚Ⓨ⏕ࡍࡿ▱ⓗ㈈⏘ᶒࡢྲྀᢅ࠸࡟ࡘ࠸࡚ࡣࠊูῧ㸱ࠕ⏘ᴗ㈈⏘ᶒ≉⣙᮲㡯ࠖ࡟ᐃࡵࡽࢀࡓ࡜࠾ࡾ࡜ࡍࡿࠋ2) ᢏ⾡᝟ሗࡢ㛤♧ไ㝈ཷὀ⪅ࡣࠊᮏዎ⣙ࢆᐇ᪋ࡍࡿࡇ࡜࡟ࡼࡾᚓࡓᢏ⾡᝟ሗࢆ➨୕⪅࡟ᑐࡋ࡚㛤♧ࡋࡼ࠺࡜ࡍࡿ࡜ࡁࡣࠊ࠶ࡽ࠿ࡌࡵ᭩㠃࡟ࡼࡾ㔞◊ࡢᢎㄆࢆᚓ࡞ࡅࢀࡤ࡞ࡽ࡞࠸ࡶࡢ࡜ࡍࡿࠋ㔞◊ࡀᮏዎ⣙࡟㛵ࡋࠊࡑࡢ┠ⓗࢆ㐩ᡂࡍࡿࡓࡵཷὀ⪅ࡢಖ᭷ࡍࡿᢏ⾡᝟ሗࢆ஢▱ࡍࡿᚲせࡀ⏕ࡌࡓሙྜࡣࠊ㔞◊࡜ཷὀ⪅࡟ࡼࡿ༠㆟ࡢୖࠊỴᐃࡍࡿࡶࡢ࡜ࡍࡿࠋ3) ᡂᯝࡢබ㛤ཷὀ⪅ࡣࠊᮏዎ⣙࡟ᇶ࡙ࡃᴗົࡢෆᐜཬࡧᡂᯝ࡟ࡘ࠸࡚ࠊⓎ⾲ⱝࡋࡃࡣබ㛤ࡋࠊཪࡣ≉ᐃࡢ➨୕⪅࡟ᥦ♧ࡋࡼ࠺࡜ࡍࡿ࡜ࡁࡣࠊ࠶ࡽ࠿ࡌࡵ᭩㠃࡟ࡼࡾ㔞◊ࡢᢎㄆࢆᚓ࡞ࡅࢀࡤ࡞ࡽ࡞࠸ࡶࡢ࡜ࡍࡿࠋ1.15 ࢢ࣮ࣜࣥ㉎ධἲࡢ᥎㐍1) ᮏዎ⣙࡟࠾࠸࡚ࠊࢢ࣮ࣜࣥ㉎ධἲ(ᅜ➼࡟ࡼࡿ⎔ቃ≀ရ➼ࡢㄪ㐩➼࡟㛵ࡍࡿἲᚊ)࡟㐺⏝ࡍࡿ⎔ቃ≀ရ(஦ົ⏝ရࠊOA ᶵჾ➼)ࡀⓎ⏕ࡍࡿሙྜࡣࠊࡇࢀࢆ᥇⏝ࡍࡿࡶࡢ࡜ࡍࡿࠋ2) ᮏ௙ᵝ࡟ᐃࡵࡿᥦฟᅗ᭩(⣡ධ༳ๅ≀)࡟ࡘ࠸࡚ࡣࠊࢢ࣮ࣜࣥ㉎ධἲࡢᇶᮏ᪉㔪࡟ᐃࡵࡿࠕ⣬㢮ࠖࡢᇶ‽ࢆ‶ࡓࡋࡓࡶࡢ࡛࠶ࡿࡇ࡜ࠋ1.16 ༠㆟ᮏ௙ᵝ᭩࡟グ㍕ࡉࢀ࡚࠸ࡿ஦㡯ཬࡧᮏ௙ᵝ᭩࡟グ㍕ࡢ࡞࠸஦㡯࡟ࡘ࠸࡚␲⩏ࡀ⏕ࡌࡓሙྜࡣࠊ㔞◊࡜༠㆟ࡢୖࠊࡑࡢỴᐃ࡟ᚑ࠺ࡶࡢ࡜ࡍࡿࠋ2.

ᢏ⾡௙ᵝ2.1 ᴗົᴫせᮏ௳࡛ࡣ 1 ḟඖ᢬ᢠᛶ☢Ẽὶయࢥ࣮ࢻ࡬ࡢࢽ࣮ࣗࢺࣥἲ࡟ࡼࡿ㝜ゎἲࡢᐇ⿦ཬࡧ࣋ࣥࢳ࣐࣮ࢡసᴗ࡟㛵ࡋࠊ௨ୗࡢᴗົࢆᐇ᪋ࡍࡿࠋ1) ௙ᵝศᯒ㸦ᡴྜࡏࢆྵࡴ㸧2) ࣉࣟࢢ࣒ࣛタィ3) ࣉࣟࢢ࣑ࣛࣥࢢ4) ࢹࣂࢵࢢ5) ࢸࢫࢺ࣭ࣛࣥ6) సᴗሗ࿌᭩➼ࡢ⣡ධ㈨ᩱࡢసᡂ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ḟඖ᢬ᢠᛶ☢Ẽὶయࢥ࣮ࢻࡢᴫせ1) ྡ⛠1DRMHD2) ࣉࣟࢢ࣑ࣛࣥࢢゝㄒFortran903) ᴫせ⌧⾜ࡢ 1 ḟඖ᢬ᢠᛶ☢Ẽὶయࢥ࣮ࢻࡣ 1 ḟඖಖᏑᆺ᢬ᢠᛶ☢Ẽὶయ᪉⛬ᘧࢆRunge-Kutta Discontinuous Galerkinἲ࡟ࡼࡾ㞳ᩓ໬ࡋ࡚࠸ࡿࠋせ⣲ᙧ≧ࡣ➼㛫㝸᱁Ꮚ࡛࠶ࡾࠊ1ḟࡢᇶᗏ㛵ᩘࢆ⏝࠸ࡓ✵㛫 2ḟ⢭ᗘࠊ✵㛫✚ศ࡟ࡣ 2ḟ⢭ᗘࡢ࢞࢘ࢫồ✚ࠊ᫬㛫✚ศ࡟ࡣ1ẁ2ḟ⢭ᗘ㝜ⓗࣝࣥࢤ࣭ࢡࢵࢱἲࢆ౑⏝ࡍࡿࠋ⛣ὶ㡯ࡢᩘ್ὶ᮰ࡣ HLLD ἲ [T. Miyoshi & K. Kusano, Journal of Computational Physics(2005)]ࠊᣑᩓἲࡢᩘ್ὶ᮰ࡣBR2ἲ [F. Bassi +, Computers & Fluids (2005)]࡟ࡼࡾᐇ⿦ࡋ࡚࠸ࡿࠋࡲࡓࠊ᫬㛫✚ศ⤊஢ᚋ࡟ᩘ್᣺ືᢚไࡢࡓࡵ࡟ Hermite WeightedEssentially Non-Oscillatory (HWENO) ἲ [J. Qiu & C.-W.

Shu,

Journal of ComputationalPhysics (2003)]ࢆ౑⏝ࡋ࡚࠸ࡿࠋࡉࡽ࡟ࠊ⌧ᅾࡢࢥ࣮ࢻ࡟ࡣHWENOἲࢆ㝖ࡃ᫬㛫Ⓨᒎ㒊ศ࡟㏆ఝࢽ࣮ࣗࢺࣥἲ࡟ᇶ࡙ࡃእ㒊཯᚟ࡀ᪤࡟ᐇ⿦ࡉࢀ࡚࠸ࡿࠋࢥ࣮ࢻ⾜ᩘࡣ2,376⾜࡛࠶ࡿࠋ2.4 せồ㡯┠2.4.1 ୍⯡஦㡯1 ḟඖ᢬ᢠᛶ☢Ẽὶయࢥ࣮ࢻᮏయࡢᨵⰋసᴗࡣࣉࣟࢢ࣑ࣛࣥࢢ࡟ࡣ Fortran ゝㄒ(Fortran90௨㝆ࡢ௙ᵝ)ࢆ⏝࠸ࡿࡇ࡜ࠋᩘ್ィ⟬ࣛ࢖ࣈࣛࣜࡢ౑⏝ࡣㄆࡵ࡞࠸ࠋ2.4.2 ࢽ࣮ࣗࢺࣥἲ࡟ࡼࡾᒁᡤ⥺ᙧ໬ࡋࡓ㞳ᩓ᪉⛬ᘧࡢᑟฟཷὀ⪅ࡣࠊ⌧⾜ࡢ1ḟඖ᢬ᢠᛶ☢Ẽὶయࢥ࣮ࢻ࡟ᐇ⿦ࡉࢀ࡚࠸ࡿ㠀⥺ᙧ㞳ᩓ᪉⛬ᘧࢆྵࡵࡓࡶࡢ࡟ಟṇࡋࠊ㠀⥺ᙧ࣑ࣜࢵࢱ࣮ࢆྵࡵ࡚ᒁᡤ⥺ᙧ໬ࡋࡓ᪉⛬ᘧࡢᑟฟࢆ⾜࠺ࠋ㝜ゎἲ࡜ࡋ࡚ࠊ1ẁ2ḟ⢭ᗘ㝜ⓗࣝࣥࢤ࣭ࢡࢵࢱἲࢆ᥇⏝ࡍࡿࠋ᫬㛫Ⓨᒎ᪉⛬ᘧࡢྑ㎶࡟ࡣࠊHWENOἲ࡟ࡼࡾ⿵ṇࢆ⾜ࡗࡓኚᩘࡀධࡿࡇ࡜࡟ὀពࡍࡿࡇ࡜㸦ูῧ㸲ཧ↷㸧ࠋ⌧⾜ࡢ1ḟඖ᢬ᢠᛶ☢Ẽὶయࢥ࣮ࢻࡢᐇ⿦࡛ࡣᩘ್ὶ᮰ࡢィ⟬࡟㞄᥋ࢭࣝࡢ᝟ሗࡋ࠿౑⏝ࡋ࡞࠸ࡓࡵࠊᚓࡽࢀࡿಀᩘ⾜ิࡣࣈࣟࢵࢡ୕㔜ᑐゅ⾜ิ࡜࡞ࡿࠋཷὀ⪅ࡣࠊ㈚୚≀ရࡢࣈࣟࢵࢡ୕㔜ᑐゅ⾜ิࢯࣝࣂ࣮ࢆ౑⏝ࡋ࡚⥺ᙧ᪉⛬ᘧࢆゎࡃࡇ࡜ࡀ࡛ࡁࡿࠋ2.4.3 ᕪศ㏆ఝࡋࡓࣖࢥࣅ⾜ิ࡟ᇶ࡙ࡃࢽ࣮ࣗࢺࣥἲࡢᐇ⿦సᴗཷὀ⪅ࡣࠊᕪศ㏆ఝ࡟ࡼࡾࣖࢥࣅ⾜ิࡢィ⟬ࢆ⾜࠸1ḟඖ᢬ᢠᛶ☢Ẽὶయࢥ࣮ࢻ࡬ࡢᐇ⿦ࢆ⾜࠺ࠋHLLDἲ➼࡟ࡼࡗ࡚ィ⟬ࡉࢀࡿᩘ್ὶ᮰࡟ᑐࡋ࡚ࣖࢥࣅ⾜ิࢆゎᯒⓗ࡟ィ⟬ࡍࡿࡇ࡜ࡣ࡛ࡁ࡞࠸ࡓࡵࠊ≀⌮㔞࡟ᚤᑠኚ఩ࢆ୚࠼࡚㠀⥺ᙧ㡯ࢆィ⟬ࡋ㸪ࡑࢀࢆኚ఩㔞࡟ࡼࡾ㝖⟬ࡍࡿࡇ࡜࡛㏆ఝⓗ࡟ࣖࢥࣅ⾜ิࢆィ⟬ࡍࡿࠋ2.4.4 ᩘ್ᐇ㦂ཷὀ⪅ࡣࠊ㔞◊ࡀᣦᐃࡍࡿࢸࢫࢺၥ㢟࡟ࡼࡿᩘ್ᐇ㦂ࢆ⾜࠺ࠋᨵⰋ๓࡜ᨵⰋᚋࡢ࢔ࣝࢦࣜࢬ࣒࡟㛵ࡍࡿィ⟬㔞ࡢẚ㍑ࢆ⾜࠺ࠋࡲࡓࠊᨵⰋᚋࡢ࢔ࣝࢦࣜࢬ࣒ࡢ᪉ࡀᩘ್Ᏻᐃᛶࡀྥୖࡋ࡚࠸ࡿࡇ࡜ࡀண᝿ࡉࢀࡿࡓࡵࠊᩘ್ᐇ㦂⤖ᯝ࡟࠾ࡅࡿὶయሙࡢศᕸ࡟ࡘ࠸࡚㆟ㄽࢆ⾜࠺ࠋᚲせࡀ࠶ࢀࡤࢹࣂࢵࢢసᴗࢆ⾜࠺ࠋ2.5 ᥦฟ㈨ᩱࡢసᡂ1) ᥦฟ㈨ᩱࡢసᡂసᴗሗ࿌᭩࡟ࡣࠊᮏ௙ᵝ࡟ᇶ࡙ࡁࠊཷὀ⪅ࡀᐇ᪋ࡋࡓసᴗෆᐜཬࡧືసヨ㦂ࡢ⤖ᯝࢆグ㏙ࡍࡿࡇ࡜ࠋ2) ᥦฟ㈨ᩱ࡟㛵ࢃࡿ㟁Ꮚࣇ࢓࢖ࣝసᴗሗ࿌᭩ࡣ⣬፹యཬࡧCD-R➼ࡢ࣓ࢹ࢕࢔፹య࡜ࡋ࡚ᥦฟࡍࡿࠋ㛤Ⓨࡋࡓࢥ࣮ࢻ୍ᘧ࡜ྜࢃࡏ࡚ࠊືసヨ㦂࡟౑⏝ࡋࡓࢥࣥࣃ࢖ࣛ࢜ࣉࢩࣙࣥࡢタᐃࣇ࢓࢖ࣝࡸධຊࣇ࢓࢖ࣝࠊࢪࣙࣈᢞධࢫࢡࣜࣉࢺ㢮ࡶྜࢃࡏ࡚⣡ධࡍࡿࡇ࡜ࠋูῧࢥࣥࣆ࣮ࣗࢱࣉࣟࢢ࣒ࣛసᡂ➼ᴗົ≉⣙᮲㡯㸦┠ⓗ≀㸧➨㸯᮲ࡇࡢዎ⣙ࡢ┠ⓗ≀ࡣࠊḟࡢྛྕࡢ୍ཪࡣ஧௨ୖࡢ⤌ࡳྜࡏ࡟ヱᙜࡍࡿࢥࣥࣆ࣮ࣗࢱࣉࣟࢢ࣒ࣛࡢⴭస≀㸦ࢹ࣮ࢱࠊࢹ࣮ࢱ࣮࣋ࢫࠊ࣐ࢽࣗ࢔ࣝཬࡧࢻ࣓࢟ࣗࣥࢸ࣮ࢩࣙࣥࢆྵࡴࠋ௨ୗྠࡌࠋ㸧ཬࡧᙜヱࢥࣥࣆ࣮ࣗࢱࣉࣟࢢ࣒ࣛ࡟ࡼࡿィ⟬⤖ᯝ࡛࠶ࡗ࡚ࠊ௙ᵝ᭩࡟ᐃࡵࡿ⠊ᅖࡢࡶࡢ࡜ࡍࡿࠋ㸧ࢥࣥࣆ࣮ࣗࢱࣉࣟࢢ࣒ࣛ㸦ࢥࣥࣆ࣮ࣗࢱࣉࣟࢢ࣒ࣛࡢタィࢆྵࡴࠋ㸧ⴭస≀㸧㔞◊ࡀᥦ౪ࡍࡿࢥࣥࣆ࣮ࣗࢱࣉࣟࢢ࣒ࣛࡢⴭస≀࡟ࡼࡾᚓࡽࢀࡓィ⟬⤖ᯝ㸧ཷὀ⪅ࡀᡤ᭷ࡍࡿࢥࣥࣆ࣮ࣗࢱࣉࣟࢢ࣒ࣛࡢⴭస≀ཬࡧࡇࢀ࡟ࡼࡾᚓࡽࢀࡓィ⟬⤖ᯝ㸦ᶒ฼⩏ົࡢㆡΏ➼㸧➨㸰᮲ ཷὀ⪅ࡣࠊࡇࡢዎ⣙࡟ᇶ࡙ࡃᶒ฼ཪࡣ⩏ົࢆࠊ➨୕⪅࡟ㆡΏࡋࠊⱝࡋࡃࡣᢎ⥅ࡉࡏࠊཪࡣᢸಖࡑࡢ௚ࡢ┠ⓗ࡟౪ࡋ࡚ࡣ࡞ࡽ࡞࠸ࠋࡓࡔࡋࠊ࠶ࡽ࠿ࡌࡵ᭩㠃࡟ࡼࡾ㔞◊ࡢᢎㄆࢆᚓࡓሙྜࡣࠊࡇࡢ㝈ࡾ࡛ࡣ࡞࠸ࠋ㸦ᶒ฼ࡢᖐᒓ➼㸧➨㸱᮲ ࡇࡢᴗົ࡟ࡼࡾసᡂࡉࢀࡓ┠ⓗ≀㸦➨᮲ྛྕ࡟ᥖࡆࡿࡶࡢࢆ࠸࠺ࠋ௨ୗྠࡌࠋ㸧࡟ಀࡿⴭసᶒࡑࡢ௚ࡇࡢ┠ⓗ≀ࡢ౑⏝ࠊ཰┈ཬࡧฎศ㸦」〇ࠊ⩻ヂࠊ⩻᱌ࠊኚ᭦ࠊㆡΏ࣭㈚୚ཬࡧ஧ḟⓗⴭస≀ࡢ฼⏝ࢆྵࡴࠋ㸧࡟㛵ࡍࡿ୍ษࡢᶒ฼ࡣ㔞◊࡟ᖐᒓࡍࡿࡶࡢ࡜ࡍࡿࠋࡓࡔࡋࠊᮏዎ⣙㐙⾜ࡢࡓࡵ࡟౑⏝ࡍࡿࣉࣟࢢ࣒ࣛ➼ࡢ࠺ࡕࠊᮏዎ⣙⥾⤖௨๓࠿ࡽࠊཷὀ⪅ࡀᡤ᭷ࡍࡿࡶࡢ࡟ࡘ࠸࡚ࡣࠊࡑࡢⴭసᶒࡣཷὀ⪅࡟ᖐᒓࡍࡿࡶࡢ࡜ࡍࡿࠋ㸦Ặྡࡢ⾲♧ࡢไ㝈㸧➨㸲᮲ ཷὀ⪅ࡣࠊ➨᮲࡟つᐃࡍࡿⴭస≀࡟ⴭస⪅Ặྡࢆ⾲♧ࡋ࡞࠸ࡶࡢ࡜ࡍࡿࠋ㸦➨୕⪅ࡢᶒ฼ࡢಖㆤ㸧➨㸳᮲ཷὀ⪅ࡣࠊࡇࡢᴗົࡢᐇ᪋࡟㛵ࡋ➨୕⪅㸦ⴭస⪅ࢆྵࡴࠋ㸧ࡢⴭసᶒࡑࡢ௚ࡢᶒ฼ࢆ౵ᐖࡍࡿࡇ࡜ࡢ࡞࠸ࡼ࠺ᚲせ࡞ᥐ⨨ࢆ⮬ࡽࡢ㈐௵࡟࠾࠸࡚ㅮࡌ࡞ࡅࢀࡤ࡞ࡽ࡞࠸ࠋ㸦ᢏ⾡᝟ሗ㸧➨㸴᮲ 㔞◊ࡀࠊࡇࡢᴗົࡢᐇ᪋࡟㛵ࡋࠊཷὀ⪅ࡢಖ᭷ࡍࡿᢏ⾡᝟ሗࢆ▱ࡿᚲせࡀ⏕ࡌࡓሙྜ࡟ࡣࠊཷὀ⪅ࡣࠊࡇࡢዎ⣙ࡢᴗົ࡟ᚲせ࡞⠊ᅖෆ࡟࠾࠸࡚ᙜヱᢏ⾡᝟ሗࢆ㔞◊࡟↓ൾ࡛ᥦ౪ࡋ࡞ࡅࢀࡤ࡞ࡽ࡞࠸ࠋ㸰 㔞◊ࡣࠊཷὀ⪅࠿ࡽࡢ᭩㠃࡟ࡼࡿ஦๓ࡢྠពࢆᚓࡓሙྜࢆ㝖ࡁ๓㡯࡟ࡼࡾ▱ࡾᚓࡓᢏ⾡᝟ሗࢆ➨୕⪅࡟ᥦ౪ࡋ࡞࠸ࡶࡢ࡜ࡍࡿࠋ㸦ࣉࣟࢢ࣒ࣛ㛤Ⓨ࡟ᚲせ࡞ᢏ⾡᝟ሗ㸧➨㸵᮲ 㔞◊ࡣࠊ௙ᵝ᭩࡟ᐃࡵࡿ࡜ࡇࢁ࡟ࡼࡾࠊཷὀ⪅ࡀࡇࡢᴗົࡢᐇ᪋࡟ᚲせ࡞ィ⟬ࢥ࣮ࢻࡑࡢ௚ᚲせ࡞ᢏ⾡᝟ሗࢆཷὀ⪅࡟౑⏝ࡉࡏࡿࡇ࡜ࡀ࠶ࡿࠋ㸦༴㝤㈇ᢸ㸧➨㸶᮲ ࡇࡢᴗົࡢ᏶஢๓࡟┠ⓗ≀࡟ࡘ࠸࡚⏕ࡌࡓᦆᐖࡑࡢ௚ࡇࡢᴗົࡢᐇ᪋࡟㛵ࡋ࡚⏕ࡌࡓᦆᐖࡣࠊཷὀ⪅ࡢ㈇ᢸ࡜ࡍࡿࠋࡓࡔࡋࠊࡑࡢᦆᐖࡢ࠺ࡕཷὀ⪅ࡢ㈐ࡵ࡟ᖐࡋ㞴࠸஦⏤࡟ࡼࡾ⏕ࡌࡓࡶࡢ࡟ࡘ࠸࡚ࡣࠊࡇࡢ㝈ࡾ࡛࡞࠸ࠋูῧ㸰ᮏዎ⣙࡟࠾࠸࡚㑂Ᏺࡍ࡭ࡁࠕ᝟ሗࢭ࢟ࣗࣜࢸ࢕ࡢ☜ಖࠖ࡟㛵ࡍࡿ஦㡯㸯㸬ཷὀ⪅ࡣࠊዎ⣙ࡢᒚ⾜࡟㛵ࡋࠊ᝟ሗࢩࢫࢸ࣒㸦᝟ሗฎ⌮ཬࡧ㏻ಙ࡟㛵ࢃࡿࢩࢫࢸ࣒࡛࠶ࡗ࡚ࠊࣁ࣮ࢻ࢙࢘࢔ࠊࢯࣇࢺ࢙࢘࢔ཬࡧࢿࢵࢺ࣮࣡ࢡ୪ࡧ࡟グ㘓፹య࡛ᵓᡂࡉࢀࡿࡶࡢࢆ࠸࠺ࠋ㸧ࢆ฼⏝ࡍࡿሙྜ࡟ࡣࠊ㔞◊ࡢ᝟ሗཬࡧ᝟ሗࢩࢫࢸ࣒ࢆಖㆤࡍࡿࡓࡵ࡟ࠊ᝟ሗࢩࢫࢸ࣒࠿ࡽࡢ᝟ሗ₃ὤࠊࢥࣥࣆ࣮ࣗࢱ࢘࢕ࣝࢫࡢ౵ධ➼ࡢ㜵Ṇࡑࡢ௚ᚲせ࡞ᥐ⨨ࢆㅮࡌ࡞ࡅࢀࡤ࡞ࡽ࡞࠸ࠋ㸰㸬ཷὀ⪅ࡣࠊḟࡢྛྕ࡟ᥖࡆࡿ஦㡯ࢆ㑂Ᏺࡍࡿ࡯࠿ࠊ㔞◊ࡢ᝟ሗࢭ࢟ࣗࣜࢸ࢕☜ಖࡢࡓࡵ࡟ࠊ㔞◊ࡀᚲせ࡞ᣦ♧ࢆ⾜ࡗࡓ࡜ࡁࡣࠊࡑࡢᣦ♧࡟ᚑࢃ࡞ࡅࢀࡤ࡞ࡽ࡞࠸ࠋ 1.(1) ཷὀ⪅ࡣࠊዎ⣙ࡢᴗົ࡟ᦠࢃࡿ⪅㸦௨ୗࠕᴗົᢸᙜ⪅ࠖ࡜࠸࠺ࠋ㸧ࢆ≉ᐃࡋࠊࡑࢀ௨እࡢ⪅࡟సᴗࢆࡉࡏ࡚ࡣ࡞ࡽ࡞࠸ࠋ 1.(2) ཷὀ⪅ࡣࠊዎ⣙࡟㛵ࡋ࡚▱ࡾᚓࡓ᝟ሗ㸦㔞◊࡟ᘬࡁΏࡍ࡭ࡁࢥࣥࣆ࣮ࣗࢱࣉࣟࢢ࣒ࣛⴭస≀ཬࡧィ⟬⤖ᯝࢆྵࡴࠋ௨ୗྠࡌࠋ㸧ࢆྲྀࡾᢅ࠺᝟ሗࢩࢫࢸ࣒࡟ࡘ࠸࡚ࠊᴗົᢸᙜ⪅௨እࡀᙜヱ᝟ሗ࡟࢔ࢡࢭࢫྍ⬟࡜࡞ࡽ࡞࠸ࡼ࠺㐺ษ࡟࢔ࢡࢭࢫไ㝈ࢆ⾜࠺ࡇ࡜ࠋ 1.(3) ཷὀ⪅ࡣࠊዎ⣙࡟㛵ࡋ࡚▱ࡾᚓࡓ᝟ሗࢆྲྀࡾᢅ࠺᝟ሗࢩࢫࢸ࣒࡟ࡘ࠸࡚ࠊ࢘࢕ࣝࢫᑐ⟇ࢶ࣮ࣝཬࡧࣇ࢓࢖࢔࢛࣮࢘ࣝᶵ⬟ࡢᑟධࠊࢭ࢟ࣗࣜࢸ࢕ࣃࢵࢳࡢ㐺⏝➼㐺ษ࡞᝟ሗࢭ࢟ࣗࣜࢸ࢕ᑐ⟇ࢆᐇ᪋ࡍࡿࡇ࡜ࠋ 1.(4) ཷὀ⪅ࡣࠊP2P ࣇ࢓࢖ࣝ஺᥮ࢯࣇࢺ࢙࢘࢔㸦WinnyࠊWinMXࠊKaZaaࠊShare ➼㸧ཬࡧ SoftEtherࢆᑟධࡋࡓ᝟ሗࢩࢫࢸ࣒࡟࠾࠸࡚ࠊዎ⣙࡟㛵ࡋ࡚▱ࡾᚓࡓ᝟ሗࢆྲྀࡾᢅࡗ࡚ࡣ࡞ࡽ࡞࠸ࠋ 1.

(5) ཷὀ⪅ࡣࠊ㔞◊ࡢᢎㅙࡢ࡞࠸㝈ࡾࠊዎ⣙࡟㛵ࡋ࡚▱ࡾᚓࡓ᝟ሗࢆ㔞◊ཪࡣཷὀ⪅ࡢ᝟ሗࢩࢫࢸ࣒௨እࡢ᝟ሗࢩࢫࢸ࣒㸦ᴗົᢸᙜ⪅ࡀᡤ᭷ࡍࡿࣃࢯࢥࣥ➼㸧࡟࠾࠸࡚ྲྀࡾᢅࡗ࡚ࡣ࡞ࡽ࡞࠸ࠋ 1.(6) ཷὀ⪅ࡣࠊጤ௵ࢆࡋཪࡣୗㄳ㈇ࢆࡉࡏࡓሙྜࡣࠊᙜヱጤ௵ཪࡣୗㄳ㈇ࢆཷࡅࡓ⪅ࡢዎ⣙࡟㛵ࡍࡿ⾜Ⅽ࡟ࡘ࠸࡚ࠊ㔞◊࡟ᑐࡋ඲࡚ࡢ㈐௵ࢆ㈇࠺࡜࡜ࡶ࡟ࠊᙜヱጤ௵ཪࡣୗㄳ㈇ࢆཷࡅࡓ⪅࡟ᑐࡋ࡚ࠊ᝟ሗࢭ࢟ࣗࣜࢸ࢕ࡢ☜ಖ࡟ࡘ࠸࡚ᚲせ࡞ᥐ⨨ࢆㅮࡎࡿࡼ࠺࡟ດࡵ࡞ࡅࢀࡤ࡞ࡽ࡞࠸ࠋ 1.(7) ཷὀ⪅ࡣࠊ㔞◊ࡀồࡵࡓሙྜ࡟ࡣࠊ᝟ሗࢭ࢟ࣗࣜࢸ࢕ᑐ⟇ࡢᐇ᪋≧ἣ࡟ࡘ࠸࡚ࡢ┘ᰝࢆཷࡅධࢀࠊࡇࢀ࡟༠ຊࡍࡿࡇ࡜ࠋ 1.(8) ཷὀ⪅ࡣࠊ㔞◊ࡢᥦ౪ࡋࡓ᝟ሗ୪ࡧ࡟ཷὀ⪅ཬࡧጤ௵ཪࡣୗㄳ㈇ࢆཷࡅࡓ⪅ࡀዎ⣙ᴗົࡢࡓࡵ࡟཰㞟ࡋࡓ᝟ሗ࡟ࡘ࠸࡚ࠊ⅏ᐖࠊ⣮ኻࠊ◚ቯࠊᨵࡊࢇࠊࡁᦆࠊ₃࠼࠸ࠊࢥࣥࣆ࣮ࣗࢱ࢘࢕ࣝࢫ࡟ࡼࡿ⿕ᐖࠊ୙ṇ࡞฼⏝ࠊ୙ṇ࢔ࢡࢭࢫࡑࡢ௚ࡢ஦ᨾࡀⓎ⏕ࠊཪࡣ⏕ࡎࡿ࠾ࡑࢀࡢ࠶ࡿࡇ࡜ࢆ▱ࡗࡓሙྜࡣࠊࡓࡔࡕ࡟㔞◊࡟ሗ࿌ࡋࠊ㔞◊ࡢᣦ♧࡟ᚑ࠺ࡶࡢ࡜ࡍࡿࠋዎ⣙ࡢ⤊஢ᚋ࡟࠾࠸࡚ࡶࠊྠᵝ࡜ࡍࡿࠋ࡞࠾ࠊ㔞◊ࡢධᮐ࡟ཧຍࡍࡿሙྜࠊࡲࡓࡣ㔞◊࠿ࡽࡢぢ✚౫㢗ࢆཷࡅࡿሙྜ࡟ࡶࠊୖグ஦㡯ࢆ㑂Ᏺࡋ࡚࠸ࡓࡔࡁࡲࡍࠋ௨ ୖูῧ㸱⏘ᴗ㈈⏘ᶒ≉⣙᮲㡯㸦ཷὀ⪅ࡀ༢⊂࡛⾜ࡗࡓⓎ᫂➼ࡢ⏘ᴗ㈈⏘ᶒࡢᖐᒓ㸧➨㸯᮲ ཷὀ⪅ࡣࠊᮏዎ⣙࡟㛵ࡋ࡚ࠊཷὀ⪅ࡀ༢⊂࡛࡞ࡋࡓⓎ᫂ཪࡣ⪃᱌㸦௨ୗࠕⓎ᫂➼ࠖ࡜࠸࠺ࠋ㸧࡟ᑐࡍࡿ≉チᶒࠊᐇ⏝᪂᱌ᶒཪࡣព໶ᶒ㸦௨ୗࠕ≉チᶒ➼ࠖ࡜࠸࠺ࠋ㸧ࢆྲྀᚓࡍࡿሙྜࡣࠊ༢⊂࡛ฟ㢪࡛ࡁࡿࡶࡢ࡜ࡍࡿࠋࡓࡔࡋࠊฟ㢪ࡍࡿ࡜ࡁࡣ࠶ࡽ࠿ࡌࡵฟ㢪࡟㝿ࡋ࡚ᥦฟࡍ࡭ࡁ᭩㢮ࡢ෗ࡋࢆῧ࠼࡚㔞◊࡟㏻▱ࡍࡿࡶࡢ࡜ࡍࡿࠋ㸦ཷὀ⪅ࡀ༢⊂࡛⾜ࡗࡓⓎ᫂➼ࡢ≉チᶒ➼ࡢㆡΏ➼㸧➨㸰᮲ ཷὀ⪅ࡣࠊཷὀ⪅ࡀ๓᮲ࡢ≉チᶒ➼ࢆ㔞◊௨እࡢ➨୕⪅࡟ㆡΏཪࡣᐇ᪋チㅙࡍࡿሙྜ࡟ࡣࠊᮏྲྀᢅ࠸ࡢྛ᮲㡯ࡢつᐃࡢ㐺⏝࡟ᨭ㞀ࢆ୚࠼࡞࠸ࡼ࠺ᙜヱ➨୕⪅࡜⣙ᐃࡋ࡞ࡅࢀࡤ࡞ࡽ࡞࠸ࠋ㸦ཷὀ⪅ࡀ༢⊂࡛⾜ࡗࡓⓎ᫂➼ࡢ≉チᶒ➼ࡢᐇ᪋チㅙ㸧➨㸱᮲ 㔞◊ࡣࠊ➨㸯᮲ࡢⓎ᫂➼࡟ᑐࡍࡿ≉チᶒ➼ࢆ↓ൾ࡛⮬ࡽヨ㦂ཪࡣ◊✲ࡢࡓࡵ࡟ᐇ᪋ࡍࡿࡇ࡜ࡀ࡛ࡁࡿࠋ㔞◊ࡀ㔞◊ࡢࡓࡵ࡟ཷὀ⪅௨እࡢ➨୕⪅࡟〇సࡉࡏࠊཪࡣᴗົࢆ௦⾜ࡍࡿ➨୕⪅࡟෌ᐇ᪋ᶒࢆチㅙࡍࡿሙྜࡣࠊཷὀ⪅ࡢᢎㅙࢆᚓࡓୖ࡛チㅙࡍࡿࡶࡢ࡜ࡋࠊࡑࡢᐇ᪋᮲௳➼ࡣ㔞◊ࠊཷὀ⪅༠㆟ࡢୖࠊỴᐃࡍࡿࠋ㸦㔞◊ཬࡧཷὀ⪅ࡀඹྠ࡛⾜ࡗࡓⓎ᫂➼ࡢ≉チᶒ➼ࡢᖐᒓཬࡧ⟶⌮㸧➨㸲᮲ 㔞◊ཬࡧཷὀ⪅ࡣࠊᮏዎ⣙࡟㛵ࡋ࡚ඹྠ࡛࡞ࡋࡓⓎ᫂➼࡟ᑐࡍࡿ≉チᶒ➼ࢆྲྀᚓࡍࡿሙྜࡣࠊඹྠฟ㢪ዎ⣙ࢆ⥾⤖ࡋࠊඹྠ࡛ฟ㢪ࡍࡿࡶࡢ࡜ࡋࠊฟ㢪ࡢࡓࡵࡢ㈝⏝ࡣࠊ㔞◊ࠊཷὀ⪅ࡢᣢศ࡟ẚ౛ࡋ࡚㈇ᢸࡍࡿࡶࡢ࡜ࡍࡿࠋ㸦㔞◊ཬࡧཷὀ⪅ࡀඹྠ࡛⾜ࡗࡓⓎ᫂➼ࡢ≉チᶒ➼ࡢᐇ᪋㸧➨㸳᮲ 㔞◊ࡣࠊඹྠ࡛⾜ࡗࡓⓎ᫂➼ࢆヨ㦂ཪࡣ◊✲௨እࡢ┠ⓗ࡟ᐇ᪋ࡋ࡞࠸ࡶࡢ࡜ࡍࡿࠋࡓࡔࡋࠊ㔞◊ࡣ㔞◊ࡢࡓࡵ࡟ཷὀ⪅௨እࡢ➨୕⪅࡟〇సࡉࡏࠊཪࡣᴗົࢆ௦⾜ࡍࡿ➨୕⪅࡟ᐇ᪋チㅙࡍࡿሙྜࡣࠊ↓ൾ࡟࡚ᙜヱ➨୕⪅࡟ᐇ᪋チㅙࡍࡿࡇ࡜ࡀ࡛ࡁࡿࡶࡢ࡜ࡍࡿࠋ㸰ཷὀ⪅ࡀ๓㡯ࡢⓎ᫂➼࡟ࡘ࠸࡚⮬ࡽၟᴗⓗᐇ᪋ࢆࡍࡿ࡜ࡁࡣࠊ㔞◊ࡀ⮬ࡽၟᴗⓗᐇ᪋ࢆࡋ࡞࠸ࡇ࡜࡟㚷ࡳࠊཷὀ⪅ࡢၟᴗⓗᐇ᪋ࡢィ⏬ࢆ຺᱌ࡋࠊ஦๓࡟ᐇ᪋ᩱ➼࡟ࡘ࠸࡚㔞◊ࠊཷὀ⪅༠㆟ࡢୖࠊู㏵ᐇ᪋ዎ⣙ࢆ⥾⤖ࡍࡿࡶࡢ࡜ࡍࡿࠋ㸦⛎ᐦࡢಖᣢ㸧➨㸴᮲ 㔞◊ཬࡧཷὀ⪅ࡣࠊ➨㸯᮲ཬࡧ➨㸲᮲ࡢⓎ᫂➼ࡢෆᐜࢆฟ㢪࡟ࡼࡾෆᐜࡀබ㛤ࡉࢀࡿ᪥ࡲ࡛௚࡟₃࠼࠸ࡋ࡚ࡣ࡞ࡽ࡞࠸ࠋࡓࡔࡋࠊ࠶ࡽ࠿ࡌࡵ᭩㠃࡟ࡼࡾฟ㢪ࢆ⾜ࡗࡓ⪅ࡢ஢ゎࢆᚓࡓሙྜࡣࡇࡢ㝈ࡾ࡛ࡣ࡞࠸ࠋ㸦ጤ௵࣭ୗㄳ㈇㸧➨㸵᮲ ཷὀ⪅ࡣࠊᮏዎ⣙ࡢ඲㒊ཪࡣ୍㒊ࢆ➨୕⪅࡟ጤ௵ࡋࠊཪࡣㄳࡅ㈇ࢃࡏࡓሙྜ࡟࠾࠸࡚ࡣࠊࡑࡢ➨୕⪅࡟ᑐࡋ࡚ࠊᮏྲྀᢅ࠸ࡢྛ᮲㡯ࡢつᐃࢆ‽⏝ࡍࡿࡶࡢ࡜ࡋࠊཷὀ⪅ࡣࡇࡢࡓࡵ࡟ᚲせ࡞ᥐ⨨ࢆㅮࡌ࡞ࡅࢀࡤ࡞ࡽ࡞࠸ࠋ㸰ཷὀ⪅ࡣࠊ๓㡯ࡢᙜヱ➨୕⪅ࡀᮏྲྀᢅ࠸࡟ᐃࡵࡿ஦㡯࡟㐪཯ࡋࡓሙྜ࡟ࡣࠊ㔞◊࡟ᑐࡋ඲࡚ࡢ㈐௵ࢆ㈇࠺ࡶࡢ࡜ࡍࡿࠋ㸦༠㆟㸧➨㸶᮲ ➨㸯᮲ཬࡧ➨㸲᮲ࡢሙྜ࡟࠾࠸࡚ࠊ༢⊂ⱝࡋࡃࡣඹྠࡢ༊ูཪࡣඹྠࡢ⠊ᅖ➼࡟ࡘ࠸࡚␲⩏ࡀ⏕ࡌࡓ࡜ࡁࡣࠊ㔞◊ࠊཷὀ⪅༠㆟ࡋ࡚ᐃࡵࡿࡶࡢ࡜ࡍࡿࠋ㸦᭷ຠᮇ㛫㸧➨㸷᮲ ᮏྲྀᢅ࠸ࡢ᭷ຠᮇ㝈ࡣࠊዎ⣙⥾⤖ࡢ᪥࠿ࡽᙜヱ≉チᶒ➼ࡢᾘ⁛ࡍࡿ᪥ࡲ࡛࡜ࡍࡿࠋ௨ୖ 別添４多次元方程式に対するニュートン法次のような非線形システムを考える．F(Un+1)=0, (1)ここで，非線形方程式F及び数値解Un+1 の次元数は等しいものとし，nは時間ステップを表す．ニュートン法では次に示すように，局所線形化を行うことで反復的にUn+1 を得る．F(Un+1,k)+∂F∂Un+1(Un+1,k)(Un+1,k+1 −Un+1,k)=0, (2)ここで，k はニュートン法の内部ステップを表す．したがって，ヤコビ行列∂F/∂Un+1 を係数行列とする線形方程式を解くことで，Un+1 の修正量Un+1,k+1 −Un+1,k を得ることができる．また，近似リーマン解法のような非線形なアルゴリズムを使用する場合，∂F/∂Un+1 の関数形を解析的に得ることができない．そこで，∂F∂Un+1i=F(.,Un+1i +ε,.)−F(.,Un+1i −ε, . . . )2ε+O(ε2), (3)と差分近似することで計算を行う．ここで，εは微小な定数である．不連続Galerkin法では隣接セルの情報のみを使用してステンシルを構成することから，ヤコビ行列∂F/∂Un+1 は疎行列となり，計算量はO(N)となる．現在のコードでは，次の形式に従って陰解法の実装が行われている．非線形磁気流体方程式は，Un+1 −UnΔt=G(Un+1), (4)非線形リミッターは，Un+1 =H(Un+1). (5)ここで，Un+1 は時間発展後の数値解，Un+1 はフィルタリング処理後の数値解を表す．すなわち，ニュートン法による外部反復は式(4) のみに対して行われ，式(5) は陽的に計算可能であるためニュートン法に組み込まれていない．今回の業務では非線形磁気流体方程式及び非線形リミッターを次のように改修する．Un+1 −UnΔt=G(Un+1), (6)Un+1 =H(Un+1). (7)すなわち，時間発展方程式の右辺にUn+1 が含まれるため，改修後のコードでは非線形リミッター部分もまたニュートン法の中に組み込まれることになる．1