Table Of ContentMoshe R. Heller
(Editor)
ASIMUTH
Nuclear Simulation
Second European Nuclear Simulation Symposium
Schliersee, October 1990 - Proceedings
With 74 Figures
Springer-Verlag
Berlin Heidelberg NewY ork
London Paris Tokyo Hong Kong Barcelona
Editor:
Moshe R. Heller
ASIMUTH -Applied Simulation Technology GmbH
Planegger StraBe 47
8000 Miinchen 60
ISBN-13:978-3-642-84281-8 e-ISBN-13:978-3-642-84279-5
DOl: 10.1007/978-3-642-84279-5
This work is subject to copyrighLAll rights are reserved, whetherthe whole or part oft he material is concerned,
specifically the rights oftranslation, reprinting, re-use of illustrations, recitation, broadcasting,reproduction on
microfilms or in other ways, and storage in data banks. Duplication of this publication or parts thereof is only
permitted under the provisions of the German Copyright Law ofSeptember9, 1965, in its current version and a
copyright fee must always be paid. Violations fall under the prosecution act of the German Copyright Law.
© Springer-Verlag Berlin, Heidelberg 1990
Softcover reprint of the hardcover 1s t edition 1990
ASIMUTH is a registered trademark of Applied Simulation Technology GmbH.
The use of registered names, trademarks, etc. in this publication does not imply, even in the absence ofa specific
statement, that such names are exempt from the relevant protective laws and regulations and therefore free for
general use.
216113020-543210 -Printed on acid-free paper
Preface
Welcome to Bavaria - West Germany to the SECOND EUROPEAN NUCLEAR SIMULATION
SYMPOSIUM. That workshop-type meeting is a follow-up of the first Symposium
which took place also at Schliersee in October 1987.
The objective of gatheri ng together is to cover most of the aspects of
Nuclear Mathematical Modelling and Simulation in theory and practice
to promote the exchange of knowledge and experience between different
national and international research groups in that field taking into con
sideration that Nuclear Safety in "United Europe" turned to be a main
issue.
Furthermore, there is to strengthen the international contact between
deve lopers and users of mode 11 i ng and simu 1a ti on techniques consi deri ng
the "open" European Community borders of 1992.
This triennial International Conference jointly promoted by ASIMUTH -
Applied Simulation Technology and the Gesellschaft fUr Reaktorsicherheit
(GRS) created an interest to publish new projects including their new
results. A large number of contributed papers has been strictly examined
and selected by the editorial committee to guarantee a high international
technical standard.
The book contai ns the accepted papers whi ch wi 11 be presented at the
Symposium. The papers have been classified according to the following
topics:
VI
1. NUCLEAR POWER PLANT SIMULATION
2. SOFTWARE TOOLS
3. ARTIFICIAL INTELLIGENCE
4. COMPUTATIONAL FLUID DYNAMICS
5. NEUTRON KINETICS
6. PLANT ANALYSER
7. COMPUTER BASED EDUCATION
8. NUCL!:AR \JASTE
Authors from 10 countries coming from three continents will meet at the
Symposium. Nuclear Community, Computer Systems Industries, Universities
and Research and Development Institutes, so that a broad spectrum of
simulation know-how is covered: Theory and Applications, Hardware and
Software, Research and Development.
The editor is grateful to the authors for making possible the publication
of this book, and especially to Mr. von Hagen and Ms. Raufelder of
Springer-Verlag for the confidence in ASIMUTH's ability and the excellent
new simulation volume respectively. Without the contributions of all
these people, the conference could not have been materialized.
Munich, October 1990 Moshe R. Heller
ASIMUTH
Applied Simulation
Technology GmbH
Table of Contents
PREFACE
NUCLEAR POWER PLANT SIMULATION SESSION
Experience with Simulation of Nuclear Systems
on Parallel Processing Computer Systems
(U. Graf, W.F. Werner) .•....••.....•.•.••...••....••..••.•..•.••.•. 3
A Pump Model for Use in Real-Time Simulation
(A.C. Barroso, L.V. Loureiro) ...•••.•...••..•.....•..•.•.•.....•.•• 15
Hydraulic Network Modelling for Real-Time Power Plant
Simulation with Computer Aided Code Generation
(R. Bakker, E. GroBe-Dunker, P. Leishman) ..•..•....•••..••.....•.•. 25
SOFTWARE TOOLS SESSION
SIPA, a Training and Engineering Simulator and AGLAE,
a Flexible Workshop for Model Generation
(C. Mathey, A. Valembois) ..•.......••..••..........••...•...•...... 41
An Environment for Parallel Structuring
of Fortran Programs
(K. Sridharan, M. McShea et al.) .•.•.••••.•...••...•..•.••.••..•.•. 53
Advanced Modular Simulation Techniques for
Nuclear Power Plants
(D. Saphier) •..••..••..••....•....•....•....•....•..•..•.•..•••.... 75
Practical Experience with Software Tools to Assess
and Improve the Quality of EXisting Nuclear Analysis
and Safety Codes
(N.H. Marshall, E.S. Marwill, S.D. Matthews, B.J. Stacey) 92
Nuclear Core and Power Plant Simulation on
High Performance Parallel Computer Systems
(R. Muller, R. Boer, H. Finnemann) .•.••••...•.••..•••••.•••..•.•... 104
VIII
ARTIFICIAL INTELLIGENCE SESSION
Knowledge-Based Systems to Support Dynamic
Process Simulation
(A. Lehmann) •.•.•••••. •. .•.• •. ..•.•. .. .. .. •.. . .. .••. ..•. .. .. .. .. ... 119
COMPUTATIONAL FLUID DYNAMICS SESSION
Two Phase Flow Analysis Capability
of Advanced Computer Codes
(W. Kohler. M. Schindler) 133
Finite Element Weak Statement CFD Algorithms
for Fluid-Thermal System Analysis
(A.J. Baker) .•.....••...•.....•...•.............•...•.•.....•.•...• 142
NEUTRON KINETICS SESSION
Real-Time Advanced Nuclear Reactor Core Model
(J. Koclas. F. Friedman. C. Paquette. P. Vivier) .........•......... 159
The Validity of the Point Kinetics Model
During Reactor Start-Up
(F.A. Tarantino. R.P. Jacqmin. A.F. Henry) .•...•........•....•.•... 170
Axial Flux Difference on Seijing Nuclear
Power Plant Simulator
(Z. Fuyu. L. Jingyu) 182
PLANT ANALYSER SESSION
APROS Nuclear Plant Analyser
(E.K. Puska. M. Saarinen. K. Porkholm) .••..••....•.•......•........ 193
COMPUTER BASED EDUCATION SESSION
Pressurised Water Reactor Simulation
in the Training Environment
(A.G. Wills) ••...••........•••.............••••....•.•............• 207
A Real-Time Core Model for
Nuclear Power Plant Simulators
(D. Lee. J.-A. Gyllander) .......................................... 220
NUCLEAR WASTE SESSION
Ground Water Flow Analysis of Potential
Low Level Radioactive Waste Disposal Sites
Using Electrical Circuit Analogies
(G.H. Miley. K.J. Kuelske) ....•.•..•..•••.•.•...•.•............•... 229
Nuclear Power Plant
Simulation Session
Experience with Simulation of Nuclear Systems
on Parallel Processing Computer Systems
U. GRAF, W.F. WERNER
Gesellschaft fur Reaktorsicherheit (GRS)mbH
Forschungsgelande
D - 8046 Garching
INTRODUCTION
Computer systems with parallel-processing capabilities are increasingly
being used in technical and scientific applications. Very crudely, two
types of computer architecture may be distinguished:
- systems consisting of several (usually 2-8) high performance
computers with a shared memory. No physical exchange of data
between memory areas allocated to the individual computer is neces
sary
- systems consisting of many (usually >16) computers which perform
the calculations using their own distributed memories. Physical exchange
of data between the local memories is required in such systems.
The paper describes some experience gained by running a fluid-dynamics
module (FLUBOX [5]) on both types of machines. The results can be considered
typical of many other applications involving the solution elliptic, hyper-
bolic and parabolic partial differential equations.
1. BASIC EQUATIONS
Many multi-dimensional simulation problems are described in terms of a
system of partial differential equations
4
=
~~ + L(u) R(u). with (1)
( ) =
L u lm(G. ( u ) aaiu: " + Ki ( u ) aax2u2 • )
i=l 1 1 1
R(u) = (Rl(u) •...• Rs(u»T is free of dervatives
G. and K. are matrix valued functions.
1 1
Initial conditions u(x, TO) = u(x) are prescribed at t = to. Boundary
conditions u(x, t) = g(t), t & [to' t1] are prescribed for x & 6 V.
The spatial semidiscretisation of the system of partial differential
equations (e.g. by a method as described in [1]) leads to initial value
problems (IVP) for ordinary differential equations (ODE):
=
~ f (t.u(t»
dt
=
u (to) Uo (2)
o ~ t & [to. t1]; u & C1 [to.tl]; u(t). ~~ & RN
For the solution of equation (2) (which is nonlinear since f is nonlinear)
the following discrete and linearized formula is considered
(3)
