日 P Outer 100p  N  ‑300*

100  0.46 

1.

0  2 > < 1 0 ‑6  1x 10‑5  1

1.

9  70  0.23 

1.l 

4 > < 1 0 ‑6  1 > < 1 0 ‑5  : 1 . 4.2 

0.086 

1.

4  2 > < 1 0 ‑6  1 > < 1 0 ‑ ; : '   16.3 

* Nu

曲

er ロ f iteration of probab11ity

皿

alys1s.

( 2 1  

収 束 判 定 因 子 と 圧 力 損 失 の 有 効 数 字

解を収束さぜるためには初速因子が重要であるとと怯(1)で指摘した通りであるが，まえ.とれ」ζ対して収束

判定因子を変えてやるのも一法である。計算時間との関連で結当然収束判定因子姐厳し〈をい方がよいが.し

か し 計 算 精 度 上 そ う む や み

κ

甘くするととはでき在い。本コート・で'/i.一 般

K.

内側ループの判定因子同舛剣

Jレープ￠それに比しやや厳しくして台〈方首都の'反束あるい同精度上よいようである。とにで位計算澗匿に捜漣し

‑3 9  ‑

J A E R 1 - M 7 0 5 9

* h ft £JR*ffiEf§ J ciEtm

5x10"

1x10"

5x10

^{i - 5}

»> l x l 0 ~

⁵

5xl0~

⁶

1x10 ^{- 6 _L JL J_}

3 4 5 S i g n i f i c a n t d i g i t

Fig. 4.12 Relation between convergence f a c t o r and s i g n i f i c a n t d i g i t in c a l c u l a t i o n of pressure l o s s

* i S I S H i K &v>X&%K<DBctlt&9i*

^:

fttfr§ t T-Ht L X^ &*» t ^ 5 *»»^©*'MfflRo

40

JAERI‑M  7059 

て.本モデル

K

合いて計算された運転負荷主主

40

婦の場合の解析例

K

みられる収束判定限「と!王対負

.

，  

失.1

P .

の有効数字との関係を，参考として下図

K

示す。収束状況が極めて急激在ため目ほとんどの 場合，図中の最小ラインを越えて収束している。つまり，

全流路

K

ついての圧力績失は.ほとんどのものが有効数字

5

桁目ぐらいまで一致している。

1  x  10‑

⁵の収束判定悶子の場合でも.

5x10‑

匂

•

lx10‑

匂

﹁トtlトlト1ト E 3 F b   e

n‑

u A U  

' i '

ム

X V

︽

唱4 E J

5x10

四

5

陣 叫 倫 相 O山田U 咽 刷 副 Uロ叫岨

μ ωh F

ロ

o u

7  6 

5 

」 3  4  1  2 

1x10‑6 

Signific

田

t digit 

Re

1ation between convergence  factor and  significant digit  in calculation of pressure  10ss 

Fig.  4.12 

ある判定因子1'1::1>‑いて各自

E

飴の底力鍋失が向桁自主で一致してのるかという有効窓字の畳小桁散。

‑40‑

JAERI‑M  7059 

て.本モデル K 合いて計算された運転負荷率 40~ の場合の解析例』てみられる収束判定応「と!王対負

0} 

失.1

P .

の有効数字との関係を，参考として下図

K

示す。収束状況が極めて急激在ため目ほとんどの 場合，図中の最小ライYを越えて収束している。つまり，

全流路

K

ついての圧力損失は.ほとんどのものが有効数字

5

桁目ぐらいまで一致している。

1  x  1 0 ‑ .

の収束判定因子の場合でも.

5x

lO‑匂

•

1x10‑

匂

﹁トtlトlト1ト r 3 p b   e

A‑

U A U  

噌 ム 噌

4

x x  

旬 ム

EJ

5x10

^四

5

障ω

‑ H

叫UO 咽 刷 副 Uロ国岨何世hFロ

ou

7  6 

5 

」 3  4  1  2 

1x10‑6 

Signific

血

t digit 

Re

1ation between convergence  factor and  significant 

d1

git  in calculation of pres8ure 

108S 

Fig.  4.12 

ある判定因子 K主砕いて各流飴の圧力鍋失が向桁自主で一致してのるかという有効窓字の量小桁散。

‑40‑

J A E R I - M 7 0 5 9

h t ft £

