> UWT7 5bjbjUU .R7|7|11ltttt 8T / dd:.......$:0 Z2. .ttw.t ..K!+4 a.Xp_ׁ` ->u.$.0/#.>22u. ttttAER Benchmark Solution Sheet
1. Test ID: AER-DYN-003
2. a, Solution Submitted by: Ulrich Rohde, Forschungszentrum Rossendorf,
Institute of Safety Research, Germany
Date: 28.07.2000
b, Reviewed by: P. Siltanen, Fortum Engineering Ltd
Date: 27.10.2000
c, Accepted by: (name)
Date:
3. Code or Program Applied: DYN3D/H1.1
4. Short description of the Code:
The code DYN3D/H1.1 consists of the 3 - dimensional neutron kinetic model and the thermohydraulic module FLOCAL. The neutron kinetics of DYN3D/H1.1 is calculated by using a nodal expansion method (NEM) for hexagonal geometry. The developed method solves the neutron diffusion equation for two energy groups. Stationary state and transient behaviour can be calculated. The code FLOCAL consisting of a two-phase coolant flow model, a fuel rod model and a heat transfer regime map up to superheated steam is coupled with neutron kinetics by the neutron physical constants. One coolant channel per fuel assembly and additional hot channels can be considered.
5. Known Approximations:
Neutron Kinetics
Neutron diffusion theory
Two group theory
Nodewise homogenised cross sections
Thermal Hydraulics
One-dimensional four equation model for two phase coolant flow (momentum equation, energy equation and mass balance of mixture, mass balance of the vapour phase)
Constitutive equations for pressure drop, evaporation and condensation, phase slip
Heat transfer regime map for fuel to coolant heat transfer
Radial heat conduction equation in fuel pin
Feedback
Calculation of neutron cross section by using libraries or input data
6. Mathematical Model:
Neutron Kinetics
The 3-dimensional neutron kinetic model is based on the solution of the 3-dimensional 2group neutron diffusion equation by a nodal expansion method which is specific for the geometry of fuel assemblies [1,2]. It is assumed that the macroscopic cross sections are spatially constant in a node being a part of the hexagonal fuel assembly. The stationary diffusion equation in the node is solved by factorising the space dependency of neutron fluxes in the radial plane and the axial direction. A 2-dimensional diffusion equation in the radial plane and a 1-dimensional equation in axial direction are obtained. The two equations are coupled by the transversal bucklings. In the hexagonal plane the fluxes are expanded by using Bessel functions being the solutions of the Helmholtz equation. The low order coefficients are expressed by the node averaged fluxes and the incoming partial currents averaged over the interface of the hexagon. In this way, the outgoing partial currents at the interfaces are given by the node fluxes and the incoming partial currents. The matrix elements of these relations depend on the transversal buckling and the eigenvalue keff. The 1-dimensional equation in axial direction is solved by a polynomial expansion up to the fourth order. The outgoing partial currents in axial direction are given by the averaged fluxes, incoming partial currents and higher order coefficients. The equations for the 3rd and 4th order polynomials are obtained by Galerkin weighting. The outgoing partial currents at a node interface are the incoming currents in the neighbouring nodes. The steady state diffusion equation is solved by an inner and outer iteration process. The outer iterations are the fission source iterations accelerated by a Chebychev extrapolation scheme. A small number of inner iterations (3 - 5) are sufficient for the convergence. During the outer iteration process the matrix elements are recalculated few times (3-5).
Concerning the time integration over the neutronic time step an implicit difference scheme with exponential transformation is used. The exponents in each node are calculated from the previous time step or during the iteration process. For the calculation of matrix elements describing the relation between partial currents and averaged fluxes it is assumed that the time behaviour of the neutron fluxes in the nodes is exponential and the local variation of the source of delayed neutrons is proportional to the source of prompt neutrons. These assumptions allow the same treatment of diffusion equation in the nodes as in the steady state. In the iteration process we have to solve an inhomogeneous problem. Similar methods as used for the steady state are applied.
Thermal Hydraulics
The thermohydraulic model of the reactor core and the fuel rod model are implemented in the module FLOCAL [3] being a part of DYN3D. The reactor core is modelled by parallel cooling channels which can describe one or more fuel elements. The parallel channels are coupled hydraulically by the condition of equal pressure drop over all core channels. Additionally, so called hot channels can be considered for the investigation of hot spots and uncertainties in power density, coolant temperature or mass flow rate. Thermohydraulic boundary conditions for the core like coolant inlet temperature, pressure, coolant mass flow rate or pressure drop must be given as input for DYN3D.
Mixing of coolant from different loops before entering the core can be modelled by applying several options. Homogeneous mixing can be assumed for each reactor type and number of loops. For VVER-440 type reactors, a mixing model for the downcomer and the lower plenum is implemented in the code. The model is based on the analytical solution of the Navier-Stokes equations in the potential flow approximation in three-dimensional cylindrical geometry and the diffusion equation for heat transport or soluble poison. Turbulent Peclet numbers for the downcomer and the lower plenum are parameters of the model, which are used for a best fit adaptation to experimental results. The mixing model represents an interface between the cold legs of the primary loop and the core inlet.
The module FLOCAL comprises
- a one- or two-phase coolant flow model on the basis of four differential balance equations for mass, energy and momentum of the two-phase mixture and the mass balance for the vapour phase allowing the description of thermodynamic non-equilibrium between the phases,
- a heat transfer regime map from one-phase liquid up to post-critical heat transfer regimes and superheated steam,
- a fuel rod model for the calculation of fuel and cladding temperatures and the determination of some parameters for fuel rod failure estimation.
The two-phase flow model is closed by constitutive laws for heat, mass and momentum transfer, e.g. vapour generation at the heated walls, condensation in the subcooled liquid, phase slip ratio, pressure drop at single flow resistance's and due to friction along the flow channels as well as heat transfer correlations. Different packages of water and steam thermophysical properties presentation can be used.
The heat transfer regime map which is implemented in FLOCAL ranges from one-phase liquid convection up to superheated steam. The occurrence of heat transfer crisis is stated by different correlation's for the critical heat flux. The transient boiling region is described by the KIRCHNER and GRIFFITH interpolation for the heat flux. In the stable post-crisis region for inverted annular or dispersed flow the GROENEVELD - DELORME or a modified BROMLEY correlation's are used. After full evaporation of coolant, heat transfer to superheated steam is estimated by a forced convection correlation [3,4].
Special emphasis is put on adopting the model to the conditions of RIA accidents where a combination of high heat fluxes with high degree of coolant subcooling is typical and thermodynamic non-equilibrium effects are important. A correlation for the Leidenfrost or minimum stable film boiling point TMSFB takes into account the influence of fluid subcooling. In the post-crisis region, a special non-equilibrium correction to the heat transfer coefficient due to subcooled liquid is used.
For the estimation of fuel and cladding temperatures, the heat conduction equation in one-dimensional radial geometry is solved. In the gas gap between fuel and cladding the heat transfer components due to conduction in the gas, radiation and fuel- cladding contact are considered. The idea of the model is that the parameters for the stationary reference state (geometrical gap width , gas pressure and composition) are given and have to be obtained from detailed fuel rod behaviour codes. The changes during the instationary process are estimated by the model.
A thermomechanical model of the fuel and cladding behaviour is implemented into the code. The aim of this model is the estimation of gas gap conductance behaviour for a realistic temperature calculation being important especially for RIA calculations. The thermomechanical model is based on following assumptions:
- one-dimensional modelling of mechanics in radial direction,
- simple superposition of radial thermal, elastic and plastic deformations without axial coupling,
- elastic deformation of the fuel is taken into account only in the case of fuel-cladding contact, plastic deformations of the fuel are not considered,
- cladding is described in the thin shell approximation.
A deterministic model of fuel rod failure during accidents is not included in FLOCAL, but some parameters for the diagnostics of possible fuel rod failure are given, that is:
- fuel enthalpy for each axial node of the rod,
- cladding oxide thickness,
- signalisation of possible cladding rupture, when the cladding stress is positive (inner pressure is larger than outer pressure) and exceeds the yield point.
For testing and validation of the fuel rod and heat transfer model in the code FLOCAL experiments on fuel rod behaviour under RIA conditions were calculated [4].
7. Features of Techniques Used:
Regarding the specification, the code had to be modified for the calculations:
Correlations for thermal conductivity and thermal capacity of fuel and cladding as well as the heat transfer coefficient in the fuel-cladding gas gap given in the benchmark definition are used.
Recommended formulaes for the calculation of pressure losses are used.
Own correlations for heat transfer from fuel to coolant, phase slip ratio, two-phase friction multiplyer and an own evaporation and condensation model are used.
Concerning DNB and post-DNB phenomena, the recommended Gidropress correlation for the DNB ratio and the Groeneveld correlation for film boiling heat transfer are used, but own models for rewetting and cladding oxidation are applied.
The solution provided was calculated with a modified version of the code implementing the modifications described above. Additionally, a solution was obtained using the standard version of the code without modifications [5]. Particularly, the own gas gap model taking into account heat transfer due to gas conduction, radiation and contact conductance was used.
8. Computer, Operational System: SUN workstation, Solaris
9. References:
[1] U. GRUNDMANN and U. ROHDE, "DYN3D/M2 - a Code for Calculation of ReactivityTransients in Cores with Hexagonal Geometry", Proc. of IAEA Technical Committee Meeting on Reactivity Initiated Accidents, Vienna (1989)
Rep. FZR 93-01, Research Centre Rossendorf (1993)
[2] U. GRUNDMANN, "A NEM for Solving Time-Dependent 3-Dimensional Diffusion Equation for Hexagonal Geometry", Proc. International Conference on the Physics of Reactors PHYSOR'90, Marseille (1990)
[3] U. ROHDE, "Modelling of Fuel Rod Behaviour and Heat Transfer in the Code FLOCAL for Reactivity Accident Analysis of Reactor Cores", 1st Baltic Heat Transfer Conference, Gotemborg, (1991), published in: Transport Processes in Engineering, 2: Elsevier Publ., Amsterdam, (1992)
[4] U. ROHDE, The modelling of fuel rod behaviour under RIA conditions in the code DYN3D, 8th Symposium of AER, Bystrice nad Pernstejnem (CZ), September 21-25, 1998, Proc. Pp. 721 - 740
[5] U.GRUNDMANN, U. ROHDE, Comparative Analysis of the Third Three-dimensional Dynamic AER Benchmark Problem with the Help of the Code DYN3D, 5th Symposium of AER, Dobogok (Hungary), October 15-19, 1995, Proc. pp. 329-344
10. Results:
The requested results:
time functions
axial distributions
power distributions
hot channel time functions
hot channel distributions
are included in the ASCII-file dyn003_solfzr.txt. The comparisons with other solutions are presented in the specification document DYN-003.
11. Comparison to Recommended Solution:
No reference solution does exist so far. The comparisons with solutions of other nodal codes are presented in the specification document DYN-003.
PAGE 1
PAGE 5 of NUMPAGES 5
*9U+3AC !<MiV
!
i
.068@ST`"e"R*r*Y-...//11h2j2\3^333Ԥ>*OJQJmH sH OJQJmH sH H*OJQJmH sH H*OJQJmH sH OJQJmH sH OJQJhmH sH >*CJOJQJ5CJCJhnHtHCJOJQJ5CJOJQJ5CJ OJQJ<893BCjk !$L^L`a$$
h8 xHa$
h8 xHd L^L`
/@&
/555!;<MOizwjO$
&F
h8 xHda$$L^L`a$!$
&F
"p
h78 xHd`a$!$
&F
"p
h8 xHd^a$$
&F
p7`a$^
&F
h8 xHd
h8 xHd V
!
g
|V%$
&F
"p
h8 xH7d^7`a$$
&F
h8 xHda$$
h78 xHd`a$"$
&F
#hS8 xHd`a$$
&F
h^a$$
h8 xHda$ g
h
i
=>?@ST#$45$V^`Va$$a$
&F
h8 xHd
h8 xHd>?3!4!##R%S%&&&&0'1'''(((((($7^7`a$$a$ ^`$a$$V^`Va$( )
)))Q*R*t*u***++o,X-Y-$
&F
8 xHda$$
&F
8 xHda$$
h8 xHda$
h8 xHd$a$$7^7`a$Y-.../////#0$0002'
8L".9
EP\hs 0^`0$a$
h8 xHd1
W`0p@P !$`'0*-/2p5@8;=@CPF I h8 xHd
22223333333
4%4?4444$
&Fa$$a$
h8 xHd'
8L".9
EP\hs 0^`0*$
8L".9
EP\hs 0^`0a$3344555555555555555555555555555͗CJCJ0JCJmHnHuj0JCJUCJOJQJhmHnHtHu jCJOJQJUhnHtHCJOJQJhnHtH0JCJ0JmHnHu0J
j0JU55OJQJmH sH OJQJmH sH CJOJQJ444555555555555555555$a$&`#$5$
Y`0p@P !$`'0*-/2p5@8;=@CPF IdOa$$a$5555555555555,&P /R . A!"#$%
i0@0Normal_HmH sH tH <@< Heading 1$@&5CJ OJQJ88 Heading 2$@&CJOJQJLL Heading 3$L@&^L`>*CJOJQJ<< Heading 4$@&>*CJOJQJ Heading 5p$$
Y`0p@P !$`'0*-/2p5@8;=@CPF IdO@&a$>*CJOJQJ<< Heading 6$$dO@&a$CJf@f Heading 76$
h8 xHd@&^CJOJQJ<A@<Default Paragraph Font2B@2 Body TextCJOJQJCBody Text Indentt$
Y`0p@P !$`'0*-/2p5@8;=@CPF I0^`0a$CJOJQJ,,Header
p#, @",Footer
p#&)@1&Page NumberP@BBody Text 2j$
Y`0p@P !$`'0*-/2p5@8;=@CPF IdO1a$CJOJQJ8YR8Document Map-D OJQJ$b$_HmH sH tH QrBody Text 3j$
Y`0p@P !$`'0*-/2p5@8;=@CPF IdOa$CJ`R`Body Text Indent 2:^`:CJOJQJmHsHu@O@VorgabetextCJhmH sH tHuPOPAufzaehl$7^7`a$CJhmH sH tHuHOHText Abstand$a$CJhmH sH tHuLOL
references$
Ya$CJhmH sH tHu
1
1R893BCjk !;<MOizV !ghi=>?@ST#$45>?34R!S!""""0#1###$$$$$$ %
%%%Q&R&t&u&&&''o(X)Y)***+++++#,$,,,....///////
0%0?00000011111111111 hh OO OOh 00@@@0@@@
000335&!g
(Y-2455 !"#$%'(5 (*3!t!T672$t1/`R9X_242$BY+#UwV/42$o#T"42$r+EDk]NUR)w42$y(G:<:,,}42$d;j"oytyt42$tE&woH-O42$Ԫ6pcWFtV+&42$ȟ)/>?.M͘)42$XΰGJsW.42$
eM!DLޣM442$}[bqv8<42$ڻ:G;'#&bܝ 42$ae#洣&,Bb$B42$dK6]Lcm!#nf4@60(
B
S ?1{'',(6(((((}......q/y/1111<=FJ#/L Q
+
_q
rz*+78C%_k !P!""""3#:###$$$$%%'(//////
0
0%0(0?0B01111333333333333333333333333333333333333333rohde;C:\TEMP\AutoWiederherstellen-Speicherung von DYN3D_AER3.asdrohde;C:\TEMP\AutoWiederherstellen-Speicherung von DYN3D_AER3.asdrohde;C:\TEMP\AutoWiederherstellen-Speicherung von DYN3D_AER3.asdrohde;C:\TEMP\AutoWiederherstellen-Speicherung von DYN3D_AER3.asdrohde;C:\TEMP\AutoWiederherstellen-Speicherung von DYN3D_AER3.asdrohde;C:\TEMP\AutoWiederherstellen-Speicherung von DYN3D_AER3.asdrohde.\\pluesch\workd\text\dyn\applic\DYN3D_AER3.docPertti Siltanen+C:\TEMP\AutoRecovery save of DYN3D_AER3.asdPertti SiltanenD:\AERbench\DYN3D_AER3.doc
Imre SzilagyiJ:\aekihtml\dyn003_solfzr.doc(ijg
$a
㚮)"jUEJ 0'!*LG$ܑG I(jc|W(\|+vZo6g/|rT3jRqb:?^=*QGR4eN0:P1j,3Sj]Z F#>n@B+T[oihur}u*.L^vvZ-vvZVD,z
|v{^`OJQJo(-hh^h`OJQJo(-r^`rOJQJo(hh^h`OJQJo(hh^h`OJQJo(hh^h`OJQJo(Vh^`VOJQJo(Vh^`VOJQJo(`h^``OJQJo(Th^T`OJQJo(psp^p`sOJQJo(psp^p`sOJQJo(Vh^`VOJQJo(hh^h`OJQJo(Vh^`VOJQJo(Vh^`VOJQJo(0^`0OJQJo(-hh^h`OJQJo(-Th^T`OJQJo(^`OJQJo(Vh^`VOJQJo(^`OJQJo(hh^h`OJQJo(Vh^`VOJQJo(psp^p`sOJQJo(psp^p`sOJQJo(Th^T`OJQJo(^`OJQJo(eNg
:P)"jT3>nLG$?^=(}uj]Za
+T[oQG|W(
|o6g/-v\|+0'!,3SG I(VD,zqb:UEJ ur.L^v 1111@p0p0ܕp0p01@@UnknownGz Times New Roman5Symbol3&z Arial3z Times5#Tahoma"qJ&KffG&*(W##20d,2t2!AER Benchmark Specification Sheetgrundman
Imre SzilagyiOh+'0 (
DP\
ht|"AER Benchmark Specification SheetosER grundmanmarrunNormal.dotrImre SzilagyiSp3reMicrosoft Word 9.0f@0@*@@5VB@v`*(՜.+,0hp|
FZRW,2
"AER Benchmark Specification SheetTitle
!"#$%&'()+,-./0123456789:;<=>?@ABCEFGHIJKMNOPQRSVRoot Entry F`X1Table*2WordDocument.RSummaryInformation(DDocumentSummaryInformation8LCompObjjObjectPool``
FMicrosoft Word Document
MSWordDocWord.Document.89q