RECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS
|
|
- Asher Richards
- 6 years ago
- Views:
Transcription
1 RECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS Dean Wang April 30, Nuclear Reactor Physics
2 Outline 2 Introduction and Background Coupled T-H/Neutronics Safety Analysis Numerical schemes for T-H transport equations and neutron kinetics diffusion equations Recent Developments: Implementation and Assessments of High-Resolution Numerical Methods in TRACE High-resolution schemes Flux Limiters Implementation and assessments Computational Challenges in Reactor Analysis
3 History of TRACE/PARCS Development 3 TRACE is the latest best-estimate reactor system thermal-hydraulics analysis code developed by NRC for analyzing transient and steady-state neutronicthermal-hydraulic behavior in light water reactors. PARCS is a three-dimensional reactor core neutron kinetics code which solves the steady-state and time-dependent, multi-group neutron diffusion equations. I Joined ISL 7/7/03 V5.0-stable branch Branch code tree for assessment 6/7/06 V5.0 Release 7/31/07 V5.0p1-stable branch Peer Review 8/07 8/08 V5.0p1 Release 10/17/08 V /4/10 PARCS 3.0 V5.0p2 Release 6/2/10 I Joined ORNL 12/5/10 V5.0p3 Release 5/11/12 V5.0p3ho 11/2012 V5.0p4 with ho Release 8/2014 I Joined UML 9/14 Main Development trunk
4 Coupled T-H/Neutronics 4 Radial Mapping Time-Step: n TRACE PARCS Axial Mapping Time-Step: n+1 TRACE PARCS Thermal-hydraulics (TRACE) Two-fluid, two-phase 6 field equations Space: Finite volume, 1st order upwind differencing Time: SETS or Semi-Implicit. Semi- Implicit is subject to the Courant limit, while SETS is more diffusive Neutronics (PARCS) Two-group nodal diffusion kinetics equations Space: CMFD with local two-node correction Time: the theta method with exponential transformation Coupling between T-H/Neutronics (TRACE/ PARCS) Space: radial and axial Time: explicit
5 Cross Sections for Neutronics 5 D. Wang, et al., Cross Section Generation Guidelines for TRACE-PARCS, NUREG/CR-7164, June 2013.
6 Reactor Modeling and Simulation Physical System Physical Laws SS and Transients LOCA AOO/ATWS Stability Mathematical Model Variables Field Equations Constitutive Equations Initial and Boundary Conditions Numerical Solution Numerical Methods Implicit or explicit First-order or highorder Computational Algorithms = 1 + 2
7 Two-Phase Flow Transport Equations 7 Mass: Momentum: Energy:
8 1 st -Order Upwind Finite-Volume Method 8 j-1 j j+1 j-3/2 j-1/2 j +1/2 u Upwind: Local Truncation Error:
9 Standard 2 nd -Order Schemes 9 They don t work well for two-phase flow!
10 Numerical Diffusion 10 Most reactor system analysis codes use the 1 st -order upwind scheme such as TRACE, RELAP, COBRA, etc. While very robust, 1 st -order upwinding leads to excessive numerical diffusion (damping). 2 nd -order methods can effectively reduce numerical diffusion, but often produce spurious oscillations and slow down convergence It has to be carefully treated for reactor analysis, particularly for BWR stability analysis and boron tracking
11 Fix for Standard 2 nd -Order Schemes 11 Flux Limiters:
12 Why These Flux Limiters Work? TVD Region and Monotonicity Limiter Functions MUSCL1.5 OSPRE Van Albada ENO 2 nd -order TVD Region 1.5 Φ(r) φ(r)=1 TVD Region r Theorem: Any TVD numerical method is monotonicity preserving
13 Desired Properties of Flux Limiters 13 Φ(1) = 1. This is a necessary requirement for 2 nd - order accuracy on smooth solu=ons Φ(r)/r = Φ(1/r). This symmetric property ensures that a flux-limiter has the same actions on forward and backward gradients Φ(r) is located in the 2 nd - order TVD region Φ(r) MUSCL1.5 OSPRE Van Albada ENO Limiter Functions r 2 nd -order TVD Region 13
14 14 Linear Transport Equation
15 Why Are These Limiters Important? 15 Provide a solution to the important and long-lasting numerical issue in the 1d system analysis codes Effectively reduce numerical damping without incurring spurious oscillations. They have a promising application in BWR stability analysis. We used to develop a non-uniform nodalization to mitigate numerical diffusion. Higher numerical accuracy can be achieved with these flux limiters on a uniform coarse nodalization. We initially introduced 2 nd -order limiter schemes into TRACE (ENO and others) Relatively simple algorithm Easy implementation within TRACE Very promising results
16 Implementation of Flux-Limiters in TRACE 16 Central Difference with flux limiters for the convection terms in the mass and energy conservation equations. The momentum equation in TRACE is in nonconservative form, in which 2 nd -order central difference is used for the convection term. Time integration scheme is semi-implicit
17 Assessment of Flux Limiters Convergence Performance 17 CPU time taken to achieve the steady-state convergence. The time-step size on the top of each column is the maximum time-step size allowed to achieve the SS convergence. 600 VA (and MUSCL) has very comparable performance with the upwind method in terms of CPU time. CPU time (s) Van Albada Upwind Time-Step Size (s)
18 18 Applications of Flux Limiter Schemes Ringhals Cycle 14 Coupled TRACE/PARCS BWR plant model 8 stability test points of Ringhals Cycle 14 Upwind largely underpredicts decay ratio of core-wide and out-of-phase oscillations 2 nd -order flux limiter VA has improved the overall accuracy by 7% TRACE DR TRACE DR TRACE V5.0p3 - Upwind Core-Wide DR Out-Of-Phase DR Uncertainty ±14% Benchmark DR TRACE V5.0p3 -VA Core-Wide DR Out-Of-Phase DR Uncertainty ±7% Benchmark DR P09
19 Oskarshamn Feedwater Transient Event 19 Oskarshamn-2 experienced a stability event on February 25, A loss of feedwater preheaters and control system logic failure resulted in high feedwater flow and low temperature Then, automatic power and flow control moved the reactor into a low flow - high power regime Very challenging coupled T-H/Neutronics problem The 2 nd -order flux limiter shows improvement in the prediction of nonlinear power oscillations with DR > 1.
20 Summary and Recommendations 20 A number of flux limiters were selected and implemented in TRACE for assessments. In this study, assessments were mainly focused on evaluating their performance for BWR stability analysis. VA, OSPRE, and MUSCL seem to work best in terms of numerical stability, accuracy, and convergence. Good numerical accuracy on a coarse nodalization can be achieved with these 2 nd -order flux limiter schemes. It is worth noting that 2 nd -order accuracy in both time and space can be achieved when using the Lax-Wendroff scheme with these flux limiters. These high-resolution methods are now officially released with TRACE V5.0p4.
21 2 nd -order Accurate both in Space and 21 Time where D. Wang, et al., "Solving Two-Phase Flow Transport Equations Using the Lax-Wendroff Scheme," Trans. AM. Nucl. Soc., 2015
22 Challenges in Computational Reactor Safety Analysis 22 Coupling between T-H and Neutronics Speedup of calculations Beyond CMFD Beyond 1D T-H
23 Implicit Coupling between T-H and Neutronics? 23 Fully implicit coupling Tightly coupled T-H, fuel heat conduction, and Neutronics equations Nonlinear system: need local linearization Jacobian matrix, Challenges Large matrix Complex coupling terms Storage Larger time-step for speedup? higher order time integration is needed J. K. Watson, PhD thesis, 2010 Is it really worth going implicit coupling?
24 Speedup of Reactor Safety Analysis CPU Time (s) Speedup # of CHANs in TRACE Parallelize XSEC and BiCGSTAB Modules Using OpenMP Peach Bottom TT2 (2s transient) BiCGSTAB XSEC Numbers of Cores TRACE PARCS - CMFD PARCS - Nodal PARCS - T/H PARCS - XS PARCS - Pin Power Coupled T-H/Neutronics computation is time-consuming BWR stability analysis AOO/ATWS Station blackout Speedup Relax time-step using high-order implicit time integration? n Not very promising SSP n Physical time scales, e.g. BWR flow/ power oscillations n Advanced iterative solvers for large linear systems Locally adaptive time stepping Parallelization and GPU computing
25 GPU Computing 25 NVIDIA Quadro K4200 CUDA Cores: 1344 GPU Memory: 4GB Memory Bandwidth: 173GB/s
26 Beyond CMFD? 26 Now the CMFD is widely used for solving core neutron diffusion equations or as a preconditioner for neutron transport calculations. However, the CMFD may encounter numerical stability problems and fails to converge. Solution: high-order Galerkin spectral element methods may provide viable approaches for solving neutron kinetics problems with better accuracy and computational efficiency.
27 Beyond 1D T-H and Low-Order CFD 27 SMRs, SFR and FHRs largely rely on natural circulation to remove decay heat Thermal stratification and mixing requires 3D modeling and simulation capability However, large scales and local complexities of the reactor system pose a paramount challenge for low-order CFD analysis, In particular, long-term transient problems. Solution: highly efficient CFD solvers
28 28 DOE CASL Advanced M&S
29 29 DOE NEAMS
30 What We Are Working on Now and 30 Plan to Do 2 nd -order numerical methods in BOTH time and space for TRACE V5.0p4. Implementation of high-resolution schemes in COBRA-TF (CASL). Locally adaptive time stepping schemes for twophase flow simulation. Development of new acceleration schemes for neutron transport calculations.
31 References 31 D. Wang, Reduce Numerical Diffusion in TRACE Using the High-Resolution Numerical Method ENO, Trans. AM. Nucl. Soc.107, 2012 D. Wang, et al., Implementation and Assessment of High-Resolution Numerical Methods in TRACE, Nuclear Engineering and Design 263 (2013) D. Wang, et al., "Solving Two-Phase Flow Transport Equations Using the Lax-Wendroff Scheme," Trans. AM. Nucl. Soc., 2015 N.P. Waterson, H. Deconinck, Design Principles for Bounded Higher-Order Convection Schemes a Unified Approach, Journal of Computational Physics, 224 (2007), P.K. Sweby, High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws, SIAM J. Num. Anal., 21 (1984), pp C.W. SHU and S. OSHER, Efficient Implementation of Essentially Non-oscillatory Shock Capturing Scheme, II, J. Comput. Phys., 83: 32-78, R.J. LeVeque, Numerical Methods for Conservation Laws, Lectures in Mathematics, ETH Zurich, Birkhauser, R.J. LeVeque, Finite Volume Methods for Hyperbolic Problems, Cambridge Texts in Applied Mathematics, M.O. Deville, et al., High-Order Methods for Incompressible Fluid Flow, Cambridge University Press, C. Hirsch, Numerical Computation of Internal and External Flows, 2 nd Edition, Wiley, J.S. Hesthaven, et al., Spectral Methods for Time-Dependent Problems, Cambridge J.S. Hesthaven, T. Warburton, Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications, Springer 2008.
32 32 Thank You!
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods Xinghui Zhong 12/09/ 2009 Outline 1 Introduction Why SSP methods Idea History/main reference 2 Explicit SSP Runge-Kutta Methods
More informationA Fifth Order Flux Implicit WENO Method
A Fifth Order Flux Implicit WENO Method Sigal Gottlieb and Julia S. Mullen and Steven J. Ruuth April 3, 25 Keywords: implicit, weighted essentially non-oscillatory, time-discretizations. Abstract The weighted
More informationTRACE/SIMULATE-3K Analysis of the NEA/OECD Oskarshamn-2 Stability Benchmark
TRACE/SIMULATE-3K Analysis of the NEA/OECD Oskarshamn-2 Stability Benchmark ABSTRACT Abdelhamid Dokhane, Omar Zerkak, and Hakim Ferroukhi Paul Scherrer Institute CH-5232 Villigen, Switzerland Abdelhamid.dokhane@psi.ch;
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
More informationAn efficient numerical method for hydraulic transient computations M. Ciccotelli," S. Sello," P. Molmaro& " CISE Innovative Technology, Segrate, Italy
An efficient numerical method for hydraulic transient computations M. Ciccotelli," S. Sello," P. Molmaro& " CISE Innovative Technology, Segrate, Italy Abstract The aim of this paper is to present a new
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationROYAL INSTITUTE OF TECHNOLOGY SENSITIVITY AND UNCERTAINTY ANALYSIS OF BWR STABILITY
ROYAL INSTITUTE OF TECHNOLOGY SENSITIVITY AND UNCERTAINTY ANALYSIS OF BWR STABILITY IVAN GAJEV LICENTIATE THESIS STOCKHOLM, SWEDEN, 2010 Abstract Best Estimate codes are used for licensing, but with conservative
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 informationInvestigation of Limit Cycle Behavior in BWRs with Time-Domain Analysis. Aaron John Wysocki
Investigation of Limit Cycle Behavior in BWRs with Time-Domain Analysis by Aaron John Wysocki A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy
More informationComputational and Experimental Benchmarking for Transient Fuel Testing: Neutronics Tasks
Computational and Experimental Benchmarking for Transient Fuel Testing: Neutronics Tasks T. Downar W. Martin University of Michigan C. Lee Argonne National Laboratory November 19, 2015 Objective of Neutronics
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationHigh-Order Finite Difference Nodal Method for Neutron Diffusion Equation
Journal of NUCLEAR SCIENCE and TECHNOLOGY, 28[4], pp. 285~292 (April 1991) 285 High-Order Finite Difference Nodal Method for Neutron Diffusion Equation Kazuo AZEKURA and Kunitoshi KURIHARA Energy Research
More informationSensitivity Analyses of the Peach Bottom Turbine Trip 2 Experiment
International Conference Nuclear Energy for New Europe 2003 Portorož, Slovenia, September 8-11, 2003 http://www.drustvo-js.si/port2003 Sensitivity Analyses of the Peach Bottom Turbine Trip 2 Experiment
More informationComparison of (Some) Algorithms for Edge Gyrokinetics
Comparison of (Some) Algorithms for Edge Gyrokinetics Greg (G.W.) Hammett & Luc (J. L.) Peterson (PPPL) Gyrokinetic Turbulence Workshop, Wolfgang Pauli Institute, 15-19 Sep. 2008 w3.pppl.gov/~hammett Acknowledgments:
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 informationANALYSIS OF THE OECD PEACH BOTTOM TURBINE TRIP 2 TRANSIENT BENCHMARK WITH THE COUPLED NEUTRONIC AND THERMAL-HYDRAULICS CODE TRAC-M/PARCS
ANALYSIS OF THE OECD PEACH BOTTOM TURBINE TRIP 2 TRANSIENT BENCHMARK WITH THE COUPLED NEUTRONIC AND THERMAL-HYDRAULICS CODE TRAC-M/PARCS Deokjung Lee and Thomas J. Downar School of Nuclear Engineering
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 informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More informationTOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES
MATHEMATICS OF COMPUTATION Volume 67 Number 221 January 1998 Pages 73 85 S 0025-5718(98)00913-2 TOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES SIGAL GOTTLIEB AND CHI-WANG SHU Abstract. In this paper we
More informationVISCOUS FLUX LIMITERS
VISCOUS FLUX LIMITERS E. F. Toro Department of Aerospace Science College of Aeronautics Cranfield Institute of Technology Cranfield, Beds MK43 OAL England. Abstract We present Numerical Viscosity Functions,
More informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationNUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS JOHN A. TRANGENSTEIN Department of Mathematics, Duke University Durham, NC 27708-0320 Ш CAMBRIDGE ЩР UNIVERSITY PRESS Contents 1 Introduction
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 informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationRattlesnake, MAMMOTH and Research in Support of TREAT Kinetics Calculations
Rattlesnake, MAMMOTH and Research in Support of TREAT Kinetics Calculations www.inl.gov DOE NEUP-IRP Meeting University of Michigan May 24, 2016 TREAT s mission is to deliver transient energy deposition
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
More informationWhole Core Pin-by-Pin Coupled Neutronic-Thermal-hydraulic Steady state and Transient Calculations using COBAYA3 code
Whole Core Pin-by-Pin Coupled Neutronic-Thermal-hydraulic Steady state and Transient Calculations using COBAYA3 code J. Jiménez, J.J. Herrero, D. Cuervo and J.M. Aragonés Departamento de Ingeniería Nuclear
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 informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationCross Section Generation Guidelines for TRACE- PARCS
NUREG/CR-764 ORNL/TM-202/58 Cross Section Generation Guidelines for TRACE- PARCS Office of Nuclear Regulatory Research AVAILABILITY OF REFERENCE MATERIALS IN NRC PUBLICATIONS NRC Reference Material As
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 information3.4. Monotonicity of Advection Schemes
3.4. Monotonicity of Advection Schemes 3.4.1. Concept of Monotonicity When numerical schemes are used to advect a monotonic function, e.g., a monotonically decreasing function of x, the numerical solutions
More informationARTICLE IN PRESS Mathematical and Computer Modelling ( )
Mathematical and Computer Modelling Contents lists available at ScienceDirect Mathematical and Computer Modelling ournal homepage: wwwelseviercom/locate/mcm Total variation diminishing nonstandard finite
More informationEFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL
EFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL Vahid SHOKRI 1*,Kazem ESMAEILI 2 1,2 Department of Mechanical Engineering, Sari Branch, Islamic Azad
More informationModelling and numerical methods for the diffusion of impurities in a gas
INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi
More informationNon-linear Methods for Scalar Equations
Non-linear Methods for Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 3, 04 / 56 Abstract
More informationLocal Truncation Error Analysis of Sn Discretizations
Research Note, Dean Wang, 28206 Local Truncation Error Analysis of Sn Discretizations D Monoenergetic Sn Equation The D monoenergetic Sn equation assuming isotropic scattering source is written as #$ μ
More informationTitle: Development of a multi-physics, multi-scale coupled simulation system for LWR safety analysis
Title: Development of a multi-physics, multi-scale coupled simulation system for LWR safety analysis Author: Yann Périn Organisation: GRS Introduction In a nuclear reactor core, different fields of physics
More informationSMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES
Proceedings of ALGORITMY 2016 pp. 155 164 SMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES TAMER H. M. A. KASEM AND FRANÇOIS G. SCHMITT Abstract. The weighted essentially non-oscillatory
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report March2008
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 08-08 March2008 The complexe flux scheme for spherically symmetrie conservation laws by J.H.M. ten Thije Boonkkamp,
More informationSome notes about PDEs. -Bill Green Nov. 2015
Some notes about PDEs -Bill Green Nov. 2015 Partial differential equations (PDEs) are all BVPs, with the same issues about specifying boundary conditions etc. Because they are multi-dimensional, they can
More informationPractical Combustion Kinetics with CUDA
Funded by: U.S. Department of Energy Vehicle Technologies Program Program Manager: Gurpreet Singh & Leo Breton Practical Combustion Kinetics with CUDA GPU Technology Conference March 20, 2015 Russell Whitesides
More informationComputational Engineering
Coordinating unit: 205 - ESEIAAT - Terrassa School of Industrial, Aerospace and Audiovisual Engineering Teaching unit: 220 - ETSEIAT - Terrassa School of Industrial and Aeronautical Engineering Academic
More informationA Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws
A Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws Mehdi Dehghan, Rooholah Jazlanian Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University
More informationA method for avoiding the acoustic time step restriction in compressible flow
A method for avoiding the acoustic time step restriction in compressible flow Nipun Kwatra Jonathan Su Jón T. Grétarsson Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA 9435 Abstract
More informationPALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations
2013 SIAM Conference On Computational Science and Engineering Boston, 27 th February 2013 PALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations U. Villa,
More informationLecture 18 Neutron Kinetics Equations
24.505 Lecture 18 Neutron Kinetics Equations Prof. Dean Wang For a nuclear reactor to operate at a constant power level, the rate of neutron production via fission reactions should be exactly balanced
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 informationSecond Order Positive Schemes by means of Flux Limiters for the Advection Equation
Second Order Positive Schemes by means of Flux Limiters for the Advection Equation Riccardo Fazio and Alessandra Jannelli Abstract In this paper, we study first and second order positive numerical methods
More informationCOMPARATIVE ANALYSIS OF WWER-440 REACTOR CORE WITH PARCS/HELIOS AND PARCS/SERPENT CODES
COMPARATIVE ANALYSIS OF WWER-440 REACTOR CORE WITH PARCS/HELIOS AND PARCS/SERPENT CODES S. Bznuni, A. Amirjanyan, N. Baghdasaryan Nuclear and Radiation Safety Center Yerevan, Armenia Email: s.bznuni@nrsc.am
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 informationOn High Order Strong Stability Preserving Runge Kutta and Multi Step Time Discretizations
Journal of Scientific Computing, Vol. 5, Nos. /, November 005 ( 005) DOI: 0.007/s095-00-65-5 On High Order Strong Stability Preserving Runge Kutta and Multi Step Time Discretizations Sigal Gottlieb Received
More informationThe Use of Serpent 2 in Support of Modeling of the Transient Test Reactor at Idaho National Laboratory
The Use of Serpent 2 in Support of Modeling of the Transient Test Reactor at Idaho National Laboratory Sixth International Serpent User s Group Meeting Politecnico di Milano, Milan, Italy 26-29 September,
More informationStrong stability preserving high order time discretizations.
SIGAL GOTTLIEB Mathematics Department 285 Old Westport Road North Dartmouth, MA 02747 sgottlieb@umassd.edu 322 Cole Avenue Providence, RI 02906 Phone: (401) 751-9416 sigalgottlieb@yahoo.com Current Research
More informationThermal Analysis Contents - 1
Thermal Analysis Contents - 1 TABLE OF CONTENTS 1 THERMAL ANALYSIS 1.1 Introduction... 1-1 1.2 Mathematical Model Description... 1-3 1.2.1 Conventions and Definitions... 1-3 1.2.2 Conduction... 1-4 1.2.2.1
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
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 informationLarge Time Step Scheme Behaviour with Different Entropy Fix
Proceedings of the Paistan Academy of Sciences: A. Physical and Computational Sciences 53 (): 13 4 (016) Copyright Paistan Academy of Sciences ISSN: 0377-969 (print), 306-1448 (online) Paistan Academy
More informationPOSITIVITY PROPERTY OF SECOND-ORDER FLUX-SPLITTING SCHEMES FOR THE COMPRESSIBLE EULER EQUATIONS. Cheng Wang. Jian-Guo Liu
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 3 Number May003 pp. 0 8 POSITIVITY PROPERTY OF SECOND-ORDER FLUX-SPLITTING SCHEMES FOR THE COMPRESSIBLE EULER EQUATIONS
More information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
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 informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationFlood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum
Archives of Hydro-Engineering and Environmental Mechanics Vol. 56 (29), No. 3 4, pp. 121 137 IBW PAN, ISSN 1231 3726 Flood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum Dariusz
More informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationClick to edit Master title style
Automated calculation sequence for group constant generation in Serpent 4th International Serpent UGM, Cambridge, UK, Sept. 17-19, 014 Jaakko Leppänen VTT Technical Research Center of Finland Click to
More informationGODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS
GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS ARTURO S. LEON Dept. of Civil and Envir. Engng., Univ. of Illinois at Urbana-Champaign, 2519 Hydrosystems Lab., MC-250. 205 North Mathews Av., Urbana,
More informationA New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods
A New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods Raymond J. Spiteri Steven J. Ruuth Technical Report CS-- May 6, Faculty of Computer Science 65 University Ave.,
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationAdaptive TVD-RK Discontinuous Galerkin Algorithms for Shallow Water Equations
Adaptive TVD-RK Discontinuous Galerkin Algorithms for Shallow Water Equations Thida Pongsanguansin, Khamron Mekchay and Montri Maleewong Abstract The adaptive Discontinuous Galerkin DG method for solving
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 informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationEntropic Schemes for Conservation Laws
CONSTRUCTVE FUNCTON THEORY, Varna 2002 (B. Bojanov, Ed.), DARBA, Sofia, 2002, pp. 1-6. Entropic Schemes for Conservation Laws Bojan Popov A new class of Godunov-type numerical methods (called here entropic)
More informationCross Section Generation Strategy for High Conversion Light Water Reactors Bryan Herman and Eugene Shwageraus
Cross Section Generation Strategy for High Conversion Light Water Reactors Bryan Herman and Eugene Shwageraus 1 Department of Nuclear Science and Engineering Massachusetts Institute of Technology 77 Massachusetts
More informationApplication of System Codes to Void Fraction Prediction in Heated Vertical Subchannels
Application of System Codes to Void Fraction Prediction in Heated Vertical Subchannels Taewan Kim Incheon National University, 119 Academy-ro, Yeonsu-gu, Incheon 22012, Republic of Korea. Orcid: 0000-0001-9449-7502
More informationStrong Stability Preserving Properties of Runge Kutta Time Discretization Methods for Linear Constant Coefficient Operators
Journal of Scientific Computing, Vol. 8, No., February 3 ( 3) Strong Stability Preserving Properties of Runge Kutta Time Discretization Methods for Linear Constant Coefficient Operators Sigal Gottlieb
More informationInternational Engineering Research Journal
Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationProcess Modelling. Table of Contents
Process Modelling 1 Process Modelling prepared by Wm. J. Garland, Professor, Department of Engineering Physics, McMaster University, Hamilton, Ontario, Canada More about this document Summary: The general
More informationFaster Kinetics: Accelerate Your Finite-Rate Combustion Simulation with GPUs
Faster Kinetics: Accelerate Your Finite-Rate Combustion Simulation with GPUs Christopher P. Stone, Ph.D. Computational Science and Engineering, LLC Kyle Niemeyer, Ph.D. Oregon State University 2 Outline
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 informationANALYSIS OF THE OECD MAIN STEAM LINE BREAK BENCHMARK PROBLEM USING THE REFINED CORE THERMAL-HYDRAULIC NODALIZATION FEATURE OF THE MARS0MASTER CODE
ANALYSIS OF THE OECD MAIN STEAM LINE BREAK BENCHMARK PROBLEM USING THE REFINED CORE THERMAL-HYDRAULIC NODALIZATION FEATURE OF THE MARS0MASTER CODE THERMAL HYDRAULICS KEYWORDS: MARS/MASTER code, coupled
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 informationREACTOR PHYSICS FOR NON-NUCLEAR ENGINEERS
International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering (M&C 2011) Rio de Janeiro, RJ, Brazil, May 8-12, 2011, on CD-ROM, Latin American Section (LAS)
More informationChapter 1. Introduction
Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear
More informationA Hybrid Deterministic / Stochastic Calculation Model for Transient Analysis
A Hybrid Deterministic / Stochastic Calculation Model for Transient Analysis A. Aures 1,2, A. Pautz 2, K. Velkov 1, W. Zwermann 1 1 Gesellschaft für Anlagen- und Reaktorsicherheit (GRS) ggmbh Boltzmannstraße
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationApplication of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms
Future Generation Computer Systems () 65 7 Application of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms R. Naidoo a,b, S. Baboolal b, a Department
More informationA Newton-Krylov Solution to the Coupled Neutronics-Porous Medium Equations
A Newton-Krylov Solution to the Coupled Neutronics-Porous Medium Equations by Andrew M. Ward A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy Nuclear
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationAnti-diffusive finite difference WENO methods for shallow water with. transport of pollutant
Anti-diffusive finite difference WENO methods for shallow water with transport of pollutant Zhengfu Xu 1 and Chi-Wang Shu 2 Dedicated to Professor Qun Lin on the occasion of his 70th birthday Abstract
More information