^ ^ - c m ^ t i - X ^ i , DREAM- FLOW => - Ktt, &Uti xjp ( f;- . -f > • -• c ,- ? fr4*

$ £ £ tn^ff-5«SfifiiffiOp^Jffflii frf$#f©£B&£ L t ^ t , Si«to!c©#£iMS!(CSM I T i l *

4 © r i> h „

i L t , IJSfffi#»Otf^oSKff5C.rt^«i(*gtfa:^i'A,4'RFW-:flfliSr*-fi>i^5 t i t f t J .

»f-lJirtffi'A;K L t n *> t o t * 6 ,

l l L i r f f * 5 £tHffiffl.b©'*i]ft#(Cl+»B$R5i©M*^tilfflWT4^

⁰

S * t t . ffi.H.*!>m

mm&umK*)i-z,mi&ftc «»«•*) a^t^ss*< ioo@y±300fagKotiuii- ®»*

^ T © $ # f # l 7 > i b J t 5 I L T ^ o o $)Mi, H & « K t t t t § t l t t ® a © ^ £ g * «'h@tSK ret * t f 3 0

~5oisi)Ti+»srrr«)-ot:4,1^^5 5o - t u , isiiffffi*^* t>^,^ow-^i&^mmsmom

*©t£g©ffffil«!?tt, i W K 6 & ^ t t - 2 ^ | c r e ; i ; - t E K l i i i £ f f i 4 i'ogtft&^ffifcffl^-rfFlffi-*

M J : i A t * 5 5o

DREAM-FLOW = - Kti, SJ»«©lSiSS*^IO

^;

flffl®*(-©lt#^ra

^,

)f«:>«-t4S0ft^* t . 7 t i ^ t f S 4 * ( C * L ^ * § W i i f i * ' ^ l l l K i P © #

5

' L - ^ * b n 4 X 5 4 » £ < ©8M&*fe44iF<fr

£(**^TiK§StL7cW#f;!5;fj4;i4^<> t o j k f c . * V 7 ^ f « f t T f i 4 , IttfLti&tflK l=] LZ I 5 4JUS4 l o K f t ^ 5 - a r - t > i * 4 > 4 ^ t t ^ * K 4 J E B t t t f f i i - 4 t * L T ^ * © « E I » » 4 r

^ • S - ^ & S T } ; * - ^ L r t ^ o t , Wf,nfeiKttfi<jeH©ie*fl'£.5K:i£;i;T:fffFffi-tSCti 4 6 5 * , e i o ^ f c o ^ - r w : . 3-SL *=> - Kt#gfi?j®ia*^iHfe]!P»cfc^siP-Q©!!*ffifi)Stt K f i J f f l L - r i ^ < ^ 5 l t o - t « S 7 t ^ t ^ ^ T ^ i o 4 7t. 5 b K i ! ) # < © ) S * f * i i t o , IHHIRK fiPL-7t»Wflctf§;l, * W * f c t f a f e « 4 * S L - t ^ < o t t T * j .

S B . * = - Ktt±(cieL7tJ:5K#iW«E#^5l^r©«

s

'U©*!Sfti!)S:l+K:if0fflLTi!i><^

£?*•&#. ^&, c o i o g ^ t r a L t , ±mm#M&kinKi>ik<Dmn&.inmLtt

i ) * 8 « t c i - f t S » f c i t * © - W T S

ii) * 8 ? f t f i . » e a $ J i f t o i f & f ) -

^,

7 > - ^©^515 4 1

-JAEHI‑M  7059 

あ と が き

本報告書で記述されている

DREAM‑FLOW

コードは.高温ガス炉〔ピJ ・イJ ・フロックをが 料とする炉心)の冷却流路を多並列

l i f L

路とした炉心の熱流動解析主として日

I i d

のおif:F干伽守行うた めのものである。

との計算コ ドは，流れの遷移域

v c

;t'ける摩擦係数のデ タのばらつきに

1 1

目した各iiiL時，")流

1 i 1

i誤差ととれ

v c i

'l'う燃料品よ度の誤，‑r

J f

佃]とを解析の主眼とLている。遷移城の'‑"樗係数

v c

i"ilてはあ る統計的1"4'買を仮定し，その解相l子仏としてモンテカノレロ注を適用する。そ Lて 専 ム71プヒ統計的 1

"

1 

