JP4323901B2 - 沸騰水型原子炉のボイド率 - Google Patents
沸騰水型原子炉のボイド率 Download PDFInfo
- Publication number
- JP4323901B2 JP4323901B2 JP2003314910A JP2003314910A JP4323901B2 JP 4323901 B2 JP4323901 B2 JP 4323901B2 JP 2003314910 A JP2003314910 A JP 2003314910A JP 2003314910 A JP2003314910 A JP 2003314910A JP 4323901 B2 JP4323901 B2 JP 4323901B2
- Authority
- JP
- Japan
- Prior art keywords
- steam
- saturated
- water
- vlh
- velocity
- 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.)
- Expired - Fee Related
Links
Images
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)
Description
本発明は特に沸騰水型原子炉のボイド率の相関式に関する。
沸騰水型原子炉は、原子炉内の核燃料で発生した熱を冷却材である水に伝達し原子炉内で液体の水を沸騰させ蒸気を発生する。蒸気はタービンに導かれ電気を発生する。沸騰水型原子炉を中核とした沸騰水型原子力発電所は、火力発電所同様に多くの機器系統から構成されている。運転には多くの計測、制御装置が必要である。大量のデータを効率的に取り扱う計算機が広く利用されている。
図１は沸騰水型原子力発電所における従来の計算機の主要な性能計算項目である。計測値を元にして、プラント性能計算と炉心性能計算を実施し表示、記録する。
プラント性能計算は、火力発電所で計算されているものと類似している。給水加熱器性能等の原子力発電所を構成している諸機器の性能、熱収支等の計算を実施し表示、記録する。
炉心性能計算は、原子炉運転を安全に管理するために原子炉熱出力等を計算し表示、記録する。核燃料であるウランやプルトニウムの中性子との反応は中性子の速度によって変わり、中性子の速度は減速材である水の量により変わる。沸騰水型原子炉では図２に示すように核燃料を内包する核燃料棒（１）の下部から流入した液体の水（１１）は核燃料棒（１）の上部の冷却材出口に行くにつれ核燃料から熱を受けて液体の一部は蒸気（１２）になる。液体である水（１１）と気体である蒸気（１２）が共存して流れている二相流の流路断面における単位高さ当たりの（蒸気が占める面積）÷（蒸気と液体である水が占める流路断面積）をボイド率と呼びαと表されている。ボイド率αは核燃料棒（１）の下部ではゼロであり核燃料棒（１）の上部では７０％にもなり、高さｚが増加するにつれ大きくなるからボイド率分布α(z)を持っている。ボイド率の違いは水の量の違いであるから核燃料と中性子との反応の結果生じる出力に変化をもたらし、逆に出力の変化はボイド率に違いをもたらす。冷却材流量が多く冷却材出口側も水だけであれば、上下端からの中性子漏洩により出力Pの高さｚに関する出力分布P(z)は中央が最大になる上下対称なコサイン分布状である。しかし、冷却材入口側である核燃料棒（１）の下部ではボイド率がゼロであるため核燃料の反応が活発であるが、冷却材出口側である核燃料棒（１）の上部ではボイド率が高いため核燃料の反応が不活発である。したがって、出力分布P(z)の最大は核燃料棒（１）の中央から下方に移動する。核燃料棒最高温度は当然核燃料棒（１）の下方で生じる。
このような挙動をする原子炉の運転では、核燃料棒最高温度等が設定値と比べて余裕があることを確認しながら運転する必要がある。また経済的な原子炉を設計するには、核燃料棒最高温度等の安全性にかかわる値が設定値内で、炉心の大きさをできるだけ小さくすることが必要である。そのためには、炉心内出力分布を平坦化して、最大出力と平均出力の比を小さくするようにする必要がある。
このように炉心全体でボイド率に分布の存在する沸騰水型原子炉では、ボイド率分布α(z)と出力分布P(z)の計算は取分け重要である。
図3は従来の計算機に内蔵されている従来のボイド率式を組み込んだソフトウエアを示すフローチャートである。
ステップ１
二相流の全質量流量GTは冷却材循環ポンプ性能等により決まる一定の値とし、高さゼロからの総熱出力Qを一定の値とする。出力分布はどの高さｚでも一定と仮定する。
ステップ２
高さzでの出力分布P0(z)が求まる。
ステップ３
高さzでの熱出力分布ｑ(z)=Q×P0(z) が求まる。
ステップ４
二相流の流路断面に含まれる蒸気流量を表す蒸気含有率χの高さｚでの蒸気含有率分布χ(z)は、水が下部から高さｚ迄流れる間に核燃料から得た全熱量にほぼ比例するからGTと熱出力分布ｑ(z)とで求まる。流路断面を通る単位時間当たり高さzでの質量流量において、飽和蒸気の質量流量Gg(z)、飽和水の質量流量Gl(z)から、流路断面を通る単位時間当たりの二相流の全質量流量GTはGT＝Gg(z)＋Gl(z)であるから、蒸気含有率はχ(z)＝Gg(z)/GTと表される。
ステップ５
χ(z)に対する従来のボイド率分布α(z)の相関式。飽和蒸気の速度Vgと飽和水の速度Vlとに差があるため、すべり比R＝Vg/Vlを実験により定め、この定数を使って、
α(z)を計算する。
ステップ６
高さｚでの核燃料の組成が決まっていると、水の組成はボイド率分布α(z)により決められるため、高さｚ位置毎の物質反応断面積が求まる。
ステップ7
拡散方程式を解いて出力分布P1(z)を計算する。
ステップ8
新たに計算された出力分布P1(z)と前回の出力分布P0(z)を比較して所定の誤差d1内でなければ、このP1(z)をP0(z)と置き換えてχ(z)を計算し再度出力分布P1(z)を計算する。新たに計算された出力分布P1(z)と前回の出力分布P0(z)を比較して所定の誤差d1内になったらそれが最終出力分布であるとする（非特許文献１）。
ステップ9
核燃料棒（１）の健全性の指標となる核燃料棒（１）の温度や限界出力比等を計算表示する。核燃料棒（１）の温度は熱伝達率から計算される。
核燃料棒（１）を多数本結束しチャンネルボックスで包囲してなる核燃料集合体（特許文献１）が多数ある場合は、ボイド率と抵抗係数との関係式等から核燃料集合体毎の圧力降下を計算し核燃料集合体毎の全質量流量を配分する繰り返し計算する部分がある。
図１は沸騰水型原子力発電所における従来の計算機の主要な性能計算項目である。計測値を元にして、プラント性能計算と炉心性能計算を実施し表示、記録する。
プラント性能計算は、火力発電所で計算されているものと類似している。給水加熱器性能等の原子力発電所を構成している諸機器の性能、熱収支等の計算を実施し表示、記録する。
炉心性能計算は、原子炉運転を安全に管理するために原子炉熱出力等を計算し表示、記録する。核燃料であるウランやプルトニウムの中性子との反応は中性子の速度によって変わり、中性子の速度は減速材である水の量により変わる。沸騰水型原子炉では図２に示すように核燃料を内包する核燃料棒（１）の下部から流入した液体の水（１１）は核燃料棒（１）の上部の冷却材出口に行くにつれ核燃料から熱を受けて液体の一部は蒸気（１２）になる。液体である水（１１）と気体である蒸気（１２）が共存して流れている二相流の流路断面における単位高さ当たりの（蒸気が占める面積）÷（蒸気と液体である水が占める流路断面積）をボイド率と呼びαと表されている。ボイド率αは核燃料棒（１）の下部ではゼロであり核燃料棒（１）の上部では７０％にもなり、高さｚが増加するにつれ大きくなるからボイド率分布α(z)を持っている。ボイド率の違いは水の量の違いであるから核燃料と中性子との反応の結果生じる出力に変化をもたらし、逆に出力の変化はボイド率に違いをもたらす。冷却材流量が多く冷却材出口側も水だけであれば、上下端からの中性子漏洩により出力Pの高さｚに関する出力分布P(z)は中央が最大になる上下対称なコサイン分布状である。しかし、冷却材入口側である核燃料棒（１）の下部ではボイド率がゼロであるため核燃料の反応が活発であるが、冷却材出口側である核燃料棒（１）の上部ではボイド率が高いため核燃料の反応が不活発である。したがって、出力分布P(z)の最大は核燃料棒（１）の中央から下方に移動する。核燃料棒最高温度は当然核燃料棒（１）の下方で生じる。
このような挙動をする原子炉の運転では、核燃料棒最高温度等が設定値と比べて余裕があることを確認しながら運転する必要がある。また経済的な原子炉を設計するには、核燃料棒最高温度等の安全性にかかわる値が設定値内で、炉心の大きさをできるだけ小さくすることが必要である。そのためには、炉心内出力分布を平坦化して、最大出力と平均出力の比を小さくするようにする必要がある。
このように炉心全体でボイド率に分布の存在する沸騰水型原子炉では、ボイド率分布α(z)と出力分布P(z)の計算は取分け重要である。
図3は従来の計算機に内蔵されている従来のボイド率式を組み込んだソフトウエアを示すフローチャートである。
ステップ１
二相流の全質量流量GTは冷却材循環ポンプ性能等により決まる一定の値とし、高さゼロからの総熱出力Qを一定の値とする。出力分布はどの高さｚでも一定と仮定する。
ステップ２
高さzでの出力分布P0(z)が求まる。
ステップ３
高さzでの熱出力分布ｑ(z)=Q×P0(z) が求まる。
ステップ４
二相流の流路断面に含まれる蒸気流量を表す蒸気含有率χの高さｚでの蒸気含有率分布χ(z)は、水が下部から高さｚ迄流れる間に核燃料から得た全熱量にほぼ比例するからGTと熱出力分布ｑ(z)とで求まる。流路断面を通る単位時間当たり高さzでの質量流量において、飽和蒸気の質量流量Gg(z)、飽和水の質量流量Gl(z)から、流路断面を通る単位時間当たりの二相流の全質量流量GTはGT＝Gg(z)＋Gl(z)であるから、蒸気含有率はχ(z)＝Gg(z)/GTと表される。
ステップ５
χ(z)に対する従来のボイド率分布α(z)の相関式。飽和蒸気の速度Vgと飽和水の速度Vlとに差があるため、すべり比R＝Vg/Vlを実験により定め、この定数を使って、
α(z)を計算する。
ステップ６
高さｚでの核燃料の組成が決まっていると、水の組成はボイド率分布α(z)により決められるため、高さｚ位置毎の物質反応断面積が求まる。
ステップ7
拡散方程式を解いて出力分布P1(z)を計算する。
ステップ8
新たに計算された出力分布P1(z)と前回の出力分布P0(z)を比較して所定の誤差d1内でなければ、このP1(z)をP0(z)と置き換えてχ(z)を計算し再度出力分布P1(z)を計算する。新たに計算された出力分布P1(z)と前回の出力分布P0(z)を比較して所定の誤差d1内になったらそれが最終出力分布であるとする（非特許文献１）。
ステップ9
核燃料棒（１）の健全性の指標となる核燃料棒（１）の温度や限界出力比等を計算表示する。核燃料棒（１）の温度は熱伝達率から計算される。
核燃料棒（１）を多数本結束しチャンネルボックスで包囲してなる核燃料集合体（特許文献１）が多数ある場合は、ボイド率と抵抗係数との関係式等から核燃料集合体毎の圧力降下を計算し核燃料集合体毎の全質量流量を配分する繰り返し計算する部分がある。
ボイド率αは、蒸気含有率χと、飽和蒸気の速度Vgと飽和水の速度Vlの比であるすべり比R＝Vg/Vlと、飽和蒸気密度ρgと飽和水密度ρlの比であるM=ρg/ρlから、
α=1/( 1- R×M×(1 - 1/χ) ) (非特許文献２)
と表される。Mは運転圧力が一定であるから蒸気表から求められ、χはGTと当該位置までの総出力から求まる。しかし、すべり比Rは今までαとχに関する実験により定数を定めていた。したがって、運転状態が大きく変わりボイド率が大きくなった場合には新たな実験をしてすべり比Rを求めねば適切なボイド率αが求まらないし、式そのものが適用できるとは実証されていない。したがって、出力分布が正しく求まるとは限らない。
沸騰水型原子炉の最近の動向として、ウランとプルトニウムの混合酸化物MOX燃料を導入しようとしている。プルトニウムを有効に利用するためには中性子速度を落とさないで燃焼させるほうがよいとされている。そのために、ボイド率αが高くなるような運転を想定した原子炉が望ましいとされている。ボイド率αの適切な求め方が求められている。
α=1/( 1- R×M×(1 - 1/χ) ) (非特許文献２)
と表される。Mは運転圧力が一定であるから蒸気表から求められ、χはGTと当該位置までの総出力から求まる。しかし、すべり比Rは今までαとχに関する実験により定数を定めていた。したがって、運転状態が大きく変わりボイド率が大きくなった場合には新たな実験をしてすべり比Rを求めねば適切なボイド率αが求まらないし、式そのものが適用できるとは実証されていない。したがって、出力分布が正しく求まるとは限らない。
沸騰水型原子炉の最近の動向として、ウランとプルトニウムの混合酸化物MOX燃料を導入しようとしている。プルトニウムを有効に利用するためには中性子速度を落とさないで燃焼させるほうがよいとされている。そのために、ボイド率αが高くなるような運転を想定した原子炉が望ましいとされている。ボイド率αの適切な求め方が求められている。
上記課題に鑑み、蒸気含有率χに関するボイド率αの新たな相関式を作成した。
以下に、説明に使われる記号の一覧を示した。
ρl：飽和水密度。規定されている運転圧力と蒸気表から求まる。
ρg：飽和蒸気密度。規定されている運転圧力と蒸気表から求まる。
M：ρgとρlの比。M=ρg/ρl。
A：流路断面積。
GT：流路断面を通る単位時間当たりの二相流の全質量流量。ポンプ性能等で求まっている液体の水の流入質量流量に等しい。
Vlh：飽和水初速度。Vlh＝GT / (ρl×A）と求められる。
Vg100：二相流が蒸気１００％になった時の飽和蒸気の速度。Vg100＝Vlh / M。
Gg：飽和蒸気の質量流量。
Gl：飽和水の質量流量。
χ：蒸気含有率。χ＝Gg/GT。
α：ボイド率。流路断面を通る単位時間当たりの二相流における単位高さ当たり、（蒸気が占める面積）÷（蒸気と液体である水が占める流路断面積）。
Vg：飽和蒸気の速度。
Vl：飽和水の速度。
ｋ：飽和蒸気の速度Vgにおけるαに対するχの重み。最適値を50とする正の値。
記号一覧終了。
E1=( M-1 ) / ( k+1 ) とし、重みｋをパラメータとしてχが0から1.0の範囲で本発明のボイド率αを、
α＝ ( ( M - E1×χ×k ) - ( ( M - E1×χ×k )2 - 4×χ×E1 )1/2 ）/ ( 2×E1 )
とする。
飽和蒸気の速度Vgをαとχ両方に依存するとして、
Vg＝Vlh＋（（α＋ｋ×χ）/（ｋ＋１））×（Vg100 - Vlh）
とする。
飽和水の速度Vlを、
Vl = ( Vlh - M×Vg×α ) / ( 1 - α )
とする。αが１では無限大になるがαが0.9999程度であれば有限な値になる。
飽和蒸気の質量流量Ggと飽和水の質量流量Glを、
Gg = ρg×Vg×α×Ａ
Gl = ρl×Vl×（1-α）×A
とする。
図４は計算機に内蔵せる、本発明のボイド率式を組み込んだ本発明のソフトウエアを示すフローチャートである。図３に示す従来との違いは、計算ステップ５では以下のようにした。
E1=( M-1 ) / ( k+1 ) とし、重みｋをパラメータとしてχが0から1.0の範囲で蒸気含有率分布χ(z)から、二相流のボイド率分布α(z)を、
α(z)＝( ( M - E1×χ(z)×k ) - ( ( M - E1×χ(z)×k )2 - 4×χ(z)×E1 )1/2 ）/ ( 2×E1 )
とした。
飽和蒸気の速度分布Vg(z)を、
Vg(z)＝Vlh＋（（α(z)＋ｋ×χ(z)）/（ｋ＋１））×（Vg100-Vlh）
とした。
ステップ９では以下のようにした。最終結果のα(z)とχ(z)から、
飽和蒸気の速度分布Vg(z)を、
Vg(z)＝Vlh＋（（α(z)＋ｋ×χ(z)）/（ｋ＋１））×（Vg100-Vlh）
飽和水の速度分布Vl(z)を、
Vl(z) = ( Vlh - M×Vg(z)×α(z) ) / ( 1 - α(z) )
とした。αが１では無限大になるがαが0.9999程度であれば有限な値になる。
飽和蒸気の質量流量分布Gg(z)と飽和水の質量流量分布Gl(z)を、
Gg(z) = ρg×Vg(z)×α(z)×Ａ
Gl = ρl×Vl(z)×（1-α(z) ）×A
とした。
核燃料棒（１）の健全性の指標となる核燃料棒（１）の温度や限界出力比等は従来通りに限界出力比を表示するか、飽和水の速度分布Vl(z)と( 1 - α(z) )から変換した水膜厚さとによる実験から限界熱流束に関わる式を定め最小限界熱流束比（MCHFR）を求め表示する。核燃料棒（１）の温度は熱伝達率から計算される（非特許文献３）。
核燃料集合体が多数ある場合は従来同様圧力降下を従来通りに、α(z)と抵抗係数との関係式等から求めるか、飽和水の速度分布Vl(z)から水のみの単相流から求めることにより各核燃料集合体毎の全質量流量GTも繰り返し計算により求める。
図5は、蒸気含有率χに対する本発明のボイド率αの相関式による、重みｋをパラメータとしたボイド率αの計算結果である。k=1000は従来のボイド率の相関式ですべり比が1.0に相当する。k=50は従来の相関式ですべり比が2.0に相当する（非特許文献４）。k=0はα増加による蒸気流路断面積の増加による蒸気流動抵抗の減少により蒸気流速が速くなると想定した場合の計算結果である。
図６，７，８は上記本発明の計算式による、運転圧力が70気圧、飽和水初速度Vlhが1m/sの場合の、重みｋによる蒸気含有率χ対飽和蒸気の速度Vg、飽和水の速度Vlの挙動を示す。縦軸速度の単位はm/sである。χが大きくなればVlもVgも速くなるため除熱効果は高いように見受けられる。
図９は、図1に示す沸騰水型原子力発電所における計測値を元にした従来の計算機での主要な性能計算項目（非特許文献５）のうち、炉心性能計算を図4に示した本発明のボイド率の計算式を組み込んだ本発明のソフトウエアを内蔵した本発明の計算機の主要な性能計算項目である。
：昭61-37591、「核燃料集合体」
：理工学社、著者三神他「核燃料管理の方法と解析」２３１頁。
：著者三神他「核燃料管理の方法と解析」２２０頁。
：日本機械学会「伝熱工学資料」。
：コロナ社、赤川著「気液二相流」52頁。横軸は対数目盛。
：オーム社、浅田他監修「原子力ハンドブック」215頁。
以下に、説明に使われる記号の一覧を示した。
ρl：飽和水密度。規定されている運転圧力と蒸気表から求まる。
ρg：飽和蒸気密度。規定されている運転圧力と蒸気表から求まる。
M：ρgとρlの比。M=ρg/ρl。
A：流路断面積。
GT：流路断面を通る単位時間当たりの二相流の全質量流量。ポンプ性能等で求まっている液体の水の流入質量流量に等しい。
Vlh：飽和水初速度。Vlh＝GT / (ρl×A）と求められる。
Vg100：二相流が蒸気１００％になった時の飽和蒸気の速度。Vg100＝Vlh / M。
Gg：飽和蒸気の質量流量。
Gl：飽和水の質量流量。
χ：蒸気含有率。χ＝Gg/GT。
α：ボイド率。流路断面を通る単位時間当たりの二相流における単位高さ当たり、（蒸気が占める面積）÷（蒸気と液体である水が占める流路断面積）。
Vg：飽和蒸気の速度。
Vl：飽和水の速度。
ｋ：飽和蒸気の速度Vgにおけるαに対するχの重み。最適値を50とする正の値。
記号一覧終了。
E1=( M-1 ) / ( k+1 ) とし、重みｋをパラメータとしてχが0から1.0の範囲で本発明のボイド率αを、
α＝ ( ( M - E1×χ×k ) - ( ( M - E1×χ×k )2 - 4×χ×E1 )1/2 ）/ ( 2×E1 )
とする。
飽和蒸気の速度Vgをαとχ両方に依存するとして、
Vg＝Vlh＋（（α＋ｋ×χ）/（ｋ＋１））×（Vg100 - Vlh）
とする。
飽和水の速度Vlを、
Vl = ( Vlh - M×Vg×α ) / ( 1 - α )
とする。αが１では無限大になるがαが0.9999程度であれば有限な値になる。
飽和蒸気の質量流量Ggと飽和水の質量流量Glを、
Gg = ρg×Vg×α×Ａ
Gl = ρl×Vl×（1-α）×A
とする。
図４は計算機に内蔵せる、本発明のボイド率式を組み込んだ本発明のソフトウエアを示すフローチャートである。図３に示す従来との違いは、計算ステップ５では以下のようにした。
E1=( M-1 ) / ( k+1 ) とし、重みｋをパラメータとしてχが0から1.0の範囲で蒸気含有率分布χ(z)から、二相流のボイド率分布α(z)を、
α(z)＝( ( M - E1×χ(z)×k ) - ( ( M - E1×χ(z)×k )2 - 4×χ(z)×E1 )1/2 ）/ ( 2×E1 )
とした。
飽和蒸気の速度分布Vg(z)を、
Vg(z)＝Vlh＋（（α(z)＋ｋ×χ(z)）/（ｋ＋１））×（Vg100-Vlh）
とした。
ステップ９では以下のようにした。最終結果のα(z)とχ(z)から、
飽和蒸気の速度分布Vg(z)を、
Vg(z)＝Vlh＋（（α(z)＋ｋ×χ(z)）/（ｋ＋１））×（Vg100-Vlh）
飽和水の速度分布Vl(z)を、
Vl(z) = ( Vlh - M×Vg(z)×α(z) ) / ( 1 - α(z) )
とした。αが１では無限大になるがαが0.9999程度であれば有限な値になる。
飽和蒸気の質量流量分布Gg(z)と飽和水の質量流量分布Gl(z)を、
Gg(z) = ρg×Vg(z)×α(z)×Ａ
Gl = ρl×Vl(z)×（1-α(z) ）×A
とした。
核燃料棒（１）の健全性の指標となる核燃料棒（１）の温度や限界出力比等は従来通りに限界出力比を表示するか、飽和水の速度分布Vl(z)と( 1 - α(z) )から変換した水膜厚さとによる実験から限界熱流束に関わる式を定め最小限界熱流束比（MCHFR）を求め表示する。核燃料棒（１）の温度は熱伝達率から計算される（非特許文献３）。
核燃料集合体が多数ある場合は従来同様圧力降下を従来通りに、α(z)と抵抗係数との関係式等から求めるか、飽和水の速度分布Vl(z)から水のみの単相流から求めることにより各核燃料集合体毎の全質量流量GTも繰り返し計算により求める。
図5は、蒸気含有率χに対する本発明のボイド率αの相関式による、重みｋをパラメータとしたボイド率αの計算結果である。k=1000は従来のボイド率の相関式ですべり比が1.0に相当する。k=50は従来の相関式ですべり比が2.0に相当する（非特許文献４）。k=0はα増加による蒸気流路断面積の増加による蒸気流動抵抗の減少により蒸気流速が速くなると想定した場合の計算結果である。
図６，７，８は上記本発明の計算式による、運転圧力が70気圧、飽和水初速度Vlhが1m/sの場合の、重みｋによる蒸気含有率χ対飽和蒸気の速度Vg、飽和水の速度Vlの挙動を示す。縦軸速度の単位はm/sである。χが大きくなればVlもVgも速くなるため除熱効果は高いように見受けられる。
図９は、図1に示す沸騰水型原子力発電所における計測値を元にした従来の計算機での主要な性能計算項目（非特許文献５）のうち、炉心性能計算を図4に示した本発明のボイド率の計算式を組み込んだ本発明のソフトウエアを内蔵した本発明の計算機の主要な性能計算項目である。
近年、沸騰水型原子炉において、従来のウラン核燃料同様に、流路断面積が広く運転時の平均ボイド率が40%程度と低くして減速材たる水を十分にして、低速な熱中性子によりプルトニウムを燃焼させることが考えられている。しかし、プルトニウムは高速な高速中性子で燃焼させると効率が良い。そこで、流路断面積が狭く運転時の平均ボイド率を60%程度に高め減速材たる水が少ない低減速炉が検討されている。
従来のボイド率相関式は、原子炉実運転時と同程度の出口蒸気含有率χが0.2近傍での実験により約2.0と定めた、すべり比R＝Vg/Vlを用いていた。低減速炉のように平均ボイド率が60%にもなるとχが0.5以上にもなり新たに大掛かりな実験をしてすべり比Rを求めねばならない。その上Rを一定としても良いのかどうか不明であった。また、VgやVlは顕には表されていなかった。
本発明のボイド率相関式では、１００％蒸気での蒸気の理論的速度を考慮しているため、χが0から1.0の範囲で正確に、αは元より飽和蒸気の速度Vg、飽和水の速度Vl、飽和蒸気の質量流量Gg、飽和水の質量流量Glを求めることができる。その結果出力分布P(z)も精度良く求めることができるため、プルトニウムを効率よく燃焼させることができて、経済性向上に役立つ低減速炉を精度よく設計することができる。VgやVlの値は熱水力に関する種々の物理的解釈に役立ち、新たな原子炉設計や高度な原子炉運転に役立つ。
特に、Vlが顕に表されるため核燃料棒（１）の表面からの除熱が精度良く計算されるため核燃料棒（１）の健全性に関して精度良く検討できるため、低減速炉の経済性向上に役立つ。
原子炉の運転に関しても計算精度が向上するため、信頼性が増し原子炉運転員の負担が減るため人員減が可能となり管理コストが下がり、結果として発電コスト低下になる。
従来のボイド率相関式は、原子炉実運転時と同程度の出口蒸気含有率χが0.2近傍での実験により約2.0と定めた、すべり比R＝Vg/Vlを用いていた。低減速炉のように平均ボイド率が60%にもなるとχが0.5以上にもなり新たに大掛かりな実験をしてすべり比Rを求めねばならない。その上Rを一定としても良いのかどうか不明であった。また、VgやVlは顕には表されていなかった。
本発明のボイド率相関式では、１００％蒸気での蒸気の理論的速度を考慮しているため、χが0から1.0の範囲で正確に、αは元より飽和蒸気の速度Vg、飽和水の速度Vl、飽和蒸気の質量流量Gg、飽和水の質量流量Glを求めることができる。その結果出力分布P(z)も精度良く求めることができるため、プルトニウムを効率よく燃焼させることができて、経済性向上に役立つ低減速炉を精度よく設計することができる。VgやVlの値は熱水力に関する種々の物理的解釈に役立ち、新たな原子炉設計や高度な原子炉運転に役立つ。
特に、Vlが顕に表されるため核燃料棒（１）の表面からの除熱が精度良く計算されるため核燃料棒（１）の健全性に関して精度良く検討できるため、低減速炉の経済性向上に役立つ。
原子炉の運転に関しても計算精度が向上するため、信頼性が増し原子炉運転員の負担が減るため人員減が可能となり管理コストが下がり、結果として発電コスト低下になる。
以下に、説明に使われる記号の一覧を示した。
ρlin ：未飽和水密度。
ρl：飽和水密度。規定されている運転圧力と蒸気表から求まる。
ρg：飽和蒸気密度。規定されている運転圧力と蒸気表から求まる。
M：ρgとρlの比。M=ρg/ρl。
A：流路断面積。
GT：流路断面を通る単位時間当たりの二相流の全質量流量。ポンプ性能等で求まっている液体の水の流入質量流量に等しい。
Vlin：未飽和液体の水の流入速度。
Vlh：飽和水初速度。Vlh＝GT / (ρl×A）と求められる。
Vg100：二相流が蒸気１００％になった時の飽和蒸気の速度。Vg100＝Vlh / M。
GH ：流路断面を通る単位時間当たりの二相流が蒸気１００％になった時の飽和蒸気の質量流量。GH＝ρg×Vg100×A。
Gg：流路断面を通る単位時間当たりの飽和蒸気の質量流量。
Gl：流路断面を通る単位時間当たりの飽和水の質量流量。
χ：蒸気含有率。χ＝Gg/GT。
α：ボイド率。流路断面を通る単位時間当たりの二相流の流路断面における単位高さ当たり、（蒸気が占める面積）÷（蒸気と液体である水が占める流路断面積）。
Vg：飽和蒸気の速度。
Vl：飽和水の速度。
ｋ：飽和蒸気の速度Vgにおけるαに対するχの重み。
記号一覧終了。
規定されている運転圧力と蒸気表から飽和水密度ρl 、飽和蒸気密度ρgが求まる。
二相流の全質量流量GTはポンプ性能等で求まっている液体の水の流入質量流量に等しく流路断面積Aは一定の値で決まっているから、液体の水の流入速度をVlin とし液体の水が飽和温度まで上昇したときの飽和水初速度をVlhすると、
GT=ρlin×Vlin×A＝ρl×Vlh×Aと表せるから飽和水初速度は、
Vlh＝GT / (ρl×A）と求められる。
二相流が蒸気１００％になった時の飽和蒸気の速度をVg100その質量流量をGHとすると、
GH＝ρg×Vg100×Aと表せる。質量保存則からGH =GTであるから、
ρg×Vg100×A＝ρl×Vlh×AよりVg100＝Vlh×ρl/ρg＝Vlh / MであるからVg100は求められる。
飽和蒸気の質量流量GgはGg＝ρg×Vg×α×Ａであるから蒸気含有率χは、今までは、
χ＝Gg/GT＝ρg×Vg×α×Ａ/（ρg×Vg×α×Ａ＋ρl×Vl×（1-α）×A）から変形してα=1/( 1-R×M×(1-1/χ) )としていた。本発明では、
χ＝Gg/GT＝ρg×Vg×α×Ａ/（ρl×Vlh×A）＝（ρg/ρl）×α×（Vg/ Vlh）と表せるからχ＝M×α×（Vg/ Vlh）と表せる。
Vgは今まであらわには表されていなかったが第一の近似として、
Vg＝Vlh＋χ×（Vg100 - Vlh）＝Vlh＋χ×（Vlh / M -Vlh）とχに関する線形近似とした。蒸気は燃料棒（１）を覆う液の表面から蒸発する近似とした。核燃料棒(1）の下部で飽和蒸気の初速度を飽和水初速度と同じVg＝Vlhとし、核燃料棒(1）の上部でχが100％蒸気になるとVg＝Vg100になるとした。
Vg＝Vlh＋χ×（Vlh / M -Vlh）＝Vlh×（１-χ×(1-1/ M) ）と表せるから、
χ＝M×α×（Vg/ Vlh）＝M×α×Vlh×（１-χ×(1-1/ M) ）/ Vlh
χ＝M×α×（１-χ×(1-1/ M) ）。これを変形すると、
α=1/( 1 - M×( 1 - 1/χ ) )となる。これは、従来のボイド率に関する相関式において飽和蒸気の速度Vgと飽和水の速度Vlの比であるすべり比R＝Vg/Vl＝１としたことになる。核燃料棒(1）の下部でχがゼロの時α＝0で、χが1になるとα＝1となり、χが0.5になるとαは約0.96である。
第二の近似として、
Vg＝Vlh＋α×（Vg100-Vlh）＝Vlh×（１-α×(1-1/ M) ）とαに関する線形近似とした。蒸気は蒸気通路が広くなれば抵抗が少なくなり速度を増すとした近似とした。すると蒸気含有率χは、
χ＝M×α×（Vg/ Vlh）＝M×α×Vlh×（１-α×(1-1/ M) ）/ Vlh
χ＝M×α×（１-α×(1-1/ M) ）。これを変形すると、
α＝( M - ( M2 - 4×χ×( M - 1 ) )1/2 ）/ ( 2×( M - 1) )となる。これは核燃料棒(1）の下部でχがゼロの時α＝0で、χが1になるとM<1だからα＝1となり、χが0.5になるとαは約0.7である。
実際の飽和蒸気の速度Vgは、上記２つの近似を同時に含むと想定できるため重みをｋとしたαとχの線形近似とした。
Vg＝Vlh＋（（α＋ｋ×χ）/（ｋ＋１））×（Vg100-Vlh）これを変形すると、
Vg＝Vlh×( 1+ （（α＋ｋ×χ）/（ｋ＋１））×（ (1/M) -1 ）)
上記と同様に、
χ＝M×α×（Vg/ Vlh）だからE1=( M-1 ) / ( k+1 ) と置くと、
α＝ ( ( M - E1×χ×k ) - ( ( M - E1×χ×k )2 - 4×χ×E1 )1/2 ）/ ( 2×E1 )
となる。
飽和水の速度VlはGg+Gl=GTより、
ρg×Vg×α×Ａ＋ρl×Vl×（1-α）×A =ρl×Vlh×A だから、
Vl = ( Vlh - M×Vg×α ) / ( 1 - α ) と求められる。
したがって、飽和蒸気の質量流量Ggと飽和水の質量流量Glは、
Gg = ρg×Vg×α×Ａ
Gl = ρl×Vl×（1-α）×A
と求められ、和は常にGTになる。
図5に示す本発明のボイド率αの相関式においてk=50の場合は、従来のボイド率αの相関式において蒸気含有率χが0.3以下の運転時ボイド率が低い原子炉において比較的良く合うとされているすべり比が2.0の場合に相当する。本発明のボイド率αの相関式は、蒸気含有率が低い状態でも従来のボイド率αの相関式と同等の値を示し、高蒸気含有率でも精度の良い相関式であることは１００％蒸気になった時の飽和蒸気の速度Vgが理論値Vg100になることが保証されていることからわかる。
図６，７，８は本発明のボイド率αの相関式において同時に求められる飽和蒸気の速度と飽和水の速度のk=0,50,1000の場合の挙動である。蒸気含有率χが1.0の時飽和蒸気の速度Vgは理論値になり、χが0.3以上になっても精度が保証されている。ボイド率αの相関式が保証されているということは、出力分布のボイド率に依存する計算も保証されるということであるから、原子炉挙動を高精度に解析できることが期待できる。したがって、運転状態の監視に高い信頼が持てる。計算された出力分布と炉心内核計側装置との比較において差が生じた場合においてもkの値を修正することにより精度が向上する。
プルトニウムを沸騰水型原子炉に適用する場合、ボイド率αが高い状態を想定した低減速炉の設計が望ましく、設計計算の際、ボイド率αの高い時にも高精度である式を使うことが必須である。
本発明のボイド率αの相関式において同時に求められる飽和蒸気の速度と飽和水の速度は、核燃料棒（１）の温度を求める際重要な熱伝達率を検討するのに使われる。核燃料棒（１）表面は液体の水で覆われているとして熱伝達率が求められる。ボイド率が1.0に近い場合でも高速な飽和蒸気と飽和水により除熱できることが計算可能になったため高度な炉心設計が可能となる。
また、多数の核燃料集合体を持った原子炉において、流速が解れば流動抵抗が解り流量分布も正確に求めることができて、炉心設計計算精度の向上につながる。
本発明のボイド率αの相関式中のパラメータkは、予め実験により定めることができる。即ち、想定した核燃料棒（１）の材料、配列形状、寸法を模擬した熱流動実験
によりkを定めることができる。
ρlin ：未飽和水密度。
ρl：飽和水密度。規定されている運転圧力と蒸気表から求まる。
ρg：飽和蒸気密度。規定されている運転圧力と蒸気表から求まる。
M：ρgとρlの比。M=ρg/ρl。
A：流路断面積。
GT：流路断面を通る単位時間当たりの二相流の全質量流量。ポンプ性能等で求まっている液体の水の流入質量流量に等しい。
Vlin：未飽和液体の水の流入速度。
Vlh：飽和水初速度。Vlh＝GT / (ρl×A）と求められる。
Vg100：二相流が蒸気１００％になった時の飽和蒸気の速度。Vg100＝Vlh / M。
GH ：流路断面を通る単位時間当たりの二相流が蒸気１００％になった時の飽和蒸気の質量流量。GH＝ρg×Vg100×A。
Gg：流路断面を通る単位時間当たりの飽和蒸気の質量流量。
Gl：流路断面を通る単位時間当たりの飽和水の質量流量。
χ：蒸気含有率。χ＝Gg/GT。
α：ボイド率。流路断面を通る単位時間当たりの二相流の流路断面における単位高さ当たり、（蒸気が占める面積）÷（蒸気と液体である水が占める流路断面積）。
Vg：飽和蒸気の速度。
Vl：飽和水の速度。
ｋ：飽和蒸気の速度Vgにおけるαに対するχの重み。
記号一覧終了。
規定されている運転圧力と蒸気表から飽和水密度ρl 、飽和蒸気密度ρgが求まる。
二相流の全質量流量GTはポンプ性能等で求まっている液体の水の流入質量流量に等しく流路断面積Aは一定の値で決まっているから、液体の水の流入速度をVlin とし液体の水が飽和温度まで上昇したときの飽和水初速度をVlhすると、
GT=ρlin×Vlin×A＝ρl×Vlh×Aと表せるから飽和水初速度は、
Vlh＝GT / (ρl×A）と求められる。
二相流が蒸気１００％になった時の飽和蒸気の速度をVg100その質量流量をGHとすると、
GH＝ρg×Vg100×Aと表せる。質量保存則からGH =GTであるから、
ρg×Vg100×A＝ρl×Vlh×AよりVg100＝Vlh×ρl/ρg＝Vlh / MであるからVg100は求められる。
飽和蒸気の質量流量GgはGg＝ρg×Vg×α×Ａであるから蒸気含有率χは、今までは、
χ＝Gg/GT＝ρg×Vg×α×Ａ/（ρg×Vg×α×Ａ＋ρl×Vl×（1-α）×A）から変形してα=1/( 1-R×M×(1-1/χ) )としていた。本発明では、
χ＝Gg/GT＝ρg×Vg×α×Ａ/（ρl×Vlh×A）＝（ρg/ρl）×α×（Vg/ Vlh）と表せるからχ＝M×α×（Vg/ Vlh）と表せる。
Vgは今まであらわには表されていなかったが第一の近似として、
Vg＝Vlh＋χ×（Vg100 - Vlh）＝Vlh＋χ×（Vlh / M -Vlh）とχに関する線形近似とした。蒸気は燃料棒（１）を覆う液の表面から蒸発する近似とした。核燃料棒(1）の下部で飽和蒸気の初速度を飽和水初速度と同じVg＝Vlhとし、核燃料棒(1）の上部でχが100％蒸気になるとVg＝Vg100になるとした。
Vg＝Vlh＋χ×（Vlh / M -Vlh）＝Vlh×（１-χ×(1-1/ M) ）と表せるから、
χ＝M×α×（Vg/ Vlh）＝M×α×Vlh×（１-χ×(1-1/ M) ）/ Vlh
χ＝M×α×（１-χ×(1-1/ M) ）。これを変形すると、
α=1/( 1 - M×( 1 - 1/χ ) )となる。これは、従来のボイド率に関する相関式において飽和蒸気の速度Vgと飽和水の速度Vlの比であるすべり比R＝Vg/Vl＝１としたことになる。核燃料棒(1）の下部でχがゼロの時α＝0で、χが1になるとα＝1となり、χが0.5になるとαは約0.96である。
第二の近似として、
Vg＝Vlh＋α×（Vg100-Vlh）＝Vlh×（１-α×(1-1/ M) ）とαに関する線形近似とした。蒸気は蒸気通路が広くなれば抵抗が少なくなり速度を増すとした近似とした。すると蒸気含有率χは、
χ＝M×α×（Vg/ Vlh）＝M×α×Vlh×（１-α×(1-1/ M) ）/ Vlh
χ＝M×α×（１-α×(1-1/ M) ）。これを変形すると、
α＝( M - ( M2 - 4×χ×( M - 1 ) )1/2 ）/ ( 2×( M - 1) )となる。これは核燃料棒(1）の下部でχがゼロの時α＝0で、χが1になるとM<1だからα＝1となり、χが0.5になるとαは約0.7である。
実際の飽和蒸気の速度Vgは、上記２つの近似を同時に含むと想定できるため重みをｋとしたαとχの線形近似とした。
Vg＝Vlh＋（（α＋ｋ×χ）/（ｋ＋１））×（Vg100-Vlh）これを変形すると、
Vg＝Vlh×( 1+ （（α＋ｋ×χ）/（ｋ＋１））×（ (1/M) -1 ）)
上記と同様に、
χ＝M×α×（Vg/ Vlh）だからE1=( M-1 ) / ( k+1 ) と置くと、
α＝ ( ( M - E1×χ×k ) - ( ( M - E1×χ×k )2 - 4×χ×E1 )1/2 ）/ ( 2×E1 )
となる。
飽和水の速度VlはGg+Gl=GTより、
ρg×Vg×α×Ａ＋ρl×Vl×（1-α）×A =ρl×Vlh×A だから、
Vl = ( Vlh - M×Vg×α ) / ( 1 - α ) と求められる。
したがって、飽和蒸気の質量流量Ggと飽和水の質量流量Glは、
Gg = ρg×Vg×α×Ａ
Gl = ρl×Vl×（1-α）×A
と求められ、和は常にGTになる。
図5に示す本発明のボイド率αの相関式においてk=50の場合は、従来のボイド率αの相関式において蒸気含有率χが0.3以下の運転時ボイド率が低い原子炉において比較的良く合うとされているすべり比が2.0の場合に相当する。本発明のボイド率αの相関式は、蒸気含有率が低い状態でも従来のボイド率αの相関式と同等の値を示し、高蒸気含有率でも精度の良い相関式であることは１００％蒸気になった時の飽和蒸気の速度Vgが理論値Vg100になることが保証されていることからわかる。
図６，７，８は本発明のボイド率αの相関式において同時に求められる飽和蒸気の速度と飽和水の速度のk=0,50,1000の場合の挙動である。蒸気含有率χが1.0の時飽和蒸気の速度Vgは理論値になり、χが0.3以上になっても精度が保証されている。ボイド率αの相関式が保証されているということは、出力分布のボイド率に依存する計算も保証されるということであるから、原子炉挙動を高精度に解析できることが期待できる。したがって、運転状態の監視に高い信頼が持てる。計算された出力分布と炉心内核計側装置との比較において差が生じた場合においてもkの値を修正することにより精度が向上する。
プルトニウムを沸騰水型原子炉に適用する場合、ボイド率αが高い状態を想定した低減速炉の設計が望ましく、設計計算の際、ボイド率αの高い時にも高精度である式を使うことが必須である。
本発明のボイド率αの相関式において同時に求められる飽和蒸気の速度と飽和水の速度は、核燃料棒（１）の温度を求める際重要な熱伝達率を検討するのに使われる。核燃料棒（１）表面は液体の水で覆われているとして熱伝達率が求められる。ボイド率が1.0に近い場合でも高速な飽和蒸気と飽和水により除熱できることが計算可能になったため高度な炉心設計が可能となる。
また、多数の核燃料集合体を持った原子炉において、流速が解れば流動抵抗が解り流量分布も正確に求めることができて、炉心設計計算精度の向上につながる。
本発明のボイド率αの相関式中のパラメータkは、予め実験により定めることができる。即ち、想定した核燃料棒（１）の材料、配列形状、寸法を模擬した熱流動実験
によりkを定めることができる。
沸騰水型原子力発電所は日本では約２０基以上建設され稼動している。核燃料の利用が進みプルトニウムが大量に蓄積されている。無駄でもよいから無理やり消耗させようとするのではなく効率良く燃焼させるには、ボイド率が高い状態で運転する低減速原子炉が良いといわれている。沸騰水型原子炉は元々ボイドを発生させているため、ボイド率の程度の問題であるから検討が進んでいる。本発明のボイド率に関する式が適用されることにより、非常に高いボイド率を想定した原子炉を設計することができ、プルトニウムの効率的運用を通して発電コストの低減に寄与できる。
ボイド率の問題は沸騰水型原子炉だけでなく、加圧水型原子力発電所の蒸気発生器にもでてくる。また、火力発電所においてもボイラーで発生した熱を水管で蒸気にすることによりタービンを回転させるわけであるから、ボイド率の問題は重要である。
ボイド率の問題は沸騰水型原子炉だけでなく、加圧水型原子力発電所の蒸気発生器にもでてくる。また、火力発電所においてもボイラーで発生した熱を水管で蒸気にすることによりタービンを回転させるわけであるから、ボイド率の問題は重要である。
１ 沸騰水型原子炉の核燃料棒
１１ 沸騰水型原子炉の二相流の液体である水
１２ 沸騰水型原子炉の二相流の気体である蒸気
１１ 沸騰水型原子炉の二相流の液体である水
１２ 沸騰水型原子炉の二相流の気体である蒸気
Claims (3)
- 水と蒸気とからなる二相流において、以下の記号を使って、
ρl：飽和水密度。
ρg：飽和蒸気密度。
M：ρgとρlの比。M=ρg/ρl。
A：流路断面積。
GT：流路断面を通る単位時間当たりの二相流の全質量流量。
Vlh：流路断面を通る単位時間当たりの飽和水初速度。Vlh＝GT / (ρl×A）。
Vg100：流路断面を通る単位時間当たりの二相流が蒸気１００％になった時の飽和蒸気の速度。Vg100＝Vlh / M。
Gg：流路断面を通る単位時間当たりの飽和蒸気の質量流量。
Gl：流路断面を通る単位時間当たりの飽和水の質量流量。
χ：蒸気含有率。χ＝Gg/GT。
α：ボイド率。
Vg：流路断面を通る単位時間当たりの飽和蒸気の速度。
Vl：流路断面を通る単位時間当たりの飽和水の速度。
ｋ：飽和蒸気の速度Vgにおけるαに対するχの重み。
記号一覧終了。
χに対するボイド率αの相関式を、E1=( M-1 ) / ( k+1 )として、
α＝ ( ( M - E1×χ×k ) - ( ( M - E1×χ×k )2 - 4×χ×E1 )1/2 ）/ ( 2×E1 )
とし、
飽和蒸気の速度Vgを、
Vg＝Vlh＋（（α＋ｋ×χ）/（ｋ＋１））×（Vg100-Vlh）
とし、
飽和水の速度Vlを、
Vl = ( Vlh - M×Vg×α ) / ( 1 - α )
とし、
飽和蒸気の質量流量Ggと飽和水の質量流量Glを、
Gg = ρg×Vg×α×Ａ
Gl = ρl×Vl×（1-α）×A
として導出する方法。 - 沸騰水型原子力発電所における炉心性能計算において、水と蒸気とからなる二相流に関わる式として請求項１における式を組み込んだことを特徴とするソフトウエア及び当該ソフトウエアを内蔵する計算機。
- 水と蒸気とからなる二相流に関わる式として請求項１における式を使ったことを特徴とする沸騰水型原子炉の炉心設計方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003314910A JP4323901B2 (ja) | 2003-09-08 | 2003-09-08 | 沸騰水型原子炉のボイド率 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003314910A JP4323901B2 (ja) | 2003-09-08 | 2003-09-08 | 沸騰水型原子炉のボイド率 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2005083842A JP2005083842A (ja) | 2005-03-31 |
| JP4323901B2 true JP4323901B2 (ja) | 2009-09-02 |
Family
ID=34415325
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2003314910A Expired - Fee Related JP4323901B2 (ja) | 2003-09-08 | 2003-09-08 | 沸騰水型原子炉のボイド率 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4323901B2 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5661274B2 (ja) * | 2009-12-04 | 2015-01-28 | 白川 利久 | 沸騰水型原子炉の核燃料棒温度の導出方法 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS52144595A (en) * | 1976-05-28 | 1977-12-01 | Nippon Atom Ind Group Co Ltd | Measuring method and device of void distribution |
| JPS5413899A (en) * | 1977-07-04 | 1979-02-01 | Toshiba Corp | Nuclear fuel assembly |
| JPS63173921A (ja) * | 1987-01-14 | 1988-07-18 | Nippon Atom Ind Group Co Ltd | 蒸気クオリテイ測定装置 |
| JPH0713664B2 (ja) * | 1989-04-26 | 1995-02-15 | 株式会社日立製作所 | 燃料集合体および燃料スペーサ |
| JP2777004B2 (ja) * | 1991-11-11 | 1998-07-16 | 東京電力株式会社 | 原子炉内流量の測定方法およびその測定システム |
| US5600073A (en) * | 1994-11-02 | 1997-02-04 | Foster-Miller, Inc. | Method and system for analyzing a two phase flow |
| JP2866021B2 (ja) * | 1995-01-27 | 1999-03-08 | 石油公団 | タ−ビン型流量計による２相流体の流量計測方法 |
| JPH11264887A (ja) * | 1998-03-17 | 1999-09-28 | Toshiba Corp | 原子炉核計装システム、このシステムを備えた原子炉出力分布監視システムおよび原子炉出力分布監視方法 |
-
2003
- 2003-09-08 JP JP2003314910A patent/JP4323901B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2005083842A (ja) | 2005-03-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yeoh | Thermal hydraulic considerations of nuclear reactor systems: Past, present and future challenges | |
| KR20140063564A (ko) | 핵 원자로 노심 내에서 연료봉 파워 분포를 모델링하기 위한 방법론 | |
| JP4854654B2 (ja) | 炉心性能計算装置 | |
| Gui et al. | Experimental investigation on quenching behavior during reflooding in a 3× 3 dual-cooled annular rod bundle | |
| JP4323901B2 (ja) | 沸騰水型原子炉のボイド率 | |
| Moon et al. | An experimental study on the critical heat flux for low flow of water in a non-uniformly heated vertical rod bundle over a wide range of pressure conditions | |
| Sanchez-Torrijos et al. | A review on the thermo-mechanical modelling needs in system codes | |
| Zhu et al. | Subchannel analysis of critical heat flux in radially Non-uniform heated rod assembly at low flow rate | |
| Svanholm et al. | Halden reactors IFA-511.2 and IFA-54x: Experimental series under adverse core cooling conditions | |
| D'Auria | The BWR stability issue | |
| JP4707826B2 (ja) | 沸騰水型原子炉の監視制御装置 | |
| Jecong et al. | Safety Analysis of 250-kW Philippine Research Reactor-1 Thermal-hydraulics under Steady-state Operations Using MARS-KS Code. | |
| Guo et al. | SAC-CFR computer code verification with Experimental Breeder Reactor II loss-of-primary-flow-without-scram tests data | |
| JP5661274B2 (ja) | 沸騰水型原子炉の核燃料棒温度の導出方法 | |
| Ma et al. | Numerical study on thermal shock of central measuring shroud with fluid–solid thermal coupling in the “Monju” reactor | |
| Salko Jr et al. | Subchannel Methods Development for Modeling of Light Water Reactors at Oak Ridge National Laboratory | |
| Moon et al. | Experimental study on the low flow CHF in vertical 3x3 rod bundle with non-uniform axial heat flux distribution | |
| Carver et al. | INSTABILITIES IN FLOW THROUGH MULTIROD BUNDLES AND THEIR EFFECT ON DRYOUT. | |
| Oh et al. | Preliminary simulation of the KAERI rods irradiated in IFA-790 using FRAPCON4. 0P1 | |
| Mclntyre et al. | CONTROLLED TRANSIENT CHF TESTS IN A 5× 5 ROD BUNDLE UNDER LOSS-OF-COOLANT ACCIDENT CONDITIONS | |
| Ivankov | Analysis of PWR-PACTEL Small Break LOCA experiment using TRACE | |
| Liang et al. | Periodic remaining life evaluation program of PWR pressurizer surge line accounting for thermal stratification effect | |
| McCREERY et al. | Primary pump power as a measure of fluid density during bubbly two-phase flow.[PWR] | |
| Hwang et al. | Critical heat flux in square-and nonsquare-array rod bundles for advanced light water reactors | |
| Middleton et al. | FACTAR 2.0 code validation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060529 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20081217 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090421 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090424 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20090602 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20090605 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120612 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120612 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150612 Year of fee payment: 6 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |