Analyses of the V1000CT-1 benchmark with the DYN3D/ATHLET and DYN3D/RELAP coupled code systems including a coolant mixing model validated against CFD calculations

Plant-measured data provided within the specification of the OECD/NEA VVER-1000 coolant transient benchmark (V1000CT) were used to validate the DYN3D/RELAP5 and DYN3D/ATHLET coupled code systems. Phase 1 of the benchmark (V1000CT-1) refers to the MCP (main coolant pump) switching on experiment conducted in the frame of the plant-commissioning activities at the Kozloduy NPP Unit 6 in Bulgaria. The experiment was started at the beginning of cycle (BOC) with average core expose of 30.7 effective full power days (EFPD), when the reactor power was at 27.5% of the nominal level and three out of four MCPs were operating. The transient is characterized by a rapid increase in the primary coolant flow through the core and, as a consequence, a decrease of the space-dependent core inlet temperature. Both DYN3D/RELAP5 and DYN3D/ATHLET analyses were based on the same reactor model, including identical MCP characteristics, boundary conditions, benchmark-specified nuclear data library and nearly identical nodalization schemes. For an adequate modelling of the redistribution of the coolant flow in the reactor pressure vessel during the transient a simplified mixing model for the DYN3D/ATHLET code was developed and validated against a computational fluid dynamics calculation.

The results of both coupled code calculations are in good agreement with the available experimental data. The discrepancies between experimental data and the results of both coupled code calculations do not exceed the accuracy of the measurement data. This concerns the initial steady-state data as well as the time histories during the transient. In addition to the validation of the coupled code systems against measured data, a code-to-code comparison between simulation results has been performed to evaluate relevant thermal hydraulic models of the system codes RELAP5 and ATHLET and to explain differences between the calculation results.     

Validation of 3D neutronic-thermalhydraulic coupled codes RELAP5/PARCSv2.7 and TRACEv5.0P3/PARCSv3.0 against a PWR control rod drop transient

In nuclear safety field, neutronic and thermalhydraulic codes performance is an important issue. New capabilities implementation, as well as models and tools improvements are a significant part of the community effort in looking for better nuclear power plants (NPP) designs. A procedure to analyze the PWR response to local deviations on neutronic or thermalhydraulic parameters is being developed. This procedure includes the simulation of Incore and Excore neutron flux detectors signals. A control rod drop real plant transient is used to validate the used codes and their new capabilities. Cross-section data are obtained by means of the SIMTAB methodology. Detailed thermalhydraulic models were developed: RELAP5 and TRACE models simulate three different azimuthal zones. Besides, TRACE model is performed with a fully three-dimensional core, thus, the cross-flow can be obtained. A Cartesian vessel represents the fuel assemblies and a cylindrical vessel the bypass and downcomer. Simulated detectors signals are obtained and compared with the real data collected during a control rod drop trial at a PWR NPP and also with data obtained with SIMULATE-3K code.     

Analysis of the hot gas flow in the outlet plenum of the very high temperature reactor using coupled RELAP5-3D system code and a CFD code

The very high temperature reactor (VHTR) system behavior should be predicted during normal operating conditions and postulated accident conditions. The plant accident scenario and the passive safety behavior should be accurately predicted. Uncertainties in passive safety behavior could have large effects on the resulting system characteristics. Due to these performance issues in the VHTR, there is a need for development, testing and validation of design tools to demonstrate the feasibility of the design concepts and guide the improvement of the plant components. One of the identified design issues for the gas-cooled reactor is the thermal mixing of the coolant exiting the core into the outlet plenum. Incomplete thermal mixing may give rise to thermal stresses in the downstream components. To provide flow details, the analysis presented in this paper was performed by coupling a VHTR model generated in a thermal hydraulic systems code to a computational fluid dynamics (CFD) outlet plenum model. The outlet conditions obtained from the systems code VHTR model provide the inlet boundary conditions to the CFD outlet plenum model. By coupling the two codes in this manner, the important three-dimensional flow effects in the outlet plenum are well modeled while avoiding modeling the entire reactor with a computationally expensive CFD code. The values of pressure, mass flow rate and temperature across the coupled boundary showed differences of less than 5% in every location except for one channel. The coupling auxiliary program used in this analysis can be applied to many different cases requiring detailed three-dimensional modeling in a small portion of the domain.     