質(流婦や燃料品itt街似度

v c

X".jして平均値なり練準偏差在り

1

によってとれらの誠差評価を行在う ものである。

、う然これらの統計Ii

H

^{てよ‑，てuT}t

l i L

在り.lfut皮;まりの悩j.

I i : i

汗・価制‑r.~-う場合. ~J の進本的在!自i閣 と

L

て.

. 1

1(.燃係般のばらつき

K

什うとれらの統計随がどん在特例分布を示すかというととがある。

これは初 4 帝ーの，m司例ヰ今までハ仙の計調例 K よって派鼠在どの特例分布は IE 規却~l持t'hiU畳分布) r 近似され仰る ζ とが灯[11> られているのそれ tl. 日目量や燃料策，~IAà/B:の統計分布が1現!愉的 K つられ た1I規Ith産自を1分 よ <t.:. l~. しているとと.さらに.紛:il1l量 x' を検定Jl\，準とする x' 一機~~てよ円てそ の判断をt

官吏 K

しているものである。

書 l

，2 の問題1':1.計算繰~しI，Jj数L打切り

1 " 1

歓)と統計最 t特に平均値と僚

i t

偏去二

1

との関遂であ る。当然繰返し￨吋数が僧大すれば得られた統計的性質の信頼度は大きく在るが.余り多〈の富十型車の 繰返しを行在うこと同使用上の情IJ約特に計算時間との関連から実用的で在い。我々は.流量左台燃 料最高温度

K

対する統計分布(度数分布)が比較的落着<

!  00

回 以 上

300

回 程 度 の 繰 返 し 回 数 を 今までの解析例から推奨している。勿論，一般的Kは統計的処理の可能季最小回訟(たとえば

30

~50 自)で計算を打切ってもよいだろう。そして，誤差評価を行在う場合の平均値や標準偏差の偲 々の怯質の信頼限界は，目的

κ

^{あるい同必要}K応じて区間推定法在どの統計的手法を用いて評価す ればよいであろう。

DREAM‑FLOW

コードは，計算機の記憶容量品、よび実用面か与の計算時間に対する制約等から，

たとえば第

4

章

K

示した多目的高温ガス実験炉の炉心

K

みられるよう全数多〈の流路から在る炉心 全体をすべて網緩した解析が行なえ在い。とのため.オリアイス領域内で等価在，たとえば熱的K

1 0 J

じよう左流路を

l

つ

K

代表させるとかあるいは代表的念流路を抽出するとかして全体の流賂叡を 減少さぜる必要がある。したがって，得られた統計的性質の結果は必要

K

応じて再評価するととと なろうが，この点については.今後.本コードを多目的高温ガス実験炉(‑(;1;‑ける炉心の熱流動設計 K利用していく際確めて置きたいと考えている。また.さらKより多〈の解析を進め， よP実 際K 即 し た 解 析

K

備え.必要あらば改善を施してい〈つもりである。

当面.本コード11上K記したようK多目的高温ガス実験炉の炉心の熱流動設計K利用してゆ〈考 えであるが.今後.との穫の誤差評価K対して.上記摩擦係数以外にも次のようま外乱を考慮した コードの開発念いしは改良を行い.より一般的左誤差評価コード

K

改善してゆきたいと考えている。

i 

)遷移域(‑(;1;‑ける熱伝達率の刊椛さ

i i

 

)摩擦係数.熱伝達率以外の工学的パラメ タの不獲さ

‑4 1   ‑

JAEHI‑M  7059 

あ と 刃 : き

本報告書で記述されている

DREAM‑FLOW

コードは.高温ガス炉〔ピJ ・イJ ・フロックをが 料とする炉心)の冷却流路を多並列

I n L

路とした炉心の熱流動前析主として日

I i d

の設ir':F干伽守行うた めのものである。

との計算コ ドは，流れの遷移域

v c

;t'ける摩擦係数のデ タのばらつきに

1 1

目した各ii¥L路，")流

1 i 1

i誤差ととれ

v c i f

う燃料温度の誤，‑r::

J f

価とを解析の主眼としている。遷移城の""樗係数

v n i l

てはあ る統計的1"4'質を仮定し，その解相

i

1二仏としてモンテカノレロ拡を適用する。そして，得ム

7 1

た統計的 1

"

1 

質(流

1

武や燃料品il:，s'i制度

K

対して平均値在り僚準偏差走り

1

によってζれらのiis差評価を行在う ものである。

吋然これらの統計Ii

H

てよ‑，てu1

t l i l

在り.

! ! ; t

皮在りの総

1 ] i

汗・価射‑r.~-う場合. æJ の~本的吐問題 と

L

て.

. 1

1(.線係数のばらつき

K

什うとれらの統計値がどん在特何分布を示すかというととがある。

これは初

4

阜の，

m

司例ヰ今吐でハ仙の計調例

K

よって陀匝在どの特科分;f[iは正規却J特t'

hiU

畳分布) r 近似され仰る ζ とが灯I められているのそれ tl. 流量や燃料照，~

i ! . l J B ' .

の統計分和が1司!愉的

Kfi

られ た11.1乳Ith似を1分 よ <t.:. l~. しているとと.きら VC. 統計 i量 x' を検定)!\，準とする x '--1喰~~てよ司てそ の判断を

C 1 C k . v c

しているものであるの

理 l

，

2 の問題1';1.計算繰~し

I ! J j

数

L

打切り

I I J ] { ' x

)と統計鼠，

持K [

平均値と標章偏長

1

との関遂であ る。当然繰返し￨口]数が噌大すれば得られた統計的性質の信頼度は大き〈在るが.余り多〈の計型車の 繰返しを行在うこと同使用上の借IJ約特

K

計算時間との関連から実用的で在い。我々は.流量左

h

燃 料最高温度

K

対する統計分布(度数分布)が比較的落着<

!  00

回以上

300

回程度の繰返し回数を 今までの解析例から推奨している。勿論，一般的

K

は統計的処理の可能季最小田叡(たとえば

30

~5D 回)で計算を打切ってもよいだろう。そして，誤差評価を行左う場合の平均値や標準偏差の個 々の但質の信頼限界は，目的Kあるいは必要K応じて区間推定法在どの統計的手法を用いて評価す ればよいであろう。

DREAM‑FLOW

コードは，計算機の記憶容量品、よび実用面かトの計算時間に対する制約等から，

たとえば第

4

章

K

示した多目的高温ガス実験炉の炉心

K

みられるよう在数多〈の流路から在る炉心 全体をすべて網緩した解析が行吐え吐い。とのため.オリアイス領域内で等価在，たとえば熱的K

I c J

じよう在流路を

l

つに代表させるとかあるいは代表的左流路を抽出するとかして全体の流路数を 減少させる必要がある。したがって.得られた統計的性質の結果は必要

K

応じて再評価するζとと 左ろうが，この点については.今後.本コードを多目的高温ガス実験炉

K

;1;'ける炉心の熱流動設計

K

利用してい〈際確めて置きたいと考えている。また.きら

K

よ

h

多〈の解析を進め，

ょ P

実 際

K

即した解析

K

備え.必要あらば改善を施してい〈つも

P

である。

当面.本コードは上

K

記したよう

K

多目的高温ガス実験炉の炉心の熱流動設計

K

利用してゆ〈考 えであるが.今後.との穫の誤差評価

K

対して.上記摩擦係数以外

K

も次のよう在外乱を考慮した コードの開発在いしは改良を行い.より一般的在誤差評価コード

K

改善してゆきたいと考えている。

i 

)遷移域VC;j;.ける熱伝達率の刊配さ

i i

 

)摩擦係数.熱伝達率以外の工学的パラメ タの不確さ

