JPS6086489A - 原子炉炉心シミユレ−シヨン方式 - Google Patents
原子炉炉心シミユレ−シヨン方式Info
- Publication number
- JPS6086489A JPS6086489A JP58194156A JP19415683A JPS6086489A JP S6086489 A JPS6086489 A JP S6086489A JP 58194156 A JP58194156 A JP 58194156A JP 19415683 A JP19415683 A JP 19415683A JP S6086489 A JPS6086489 A JP S6086489A
- Authority
- JP
- Japan
- Prior art keywords
- output
- control rod
- value
- notches
- current
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔発明の利用分野〕
本発明は、原子炉炉心内局所出力領域モニタ出ノ月ｌ［
ｆｔ−制御棒バタンに対応して速やかに模擬計算によシ
得られる、原子炉炉心シミュレーション方式に関するも
のである。
ｆｔ−制御棒バタンに対応して速やかに模擬計算によシ
得られる、原子炉炉心シミュレーション方式に関するも
のである。
従来の原子炉炉心シミュレーション方式〇よ、主に、核
熱水力モデルを使用している。この方式によれば、演算
方法は行列演算を主としたくシかえし算法になっている
ので、計算時間が長く、実時間でプロセスｊＩｋを模擬
しその模擬出力値葡求めることは困−でめった。この理
由によシ、原子炉の運転状態を実時間で模擬した運転−
１１練用シミユレータを構成することができなかった。
熱水力モデルを使用している。この方式によれば、演算
方法は行列演算を主としたくシかえし算法になっている
ので、計算時間が長く、実時間でプロセスｊＩｋを模擬
しその模擬出力値葡求めることは困−でめった。この理
由によシ、原子炉の運転状態を実時間で模擬した運転−
１１練用シミユレータを構成することができなかった。
本発明の目的は、原子炉炉心内に配置された局所出力領
域モ＝り（Ｉ、ｏｃａｌ ｐｏｗｅｒ几ａｎｔｅＭｏｎ
ｉｔｏｒ、ＬＰＲＭと以下称す）の出力値を、現状の制
御棒パタ／に対応して実時間で速やかに演算によシ原子
炉炉心シミュレーション方式を提供することにある。
域モ＝り（Ｉ、ｏｃａｌ ｐｏｗｅｒ几ａｎｔｅＭｏｎ
ｉｔｏｒ、ＬＰＲＭと以下称す）の出力値を、現状の制
御棒パタ／に対応して実時間で速やかに演算によシ原子
炉炉心シミュレーション方式を提供することにある。
上記の目的のため本発明は、制御棒の単位長引抜による
ＬＰ几Ｍ位置でのＬＰｆＬＭ出力値の変化率を用い、Ｌ
Ｐ几１１−囲む限定した煩域の制御棒のパタンの変化を
制御棒の引抜／そう人長さとその置さ方間位置で検出し
、ＬＰＲＭの出力値に変化率を引抜／そう人長さに応じ
重ね合せて乗じ、更にＬＰＲＭを囲む限定した領域ごと
に、距離による補正のために係数を乗じた直をもって、
現状のＬＰＲＭ出力１直とする方式とした。
ＬＰ几Ｍ位置でのＬＰｆＬＭ出力値の変化率を用い、Ｌ
Ｐ几１１−囲む限定した煩域の制御棒のパタンの変化を
制御棒の引抜／そう人長さとその置さ方間位置で検出し
、ＬＰＲＭの出力値に変化率を引抜／そう人長さに応じ
重ね合せて乗じ、更にＬＰＲＭを囲む限定した領域ごと
に、距離による補正のために係数を乗じた直をもって、
現状のＬＰＲＭ出力１直とする方式とした。
以下、本発明を第１図から第 図によシ説明する。
第１図は、本発明の原子炉炉心シミュレーション方式に
よ、！７ＬＰＲＭ出力値を演算・出力する処理手順の流
れ図でおる。
よ、！７ＬＰＲＭ出力値を演算・出力する処理手順の流
れ図でおる。
まず、代表的な出力における制御棒パタ／と、その出力
でのＬＰＲＭ出力値金参照値として入力する。Ｉｌｏｏ
ＭＷ、級ＢＷＲにおける購御棒総数は１８５本である。
でのＬＰＲＭ出力値金参照値として入力する。Ｉｌｏｏ
ＭＷ、級ＢＷＲにおける購御棒総数は１８５本である。
ＬＰＲＭは制御棒のそう人されない燃料チャンネル間の
、燃料チャンネル４体に一本が割当てられる位置（以下
ストリング位置と呼ぶ）４３ケ所に、高さ方向４点に配
置されている。高さ方向に、炉心を２４分割したノード
の、となシ合う上下ノード間の高式にＬＰＲＭは置かれ
る。下から順に第３ノードと第４ノードの間に最下位置
のＬＰＲＭがおかれ、第９ノードの間に下から３番目の
ＬＰＲＭが、そして、第２１ノードと第２２ノードの間
に最上位置のＬ ｌ）几Ｍが置かれている。
、燃料チャンネル４体に一本が割当てられる位置（以下
ストリング位置と呼ぶ）４３ケ所に、高さ方向４点に配
置されている。高さ方向に、炉心を２４分割したノード
の、となシ合う上下ノード間の高式にＬＰＲＭは置かれ
る。下から順に第３ノードと第４ノードの間に最下位置
のＬＰＲＭがおかれ、第９ノードの間に下から３番目の
ＬＰＲＭが、そして、第２１ノードと第２２ノードの間
に最上位置のＬ ｌ）几Ｍが置かれている。
制御棒の引抜ノツチ数は、完全そう人時に０、完全引抜
時に４８となるように引抜の有効炉心内憂を４８分割し
て表わし、制御棒総数１８５個用意する。
時に４８となるように引抜の有効炉心内憂を４８分割し
て表わし、制御棒総数１８５個用意する。
ＬＰＲＭの出力値は、相対値で表わし、炉心の４３本の
ストリング位置の高さ方向４点の１７２個用意する。
ストリング位置の高さ方向４点の１７２個用意する。
以上の、制御棒引抜ノツチ数に関する数値と、ＬＰ几Ｍ
出力値の参照値を、代表的な出力レベルに対応させ入力
する。以下の実施例では出力レベルを、２５％、５０％
、７５％、１００％とする。
出力値の参照値を、代表的な出力レベルに対応させ入力
する。以下の実施例では出力レベルを、２５％、５０％
、７５％、１００％とする。
炉心シュミレータとして、本発明では任意の制御棒パタ
ンを入力すれば、ＬＰＲＭ出力値をＬＰ几Ｍ出力値の参
照値を用いて、後述する方式による内そう計算によシ得
る。そのため、次のステップでは、制御棒バタン金入力
して、その引抜ノツチ数金計算して判別する。
ンを入力すれば、ＬＰＲＭ出力値をＬＰ几Ｍ出力値の参
照値を用いて、後述する方式による内そう計算によシ得
る。そのため、次のステップでは、制御棒バタン金入力
して、その引抜ノツチ数金計算して判別する。
（１）もし、引抜ノツチ数が、出力レベル２５％。
５０％、７５％、１００％といずれかでの制御棒引抜ノ
ツチ数であれば、その出力レベルでのＬＰ几Ｍ出力値の
参照＋ｉを規格化した値に、出力レベル（％）をｌＯで
割った値忙かけて、ＬＰ几Ｍ出力値とする。更に、 （２）もし引抜ノツチ総数が、２５％、５０九、７５％
、１００％の出力レベルの中のとなシ甘うレベルでの引
抜ノツチ総数の間の値ｔ−取シ、シかも、この引抜ノツ
チ総、数が、シミュレータの起動後はじめて２５％、５
０％、７５％、１００％での引抜ノツチ総数のいずれか
の、間の匝を取った場合、（３）もし、引抜ノツチ数が
、２５３１６．５０％、７５％、１００％の出力レベル
の中のとなシ会うレベルでの引抜ノツチ総数の間の値を
取シ、しかも、シミュレータの起動後のＬＰＢＭ出力値
計ｎ時に同じ引抜ノツチ総数間の値をこれまでに取った
ことがめる場合、に分けてＬＰＩ（Ｍ出力値ｅ ｉ４’
ｌ詳すゐ。
ツチ数であれば、その出力レベルでのＬＰ几Ｍ出力値の
参照＋ｉを規格化した値に、出力レベル（％）をｌＯで
割った値忙かけて、ＬＰ几Ｍ出力値とする。更に、 （２）もし引抜ノツチ総数が、２５％、５０九、７５％
、１００％の出力レベルの中のとなシ甘うレベルでの引
抜ノツチ総数の間の値ｔ−取シ、シかも、この引抜ノツ
チ総、数が、シミュレータの起動後はじめて２５％、５
０％、７５％、１００％での引抜ノツチ総数のいずれか
の、間の匝を取った場合、（３）もし、引抜ノツチ数が
、２５３１６．５０％、７５％、１００％の出力レベル
の中のとなシ会うレベルでの引抜ノツチ総数の間の値を
取シ、しかも、シミュレータの起動後のＬＰＢＭ出力値
計ｎ時に同じ引抜ノツチ総数間の値をこれまでに取った
ことがめる場合、に分けてＬＰＩ（Ｍ出力値ｅ ｉ４’
ｌ詳すゐ。
以下では、上記（２）・（３）での処理の手順の笑ｈｔ
ｉＩ例を説明する。
ｉＩ例を説明する。
上記（２）の例えでは、次の処理手順を実行する。
制ｒＩ＆ｔ１棹引抜総ノツチ数Ｎが、１ヒとえは、５０
夕Ｃ出力レベルでの制御棒パタンノツチ数ＮＳＯと７５
％出力レベルでの制御棒バタンノツチＭ Ｎｙ ａの間
にるる場合、すなわち、 Ｎｓ ｏ ＜ Ｎ ＜ Ｎｔ ａとなる場合、と考える
。後述する方式による内そう計算により、まず５０％出
力レベルでのＬＰ几Ｍ出力値から得た、現制御棒総ノツ
チ数Ｎに対応する制御棒バタンでのＬＰＲＭ出力値を相
対値も。とじてめる。同様に、７５％出力レベルでのＬ
Ｐ几Ｍ出力値から、ＬＰＡＭ出力値全相対値Ｒ７！とし
てめる。Ｂ、。。
夕Ｃ出力レベルでの制御棒パタンノツチ数ＮＳＯと７５
％出力レベルでの制御棒バタンノツチＭ Ｎｙ ａの間
にるる場合、すなわち、 Ｎｓ ｏ ＜ Ｎ ＜ Ｎｔ ａとなる場合、と考える
。後述する方式による内そう計算により、まず５０％出
力レベルでのＬＰ几Ｍ出力値から得た、現制御棒総ノツ
チ数Ｎに対応する制御棒バタンでのＬＰＲＭ出力値を相
対値も。とじてめる。同様に、７５％出力レベルでのＬ
Ｐ几Ｍ出力値から、ＬＰＡＭ出力値全相対値Ｒ７！とし
てめる。Ｂ、。。
几７Ｂは、１７２個のＬＰＲＭについて計算される。
しかるに、現制御棒総ノツチ数Ｎに対応する制御棒パタ
ン時のＬＰ几Ｍ出力値は、まず相対値几を、によ請求め
、この相対値Ｒｅ、１７２個のＬＰＲＭについての平均
値によシ割った値に、現制御棒総ノンチ故Ｎにより定ま
る出力レベル係数Ｆをかけた値とする。出力レベル係数
Ｆは、ＬＰ１’ＬＭ出力値の規格化後の平均値が、２５
％出力レベルで１５．１００％出力レベルで１０．０と
なるように定め、次式で与える。
ン時のＬＰ几Ｍ出力値は、まず相対値几を、によ請求め
、この相対値Ｒｅ、１７２個のＬＰＲＭについての平均
値によシ割った値に、現制御棒総ノンチ故Ｎにより定ま
る出力レベル係数Ｆをかけた値とする。出力レベル係数
Ｆは、ＬＰ１’ＬＭ出力値の規格化後の平均値が、２５
％出力レベルで１５．１００％出力レベルで１０．０と
なるように定め、次式で与える。
１２、Ｎ１００ｉｌｌ：１００％出力レベルでの制御棒
引接線ノツチ数を表わす。
引接線ノツチ数を表わす。
前記（３）の例では、次の処理手順を実行する。前記（
３）の例は、すでに、５０％と７５％出力レベルを与え
る制御棒バタンとＬＰＲＭ出ノ月直とから、前記（２）
の場合を経て２つのＬＰ几Ｍ出力値の相対値Ｒ５゜、几
７１が、得られてい／）場合でるる。制御棒が、制御棒
引抜総ノンテ数Ｎをこえ、さらに引抜かれ、しかも、 Ｎ、。＜Ｎ’＜Ｎｔｓ を満たす制御棒引抜総ノンチ数Ｎ′になった時、前記（
２）の手順でＮをＮ’とおきかえて処理するかわシに、
几う。、几７Ｉｌから、後述する方式による内そう計算
によシ、制御棒引抜総ノンチ赦Ｎに対応する制御棒バタ
ンからＮ′に対応する制御棒バタンに変化した後のＬＰ
几Ｍ出力値の相対値をめる。Ｒ５゜からめられｆｃＭｉ
ｒｆｃな相対値Ｒ５゜′。
３）の例は、すでに、５０％と７５％出力レベルを与え
る制御棒バタンとＬＰＲＭ出ノ月直とから、前記（２）
の場合を経て２つのＬＰ几Ｍ出力値の相対値Ｒ５゜、几
７１が、得られてい／）場合でるる。制御棒が、制御棒
引抜総ノンテ数Ｎをこえ、さらに引抜かれ、しかも、 Ｎ、。＜Ｎ’＜Ｎｔｓ を満たす制御棒引抜総ノンチ数Ｎ′になった時、前記（
２）の手順でＮをＮ’とおきかえて処理するかわシに、
几う。、几７Ｉｌから、後述する方式による内そう計算
によシ、制御棒引抜総ノンチ赦Ｎに対応する制御棒バタ
ンからＮ′に対応する制御棒バタンに変化した後のＬＰ
几Ｍ出力値の相対値をめる。Ｒ５゜からめられｆｃＭｉ
ｒｆｃな相対値Ｒ５゜′。
几７．からめられた新たな相対値ル。′を用いて、現制
御棒総ノツチ数Ｎ′に対応する制御棒バタン時のＬＰＲ
Ｍ出力直出力対値几′は、 によシ氷めゐ。前記（２）の場合と同様に、各ＬＰ几Ｍ
出力１直は、Ｎ′をＮ′の１７５個のＬＰＲＭについて
の平均値によシ割りた値に、出力レベル係数Ｆ′をかけ
た値とする。出力レベル係数Ｆ′は、この場合、 によシ与える。
御棒総ノツチ数Ｎ′に対応する制御棒バタン時のＬＰＲ
Ｍ出力直出力対値几′は、 によシ氷めゐ。前記（２）の場合と同様に、各ＬＰ几Ｍ
出力１直は、Ｎ′をＮ′の１７５個のＬＰＲＭについて
の平均値によシ割りた値に、出力レベル係数Ｆ′をかけ
た値とする。出力レベル係数Ｆ′は、この場合、 によシ与える。
以上の、（３）の場合は、制御棒引接線ノツチ数Ｎ ’
が、Ｎ、。（Ｎ’ ＜ Ｎｕｌｌを満たすＮ′を継続し
てとゐような制御棒バタンか入力として与えられる間、
〈りかえして実行される処理手順である。
が、Ｎ、。（Ｎ’ ＜ Ｎｕｌｌを満たすＮ′を継続し
てとゐような制御棒バタンか入力として与えられる間、
〈りかえして実行される処理手順である。
更に、制御棒引抜が進み、Ｎ′がＮ７．に等しくなった
時は、前記（１）の場合の処理手順を実行し、Ｎ、。＜
Ｎ ’ ＜ Ｎｔ ｓ を満たすＮ′の次に、 Ｎ□＜ Ｎ” ＜ Ｎｌ。。
時は、前記（１）の場合の処理手順を実行し、Ｎ、。＜
Ｎ ’ ＜ Ｎｔ ｓ を満たすＮ′の次に、 Ｎ□＜ Ｎ” ＜ Ｎｌ。。
を満たすＮ”を与える制御棒バタンか人力されると、前
記（２）の場合に該当し、その処理手順を実行する。
記（２）の場合に該当し、その処理手順を実行する。
以上の実施例に２ける、ＬＰＲＭ出力値の内そう計算方
式について、以下説明する。以下説明するＬＰＲＭ出力
値の内そう計算方式とは、制御棒バタンとその制御棒パ
タン時のＬＰＲＭ出力値の相対値をもとに、制御棒のそ
う入拳引抜操作後に変化した制御棒バタンとなった時の
ＬＰＲＭ出力値の相対値をめる計算方式である。
式について、以下説明する。以下説明するＬＰＲＭ出力
値の内そう計算方式とは、制御棒バタンとその制御棒パ
タン時のＬＰＲＭ出力値の相対値をもとに、制御棒のそ
う入拳引抜操作後に変化した制御棒バタンとなった時の
ＬＰＲＭ出力値の相対値をめる計算方式である。
第２図に、ＬＰ几Ｍ出力値の変化率を与える制御棒引抜
位置と、ＬＰＲＭ位置の位置関係を示す。
位置と、ＬＰＲＭ位置の位置関係を示す。
ＬＸ及至Ｌ６は、炉心を高さ方向２４分割したノード間
位置を示し、Ｂ６が炉心上方に位置する。
位置を示し、Ｂ６が炉心上方に位置する。
Ｂ１及至Ｂ７は、分割され九ノードを示すとともに、２
ノツチ引抜・そう大分の制御棒部分を示す。
ノツチ引抜・そう大分の制御棒部分を示す。
１ノードは、制御棒２ノツチに相当する。いま、制御棒
先端がノードＢ４とＢ５の間から、ノードＲ３と８４の
間に１ノード移動した場合を考える。
先端がノードＢ４とＢ５の間から、ノードＲ３と８４の
間に１ノード移動した場合を考える。
中性子吸収体がノードＢ４から消失したことによシ、こ
のノードをはさむ上下６点のノード間位置Ｌ１からＲ６
各点にもし局所出力績域モニタを設・けた時のその出力
値の変化率は、ＬｌではＲ１゜Ｒ２では几６．Ｌ３では
Ｒｓ、Ｒ４ではＲ４゜Ｒ５ではＲ，、、Ｒ６ではＲ６と
与える。この変化率Ｒ１及至Ｒ０を用いて、複数ノード
数にわたる制御棒の引抜・そう人によるＬＰ几Ｍ出力値
の変化４をめる行列ＡＩｊを次のように定める。
のノードをはさむ上下６点のノード間位置Ｌ１からＲ６
各点にもし局所出力績域モニタを設・けた時のその出力
値の変化率は、ＬｌではＲ１゜Ｒ２では几６．Ｌ３では
Ｒｓ、Ｒ４ではＲ４゜Ｒ５ではＲ，、、Ｒ６ではＲ６と
与える。この変化率Ｒ１及至Ｒ０を用いて、複数ノード
数にわたる制御棒の引抜・そう人によるＬＰ几Ｍ出力値
の変化４をめる行列ＡＩｊを次のように定める。
ＡＩＪは、７行７列の正方行列で、その１６１列要素の
値ａＩｊは、 ａｌｌ ＝＝Ｑ、 Ｏｊ＝ｊ ａｕ ＝Ｒｔ ｊ＝ｉ＋ｌ、 ｉ＝１．２．３．４．５
．６ａＢ−几ＩＸ几１．１ ＋ Ｊ−’＋２Ｉｉ＝Ｌ
２＋ ３．’ｓ ５ａＩＨ＝Ｒｔ ＸＲ，＋ ＋３ Ｘ
Ｒｔ ＋ｔ ＋ Ｊ＝’ ＋ｄ、ｉ＝ｘ、 ２゜３．４ ａ ＋ １ ＝ｋＬ＋ ＸＲｔ ＋ 、 Ｘｆｔ量ヤ、
×土し＋＋ｓ 、 ｊ＝ｉ＋４、ｉ＝ｉ、 ２． ｓ ａ ＩＩ ＝ｆｌしｔ ＸＲｔ ＋ 、 ＸＪコＬＩ＋
ｔ ＸＲｔ ＋ Ｓ ×Ｌもｌ ＋４ ＋ ３ ＝’
＋５ ＋ｉ＝ｉ、２ ａｓ ？ ＝Ｒ１’ｌ ＸＲ，×几、×几、×几、ＸＩ
Ｌ。
値ａＩｊは、 ａｌｌ ＝＝Ｑ、 Ｏｊ＝ｊ ａｕ ＝Ｒｔ ｊ＝ｉ＋ｌ、 ｉ＝１．２．３．４．５
．６ａＢ−几ＩＸ几１．１ ＋ Ｊ−’＋２Ｉｉ＝Ｌ
２＋ ３．’ｓ ５ａＩＨ＝Ｒｔ ＸＲ，＋ ＋３ Ｘ
Ｒｔ ＋ｔ ＋ Ｊ＝’ ＋ｄ、ｉ＝ｘ、 ２゜３．４ ａ ＋ １ ＝ｋＬ＋ ＸＲｔ ＋ 、 Ｘｆｔ量ヤ、
×土し＋＋ｓ 、 ｊ＝ｉ＋４、ｉ＝ｉ、 ２． ｓ ａ ＩＩ ＝ｆｌしｔ ＸＲｔ ＋ 、 ＸＪコＬＩ＋
ｔ ＸＲｔ ＋ Ｓ ×Ｌもｌ ＋４ ＋ ３ ＝’
＋５ ＋ｉ＝ｉ、２ ａｓ ？ ＝Ｒ１’ｌ ＸＲ，×几、×几、×几、ＸＩ
Ｌ。
ａｌｌ ＝１．０／ａｎ
とする。
この行列ｋ１１を用いて、下から工番目と（ＩＩ１）番
目のノードにあるＬＰｉＬＭの出力値の変化率は、その
ＬＰＲＭを含むＬＰ几Ｍストリングの周囲の制御棒の先
端が下からｊ着目のノードからに番目のノードに移動し
た時、次の様に計算する。
目のノードにあるＬＰｉＬＭの出力値の変化率は、その
ＬＰＲＭを含むＬＰ几Ｍストリングの周囲の制御棒の先
端が下からｊ着目のノードからに番目のノードに移動し
た時、次の様に計算する。
変化率＝ Ａ Ｉ Ｊ
ただし １＝Ｉ−Ｊ＋４
ｊＩＩ −に＋ ４
ただし、ｉ２７の時ｉ＝７、Ｊ２７の時ｊ＝７第３図に
、一本の局所出力憤城モニタストリ／グ（ＬＰＲＭスト
リング）と、それを囲む１６本の制御棒の位置関係を示
す。
、一本の局所出力憤城モニタストリ／グ（ＬＰＲＭスト
リング）と、それを囲む１６本の制御棒の位置関係を示
す。
ＬＰＲＭストリングに最も近い制御棒２のバタン変化に
伴うＬＰＲＭ出力値の変化車積（前記のＲ８からＲ６）
に比べて、制御ａ２ａのバタン変化に伴うＬＰＲ，Ｍ出
力値の変化率１１ｘは、小さくとることにし、 Ｒ＋’ −（Ｒ＋ １．０ ） ＸＣＥ＋１．０により
足める。ただし、０≦α（ｉ、 。
伴うＬＰＲＭ出力値の変化車積（前記のＲ８からＲ６）
に比べて、制御ａ２ａのバタン変化に伴うＬＰＲ，Ｍ出
力値の変化率１１ｘは、小さくとることにし、 Ｒ＋’ −（Ｒ＋ １．０ ） ＸＣＥ＋１．０により
足める。ただし、０≦α（ｉ、 。
更に、？ｌ１ＩＩ御＠２ｂのバタン変化に伴う変化車積
″も Ｂ（’ ＝ （ｆＬ ｒ ｌ、Ｕ ） Ｘβ十ｉ、。
″も Ｂ（’ ＝ （ｆＬ ｒ ｌ、Ｕ ） Ｘβ十ｉ、。
により定め、β〈αとして距離による効果を反映させる
。
。
以上で説明したＬＰ几Ｍ出カ値の変化率をめる行列ＡＨ
を用いて、制御バタン変化後のＬｌ（Ｍ出力値は、次の
手順でめる。
を用いて、制御バタン変化後のＬｌ（Ｍ出力値は、次の
手順でめる。
（１） お心制御棒に注目し１本毎に、バタンの変化を
制＠欅先端ノードの位Ｉｔ変化としてめる。
制＠欅先端ノードの位Ｉｔ変化としてめる。
（２）制御棒に最も近いＬ Ｐ ＲＭストリング、次に
近いＬＰＲＭストリング、次に近いＬＰＲＭストリング
の谷ＬＰＲＭ尚さでの変化率を行列ＡＩＪによ請求める
。
近いＬＰＲＭストリング、次に近いＬＰＲＭストリング
の谷ＬＰＲＭ尚さでの変化率を行列ＡＩＪによ請求める
。
（３）求められた、変化率を対応するＬＰｋＬＭ出方１
直にかける。これにより注目し７辷制御褌のバタ／（引
抜）ツチ赦）変化による、周囲１．ＰＲＭ出力値への効
果が計算される。
直にかける。これにより注目し７辷制御褌のバタ／（引
抜）ツチ赦）変化による、周囲１．ＰＲＭ出力値への効
果が計算される。
（４）注目する制御棒をかえる。すべての制御棒につい
て、計算終了したら終了。
て、計算終了したら終了。
第４図に、本発明による原子炉炉心シミュレーション方
式で模擬した、ＬＰＢ、Ｍ出力値の制御棒引抜に伴う推
移を示す。ＬＰＲＭストリングは、炉心中央付近とし、
２５％出レ出歩ベル５０％出力レベルへ出力上昇する過
程を示す。制御棒を、このＬＰ几ＭストＩＪングのまわ
りで、第５図に示した制御棒引抜手順で引抜いた時、Ｌ
Ｐ几Ｍ出力値は、上のＬＰＲＭからａ、ｂ、ｃ、ｄで示
す様に変化する。なお、ｅは出力レベルをｌＯで割った
ものでＬＰＲＭ出力値の平均値でその変化は、ＬＰ几Ｍ
出力１直の増加率となる。
式で模擬した、ＬＰＢ、Ｍ出力値の制御棒引抜に伴う推
移を示す。ＬＰＲＭストリングは、炉心中央付近とし、
２５％出レ出歩ベル５０％出力レベルへ出力上昇する過
程を示す。制御棒を、このＬＰ几ＭストＩＪングのまわ
りで、第５図に示した制御棒引抜手順で引抜いた時、Ｌ
Ｐ几Ｍ出力値は、上のＬＰＲＭからａ、ｂ、ｃ、ｄで示
す様に変化する。なお、ｅは出力レベルをｌＯで割った
ものでＬＰＲＭ出力値の平均値でその変化は、ＬＰ几Ｍ
出力１直の増加率となる。
以上、本実施例によれば、制御棒操作後のＬＰ几Ｍ出力
値を、あらかじめたくわえた代表的出力レベルでの制御
棒パタ／としＬＰｉ（Ｍ出力値（相対１［）を参照し、
出ノ月直の変化率をめる行列をよみとシ基準となるＬＰ
ＲＭ出力値にかけることで計算することができる。この
計算法は、＜）かえし行列算法、特殊関数、級数等の反
復算出法を用層いていないので、上記の処理は、大型計
ＪＥｄを用いると関数演算時ｆｆ１ｌ程度の演算処理時
間で行なうことができる効果がある。
値を、あらかじめたくわえた代表的出力レベルでの制御
棒パタ／としＬＰｉ（Ｍ出力値（相対１［）を参照し、
出ノ月直の変化率をめる行列をよみとシ基準となるＬＰ
ＲＭ出力値にかけることで計算することができる。この
計算法は、＜）かえし行列算法、特殊関数、級数等の反
復算出法を用層いていないので、上記の処理は、大型計
ＪＥｄを用いると関数演算時ｆｆ１ｌ程度の演算処理時
間で行なうことができる効果がある。
〔発明の効果」
本発明によれば、原子炉炉心シュミレーションにおいて
、原子炉の起動・出力運転・停止の各状態Ｖご２いて圧
怠の制御棒パタンを入力として、′速ヤかに、ＬＰＲＭ
出刃値を模擬した計算結果が得られるという効果がある
。
、原子炉の起動・出力運転・停止の各状態Ｖご２いて圧
怠の制御棒パタンを入力として、′速ヤかに、ＬＰＲＭ
出刃値を模擬した計算結果が得られるという効果がある
。
第１図は本発明による原子炉炉心シミュレーション方式
の処理手順の一実施例を示す図、第２図は局所出力領域
モニタと、その出力値の変化率を与える制鍔揮位置との
位置関係を示す図、第３図は局所出力領域モニタストリ
ングを含む制御棒・燃料チャネルの位置関係を示す図、
第４図及第５図は本発明による原子炉炉心シミュレーシ
ョン方式により模擬された局所出力領域モニタ出力値の
制御棒引抜にともなう推移の例を示す図である。 ｌ・・・局所出力領域モニタストリング、２・・・制御
棒、３・・・燃料チャネル。 代理人 弁理士 高橋明夫 耐 ｌ 図 第 ｚｍ 稟 ３ 図 別寸卯棒引′Ｌ鋭・ノッ＋氏
の処理手順の一実施例を示す図、第２図は局所出力領域
モニタと、その出力値の変化率を与える制鍔揮位置との
位置関係を示す図、第３図は局所出力領域モニタストリ
ングを含む制御棒・燃料チャネルの位置関係を示す図、
第４図及第５図は本発明による原子炉炉心シミュレーシ
ョン方式により模擬された局所出力領域モニタ出力値の
制御棒引抜にともなう推移の例を示す図である。 ｌ・・・局所出力領域モニタストリング、２・・・制御
棒、３・・・燃料チャネル。 代理人 弁理士 高橋明夫 耐 ｌ 図 第 ｚｍ 稟 ３ 図 別寸卯棒引′Ｌ鋭・ノッ＋氏
Claims (1)
- 【特許請求の範囲】 １、実際の原子炉を模擬し、起動・出力運転・停止の状
態に応じた原子炉内の局所出力領域モニタの出力＋１ｔ
を演算し、その数１ｉｉ’に原子炉内中性子束分布の模
擬測定結果として得る、原子炉炉心シミュレーション方
式に２いて、代表的出力レベルでの制御棒パタン、引抜
線ノツチ数と局部出力領域モニタからの出力値ｔデータ
として電子計算機内の記憶装置に記憶しておき、現状の
制御棒パタンを入力すると、最も近い制御棒パタンの出
力レベルを現状の出力レベルの上下２点選び、各出力レ
ベルに対応する制御棒バタンと現状の制御棒パタ／を比
較し、制御棒の単位長引抜令そう人による局所出力領域
モニタの出力値の変化率を用いて下の出力レベルでの局
所出力領域モニタ出力値から計算された現状の局所出力
領域モニタ出力値の推定値と、同様に上の出力レベルで
の局所出力領域モニタ出力値から計算された現状の局所
出力機械モニタ出力値の推に値の両者の平均【、上の出
力レベルでの制御棒線ノツチ数と現状のノツチ数の差と
下の出力レベルでの制御棒線ノツチ数と現状のノツチ数
の差を各々前者、後者の推定呟の重みにしてＩｔｔ算し
た値をもって、現状の制御棒／くタンでの局所出力領域
モニタ出力値とすゐことを特徴とする原子炉炉心シミュ
レーション方式。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58194156A JPS6086489A (ja) | 1983-10-19 | 1983-10-19 | 原子炉炉心シミユレ−シヨン方式 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58194156A JPS6086489A (ja) | 1983-10-19 | 1983-10-19 | 原子炉炉心シミユレ−シヨン方式 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPS6086489A true JPS6086489A (ja) | 1985-05-16 |
Family
ID=16319849
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP58194156A Pending JPS6086489A (ja) | 1983-10-19 | 1983-10-19 | 原子炉炉心シミユレ−シヨン方式 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS6086489A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8739345B2 (en) | 2009-03-31 | 2014-06-03 | Industrial Brushware Limited | Down hole cleaning tool |
-
1983
- 1983-10-19 JP JP58194156A patent/JPS6086489A/ja active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8739345B2 (en) | 2009-03-31 | 2014-06-03 | Industrial Brushware Limited | Down hole cleaning tool |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110402467A (zh) | 用于为核反应堆建模的系统和方法 | |
| US20060149514A1 (en) | Method of determining a fresh fuel bundle design for a core of a nuclear reactor | |
| US20020122521A1 (en) | Determination of operating limit minimum critical power ratio | |
| US8762115B2 (en) | System and method for evaluating nuclear reactor fueling plan | |
| CN110737870B (zh) | 一种用于在gpu上合并舒尔矩阵的方法及装置 | |
| JPS6086489A (ja) | 原子炉炉心シミユレ−シヨン方式 | |
| Kang et al. | Coupling Gothic and RELAP5 for Passive Containment Cooling Modeling | |
| Martínez-Quiroga et al. | NIRK3D and 3DKIN: General Description and current status of the new 3D kinetics capabilities of RELAP5/SCDAPSIM/MOD4. 0 | |
| Arafa et al. | Development of an educational nuclear research reactor simulator | |
| JP2023036106A (ja) | 運転訓練装置、及び、運転訓練方法 | |
| JP2557462B2 (ja) | 原子炉の燃料配置作成装置 | |
| Galloway | Boiling Water Reactor Core Simulation with Generalized Isotopic Inventory Tracking for Actinide Management | |
| JP2696049B2 (ja) | 原子炉炉心核特性模擬装置 | |
| JP7470659B2 (ja) | 評価方法、評価システム及びプログラム | |
| JPS6151208A (ja) | 予測シミユレ−タ | |
| Wozny et al. | Hybrid computer simulation of xenon spatial oscillations | |
| CN115544732A (zh) | 适用于轻水堆的衰变热等效简化方法、存储介质及模拟机 | |
| Zoino et al. | Simulation tools and approaches for the compliance with performance-based ECCS cladding acceptance criteria (10 CFR 50.46 C) | |
| JPS63195593A (ja) | 炉心性能計算装置 | |
| Hall | A simulation model to aid in the design and tuning of hierarchical databases | |
| Boucheron et al. | MELCOR analysis of the Three Mile Island Unit 2 accident | |
| Aybar et al. | Dynamic probabilistic reliability and safety assessment: an application to IS-BWR | |
| Armstrong | Evaluating the Feasibility of MCNP/FDTD EMP Calculations at Altitudes above 20km | |
| Kobayash et al. | Optimization of refueling scheme of fast reactor core by means of man machine interactive method | |
| CN117971375A (zh) | 干涉服务分离的核电厂仿真系统 |