Effecting the condition of quasistationary thermonuclear burning in tokamaks with a high plasma density

Conclusion The zero-dimensional analysis and the numerical calculations have shown that stable conditions of the thermonuclear burning in a tokamak with a high plasma density are achieved only when _E depends strongly upon the plasma temperature (_E – ^m, where m<–2). Since in the experiments the retention laws have a significantly weaker temperature dependence, special measures are required for obtaining stable burning. Such conditions are reached by introducing feedbacks over, say, the plasma temperature and the rate of plasma-temperature change. It was shown in the present work that the corresponding regulation parameters are on the 0.3–0.5 level. Feedback can be obtained by, say, relaying operation of the source of additional heating.Translated from Atomnaya Énergiya, Vol. 63, No. 2, pp. 147–148, August, 1987.     

Fluid simulation of the effect of a dielectric window with high temperature on plasma parameters in inductively coupled plasma

To maintain the high-density plasma source in inductively coupled plasma (ICP), very high radiofrequency power is often delivered to the antenna, which can heat the dielectric windows near the antenna to high temperature. This high temperature can modulate the plasma characteristics to a large degree. We thus study the effect of dielectric window temperature on plasma parameters in two different ICP structures based on COMSOL software. The distributions of various plasma species are examined at different dielectric window temperatures. The concentration of neutral gas is found to be largely modulated at high dielectric window temperature, which further affects the electron collision probability with neutrals and the electron temperature. However, the electron density profiles are barely affected by the dielectric window temperature, which is mainly concentrated at the center of the reactor due to the fixed power input and pressure.     

Impact of cross-section generation procedures on the simulation of the VVER-1000 pump startup experiment in the OECD/DOE/CEA V1000CT benchmark by coupled 3D thermal-hydraulics/neutron kinetics models

In the framework of joint effort between the Nuclear Energy Agency (NEA) of OECD, the United States Department of Energy (US DOE), and the Commissariat a l'Enerige Atomique (CEA), France a coupled three-dimensional (3D) thermal-hydraulics/neutron kinetics benchmark was defined. The overall objective of OECD/NEA V1000CT benchmark is to assess computer codes used in analysis of VVER-1000 reactivity transients where mixing phenomena (mass flow and temperature) in the reactor pressure vessel are complex. Original data from the Kozloduy-6 Nuclear Power Plant are available for the validation of computer codes: one experiment of pump start-up (V1000CT-1) and one experiment of steam generator isolation (V1000CT-2). Additional scenarios are defined for code-to-code comparison. As a 3D core model is necessary for a best-estimate computation of all the scenarios of the V1000CT benchmark, all participants were asked to develop their own core coupled 3D thermal-hydraulics/neutron kinetics models using the data available in the benchmark specifications and a common cross-section library. The first code-to-code comparisons based on the V1000CT-1 Exercise 2 specifications exhibited unacceptable discrepancies between two sets of results. The present paper focuses on the analysis of the observed discrepancies. The VVER-1000 3D neutron kinetics models are based on cross-section data homogenized on the assembly level. The cross-section library, provided as part of the benchmark specifications, thus consists in a set of parameterized two group cross sections representing the different assemblies and the reflectors. The origin of the observed large discrepancies was found mainly to lie in the methods used to solve the diffusion equation. The VVER reflector properties were also found to enhance discrepancies by increasing flux gradients at the core/reflector interface thus highlighting more the difficulties in some codes to handle high exponential flux gradients. This paper summarizes the different steps applied to analyze the neutronic codes and their predictions as well as the impact of cross-section generation procedures.     

Validation of the 3D code SARUMAN in a stability and quench experiment with the superconducting magnet SAFFO

Experiments on stability and quench propagation were performed on a NbTi hollow conductor test module allowing heat conduction only between adjacent layers, while the turns in a layer are thermally insulated (2-D winding). The measurements performed over a wide range of operating conditions were used as the experimental data-base for the comparison with the results of the simulations done using the quench analysis code SARUMAN. In the paper the experimental results are compared with the predictions of a 2-D version of SARUMAN, which show that the simulation is able to reproduce the general features of the quench propagation experiments.