(57)【要約】【課題】高温・高圧の沸騰水型原子炉内に発生する沸騰二相流の熱流動現象を常温・低圧の環境で再現し、沸騰二相流が燃料棒の伝熱表面に形成する液膜の挙動を常温・低圧の環境で模擬するとともに、液膜の膜厚又は液膜流量の測定を可能にする。【解決手段】沸騰二相流の密度、表面張力及び粘度の物性値と実質的に同一又は同等の密度、表面張力及び粘度の物性値を常温・低圧の環境で示す二相流体を所定断面且つ常温・低圧の直線流路に通して、燃料棒の伝熱表面の液膜を常温・低圧の環境で再現して液膜の挙動を模擬し、液膜の膜厚又は液膜流量を測定する。【選択図】図６
2007316059.tif 000002
【特許請求の範囲】
【請求項１】
沸騰水型原子炉内の高温・高圧環境に存在する沸騰二相流の熱流動現象を再現する熱流動現象の模擬方法において、
前記沸騰二相流の密度、表面張力及び粘度の物性値と実質的に同一又は同等の密度、表面張力及び粘度の物性値を常温・低圧の環境で示す二相流体の流れを液体及び気体の混合流体によって形成し、
前記二相流体を所定断面且つ常温・低圧の直線流路に通して、該流路の流路壁面に液膜を形成し、高温・高圧の環境で前記沸騰二相流が前記沸騰水型原子炉の燃料棒の伝熱表面に形成する液膜を常温・低圧の環境で流路壁面に再現して該液膜の挙動を模擬し、
前記液膜の膜厚及び／又は液膜流量を測定することを特徴とする熱流動現象の模擬方法。
【請求項２】
前記液体は、エタノールであり、前記気体は、HFC134a ガスであることを特徴とする請求項１に記載の熱流動現象の模擬方法。
【請求項３】
流路方向に間隔を隔てた複数の計測点において液膜の膜厚を電気抵抗法又は光学的測定法によって計測することを特徴とする請求項１又は２に記載の熱流動現象の模擬方法。
【請求項４】
請求項１乃至３のいずれか１項に記載の模擬方法によって限界出力を最適化した燃料集合体を有することを特徴とする沸騰水型原子炉。
【請求項５】
沸騰水型原子炉内の高温・高圧環境に存在する沸騰二相流の熱流動現象を再現する熱流動現象の模擬試験装置において、
前記沸騰二相流の密度、表面張力及び粘度の物性値と実質的に同一又は同等の密度、表面張力及び粘度の物性値を常温・低圧の環境で示す二相流体が供給される常温・低圧の直線流路と、
前記二相流体を形成する液体及び気体の混合流体を前記直線流路に供給する流体供給手段と、
前記混合流体によって前記直線流路の流路壁面に形成された液膜の膜厚及び／又は液膜流量を測定する液膜計測手段とを備え、
前記沸騰水型原子炉の燃料棒の伝熱表面に形成される液膜を常温・低圧の環境で前記流路壁面に再現して前記液膜の挙動を模擬するとともに、前記液膜の膜厚及び／又は液膜流量を測定するようにしたことを特徴とする熱流動現象の模擬試験装置。
【請求項６】
前記液膜計測手段は、流路方向に所定間隔を隔てた複数の計測点において液膜の膜厚を電気抵抗法又は光学的測定法によって計測する膜厚計測装置を備えることを特徴とする請求項５に記載の模擬試験装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、熱流動現象の模擬方法及び模擬試験装置に関するものであり、より詳細には、沸騰水型原子炉内の高温・高圧環境に存在する沸騰二相流の熱流動現象を常温・低圧の環境で再現する熱流動現象の模擬方法及び模擬試験装置に関するものである。
【背景技術】
【０００２】
沸騰水型原子炉（ＢＷＲ）の原子炉圧力容器内に配置される燃料集合体においては、燃料棒は、例えば、数mm程度の間隔で正方格子状に稠密に配置され、冷却材としての水が単相流の状態で燃料集合体の下部から流入し、発熱する燃料棒の間隙を流動し、燃料棒の発熱によって沸騰しながら燃料棒の間の間隙を上昇する。蒸気及び水の沸騰二相流が燃料集合体の内部に発生するので、炉心熱設計を行う場合、この沸騰二相流の挙動を詳細に把握する必要がある。
【０００３】
燃料棒の限界出力（限界熱流束）付近では、高速の蒸気流が燃料棒廻りに生成するので、極めて薄い液膜が環状流として燃料棒の外周面を覆い、熱伝達率が非常に高い状態が燃料棒の伝熱表面に生じる。更に出力が上昇すると、燃料棒表面の液膜が蒸発して燃料棒表面が乾く液膜ドライアウト現象が発生する。このようなドライアウトが一旦発生すると、熱伝達率が急激に悪化するので、燃料棒の壁温が急上昇し、場合によっては発熱管の焼損（バーンアウト）現象が発生し得る。このため、このような限界条件における燃料棒上の液膜挙動を正確に把握することは、原子炉の熱設計を行う上で極めて重要である。
【０００４】
従来、沸騰水型原子炉の開発・設計では、開発から実証に至る過程において、試験・実験的アプローチが重視されてきた。現状では、原子炉燃料の開発・設計は、実機条件（圧力7MPa、温度285 ℃）の実証試験（電気加熱）を行い、設計変更又は条件変更を繰り返して最適化を図るという手法が採用されている。このような開発・設計手法では、高温高圧試験を繰り返し実施しなければならないことから、莫大なコスト及び時間が必要とされる。例えば、原子炉燃料の設計を最適化するには、数週間を要する高額な実証試験を数十回も繰り返し実施して改良を図る必要が生じる。このような試験・実験的アプローチは、コスト及び時間的な不利を生じさせるのみならず、燃料の大型化や、格子稠密化等の近年の傾向と関連して、試験・実験機器の容量的な問題を生じさせており、試験・実験的アプローチを採用すること自体、困難になりつつある。このため、試験・実験に依存した開発・設計から解析中心の開発・設計への移行が望まれるとともに、簡易且つ効率的に試験・実験を行うことができる試験・実験方法及び試験・実験機器の開発が要望される。
【０００５】
熱流動現象の試験・実験を簡易且つ効率的に実施することを可能にする方法として、６フッ化硫黄ガスを液体アルコールに噴霧して気液二相流を形成し、高温・高圧環境の気液二相流を低温・低圧環境において模擬する模擬試験方法が特開２００２−１８９０９６号公報に開示されている。しかしながら、この模擬試験方法は、加圧水型原子炉（ＰＷＲ型原子炉）の蒸気発生器における気液二相流の流動現象を６フッ化硫黄ガス及び液体アルコールの気液二相流によって模擬し、液滴分離装置の液滴分離性能を評価することを意図したものであるにすぎない。
【特許文献１】特開２００２−１８９０９６号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
沸騰水型原子炉（ＢＷＲ型原子炉）の液膜ドライアウト現象は、気液二相流自体の挙動や物性に関連するのみならず、気液二相流の流動によって伝熱表面に生じる液膜の挙動と密接に関連する。即ち、沸騰水型原子炉の燃料集合体内を上昇する沸騰二相流は、前述の如く、燃料棒の表面(伝熱表面)に薄い液膜を形成するが、この液膜は、燃料棒の伝熱表面において非常に激しく変動し、複雑な液膜挙動を示すと考えられるので、液膜ドライアウトの発生条件等を正確に把握するには、このような液膜の挙動を模擬する必要がある。
【０００７】
他方、解析を中心として原子炉の開発・設計手法として、サブチャンネル解析コード(NASCA)が知られており、この解析手法を改良して原子炉の汎用的解析手法を確立するための研究が、多くの研究者によって実施されている。しかし、沸騰水型原子炉の炉内には、燃料棒の除熱性能に大きく影響を与える複数の因子が存在し、これらの因子は、解析手法の確立を困難にしている。例えば、沸騰水型原子炉の炉内には、燃料棒の振動を抑制し且つ相互接触を防止するためのスペーサが燃料集合体の適所に配置されており、燃料棒表面の液膜は、スペーサの近傍において複雑に膜厚を変化させる。このため、スペーサの形状及び位置は、燃料棒の限界出力を設定する上で極めて重要な因子であるが、スペーサが燃料棒表面の液膜挙動に与える影響は、極めて予測し難く、燃料棒近傍の熱流動現象に対するスペーサの影響を定量的に評価する手法は、依然として確立されていない。
【０００８】
本発明は、このような事情に鑑みてなされたものであり、その目的とするところは、高温・高圧の沸騰水型原子炉内に発生する沸騰二相流の熱流動現象を常温・低圧の環境で再現し、沸騰二相流が燃料棒の伝熱表面に形成する液膜の挙動を常温・低圧の環境で模擬するとともに、液膜挙動の測定を可能にし、これにより、燃料棒の除熱性能に大きく影響を与える因子の定量的評価を可能にする熱流動現象の模擬方法及び模擬試験装置を提供することにある。
【課題を解決するための手段】
【０００９】
本発明は、上記目的を達成すべく、沸騰水型原子炉内の高温・高圧環境に存在する沸騰二相流の熱流動現象を再現する熱流動現象の模擬方法において、
前記沸騰二相流の密度、表面張力及び粘度の物性値と実質的に同一又は同等の密度、表面張力及び粘度の物性値を常温・低圧の環境で示す二相流体の流れを液体及び気体の混合流体によって形成し、
前記二相流体を所定断面且つ常温・低圧の直線流路に通して、該流路の流路壁面に液膜を形成し、高温・高圧の環境で前記沸騰二相流が前記沸騰水型原子炉の燃料棒の伝熱表面に形成する液膜を常温・低圧の環境で流路壁面に再現して該液膜の挙動を模擬し、
前記液膜の膜厚及び／又は液膜流量を測定することを特徴とする熱流動現象の模擬方法を提供する。
【００１０】
本発明の上記構成によれば、沸騰水型原子炉内の沸騰二相流の流れを模擬する二相流体の流れが、直線流路内に形成される。燃料棒の伝熱表面の液膜挙動が常温・低圧の環境で流路壁面に模擬され、液膜の膜厚及び／又は液膜流量が常温・低圧の環境で測定される。このような模擬方法によれば、例えば、スペーサを模擬する障害物を流路内に配置した状態で液膜挙動を模擬することができ、これにより、スペーサが液膜挙動に与える影響を分析し、サブチャンネル解析コード(NASCA)に適用可能なスペーサ増倍係数等を設定することが可能となる。
【００１１】
なお、本明細書において、「常温」は、概ね室温又は大気温度に相当する温度範囲（概ね０℃から５０℃の範囲）の温度を意味し、「低圧」は、通常の実験機器によって形成可能な範囲（概ね０．１〜１MPaの範囲）の圧力を意味する。  
【００１２】
好ましくは、上記液体は、エタノールであり、上記気体は、HFC134a ガスである。更に好ましくは、流路方向に間隔を隔てた複数の計測点における液膜の膜厚が、電気抵抗法（例えば、定電流法）又は光学的な測定法によって計測される。
【００１３】
本発明は又、上記模擬方法によって限界出力を最適化した燃料集合体を有する沸騰水型原子炉を提供する。
【００１４】
本発明は更に、沸騰水型原子炉内の高温・高圧環境に存在する沸騰二相流の熱流動現象を再現する熱流動現象の模擬試験装置において、
前記沸騰二相流の密度、表面張力及び粘度の物性値と実質的に同一又は同等の密度、表面張力及び粘度の物性値を常温・低圧の環境で示す二相流体が供給される常温・低圧の直線流路と、
前記二相流体を形成する液体及び気体の混合流体を前記直線流路に供給する流体供給手段と、
前記混合流体によって前記直線流路の流路壁面に形成された液膜の膜厚及び／又は液膜流量を測定する液膜計測手段とを備え、
前記沸騰水型原子炉の燃料棒の伝熱表面に形成される液膜を常温・低圧の環境で前記流路壁面に再現して前記液膜の挙動を模擬するとともに、前記液膜の膜厚及び／又は液膜流量を測定するようにしたことを特徴とする熱流動現象の模擬試験装置を提供する。
【００１５】
例えば、模擬試験装置は、流路方向に所定間隔を隔てた複数の計測点において液膜の膜厚を電気抵抗法又は光学的測定法によって計測する膜厚計測手段や、液膜流量計測手段を備える。
【発明の効果】
【００１６】
本発明の模擬方法及び模擬試験装置によれば、高温・高圧の沸騰水型原子炉内に発生する沸騰二相流の熱流動現象を常温・低圧の環境で再現し、沸騰二相流が燃料棒の伝熱表面に形成する液膜の挙動を常温・低圧の環境で模擬するとともに、液膜の挙動を測定することができる。このような模擬方法及び模擬試験装置においては、例えば、スペーサを模擬する障害物を流路に配置することにより、スペーサの影響等を定量的に評価することができ、従って、本発明によれば、燃料棒の除熱性能に大きく影響を与える因子を定量的に評価することが可能となる。これは、サブチャンネル解析コード(NASCA)による解析法の精度向上に大きく寄与するであろう。
【発明を実施するための最良の形態】
【００１７】
図１は、燃料集合体の構成を部分的に示す正面図及び横断面図である。
【００１８】
沸騰水型原子炉の炉心を構成する複数の燃料集合体１が、図１に部分的に示されている。燃料集合体１は、所定間隔を隔てて垂直な燃料棒２を稠密に配置した構成を有する。燃料棒間ギャップを一定に保持するための複数のスペーサ３が、燃料棒２の軸芯方向に所定間隔を隔てて配置される。冷却材としての水が単相流の状態で燃料集合体１の下部から流入し、矢印で示すように上昇する。水は、燃料棒２の発熱によって沸騰し、蒸気及び水の沸騰二相流として燃料棒２の間の間隙４を上昇する。図１（Ｂ）には、円形断面のスペーサ３が例示され、図１（Ｃ）には、六角形断面のスペーサ３が例示されている。
【００１９】
図２は、燃料棒２の管壁５を部分的に拡大して示す部分拡大断面図である。
【００２０】
上昇する沸騰二相流は、環状流として燃料棒２の外周面を覆い、薄い液膜１０が燃料棒２の表面（伝熱表面）に形成される。液膜１０は、燃料棒２の伝熱表面の熱伝達率を向上させる。大小様々な液塊（図示せず）が燃料棒２の表層近傍を頻繁に通過し、液膜１０は、非常に激しく変動し、じょう乱波１１として示す如く、複雑且つ非定常な液膜挙動を生じさせる。沸騰二相流は、スペーサ３の近傍では、更に複雑な流動様相を示す。従って、燃料棒２の伝熱表面における液膜挙動の態様及び性質は、現状では、判明していない。
【００２１】
更に出力が上昇すると、燃料棒表面の液膜が蒸発して燃料棒表面が乾く液膜ドライアウト現象（符号１２で示す）が発生する。このようなドライアウトが一旦発生すると、熱伝達率が急激に悪化するので、燃料棒２は、その壁温が急上昇し、場合によっては管壁５の焼損（バーンアウト）現象が発生し得る。このため、このような限界条件における燃料棒上の液膜挙動を正確に把握することは、原子炉の熱設計を行う上で極めて重要である。
【００２２】
図３は、スペーサ３の近傍に生じる液膜挙動を例示する部分拡大断面図である。
【００２３】
前述の如く、スペーサ３の形状・寸法及び位置が液膜１０の挙動に与える影響は、極めて複雑である。例えば、液膜１０の膜厚は、図３に示すようにスペーサ３の近傍で局部的に大きく低減し又は増大する傾向があると考えられ、スペーサ３の近傍では、液膜ドライアウト現象（符号１２）が発生することが懸念される。燃料集合体１の流路を模式的に各燃料棒２毎の小流路（サブチャンネル）に分割して燃料集合体１内の熱流動現象を解析するサブチャンネル解析が、燃料集合体１の限界出力を求めるために近年使用されているが、スペーサ３の存在は、液膜ドライアウト現象と密接に関連することから、スペーサ３の影響をサブチャンネル解析に反映することが望まれる。
【００２４】
図４には、燃料棒２を模擬した発熱管上に生じる液膜の膜厚変動特性が示されている。
【００２５】
図４に示す変動特性は、本発明者が大気圧の実験環境において従来の実験機器を用いて実施した模擬実験の実験結果である。この模擬実験では、水蒸気の見掛け流速（jG) を16.6m/s に設定し、水の見掛け流速（jL) を0.41m/s に設定している。
【００２６】
図４に示す変動特性から理解し得るように、発熱管上に形成される液膜は、非常に激しく変動しており、大小様々な液塊が発熱管近傍を頻繁に通過する複雑な液膜挙動を示す。このような実験による観察や計測等は、或る程度までは高温・高圧の条件で実施し得るかもしれない。
【００２７】
しかし、沸騰水型原子炉（実機）の一般的な設計条件では、燃料集合体１の内部環境は、圧力7MPa、温度285 ℃に達する。従って、沸騰水型原子炉における燃料集合体１の除熱性能を把握するには、このような高温・高圧の条件で燃料棒２の表面を流れる沸騰二相流を分析しなければならず、しかも、前述の如く、スペーサ４の近傍では、更に複雑な流れが形成されると考えられるところ、仮に、このような高温・高圧の実機条件を再現する電気加熱式試験機等を使用し得たとしても、試験機内に生じる沸騰二相流の流動様相を観察し又は測定することは、現在の計測技術では、極めて困難であり、高温・高圧の実機内で沸騰二相流が実際にどのような挙動を示すのかについて観察し又は測定することはできない。
【００２８】
他方、このような沸騰二相流の流動現象又は流動様相を支配する支配因子は、気相と液相のそれぞれの流体物性、殊に、密度、表面張力、粘度であると考えられる。従って、これらの物性の絶対値が、実機の沸騰二相流における気相と液相のそれぞれの流体物性（密度、表面張力、粘度）の絶対値と実質同一又は同等であれば、観察可能且つ測定可能な常温・低圧条件の環境であっても、実機の燃料集合体１に発生する沸騰二相流の流動様相を模擬し、高温・高圧の実機内に発生する熱流動現象を擬似的に再現し得ると考えられる。
【００２９】
図５（Ａ）には、大気圧及び実機圧力（70atm)における水の諸物性が夫々示されている。
【００３０】
図５（Ａ）に物性値として示すように、水の大気圧下の飽和温度（100 ℃）における諸物性は、高温・高圧（圧力7MPa、温度285 ℃）下における水の諸物性とは、大きく相違する。従って、観察可能且つ測定可能な常温・常圧条件の環境で水を用いて沸騰二相流を形成したとしても、実機内の熱流動現象を擬似的に再現することはできず、仮にこれを行ったとしても、燃料集合体１における除熱性能の予測精度は低く、現状では実用の域に達しない。
【００３１】
図５（Ｂ）には、実機条件（70atm)における水の諸物性と、エタノール及びHFC134a ガスの諸物性が、物性値として示されている。
【００３２】
常温・低圧の環境において気体及び液体を混合して二相流体を形成した場合、両者の化学反応が進行しないと仮定すると、適切に選択して組み合わせた気体及び液体の二相流体は、気体及び液体の夫々の諸物性のうち二相流体の性質を実質的に支配する物性を適当に利用することにより、高温・高圧の沸騰二相流体の挙動を疑似することができる。本実施形態においては、沸騰水型原子炉の高温・高圧（圧力7MPa、温度285 ℃）の環境で発生する沸騰二相流を常温・低圧の環境で再現するための液体及び気体として、エタノール及びHFC134a ガスが採用される。
【００３３】
液体模擬物質としてエタノールを用い、蒸気模擬物質としてHFC134a ガスを用いた場合、これら二つの物質を混合した二相流体の密度、表面張力及び粘度は、高温・高圧（圧力7MPa、温度285 ℃）の環境に発生する沸騰二相流の密度、表面張力及び粘度と類似する。従って、エタノール及びHFC134a ガスを混合した二相流体によって上述の沸騰二相流を疑似することにより、高温高圧条件（圧力7MPa、温度285 ℃）の気液二相流動現象を観察可能且つ計測可能な常温・低圧の条件で再現することができる。
【００３４】
即ち、高温・高圧の環境に存在する沸騰二相流の流れの特性を支配する物性を特定し、この物性の物性値と同等又は同程度の物性値を常温・低圧の環境で保有する気体及び液体を適切に組合せた二相流体を形成し、この二相流体を擬似的な沸騰二相流として使用し、これにより、高温高圧条件の沸騰二相流の液膜挙動を観察可能且つ計測可能な常温・低圧の環境で再現することができる。
【００３５】
このような方法によれば、高温・高圧の沸騰二相流を用いた装置の設計・開発に要するコストを大幅に削減するとともに、開発期間をかなり短縮することができる。殊に、原子炉は、発電量が非常に大きいことから、限界出力を僅か１％向上させることができれば、莫大な利益につながるので、本発明の模擬方法を原子力の設計・開発に適用し、限界出力を高精度に設定して燃料集合体を最適化し、１％程度の性能向上を図ることができた場合、電力会社が所有する100 万kW級の原子炉では、1 万kW程度、更に発電することができ、これを現在の物価で経済効果として試算すると、１ヶ月で約1.5 億円（22円／1kWhとして試算）、１年で18億円の経済効果が得られる。
【００３６】
常温・低圧、例えば、大気温（又は室温）且つ大気圧、或いは、通常の実験機器で容易に形成可能な温度及び圧力の環境の下で上記二相流体（エタノール及びHFC134a ガス）を供給する流体供給手段と、この二相流体が供給される直線流路とを備えた熱流動現象の模擬試験装置を本発明に従って設計することができる。
【実施例１】
【００３７】
図６は、本発明の実施例に係る模擬試験装置の構成を概略的に示す断面図である。
【００３８】
試験装置は、液体及び気体を混合する気液混合器２０と、液体及び気体を分離する気液分離器３０と、液体及び気体の混合気が流通する管体４０とを備える。管体４０は垂直に配置される。気液混合器２０は管体４０の下端部に配置され、気液分離器３０は管体４０の上端部に配置される。
【００３９】
HFC134a ガスを供給する気体供給管２１（破線で示す）が、気液混合器２０の気体流入ポート２２に接続される。気体供給管２１は、気体給送手段（図示せず）を介して管HFC134a ガスのガス発生器（図示せず）に接続される。エタノール液を供給する液体供給管２３（破線で示す）が、気液混合器２０の液体流入ポート２４に接続される。液体供給管２３は、液体給送手段（図示せず）を介してエタノール液の貯留槽（図示せず）に接続される。
【００４０】
HFC134a ガスをガス発生器に還流させる気体回収管３１（破線で示す）が気液分離器３０の気体流出ポート３２に接続される。気体回収管３１は、気体還流手段（図示せず）を介してガス発生器に接続される。エタノール液を液槽に還流する液体回収管３３（破線で示す）が、気液分離器の液体流出ポート３４に接続される。液体回収管３３は、液体還流手段及び液膜流量検出手段（図示せず）を介してエタノール液の貯留槽に接続される。液膜流量検出手段は、液体回収管３３によって回収される液量（エタノール液量）の時間変化を検出する液量計を有する。
【００４１】
管体４０は、HFC134a及びエタノールの気液混合気Ｆが流通可能な円形断面の垂直流路４１を有する。流路４１の下端部は、気液混合器２０の流出ポート２５に接続され、流路４１の上端部は、気液分離器３０の流入ポート３５に接続される。
【００４２】
試験装置は、流路壁４２上の液膜の膜厚を計測するための膜厚計測手段を備える。空気・水系のシステムにおいて非接触で気液比及び液膜厚さを測定する手法として定電圧法及び定電流法が知られているが、本実施例においては、センサ間の干渉がなく、しかも、膜厚の減少に対して出力を増大する（従って、計測精度が高い）定電流法が採用される。膜厚計測手段は、流路方向に間隔を隔てた流路壁４２の２点間電位差を夫々測定する上流側及び下流側の電圧計（Ｖ）と、各電圧計（Ｖ）の計測値を記録するデータレコーダ（ＤＲ）とを備える。直流電源に接続された正負の電極対４５が、管体４０の上部及び下部に夫々挿入され、流路壁４２の液膜に通電する。各電圧計（Ｖ）の検知電極対４６が、流路壁４２の所定位置において流路壁４２上の液膜に接触するように配置される。管体４０は、非導電性樹脂等によって全体的に形成され、定電流電源の電力を印加する電極対４５と、液膜厚さ検知用の電極対４６のみが金属によって形成される。
【００４３】
試験装置は更に、差圧センサ（ＤＰ）を有し、差圧センサ（ＤＰ）の検知部４３が流路４１に設置される。検知部４３は、流路４１の流入部及び流出部に夫々配置され、金属管４９によって差圧センサ本体に接続される。金属管４９には、圧力信号を差圧センサ本体に伝達する液体が充填され、差圧センサ（ＤＰ）は、流路４１の流入部及び流出部における二相流体の圧力差を検出する。このような差圧検出により、原子炉ポンプの容量や必要性能の設定に必要な流体圧損データが得られる。
【００４４】
図７は、図６に示す試験装置を用いて実施した液膜挙動再現実験の実験結果を示す線図である。
【００４５】
液膜挙動再現実験において、エタノール及びHFC134a ガスの二相流体が気液混合器２０から流路４１内に連続的に給送された。二相流体は、流路４１を上昇して気液分離器３０に流出する。上下の電圧計（Ｖ）によって、計測点α、βにおける電位差が夫々計測され、データレコーダ（ＤＲ）に記録された。計測点α、βにおいて夫々計測された電位差の値は、下式により、流路壁４２に形成された液膜の膜厚値に変換された。
【数１】
2007316059.tif 000003
【００４６】
なお、上式における各符号は、以下の数値を夫々示す。
d    :管内径
d0  :検定棒外径
V0  :検定棒挿入時の出力電圧
I0  ：検定棒挿入時の定電流値
VTP  : 膜厚測定時の出力電圧
I  ：膜厚測定時の定電流値
【００４７】
図７（Ａ）及び図７（Ｂ）は、エタノール及びHFC134aの流量が異なる二つの実験の各実験結果を示す線図である。図７（Ａ）及び図７（Ｂ）には、計測点α、βにおける膜厚の時間変化が示されている。図７において、符号JGは、流路４１を気相が充満して流れたと仮定して求められた見かけの速度であり、符号JL は、同様に流路４１を液相が充満して流れたと仮定して求められた見かけの速度を示す。また、符号τiは，気液界面のせん断応力であり、次式で定義される。符号Δtは、上下流の膜厚データの相互相関係数の最大値として得られる２つの波形の時刻のずれ時間であり、符号vDWは、センサ間の距離（57mm）をΔtで除すことにより算定されるじょう乱波速度を示す。 
【数２】
2007316059.tif 000004
【００４８】
計測点α、βは、鉛直方向に所定距離（本例では、距離約57mm）を隔てており、流路壁４２の内壁面に形成された液膜は、波状形態をなして下流側に遷移することから、上流側の計測点αに生じた膜厚変化と類似な膜厚変化が、下流側の計測点βにおいて所定時間Δｔ後に現れる．
【００４９】
図７（Ａ）及び図７（Ｂ）には、計測点α、βにおける膜厚変化が夫々示されている。各図において、計測点βの膜厚変化と計測点αの膜厚変化とを容易に対比し得るように、計測点βの膜厚変化は、その位相を時間Δｔだけ変移させた状態で示されている。
【００５０】
図７（Ａ）及び図７（Ｂ）に示すように、実質的に同一の膜厚変化が時間差Δｔで計測点α、βに顕れており、沸騰二相流の流動現象を再現するエタノール及びHFC134a ガスの二相流体の場合においても、測定上の矛盾がないことが確認された。
【００５１】
図８は、サブチャンネル解析コード(NASCA)による解析結果と、液膜挙動再現実験の実験結果とを示す線図である。図８(A)には、気相の見かけ流速と、液膜流量との関係が示され、図８(B)には、気相の見かけ流速と、液膜厚さとの関係が示されている。
【００５２】
図８に示す如く、エタノール及びHFC134a ガスの流量に応じて、流路４１内の気体流速（JG）は増大し、流路壁４２の内壁面に形成される液膜の膜厚(tFm)は変化する。再現実験の計測値は、液膜流量（図８(A)）及び液膜厚さ（図８(B)）の双方において、サブチャンネル解析コード(NASCA)による解析結果と極めて類似する。ここに、サブチャンネル解析コード(NASCA)は、スペーサが設けられない場合には、液膜挙動を精度良く再現できる。上記模擬試験装置により測定した液膜流量及び液膜厚さの双方が、そのようなサブチャンネル解析コード(NASCA)の解析結果と極めて良く一致することは、本手法の妥当性が確認されたことを意味する。即ち、本実施例によれば、主要な物性値が高温高圧下の水及び蒸気と同程度である液体（エタノール）及び気体(HFC134a)を作動流体(実験圧力: 0.7MPa，実験温度：27℃)として使用し、高温高圧の蒸気−水系の二相流動を直線流路４１において生じさせることにより、液膜厚さ、液膜流量及び圧力損失を測定可能な高温高圧下の二相流動を模擬することができる。
【００５３】
沸騰水型原子炉（ＢＷＲ）の開発・設計においては、開発から実証に至る過程で試験・実験的アプローチが重視されてきたが、近年の沸騰水型原子炉や、次世代の沸騰水型原子炉では、燃料の大型化や、格子稠密化等の理由により、試験・実験的アプローチの採用が実験機器容量や開発コスト等との関係で困難になりつつあり、試験・実験に依存した開発・設計から解析中心の開発・設計に移行することが望まれている。しかしながら、スペーサ廻りの乱流解析や、燃料バンドルにおける気液クロスフロー効果等については、サブチャンネル解析コード(NASCA)による解析では、得られない。しかしながら、図６に示す試験装置によれば、スペーサを模擬する障害物を流路４１内に配置して前述の液膜挙動再現実験を実施し、これにより、スペーサが液膜挙動に与える影響を分析し、サブチャンネル解析コード(NASCA)に適用可能なスペーサ増倍係数等を設定することができる。
【００５４】
以上、本発明の好適な実施形態について詳細に説明したが、本発明は上記実施形態に限定されるものではなく、特許請求の範囲に記載された本発明の範囲内で種々の変形又は変更が可能である。
【００５５】
例えば、試験装置の流路断面は、必ずしも全長に亘る均等又は均一な円形断面に限定されるものではなく、流路断面を楕円形断面、多角形断面等に設計し、或いは、流路断面を適当に変化させても良い。
【００５６】
また、試験装置の管体を流通する二相流体の視覚的な観察を可能にするように管体を透明又は半透明の素材で形成することも可能である。
【００５７】
更に、上記実施例では、定電流法を用いて液膜の膜厚を測定しているが、膜厚測定に定電圧法を採用し、或いは、レーザー光によって液膜の膜厚を測定する非接触型の光学式膜厚測定装置等を使用した光学的測定法を採用しても良い。  
【産業上の利用可能性】
【００５８】
本発明の模擬方法及び模擬試験装置は、沸騰水型原子炉（ＢＷＲ型原子炉）の開発から実証に至る過程において、燃料棒の限界出力を設定するための試験又は実験に使用される。本発明の模擬方法及び模擬試験装置によって得られた液膜挙動の情報に基づいて、燃料棒表面の液膜挙動に対するスペーサの影響等の如く、燃料棒の除熱性能に大きく影響を与える因子を定量的に評価することが可能となる。これは、沸騰水型原子炉内の熱流動現象を解析する上で必要な未知データの取得を可能にし、サブチャンネル解析コード(NASCA)の精度向上に大きく寄与すると考えられ、その実用的価値は、顕著である。
【図面の簡単な説明】
【００５９】
【図１】燃料集合体の構成を部分的に示す正面図及び横断面図である。
【図２】燃料棒の管壁を部分的に拡大して示す部分拡大断面図である。
【図３】スペーサの近傍に生じる液膜挙動を例示する部分拡大断面図である。
【図４】燃料棒を模擬した発熱管上に生じる液膜の膜厚変動特性を示す線図である。
【図５】図４（Ａ）は、大気圧及び実機圧力（70atm)における水の諸物性を示す図表であり、図４（Ｂ）は、実機条件（70atm)における水の諸物性と、エタノール及びHFC134a ガスの諸物性を示す図表である。
【図６】本発明の実施例に係る模擬試験装置の構成を概略的に示す断面図である。
【図７】図６に示す試験装置を用いて実施した液膜挙動再現実験の実験結果を示す線図である。
【図８】気相の見掛け速度と液膜流量及び平均液膜厚さとの関係を示す線図である。
【符号の説明】
【００６０】
１  燃料集合体
２  燃料棒
３  スペーサ
４  間隙
５  管壁
１０  液膜
１１  じょう乱波
１２  液膜ドライアウト現象
２０  気液混合器
３０  気液分離器
４０  管体
４１  流路
４２  流路壁
【図１】
2007316059.tif 000005
【図２】
2007316059.tif 000006
【図３】
2007316059.tif 000007
【図４】
2007316059.tif 000008
【図５】
2007316059.tif 000009
【図６】
2007316059.tif 000010
【図７】
2007316059.tif 000011
【図８】
2007316059.tif 000012