A FINITE VOLUME-BASED NETWORK METHOD FOR THE PREDICTION OF HEAT, MASS AND MOMENTUM TRANSFER IN A PEBBLE BED REACTOR
|
|
- Millicent Newton
- 6 years ago
- Views:
Transcription
1 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 Bag X6001, Potchefstroom 2520, South Africa Corresponding author HJ van Antwerpen: Fax: mgihjva@puk.ac.za 1. ABSTRACT In this paper a new control volume based network scheme for the prediction of heat, mass and momentum transfer in a pebble bed reactor is presented. The discretization scheme takes detail effects such as convection, radiation and conduction heat transfer into account and also provides for several types of reactor flow elements such as pipes, annuli, risers and porous media. The network representation of the reactor is solved as an integral part of the overall thermal-fluid network for the entire power conversion unit. The new control volume based network scheme facilitates the setting up of the reactor network taking detail effects such as complex radiation and effective pebble bed-solid interface conduction into account. The model is used as a benchmark for the current element-based network scheme and it is also used to investigate the effect of more accurate radiation discretization on the temperature profile in the reactor. Excellent agreement has been found with the current element-based reactor model. 2. INTRODUCTION To design a complex thermal-fluid system such as the Pebble Bed Modular Reactor (PBMR) various levels of analysis and simulation are required, from detailed component-level CFD simulations to integrated system-level simulation. In most system-level simulation codes simple one-dimensional models represent components such as pipes, pumps and heat transfer elements. However, nuclear reactor behaviour is too complex to be accurately modelled as a simple one-dimensional element in a system simulation code. Modelling the reactor is not a simple task because of the complexity of the phenomena involved. In the pebble bed, heat is generated in kernels inside the fuel spheres. The heat is then conducted to the sphere surface where it is transferred to the helium by means of convection. In the bed itself, heat is transferred by means of conduction and interstitial radiation in the radial and axial directions. In addition to this, the nuclear reaction that produces the heat is temperature dependent, making the variables in the problem very interrelated. One approach to solve this problem is to couple a systems simulation code with a threedimensional Computational Fluid Dynamics (CFD) code that models the reactor. Classical CFD calculations are indeed very powerful and most commercially available CFD codes are capable of solving the flow in complex three-dimensional geometries. However, these capabilities require large meshes, sophisticated discretization schemes and intricate data management that make CFD calculations computationally very expensive. An additional consideration is that in the case of a pebble bed, even CFD models assume a porous medium, which causes the flow equations to simplify considerably due to the dominance of the bed resistance. Given this simplification and the fact that the cost of the CFD computational overhead remains, the use of classical CFD offers very little benefit in this case. Even if a very simplified CFD model is coupled to a system simulation code, there are still two separate solvers running simultaneously on parts of the same problem, so that information has to be 1
2 exchanged either with an explicit or semi-implicit interface. This requires considerable programming effort to maintain consistency between the two codes. An explicit interface furthermore severely limits the timestep length in transient simulations [1]. Reactor models for system simulation have previously been developed using the thermal-fluid network approach, in which a network of elements was assembled. Elements are onedimensional entities that represent either flow paths or heat transfer paths. Heat transfer paths include conduction, convection and radiation heat transfer elements. Du Toit et al[3] used such an element-based approach in Flownex[2], the code used for the PBMR plant simulation. Flownex is a systems-cfd code in which the unsteady compressible flow and energy equations are solved for an unstructured network of one-dimensional elements representing components such as pipes, valves, turbomachinery, heat exchangers or heat transfer elements. As described for this context, the fundamental mass, momentum and energy equations governing unsteady compressible flow through the reactor is discretised in nodes representing mass and volume, interconnected by elements representing heat and mass flux between nodes. Phenomena that are taken into account include convection heat transfer between all gas and solid structures as well as solid conduction. The assembled network comprising the reactor sub-network and the rest of the power conversion system is solved as one system using the implicit pressure correction method [4]. In this way, a two-dimensional solution of the reactor is obtained as part of the solution of the overall power conversion system. This has the advantage that the solution is in no way constrained by the coupling interface between a systems-cfd and a classical CFD code, as a fully implicit coupling exists between the reactor and the rest of the system. In this way, a balance is found between accuracy and calculation speed, while taking into account all phenomena that influence the reactor interaction with the rest of the system. However, this element based discretization does not lend itself to easily take detail effects like radiation or surface temperature into account. A need therefore developed to devise an alternative discretization scheme that firstly will facilitate the modelling of certain detail effects and secondly, that can serve as a benchmark for element based discretization model. This paper describes a new control volume based network approach with which the discretization process for complex geometries is eased considerably. Complex effects like pebble bed effective conduction and radiation heat transfer in cavities are easier taken into account. 3. RELATION BETWEEN THE CURRENT AND BENCHMARK REACTOR MODELS The control volume based approach serves the following three purposes: To facilitate the discretization of complex geometries, taking detail effects into account. To serve as benchmark for the element-based reactor model because the discretization approach is totally different than that of the existing model. To investigate the effect of radiation in the reactor as detail radiation is taken into account with the control volume based discretization scheme. 2
3 With the control volume based approach full control volumes are used throughout the discretization, as opposed to half control volumes at the boundaries in case of the elementbased approach. Boundary or interface values are then calculated with massless nodes connected to full control volumes having mass, as shown in FIG 1 and FIG 8. The use of massless nodes is made possible by the implicit solution algorithm, as an explicit time discretization requires nonzero node masses. Zero volume nodes Full volume nodes FIG 1. Element-based model (left) with half control volumes at boundaries and some interfaces, and control volume based model (right) with full control volumes and massless nodes at boundaries. What the control volume based network approach has in common with classical CFD control volume approach is that momentum conservation is solved at control volume interfaces (elements) while mass and energy conservation are solved at nodes. In contrast to the control volumes used in classical CFD, the control volumes in this approach are of a much larger scale, incorporating empirical correlations or analytical models to take smaller-scale detail into account such as flow channels through a solid or conduction shape factors, without losing the simplicity of the large-scale control volumes. To construct the simulation network from a control volume specification chart, five fundamental elements are used, namely gas flow elements, conduction elements, convection heat transfer elements as well as two types of radiation elements: surface elements and spatial elements. The control volume type, as well as the types of its neighbours determine the network topology. The use of these elements at control volume-level is now illustrated with some control volume type examples. 4. SOLID WITH ONE-DIMENSIONAL FLOW A control volume type that illustrate the approach very well is solid with one-dimensional flow. Solid material, with flow passages passing through it, is typically encountered at the coolant riser channels in the side reflector or at control rod channels. Because of the pipegeometry of the flow channels, the Darcy-Weisbach correlation is used for calculating pressure drop, while the Dittus-Boelter correlation or a fixed heat transfer coefficient is used in the convection heat transfer between solid and fluid nodes [2]. For the solid conduction along the flow channel direction, area is reduced with a permeability-factor, while a conduction form factor can be used to account for the flow channel in transverse conduction. In this way, geometrical detail is taken into account without adding complexity to the solution. The network topology for this control volume type is as shown in FIG 2. 3
4 Solid with flow passage Solid node Gas node Solid conduction Convection heat transfer element Pipe flow element FIG 2. Nodes and elements used to represent a solid with one-dimensional flow through it. 5. PEBBLE BED REPRESENTATION In a pebble bed zone, heat is generated inside spheres by nuclear reaction, conducted to the sphere surface and transferred to the fluid by convection. When a temperature gradient exists across the bed, heat is transferred between pebbles by means of conduction through contact resistances as well as surface radiation exchange. These complex phenomena are modeled in a similar method to the current reactor model, of which the network topology for the pebble bed was established by Greyvenstein and Van Antwerpen [5]. This topology is based on the assumption of the pebble bed as a homogeneous conducting material of which the effective conductivity can be calculated with the Zehner-Schlünder correlation [2]. Maximum fuel temperature is calculated with a single representative sphere, having the heat generation, conduction resistances and internal node masses equivalent to the sum of all spheres in a control volume. Gas flow pressure drop is calculated with the well-known Ergun equation for a pebble bed, while the convection heat transfer coefficient is calculated with the correlation by Kugeler and Schulten [2]. The network representation of the abovementioned is presented in FIG 3. Packed bed with porous flow Solid node Gas node Representative fuel sphere with internal layers Point of maximum fuel temperature Porous flow elements Convection heat Pebble bed effective conduction and transfer element interstitial radiation heat transfer Fuel sphere surface node FIG 3. Assembly of nodes and elements representing a pebble bed zone. 6. CAVITIES AND RADIATIVE HEAT TRANSFER The surfaces of all solids bordering on cavity zones, exchange heat by means of radiation, which is of particular importance as radiative heat transfer increases with temperature to the fourth power and gas cooled reactors do operate at very high temperatures. Radiation heat exchange has previously not been fully calculated because of the complexity of calculating the 4
5 geometric view factor. In the element-based reactor model, only radiative heat transfer between directly opposing vertical surfaces is taken into account, so that the effect of radiation has not yet been accurately assessed. However, in the benchmark reactor model presented here, this effect has to be ascertained, even though it adds a large number of elements. For that purpose, analytical view factor formulae for discretised surfaces was derived from the formulae Brockmann [6] presented for right-angular cavities in a cylindrical geometry. As shown in FIG 4, the radiation network in a cavity is set up by connecting a radiative surface node to a massless surface node by means of a radiative surface element. The equation for its heat flux is then: A q = σε T T 1 ε 4 4 ( s sp ) (1) Each radiative surface node in a cavity is connected to every other radiative surface node by means of a radiative spatial element, for which the heat flux formula is given by: 4 4 ( ) q = A F σ T T (2) sp, sp1 sp2 Solid Cavity Primary nodes Massless nodes Radiative spatial elements Radiative surface elements Radiative surface nodes FIG 4. Radiation network across a cavity, showing surface and spatial radiation elements. 7. DISCRETIZATION OF REACTOR AND ASSEMBLY OF ZONES With the approach of having different control volumes, giving rise to appropriate network topologies, it is possible to automatically generate complete networks for system simulation. The process is illustrated in FIG 5 for a pebble bed type reactor with a centre column. 5
6 Pebble bed type reactor core assembly 2D Axi-symmetric Control volume type specification System CFD simulation network FIG 5. The discretization process for a pebble bed type reactor with a central column. 8. HEAT GENERATION AND NEUTRONICS Heat generation is managed by distributing a total specified power with normalised distribution curves obtained from detailed neutronics codes. Temperature and fast neutron flux history is not kept track of at present, but it is taken into account by specifying an appropriate temperature-dependent conductivity according to location in the reactor assembly. 9. STEADY-STATE RESULTS Using a pebble bed type reactor operating at full power, full mass flow conditions was simulated in steady state. The figures below present solid temperature profiles in two sections through the bed, axially and radially, with pebble bed surface temperature taken as solid temperature in the pebble bed. FIG 6 shows solid temperature profiles along an axial line at the middle of the pebble bed zone. In the pebble bed zone, sphere surface temperatures are used as solid temperatures. Maximum fuel temperatures differed by less than five percent, while the largest difference on the axial temperature profile is less than ten percent. Considering the difference in discretization methods, this can be considered excellent agreement. 6
7 Solid temperature, axial 1.2 Normalised temperature Normalised height Bencmark Current FIG 6. Axial solid temperature profile, in the middle of the pebble bed zone. 1.2 Solid temperature,radial Normalised temperature Normalised radius Benchmark Current FIG 7. Radial solid temperature profile, at the middle of the pebble bed zone. The radial temperature profile in FIG 7 has a large temperature percentage difference at a normalised radius of 0.7, just next to the pebble bed in the side reflector. Referring to FIG 8, this temperature difference can be explained as a direct result of discretization difference between the two topologies. The temperature of surface nodes in the pebble bed are largely determined by the heat transfer coefficient and the heat production rate in the pebbles, so that in the case of the element-based (current) model, the pebble bed boundary node effectively fixes the side reflector inner temperature at the pebble bed surface temperature. The control volume based (benchmark) model has no heat generation (or pebble bed surface) nodes in direct contact with the side reflector, but rather has the pebble bed surface temperature separated from the side reflector by an effective conductivity element. Because of the low effective conductivity of the pebble bed compared to the side reflector, there is a large thermal resistance between the pebble bed surface temperature and the side reflector. This highlights the impact of discretization at interfaces on the temperature profile in the reactor. 7
8 Pebble bed/side reflector interface Current model: element-based discretisation Benchmark model: volume-based discretisation Pebble bed Side reflector Pebble bed Side reflector Full mass node in reflector Pebble bed/reflector lumped mass Pebble bed Massless node surface node at interface Pebble bed effective conduction FIG 8. Different network topologies at the pebble bed/side reflector interface. 10. EFFECT OF RADIATION The current model only takes radiation in vertical cavity zones into account, which are typically used to model cavities between the core assembly and radiative heat shields and between radiative heat shields and the pressure vessel. The control-volume based benchmark model has a detailed radiation model. The only effect of taking radiation into account during steady-state full power operation is that the reactor assembly at the top and bottom has lower temperatures. 1.2 Solid temperature, axial Normalised temperature Normalised height Bench,Rad Current,Rad Bench,No Rad FIG 9. Axial solid temperature profile, in the middle of the pebble bed zone. One Benchmark case without radiative heat transfer in cavities is repeated as reference, the other cases with radiation. 11. CONCLUSION With the control-volume based network approach, the reactor can be modelled in greater detail, so that it can be used to validate other reactor models. Excellent agreement has been found with the current element-based reactor model. 8
9 12. ACKNOWLEGDEMENT The authors wish to acknowledge PBMR (Pty.) Ltd and the Technical Human Resources for Industry Programme (THRIP) of the National Research Foundation (NRF) who partially funded this research. 13. NOMENCLATURE A = area F 12 = radiation view factor Q = heat transfer T = temperature ε = emissivity σ = Stefan-Boltzmann constant 14. REFERENCES 1. Aumiller, D.L., Tomlinson, E.T., Bauer, R.C., A coupled RELAP5-3D/CFD methodology with a proof-of-principle calculation. Nuclear Engineering and Design 205, p83-90, Van der Merwe, J. and Van Ravenswaay, J.P., Flownex Version 6.4 User Manual, M-Tech Industrial, Potchefstroom, South Africa, Du Toit, C.G., Greyvenstein, G.P., Rousseau, P.G. (2003) A Comprehensive reactor model for the integrated network simulation of the PBMR power plant. Proceedings of the 2003 International Congress on Advances in Nuclear Power Plants (ICAPP 03). Cordoba, Spain, 4-7 May, Greyvenstein, G.P., An implicit method for the analysis of transient flow in pipe networks, Int. J. Numer. Meth. Engrng., 53, , Greyvenstein, G.P., Van Antwerpen, H.J., Rousseau P.G. The system CFD approach applied to a pebble bed reactor core. International Journal of Nuclear Energy Science and Technology. Inaugural Issue, Brockmann, H. Analytic angle factors for the radiant interchange among the surface elements of two concentric cylinders. Int. J. Heat Mass Transfer. Vol.37, No. 7. pp
Coolant Flow and Heat Transfer in PBMR Core With CFD
Heikki Suikkanen GEN4FIN 3.10.2008 1/ 27 Coolant Flow and Heat Transfer in PBMR Core With CFD Heikki Suikkanen Lappeenranta University of Technology Department of Energy and Environmental Technology GEN4FIN
More informationVALIDATION OF THE POINT KINETIC NEUTRONIC MODEL OF THE PBMR. D Marais BEng
Faculty Engineering VALIDATION OF THE POINT KINETIC NEUTRONIC MODEL OF THE PBMR D Marais BEng Mini-dissertation submitted in partial fulfilment of the requirements for the degree Master of Engineering
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 informationDepartment of Engineering and System Science, National Tsing Hua University,
3rd International Conference on Materials Engineering, Manufacturing Technology and Control (ICMEMTC 2016) The Establishment and Application of TRACE/CFD Model for Maanshan PWR Nuclear Power Plant Yu-Ting
More informationEFFECT OF DISTRIBUTION OF VOLUMETRIC HEAT GENERATION ON MODERATOR TEMPERATURE DISTRIBUTION
EFFECT OF DISTRIBUTION OF VOLUMETRIC HEAT GENERATION ON MODERATOR TEMPERATURE DISTRIBUTION A. K. Kansal, P. Suryanarayana, N. K. Maheshwari Reactor Engineering Division, Bhabha Atomic Research Centre,
More informationDocumentation of the Solutions to the SFPE Heat Transfer Verification Cases
Documentation of the Solutions to the SFPE Heat Transfer Verification Cases Prepared by a Task Group of the SFPE Standards Making Committee on Predicting the Thermal Performance of Fire Resistive Assemblies
More informationA Specialized Thermal-hydraulic Code with Porous Media Model and SIMPLEC Algorithm for PB-FHRs
A Specialized Thermal-hydraulic Code with Porous Media Model and SIMPLEC Algorithm for PB-FHRs Yanzhi E, Yang Zou, and Hongjie Xu Shanghai Institute of Applied Physics, Chinese Academy of Sciences Jialuo
More informationTHERMAL ANALYSIS OF A SPENT FUEL TRANSPORTATION CASK
Excerpt from the Proceedings of the COMSOL Conference 2009 Bangalore THERMAL ANALYSIS OF A SPENT FUEL TRANSPORTATION CASK P. Goyal*, Vishnu Verma, R.K. Singh & A.K. Ghosh Reactor Safety Division Bhabha
More informationSIMULATION OF FLOW IN A RADIAL FLOW FIXED BED REACTOR (RFBR)
SIMULATION OF FLOW IN A RADIAL FLOW FIXED BED REACTOR (RFBR) Aqeel A. KAREERI, Habib H. ZUGHBI, *, and Habib H. AL-ALI * Ras Tanura Refinery, SAUDI ARAMCO, Saudi Arabia * Department of Chemical Engineering,
More informationThis chapter focuses on the study of the numerical approximation of threedimensional
6 CHAPTER 6: NUMERICAL OPTIMISATION OF CONJUGATE HEAT TRANSFER IN COOLING CHANNELS WITH DIFFERENT CROSS-SECTIONAL SHAPES 3, 4 6.1. INTRODUCTION This chapter focuses on the study of the numerical approximation
More informationLectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 6
Lectures on Nuclear Power Safety Lecture No 6 Title: Introduction to Thermal-Hydraulic Analysis of Nuclear Reactor Cores Department of Energy Technology KTH Spring 2005 Slide No 1 Outline of the Lecture
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 informationTHE USE OF PB-BI EUTECTIC AS THE COOLANT OF AN ACCELERATOR DRIVEN SYSTEM. Joint research Centre of the European Commission Ispra, Italy.
THE USE OF PB-BI EUTECTIC AS THE COOLANT OF AN ACCELERATOR DRIVEN SYSTEM Alberto Peña 1, Fernando Legarda 1, Harmut Wider 2, Johan Karlsson 2 1 University of the Basque Country Nuclear Engineering and
More informationAnalysis of the Cooling Design in Electrical Transformer
Analysis of the Cooling Design in Electrical Transformer Joel de Almeida Mendes E-mail: joeldealmeidamendes@hotmail.com Abstract This work presents the application of a CFD code Fluent to simulate the
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 informationThis section develops numerically and analytically the geometric optimisation of
7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically
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 informationVALIDATION AND VERIFICATION OF FLOWNEX NUCLEAR
Proceedings HTR2006: 3 rd International Topical Meeting on High Temperature Reactor Technology October 1-4, 2006, Johannesburg, South Africa C00000169 VALIDATION AND VERIFICATION OF FLOWNEX NUCLEAR JJ
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 informationC ONTENTS CHAPTER TWO HEAT CONDUCTION EQUATION 61 CHAPTER ONE BASICS OF HEAT TRANSFER 1 CHAPTER THREE STEADY HEAT CONDUCTION 127
C ONTENTS Preface xviii Nomenclature xxvi CHAPTER ONE BASICS OF HEAT TRANSFER 1 1-1 Thermodynamics and Heat Transfer 2 Application Areas of Heat Transfer 3 Historical Background 3 1-2 Engineering Heat
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 informationStudies on flow through and around a porous permeable sphere: II. Heat Transfer
Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India
More informationOECD/NEA Transient Benchmark Analysis with PARCS - THERMIX
OECD/NEA Transient Benchmark Analysis with PARCS - THERMIX Volkan Seker Thomas J. Downar OECD/NEA PBMR Workshop Paris, France June 16, 2005 Introduction Motivation of the benchmark Code-to-code comparisons.
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 informationSimplified Model of WWER-440 Fuel Assembly for ThermoHydraulic Analysis
1 Portál pre odborné publikovanie ISSN 1338-0087 Simplified Model of WWER-440 Fuel Assembly for ThermoHydraulic Analysis Jakubec Jakub Elektrotechnika 13.02.2013 This work deals with thermo-hydraulic processes
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 informationCFD STUDIES IN THE PREDICTION OF THERMAL STRIPING IN AN LMFBR
CFD STUDIES IN THE PREDICTION OF THERMAL STRIPING IN AN LMFBR K. Velusamy, K. Natesan, P. Selvaraj, P. Chellapandi, S. C. Chetal, T. Sundararajan* and S. Suyambazhahan* Nuclear Engineering Group Indira
More informationModeling melting/solidification processes in the Molten Salt Fast Reactor (MEP)
Modeling melting/solidification processes in the Molten Salt Fast Reactor (MEP) The most innovative aspects of the Molten Salt Fast Reactor (MSFR), one of the six Generation IV nuclear reactors, are that
More informationChapter 2 HEAT CONDUCTION EQUATION
Heat and Mass Transfer: Fundamentals & Applications 5th Edition in SI Units Yunus A. Çengel, Afshin J. Ghajar McGraw-Hill, 2015 Chapter 2 HEAT CONDUCTION EQUATION Mehmet Kanoglu University of Gaziantep
More informationIntroduction to Heat Transfer
FIFTH EDITION Introduction to Heat Transfer FRANK P. INCROPERA College of Engineering University ofnotre Dame DAVID P. DEWITT School of Mechanical Purdue University Engineering THEODORE L. BERGMAN Department
More informationFAULT DIAGNOSIS OF GENERATION IV NUCLEAR HTGR COMPONENTS USING THE ENTHALPY-ENTROPY GRAPH APPROACH
Vol.3(3) September 22 SOUTH AFRICAN INSTITUTE OF ELECTRICAL ENGINEERS 27 FAULT DIAGNOSIS OF GENERATION IV NUCLEAR HTGR COMPONENTS USING THE ENTHALPY-ENTROPY GRAPH APPROACH C. P. du Rand* and G. van Schoor**
More informationUtilization of two-dimensional deterministic transport methods for analysis of pebble-bed reactors
annals of NUCLEAR ENERGY Annals of Nuclear Energy 34 (2007) 396 405 www.elsevier.com/locate/anucene Utilization of two-dimensional deterministic transport methods for analysis of pebble-bed reactors Bismark
More informationLecture 28. Key words: Heat transfer, conduction, convection, radiation, furnace, heat transfer coefficient
Lecture 28 Contents Heat transfer importance Conduction Convection Free Convection Forced convection Radiation Radiation coefficient Illustration on heat transfer coefficient 1 Illustration on heat transfer
More informationBiotransport: Principles
Robert J. Roselli Kenneth R. Diller Biotransport: Principles and Applications 4 i Springer Contents Part I Fundamentals of How People Learn (HPL) 1 Introduction to HPL Methodology 3 1.1 Introduction 3
More informationA 2-D Test Problem for CFD Modeling Heat Transfer in Spent Fuel Transfer Cask Neutron Shields. G Zigh and J Solis U.S. Nuclear Regulatory Commission
CFD4NRS2010 A 2-D Test Problem for CFD Modeling Heat Transfer in Spent Fuel Transfer Cask Neutron Shields G Zigh and J Solis U.S. Nuclear Regulatory Commission Abstract JA Fort Pacific Northwest National
More informationMATLAB Solution of Flow and Heat Transfer through a Porous Cooling Channel and the Conjugate Heat Transfer in the Surrounding Wall
MATLAB Solution of Flow and Heat Transfer through a Porous Cooling Channel and the Conjugate Heat Transfer in the Surrounding Wall James Cherry, Mehmet Sözen Grand Valley State University, cherryj1@gmail.com,
More informationAN ABSTRACT OF THE THESIS OF
AN ABSTRACT OF THE THESIS OF Jordan Cox for the degree of Master of Science in Nuclear Engineering presented on September 16, 2015. Title: Scaling of the Chinese HTR-PM Reactor Design for Licensing and
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 informationIntroduction to Heat and Mass Transfer. Week 9
Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional
More informationThe Dynamical Loading of the WWER440/V213 Reactor Pressure Vessel Internals during LOCA Accident in Hot and Cold Leg of the Primary Circuit
The Dynamical Loading of the WWER440/V213 Reactor Pressure Vessel Internals during LOCA Accident in Hot and Cold Leg of the Primary Circuit ABSTRACT Peter Hermansky, Marian Krajčovič VUJE, Inc. Okružná
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 informationTHERMAL HYDRAULIC ANALYSIS IN REACTOR VESSEL INTERNALS CONSIDERING IRRADIATION HEAT
THERMAL HYDRAULIC ANALYSIS IN REACTOR VESSEL INTERNALS CONSIDERING IRRADIATION HEAT Sungje Hong, Kunwoo Yi, Jin Huh, Inyoung Im and Eunkee Kim KEPCO Engineering and Construction Company. INC. NSSS Division.
More informationDEVELOPMENT OF A COUPLED CODE SYSTEM BASED ON SPACE SAFETY ANALYSIS CODE AND RAST-K THREE-DIMENSIONAL NEUTRONICS CODE
DEVELOPMENT OF A COUPLED CODE SYSTEM BASED ON SPACE SAFETY ANALYSIS CODE AND RAST-K THREE-DIMENSIONAL NEUTRONICS CODE Seyun Kim, Eunki Lee, Yo-Han Kim and Dong-Hyuk Lee Central Research Institute, Korea
More informationPWR CONTROL ROD EJECTION ANALYSIS WITH THE MOC CODE DECART
PWR CONTROL ROD EJECTION ANALYSIS WITH THE MOC CODE DECART Mathieu Hursin and Thomas Downar University of California Berkeley, USA mhursin@nuc.berkeley.edu,downar@nuc.berkeley.edu ABSTRACT During the past
More informationIf there is convective heat transfer from outer surface to fluid maintained at T W.
Heat Transfer 1. What are the different modes of heat transfer? Explain with examples. 2. State Fourier s Law of heat conduction? Write some of their applications. 3. State the effect of variation of temperature
More informationPebble flow in a High Temperature Reactor
Pebble flow in a High Temperature Reactor Andrea Ooms Physics of Nuclear Reactors TU Delft PNR-131-2008-003 May 15, 2008 Preface This thesis represents my graduation work for the Master in Applied Physics
More informationTutorial for the heated pipe with constant fluid properties in STAR-CCM+
Tutorial for the heated pipe with constant fluid properties in STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities
More informationA First Course on Kinetics and Reaction Engineering Unit D and 3-D Tubular Reactor Models
Unit 34. 2-D and 3-D Tubular Reactor Models Overview Unit 34 describes two- and three-dimensional models for tubular reactors. One limitation of the ideal PFR model is that the temperature and composition
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 informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationDevelopments and Applications of TRACE/CFD Model of. Maanshan PWR Pressure Vessel
Developments and Applications of TRACE/CFD Model of Maanshan PWR Pressure Vessel Yu-Ting Ku 1, Yung-Shin Tseng 1, Jung-Hua Yang 1 Shao-Wen Chen 2, Jong-Rong Wang 2,3, and Chunkuan Shin 2,3 1 : Department
More informationMONTE CARLO SIMULATION OF VHTR PARTICLE FUEL WITH CHORD LENGTH SAMPLING
Joint International Topical Meeting on Mathematics & Computation and Supercomputing in Nuclear Applications (M&C + SNA 2007) Monterey, California, April 5-9, 2007, on CD-ROM, American Nuclear Society,
More informationQUESTION ANSWER. . e. Fourier number:
QUESTION 1. (0 pts) The Lumped Capacitance Method (a) List and describe the implications of the two major assumptions of the lumped capacitance method. (6 pts) (b) Define the Biot number by equations and
More informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
More informationCalculating equation coefficients
Fluid flow Calculating equation coefficients Construction Conservation Equation Surface Conservation Equation Fluid Conservation Equation needs flow estimation needs radiation and convection estimation
More informationInstability Analysis in Peach Bottom NPP Using a Whole Core Thermalhydraulic-Neutronic Model with RELAP5/PARCS v2.7
Progress in NUCLEAR SCIENCE and TECHNOLOGY, Vol., pp.10-18 (011) ARTICLE Instability Analysis in Peach Bottom NPP Using a Whole Core Thermalhydraulic-Neutronic Model with RELAP/PARCS v. Agustín ABARCA,
More informationHEAT TRANSFER. PHI Learning PfcO too1. Principles and Applications BINAY K. DUTTA. Delhi Kolkata. West Bengal Pollution Control Board
HEAT TRANSFER Principles and Applications BINAY K. DUTTA West Bengal Pollution Control Board Kolkata PHI Learning PfcO too1 Delhi-110092 2014 Contents Preface Notations ix xiii 1. Introduction 1-8 1.1
More informationReview: Conduction. Breaking News
CH EN 3453 Heat Transfer Review: Conduction Breaking News No more homework (yay!) Final project reports due today by 8:00 PM Email PDF version to report@chen3453.com Review grading rubric on Project page
More informationAnalysis and interpretation of the LIVE-L6 experiment
Analysis and interpretation of the LIVE-L6 experiment A. Palagin, A. Miassoedov, X. Gaus-Liu (KIT), M. Buck (IKE), C.T. Tran, P. Kudinov (KTH), L. Carenini (IRSN), C. Koellein, W. Luther (GRS) V. Chudanov
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationNATURAL CONVECTION HEAT TRANSFER CHARACTERISTICS OF KUR FUEL ASSEMBLY DURING LOSS OF COOLANT ACCIDENT
NATURAL CONVECTION HEAT TRANSFER CHARACTERISTICS OF KUR FUEL ASSEMBLY DURING LOSS OF COOLANT ACCIDENT Ito D*, and Saito Y Research Reactor Institute Kyoto University 2-1010 Asashiro-nishi, Kumatori, Sennan,
More informationCFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel
CFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel *1 Hüseyin Kaya, 2 Kamil Arslan 1 Bartın University, Mechanical Engineering Department, Bartın, Turkey
More informationFLOW CHARACTERIZATION WITHIN A SPHERE-PACKED BED USING PIV MEASUREMENT
FLOW CHARACTERIZATION WITHIN A SPHERE-PACKED BED USING PIV MEASUREMENT J. H.ZHANG, L.XIONG, N.X.WANG and W ZHOU Department of reactor physics, Shanghai institute of applied physics, Chinese academy of
More informationChemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017
Chemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017 Objective: Text: To introduce the basic concepts of fluid mechanics and heat transfer necessary for solution of engineering
More informationNumerical Investigation of Convective Heat Transfer in Pin Fin Type Heat Sink used for Led Application by using CFD
GRD Journals- Global Research and Development Journal for Engineering Volume 1 Issue 8 July 2016 ISSN: 2455-5703 Numerical Investigation of Convective Heat Transfer in Pin Fin Type Heat Sink used for Led
More informationSafety Analyses for Dynamical Events (SADE) SAFIR2018 Interim Seminar Ville Sahlberg
VTT TECHNICAL RESEARCH CENTRE OF FINLAND LTD Safety Analyses for Dynamical Events (SADE) SAFIR2018 Interim Seminar Ville Sahlberg SADE Project SADE produces more reliable answers to safety requirements
More informationQUALIFICATION OF A CFD CODE FOR REACTOR APPLICATIONS
QUALIFICATION OF A CFD CODE FOR REACTOR APPLICATIONS Ulrich BIEDER whole TrioCFD Team DEN-STMF, CEA, UNIVERSITÉ PARIS-SACLAY www.cea.fr SÉMINAIRE ARISTOTE, NOVEMBER 8, 2016 PAGE 1 Outline Obective: analysis
More informationChapter 5. Formulation of FEM for Unsteady Problems
Chapter 5 Formulation of FEM for Unsteady Problems Two alternatives for formulating time dependent problems are called coupled space-time formulation and semi-discrete formulation. The first one treats
More informationTutorial for the supercritical pressure pipe with STAR-CCM+
Tutorial for the supercritical pressure pipe with STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities with
More informationUsing LBM to Investigate the Effects of Solid-Porous Block in Channel
International Journal of Modern Physics and Applications Vol. 1, No. 3, 015, pp. 45-51 http://www.aiscience.org/journal/ijmpa Using LBM to Investigate the Effects of Solid-Porous Bloc in Channel Neda Janzadeh,
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationNumerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions
Advanced Computational Methods in Heat Transfer X 25 Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions F. Selimovic & B. Sundén
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationThe thermal conductance of collection tubes in geothermal energy systems
Advanced Computational Methods and Experiments in Heat Transfer XIII 83 The thermal conductance of collection tubes in geothermal energy systems R. L. Frederick & A. Zenteno Departamento de Ingeniería
More informationThree-dimensional RAMA Fluence Methodology Benchmarking. TransWare Enterprises Inc., 5450 Thornwood Dr., Suite M, San Jose, CA
Three-dimensional RAMA Fluence Methodology Benchmarking Steven P. Baker * 1, Robert G. Carter 2, Kenneth E. Watkins 1, Dean B. Jones 1 1 TransWare Enterprises Inc., 5450 Thornwood Dr., Suite M, San Jose,
More informationEntropy 2011, 13, ; doi: /e OPEN ACCESS. Entropy Generation at Natural Convection in an Inclined Rectangular Cavity
Entropy 011, 13, 100-1033; doi:10.3390/e1305100 OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Article Entropy Generation at Natural Convection in an Inclined Rectangular Cavity Mounir
More informationLecture 30 Review of Fluid Flow and Heat Transfer
Objectives In this lecture you will learn the following We shall summarise the principles used in fluid mechanics and heat transfer. It is assumed that the student has already been exposed to courses in
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 informationA NUMERICAL APPROACH FOR ESTIMATING THE ENTROPY GENERATION IN FLAT HEAT PIPES
A NUMERICAL APPROACH FOR ESTIMATING THE ENTROPY GENERATION IN FLAT HEAT PIPES Dr. Mahesh Kumar. P Department of Mechanical Engineering Govt College of Engineering, Kannur Parassinikkadavu (P.O), Kannur,
More informationExact Solution of an MHD Natural Convection Flow in Vertical Concentric Annulus with Heat Absorption
International Journal of Fluid Mechanics & Thermal Sciences 217; 3(5): 52-61 http://www.sciencepublishinggroup.com/j/ijfmts doi: 1.11648/j.ijfmts.21735.12 ISSN: 2469-815 (Print); ISSN: 2469-8113 (Online)
More informationStudy of the End Flux Peaking for the CANDU Fuel Bundle Types by Transport Methods
International Conference Nuclear Energy for New Europe 2005 Bled, Slovenia, September 5-8, 2005 Study of the End Flux Peaking for the CANDU Fuel Bundle Types by Transport Methods ABSTRACT Victoria Balaceanu,
More informationMNRSIM: An Interactive Visual Model which Links Thermal Hydraulics, Neutron Production and other Phenomena. D. Gilbert, Wm. J. Garland, T.
MNRSIM: An Interactive Visual Model which Links Thermal Hydraulics, Neutron Production and other Phenomena D. Gilbert, Wm. J. Garland, T. Ha McMaster University, Hamilton Ontario, CANADA 6th International
More informationSafety Analysis of Loss of Flow Transients in a Typical Research Reactor by RELAP5/MOD3.3
International Conference Nuclear Energy for New Europe 23 Portorož, Slovenia, September 8-11, 23 http://www.drustvo-js.si/port23 Safety Analysis of Loss of Flow Transients in a Typical Research Reactor
More informationmidas NFX 2015 R1 Release Note
Total Solution for True Analysis-driven Design midas NFX 2015 R1 Release Note 1 midas NFX R E L E A S E N O T E 2 0 1 5 R 1 Accurate whenever you Need Be more efficient by focusing on the right details
More informationTdyn-CFD+HT - Validation Case 9
Two-dimensional heat conduction with heat generation Version 14.0.0 Compass Ingeniería y Sistemas http://www.compassis.com Tel.: +34 932 181 989 - Fax.: +34 933 969 746 - E: info@compassis.com - C/ Tuset
More informationMathematical Investigation and CFD Simulation of Monolith Reactors: Catalytic Combustion of Methane
Excerpt from the Proceedings of the COMSOL Conference 8 Hannover Mathematical Investigation and CFD Simulation of Monolith Reactors: Catalytic Combustion of Methane Maryam Ghadrdan *,, Hamid Mehdizadeh
More informationHEAT TRANSFER COEFFICIENT CHARACTERIZATION AT THE SOLAR COLLECTOR WALL-FLUID INTERFACE
SASEC15 Third Southern African Solar Energy Conference 11 13 May 15 Kruger National Park, South Africa HEAT TRANSFER COEFFICIENT CHARACTERIZATION AT THE SOLAR COLLECTOR WALL-FLUID INTERFACE Mébarki Ghazali*
More informationSteady-State and Transient Neutronic and Thermal-hydraulic Analysis of ETDR using the FAST code system
Steady-State and Transient Neutronic and Thermal-hydraulic Analysis of ETDR using the FAST code system Sandro Pelloni, Evaldas Bubelis and Paul Coddington Laboratory for Reactor Physics and Systems Behaviour,
More informationModeling and analysis of a depressurized loss of forced cooling event in a thorium fueled high temperature reactor
Department Radiation Science & Technology Section Nuclear Energy and Radiation Applications Faculty of Applied Sciences Delft University of Technology Mekelweg 15, 2629 JB Delft The Netherlands Modeling
More information1D-3D COUPLED SIMULATION OF THE FUEL INJECTION INSIDE A HIGH PERFORMANCE ENGINE FOR MOTORSPORT APPLICATION: SPRAY TARGETING AND INJECTION TIMING
1D-3D COUPLED SIMULATION OF THE FUEL INJECTION INSIDE A HIGH PERFORMANCE ENGINE FOR MOTORSPORT APPLICATION: SPRAY TARGETING AND INJECTION TIMING M. Fiocco, D. Borghesi- Mahindra Racing S.P.A. Outline Introduction
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 informationPebble Bed Heat Transfer Particle-to-Fluid Heat Convection
Pebble Bed Heat Transfer Particle-to-Fluid Heat Convection Raluca Scarlat Thermal Hydraulics Laboratory Department of Nuclear Engineering University of California, Berkeley Group Meeting 26 February 2009
More informationCOURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour. Basic Equations in fluid Dynamics
COURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour Basic Equations in fluid Dynamics Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET 1 Description of Fluid
More informationEVALUATION OF PWR AND BWR CALCULATIONAL BENCHMARKS FROM NUREG/CR-6115 USING THE TRANSFX NUCLEAR ANALYSIS SOFTWARE
ANS MC2015 - Joint International Conference on Mathematics and Computation (M&C), Supercomputing in Nuclear Applications (SNA) and the Monte Carlo (MC) Method Nashville, Tennessee April 19 23, 2015, on
More informationThermal Systems. What and How? Physical Mechanisms and Rate Equations Conservation of Energy Requirement Control Volume Surface Energy Balance
Introduction to Heat Transfer What and How? Physical Mechanisms and Rate Equations Conservation of Energy Requirement Control Volume Surface Energy Balance Thermal Resistance Thermal Capacitance Thermal
More informationDEPARTMENT OF NUCLEAR ENGINEERING GRADUATE SCHOOL OF NUCLEAR AND ALLIED SCIENCES UNIVERSITY OF GHANA, LEGON. B.Sc. (KNUST), 2001
TH ER M A L ANALYSIS OF A PEBBLE BED H IG H TEM PER A TU R E GAS COOLED R EA C TO R A Thesis submitted to the: DEPARTMENT OF NUCLEAR ENGINEERING GRADUATE SCHOOL OF NUCLEAR AND ALLIED SCIENCES UNIVERSITY
More informationThermal resistance model for CSP central receivers
Thermal resistance model for CSP central receivers O. A. J. de Meyer, F. Dinter, and S. Govender Citation: AIP Conference Proceedings 1734, 030010 (2016); doi: 10.1063/1.4949062 View online: http://dx.doi.org/10.1063/1.4949062
More informationThermal Field in a NMR Cryostat. Annunziata D Orazio Agostini Chiara Simone Fiacco
Thermal Field in a NMR Cryostat Annunziata D Orazio Agostini Chiara Simone Fiacco Overall Objective of Research Program The main objective of the present work was to study the thermal field inside the
More informationNonlinear Iterative Solution of the Neutron Transport Equation
Nonlinear Iterative Solution of the Neutron Transport Equation Emiliano Masiello Commissariat à l Energie Atomique de Saclay /DANS//SERMA/LTSD emiliano.masiello@cea.fr 1/37 Outline - motivations and framework
More informationChapter 10 Flow in Conduits
Chapter 10 Flow in Conduits 10.1 Classifying Flow Laminar Flow and Turbulent Flow Laminar flow Unpredictable Turbulent flow Near entrance: undeveloped developing flow In developing flow, the wall shear
More information