‑4 1   ‑

J A E K 1 - M 7 0 5 9

-J y a yimt. t±, --')<•> A o W M K ^ ' T W I s W ^ ' s l o ^ & i f i J t i t ^ t f unit i 0 * ( f

42

-JAEIiJ ‑M  7059 

本コ ドを作成するjt(当 1)， 事目的炉設計研究翁の荒舛長利氏には開発の初期の段階で千イスカ

ッ シg ; /を頂き，また.ヘリウムの物件

i (

直

K

ついては同研究家の瀬谷東光氏が存呼ず れたものを使

用させて頂いた。 ζ こjt( ~ê して深甚在る謝意を表します。

‑42‑

JAERJ‑M  7059 

本コ ドを作成するK当1)， 宰目的炉設計研究箪の荒井長利氏には開発の初期の段陪で干イスカ ッショ Y を頂き，また.ヘリウムの物件 (i直 K ついては同研究事の瀬谷東光氏が~呼 fれたものを使 用させて頂いた。 ζ乙jI[~t.:して深甚在る謝意を表します。

‑4 2  ‑

ドキュメント内 ミ_~) -l--^;'r ~l ( }K.i. H «'. i' OW ~~ ">i!i~ ウ J - JAERI.M JAERI-M v,t:-- ~としご. hit':. ^'cffli:^jf] L :^«>.~'t~!;ヘ W; ヘ,'}て A て j- へ f 代 111~~にど j, (ページ 44-50)