(57)【要約】【課題】制御棒が部分的に挿入されたメッシュが存在する三次元炉心の核特性シミュレーションプログラムにおいて、従来より制御棒の反応度効果を精度良く評価できる方法の提供。【解決手段】炉心の三次元分割メッシュに基づく炉心核特性シミュレーションプログラムを用いて、制御棒の反応度効果を評価する方法において、制御棒部分挿入メッシュに対する巨視的中性子断面積を当該メッシュの制御棒挿入量及び中性子束による重み付けで平均化したものを用いて制御棒反応度効果を評価する際に、前記平均化のための重み付け係数として、燃料集合体を対象とする集合体核特性評価プログラムで求められた燃料棒及び非燃料棒における中性子エネルギースペクトルを用いて制御棒案内管セルにおける制御棒挿入時と制御棒非挿入時のそれぞれについて求めた第１群の中性子エネルギースペクトルφ１  と第２群以下の中性子エネルギースペクトルφｇ  との比φｇ／φ１を用いる。【選択図】      なし【特許請求の範囲】
【請求項１】
炉心の三次元分割メッシュに基づく炉心核特性シミュレーションプログラムを用い、各燃料棒をメッシュに割り当てると共に、制御棒部分挿入メッシュに対する巨視的中性子断面積を当該メッシュの制御棒挿入量及び中性子束による重み付けで平均化して三次元炉心の核的挙動、特に制御棒の反応度効果を評価する方法において、
前記平均化のための重み付け係数として、燃料集合体を対象とする集合体核特性評価プログラムで求められた燃料棒及び非燃料棒における中性子エネルギースペクトルを用いて制御棒案内管セルにおける制御棒挿入時と制御棒非挿入時のそれぞれについて求めた第１群の中性子エネルギースペクトルφ１  と第２群以下の中性子エネルギースペクトルφｇ  との比φｇ／φ１を用いることを特徴とする炉心制御棒反応度効果評価方法。
【請求項２】
前記制御棒部分挿入メッシュに対する巨視的中性子断面積の平均化を以下の数１に基づいて演算することを特徴とする請求項１に記載の炉心制御棒反応度効果評価方法。
【数１】
2004226245.tif 000002
ただし、Σａｖｅ  はメッシュ平均巨視的断面積、ΣＲｏｄ−ｉｎは制御棒挿入部巨視的断面積、ΣＲｏｄ−ｏｕｔ  は制御棒非挿入部巨視的断面積、ＶＲｏｄ−ｉｎは制御棒挿入部体積、ＶＲｏｄ−ｏｕｔ  は制御棒非挿入部体積、φｇＲｏｄ−ｉｎ  は制御棒挿入時の制御棒案内管セルにおけるｇ群の中性子スペクトル、φ１Ｒｏｄ−ｉｎ  は制御棒挿入時の制御棒案内管セルにおける１群の中性子スペクトル、φｇＲｏｄ−ｏｕｔは制御棒非挿入時の制御棒案内管セルにおけるｇ群の中性子スペクトル、φ１Ｒｏｄ−ｏｕｔは制御棒非挿入時の制御棒案内管セルにおける１群の中性子スペクトル、である。
【請求項３】
前記三次元分割メッシュの単位メッシュ幅寸法をほぼ燃料棒１本分の径寸法に設定することを特徴とする請求項１に記載の炉心制御棒反応度効果評価方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、原子炉制御棒の反応度効果の評価方法に関するものである。
【０００２】
【従来の技術】
原子炉では、熱の発生源である核分裂の増減をコントロールするために、核分裂に寄与する中性子を吸収する物質から構成されている制御棒を用いる。従って、炉心内に制御棒を挿入すると、中性子が吸収されて核分裂が減少し、原子炉の熱出力が低下することとなり、逆に炉心内から制御棒を引く抜くと、中性子の吸収が少なくなるため核分裂が増加し、原子炉の熱出力が増加することとなる。
【０００３】
このような制御棒の炉心内への挿入及び引き抜きによって原子炉に与えられる影響は一般に反応度の添加と呼ばれている。制御棒の引き抜きによる原子炉の熱出力増加は、正の反応度添加と称され、制御棒の挿入による原子炉の熱出力減少は、負の反応度添加と称されている。
【０００４】
反応度の添加量（ρ）は、制御棒を移動させる前の炉心の実効増倍率Ｋ０ｅｆｆと制御棒を移動させた後の炉心の実効増倍率Ｋ１ｅｆｆとから、反応度添加量（ρ）＝（Ｋ１ｅｆｆ−Ｋ０ｅｆｆ）／Ｋ１ｅｆｆ・Ｋ０ｅｆｆの式から計算することができる。
【０００５】
制御棒は、車におけるアクセルとブレーキの役割を果たしているため、挿入もしくは引き抜きによる反応度添加量を正確に予測することは極めて重要である。そこで、原子炉の炉心設計においては、複数組ある制御棒の反応度添加量を予測している。
【０００６】
従来から、制御棒の反応度添加量は、原子炉の炉心の核的な振る舞いをシミュレートする計算プログラムによって計算されている。原子炉の炉心は三次元形状であることから、核的な振る舞いのシミュレーションも実際の炉心を模擬して三次元で行われるのが一般的である。
【０００７】
炉心核特性シミュレーションを行うプログラムでは、炉心はメッシュと呼ばれる小領域に分割されて核的挙動が評価されるが、制御棒の反応度効果もこのメッシュ分割されたプログラムで評価される。具体的には、各メッシュに、中性子の反応のしやすさを表すパラメータである巨視的断面積を割り当てて計算を行い、制御棒の反応度添加量が求められる。
【０００８】
計算方法によってメッシュの大きさは様々であるが、例えば加圧水型軽水炉の場合、半径方向には１〜１０ｃｍ程度、軸方向には１０〜２０ｃｍ程度の幅で分割が行われている。各メッシュでは、温度、積算熱出力、燃料タイプ、制御棒挿入の有無などが異なるため、各メッシュ毎にそれぞれ対応する巨視的断面積が割り当てられる。
【０００９】
しかしながら、炉心の上部又は下部から炉心に挿入される制御棒は、全引き抜き位置から全挿入位置までの間を連続的に移動するため、場合によっては制御棒の先端位置が軸方向のメッシュ境界に一致せず、一つのメッシュ内部、即ち当該メッシュの上下軸方向の途中に位置してしまうことがある。炉心の核特性シミュレーションプログラムでは、一つのメッシュに対しては１セットの巨視的断面積しか割り当てないため、このような１メッシュ内部に制御棒先端が位置するような場合はそのまま取り扱うことができない。
【００１０】
即ち、制御棒が部分的に挿入されたメッシュの場合、制御棒が挿入されている領域と挿入されてない領域とに対応する２セットの巨視的断面積が必要となるのに対して、炉心核特性シミュレーションプログラムでは、１セットの巨視的断面積しか使えないため、特別な取り扱いが必要となる。この問題を解決するために大きく分けて二種類の方法が提案されている。
【００１１】
第１の方法は、制御棒先端部で新たにメッシュ分割を行い、制御棒が部分的に挿入されたメッシュを、制御棒が完全に挿入された領域のみからなるメッシュと、制御棒が完全に非挿入である領域からなるメッシュとに分割する方法である。しかし、この方法では、プログラムの構造が複雑になり、反応度添加量の計算に手間がかかるため、結果としてプログラムの拡張・保守などの面で負担となるなどの問題が生じる可能性がある。
【００１２】
第２の方法は、制御棒が部分的に挿入されたメッシュ内の巨視的断面積を何らかの方法で平均化するものである。平均化を行うことでメッシュの分割を避けることができるため、プログラムの構造を簡素化することが可能となる。平均化の方法としては、さらにいくつかの方法に分類されている。
【００１３】
例えば、制御棒が部分的に挿入されたメッシュの巨視的断面積を、以下の（２）式に示すように制御棒の挿入量を重みとして平均化する方法がある。
【００１４】
Σａｖｅ  ＝（ΣＲｏｄ−ｉｎＶＲｏｄ−ｉｎ＋ΣＲｏｄ−ｏｕｔＶＲｏｄ−
ｏｕｔ）／（ＶＲｏｄ−ｉｎ＋ＶＲｏｄ−ｏｕｔ  ）                        ・・・（２）式
ここでは、Σａｖｅ  ：メッシュ平均巨視的断面積、ΣＲｏｄ−ｉｎ：制御棒挿入部巨視的断面積、ΣＲｏｄ−ｏｕｔ  ：制御棒非挿入部巨視的断面積、ＶＲｏｄ−ｉｎ：制御棒挿入部体積、ＶＲｏｄ−ｏｕｔ  ：制御棒非挿入部体積、である。
【００１５】
しかしこの方法では、メッシュに部分的に挿入された制御棒の反応度効果を過大評価し、精度良く計算できないことが知られている。そこで、以下の（３）式に示すように、制御棒が部分的に挿入されたメッシュを、制御棒の挿入量および中性子束を重みとして平均化する方法がある。
【００１６】
Σａｖｅ  ＝（ΣＲｏｄ−ｉｎＶＲｏｄ−ｉｎφＲｏｄ−ｉｎ＋ΣＲｏｄ−ｏ
ｕｔＶＲｏｄ−ｏｕｔφＲｏｄ−ｏｕｔ  ）／（ＶＲｏｄ−ｉｎφＲｏｄ−ｉｎ
＋ＶＲｏｄ−ｏｕｔφＲｏｄ−ｏｕｔ）              ・・・（３）式
ここでは、Σａｖｅ  ：メッシュ平均巨視的断面積、ΣＲｏｄ−ｉｎ：制御棒挿入部巨視的断面積、ΣＲｏｄ−ｏｕｔ  ：制御棒非挿入部巨視的断面積、ＶＲｏｄ−ｉｎ：制御棒挿入部体積、ＶＲｏｄ−ｏｕｔ  ：制御棒非挿入部体積、φＲｏｄ−ｉｎ：制御棒挿入部中性子束、φＲｏｄ−ｏｕｔ：制御棒非挿入部中性子束、である。
【００１７】
後者の平均化法を用いることで、部分的に挿入された制御棒の効果を精度良く評価可能になることが知られている。そこでこの方法が炉心の核特性シミュレーションで広く用いられている。ただし、この平均化法を用いる場合、制御棒挿入部と非挿入部の中性子束を求める必要がある。
【００１８】
一方、多くの原子炉では、直径１ｃｍ程度の燃料棒を多数束ねた集合体を単位として燃料交換を行っており、この燃料集合体の径方向の幅寸法は、例えば沸騰水型軽水炉用燃料集合体で約１０ｃｍ、加圧水型軽水炉用燃料集合体で約２０ｃｍである。
【００１９】
これに対して、広く使われている炉心の核特性シミュレーションプログラムでは、径方向のメッシュ分割幅が計算時間の観点から１０ｃｍ程度となっている（例えば、非特許文献１参照）。このようにメッシュの径方向分割幅を１０ｃｍ程度とすると、当然ながら個々の燃料棒をそれぞれ１メッシュ毎に割り当てて直接扱うことができず、燃料集合体の中を平均化して取り扱う必要が生じる。沸騰水型軽水炉では、一つの燃料集合体を径方向１メッシュに、加圧水型軽水炉では一つの燃料集合体を径方向２×２の４メッシュに分割することが通常行われている。
【００２０】
なお、前記炉心の核特性シミュレーションプログラムでは、制御棒挿入部と非挿入部の各中性子束分布を求める際、制御棒が部分的に挿入されたメッシュに注目する。即ち、図１に示すように、当該メッシュを取り出し、制御棒が挿入された部分と挿入されていない部分との軸方向二領域に分けた上で、各メッシュ領域内の中性子束を独立した一次元問題として求め、得られた中性子束を前記（３）式に利用している。
【００２１】
しかし近年、より精度良く炉心の核特性をシミュレーションするため、１メッシュの径方向幅寸法を１ｃｍ程度とし、従来のように燃料集合体内の燃料棒を平均化せず、個々の燃料棒を直接取り扱うことが試みられている（例えば、非特許文献２参照）。このような試みは、特に加圧水型軽水炉において盛んに行われている。
【００２２】
【非特許文献１】
Ｋ．Ｓ．Ｓｍｉｔｈ，”Ａｓｓｅｍｂｌｙ  Ｈｏｍｏｇｅｎｉｚａｔｉｏｎ  Ｔｅｃｈｎｉｑｕｅ  ｆｏｒ  Ｌｉｇｈｔ  Ｗａｔｅｒ  Ｒｅａｃｔｏｒ  Ａｎａｌｙｓｉｓ，”  Ｐｒｏｇｒｅｓｓ  Ｎｕｃｌｅａｒ  Ｅｎｅｒｇｙ，１７，３０３（１９８０）．
【非特許文献２】
Ｍ．Ｔａｔｓｕｍｉ，Ａ．Ｙａｍａｍｏｔｏ，”Ｏｂｊｅｃｔ−Ｏｒｉｅｎｔｅｄ  Ｔｈｒｅｅ−Ｄｉｍｅｎｓｉｏｎａｌ  Ｆｉｎｅ−Ｍｅｓｈ  Ｔｒａｎｓｐｏｒｔ  Ｃａｌｃｕｌａｔｉｏｎ  ｏｎ  Ｐａｒａｌｌｅｌ／Ｄｉｓｔｒｉｂｕｔｅｄ  Ｅｎｖｉｒｏｎｍｅｎｔｓ  ｆｏｒ  Ａｄｖａｎｃｅｄ  Ｒｅａｃｔｏｒ  Ｃｏｒｅ  Ａｎａｌｙｓｉｓ，”  Ｎｕｃｌｅａｒ  Ｓｃｉｅｎｃｅ  ａｎｄ  Ｅｎｇｉｎｅｅｒｉｎｇ，１４１，１９０（２００２）．
【００２３】
【発明が解決しようとする課題】
しかしながら、前記のように１メッシュの径方向幅寸法を１０ｃｍ程度として取り扱っている従来のプログラムにおいて、制御棒部分挿入メッシュを制御棒挿入部と非挿入部との二領域に分けて一次元問題として取り扱うことで精度良い評価が可能であったのは、１メッシュの径方向幅寸法が比較的大きく、集合体内部が平均化されたためであり、このような方法は、個々の燃料棒を直接扱う場合には直接適用できない。
【００２４】
これは、個々の燃料棒を直接取り扱う場合、１メッシュの径方向幅寸法が１ｃｍ程度、軸方向の長さ寸法が１０〜２０ｃｍとなり、当該メッシュにつき横（幅）方向からの中性子の流入・流出の影響が大きくなることから、これまでと同様の軸方向一次元での取り扱いでは、精度良い解析は困難であるためであり、実際、１メッシュの幅寸法を約１ｃｍとする設定で中性子の流入・流出を考慮しない一次元計算では、（３）式で用いる重みの中性子束を正確に評価できないことが確認されている。
【００２５】
横方向からの中性子の流入・流出を精度良く評価するためには、図２に示すように隣接燃料を考慮した三次元体系で中性子束の分布を計算する必要が生じる。即ち、従来のプログラムを用いた方法では、軸方向一次元計算で精度良い評価が行えたのに比べて、１メッシュの径方向の幅を１ｃｍ程度にして個々の燃料棒が割り当てられる設定にした場合には三次元計算が必要となるため、計算時間が多大になるだけでなく、そのためのプログラムが必要となり、作業が煩雑になる。
【００２６】
本発明の目的は、上記問題点に鑑み、制御棒が部分的に挿入されたメッシュが存在する三次元炉心の核特性シミュレーションプログラムにおいて、従来より制御棒の反応度効果を精度良く評価することのできる方法を提供することにある。
【００２７】
【課題を解決するための手段】
上記目的を達成するため、請求項１に記載の発明に係る炉心制御棒反応度効果評価方法は、炉心の三次元分割メッシュに基づく炉心核特性シミュレーションプログラムを用い、各燃料棒をメッシュに割り当てると共に、制御棒部分挿入メッシュに対する巨視的中性子断面積を当該メッシュの制御棒挿入量及び中性子束による重み付けで平均化して三次元炉心の核的挙動、特に制御棒の反応度効果を評価する方法において、前記平均化のための重み付け係数として、燃料集合体を対象とする集合体核特性評価プログラムで求められた燃料棒及び非燃料棒における中性子エネルギースペクトルを用いて制御棒案内管セルにおける制御棒挿入時と制御棒非挿入時のそれぞれについて求めた第１群の中性子エネルギースペクトルφ１  と第２群以下の中性子エネルギースペクトルφｇ  との比φｇ／φ１を用いるものである。
【００２８】
また、請求項２に記載の発明に係る炉心制御棒反応度効果評価方法は、請求項１に記載の炉心制御棒反応度効果評価方法において、前記制御棒部分挿入メッシュに対する巨視的中性子断面積の平均化を数１に基づいて演算するものである。
【数１】
Σａｖｅ  ＝｛ΣＲｏｄ−ｉｎＶＲｏｄ−ｉｎ（φｇＲｏｄ−ｉｎ／φ１Ｒｏｄ−ｉｎ）＋ΣＲｏｄ−ｏｕｔＶＲｏｄ−ｏｕｔ（φｇＲｏｄ−ｏｕｔ／φ１Ｒｏｄ−ｏｕｔ  ）｝／｛ＶＲｏｄ−ｉｎ（φｇＲｏｄ−ｉｎ／φ１Ｒｏｄ−ｉｎ）＋ＶＲｏｄ−ｏｕｔ（φｇＲｏｄ−ｏｕｔ／φ１Ｒｏｄ−ｏｕｔ  ）｝
ただし、Σａｖｅ  はメッシュ平均巨視的断面積、ΣＲｏｄ−ｉｎは制御棒挿入部巨視的断面積、ΣＲｏｄ−ｏｕｔ  は制御棒非挿入部巨視的断面積、ＶＲｏｄ−ｉｎは制御棒挿入部体積、ＶＲｏｄ−ｏｕｔ  は制御棒非挿入部体積、φｇＲｏｄ−ｉｎ  は制御棒挿入時の制御棒案内管セルにおけるｇ群の中性子スペクトル、φ１Ｒｏｄ−ｉｎ  は制御棒挿入時の制御棒案内管セルにおける１群の中性子スペクトル、φｇＲｏｄ−ｏｕｔは制御棒非挿入時の制御棒案内管セルにおけるｇ群の中性子スペクトル、φ１Ｒｏｄ−ｏｕｔは制御棒非挿入時の制御棒案内管セルにおける１群の中性子スペクトル、である。
【００２９】
また請求項３に記載の発明に係る炉心制御反応度効果評価方法は、請求項１に記載の炉心制御棒反応度効果評価方法において、前記三次元分割メッシュの単位メッシュ幅寸法をほぼ燃料棒１本分の径寸法に設定するものである。
【００３０】
【発明の実施の形態】
本発明は、炉心の三次元分割メッシュに基づく核特性シミュレーションプログラムを用いて核的挙動を評価する際に、各メッシュに対応する巨視的断面積を割り当てて制御棒の反応度添加量を計算し、制御棒反応度効果を評価するものであり、径方向（炉心軸方向と直交する方向）に関して各燃料棒をメッシュに割り当てて個々の燃料棒を直接取り扱う設定においても、制御棒反応度効果を精度良く評価するためのものである。
【００３１】
即ち、個々の燃料棒を扱う設定での核特性シミュレーションプログラムにおいては、制御棒部分挿入メッシュに対する巨視的断面積として当該メッシュの制御棒挿入量及び中性子束による重み付けで平均化したものを用いる際に、従来はこの平均化巨視的中性子断面積を精度良く求めるのが困難であったのに対して、炉心の核的挙動を計算するプログラムでは燃料集合体を対象とする集合体核特性評価プログラムにより計算された巨視的断面積が使用されており、この集合体核特性評価プログラムでは全ての燃料棒および非燃料棒における中性子のエネルギー分布（中性子のエネルギースペクトル）を求めていることに着目し、本発明においては、この中性子エネルギースペクトルを前記重み係数に利用することによって、制御棒部分挿入メッシュの巨視的断面積の平均化を精度良く行うことができることを見出し、従って制御棒反応度添加効果の評価を精度良く行うことを可能とするものである。
【００３２】
一般に、炉心の核的挙動を計算するプログラムでは、中性子のエネルギー群は２群ないしそれより多い群数で計算され、１群目は高速群となっている。高速群では中性子の飛行距離が長いため、メッシュ内で中性子束の分布が平坦であると仮定できる場合が多い。そこで、制御棒が部分的に挿入されたメッシュにおいても、第１群の中性子束は平坦である、即ちメッシュ内で分布がないと仮定する。さらに、制御棒が挿入された部分と挿入されていない部分とのそれぞれにおいて、中性子エネルギースペクトルが集合体核特性評価プログラムで求められたものと大きく異なっていないと仮定する。
【００３３】
これら二つの仮定をおくと、前記（３）式で重みに用いている中性子束φＲｏｄ−ｉｎ、及びφＲｏｄ−ｏｕｔ  は、それぞれ集合体核特性評価プログラムの計算結果から求められるφｇＲｏｄ−ｉｎ／φ１Ｒｏｄ−ｉｎ・・・（４）式、及びφｇＲｏｄ−ｏｕｔ／φ１Ｒｏｄ−ｏｕｔ・・・（５）式、の重み係数で精度良く置き換え（数１を得る）ることができる。
【００３４】
ここで、φｇＲｏｄ−ｉｎ  は制御棒挿入時の制御棒案内管セルにおけるｇ群（２群以下）の中性子スペクトル、φ１Ｒｏｄ−ｉｎ  は制御棒挿入時の制御棒案内管セルにおける１群の中性子スペクトル、また、φｇＲｏｄ−ｏｕｔは制御棒未挿入時の制御棒案内管セルにおけるｇ群の中性子スペクトル、φ１Ｒｏｄ−ｏｕｔは制御棒未挿入時の制御棒案内管セルにおける１群の中性子スペクトル、である。
【００３５】
以上の（４）式および（５）式で求められる第１群の中性子エネルギースペクトルφ１と第２群以下（ｇ群）の中性子エネルギースペクトルφｇとの比φｇ／φ１を重み係数に用いると、第１群に対する巨視的断面積は体積平均となり、２群以下（ｇ群）は、制御棒の有無による中性子束のエネルギースペクトルの変化を考慮したものとなる。（４）式、（５）式は、核特性評価で用いるスペクトルインデックスの逆数となっているため、このスペクトルインデックス逆数を制御棒部分挿入メッシュ内の巨視的断面積の平均化における重み付け係数として用いる方法を逆スペクトルインデックス重み法と呼ぶことができる。
【００３６】
このような逆スペクトルインデックス重み法を用いれば、径方向に関して各燃料棒を１メッシュ毎に割り当てて個々の燃料棒を直接扱う設定での炉心核特性シミュレーションプログラムにおいて、（３）式で表した制御棒部分挿入メッシュの巨視的断面積の平均化の重み係数としての中性子束を求めるために多大な計算時間が掛かる三次元計算をおこなう必要なく、精度良く制御棒の反応度効果の評価を行うことができる。
【００３７】
なお、個々の燃料棒を直接取り扱うための炉心の三次元メッシュ分割は、分割メッシュの単位メッシュ幅寸法をほぼ燃料棒１本分の径寸法に設定すれば良い。
【００３８】
一方、上記の本発明における逆スペクトルインデックス重み法は、径方向に関して１燃料棒を１メッシュ毎に割り当て個々の燃料棒を直接取り扱う炉心核特性シミュレーションプログラムによる評価の場合に限らず、軸方向に異なる組成の燃料が存在する軸方向非均質燃料に対しても適用可能である。即ち、組成の境界がメッシュ境界に一致していなくても制御棒の場合と同様の考え方で本方法を適用できる。
【００３９】
従って、従来の径方向に関して１燃料集合体を１メッシュに割り当てるプログラムにおいても、この逆スペクトルインデックス重み法を用いれば、制御棒部分挿入メッシュを挿入部と非挿入部との二領域に分けてそれぞれ行っていた一次元計算までも必要なく、メッシュ平均巨視的断面積を簡便に精度良く求めることができ、制御棒の反応度効果を精度良く評価することが可能となる。
【００４０】
【実施例】
本発明の一実施例として、燃料集合体３行３列の小型炉心を想定し、制御棒の反応度添加量の計算を行った場合を以下に説明する。本実施例においては、使用燃料は、加圧水型軽水炉で用いられている１７×１７型のものであり、濃縮度が２ｗｔ％のものを想定した。
【００４１】
また制御棒は中央の集合体に炉心上部から挿入されるものとした。１燃料集合体のサイズは径方向が約２０ｃｍ、軸方向が約３６０ｃｍである。メッシュ分割は、径方向に関して単位メッシュ幅寸法をほぼ燃料棒単位に一致（約１ｃｍ）させ、軸方向に関しては等分割の２０メッシュ（１メッシュ約１８ｃｍ）とした。
【００４２】
以上の炉心設計において、制御棒の挿入は、その先端が炉心上部から６メッシュ目内に位置するものとし、この制御棒が部分挿入された第６メッシュについて前記数１の式を用いて巨視的断面積の平均化処理を行った上で反応度添加量を求めた。
【００４３】
ここで、本実施例の比較対照として、同じ第６メッシュに関して、従来の体積重み法による前記（２）式を用いて巨視的断面積の平均化処理を行った上で反応度添加量も求め、それぞれの結果を図３に示した。なお、リファレンスとして、当該メッシュを制御棒の先端部位置で制御棒挿入領域と非挿入領域とに分け（このとき軸方向メッシュ数は２１個となる）て計算を行った結果に基づいても反応度添加量を求め、これも同様に図３に示した。この場合、対象メッシュ内における制御棒の平均化を行っておらず、最も正確な解が与えられる。
【００４４】
本実施例における、メッシュ平均巨視的断面積の解、およびこれを割り当てて反応度添加量を得るまでの演算は、図４のフローチャートに示す手順の通りである。
【００４５】
図３に示す結果から明らかなように、図中白四角（□）で表した本実施例における逆スペクトルインデックス重み法による巨視的断面積の平均化に基づく制御棒の反応度添加量の解（図中縦軸に反応度効果Δｋ／ｋとして示す）は、図中黒三角（▲）で表したリファレンス解とほぼ一致しており、本発明の方法によれば、制御棒の反応度効果評価は、多大な計算時間を必要とせず、簡便でありながらも高精度なものであることが確認できた。
【００４６】
なお、炉心の核特性プログラムでは、メッシュ内の中性子束分布を多項式などの関数で展開することがよく行われる。以上の実施例においては、第１群の中性子エネルギースペクトルと第２群以下の中性子エネルギースペクトルとの比を、（４）式及び（５）式において第１群の中性子束を平坦であると仮定して導出したが、中性子束分布を多項式などの関数で展開している場合、制御棒挿入領域と非挿入領域とで展開関数を積分することで制御棒部分挿入メッシュ内の制御棒挿入領域と非挿入領域における第１群の中性子束の変化割合を計算して重み係数に反映することも可能である。
【００４７】
また、注目するメッシュと軸方向の上下に隣接するメッシュの平均中性子束を多項式の関数でフィッティングし、得られた多項式を用いて前記メッシュ内の制御棒挿入領域と非挿入領域における第１群の中性子束の変化割合を計算して重み係数に反映することもできる。さらに、Ｊ．Ｃ．Ｇｅｈｉｎ，”Ａ  Ｑｕａｓｉ  Ｓｔａｔｉｃ  ＰｏｌｙｎｏｍｉａｌＮｏｄａｌ  Ｍｅｔｈｏｄ  ｆｏｒ  Ｎｕｃｌｅｓｒ  Ｒｅａｃｔｏｒ  Ａｎａｌｙｓｉｓ，”  Ｐｈ．Ｄ．Ｔｈｅｓｉｓ，ＭＩＴ，ＵＳＡ，ＰＰ．５９−６０，（１９９２）．  に示す方法で制御棒挿入領域と非挿入領域における第１群の中性子束の変化割合を計算し、重み係数に反映することもできる。
【００４８】
さらに、注目するメッシュの平均中性子束と、上下端における表面中性子束を多項式等の関数でフィッティングし、得られた多項式を用いて前記メッシュ内の制御棒挿入領域と非挿入領域における第１群の中性子束の変化割合を計算し、重み係数に反映することもできる。
【００４９】
なお、本発明は、制御棒が部分的に挿入されたメッシュに関して考察したものであるが、本発明の考え方は、メッシュ内に制御棒の境界だけでなく、他の物質境界が存在する場合への応用も可能であると考えられる。例えば、沸騰水型軽水炉では燃料集合体の軸方向に燃料濃縮度や可燃性毒物濃度などの分布がつけられており、燃料濃縮度の境界がメッシュ境界と一致しないでメッシュ内部に位置する場合などにも、炉心の核特性評価プログラムにおいて、中性子エネルギースペクトル比を重み係数に用いることで、付加的な計算を行うことなく精度よく炉心特性の評価を行える可能性がある。
【００５０】
【発明の効果】
以上説明したとおり、本発明の方法によれば、制御棒が部分的に挿入されてたメッシュが存在する炉心の核特性評価プログラムにおいても、時間の掛かる付加的な計算を行うことなく、簡便に精度の良い制御棒の反応度効果を評価することができるという効果がある。
【図面の簡単な説明】
【図１】従来の三次元分割メッシュに基づく炉心核特性評価プログラムにおける制御棒部分挿入メッシュの取り扱いを説明する模式図である。
【図２】詳細なメッシュ分割に基づく炉心核特性評価プログラムにおける制御棒部分挿入メッシュの取り扱いを説明する模式図である。
【図３】本発明の一実施例による小型炉心における制御棒部分挿入メッシュの制御棒の反応度価値の計算結果を示すグラフ図であり、横軸の当該メッシュ下端からの距離に応じた反応度効果（Δｋ／ｋ）を縦軸に表したものである。
【図４】本実施例における制御棒の反応度価値を算出するための手順を示すフローチャート図である。
【図１】
2004226245.tif 000003
【図２】
2004226245.tif 000004
【図３】
2004226245.tif 000005
【図４】
2004226245.tif 000006