FIRST RESULTS ON THE EFFECT OF FUEL-CLADDING ECCENTRICITY
|
|
- Corey Bryan
- 5 years ago
- Views:
Transcription
1 FIRST RESULTS ON THE EFFECT OF FUEL-CLADDING ECCENTRICITY I. Panka, A. Keresztúri Hungarian Academy of Sciences KFKI Atomic Energy Research Institute, Reactor Analysis Department H-155 Budapest, P.O.BOX 49, Hungary ABSTRACT In the traditional fuel-behaviour or hot channel calculations it is assumed that the fuel pellet is centered within the clad. However, in the real life the pellet could be positioned asymmetrically within the clad, which leads to asymmetric gap conductance and therefore it is worthwhile to investigate the magnitude of the effect on maximal fuel temperature and surface heat flux. In this paper our first experiences are presented on this topic. 1. INTRODUCTION Traditionally, the effect of fuel-cladding eccentricity is neglected or it is not taken into account properly in the fuel-behaviour or hot channel calculation. Therefore we have decided to deal with this problem. It must be mentioned, that only the radial effects are considered, the axial heat conduction is neglected in this investigations. First, we were intended to solve the problem with available traditional codes, like FRAPT-6, FRAPTRAN, but these codes were not able to calculate the problem. Another attempt was to use an older code TOODEE, but there were some problems in adopting this code. Finally we have decided to elaborate a new code. In the first part of the papers we show briefly this code. It follows its verification and some results will be presented.. THE CODE First of all we define the problem to be solved. In Fig. 1, the eccentric geometry can be seen (everything in cm), where R 1 is the clad inner radius, R is the radius of the UO pellet, d is the fuel pellet displacement (the measure of eccentricity, d=0.007 cm, which is considered as a maximal conservative value), R UO is the distance of UO surface from the center of new coordinate system, dgap is the fuel-clad radial gap.
2 R = dgap alfa R UO R 1 = d = dclad = Figure 1: The eccentric fuel-cladding geometry In this case, the azimuthally depending gap size can be calculated as dgap = R - R = R - dcos( α) R d sin ( ). (1) 1 UO 1 θ If the nominal gap size is dgap0 = (R1 - R ) = << R UO, an it is true in our case, then the size of the gap can be approximated by dgap = dgap0 - d * cos( θ ) () A hypothetical (d=0.007 cm) angular depended gap size can be seen in Fig..
3 dgap dgap (cm) Theta (grid points) Figure : The angular dependence of the gap size in case of d=0.007 cm, azimuthally 40 points are considered and here the gap is given in mm. Our task is to solve the time independent heat conduction equation for the entire pin: T r 1 T + r r 1 + r T θ,,, q& + = 0, (3) λ,,, where q& : q.,,, is not zero only in the region of pellet and λ is different for the pellet, for the gap and for the clad and λ is depending on the temperature in general case. It must be mentioned that in the gap region only the heat conduction is considered and we are not dealing with the thermal contact resistance and with the radiation, these effects are firstly neglected. The code solves equation (3) with the following assumptions: We use standard finite difference scheme for the discretization of eq. 3. In any material, many grid points can be given in both direction Generally, the λ coefficients are temperature dependent, and we used VVER related date, but in case of gap we considered only He gas and we used MATPRO data. It is possible to give constant λ values, it is important for verification The heat generation can be given in radial and azimuthal directions The clad outer surface temperature is fixed, in case of verification we used very bigλ for clad (see Chapter 3) The sparse linear system are solved by Gauss-Seidel iteration algorithm There is an other iteration cycle because of the temperature dependence of the λ coefficients. Of course, if we have the temperature field we can calculate other quantities, like heat flux, etc.
4 3. VERIFICATION OF In the literature, we found a theoretical investigation [1] for solving the following problem using some simplifications. These simplifications were: Neglecting the heat conduction in the clad, the temperature at inner surface of the clad is given: T 0 In the pellet region constant λ and heat generation is considered In case of gap: gap conductance is used, but calculating the gap size approximated theta dependence is considered: = h n d 1 d h( d, θ ) 1 + cos( θ ) + cos(θ 1 ( / 0) 0 0 ), (4) d dgap dgap dgap where h n is the nominal gap conductance. Our investigation pointed out that this approximation can lead different results as we used the correct gap size (see later and Fig. 3) /dgap (1/cm) 900 1/dgap (correct) 1/dgap_approximated average 1/dgap nominal (1 Dimensional): 1/dgap Theta (grid points) Figure 3: The angular dependence of the 1/(gap size) in case of d=0.007 cm, azimuthally 40 points are considered, and these values are proportional to the gap conductance values. The theta dependence is taken into account by the following function:
5 ,,, q& F( r, θ ) = T ( r, θ ) T0 1 4λ and r uo R uo n, (5),,, q r F( r, θ ) = & an cos( nθ ) 4λ uo n= 0 R. (6) uo Note, that here R UO is the radius of UO, earlier it was signed as R. For the determination of the a n coefficients, an infinite set of linear equations can be written, and the solution is based on stopping at a finite n. This means, that this solution is quasi theoretical. In the following some results will be given. First, for the one dimensional case, and it follows the two dimensional case. In all cases, two sub-cases are considered: pellet with central hole and pellet without central hole. In case of eccentricity, the results will be given for the cases: using a gap size coming from an approximated gap size (see equation (4)) or using the correct gap size, the results will be different. Unfortunately, for the two dimensional case it is not given formula for the case pellet with central hole. In this case, the results can be compared to the results of the case: using pellet without central hole, but the central hole is taken into account by such a material which have very small lambda and in this region there is no heat generation. In all cases, we used constant parameters:,,, W q & =.13 (relating for UO 3 (R uo ) only), cm W λ gap = , cmk λ = 0. W uo 0369 cmk o T = 317 C, clad inner surface temperature. 0
6 Case 1: One dimension. Pellet without central hole: Theoretical Figure 4: Fuel temperature in Case 1 Case : One dimension. Pellet with central hole: Theoretical Figure 5: Fuel temperature in Case As it can be seen, in 1-dimensional case the agreement is perfect.
7 Case 3: Two-dimensional problem. Pellet without central hole: Figure 6: Angular depended fuel temperature in Case 3, theoretical one Figure 7: Angular depended fuel temperature in Case 3, using approx. gap size
8 Figure 8: Angular depended fuel temperature in Case 3, using correct gap size Figs. 6-8 shows qualitatively good agreements between the calculations, but if we use the correct gap size the maximal and minimal values are different as if we used the approximated gap size (see Eq.( 4) and Fig. 3) Table 1: Temperatures of the fuel pellet, without central hole T, center T, outer surface of UO -dim, approx. gap size, theoretical -dim, approx. gap size, -dim, correct gap size, Max. Min. Max. Min (+0. %) (+5.5 %) 1-dim (+7.4 %) 75.1 (+0. %) (+5.5 %) (+7.4 %) (-.4 %) (-6. %) (+3.7 %) (-1.6 %) (+14.6 %) (-4.1 %) This table shows the maximal and the minimal temperature values at the fuel center and at the outer surface of the UO pellet. From the point of view of verification the bold one is important: at the center there is good agreement, but at the outer surface the difference is approx. %. This value is not large, and this difference can be explained by a various reason (e.g. numerical approximations, etc). However, it can be a possible reason that the theta dependence is not taken into account properly by Eq. (5). The problem is that the theta dependence of the gap size is very strange if we use approximated gap size given in Eq. (4), see Fig 3. Because in Eq (5) only cosine functions are used it is not enough to describe the theta dependence, consequently we are even needed to use sinus function.
9 Note, that we performed calculations for the case if cos(θ) is not considered in Eq. (5), and in this case the differences were below 0.3 %. The other results (given in Table 1) show how different values can be get using different approximations. E.g. in case of using the correct gap size the maximal and minimal temperate values are different as if we used approximated gap size. The maximal temperature of fuel is smaller in all cases than the values coming from the one-dimensional case. Another observation can be that the angular dependence of the temperature at UO outer surface is very big, which means more than 60 % difference to the average in heat flux. Case 4: Two-dimensional problem. Pellet with central hole: Figure 9: Angular depended fuel temperature in Case 4, using correct gap size Figure 10: Angular depended fuel temperature in Case 4, using approx. gap size
10 Figure 11: Angular depended fuel temperature in Case 4, using approx. gap size, and a material with very small lambda and there is no heat generation in the central hole Table : Temperatures of the fuel pellet, with central hole T, center T, outer surface of UO -dim, approx. gap size, -dim, approx. gap size, hole filled with material having small lambda -dim, correct gap size, Max. Min. Max. Min (+0.01 %) (+7. %) 1-dim (+4.1 %) (-0.03 %) (+3.7 %) (+11.1 %) (+0.0 %) (+16.5 %) (-.6 %) 356. (+0.0 %) 34.7 (-3.8 %) (+35.3 %) In this case, similar observations can be get as in case 3, but now we were not able to used theoretical considerations. The verification is based on the results of the previous case. 3. RESULTS FOR THE ENTIRE PIN In this section the entire fuel pin will be calculated. First we use constant values and than we show the temperature dependent caser. Firstly, λ clad =0.199 W/cmK T0= o C (at outer surface of clad)
11 Case 5: Two-dimensional problem. Pellet with central hole, constant lambda values, calculation of the full fuel pin Figure 1: Angular depended and 1-dimensional fuel temperatures in Case 5, using correct gap size Relative heat flux (-) Theta (grid points) Figure 13: Angular depended relative heat flux in Case 5, using correct gap size.
12 Case 6: Two-dimensional problem. Pellet with central hole, temperature dependent lambda values, calculation of the full fuel pin Figure 14: Angular depended fuel temperature in Case 6, using correct gap size Relative heat flux (-) Heat flux, iter No. 1 Heat flux, iter No. 9 Heat flux, iter No. 10 Constant lambda values Theta (grid points) Figure 15: Angular depended relative heat flux in Case 6, using correct gap size. From the point of view of DNBR calculations, the most important results can be seen in Fig. 15. This figure shows the angular dependent relative heat flux at the clad outer surface can be more than 1.5 (approx. 1.7 in this case) corresponding to the average relative heat flux. Additionally, the maximal fuel temperature is grater than in 1-dimensional case, but tis difference is not so big: approx. %, see Fig 14.
13 4. SUMMARY A new code, was presented with its verification and some results were showed. The results pointed out that the angular dependence of the gap size can be important, especially in case of the heat flux calculation at the outer surface of the clad: the angular dependent relative heat flux at the clad outer surface can be more than 1.5 (approx. 1.7 in our case) corresponding to the average relative heat flux. In the future, we are going to continue this work. We plan to take into account the uncertainties of the angular dependence of the gap sizes, considering a hydraulic module, etc. Another task can be the extension of this work for the time dependent case, which can be important in case of Reactivity Initiated Transients. REFERENCES [1] O. McNARY and T.H. BAUER, The effect of asymmetric fuel-clad gap conductance on fuel pin thermal performance, NED, 63 (1981) 39-46
Coupling of thermal-mechanics and thermalhydraulics codes for the hot channel analysis of RIA events First steps in AEKI toward multiphysics
Coupling of thermal-mechanics and thermalhydraulics codes for the hot channel analysis of RIA events First steps in AEKI toward multiphysics A. Keresztúri, I. Panka, A. Molnár KFKI Atomic Energy Research
More informationAPPLICATION OF THE COUPLED THREE DIMENSIONAL THERMAL- HYDRAULICS AND NEUTRON KINETICS MODELS TO PWR STEAM LINE BREAK ANALYSIS
APPLICATION OF THE COUPLED THREE DIMENSIONAL THERMAL- HYDRAULICS AND NEUTRON KINETICS MODELS TO PWR STEAM LINE BREAK ANALYSIS Michel GONNET and Michel CANAC FRAMATOME Tour Framatome. Cedex 16, Paris-La
More informationSUB-CHAPTER D.1. SUMMARY DESCRIPTION
PAGE : 1 / 12 CHAPTER D. REACTOR AND CORE SUB-CHAPTER D.1. SUMMARY DESCRIPTION Chapter D describes the nuclear, hydraulic and thermal characteristics of the reactor, the proposals made at the present stage
More informationX. Neutron and Power Distribution
X. Neutron and Power Distribution X.1. Distribution of the Neutron Flux in the Reactor In order for the power generated by the fission reactions to be maintained at a constant level, the fission rate must
More informationE. not enough information given to decide
Q22.1 A spherical Gaussian surface (#1) encloses and is centered on a point charge +q. A second spherical Gaussian surface (#2) of the same size also encloses the charge but is not centered on it. Compared
More informationA NEW LUMPED THERMAL RESISTANCE HEAT TRANSFER MODEL FOR FUEL PIN STRUCTURE
A NEW LUMPED THERMAL RESISTANCE HEAT TRANSFER MODEL FOR FUEL PIN STRUCTURE Shisheng Wang, Andrei Rineiski, Liancheng Guo Institute for Nuclear and Energy Technologies (IKET) Karlsruhe Institute of Technology
More information22.06 ENGINEERING OF NUCLEAR SYSTEMS OPEN BOOK FINAL EXAM 3 HOURS
22.6 ENGINEERING OF NUCLEAR SYSTEMS OPEN BOOK FINAL EXAM 3 HOURS Short Questions (1% each) a) The specific power in a UO 2 pellet of a certain LWR is q"'=2 W/cm 3. The fuel 235 U enrichment is 4 % by weight.
More informationTHREE-DIMENSIONAL INTEGRAL NEUTRON TRANSPORT CELL CALCULATIONS FOR THE DETERMINATION OF MEAN CELL CROSS SECTIONS
THREE-DIMENSIONAL INTEGRAL NEUTRON TRANSPORT CELL CALCULATIONS FOR THE DETERMINATION OF MEAN CELL CROSS SECTIONS Carsten Beckert 1. Introduction To calculate the neutron transport in a reactor, it is often
More informationComparison of Silicon Carbide and Zircaloy4 Cladding during LBLOCA
Comparison of Silicon Carbide and Zircaloy4 Cladding during LBLOCA Prepared By: Kwangwon Ahn Prepared For: 22.314 Prepared On: December 7 th, 2006 Comparison of Silicon Carbide and Zircaloy4 Cladding during
More informationDesign constraints Maximum clad temperature, linear power rating
Design constraints Maximum clad temperature, linear power rating K.S. Rajan Professor, School of Chemical & Biotechnology SASTRA University Joint Initiative of IITs and IISc Funded by MHRD Page 1 of 7
More informationCALCULATION OF TEMPERATURE REACTIVITY COEFFICIENTS IN KRITZ-2 CRITICAL EXPERIMENTS USING WIMS ABSTRACT
CALCULATION OF TEMPERATURE REACTIVITY COEFFICIENTS IN KRITZ-2 CRITICAL EXPERIMENTS USING WIMS D J Powney AEA Technology, Nuclear Science, Winfrith Technology Centre, Dorchester, Dorset DT2 8DH United Kingdom
More informationA PWR HOT-ROD MODEL: RELAP5/MOD3.2.2Y AS A SUBCHANNEL CODE I.C. KIRSTEN (1), G.R. KIMBER (2), R. PAGE (3), J.R. JONES (1) ABSTRACT
FR0200515 9 lh International Conference on Nuclear Engineering, ICONE-9 8-12 April 2001, Nice, France A PWR HOT-ROD MODEL: RELAP5/MOD3.2.2Y AS A SUBCHANNEL CODE I.C. KIRSTEN (1), G.R. KIMBER (2), R. PAGE
More informationNEW FUEL IN MARIA RESEARCH REACTOR, PROVIDING BETTER CONDITIONS FOR IRRADIATION IN THE FAST NEUTRON SPECTRUM.
NEW FUEL IN MARIA RESEARCH REACTOR, PROVIDING BETTER CONDITIONS FOR IRRADIATION IN THE FAST NEUTRON SPECTRUM. M. LIPKA National Centre for Nuclear Research Andrzeja Sołtana 7, 05-400 Otwock-Świerk, Poland
More informationStudy of Control rod worth in the TMSR
Nuclear Science and Techniques 24 (2013) 010601 Study of Control rod worth in the TMSR ZHOU Xuemei * LIU Guimin 1 Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China
More informationULOF Accident Analysis for 300 MWt Pb-Bi Coolled MOX Fuelled SPINNOR Reactor
ULOF Accident Analysis for 300 MWt Pb-Bi Coolled MOX Fuelled SPINNOR Reactor Ade afar Abdullah Electrical Engineering Department, Faculty of Technology and Vocational Education Indonesia University of
More informationFuel - Coolant Heat Transfer
Heat Transfer 5-1 Chapter 5 Fuel - Coolant Heat Transfer 5.1 Introduction The interface between the fuel and the coolant is centrally important to reactor design since it is here that the limit to power
More informationPOST SERVICE INVESTIGATIONS OF VVER-440 RPV STEEL FROM NPP GREIFSWALD
POST SERVICE INVESTIGATIONS OF VVER-440 RPV STEEL FROM NPP GREIFSWALD 1. Introduction Bertram Boehmer, Juergen Boehmert, and Udo Rindelhardt Pressure vessel integrity assessment after long-term service
More informationPhD Qualifying Exam Nuclear Engineering Program. Part 1 Core Courses
PhD Qualifying Exam Nuclear Engineering Program Part 1 Core Courses 9:00 am 12:00 noon, November 19, 2016 (1) Nuclear Reactor Analysis During the startup of a one-region, homogeneous slab reactor of size
More informationFuel BurnupCalculations and Uncertainties
Fuel BurnupCalculations and Uncertainties Outline Review lattice physics methods Different approaches to burnup predictions Linkage to fuel safety criteria Sources of uncertainty Survey of available codes
More informationTHERMAL HYDRAULIC REACTOR CORE CALCULATIONS BASED ON COUPLING THE CFD CODE ANSYS CFX WITH THE 3D NEUTRON KINETIC CORE MODEL DYN3D
THERMAL HYDRAULIC REACTOR CORE CALCULATIONS BASED ON COUPLING THE CFD CODE ANSYS CFX WITH THE 3D NEUTRON KINETIC CORE MODEL DYN3D A. Grahn, S. Kliem, U. Rohde Forschungszentrum Dresden-Rossendorf, Institute
More informationKord Smith Art DiGiovine Dan Hagrman Scott Palmtag. Studsvik Scandpower. CASMO User s Group May 2003
Kord Smith Art DiGiovine Dan Hagrman Scott Palmtag CASMO User s Group May 2003 TFU-related data is required input for: -CASMO-4 - SIMULATE-3 - SIMULATE-3K and SIMULATE-3R (implicit in XIMAGE and GARDEL)
More informationIdaho National Laboratory Reactor Analysis Applications of the Serpent Lattice Physics Code
Idaho National Laboratory Reactor Analysis Applications of the Serpent Lattice Physics Code By Frederick N. Gleicher II, Javier Ortensi, Benjamin Baker, and Mark DeHart Outline Intra-Pin Power and Flux
More informationDETAILED CORE DESIGN AND FLOW COOLANT CONDITIONS FOR NEUTRON FLUX MAXIMIZATION IN RESEARCH REACTORS
Proceedings of ICONE14 14th International Conference on Nuclear Engineering June 17-21, 2006, Miami, FL, USA ICONE14-89547 DETAILED CORE DESIGN AND FLOW COOLANT CONDITIONS FOR NEUTRON FLUX MAXIMIZATION
More informationModernization of Cross Section Library for VVER-1000 Type Reactors Internals and Pressure Vessel Dosimetry
EPJ Web of Conferences 106, 04014 (2016) DOI: 10.1051/epjconf/201610604014 C Owned by the authors, published by EDP Sciences, 2016 Modernization of Cross Section Library for VVER-1000 Type Reactors Internals
More informationCritical Experiment Analyses by CHAPLET-3D Code in Two- and Three-Dimensional Core Models
Journal of NUCLEAR SCIENCE and TECHNOLOGY, Vol. 42, No. 1, p. 101 108 (January 2005) TECHNICAL REPORT Critical Experiment Analyses by CHAPLET-3D Code in Two- and Three-Dimensional Core Models Shinya KOSAKA
More informationVERIFICATION OF A REACTOR PHYSICS CALCULATION SCHEME FOR THE CROCUS REACTOR. Paul Scherrer Institut (PSI) CH-5232 Villigen-PSI 2
VERIFICATION OF A REACTOR PHYSICS CALCULATION SCHEME FOR THE CROCUS REACTOR M. Hursin 1,*, D. Siefman 2, A. Rais 2, G. Girardin 2 and A. Pautz 1,2 1 Paul Scherrer Institut (PSI) CH-5232 Villigen-PSI 2
More information(Refer Slide Time: 01:17)
Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay Lecture No. 7 Heat Conduction 4 Today we are going to look at some one dimensional
More informationCurrent Developments of the VVER Core Analysis Code KARATE-440
Current Developments of the VVER Core Analysis Code KARATE-440 György Hegyi Hungarian Academy of Sciences Centre for Energy Research, Budapest, Hungary Reactor Analysis Department Konkoly Thege Miklós
More informationFuel Element Burnup Determination in HEU - LEU Mixed TRIGA Research Reactor Core
Fuel Element Burnup Determination in HEU - LEU Mixed TRIGA Research Reactor Core Tomaž Žagar, Matjaž Ravnik Institute "Jožef Stefan", Jamova 39, Ljubljana, Slovenia Tomaz.Zagar@ijs.si Abstract This paper
More informationExtension of the MCBEND Monte Carlo Code to Perform Adjoint Calculations using Point Energy Data
Extension of the MCBEND Monte Carlo Code to Perform Adjoint Calculations using Point Energy Data Malcolm Grimstone Abstract In radiation transport calculations there are many situations where the adjoint
More informationNumerical simulations of hardening and cooling of the early-age concrete
Numerical simulations of hardening and cooling of the early-age concrete Grzegorz Knor IPPT PAN, 5b Pawińskiego str., 02-106 Warsaw, Poland gknor@ippt.pan.pl received 26 September 2013, in final form 30
More informationModeling and Identification of a Nuclear Reactor with Temperature Effects and Xenon Poisoning
Modeling and Identification of a uclear Reactor with Temperature Effects and Xenon Poisoning Attila Gábor, Csaba Fazekas, Gábor Szederkényi and Katalin M. Hangos Process Control Research Group, Computer
More informationMAGNETIC MEASUREMENTS ON THE GUN FOCUSING SOLENOID
SPARC-MM-05/003 20 September 2005 MAGNETIC MEASUREMENTS ON THE GUN FOCUSING SOLENOID B. Bolli, S. Ceravolo,M. Esposito, F. Iungo, M. Paris, M.Preger, C. Sanelli, F.Sardone, F. Sgamma, M. Troiani, P.Musumeci,
More informationModified lumped model for Transient heat conduction in spherical shape
American International Journal of Research in Science, Technology, Engineering & Mathematics Available online at http://www.iasir.net ISSN (Print): 2328-3491, ISSN (Online): 2328-3580, ISSN (CD-ROM): 2328-3629
More informationDOPPLER COEFFICIENT OF REACTIVITY BENCHMARK CALCULATIONS FOR DIFFERENT ENRICHMENTS OF UO 2
Supercomputing in Nuclear Applications (M&C + SNA 2007) Monterey, California, April 15-19, 2007, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2007) DOPPLER COEFFICIENT OF REACTIVITY BENCHMARK
More information7 Mathematical Methods 7.6 Insulation (10 units)
7 Mathematical Methods 7.6 Insulation (10 units) There are no prerequisites for this project. 1 Introduction When sheets of plastic and of other insulating materials are used in the construction of building
More information«CALCULATION OF ISOTOPE BURN-UP AND CHANGE IN EFFICIENCY OF ABSORBING ELEMENTS OF WWER-1000 CONTROL AND PROTECTION SYSTEM DURING BURN-UP».
«CALCULATION OF ISOTOPE BURN-UP AND CHANGE IN EFFICIENCY OF ABSORBING ELEMENTS OF WWER-1000 CONTROL AND PROTECTION SYSTEM DURING BURN-UP». O.A. Timofeeva, K.U. Kurakin FSUE EDO «GIDROPRESS», Podolsk, Russia
More informationBERYLLIUM IMPREGNATION OF URANIUM FUEL: THERMAL MODELING OF CYLINDRICAL OBJECTS FOR EFFICIENCY EVALUATION
BERYLLIUM IMPREGNATION OF URANIUM FUEL: THERMAL MODELING OF CYLINDRICAL OBJECTS FOR EFFICIENCY EVALUATION A Senior Scholars Thesis by NICHOLAS MORGAN LYNN Submitted to the Office of Undergraduate Research
More informationdt Now we will look at the E&M force on moving charges to explore the momentum conservation law in E&M.
. Momentum Conservation.. Momentum in mechanics In classical mechanics p = m v and nd Newton s law d p F = dt If m is constant with time d v F = m = m a dt Now we will look at the &M force on moving charges
More informationMechanics, Heat, Oscillations and Waves Prof. V. Balakrishnan Department of Physics Indian Institute of Technology, Madras
Mechanics, Heat, Oscillations and Waves Prof. V. Balakrishnan Department of Physics Indian Institute of Technology, Madras Lecture 08 Vectors in a Plane, Scalars & Pseudoscalers Let us continue today with
More informationLectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 7
ectures on Nuclear Power Safety ecture No 7 itle: hermal-hydraulic nalysis of Single-Phase lows in Heated hannels Department of Energy echnology KH Spring 005 Slide No Outline of the ecture lad-oolant
More informationCONFORMITY BETWEEN LR0 MOCK UPS AND VVERS NPP PRV ATTENUATION
CONFORMITY BETWEEN LR MOCK UPS AND VVERS NPP PRV ATTENUATION D. Kirilova, K. Ilieva, S. Belousov Institute for Nuclear Research and Nuclear Energy, Bulgaria Email address of main author: desi.kirilova@gmail.com
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science Electric Machines
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.685 Electric Machines Problem Set 10 Issued November 11, 2013 Due November 20, 2013 Problem 1: Permanent
More informationROTATING RING. Volume of small element = Rdθbt if weight density of ring = ρ weight of small element = ρrbtdθ. Figure 1 Rotating ring
ROTATIONAL STRESSES INTRODUCTION High centrifugal forces are developed in machine components rotating at a high angular speed of the order of 100 to 500 revolutions per second (rps). High centrifugal force
More informationEngineering Notebook
Engineering Notebook miscellaneous problems and solutions compiled by C Bond Vol1 no9 Actuator with zero bend angle Rotary Actuator Equations Given two points at which the track radius and gap skew angle
More informationOECD-PBMR Benchmark Thermal Hydraulics J.B.M. de Haas (NRG) The Netherlands. OECD Paris June, 2005
OECD-PBMR Benchmark Thermal Hydraulics J.B.M. de Haas (NRG) The Netherlands OECD Paris 16-17 June, 2005 PBMR Codes and Methods At NRG the PBMR has been modelled in the reactor code PANTHERMIX, a combination
More informationGrid supports design for dual-cooled fuel rods
Int. Jnl. of Multiphysics Volume 5 Number 4 011 99 Grid supports design for dual-cooled fuel rods JaeYong Kim, YoungHo Lee, KyungHo Yoon and HyungKyu Kim Innovative Nuclear Fuel Division, Korea Atomic
More informationConduction and radiation Prof. C. Balaji Department of Mechanical Engineering Indian Institute of Technology, Madras
Conduction and radiation Prof. C. Balaji Department of Mechanical Engineering Indian Institute of Technology, Madras Module No. #01 Lecture No. #03 Solid angle, spectral radiation intensity In yesterday
More informationLecture D Steady State Heat Conduction in Cylindrical Geometry
Conduction and Convection Heat Transfer Prof. S.K. Som Prof. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology Kharagpur Lecture - 08 1D Steady State Heat Conduction
More informationChapter 5. Effects of Photonic Crystal Band Gap on Rotation and Deformation of Hollow Te Rods in Triangular Lattice
Chapter 5 Effects of Photonic Crystal Band Gap on Rotation and Deformation of Hollow Te Rods in Triangular Lattice In chapter 3 and 4, we have demonstrated that the deformed rods, rotational rods and perturbation
More informationBiaxial Analysis of General Shaped Base Plates
Biaxial Analysis of General Shaped Base Plates R. GONZALO ORELLANA 1 Summary: A linear model is used for the contact stresses calculation between a steel base plate and a concrete foundation. It is also
More informationCHAPTER 8 CONSERVATION LAWS
CHAPTER 8 CONSERVATION LAWS Outlines 1. Charge and Energy 2. The Poynting s Theorem 3. Momentum 4. Angular Momentum 2 Conservation of charge and energy The net amount of charges in a volume V is given
More informationThere are four questions on this exam. All are of equal value but not necessarily of equal difficulty. Answer all four questions
Electricity and Magnetism (PHYS2016) Second Semester Final Exam 2015 There are four questions on this exam. All are of equal value but not necessarily of equal difficulty. Answer all four questions You
More informationFigure 22.1 Unflattened Flux Distribution
22 Neutron Flux Control If nothing were done to flatten the flux in our reactors, it would look something like Figure 22.1. The flux would be a maximum in the efuel of the reactor (where neutrons are moving
More informationComparison of assessment of neutron fluence affecting VVER 440 reactor pressure vessel using DORT and TORT codes
Comparison of assessment of neutron fluence affecting VVER 440 reactor pressure vessel using DORT and TORT codes P. Montero Department of Neutronics, Research Center Rez, Cz International Conference on
More informationAnalysis of a Double Pipe Heat Exchanger Performance by Use of Porous Baffles and Nanofluids
Analysis of a Double Pipe Heat Exchanger Performance by Use of Porous Baffles and Nanofluids N. Targui, H. Kahalerras Abstract The present work is a numerical simulation of nanofluids flow in a double
More informationStudy of boron dilution phenomenon in the core and fuel assemblies of Bushehr VVER-1000 reactor in normal operating conditions
NUKLEONIKA 2010;55(3:323 330 ORIGINAL PAPER Study of boron dilution phenomenon in the core and fuel assemblies of Bushehr VVER-1000 reactor in normal operating conditions Yashar Rahmani, Ehsan Zarifi,
More informationANALYSIS OF THE OECD MSLB BENCHMARK WITH THE COUPLED NEUTRONIC AND THERMAL-HYDRAULICS CODE RELAP5/PARCS
ANALYSIS OF THE OECD MSLB BENCHMARK WITH THE COUPLED NEUTRONIC AND THERMAL-HYDRAULICS CODE RELAP5/PARCS T. Kozlowski, R. M. Miller, T. Downar School of Nuclear Engineering Purdue University United States
More informationNEUTRONIC ANALYSIS STUDIES OF THE SPALLATION TARGET WINDOW FOR A GAS COOLED ADS CONCEPT.
NEUTRONIC ANALYSIS STUDIES OF THE SPALLATION TARGET WINDOW FOR A GAS COOLED ADS CONCEPT. A. Abánades, A. Blanco, A. Burgos, S. Cuesta, P.T. León, J. M. Martínez-Val, M. Perlado Universidad Politecnica
More informationNEUTRON AND GAMMA FLUENCE AND RADIATION DAMAGE PARAMETERS OF EX-CORE COMPONENTS OF RUSSIAN AND GERMAN LIGHT WATER REACTORS
NEUTRON AND GAMMA FLUENCE AND RADIATION DAMAGE PARAMETERS OF EX-CORE COMPONENTS OF RUSSIAN AND GERMAN LIGHT WATER REACTORS Bertram Boehmer, Joerg Konheiser, Klaus Noack, Anatolij Rogov, Gennady Borodkin
More informationChapter 22 Gauss s Law. Copyright 2009 Pearson Education, Inc.
Chapter 22 Gauss s Law 22-1 Electric Flux Electric flux: Electric flux through an area is proportional to the total number of field lines crossing the area. 22-1 Electric Flux Example 22-1: Electric flux.
More informationREFLECTOR FEEDBACK COEFFICIENT: NEUTRONIC CONSIDERATIONS IN MTR-TYPE RESEARCH REACTORS ABSTRACT
REFLECTOR FEEDBACK COEFFICIENT: NEUTRONIC CONSIDERATIONS IN MTR-TYPE RESEARCH REACTORS L. MANIFACIER, L. CHABERT, M. BOYARD TechnicAtome, CS 50497, 13593 Aix-en-Provence Cedex 3, France ABSTRACT Having
More informationANALYTICAL COMPUTATION OF RELUCTANCE SYN- CHRONOUS MACHINE INDUCTANCES UNDER DIF- FERENT ECCENTRICITY FAULTS
Progress In Electromagnetics Research M, Vol. 24, 29 44, 2012 ANALYTICAL COMPUTATION OF RELUCTANCE SYN- CHRONOUS MACHINE INDUCTANCES UNDER DIF- FERENT ECCENTRICITY FAULTS H. Akbari * Department of Electrical
More informationSteam Generator Tubing Inspection
Steam Generator Tubing Inspection Analytical Determination of Critical Flaw Dimensions in Steam Generator Tubing I. Kadenko, N. Sakhno, R. Yermolenko, Nondestructive Examination Training and Certification
More informationData analysis with uncertainty evaluation for the Ignalina NPP RBMK 1500 gas-gap closure evaluation
Data analysis with uncertainty evaluation for the Ignalina NPP RBMK 1500 gas-gap closure evaluation M. Liaukonisl &J. Augutis2 Department ofmathematics and Statistics, Vytautas Magnus University, Lithuania.
More informationHydrogen and Helium Edge-Plasmas
Hydrogen and Helium Edge-Plasmas Comparison of high and low recycling T.D. Rognlien and M.E. Rensink Lawrence Livermore National Lab Presented at the ALPS/APEX Meeting Argonne National Lab May 8-12, 2
More informationSIMULATION OF LEAKING FUEL RODS
SIMULATION OF LEAKING FUEL RODS Zoltán Hózer KFKI Atomic Energy Research Institute P.O.B. 49, H-1525 Budapest, Hungary Hozer@sunserv.kfki.hu Abstract The behaviour of failed fuel rods includes several
More informationTask 3 Desired Stakeholder Outcomes
Task 3 Desired Stakeholder Outcomes Colby Jensen IRP Kickoff Meeting, Nov 19-20, 2015 Instrumentation Overview Three general levels of core instrumentation: Reactor control and operation Additional reactor
More informationCalculation of Spatial Weighting Functions for Ex-Core Detectors of VVER-440 Reactors by Monte Carlo Method
International Conference Nuclear Energy for New Europe 2003 Portorož, Slovenia, September 8-11, 2003 http://www.drustvo-js.si/port2003 Calculation of Spatial Weighting Functions for Ex-Core Detectors of
More informationChapter 2 HEAT CONDUCTION EQUATION
Heat and Mass Transfer: Fundamentals & Applications Fourth Edition Yunus A. Cengel, Afshin J. Ghajar McGraw-Hill, 2011 Chapter 2 HEAT CONDUCTION EQUATION Mehmet Kanoglu University of Gaziantep Copyright
More informationLesson 9: Multiplying Media (Reactors)
Lesson 9: Multiplying Media (Reactors) Laboratory for Reactor Physics and Systems Behaviour Multiplication Factors Reactor Equation for a Bare, Homogeneous Reactor Geometrical, Material Buckling Spherical,
More informationSodium void coefficient map by Serpent
Wir schaffen Wissen heute für morgen Paul Scheer Institut 4th Annual Serpent Users Group Meeting (Cambridge Sept 17-19, 2014): Sandro Pelloni Sodium void coefficient map by Serpent PSI, 8. September 2014
More informationCFX SIMULATION OF A HORIZONTAL HEATER RODS TEST
CFX SIMULATION OF A HORIZONTAL HEATER RODS TEST Hyoung Tae Kim, Bo Wook Rhee, Joo Hwan Park Korea Atomic Energy Research Institute 150 Dukjin-Dong, Yusong-Gu, Daejon 305-353, Korea kht@kaeri.re.kr Abstract
More informationBurn Stabilization of a Tokamak Power Plant with On-Line Estimations of Energy and Particle Confinement Times
Burn Stabilization of a Tokamak Power Plant with On-Line Estimations of Energy and Particle Confinement Times Javier E. Vitela vitela@nucleares.unam.mx Instituto de Ciencias Nucleares Universidad Nacional
More informationLight source thermal analysis I Sodium high-pressure discharge Lamp
1 Portál pre odborné publikovanie ISSN 1338-0087 Light source thermal analysis I Sodium high-pressure discharge Lamp Fric Róbert Elektrotechnika 24.11.2010 This is the first chapter from the series of
More informationSimple examples of MHD equilibria
Department of Physics Seminar. grade: Nuclear engineering Simple examples of MHD equilibria Author: Ingrid Vavtar Mentor: prof. ddr. Tomaž Gyergyek Ljubljana, 017 Summary: In this seminar paper I will
More information1. INTRODUCTION 2. EAEA EXISTING CAPABILITIES AND FACILITIES
EGYPT FINAL REPORT FOR THE CRP ON DEVELOPING TECHNIQUES FOR SMALL- SCALE, INDIGENOUS PRODUCTION OF MO-99 USING LOW- ENRICHED URANIUM (LEU) OR NEUTRON ACTIVATION 1. INTRODUCTION The Egypt country report
More informationDEVELOPMENT OF HIGH-FIDELITY MULTI- PHYSICS SYSTEM FOR LIGHT WATER REACTOR ANALYSIS
The Pennsylvania State University The Graduate School College of Engineering DEVELOPMENT OF HIGH-FIDELITY MULTI- PHYSICS SYSTEM FOR LIGHT WATER REACTOR ANALYSIS A Dissertation in Nuclear Engineering by
More informationAdvanced Heavy Water Reactor. Amit Thakur Reactor Physics Design Division Bhabha Atomic Research Centre, INDIA
Advanced Heavy Water Reactor Amit Thakur Reactor Physics Design Division Bhabha Atomic Research Centre, INDIA Design objectives of AHWR The Advanced Heavy Water Reactor (AHWR) is a unique reactor designed
More informationABSTRACT 1 INTRODUCTION
A NODAL SP 3 APPROACH FOR REACTORS WITH HEXAGONAL FUEL ASSEMBLIES S. Duerigen, U. Grundmann, S. Mittag, B. Merk, S. Kliem Forschungszentrum Dresden-Rossendorf e.v. Institute of Safety Research P.O. Box
More informationUSA HTR NEUTRONIC CHARACTERIZATION OF THE SAFARI-1 MATERIAL TESTING REACTOR
Proceedings of HTR2008 4 th International Topical Meeting on High Temperature Reactors September 28-October 1, 2008, Washington, D.C, USA HTR2008-58155 NEUTRONIC CHARACTERIZATION OF THE SAFARI-1 MATERIAL
More informationA FINITE VOLUME-BASED NETWORK METHOD FOR THE PREDICTION OF HEAT, MASS AND MOMENTUM TRANSFER IN A PEBBLE BED REACTOR
A FINITE VOLUME-BASED NETWORK METHOD FOR THE PREDICTION OF HEAT, MASS AND MOMENTUM TRANSFER IN A PEBBLE BED REACTOR GP Greyvenstein and HJ van Antwerpen Energy Systems Research North-West University, Private
More informationThis paper is about a new method for analyzing the distribution property of nonuniform stochastic samples.
This paper is about a new method for analyzing the distribution property of nonuniform stochastic samples. Stochastic sampling is a fundamental component in many graphics applications. Examples include
More informationGauss s Law. Chapter 22. Electric Flux Gauss s Law: Definition. Applications of Gauss s Law
Electric Flux Gauss s Law: Definition Chapter 22 Gauss s Law Applications of Gauss s Law Uniform Charged Sphere Infinite Line of Charge Infinite Sheet of Charge Two infinite sheets of charge Phys 2435:
More informationExperimental Studies of Active Temperature Control in Solid Breeder Blankets
Experimental Studies of Active Temperature Control in Solid Breeder Blankets M. S. Tillack, A. R. Raffray, A. Y. Ying, M. A. Abdou, and P. Huemer Mechanical, Aerospace and Nuclear Engineering Department
More informationCoupled Neutronics Thermalhydraulics LC)CA Analysis
Coupled Neutronics Thermalhydraulics LC)CA Analysis B.Rouben, Manager Reactor Core Physics Branch,AE,CL Presented at Chulalongkorn University Bangkok, Thailand 9 1997 December RFSP -Reactor Fuelling Simulation
More informationResearch Article Analysis of NEA-NSC PWR Uncontrolled Control Rod Withdrawal at Zero Power Benchmark Cases with NODAL3 Code
Hindawi Science and Technology of Nuclear Installations Volume 2017, Article ID 5151890, 8 pages https://doi.org/10.1155/2017/5151890 Research Article Analysis of NEA-NSC PWR Uncontrolled Control Rod Withdrawal
More informationDistribution of induced charge
E&M Lecture 10 Topics: (1)Distribution of induced charge on conducting plate (2)Total surface induced charge on plate (3)Point charge near grounded conducting sphere (4)Point charge near floating conducting
More informationAdvanced Structural Analysis EGF Cylinders Under Pressure
Advanced Structural Analysis EGF316 4. Cylinders Under Pressure 4.1 Introduction When a cylinder is subjected to pressure, three mutually perpendicular principal stresses will be set up within the walls
More informationMeasurement of electric potential fields
Measurement of electric potential fields Matthew Krupcale, Oliver Ernst Department of Physics, Case Western Reserve University, Cleveland Ohio, 44106-7079 18 November 2012 Abstract In electrostatics, Laplace
More informationSensitivity and Uncertainty Analysis Methodologies for Fast Reactor Physics and Design at JAEA
Sensitivity and Uncertainty Analysis Methodologies for Fast Reactor Physics and Design at JAEA Kick off meeting of NEA Expert Group on Uncertainty Analysis for Criticality Safety Assessment IRSN, France
More information2 Navier-Stokes Equations
1 Integral analysis 1. Water enters a pipe bend horizontally with a uniform velocity, u 1 = 5 m/s. The pipe is bended at 90 so that the water leaves it vertically downwards. The input diameter d 1 = 0.1
More informationAdvanced Multi-Physics Modeling & Simulation Efforts for Fast Reactors
Advanced Multi-Physics Modeling & Simulation Efforts for Fast Reactors J.W. Thomas Argonne National Laboratory IAEA Technical Meeting on Priorities in Modeling & Simulation in Fast Reactor Systems April
More informationThe total luminosity of a disk with the viscous dissipation rate D(R) is
Chapter 10 Advanced Accretion Disks The total luminosity of a disk with the viscous dissipation rate D(R) is L disk = 2π D(R)RdR = 1 R 2 GM Ṁ. (10.1) R The disk luminosity is half of the total accretion
More informationIntroduction to Heat and Mass Transfer. Week 7
Introduction to Heat and Mass Transfer Week 7 Example Solution Technique Using either finite difference method or finite volume method, we end up with a set of simultaneous algebraic equations in terms
More informationChapter 5 MATHEMATICAL MODELING OF THE EVACATED SOLAR COLLECTOR. 5.1 Thermal Model of Solar Collector System
Chapter 5 MATHEMATICAL MODELING OF THE EVACATED SOLAR COLLECTOR This chapter deals with analytical method of finding out the collector outlet working fluid temperature. A dynamic model of the solar collector
More informationLesson 6: Diffusion Theory (cf. Transport), Applications
Lesson 6: Diffusion Theory (cf. Transport), Applications Transport Equation Diffusion Theory as Special Case Multi-zone Problems (Passive Media) Self-shielding Effects Diffusion Kernels Typical Values
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More informationVHTR Thermal Fluids: Issues and Phenomena
VHTR Thermal Fluids: Issues and Phenomena www.inl.gov Technical workshop at PHYSOR 2012: Advanced Reactor Concepts April 15, 2012 Knoxville, TN Gerhard Strydom Idaho National Laboratory (INL) Overview
More informationSteam Condensation Induced Water Hammer Phenomena, a theoretical study
Steam Condensation Induced Water Hammer Phenomena, a theoretical study Imre Ferenc Barna and György Ézsöl Hungarian Academy of Sciences, KFKI Atomic Energy Research Institute(AEKI) Thermohydraulic Laboratory
More information