KFKI-1989-13/G
MARÓTI L., TROSZTEL I.
FORRÓPONT ANALÍZISRE ALAPOZOTT ÜZEMVITEL A PAKSI ATOMERŐMŰ VVER-440 REAKTORÁBAN III. RÉSZ.
Hungarian Academy of Sciences CENTRAL RESEARCH INSTITUTE FOR PHYSICS BUDAPEST
KFKI-1989 13/G PREPRINT
FORRÓPONT ANALÍZISRE ALAPOZOTT ÜZEMVITEL A PAKSI ATOMERŐMŰ VVER-440 REAKTORÁBAN III. RÉSZ (A MUNKA AZ OKKFT G/11 1. ALPROGRAM 1.11.01 FELADATÁNAK TELJESÍTÉSE KERETÉBEN KÉSZÜLT) MARÓTI L., TROSZTEL I. Központi Fizikai Kutató Intézet 1525 Budapest 114, Pf. 49
HU ISSN 0368 5330
Maróti L., Trotztet I.: Forrópon! analízisre alapozott üzemvitel a Paksi Atomerőmű WER 440 reaktorában III rósz KFKI 1989 13/Q KIVONAT A riport extrém határértékek számítására alkalmas módszereket ismertet Fenti módszerekkel olyan paraméterek határozhatók meg. amelyek a biztonsági számitások alapvető input adatai A bemutatott módszereket a WER 440 típusra alkalmazzuk, s a kapott eredmények lehetővé teszik abszolút konzervatív számítások elvégzését
Л. Maроти, И. Тростел: Режимы работы на реакторгх ВВЭР-МО АЭС Пакш, основанные HŐ анализе горячей точки. Часть I I I . KFKI-1983-13/G АННОТАЦИЯ В отчете даны методы расчета экстремных граничных условий. С помощью этих методов можно определить параметры, являющиеся основными входными данными для расчетов безопасности. Представленные методы были использованы для реакторов ти па ВВЭР-МО; полученные результаты дают возможность проведения абсолютно консер вативных расчетов.
L. Maróti, I. Trotztet.: Operational procedure based on hot spot analysis at the WWER 440 type block of Paks Nuclear Power Plant Part III KFKI 1989 13/Q ABSTRACT The paper describes some methods that can be used for the determination of extreme marginal values of parameters that are basic input data for reactor safety calculations The methods presented in the paper are applied for the WWER 440 type and the results make It possible to carry out totally conservative calculations
- 1-
I. Forrópont analízist ti intő hidraulikai határértékek a Paksi Atomerőműben
Bevezetés Л VVER-440 típusú atomerőim Den a teljesítmény korlátját a főkeringető szivattyú pillanatszerű beékelődése esetén lejátszódó tranziens folyamán a forráskrízis elkerülésének feltétele jelen ti. Amennyiben az aktív zónát a fenti kritérium alapján akarjuk üzemeltetni, akkor rendkívül fontos az üzemzavari folyamatot szignifikánsan érintő hidraulikai jellemzők lehetséges legked vezőtlenebb értékeinek ismerete, elsősorban azon jellemzők ese tében, melyek ismerete bizonytalan. Szükséges ez a konzervati vizmus azért, hogy a reaktor biztonságos működését mindenképpen biztosítsuk. Jelen dolgozat az aktív zóna hűtőközogforgalmnt érintő legfonto sabb hidraulikai jellemzők szélső értékeit és ezek hatását ismerteti.
Áramlás-visszafordulás a kieső hurokágban Az 1 FKSZ beékelődése típusú üzemzavar esetén a védelmi rendszer működésének feltételezett késése a teljesítmény csökkenésének kezdő időpontját illetően ДТ = 3 s a Paksi Atomerőműben. Mint hogy a kieső hurokág szeparálását végző tolózárak működtetése is késik és zárási sebességük rendkívül alacsony, a hurokág a ДТ = 3 s időtartam alatt névleges teljesítménnyel üzemelő aktív zóna hűtőközegforgalma szempontjából bypass-ként működik. A kie ső hurokágban fordított irányban áramló hűtőközeg az aktív zónc. forgalmát csökkenti és hűtését rontja. Rendkívül fontos tehát annak ismerete, hogy az áramlásvisszafordulás a kritikus ДТ=3 snál rövidebb idő alatt következik-e be. Ha igen, akkor, amint azt korábbi vizsgálatok is tanúsítják, a DNBR . fellépése a
- 2 -
szivattyú beékelődését követő
3 s elteltével következik be.
Az áramlás visszafordulásához szükséges idő becsléséhez szükséges a bypass-ként működő, kiesett hurokág visszforgalmának ismerete. A Paksi Atomerőmű 1. blokkján végzett üzembehelyezési mérések alapján ennek a forgalomnak és a kieső hurokágban a visszáramlással szemben mutatkozó szivattyú ellenállásnak a meghatározá sára vizsgálatokat végeztünk és az eredményeket az [1] riportban összefoglaltuk. Az eredmények fontosságát az alábbi számítások megvilágítják. 1 FKSZ beékelődésekor a kieső hurokágban az [1] irodalom szerint a recirkulációs forgalom
Q
3
b y
= 1990 m /ó
amely forgalom akkor jön létre, ha a szivattyún kívüli hurokele mek ellenállástényezője a fordított irányú áramlás esetén mege gyezik a korábbival, a szivattyú hidraulikai ellenállása pedig С ^ 60
Az áramlásvisszafordulás időtartamának becslésénél a súrlódási tagot elhanyagoltuk, mivel az áramlás leállása alatt a súrlódás az áramlás leállást gyorsítja. A visszaáramlás jelentősen kisebb forgalmához tartozó súrlódás, ami a recirkuláció kialakulását gátolja, ugyanakkor legfeljebb kompenzálja az ellentétes irányú elhanyagolást. Ilyen feltételek mellett az áramlás visszafordulás időtartamának meghatározására a következő egyenlet írható fel:
2
* d* PRe* T "
A«
A
- m Ül * dT -' ,* - (vvissza . ) ^ _ v előre m
,
л /
- 3 -
2 p, -~ 1 4 d
ahol
m =
L. . hurok
azaz az a folyadéktömeg, amit a hurokág magába foglal.
v = előre
névleges 6
v... = vissza
Q
4 ,2 a n
by ,2 •* а и
3600
3600
AT értékét az /1/ egyenletből az alábbi kiinduló adatokkal határoztuk meg: 3
Qnévleges , , =
41760 m /ó
d = 500 mm P
= 746 kg/m
x
4 4
Чгагок* '
3
1
m
é s 5 hurok működésekor ( l á s d Ap
R e
=2,7
Ц])
bar
Fenti adatokkal a keresett visszafordulási idő: AT = 1,92 s A AT* 1,92 s igen jó egyezésben van az ÜMBJ-ben [2] közölt szovjet görbével, amely szerint ez az idő valamivel kisebb, mint két má sodperc . Vizsgálatunk célja annak megállapítása volt, hogy a visszafordu lás a beékelődést követő 3 s időtartamon belül létrejön-e. Azon ban feltétlenül megemlítendő, hogy az OKKPT A/11. összefoglalója-
- 4 -
ban [3] a kieső szivattyú hidraulikai ellenállástényezőjére С = 15,8 értéket adnak meg. Ha ezzel az értékkel számításainkat megismételjük, akkor
és
3
Q.
= 3878 m /ó
ЛТ
= 2,33 s
adódik. Ezek az adatok korábbi megállapításunkat nem érintik jelentősen, ugyanakkor fel kell hívnunk a figyelmet arra, hogy ilyen bypass forgalom mellett a reaktoron átáramló vízmennyiség öt szivattyú és egy kieső hurokág esetén
Qr. * névleges Q ,. Re,5с = 0,782 ' v
w
amely érték ellentmondásban van a korábbi vizsgálatainkban figye lembevett, az UMBJ-ben közölt
C
Re,5
=
°'
8 1 5
Q
* névleges
értékkel. Az üzembehelyezési mérések alapján meghatározott Q. adat felhasználásával kapott
Q
Re,5
"
°'
8 2 7
Q
* névleges
az ÜMBJ adataival kielégítő egyezést mutat.
3 « 1990 m /ó
- 5 -
A kazetták belépésénél elhelyezett szűkítők hidraulikai ellenállása Az aktív zónába betöltött kazetták alatt elhelyezkedő szűkítő elemek hidraulikai ellenállása lényegesen meghaladja az aktív szakasz veszteségtényezőjét. Ez azért van, mert ilyen módon kerülhető el esetleges instabilitások fellépése a hutoköze;gáramlásban, elsősorban olyan üzemzavari viszonyok között, ami kor kétfázisú áramlás kialakulására kerül sor. Ez az eset áll elő 1 FKSZ beékelődésekor is. A hűtőközegforgalom konzervatív megítéléséhez szükséges tehát a szóbanforgó, hidraulikai vesz teségek tekintetében rendkívül összetett, szűkítőelemek ellen állástényezője minimumát megbecsülni. £ -и Ч-А minimális értéke S Z UK1tu
esetén okozza ugyanis a kialakuló kétfázisú áramlás a legnagyobb forgalomcsökkenést a kritikus kazettában. A szűkítő ellenállástényezőjének legkisebb értéke úgy becsülhe tő, ha például egyetlen éles átmenettel közelítjük az egyébként igen bonyolult geometriát. Mérőperem típusú éles átmenetek vesz teségének számítására a [4] irodalom a következő összefüggést ajánlja: 6p •--• Ap (1 - 1,1 m) ahol
6p
Д
Р
az átmenet maradó vesztesége
=
P 2.1 2 o ~2mV
A m = -r^ Aо A m A о
(
.. 1 ]
a szűkítés mértéke
a legszűkebb keresztmetszet ' az alsó keverőtérből egy kazettára jutó keresztmetszet
- 6 -
Amennyiben a kapott ellenállástényezőt a zónában uralkodó se bességre akarjuk vonatkoztatni, akkor át kell térnünk a zóna sebességre. Aоvо = A,kaz v,kaz felhasználásával kapjuk, hogy 2 6
P
=
v
! kaz "
2
ilZ
)
l~
1
1
1
l " ' *» m
. , M„ = A,. ahol —Akaz г— о és így a keresett ellenállástényező
2
'szűkító^ ^* m
1
• Ы
»I
A szűkitőelem legkisebb mérete ismereteink szerint 50 mm átmérőjű. Az alsó keverótér átmérőjére 3000 mm-t felvéve Szűkítő
=
1 7
'
6
a d ó d i k
'
A [3] OKKFT összefoglalóban közölt csehszlovák adat a szűkítőre С = 25,4 Mint korábban említettük, az alkalmazott számítási módszer a mi nimális érték becslésére használható. Ennek figyelembevételével a 25,4 számérték igen jónak tűnik, mivel a valóságos adat meg kell, hogy haladja az általunk számítottat, azonban túlzottan nagy eltérés a számítást kétségessé tette volna.
- 7-
A kazetták közötti rés hűtőközeg forgalmának becslése A teljesítménykorlátozó üzemzavar forgalmának konzervatív fel vételéhez szükséges tudnunk a kazetták közötti résben, hűtés nélkül, eláramló folyadék mennyiségének maximumát. A kazetták közötti résbe a hűtőközeg a fűtőelemkötegek oldalfalán elhelyezett nyílásokon át kerül, majd a falak között zavartalanul áramlik felfelé a kilépésig. A legnagyobb forgalmat akkor kap juk, ha a résbe történő belépés, azaz a nyílások - egyébként valóságosan is igen kicsi - hidraulikai ellenállását elhanya goljuk. Ekkor a résben létrejött nyomásesésre írható, hogy H rés de , res
r
p 2
2 res
A kazettán belül ugyanakkora nyomásesés hajtja át a hűtőközeget és írható, hogy
Д
Р
=
+
< kaz df~ X
n
s
£ P
sp
+
£
ki
+
£
v
be> f kaz
kaz Az aktív zónában a kazetták közötti osztás T. = 147 mm. A kazetták külső laptávolságára 144,2 mm-t felvéve a rés egyen értékű átmérőjére kapjuk, hogy d
= 5,6 mm rés
míg a kazetta egyenértékű átmérője d e
= 8 , 5 mm. kaz
A [3] irodalomból vett adatok alapján, azaz
'be - °' £
4
ki • °'
c
e p
-°'
5
2 7
-
8
-
értékek felhasználásával, üzemi viszonyok C
összes *
7
°'
*
X
feltételezésével
H
kaz
d
'kaz
Mivel a vizsgálatok Reynolds szám tartományában v d
-0,2
A. = áll * (
a kialakuló sebességek viszonyára azt kapjuk, hogy
e
x
rés ' 1,7 ( j - ^ ) kaz
'rés 'kaz
2
1
1 1,8 = 1,02
e
Következésképpen a rés hűtőközegforgalma lényegében a kereszt metszetek arányában alakul, azaz
'rés w
A
rés
kaz
0,0802
kaz
és ezzel a zóna bypass forgalom maximumára azt kapjuk, hogy
0»
Q
'rés rée kaz
0,0743
+Q
Az a körUlmény, hogy számításaiban a GIDROPRESS
- 9-
A zónahűtés gőzfejlődés okozta romlásának határértéke Miután a hidraulikai adatok szélsőséges értékeit meghatároztuk, lehetőség van annak kiszámítására, hogy a legterheltebb kötsg milyen mértékben "dobja le" magáról a forgalmat a gőzfejlődés következtében. Ehhez ki kell számítanunk a legterheltebb köteg ben a legkedvezőtlenebb körülmények között létrejövő gőztartal mat a következő egyenletek alapján: Д
N.kaz,max = n,V
Кq К Ч
1,12 =
Kck Z
=
1 , 3
34
1 0 9
9 —
Ь 3 5 • 1,15 • 1,12 = 6,851 • 10
6
W
Itt az 1,12 szorzótényező 10 % többletteljesítményt tételez fel, amit 2 % hibával mérünk. Ugyanakkor irható, hogy N.kaz = Gm m .
(i. ki. - i. be)
ahol r
max
G . -ot akkor kapjuk, ha a korábban kiszámított maximális hurokmin ^ recirkulációt és maximális zónabypass áramlást tételezzük fel. J
Ekkor G
. = 11760 JJbe min
G
min
349 =
1 8
'
7 4
0 > 9 2 6
л
0 7 8 2
3600 k
*
/ s
p= 122,6 bar és T. = 267 + 2 °C felvételével azt kapjuk, hogy be xmax = 0,035
-10-
Végül az x - h o z tartozó kétfázisú szorzóra a [6]-ban közölt m =
Шал
formalizmus felhsználásával Ф
2
= 1 , 7 0 adódik.
°ki A forgalomcsökkenés mértékének megállapítására a következő összefüggést használjuk f e l : Д
+
v
Р - < CA 4» 5 o -
( C
Ч в> I í
+
с
A
v
2 azaz a normál - gőzt nem tartalmazó - kötegekhez képest a Ф -et létrehozó gőztártalom hatására forgalomcsökkenés lép fel. A csökkenés mértékének konzervatív becslésére feltételezzük, hogy 2 Ф a kilépési értékre a belépéstől közvetlenül lineárisan nő, azaz 12 *o
1+
_
*r.k<
_ X
'
3 5
Az ellenállástényezők értékeire az alábbi, ugyancsak konzervatív értékeket használtuk: C
C
A - C
+ C
8 Z Ű k í t 6
B "*• Ie"
+
n
Ц р sp
+
b e
C
= 17,6 + 0,4 = 18,0
ki '
4
'
4 6
+
2
'
7 +
°'
5=
7
'
6 6
A csősurlódási tényező kiszámítására a korábban is használt 1 8 4
х - ^ ЦR e
0
'
2
formulát alkalmazva a forgalomcsökkenésre azt kapjuk, hogy v.
Qj
TTо
" T o w
* °' ' 963
Tehát a csökkenés mértéke nem haladja meg a 4 %-ot,
- 11 -
II.
A biztonsági korlátok számítása a COCONT kódban
A biztonsági limitek számításához a COCONT programot el kellett látni a fűtőelem maximális hőmérsékletét meghatározó algoritmus sal és egy olyan módszerrel, ami az 1 FKSZ beszorulásos tranzi enst kvázistacionárius közelítésben biztosan konzervatív módon kezeli. A maximális üzemanyaghómórséklet számítása A fűtőelem hőmérsékleteloszlásának számítása a hővezetési egyen let segítségével az 1. ábra jelöléseivel történhet.
burkolat UO, rés
központi furat
1. ábra
- 12 -
Névleges üzemi állapotban a felületi hőfluxus a teljesítményből meghatározható q
3
= N
/ dnH
A hűtők3-?egoldali COCONT számításból ismeretes a hőátadási tényező és folyadékhőmérséklet értéke, így a burkolathőmérséklet kiszámít ható 3 3 a foly q
A továbbiakban megoldandó hóvezetési egyenlet a következő [7]: I r ahol
d dr Q
dT dr = q
v
v
v
|.
3
A fenti egyenlet megoldása л burkolatra és a résre hőforrásmentes eset feltételezésével a következő eredményre vezet: q T
=
2
T
+
3
1
r
3 r
X burkolat
J
l n
3
A
2
r~ 3
J
r
és q T
i
= T
2 ТГ~
1
+
l
Végül az r
\és
r
3
r
2
x
T
l n l
r
2
helyen fellépő maximális hőmérséklet a [7J irodalom
ból vett
függvény felhasználásával oldva meg a hővezetési egyenletet, a következő formulával számítható: Q
T
max *
, B
+
T
l
, l e x
r
v» P TS l -
rl ni n
í ° 1l «7" 4--!,],
F^
+ +
О
- 13 -
A T számításához a COCONT-ban az alábbi adatokat használjuk max fel geometria
r = 0.7 mm r. = 3.765 mm ij 3.86 mm r, = 4.55 mm =
hővezetés X
=
UO
B+T
w
m
c
a h o 1
í / '° í
A
2
burkolat * X
rés
=
2 0
-
2
*He * ° '
W / m
2 4 2
'° w / m
= 3172,9 W/m В = 287,6 °C
C
'°
c
Névleges üzemállapot esetén a legterheltebb fűtőelem maximális hőmérsékletére T„_ = 1611,0 °C adódik. л
U\J л f ТПЭХ
1 FKSZ beékelődéses üzemzavar kvázistacionárius közelítése Az 1 FKSZ pillanatszerű megállásakor lejátszódó üzemzavari folyamat alatt a legfontosabb jellemzők változását a 2. ábra mutatja. A teljesítmény-leszabályozás késése 3 s, majd a teljesítmény az ábrán látható módon átáll az 5 hurokkal történő üzem ese tén megengedett értékre. A hűtőkőzegforgalom a hűtőközeg inerciája következtében az ábrán látható módon változik. A 2. ábra alapján megállapítható, hogy a T - 3 s időpontban & folyamat igen közel van az állandósult állapothoz, hiszen a paraméterek időszerinti deriváltjai közel zérus értékűek. Ugyanakkor a DNBR . legkisebb értéke mindig ebben az időpont ban lép fel, hiszen a leromlott hűtés mellett ekkor legnagyobb a betáplált teljesítmény az akkumulált hőt is beleértve.
-
юо]
/
14 -
Л ТЦНВмЛм KCRCSZTlX X n A k Ä O ' c A l M l W u s i k
'
«о
I--!•-
со
I - - - l - f , ,i
N^BCAKII
J,
Г
-J
fc Ol
IllkS'fc- ГОМПЛЦОМ
»
_..XJ -Х_±_Х «o~
»о
"во
юо Ü >
MO
' MO
I
I
I
«0
MO
«O
B&(3ECJ
Ш _И
I
вб
1 £BC)
DO
2. ábra Az egyik fSkeringető szivattyú (GCN) beékelődése a 6 üzemelő szivattyú közül
ÍS3C)
- 15 -
Ennek következtében kézenfekvőnek látszott kvázistacionárius közelítés alkalmazása, de meg kellett vizsgálni, hogy a le hetséges közelítések közül melyik vezet konzervatív, vagy konzervatívabb eredményre. A probléma vizsgálata során arra a következtetésre jutottunk, hogy két úton közelíthetjük
a valóságos viszonyokat.
Az egyik esetben a hűtés pillanatszerű leromlását feltételezve azt kapjuk, hogy a hűtéssel el nem vitt teljesítmény a T=3 s időpontig növeli az üzemcnyag hőmérsékletét és ennek eredménye képper. ebben az időpillanatban a lokális hőfluxus nagyobb lesz. Ugyanakkor azonban a hűtőközegbe kevesebb hő jutván a gőztar talom csökkenni fog. Az első effektus csökkenti, a második pedig növeli a DNBR értekét. Nem állapítható meg tehát közvet lenül a közelités hatása az eredményre.
A második lehetséges esetben feltételezzük, hogy a korábbi tel jesítmény továbbra
is eltávozik az üzemanyagból. A hőfluxus
tehát a T=3 s időpontig állandó marad. A lecsökkent forgalmú hűtőközegbe jutó nagymennyiségű hő erősen megnöveli a gőztar talmat. Az eredmény most is hasonló. A lokális hőfluxus kisebb, ez növeli a DNBR-t, míg a gőztartalom-növekedés ellentétesen hat. A változatok részletes vizsgálatát a COBRA kóddal végeztük el. A valóságos eset számítása a 2. ábra görbéivel, mint inputtal készült. Az
első közelítésnél a következőképpen számítottuk
ki a kvázistacionárius futás inputját. A hűtés romlása a hőátadási tényező romlásából adódik, mivel a - áll G ° '
8
A folyadékoldali hőflusus romlása tehát 0.8
Hl , ,*•-!, 'о
о
- 16 -
A
AT= 3 s alatt el nem távolított energia ДЕ = (^o q - q,.) uj F Дх
Itt F a V nagysága.
üzemanyagtérfogathoz tartozó hőátadó felület
u o
А ДЕ energiatöbblet okozta hőmérsékletnövekedés ДЕ
ДТ = < p
Vuo
v
2
* uo
2
A megnevekedett, lokális fluxus kiszámítására közelítőleg írhatjuk, hogy g_
^
дто + дт
^o о
AT
'
mivel a hőmérsékletnövekedés mértéke nem túlzottan nagy. A 2. ábrán T = 3 s-nál a forgalomcsökkenés mértéke 0.815. Ezzel az adattal, valamint (pC )„_ = 3-10~ J/kg°C hőkapacitással p u\J2
a hőmérs-ikletnövekményre ДТ = 55.7 °C adódott, és ezzel a lokális fluxus növekedés a következő: <** - 4
r é g i
* bl2
Ezután a hűtőközegbe jutó, csökkent értékű q - és G= G * 0.815 felhasználásával egyébként változatlan adatokkal stacionárius COBRA számítást készítettünk, majd a D N B R értékét a q* meg növekedett lokális hőfluxussal átszámítottuk, s ezzel a kvázistacionárius közelítés eredményéhez jutottunk. u
min
Q
- 17 -
A második közelítés COBRA inputjainak előállítása nem igényel járulékos számításokat. Minthogy feltételezzük, hogy a 3 з időtartamig fennmaradó teljes teljesítmény a hűtőközegbe átszármaztatódik, így a fűtés a névleges q hőfluxussal törté nik és ugyanez a fluxus veendő figyelembe a DNBR . képzésénél is a lokális hőfluxus alapjául. A hűtőközegnél feltételezésünk szerint pillanatszerű forgalomcsökkenés áll be, tehát a kiin dulási teljesítmény a csökkent forgalmú hűtőközeget fűti. Ebben az esetben a COBRA futtatás tehát q = q hőfluxussal és ^o G = G * 0.815 tömegfluxussal történt. Az eredmények a következő, 1. táblázatban találhatók, ahol a változatok főbb jellemzőit foglaltuk össze. Az 1. táblázat adatai alapján megállapíthatjuk, hogy az a kö zelítő megoldás a konzervatív, amelyik a fűtőelemben tárolandó hő figyelembevételével számítja a DNBR . értékét. Ugyanakkor az eltérés a közelítő megoldások és a tranziens számítás között meglehetősen к ics i. Az eredmények egyébként logikus sorrendet alkotnak. A hűtőközeg be jutó hőmennyiség a legnagyobb a "3" esetben (változatlan teljesítmény mellett). Ennél kevesebb az "1" esetban (a pontos tranziens figyelembevételekor), és legkisebb a "2" közelítés nél, mikor is a hűtéssei el nem vitt hő az üzemanyag hőmérsék letét növeli. A lokális fluxust illetően a sorrend pontosan fordított. Várható volt t~hát, hogy attól függően, hogy a loká lis fluxus, vagy a gőztartalom a meghatározó, valamiféle mono ton változást kapunk. Az eredmények végül is a várakozást iga zolták.
1. táblázat
No.
N folyadékban
q
q
lokális
N(x)
q(x)
G(x)
N-ÄN tárolt
q+Aq
0.815
N(x=0)
q(x=0)
k v
(Kw/in )
DNBRm.m
2883.7
2.52
G(T=0)
3254.5
2.40
0.815 G(x=0)
2874.9
2.53 oo i
1 -
tranziens számítás
2 -
közelítés az akkumulált hő figyelembevételével
3 -
közelítés változatlan teljesítmény feltételezésével
- 19 -
összefoglalás A kutatási jelentésben összefoglalt módszerek és eredmények alkalmazása lehetővé teszi az aktív zóna konzervatív termohidraulikai analízisét. Az eredmények jelentősége abban rej lik, hogy ez a konzervativizmus a leginkább bizonytalan para méterek predikcióját teszi biztonságossá, biztosítva ezáltal azt, hogy az analízis eredményei mindenkor a biztonságos üzem korlátait tegyék ismertté.
- 20 -
Irodalom [1]
Ézsöl Gy. ós Maróti L.: Főkeringető szivattyú kiesés a Paksi Atomerőműэеп. KFKI-1986-87/G.
[2]
Uzembehelyezést megelőző biztonsági jelentés. Paksi Atomerőmű I. blokk. MVMT-ERÖTERV, 1982.
[3]
Szabados L.: A Paksi Atomerőmű biztonsági analízise. OKKFT A/11. MTA KFKI, 1986.
[4]
Helm L.: Ipari folyamatok műszerezése. Műszaki Könyv kiadó, 1966.
[5]
Реакторная установка - Теплогидравлический расчёт актив ной зоны в режимах повышения передача давления на реак торе. У213-ТР-1704.
[6]
Maróti L. és Trösztéi I.: A COCONT kódban alkalmazott modellek leírása. (KFKI riport, mejelenés alatt)
[7]
P.A. Lottes: Nuclear Reactor Heat Transfer. ANL-6469.
The issues of the KFKI preprint/report series are classified as follows: A. Particle and Nuclear Physics В General Relativity and Gravitation C. Cosmic Rays and Space Research D. Fusion and Plasma Physics E. Solid State Physics F. Semiconductor and Bubble Memory Physics and Technology G. Nuclear Reactor Physics and Technology H. Laboratory, Biomedical and Nuclear Reactor Electronics I. Mechanical. Precision Mechanical and Nuclear Engineering J. Analytical and Physical Chemistry K. Health Physics L. Vibration Analysis. CAD, CAM M. Hardware and Software Development Computer Applications. Programming N. Computer Design CAMAC, Compuier Controlled Measurements The complete series or issues discussing one or more of the subjects can be ordered, institutions are kindly requested to contact the KFKI Library individuals the authors
Title and classification of the issues published this year: KFKI-1989-01/D G Kocsis et al A possible method for ion temperature measurement by ion sen sitivé probes KFKI-1989-02/G L Perneczky et al Using the pressurizer spray line in order to minimize loop seal effects мп Hungarian) KFKI-1989-03/E T Csiba et al Propagation of charge density wave voltage noise along a blue bronze, Rbo <МоОз crystal KFKI-1989-04/G G Baranyai et al Experimental investigation of leakage of safety valves by means of acoustic emission detectors (in Hungarian) KFKI-1989-05/A Nguyen Ai Viet et al Can solitons exist in non linear models constructed by the non linear invariance principle 7
KFKI-1989-06/A Nguyen Ai Viet el al A non linearly invariant Skyrme type model KFKI-1989-07/A Nguyen Ai Viet et al Static properties of nucleons in a modified Skyrme model KFKM989-08/B Z Perjós Factor structure of the Tomimatsu Sato metrics KFKI-1989-09/В Z Perjés Unitary spinor methods in general relativity KFKI-1989-10/G G Baranyai et al Reflooding investigations Part I (in Hungarian) KFKI-1989-11/G L Maróti et al Description of the physical models applied in the COCONT code (in Hungarian)
KFKI-1989-12/G L Maróti et al Investigation of the COBRA and COCONT codes by comparison of calculated results, (in Hungarian) KFKI-1989-13/G L. Maróti et al Operational procedure based on hot spot analysis at the WWER 440 type block of Paks Nuclear Power Plant Part III. (In Hungarian)
Kiadja a Központi Fizikai Kutató Intézet Felelős kiadó: Gyimesi Zoltán Szakmai lektor dr Ézsöl Qyórgy Példányszám 88 Törzsszám 89 103 Készült a KFKI sokszorosító üzemében Felelős vezető Qonda Péter Budapest, 1989 február hó !
