JP2005050027A - 流動解析システムおよび流動解析方法 - Google Patents

Abstract

【課題】流動解析において、個別機器の設計がシステム全体に及ぼす影響を把握し、かつ、システム全体挙動の一部分として個別機器の境界条件を適切に与える。
【解決手段】流体系全体解析モデル１０１と、流動スケールに応じて階層化された複数の部分流動解析モデル１０２とを用いて解析する流動解析システムである。全体システム挙動を解析する手段と、全体解析結果を部分解析モデルの境界条件またはモデルパラメーターとして利用し、階層化された部分解析モデルの上位のモデルによる解析結果を下位の部分解析モデルの境界条件またはモデルパラメーターとして利用して部分解析モデルによる解析を行なう手段と、下位の部分解析モデルによる解析結果に基づいて作成される物理モデルまたは相関式モデルまたは物理パラメーターを上位の部分解析モデル及び全体解析モデルに組み込む手段と、流体系全体解析モデルによる解析結果に基づいて作成される物理モデルまたは相関式モデルまたは物理パラメーターを流体系全体解析モデルに組み込む手段と、を有する。
【選択図】 図１

Description

【０００１】
【発明の属する技術分野】
本発明は、計算機を用いた流動解析システムおよび流動解析方法に関し、特に、原子炉などの複雑流路を有する大規模流体系に適した流動解析システムおよび流動解析方法に関する。
【０００２】
【従来の技術】
従来の原子炉系解析では、図１６に示すように管路網モデルなどの比較的計算機負荷の少ないモデルで原子炉系全体システムの挙動を解析し、個別機器性能に対する流動の影響や局所的流動状態の把握のために流動解析が用いられている。原子炉内の冷却材流動挙動は個別機器の性能のみならずシステム全体の性能にも影響を及ぼすため、システム全体に亘る流動解析が望まれるが、多数の機器から構成される複雑大規模流体系である原子炉系全体の詳細流動解析は実現していない。
【０００３】
このような状況の中、複雑形状における流動解析手法が特許文献１に提案されている。特許文献１では、解析対象を微小要素に分割した集合を複数の部分集合領域に分割し、並列計算機を用いて、各部分集合領域ごとに別々の演算処理装置により時間変化に伴う流体の流動状態を計算するというものである。
【０００４】
また、大規模な非定常解析を短時間に効率的に実行する方法が特許文献２に提案されている。特許文献２は、解析対象領域を複数の微小要素に分割した解析モデルを構築し、解析モデルの解析領域内に状態量の時間変化の激しい微小要素を含む領域Ａ及び状態量の時間変化の緩やかな微小要素を含む領域Ｂを設定し、領域Ｂについては時間的な計算負荷の小さな計算手法により状態量を求めて状態量Ｂとし、得られた状態量Ｂに基づいて境界条件を設定し、領域Ａについて時間的な計算負荷の大きな計算手法により求めるという方法である。
【０００５】
更に、大規模システムを階層的に分解し、分散して設計（解析）を行なう方法が特許文献４に提案されている。特許文献４では、システムを複数の下位システムに分割して、下位システムを各々のコンピューターを利用してシミュレーションで個別に設計し、各コンピューターは通信装置等で接続され、設計に関する情報を送受する構成となっており、下位システムの設計をまとめる形で上位システムの設計が行なわれる。
【０００６】
これに類似した技術が特許文献３に提案されている。特許文献３では、遠隔地に存在する複数のシミュレーター装置をシミュレーションプラットフォーム上に接続し、これらのシミュレーター装置を複数のグループに分割して管理、実効制御するものである。また、複数のシミュレーションプログラムを同時に実行処理する方法が特許文献５に提案されている。
【０００７】
【特許文献１】
特開２００２−１７２６６８号公報
【特許文献２】
特開２００２−１５７２８６号公報
【特許文献３】
特開平１１−１６１６３０号公報
【特許文献４】
特開２００２−１５００９号公報
【特許文献５】
特開２０００−２６７８９０号公報
【０００８】
【発明が解決しようとする課題】
上述の従来の原子炉系解析においては、個別機器の性能評価と原子炉系全体システム挙動の評価がそれぞれ独立して行なわれる。そのため、個別機器間の境界条件に含まれる不確定性あるいは、流動解析モデルと管路網モデルのモデルの違いに起因する誤差あるいは解析できない現象の存在等が機器設計あるいはシステム設計に安全余裕を見込んだ“設計裕度”として取り込まれ、原子炉系システム全体として過度の設計裕度を生じている。
【０００９】
このような過度の設計裕度を適正化し、原子炉系システムの合理化を図るためには、個別機器の設計がシステム全体に及ぼす影響を正確に把握し、なおかつ、システム全体挙動の一部分として個別機器の境界条件を適切に与える必要がある。
【００１０】
これらの問題を解決するための一つの策として、原子炉系全体の詳細な流動解析モデルを作成し、一貫して解析する方法が考えられるが、複雑流路形状を有し大規模領域を対象とする原子炉系解析では計算格子数が膨大になり、計算機のメモリ容量、解析結果のデータ容量、解析時間などの観点から実現は難しい。
【００１１】
従来技術の例として提示した特許文献１では複雑形状における流動解析手法が提案されている。しかしながら、これはいわゆる並列計算における一般的な領域分割法に他ならず、仮に原子炉系全体の数値流動解析モデルにこの手法を適用した場合、１つの部分集合領域の計算時間を現実的な値すなわち数日のオーダーに抑えようとすれば、部分集合領域数は１０^４のオーダーに上る。つまり、１０^４のオーダーの演算処理装置を並列計算させることになる。従って、特許文献１の方法では、原子炉系全体の数値流動解析を現実的な時間及び演算処理装置数で行なうことは現時点で不可能である。
【００１２】
また、特許文献２では、大規模な非定常解析を短時間に効率的に実行する方法が提案されているが、これを沸騰水型原子炉系の解析に当てはめて考えた場合、沸騰水型原子炉系では二相流領域の占める割合が大きく、更に流路面積の大きな場所からその数百分の一、数千分の一の流路面積を有する複数の狭隘並列流路へ分配される部分が存在し、特許文献２の中で述べている状態量の時間変化の激しい領域Ａに相当する部分の割合が非常に大きくなる。
【００１３】
領域Ａについては、時間的な計算負荷の大きな解析手法を用いるあるいは、時間ステップを小さく取った解析を行なうことが提案されているが、その場合、領域Ａの解析時間が全体解析の中における律速となる。原子炉系全体の詳細な流動解析において、その実現を困難にしているのがまさに上述の領域Ａに相当する部分の解析負荷の大きさである。しかしながら、特許文献２では、この領域Ａの解析負荷自体を低減する提案は全くなされておらず、原子炉系全体の数値流動解析を実現するための手段とはなり得ない。
【００１４】
一方、特許文献４では、大規模システムを階層的に分解し、分散して設計（解析）を行なう方法が提案されている。具体的には大規模機械を構成するシステムを複数の下位システムに分解し、この下位システムを更にメインシステムとそれに付随するサブシステムに分類する。サブシステムは個別に設計され、それを統合する形でメインシステムが構成されるというものである。この概念は、原子炉系全体の数値流動解析を実現する上でも参考となるものである。
【００１５】
そこで、特許文献４の方法を原子炉系全体の数値流動解析に当てはめてみると、メインシステムが原子炉流動系となり、サブシステムは制御系、核特性系、流体流れ詳細となる。このサブシステムの流体流れ詳細の部分を更に、解析領域ごとにメインシステムとサブシステムに分類すると、原子炉系全体流動というメインシステムの中に炉心―上部プレナム―気水分離機流動、蒸気乾燥器―蒸気ドーム流動、原子炉水位近傍流動、給水混合部―ダウンカマ―再循環ポンプ―下部プレナム流動、などの複数のサブシステムに分かれる。
【００１６】
こうしたサブシステムの流動状態は更に下位の機器レベルの流動状態に分類できる。簡単な例では、気水分離機流動は気水分離機一体を一つのサブシステムとした時、この気水分離機が数百体集まって構成している一つのメインシステムである。原子炉系流動解析では、これらの全ての階層におけるサブシステム及びメインシステム全体が、共通の流体場を介して相互に密接に関連している。
【００１７】
このような相互作用をフィードバックさせずに、原子炉系を階層的に分解してそれぞれのサブシステムを個別に解析し、その結果だけをメインシステムに戻すのであれば、個別機器の性能評価と原子炉系全体システム挙動の評価がそれぞれ独立して行なわれている従来の原子炉系解析と変わるところはない。逆に、それぞれのサブシステム及びメインシステムで情報を相互にやり取りしながら詳細流動解析を進めるのであれば、それは、原子炉系全体の詳細流動解析を単に個別システムとして分散させて実行するものとなり、やはり解析時間及び収束性などの観点から実現は難しい。
【００１８】
このように、上位システムの詳細流動解析が下位システムの詳細流動解析の集合体として位置付けられる限り、原子炉系全体の詳細流動解析を実用化できる時間スケール及び並列演算処理装置数で実施することは現時点で不可能である。特許文献４においても、流動解析に関するメインシステムとサブシステムの関係は特に明示されておらず、この問題に対する解決策は見出せない。
【００１９】
本発明は上記従来技術の課題を解決するためになされたものであり、原子炉系などの大規模流動解析において、個別機器の設計がシステム全体に及ぼす影響を正確に把握し、なおかつ、システム全体挙動の一部分として個別機器の境界条件を適切に与えることのできる流動解析システムおよび流動解析方法を得ることを目的とする。
【００２０】
【課題を解決するための手段】
この発明は上記目的に沿うものであって、請求項１に記載の発明は、流体系全体を複数のノードに分割し、ノード間をジャンクションで接続し、該ノード内で状態量を計算し、該ジャンクションでノード間の移行量を計算する最上位の数値流動解析モデルとしての流体系全体解析モデルと、解析対象領域によって前記流体系全体を分割してかつ流動スケールに応じて階層化された複数の部分流動解析モデルとを用いて、流体系の流動を解析する流動解析システムにおいて、前記流体系全体解析モデルを用いて流体系全体のシステム挙動を解析する手段と、前記流体系全体解析モデルによる解析結果を前記複数の部分流動解析モデルの境界条件あるいはモデルパラメーターとして利用し、前記流動スケールに応じて階層化された部分流動解析モデルの上位のモデルによる解析結果をそれより下位の部分流動解析モデルの境界条件あるいはモデルパラメーターとして利用して前記複数の部分流動解析モデルによる解析を行なう手段と、下位の部分流動解析モデルを用いた解析結果に基づいて作成される物理モデルまたは相関式モデルまたは物理パラメーターを上位の部分流動解析モデル及び前記流体系全体解析モデルの一部として組み込む手段と、前記流体系全体解析モデルを用いた解析結果に基づいて作成される物理モデルまたは相関式モデルまたは物理パラメーターを前記流体系全体解析モデルの一部として組み込む手段と、を有することを特徴とする。
【００２１】
また、請求項５に記載の発明は、前記流動解析システムにおいて、前記複数の部分流動解析モデルの下位に、該部分流動解析モデルに領域として内包される局所領域の局所流動解析モデルがあって、前記流体系全体解析モデルもしくは前記部分流動解析モデルの解析領域に対する最適化問題の独立変数をｘ１、ｘ２、・・・ｘｍとし、目的関数をＦ（ｘ１、ｘ２、・・・ｘｍ）とし、独立変数ｘｍの実行可能領域をＳｍとし、それよりも下位の流動解析モデルの解析領域に対する最適化問題の独立変数をｙ１、ｙ２、・・・ｙｎとし、目的関数をＧ（ｙ１、ｙ２、・・・ｙｎ）とし、独立変数ｙｎの実行可能領域をＴｎとするとき、該解析モデルにおいて、任意のｘｍおよびｙｎは解析パラメーターとして設定自在に構成され、かつｘｍとｙｎの関係を解析によって算出可能にモデル化することによって、その関係に基づいてＳｍとＴｎのいずれか一方から他方を求め、目的関数Ｆ（ｘ１、ｘ２、・・・ｘｍ）もしくは、目的関数Ｇ（ｙ１、ｙ２、・・・ｙｎ）の最適解が得られることが可能に、前記流体系全体解析モデルおよび部分流動解析モデルが構築されていること、を特徴とする。
【００２２】
また、請求項１０に記載の発明は、流体系全体を複数のノードに分割し、ノード間をジャンクションで接続し、該ノード内で状態量を計算し、該ジャンクションでノード間の移行量を計算する最上位の流動解析モデルとしての流体系全体解析モデルと、解析対象領域によって前記流体系全体を分割してかつ流動スケールに応じて階層化された複数の部分流動解析モデルとを用いて、流体系の流動を解析する流動解析方法において、前記流体系全体解析モデルを用いて流体系全体のシステム挙動を解析するステップと、前記流体系全体解析モデルによる解析結果を前記複数の部分流動解析モデルの境界条件あるいはモデルパラメーターとして利用し、前記流動スケールに応じて階層化された部分流動解析モデルの上位のモデルによる解析結果をそれより下位の部分流動解析モデルの境界条件あるいはモデルパラメーターとして利用して前記複数の部分流動解析モデルによる解析を行なうステップと、前記下位の部分流動解析モデルを用いた解析結果に基づいて作成される物理モデルまたは相関式モデルまたは物理パラメーターを上位の部分流動解析モデル及び前記流体系全体解析モデルの一部として組み込むステップと、前記流体系全体解析モデルを用いた解析結果に基づいて作成される物理モデルまたは相関式モデルまたは物理パラメーターを前記流体系全体解析モデルの一部として組み込むステップと、を有することを特徴とする。
【００２３】
【発明の実施の形態】
図１〜図１０は本発明の大規模流動解析システムを利用して炉心への流量配分を評価する場合を例題とした、第１の実施の形態の説明図である。
解析対象である沸騰水型原子炉では、図２に示すように、原子炉圧力容器２内に複数の上向き平行流路である炉心チャンネル４が配置され、さらにこの炉心チャンネル４に平行にバイパス領域６が設けられている。各炉心チャンネル４の入口部には、炉心入口オリフィス１６が配置されていて、原子炉圧力容器２の最下部に形成された下部プレナム１８から冷却材が流入するようになっている。また、冷却材は、下部プレナム１８では液単相であるが、炉心チャンネル４を通るときに加熱されて沸騰し、気液二相流となる。
【００２４】
炉心チャンネル４およびバイパス領域６の上方には上部プレナム８が形成され、その上方に複数の気水分離器８が配置されている。気水分離器８の上方にはベッセルドーム１０が形成されている。炉心チャンネル４、バイパス領域６および上部プレナム８は筒状のシュラウド１２で囲まれ、シュラウド１２の外側にダウンカマ１４が形成されている。ダウンカマ１４の底部には複数（例えば１０機）の再循環ポンプ１６がほぼ等間隔に配置されている。再循環ポンプ２０によって、ダウンカマ１４内の冷却材が下部プレナム１８に圧入される。下部プレナム１８には、制御棒案内管２２等の管群が配置されている。
【００２５】
炉心への流量配分に影響を及ぼす要因には、次のようなものがある。
（１）周方向に均等配置された１０機の再循環ポンプの吐出流量
（２）下部プレナム領域での管群間隙での冷却材流動
（３）炉心支持板下部に設置されたビームによる流れの撹乱
（４）炉心入口オリフィス部での渦流の発生と圧力損失
（５）炉心部における流動と相変化による圧力損失
【００２６】
これらはいずれも冷却材の流動に起因するものであるが、流動が生じている空間の規模、流速変化の時間スケールは現象毎に異なっている。本実施の形態の大規模流動解析システムは、流動の空間的、時間的なスケールに合わせて、複数の流動解析プログラム（モデル）を準備している。上述の炉心への流量配分に影響を及ぼす要因に対しては、ユーザーは大規模流動解析システムにおいて次の２つの解析プログラムを使用できる。
【００２７】
（ａ）下部プレナム領域全体を流動解析対象とした下部プレナム流動解析プログラム
（ｂ）燃料集合体４バンドルの炉心入口部分を流動解析対象とした炉心入口オリフィス流動解析プログラム
【００２８】
下部プレナム流動解析プログラムは、ダウンカマの下端から炉心入口までの下部プレナム領域を多数の計算格子に分割し、ナビエ−ストークス方程式を解くことにより流速場を計算するものである。計算格子の大きさは、制御棒駆動機構のガイドチューブからなる管群の間の冷却材流動に基づいている。また、下部プレナム流動解析プログラムは、下部プレナムを３６０°にわたって計算格子化したものである。下部プレナム流動解析プログラムは１０機の再循環ポンプ一台ずつに対応する１０箇所の流入境界を有している。下部プレナム流動解析プログラムは、炉心入口のビームが入っている領域から燃料集合体領域までを、圧損を与えるための計算格子として有している。下部プレナム流動解析プログラムは燃料集合体の圧損を与える計算格子の上端に圧力境界を有している。
【００２９】
炉心入口オリフィス流動解析プログラムは、炉心入口のビームによる流動への影響が無視できる程度に離れた下部プレナム領域から、炉心入口オリフィスを通って燃料チャンネル下端までを流動解析の対象としている。炉心入口オリフィス流動解析プログラムは、オリフィスを通る冷却材の流動に基づいて計算格子の大きさを決定している。
【００３０】
炉心入口オリフィス流動解析プログラムは、下部プレナム領域の下端に流入境界を有している。炉心入口オリフィス流動解析プログラムは、燃料チャンネル下端から上部プレナム入口までの圧力損失を与えるための計算格子を有している。炉心入口オリフィス流動解析プログラムは、燃料チャンネルでの圧力損失を与えるための計算格子の上端に圧力境界を有している。
【００３１】
大規模流動解析システムは、局所的な流動解析を実施するための解析プログラムの他に、プラント全体の挙動を計算するためのプラント全体解析プログラムを有している。プラント全体解析プログラムは、プラントを構成している流体の流路をノードに分割し、ノード間をジャンクションで結合したものである。ノード内では流体の圧力、温度、密度、組成は均一である。ノード内で圧力、温度、密度、組成の変化量を計算し、ジャンクションでノード間の質量及びエネルギーの移動量を計算する。更にプラント全体解析プログラムは、プラント制御系を模擬したロジックを備えている。
【００３２】
大規模流動解析システムは、プラント全体解析プログラムと下部プレナム流動解析プログラムと炉心入口オリフィス流動解析プログラムとの間のデータの変換及び受け渡しを行なうインターフェースプログラムを有している。
【００３３】
図３に、下部プレナム流動解析プログラム境界条件とプラント全体解析プログラム出力及びインターフェースプログラムでのデータ加工の関係を示す。また、図４に、プラント全体解析プログラムと下部プレナム流動解析プログラムとの間でインターフェースプログラムにおけるノードとジャンクションを示す。
【００３４】
インターフェースプログラムは、ユーザーと大規模流動解析システムとの間のインターフェースの機能も有する。図５〜９は本実施の形態の大規模流動解析システムを利用して炉心への流量配分を評価する場合の計算フロー図である。また、図１０は、この実施の形態におけるデータの相互関係を示す。以下、図５〜９に従って説明を行なう。
【００３５】
インターフェースプログラムが大規模流動解析システムの入力画面１を端末のディスプレイ装置に表示する（ｉ１）。ユーザーは、インターフェースプログラムが表示する入力画面上で、解析対象とするプラントと流動解析を実施する領域を指定する（ｕ１）。本実施の形態では、炉心への流量配分を評価する場合を例題としているため、流動解析を実施する領域は、下部プレナム領域である。インターフェースプログラムは、ユーザーが解析対象として指定したプラントに対応したプラント全体解析プログラムＰ１をデータベースＤ１から検索する（ｉ２）。
【００３６】
インターフェースプログラムは、検索された全体解析プログラムＰ１の計算実行に必要なファイル一式をデータベースＤ１からユーザーの作業ディレクトリにコピーする（ｉ３）。また、インターフェースプログラムは、ユーザーが指定した流動解析を実施する領域に対応した流動解析プログラムをデータベースＤ１から検索する（ｉ４）。本実施の形態では、インターフェースプログラムは、下部プレナム流動解析プログラムＬＰ１を検索する。インターフェースプログラムは、検索された下部プレナム流動解析プログラムＬＰ１の計算実行に必要なファイル一式をデータベースＤ１からユーザーの作業ディレクトリにコピーする（ｉ５）。
【００３７】
ユーザーは、インターフェースプログラムが表示する入力画面上で、プラントの初期運転状態（炉心出力、再循環ポンプ流量）と過渡事象、チェックポイントの間隔を指定する（ｕ２）。ここで、チェックポイントとは、Ｐ１からＬＰ１に境界条件を受け渡す時刻であると共に、Ｐ１とＬＰ１の計算結果を比較し、Ｐ１の係数調整を行なうかどうかの判断を行なう時刻である。
【００３８】
ユーザーはあるチェックポイントから次のチェックポイントまでの時間間隔を指定する。この間隔は、一定値である場合と、時間の経過と共に変化させる場合とがある。例えば、解析対象となる時間全てにわたって、１０秒間隔のチェックポイントを指定する場合、あるいは、初期の１０秒間は１秒間隔のチェックポイントを設定し、次の１００秒間は１０秒間隔のチェックポイントを設定する場合である。
【００３９】
ユーザーは境界条件の時間変化が激しいと予想される時間帯については、チェックポイント間隔を小さく、境界条件の時間変化が緩やかな時間帯については、チェックポイント間隔を大きくとる。インターフェースプログラムは、ユーザーによって指定された初期運転状態と過渡事象及びチェックポイントの間隔に基づいて、プラント全体解析プログラムＰ１の入力ファイルＰ１＿ＩＮＰ１をユーザーの作業ディレクトリ上に作成する（ｉ６）。
【００４０】
ユーザーは、インターフェースプログラムが表示する入力画面上で、解析実行を指示する（ｕ３）。インターフェースプログラムは、ユーザーによる解析実行の指示（ｕ３）を認識して、プラント全体解析プログラムＰ１を起動する（ｉ７）。プラント全体解析プログラムＰ１は、起動後（ｐ１）、時刻ｔ_{ｐｉ− １}をタイムステップ時間Δｔｐ分だけ先に進め、ｔ_ｐｉとする（ｐ２）。時刻ｔ_ｐｉにおける各ノードの流体の圧力、温度、組成、密度と、各ジャンクションにおける流体の質量移動速度を計算する（ｐ３）。
【００４１】
プラント全体解析プログラムＰ１は、時刻ｔ_ｐｉとチェックポイントＣＰとを比較する（ｐ４）。時刻ｔ_ｐｉがチェックポイントＣＰ以上である場合、プラント全体解析プログラムＰ１は、前回のチェックポイント（ＣＰ−ΔＣＰ）の時刻から現在のチェックポイントＣＰまでの各ノードの流体の圧力、温度、組成、密度と、各ジャンクションにおける流体の質量移動速度を時系列データＰ１＿ＯＵＴｎとして出力する（ｐ５）。時刻ｔ_ｐｉがチェックポイントＣＰ未満である場合、プラント全体解析プログラムＰ１は、次のタイムステップの計算を行なう。
【００４２】
ｔ_ｐｉが計算終了時刻以上である場合（ｐ６）、プラント全体解析プログラムＰ１は、計算を終了する（ｐ７）。プラント全体解析プログラムＰ１が計算を終了すると、プロセスＩＤ（ＩＤ１）が消滅したことを、インターフェースプログラムが認識する（ｉ９）。ｔ_ｐｉが計算終了時刻未満である場合（ｐ６）、プラント全体解析プログラムＰ１は、チェックポイントＣＰの値を更新する（ｐ８）。
【００４３】
インターフェースプログラムはプラント全体解析プログラムが出力した時系列データＰ１＿ＯＵＴｎを用いて、図３に示す演算を行ない、下部プレナム流動解析プログラムの境界条件を算出する（ｉ８）。インターフェースプログラムは、下部プレナム流動解析プログラムの境界条件を、予め定められた書式に従って、ユーザーの作業ディレクトリ上のファイルＬＰ１＿ＢＤｎに書き出す（ｉ１０）。このとき、ＬＰ１＿ＢＤｎにはチェックポイントの値が含まれている。ＬＰ１＿ＢＤｎを作成する基になったファイルＰ１＿ＯＵＴｎは２つのチェックポイント（ＣＰ−ΔＣＰ）とＣＰの間の時系列データである。
【００４４】
インターフェースプログラムは、Ｐ１＿ＯＵＴｎの最初の時刻の値、すなわちＣＰ−ΔＣＰをＴｉｍｅ１としてＬＰ１＿ＢＤｎに書き出す（ｉ１０）。また、インターフェースプログラムは、Ｐ１＿ＯＵＴｎの最後の時刻の値、すなわちＣＰをＴｉｍｅ２としてＬＰ１＿ＢＤｎに書き出す（ｉ１０）。下部プレナム流動解析プログラムの境界条件のファイルを作成した後、インターフェースプログラムは、下部プレナム流動解析プログラムＬＰ１を起動する（ｉ１１）。
【００４５】
下部プレナム流動解析プログラムＬＰ１は、起動後（ｌｐ１）、時刻ｔ_{ｌｐｉ− １}をタイムステップΔ時刻ｔ_ｌｐ分だけ先に進め、ｔ_ｌｐｉとする（ｌｐ２）。下部プレナム流動解析プログラムＬＰ１は、時刻ｔ_ｐｉにおける下部プレナム領域の計算格子の圧力、温度、密度、計算格子間の流速を計算する（ｌｐ３）。
【００４６】
下部プレナム流動解析プログラムＬＰ１は、時刻ｔ_ｌｐｉとチェックポイントＣＰとを比較する（ｌｐ４）。時刻ｔ_ｌｐｉがチェックポイントＣＰ以上である場合、下部プレナム流動解析プログラムＬＰ１は、前回のチェックポイント（ＣＰ−ΔＣＰ）の時刻から現在のチェックポイントＣＰまでの下部プレナム領域の計算格子の圧力、温度、密度、計算格子間の流速を時系列データＬＰ１＿ＯＵＴｎとして出力する（ｌｐ５）。
【００４７】
時刻ｔ_ｌｐｉがチェックポイントＣＰ未満である場合、下部プレナム流動解析プログラムＬＰ１は、次のタイムステップの計算を行なう。ｔ_ｌｐｉが計算終了時刻以上である場合（ｌｐ６）、下部プレナム流動解析プログラムＬＰ１は、計算を終了する（ｌｐ７）。下部プレナム流動解析プログラムＬＰ１が計算を終了すると、プロセスＩＤ（ＩＤ２）が消滅したことを、インターフェースプログラムが認識する（ｉ１２）。
【００４８】
ｔ_ｌｐｉが計算終了時刻未満である場合（ｌｐ６）、下部プレナム流動解析プログラムＬＰ１は、インターフェースプログラムから計算中止命令（ｉ２７， ｉ３７）または計算待機命令（ｉ３４）が入らなければ、境界条件ファイルＬＰ１＿ＢＤｎの値が更新されるのを待つ（ｌｐ１３）。下部プレナム流動解析プログラムＬＰ１は境界条件ファイルＬＰ１＿ＢＤｎを読み込む。ファイルが更新されていなければ、ＣＰ＝Ｔｉｍｅ２になっている。
【００４９】
プラント全体解析プログラムＰ１が次のチェックポイントまでの出力ファイルを作成し、インターフェースプログラムが境界条件ファイルＬＰ１＿ＢＤｎを更新すると、更新前のＴｉｍｅ２の値が、更新後のＴｉｍｅ１の値になっている。
【００５０】
下部プレナム流動解析プログラムＬＰ１では、チェックポイントＣＰは旧いままであるため、ファイルが更新されるとＣＰはＴｉｍｅ１の値と等しくなる。こうすることで、下部プレナム流動解析プログラムＬＰ１は境界条件ファイルＬＰ１＿ＢＤｎが更新されたことを認識する。境界条件ファイルＬＰ１＿ＢＤｎの値が更新されたら、下部プレナム流動解析プログラムＬＰ１は、チェックポイントＣＰの値を更新し（ｌｐ１４）、次のタイムステップの計算を行なう。
【００５１】
インターフェースプログラムは、プラント全体解析プログラムＰ１のプロセスＩＤ（ＩＤ１）及び下部プレナム流動解析プログラムＬＰ１のプロセスＩＤ（ＩＤ２）の消滅を検知すると（ｉ９， ｉ１２）、計算終了のメッセージ３を端末のディスプレイ装置に表示する（ｉ１４）。プラント全体解析プログラムＰ１及び下部プレナム流動解析プログラムＬＰ１の係数が変更されていない場合、インターフェースプログラムは解析結果表示画面を表示する（ｉ２０）。
【００５２】
プラント全体解析プログラムＰ１の係数が変更された場合（ｉ１５）に、インターフェースプログラムは、プラント全体解析プログラムＰ１のバージョン名を変更する（ｉ１６）。本実施の形態においては、下部プレナム領域の圧力損失係数あるいは炉心入口オリフィス圧力損失係数が変更された場合に、インターフェースプログラムは、プラント全体解析プログラムＰ１のバージョン名を変更する（ｉ１６）。バージョン名とは、プラント全体解析プログラムＰ１の変更履歴を表したものであり、プラント全体解析プログラムＰ１の実行に必要なファイル一式を格納したディレクトリ名の一部である。
【００５３】
インターフェースソフトはバージョン名を変更したプラント全体解析プログラムＰ１をデータベースＤ１に保存する（ｉ１７）。更にインターフェースプログラムは、変更された係数の名前と変更前の値及び変更後の値、変更した年月日、誤差の比較対照となった下部プレナム流動解析プログラムのバージョン名をログファイルに書き出す（ｉ１８）。インターフェースプログラムは変更内容を書き出したファイルをデータベースＤ１に保存する（ｉ１９）。その後、インターフェースプログラムは解析結果表示画面をディスプレイ装置に表示する（ｉ２０）。
【００５４】
下部プレナム流動解析プログラムＬＰ１の係数が変更された場合（ｉ１５）に、インターフェースプログラムは、下部プレナム流動解析プログラムＬＰ１のバージョン名を変更する（ｉ１６）。本実施の形態においては、炉心入口オリフィス圧力損失係数が変更された場合に、インターフェースプログラムは、下部プレナム流動解析プログラムＬＰ１のバージョン名を変更する（ｉ１６）。後の処理はプラント全体解析プログラムＰ１の場合と同様である。
【００５５】
下部プレナム流動解析プログラムＬＰ１から出力ファイルＬＰ１＿ＯＵＴｎが出力されたら、インターフェースプログラムは、ＬＰ１＿ＯＵＴｎの値を読み込む（ｉ２１）。インターフェースプログラムは、下部プレナム流動解析プログラムＬＰ１が出力した炉心の圧損を模擬する計算格子における冷却材の鉛直方向質量流量［ｋｇ／ｓ・ｍ^２］を、燃料チャンネル一体に相当する領域に渡って積算する （ｉ２１）。
【００５６】
インターフェースプログラムは、プラント全体解析プログラムＰ１で計算された各燃料チャンネルでの冷却材流量と下部プレナム流動解析プログラムＬＰ１で計算された各燃料チャンネルでの冷却材流量との誤差を計算する （ｉ２２）。
【００５７】
インターフェースプログラムは上記冷却材流量の誤差と、ユーザーが予め設定している誤差の許容値とを比較する（ｉ２３）。誤差が許容値以内であるとインターフェースプログラムが判断した場合には、インターフェースプログラムは、炉心入口オリフィス部分の圧力損失係数のチェックを行なう（後述、図８および９のａ４以降）。この時、ｉ２３で誤差が許容値以内であると判断された場合の直後の処理として、後述する係数調整フラグをｏｆｆする処理を加えても良い。このようにすると誤差が許容値以内にある場合には係数調整を行なわず、Ｐ１とＰＬ１の解析をそれぞれ続行させることができる。
【００５８】
誤差が許容値外であるとインターフェースプログラムが判断した場合、インターフェースプログラムは、係数調整フラグがＯＮになっているかどうか判断する（ｉ２４）。ここで、係数調整フラグは、通常はＯＦＦに設定されており、ユーザーが係数調整を指示した場合にのみＯＮになる変数である。係数調整フラグがＯＦＦの場合、インターフェースプログラムは、端末のディスプレイ装置にメッセージ１とともに許容値外の誤差が生じたチャンネルの番号と誤差の値、係数の調整を実行するためのボタンを表示する（ｉ２５）。
【００５９】
インターフェースプログラムが表示する入力画面上で、上記メッセージをユーザーが確認（ｕ４）後、ユーザーが解析中断の指示（ｕ５）を行なった場合には、インターフェースプログラムはプラント全体解析プログラムＰ１及び下部プレナム流動解析プログラムＬＰ１の実行を中断する（ｉ２６）。
【００６０】
ここで、ｉ２４のステップの直前に係数フラグがＯｎになってからのｉ２４のステップに至る回数を数え、予め設定した回数より多い場合に係数フラグをＯｆｆする処理を加えても良い。このようにすると、予め設定した係数調整を行なっても誤差が許容値以内にならない場合にもＰ１およびＬＰ１の解析を続行させることができる。
【００６１】
例えばＰ１とＰＬ１の誤差を比較する対象である燃料チャンネルの本数は、ＡＢＷＲタイプのプラントの場合、８７２体もある。これら全ての燃料チャンネルの流量を許容値以内の誤差に収めようとして、係数調整を行なったとしても、収束しない場合がある。このようなときに、ユーザーは、どの燃料チャンネルでどの程度の誤差が生じているのか把握し、全燃料チャンネルのうち一部がわずかに許容値を満たさない場合であれば、解析を続行させることができる。
【００６２】
インターフェースプログラムが表示する入力画面上で、ユーザーが係数の調整を実行するためのボタンを押して、プラント全体解析プログラムＰ１の係数調整を指示（ｕ６）した場合、インターフェースプログラムは係数調整フラグをＯＮにして（ｉ２９）、係数調整画面２を表示する（ｉ３０）。ユーザーは係数調整画面２で調整する係数を指定する（ｕ７）。なお、本実施の形態においては、ユーザーが調整を指定できる係数は、下部プレナムを模擬しているノード間の圧力損失係数である。
【００６３】
インターフェースプログラムは、ユーザーが指定した係数を元の値から変化させた入力ファイルＰ１＿ＩＮＰｎを作成し（ｉ３１）、プラント全体解析プログラムＰ１を起動する（ｉ７）。中断前と同様の手順により、インターフェースプログラムは下部プレナム流動解析プログラムＬＰ１の境界条件のファイルＬＰ＿ＢＤｎを作成し、下部プレナム流動解析プログラムＬＰ１を実行する。
【００６４】
プラント全体解析プログラムＰ１による燃料チャンネル流量の計算結果Ｐ＿ＯＵＴｎと下部プレナム流動解析プログラムＬＰ１による燃料チャンネル流量の計算結果ＬＰ＿ＯＵＴｎとの誤差と許容値とを比較する。誤差が許容値外であり、係数調整フラグがＯＮであれば、インターフェースプログラムはプラント全体解析プログラムＰ１と下部プレナム流動解析プログラムＬＰ１の実行を中断する（ｉ２７）。
【００６５】
誤差が許容値内に収まるまで、インターフェースプログラムは、調整対象になっている係数を変化させて、プラント全体解析プログラムＰ１の入力ファイルＰ１＿ＩＮＰｎを作成し（ｉ３１）、プラント全体解析プログラムＰ１を起動し、下部プレナム流動解析プログラムＬＰ１の境界条件ファイルＬＰ＿ＢＤｎを作成し、下部プレナム流動解析プログラムＬＰ１を起動して、解析結果の誤差を求め、許容値と比較するという作業を繰り返す。調整対象となっている係数を変化させる方法は、一般に最適設計手法で用いられている遺伝的アルゴリズムや最急降下法などを用いる。
【００６６】
プラント全体解析プログラムＰ１と下部プレナム流動解析プログラムＬＰ１との各燃料チャンネルでの冷却材流量の計算結果が、許容値以内に収まるか、許容値外でもユーザーが解析中断を指定しない場合は、インターフェースプログラムは炉心入口オリフィスの圧力損失係数の適用範囲のチェックを行なう（ａ４以降）。
【００６７】
インターフェースプログラムは、炉心入口オリフィスを模擬しているノードの冷却材圧力、流速をプラント全体解析プログラムの出力ファイルＰ１＿ＯＵＴｎから取得する（ｉ３２） 。また、インターフェースプログラムは、炉心入口オリフィスを模擬している計算格子の冷却材圧力、流速を下部プレナム流動解析プログラムの出力ファイルＬＰ１＿ＯＵＴｎから取得する（ｉ３２）。
【００６８】
インターフェースプログラムは、炉心入口オリフィスを模擬しているノードの冷却材圧力、流速及びその変化率を炉心入口オリフィスの圧力損失係数の適用範囲と比較する（ｉ３３）。炉心入口オリフィスを模擬しているノード及び計算格子の冷却材圧力、流速及びその変化率のいずれも、炉心入口オリフィスの圧力損失係数の適用範囲を逸脱していない場合、インターフェースプログラムは、次のチェックポイントまでプラント全体解析プログラムＰ１と下部プレナム流動解析プログラムＬＰ１の実行を継続する。
【００６９】
炉心入口オリフィスを模擬しているノードまたは計算格子の冷却材圧力、流速及びその変化率のいずれかが、炉心入口オリフィスの圧力損失係数の適用範囲を逸脱した場合、インターフェースプログラムはプラント全体解析プログラムＰ１及び下部プレナム流動解析プログラムＬＰ１に待機命令を出す（ｉ３４）。インターフェースプログラムは、端末のディスプレイ装置にメッセージ２とともに適用範囲外の値が生じたオリフィスのタイプと適用範囲を逸脱した変数の名前とその値、係数の再計算を実行するためのボタンを表示する（ｉ３５）。
【００７０】
インターフェースプログラムが表示する入力画面上で、上記メッセージ２をユーザーが確認後（ｕ８）、ユーザーが解析中断の指示を行なった場合（ｕ９）には、インターフェースプログラムはプラント全体解析プログラムＰ１及び下部プレナム流動解析プログラムＬＰ１の実行を中断する（ｉ３７）。
【００７１】
インターフェースプログラムが表示する入力画面上で、ユーザーが係数の再計算を実行するためのボタンを押して、炉心入口オリフィスの圧力損失係数の再計算を指示した場合（ｕ１０）、インターフェースプログラムは係数再計算画面３を表示する（ｉ３８）。ユーザーは係数再計算画面３で計算対象を指定する（ｕ１１）。なお、この場面において、ユーザーが指定する項目は、炉心入口オリフィスのタイプ、境界条件の範囲（流速、圧力）、計算の種類（定常か非定常か）である。
【００７２】
インターフェースプログラムは、ユーザーが指定した炉心入口オリフィスのタイプに対応した炉心入口オリフィス流動解析プログラムＥＯ１をデータベースＤ１から検索する（ｉ３９）。インターフェースプログラムは、検索された炉心入口オリフィス流動解析プログラムＥＯ１の計算実行に必要なファイル一式をデータベースＤ１からユーザーの作業ディレクトリにコピーする（ｉ４０）。
【００７３】
インターフェースプログラムは、係数再計算画面３でユーザーが指定した境界条件の範囲及び計算の種類に基づいて、炉心入口オリフィス流動解析プログラムＥＯ１の入力ファイルＥＯ＿ＩＮＰ１及び境界条件ファイルＥＯ＿ＢＤ１をユーザーの作業ディレクトリ上に作成する（ｉ４１）。インターフェースプログラムは炉心入口オリフィス流動解析プログラムＥＯ１を起動する（ｉ４２）。
【００７４】
炉心入口オリフィス流動解析プログラムＥＯ１の計算が終了し、出力ファイルＥＯ１＿ＯＵＴ１が出力されると、インターフェースプログラムは炉心入口オリフィスの出力ファイルＥＯ＿ＯＵＴ１中に格納されている以下に示すデータを用いて、圧力損失係数を計算する（ｉ４３）。
【００７５】
（１）炉心入口のビームより下方の計算格子における冷却材の圧力と鉛直方向速度
（２）炉心入口オリフィス部分の計算格子より上方でかつ炉心の圧損を模擬している計算格子より下方の計算格子における冷却材の圧力
【００７６】
インターフェースプログラムは、データベースＤ１に保存されている同じ炉心入口オリフィスのタイプについて実施された既存の計算結果に、今回計算した結果を追加する（ｉ４４）。インターフェースプログラムはプラント全体解析プログラムＰ１の炉心入口オリフィスの圧力損失係数を修正し、入力ファイルＰ１＿ＩＮＰｎ＋１をユーザーの作業ディレクトリに作成する（ｉ４５）。インターフェースプログラムは、プラント全体解析プログラムＰ１を起動する（ｉ７）。後の手順は、炉心入口オリフィスの圧力損失係数が適用範囲を逸脱していない場合と同様である。
【００７７】
以上のように、本実施の形態の大規模流動解析システムを用いて、炉心への流量配分を計算した場合、プラント全体の冷却材圧力変化、再循環ポンプの運転状態、炉心の出力等はプラント全体解析プログラムが計算する。このプラント全体解析プログラムの結果を境界条件として、下部プレナム流動解析プログラムで炉心への流量配分を計算する。逆に炉心への流量配分が下部プレナム流動解析プログラムの結果と同等になるよう、プラント全体解析プログラムの圧力損失係数を調整する。調整された係数及び解析プログラムと入力データは全てデータベースに保存する。解析結果が蓄積されるに従って、プラント全体解析プログラムの精度が向上していく。
【００７８】
また、プラント全体解析プログラムで使用している係数や相関式の適用範囲外の条件が生じた場合（この実施の形態では炉心入口オリフィスの圧力損失係数を取り上げたが）、その部分だけ詳細に解析するプログラムを起動して、係数あるいは相関式の設定を行なうことができる。このように、本実施の形態の大規模流動解析システムは、詳細な流動解析の結果を反映したプラント全体解析を行なう。
【００７９】
以上説明したように、第１の実施の形態においては、着目する現象あるいはパラメーターが顕在化するスケール及び領域にフォーカスし、各スケール及び領域について、それらを許容可能な精度で評価できる解析モデルを用いる。例えば、システム挙動に着目した場合、重要となるのは系内の状態量のマクロな分布とその変化率であり、こうした現象を評価するためには、ノード・ジャンクション法による解析モデルが適している。
【００８０】
また、狭隘複雑流路での圧力損失に着目した場合、重要となるのは壁面摩擦効果及び渦流の発生などの局所的な流れの状態であり、こうした現象を評価するためには、詳細な解析格子を有した数値流動解析モデルが適している。
【００８１】
こうして個々の現象については、その現象が顕在化するスケール及び領域にフォーカスしたモデルで解析が行なわれ、その現象に直接的に関与しないスケール及び領域については、境界条件あるいはモデルパラメーターあるいは物理モデルあるいは相関式モデルとして取り扱うことで、解析精度を損なうことなく、解析時間及び計算機メモリ容量及び演算処理装置数の負荷低減を図ることができる。
【００８２】
ここで重要なことは、個々の解析が独立して存在しているのではなく、あらゆるスケール及びあらゆる領域の解析同士が互に境界条件あるいはモデルパラメーターあるいは物理モデルあるいは相関式モデルを介して相互に関わり合っている点である。こうした構成にすることによって、個別機器の設計がシステム全体に及ぼす影響を把握し、なおかつ、システム全体挙動の一部分として個別機器の境界条件を適切に与えることのできる実用的な大規模流動解析システムを得ることができる。
なお、例えば単相流領域と二相流領域を分割した解析を実行することもできる。
【００８３】
更に、数値流体解析モデルから流体力及び状態量に関する境界条件を与えられ、流体場に置かれた構造物の応力と変位を数値構造解析プログラムを用いて計算する構造解析モデルを用いることにより、原子炉システムのある事象、特に過渡事象における流体側の熱水力的な挙動を、構造解析への境界条件として与えることが可能となる。
【００８４】
図１１は複数の同一形状機器を有する流体系システムに本発明の大規模流動解析システムを適用した場合の、第２の実施の形態の説明図である。機器単体の性能Ｐは機器形状ｇと環境条件ｅ、境界条件ｂによって次のように決定できる。
Ｐ＝ｆ（ｇ， ｅ， ｂ） …（１）
【００８５】
この関係を単体流動解析モデル１０４を用いて作成する。単体流動解析モデル１０４では、数値流動解析手法を用いて機器形状ｇにおける機器内部の詳細流れを解くことにより、流動状態に影響を受ける機器の特性及び性能を種々の環境条件ｅ及び境界条件ｂの下で評価する。次にこれらの結果から機器単体性能Ｐと機器形状ｇ、環境条件ｅ、境界条件ｂの関係（１）を作成する。
【００８６】
次に（１）の単体性能を持つ機器が複数集まって構成されるコンポーネントに対し、機器の集合体流動解析モデル１０５を用いて、機器の集合体としての挙動を評価すると共に、集合体内の個々の機器への境界条件ｂ及び環境条件ｅを評価する。ここで、集合体流動解析モデルにおける単体機器性能評価に関係（１）を用いることで、解析負荷の低減を図る。
【００８７】
また、機器集合体流動解析モデル自体の環境条件Ｅ及び境界条件Ｂは、流体系全体解析１０６から与えられ、一方、関係（１）は流体系全体解析１０６における機器単体性能評価モデルとなる。更に、流体系全体解析で計算した機器単体への境界条件Ｂ’及び環境条件Ｅ’が、機器集合体流動解析モデル１０５を用いて解析した境界条件ｂ及びｅと同等になるよう、流体系全体解析の物理パラメーター、物理モデル、相関式モデルの修正または追加を行なう。
【００８８】
こうすることで、流体系全体解析において、機器の形状に依存した流動状態を考慮した性能評価と、機器集合体の挙動を考慮した境界条件及び環境条件の評価が可能になるとともに、機器単体の性能評価においても、システム挙動及び集合体としての挙動を考慮した境界条件及び環境条件の取得が可能となる。
【００８９】
以上説明したように、第２の実施の形態においては、特に原子炉系のような多数の同一形状機器の集合体を有する流体系において、個々の機器内部の詳細流動を数値流動解析によって解くとともに、その結果に基づいて作成された物理モデルまたは相関式モデルまたは物理パラメーターに基づいた解析負荷の小さいモデルで機器の集合体としての挙動を解析することにより、解析精度を損なうことなく、解析時間を大幅に短縮することが可能となる。
【００９０】
図１２は、本発明の大規模流動解析システムを、流体系全体と個別機器の相互作用を考えた最適設計に適用した場合の第３の実施の形態の説明図である。以下の最適設計への適用に関する実施の形態の説明では、流体系全体のシミュレーション結果で評価可能な炉心出力の安定性および最小限界出力比が設計基準を満足するように炉心入口部における圧損係数に関わる主要因子である炉心入口オリフィスの口径を最適化する場合を例に記述する。
【００９１】
なお、炉心出力の安定性および最小限界出力比は、炉心入口部における圧損係数が大きくなると炉心出力の安定性は向上するものの、最小限界出力比は低下し熱的に厳しくなるというトレードオフの関係にある。
【００９２】
プラント全体解析プログラムによる解析で得られる原子炉出力、原子炉圧力、原子炉水位、炉心流量などの挙動に対して解析者が設定した、例えば炉心出力の安定性の設計基準や炉心燃料の最小限界出力比の設計基準といった目的関数Ｈ１（Ｐ，Ｑ，Ｒ，・・・）、Ｈ２（Ｐ，Ｓ，Ｔ・・・）に対して、その独立変数の一つであるＰを最適化することにより、目的関数Ｈ１及びＨ２を満足する場合を考える。なお、ここでは説明のために、独立変数Ｐを機器単体１０７の性能、例えば炉心入口に複数設置されている炉心入口オリフィスの圧損係数と仮定する。
【００９３】
目的関数Ｈ１及びＨ２を満足するために、機器単体１０７の形状を見直す場合、形状変更に伴う機器単体内部での詳細流動情報については、全体システム１０９を対象とした流体系全体解析のレベルで直接評価することができない。そこで、目的関数Ｈ１及びＨ２を満足するように選ばれた流体系全体解析における独立変数Ｐを、機器単体１０７に関する単体流動解析モデルの目的関数Ｐ（ｇ，ｅ，ｂ）として新たに解析者が設定し、目的関数Ｐ（ｇ，ｅ，ｂ）が全体システムレベルでの要求をすなわち目的関数Ｈ１及びＨ２を満足するよう、その独立変数である機器形状因子ｇの最適化を行なう。 ここでの例の場合、機器形状因子ｇは炉心入口オリフィスの口径が挙げられる。機器集合体１０８全体を解析対象として、その領域内の流体挙動、および集合体内の個々の機器の環境条件及び境界条件を評価するための集合体流動解析においては、前記単体流動解析で得られる目的関数Ｐ（ｇ，ｅ，ｂ）の値すなわち炉心入口オリフィス圧損係数は、モデルパラメーターＰ（ｇ，ｅ，ｂ）すなわち集合体流動解析モデルにおける炉心入口オリフィス部の圧力損失を別途計算して付加するための定数として設定される。
【００９４】
集合体流動解析と単体流動解析及び流体系全体解析間の解析結果のやり取りは第２の実施の形態と同様である。
このような構成によれば、例えば、炉心入口オリフィス口径をパラメーターに単体流動解析でその圧損係数を求め、それを集合体解析、流体系全体解析モデルにモデルパラメーターとしてフィードバックしながら、単体流動解析、集合体解析、およびプラント全体解析を連成させて解析する。
【００９５】
これにより、プラント全体解析で評価可能な炉心出力の安定性の減幅比や炉心燃料の最小限界出力比といった目的関数Ｈ１（Ｐ，Ｑ，Ｒ，・・・）、Ｈ２（Ｐ，Ｓ，Ｔ・・・）の値が設計基準を満足する最適なオリフィス口径を求めることが可能となる。このように本発明の実施の形態の大規模流動解析システムを用いることにより、流体系全体と個別機器の相互作用を考えた最適化問題を解析的に解くことで、例えば流体系全体の流体挙動を最適化する機器レベルの形状を選定するなどの最適設計が可能になる。
【００９６】
更に、以上説明した第１〜第３の実施の形態によれば、多数のユーザーが多数の解析ケースを実施し解析の入出力データが蓄積されれば、有効利用可能なリソースが増えるために大規模な流動解析を効率よく実行することが可能となり、かつ、流動スケールに応じて下位は上位の部分領域となるように複数の解析モデルによって階層化された構成とし、下位の解析モデルの解析結果に基づいて上位の解析モデルに内在する物理モデル、相関式モデル、および物理パラメーターを更新することが可能なシステムであるために、利用可能な解析結果が増えるほど精度が高い解析モデルに更新される。
【００９７】
次に、第４の実施の形態として、図１３および図１４を用いて解析の入力データやそれをもとに解析した結果である出力データをデータサーバーにおいて集中管理するディレクトリ構造を適用した大規模流動解析システムについて説明する。ここで、本実施の形態の大規模流動解析システムを運用する計算機環境の例を図１４に示す。
【００９８】
上記第１〜第３の実施の形態として説明した原子炉系解析システムにおいては、多数の解析者が多数の解析ケースを実施することによって解析の入出力データが蓄積されれば、有効利用可能なリソースが増えるために大規模な流動解析を効率よく実行することが可能となる。しかも、流動スケールに応じて下位は上位の部分領域となるように複数の解析モデルによって階層化された構成として、下位の解析モデルの解析結果に基づいて上位の解析モデルに内在する物理モデル、相関式モデル、および物理パラメーターを更新することが可能なシステムであるために、利用可能な解析結果が増えるほど精度が高い解析モデルに更新される。
【００９９】
そこで、本実施の形態は、例えば図１４に示した計算機環境において、データサーバーで解析の入力データやそれをもとに解析した結果である出力データを集中管理するディレクトリ構造を、共有データベース領域３１１と各クライアントの作業領域３１２とに分ける。
【０１００】
解析者は自分に割り当てられたクライアント作業領域では読み書き自在にデータにアクセス可能で、共有データベース領域および他の解析者に割り当てられたクライアント作業領域には読み込みのみ許可され、かつ、管理者が認めた場合には、各クライアント領域にある解析の入力データやそれをもとに解析した結果である出力データを共有データベース領域に登録可能な状態に管理する。
【０１０１】
共有データベース領域３１１と各クライアントの作業領域３１２の両領域は、例えば原子炉の種類、プラントの初期運転状態、解析対象とするプラント事象、解析ケース、および解析モデルのように解析対象と解析条件の違いによって分類された同一の階層構造３１３とする。
【０１０２】
例えば、原子炉の種類Ｉ、初期運転状態１、事象αにおいて、解析ケース（１）３１４については解析が実施されデータが共有データベース領域３１１に登録されている場合に、クライアントＡが、同事象において、共有データベース領域３１１に登録されている解析とは異なる条件、すなわち新たな解析ケースとして、原子炉系全体モデル、グローバル領域解析モデル（ｉ）、およびローカル解析モデルａによる解析を実施したい場合、クライアントＡ作業領域３１５に新たな解析ケース（Ａ１）３１６のフォルダを作成する。
【０１０３】
そして、その下を共有データベース領域と同一のディレクトリ構成として、既存の原子炉系全体モデルの入力データ３１７の読込み３１８を行ない編集した同モデルの入力データ３１９、原子炉系全体モデルの解析結果３２０、グローバル領域解析モデル（ｉ）による解析結果３２１、およびローカル解析モデルによる解析結果３２２を保存する。
【０１０４】
ここで、グローバル領域解析モデルおよびローカル解析モデルとは第１の実施の形態として記載の流動スケールに応じて階層化された数値流体解析モデルの分類の一例であり、第１の実施の形態においては、下部プレナム流動解析モデルがグローバル領域解析モデルに、その一部を詳細な空間分解能でモデル化した炉心入口オリフィス流動解析モデルがローカル領域モデルに各々相当する。
【０１０５】
グローバル領域解析モデル（ｉ）による解析結果３２１、ないしはローカル解析モデルによる解析結果３２２が他の解析者にとって活用できるものとシステム管理者が判断した場合には、管理者権限において、クライアントＡ作業領域の解析ケース（Ａ１）以下にある原子炉系全体モデル、グローバル領域解析モデル（ｉ）、およびローカル領域解析モデルａの入力データファイル、解析結果ファイルをディレクトリ構造ごと共有データベース領域に登録３２３すなわち書き込みすることが可能な構成とする。
【０１０６】
他のクライアント作業領域、例えばクライアントＡ作業領域を利用する解析者にとってクライアントＢ作業領域３２４に活用可能な入力データもしくは解析結果がある場合には、そこからコピーしてそのままもしくは入力データに関しては編集するなどして活用することも可能に管理する。
【０１０７】
この方法によれば、各領域において、原子炉の種類、初期運転状態、プラント事象といった解析対象、およびその解析対象の解析に適用した原子炉系全体モデル、グローバル領域解析モデルやローカル領域解析モデルなどの数値流体解析モデルによって階層的に分類された同一のディレクトリ構成である。
【０１０８】
その結果、解析者は保存されている入力データファイルや解析結果ファイルのプロパティがディレクトリ名に基づいて同定可能となる。そして、管理者による共有データベース領域３１１への解析入出力データの登録、解析者による他の領域の解析入出力データを活用した解析実行を、ディレクトリ構造を辿って必要な部分をディレクトリごとコピーして編集、活用するといった画一的なプロセスで行なうことが可能となる。
【０１０９】
上記作用によって、共有データベース３１１への解析データの蓄積が促進され、かつ、これらの有効活用が促進されるために、大規模解析の作業効率の向上、および解析モデルに内在する物理モデル等の更新に基づく解析精度の向上を図ることが可能となる
図１４は本実施の形態の大規模流動解析システムを運用する計算機環境の例を示したものである。
解析者が大規模流動解析システムの端末として利用するクライアント計算機３０１と、大規模流動解析システムのサーバー計算機３０２、ストレージ機器３０３、前記プラント全体解析プログラム、数値流体解析プログラム、もしくは構造解析プログラムを有し実行する高速計算機３０５、およびネットワーク３０６から構成されている。例えば、サーバー計算機３０２はウェブサーバーとなっており、クライアント計算機３０１からは、ウェブブラウザを用いて大規模流動解析システムプログラムを利用可能に構成されている。
【０１１０】
図１５は、本発明の第５の実施の形態を示すもので、大規模流動解析システムを運用する計算機環境の例を示したものである。これは、図１４に示した計算機環境を基本構成として、分散配置された複数の高速計算機３０５ａ〜ｆによる解析実行環境を、十分に高速な通信速度を有するネットワーク３０６を介して結合したグリッドコンピューティング構成である。
【０１１１】
この実施の形態によれば、分散配置された複数の高速計算機３０５ａ〜ｆ上にある解析プログラムや高速計算機３０５ａ〜ｆ等ハードウェアのリソースを大規模解析システムを構成する解析モジュールもしくはハードウェア要素として、サーバー計算機３０２によって集中管理することが可能となり、その有効活用が容易になる。そのため、大規模解析の解析負荷を分散して解析時間の短縮を図ることができる。
【０１１２】
【発明の効果】
以上説明したように、本発明によれば、原子炉内のように小さなスケールと大きなスケールの流動現象が互いに影響し合う大規模な体系について、スケールに応じた階層構造で領域分けした複数の解析モデルで構築し、各解析モデル同士が互いに境界条件、モデルパラメーター、物理モデル、あるいは相関式を解して相互に連成させることによって、個別機器の設計がシステム全体に及ぼす影響を正確に把握し、なおかつ、システム全体挙動の中における個別機器レベルの流動現象の挙動を評価できる大規模流動解析システムを提供することができる。
【図面の簡単な説明】
【図１】本発明に係る流動解析システムの第１の実施の形態における情報の流れを示すブロック図。
【図２】図１の流動解析の対象となる原子炉における流動を示す模式図。
【図３】本発明に係る流動解析システムの第１の実施の形態における下部プレナム流動解析プログラム境界条件とプラント全体解析プログラム出力及びインターフェースプログラムでのデータ加工の関係を示す表。
【図４】本発明に係る流動解析システムの第１の実施の形態におけるプラント全体解析プログラムと下部プレナム流動解析プログラムとの間でインターフェースプログラムにおけるノードとジャンクションを示す模式図。（ａ）は下部プレナム流動解析プログラムを示す図、（ｂ）はプラント全体解析プログラムを示す図。
【図５】本発明に係る流動解析システムの第１の実施の形態における計算フロー図。
【図６】図５の計算フロー図のつづき。
【図７】図６の計算フロー図のつづき。
【図８】図７の計算フロー図のつづき。
【図９】図８の計算フロー図のつづき。
【図１０】本発明に係る流動解析システムの第１の実施の形態におけるデータの相関関係を示すブロック図。
【図１１】本発明に係る流動解析システムの第２の実施の形態における情報の流れを示すブロック図。
【図１２】本発明に係る流動解析システムの第３の実施の形態における情報の流れを示すブロック図。
【図１３】本発明に係る流動解析システムの第４の実施の形態を示すブロック図。
【図１４】本発明に係る流動解析システムの第４の実施の形態における情報の流れを示すブロック図。
【図１５】本発明に係る流動解析システムの第５の実施の形態を示すブロック図。
【図１６】従来の流動解析システムにおける情報の流れを示すブロック図。
【符号の説明】
１０１…原子炉系全体モデル、１０２…下部プレナム流動解析モデル（大領域流動解析モデル）、１０３…炉心入口オリフィス流動解析モデル（局所領域解析モデル）、１０４…単体流動解析モデル、１０５…集合体流動解析モデル、１０６…流体系全体解析モデル、１０７…機器単体、１０８…機器集合体、１０９…全体システム、３０１…計算機（クライアント機能）、３０４…サーバー機能（ウェブサーバー＆データサーバー）、３０５…高速計算機、３０６…ネットワーク、３１１…共有データベース領域、３１２…クライアント作業領域。

Claims (10)

1. 流体系全体を複数のノードに分割し、ノード間をジャンクションで接続し、該ノード内で状態量を計算し、該ジャンクションでノード間の移行量を計算する最上位の数値流動解析モデルとしての流体系全体解析モデルと、解析対象領域によって前記流体系全体を分割してかつ流動スケールに応じて階層化された複数の部分流動解析モデルとを用いて、流体系の流動を解析する流動解析システムにおいて、
   前記流体系全体解析モデルを用いて流体系全体のシステム挙動を解析する手段と、
   前記流体系全体解析モデルによる解析結果を前記複数の部分流動解析モデルの境界条件あるいはモデルパラメーターとして利用し、前記流動スケールに応じて階層化された部分流動解析モデルの上位のモデルによる解析結果をそれより下位の部分流動解析モデルの境界条件あるいはモデルパラメーターとして利用して前記複数の部分流動解析モデルによる解析を行なう手段と、
   下位の部分流動解析モデルを用いた解析結果に基づいて作成される物理モデルまたは相関式モデルまたは物理パラメーターを上位の部分流動解析モデル及び前記流体系全体解析モデルの一部として組み込む手段と、
   前記流体系全体解析モデルを用いた解析結果に基づいて作成される物理モデルまたは相関式モデルまたは物理パラメーターを前記流体系全体解析モデルの一部として組み込む手段と、を有することを特徴とする流動解析システム。
2. 請求項１に記載の流動解析システムにおいて、前記複数の部分流動解析モデルを適用する領域が、流動を記述する支配方程式の相違によって分割されていること、を特徴とする流動解析システム。
3. 請求項１または２に記載の流動解析システムにおいて、
   前記流体系全体は互いに共通の単体性能を有する複数の機器を含み、前記流動解析システムは、
   前記複数の機器の単体性能を評価するための単体流動解析モデルによる解析を行なう手段と、
   前記単体流動解析モデルによる解析の結果に基づいて作成された物理モデルまたは相関式モデルまたは物理パラメーターを、前記複数の機器の集合体としての性能を評価するための集合体解析モデルの一部として組み込む手段と、
   前記集合体解析モデルによって、前記複数機器の集合体における個々の機器の境界条件を求める手段と、を有することを特徴とする流動解析システム。
4. 請求項１ないし３のいずれかに記載の流動解析システムにおいて、
   前記部分流動解析モデルを用いた計算から流体力及び状態量に関する境界条件を与える手段と、
   前記流体系の流体場に置かれた構造物の応力及び変位を数値構造解析モデルを用いて計算する手段と、
   をさらに有することを特徴とする流動解析システム。
5. 請求項１ないし４のいずれかに記載の流動解析システムにおいて、
   前記複数の部分流動解析モデルの下位に、該部分流動解析モデルに領域として内包される局所領域の局所流動解析モデルがあって、
   前記流体系全体解析モデルもしくは前記部分流動解析モデルの解析領域に対する最適化問題の独立変数をｘ１、ｘ２、・・・ｘｍとし、目的関数をＦ（ｘ１、ｘ２、・・・ｘｍ）とし、独立変数ｘｍの実行可能領域をＳｍとし、それよりも下位の流動解析モデルの解析領域に対する最適化問題の独立変数をｙ１、ｙ２、・・・ｙｎとし、目的関数をＧ（ｙ１、ｙ２、・・・ｙｎ）とし、独立変数ｙｎの実行可能領域をＴｎとするとき、該解析モデルにおいて、任意のｘｍおよびｙｎは解析パラメータとして設定自在に構成され、かつｘｍとｙｎの関係を解析によって算出可能にモデル化することによって、その関係に基づいてＳｍとＴｎのいずれか一方から他方を求め、目的関数Ｆ（ｘ１、ｘ２、・・・ｘｍ）もしくは、目的関数Ｇ（ｙ１、ｙ２、・・・ｙｎ）の最適解が得られることが可能に、前記流体系全体解析モデルおよび部分流動解析モデルが構築されていること、
   を特徴とする流動解析システム。
6. 請求項１ないし５のいずれかに記載の流動解析システムにおいて、
   複数のクライアント機能を有する複数のクライアント計算機と、
   ネットワークを介して前記クライアント計算機に接続されてサーバー機能を有するサーバー計算機およびストレージ機器と、
   前記ネットワークを介して前記クライアント計算機、サーバー計算機およびストレージ機器に接続され、前記流体系全体解析モデルおよび部分流動解析モデル、もしくは構造解析モデルを用いた解析を実行する高速計算機と、
   を有し、
   前記各解析モデルへの入力データの生成、解析実行、解析出力データの処理および表示機能を有し、前記流動解析を実行するためのプログラム、入力データ、解析出力データ、およびその処理データを前記サーバー機能で一元管理し、前記クライアント計算機からの利用を可能とした計算機環境で動作するように構成されていること、
   を特徴とする流動解析システム。
7. 請求項６に記載の流動解析システムにおいて、前記サーバー機能はウェブサーバー機能を含み、前記各クライアント計算機はウェブブラウザを介して前記高速計算機にアクセス可能に構成されていること、を特徴とする流動解析システム。
8. 請求項６または７に記載の流動解析システムにおいて、前記サーバー機能はデータサーバー機能を含み、そのデータサーバー機能は、前記各解析モデルの入力データ、解析出力データ、およびその処理データを、共有データベース領域と各クライアント計算機の作業領域とに分けて、両領域は解析の種類によって分類された、同一の階層構造によって前記ストレージ機器上に管理された計算機環境で動作するように構成されていること、を特徴とする流動解析システム。
9. 請求項６ないし８のいずれかに記載の流動解析システムにおいて、前記高速計算機は複数の単位高速計算機をネットワークを介して結んだグリッドコンピューティング構成であって、前記サーバー計算機から前記各単位高速計算機にある複数の解析プログラムを実行する計算機環境で動作するように構成されていること、を特徴とする流動解析システム。
10. 流体系全体を複数のノードに分割し、ノード間をジャンクションで接続し、該ノード内で状態量を計算し、該ジャンクションでノード間の移行量を計算する最上位の流動解析モデルとしての流体系全体解析モデルと、解析対象領域によって前記流体系全体を分割してかつ流動スケールに応じて階層化された複数の部分流動解析モデルとを用いて、流体系の流動を解析する流動解析方法において、
    前記流体系全体解析モデルを用いて流体系全体のシステム挙動を解析するステップと、
    前記流体系全体解析モデルによる解析結果を前記複数の部分流動解析モデルの境界条件あるいはモデルパラメーターとして利用し、前記流動スケールに応じて階層化された部分流動解析モデルの上位のモデルによる解析結果をそれより下位の部分流動解析モデルの境界条件あるいはモデルパラメーターとして利用して前記複数の部分流動解析モデルによる解析を行なうステップと、
    前記下位の部分流動解析モデルを用いた解析結果に基づいて作成される物理モデルまたは相関式モデルまたは物理パラメーターを上位の部分流動解析モデル及び前記流体系全体解析モデルの一部として組み込むステップと、
    前記流体系全体解析モデルを用いた解析結果に基づいて作成される物理モデルまたは相関式モデルまたは物理パラメーターを前記流体系全体解析モデルの一部として組み込むステップと、
    を有することを特徴とする流動解析方法。

Images