反应堆退役三维辐射场实时计算及可视化

本文基于DELMIA和VIRTOOLS平台开发的反应堆退役三维仿真原型系统,提出了仿真系统、数据库和计算内核既相互独立又集成统一的三维辐射场计算和可视化技术方案。利用点核积分算法建立了三维辐射场计算模型,得到了能量的对数与转换系数的多项式拟合公式,考虑了设备屏蔽和自吸收效应。采用VS语言和SQL server软件平台编制了三维辐射场计算程序,经验证,在关键点处的辐射水平计算值与测量值的比值小于10,并嵌入了仿真系统,实现了退役场景三维辐射场的实时计算和数据更新。提出了基于行走路径的人员受照剂量计算方法,并实现了可视化显示。     

Use of MCAM in creating 3D neutronics model for ITER building

The three dimensional (3D) neutronics reference model of International Thermonuclear Experimental Reactor (ITER) only defines the tokamak machine and extends to the bio-shield. In order to meet further 3D neutronics analysis needs, it is necessary to create a 3D reference model of the ITER building. Monte Carlo Automatic Modeling Program for Radiation Transport Simulation (MCAM) was developed as a computer aided design (CAD) based bi-directional interface program between general CAD systems and Monte Carlo radiation transport simulation codes. With the help of MCAM version 4.8, the 3D neutronics model of ITER building was created based on the engineering CAD model. The calculation of the neutron flux map in ITER building during operation showed the correctness and usability of the model. This model is the first detailed ITER building 3D neutronics model and it will be made available to all international organization collaborators as a reference model.     

Simulation of bubbly flow in vertical pipes by coupling Lagrangian and Eulerian models with 3D random walks models: Validation with experimental data using multi-sensor conductivity probes and Laser Doppler Anemometry

A set of two phase flow experiments for different conditions ranging from bubbly flow to cap/slug flow have been performed under isothermal concurrent upward air–water flow conditions in a vertical column of 3 m height. Special attention in these experiments was devoted to the transition from bubbly to cap/slug flow. The interfacial velocity of the bubbles and the void fraction distribution was obtained using 2 and 4 sensors conductivity probes.Numerical simulations of these experiments for bubbly flow conditions were performed by coupling a Lagrangian code with an Eulerian one. The first one tracks the 3D motion of the individual bubbles in cylindrical coordinates (r, ?, z) inside the fluid field under the action of the following forces: buoyancy, drag, lift, wall lubrication. Also we have incorporated a 3D stochastic differential equation model to account for the random motion of the individual bubbles in the turbulent velocity field of the carrier liquid. Also we have considered the deformations undergone by the bubbles when they touch the walls of the pipe and are compressed until they rebound.The velocity and turbulence fields of the liquid phase were computed by solving the time dependent conservation equations in its Reynolds Averaged Transport Equation form (RANS). The turbulent kinetic energy k, and the dissipation rate ? transport equations were simultaneously solved using the k, epsilon model in a (r, z) grid by the finite volume method and the SIMPLER algorithm. Both Lagrangian and Eulerian calculations were performed in parallel and an iterative self-consistent method was developed. The turbulence induced by the bubbles is an important issue considered in this paper, in order to obtain good predictions of the void fraction distribution and the interfacial velocity at different gas and liquid flow conditions.     

The study of aeroball system for measuring 3D neutron flux distribution in reactor core

Aeroball system is attractive in several aspects because it can easily transport the map of neutron flux distribution to be measured from incore to outside of a reactor vessel.However,before the aeroball system is put to practical use in the heating reactor.there are four topics that have to be further studied.They are the stability of the activated positions,enhancement of signal/noise(S/N)ratio,distributed control and data-acquisition system and on-lin nbeutron flux distribution reconstruction.Besides describing the rasons for them,this paper gives out the theory,concept and solution about the first two topics and it is helptul to give the possibility to enhance the reactor-power.     

MeshCNN-based BREP to CSG conversion algorithm for 3D CAD models and its application

In the field of neutronics analysis, it is imperative to develop computer-aided modeling technology for Monte Carlo codes to address the increasing complexity of reactor core components by converting 3D CAD model(boundary representation, BREP) to MC model(constructive solid geometry, CSG). Separation-based conversion from BREP to CSG is widely used in computer-aided modeling MC codes because of its high efficiency, reliability, and easy implementation. However, the current separation-based BREP-...     

Validation of coupled neutronic/thermal-hydraulic code RELAP5-3D for RBMK-1500 reactor analysis application

This paper summarizes RELAP5-3D code validation activities carried out at the Lithuanian Energy Institute, which was performed through the modeling of RBMK-1500 specific transients taking place at Ignalina NPP. A best estimate RELAP5-3D model of the INPP RBMK-1500 reactor has been developed and validated against real plant data, as well as with the calculation results obtained using the Russian STEPAN/KOBRA code. The obtained calculation results demonstrate reasonable agreement with Ignalina NPP measured data. Behaviors of the separate MCC thermal-hydraulic parameters, as well as physical processes are predicted reasonably well to the real processes, occurring in the primary circuit of RBMK-1500 reactor. The calculated reactivity and the total reactor core power behavior in time are also in reasonable agreement with the measured plant data, which demonstrates the correct modeling of the neutronic processes taking place in RBMK-1500 reactor core. The performed validation of RELAP5-3D model of Ignalina NPP RBMK-1500 reactor allowed to improve the model, which in the future would be used for the safety substantiation calculations of RBMK-1500 reactors. Future activities are discussed.     

Development of three-dimensional reactor analysis code system for accelerator-driven system,ADS3D and its application with subcriticality adjustment mechanism

In order to perform the parametric survey for an accelerator-driven system (ADS) core with the subcriticality adjustment mechanism, a new calculation code system, ADS3D, was developed on MARBLE which is a comprehensive and versatile framework for reactor analysis. The application of ADS3D was also demonstrated on the neutronics design of ADS operated by control rod (CR) movement. Through the neutronics calculation, it was shown that the maximum proton beam current was decreased from 20.5 to 11.6 mA due to the switch from beam-operated to CR-operated core.     

R&D activities based on fast reactor cycle technologies for transmutation of TRU and LLFP by JNC

Research and development(R&D) activities on partitioning and transmutation of trans-uranium nuclides (TRU) and LLFP and future R&D program in JNC were summarized. Feasibility design studies have been conducting to investigate the characteristics of a fast reactor core with TRU and LLFP transmutation. It was reconfirmed that the fast reactor has a strong potential for transmuting TRU and LLFP, effectively. R&D for establishing partitioning process of TRU apart from the high-level radioactive wastes have been carried out. By several counter-current runs of the TRUEX process using highly active raffinates, a process flow sheet capable of selective partitioning of actinides and fission products was newly developed. JNC has settled a new R&D program concerning partitioning and transmutation of long-lived radioactive waste based on recommendation of check & review for OMEGA program performed by the Ad Hoc Committee under the Atomic Energy Commission of Japan (AEC). The R&D program is composed of the design studies and development of element technologies (fabrication, irradiation) in the “Feasibility Studies” on commercialized fast reactor system and the basic studies with experiments (nuclear data, reactor physics, fuel property, etc.) to establish database and analytical tools for the TRU- and LLFP- containing fuel and core design.     

Visualization of models of nuclear reactors for investigating and monitoring reactor behavior in the regular regime

The means for two-and three-dimensional visualization which are used in the MCU program and provide convenient methods for obtaining a graphical representation of computational models are described. An example of the application of these visualization means for solving a practical problem where it is necessary to determine the optimal number of neutrons per generation required to obtain reliable results in Monte Carlo calculations of RBMK reactors is presented. __________ Translated from Atomnaya énergiya, Vol. 104, No. 1, pp. 33–37, January, 2008.     

Modeling of void, fast power and graphite temperature reactivity coefficients measurements for the validation of RELAP5-3D RBMK-1500 reactor model

This paper deals with the modeling of RBMK-1500 specific transients taking place at Ignalina NPP: measurements of void and fast power reactivity coefficients, as well as change of graphite cooling conditions transient. The simulation of these transients was performed using RELAP5-3D code model of RBMK-1500 reactor. At the Ignalina NPP void and fast power reactivity coefficients are measured on a regular basis and based on the obtained experimental results the actual values of these reactivity coefficients are determined. Graphite temperature reactivity coefficient at the plant is determined by changing graphite cooling conditions in the reactor cavity. This type of transient is unique and important from the point of view of model validation for the gap between fuel channel and the graphite bricks. The measurement results, obtained during this transient, enabled to determine the thermal conductivity coefficient for this gap and to validate the graphite temperature reactivity feedback model. The performed validation of RELAP5-3D model of Ignalina NPP RBMK-1500 reactor allowed to improve the model, which in the future would be used for the safety substantiation calculations of RBMK-1500 reactors.     

Evaluation of D(d,n)3 He reaction neutron source models for BNCT irradiation system design

A mathematical method was developed to calculate the yield,energy spectrum and angular distribution of neutrons from D(d,n)~3He(D-D)reaction in a thick deuterium-titanium target for incident deuterons in energies lower than 1.0MeV.The data of energy spectrum and angular distribution were applied to set up the neutron source model for the beam-shaping-assembly(BSA)design of Boron-Neutron-Capture-Therapy(BNCT)using MCNP-4C code. Three cases of D-D neutron source corresponding to incident deuteron energy of 1000,400 and 150 key were inves- tigated.The neutron beam characteristics were compared with the model of a 2.45 MeV mono-energetic and isotropic neutron source using an example BSA designed for BNCT irradiation.The results show significant differences in the neutron beam characteristics,particularly the fast neutron component and fast neutron dose in air,between the non-isotropic neutron source model and the 2.5 MeV mono-euergetic and isotropic neutron source model.     

Simulation of start-up measurements with the KIKO3D code and its application for determining the scram rod worth uncertainty in the VVER-440 reactor

The reasons for large discrepancies between the computed and measured values of the efficiency of control rods observed during start-up experiments on the Russian pressurized water type VVER reactors are discussed. The numerical simulation of the measurements including the prediction of the ex-core detector signals was used to resolve the discrepancies. The time and space dependent neutron flux in the core during these measurements have been calculated by the KIKO3D nodal kinetic code. For calculating the ionization chamber signals the Green function technique has been applied. The Green functions of ionization chambers have been evaluated via solving the neutron transport equation in the reflector regions with the MCNP Monte Carlo code. The detector signals have been calculated and compared with measured ones using the inverse point kinetics transformation. Large number of asymmetric rod drop measurements (with one rod stuck) and some differential rod worth measurements from the Zero Power Physics Tests were provided by the Paks NPP for validation. The experiments cover different fuels (without and with enrichment zoning) and loading patterns. The intermediate range ionization chambers have been used during the scram measurements. The newly developed method provides fairly sufficient match of measured and calculated results. The time behavior of the detector readings observed in the measurements are described by the code in a consistent manner.As a further application the uncertainty of scram rod worth of the KARATE-440 code system was determined by static calculations and subsequent simulation of rod drop with the KIKO3D code. The calculated results were compared to measurements carried out by the Paks NPP. The uncertainty of scram rod worth is established by statistical analysis.     

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.     

下庄(336)铀矿床地质特征及主矿体三维模型构建

在充分研究下庄(336)矿床的成矿地质特征、控矿构造条件及矿体赋存形态的基础上,利用Surpac软件,实现地下矿体三维模型的构建,清晰展示主矿体三维空间的分布特点,总结了该“交点”型矿床的成矿规律.该模型的建立对下庄(336)矿床深部找矿突破具有重要的指导意义和应用价值.     

SCALE5.1程序系统中蒙特卡罗方法模块与离散纵坐标模块在乏燃料运输容器屏蔽计算中的比较分析

以某型乏燃料运输容器为计算模型,分别利用SCALE5.1程序系统中的一维离散纵标法程序和三维蒙特卡罗方法程序对运输容器进行了屏蔽计算,计算结果表明,两种方法的总当量剂量率结果相对偏差在10%以内。最后对两个模块的应用特点及差异进行了比较分析,为其在乏燃料容器屏蔽计算中的应用提供参考。     

HFETR三维堆芯燃料管理程序的开发与应用

介绍了用于高通量工程试验堆（HFETR）的三维堆芯燃料管理程序（HFM）的原理及功能,并应用HFM对HFETR堆5个临界装置实验堆芯和前3炉堆芯进行了跟踪计算。结果表明,HFM所建立的栅元计算、堆芯计算模型正确,且计算值与实验值符合良好,该程序能快速、准确地于HFETR的堆芯燃料管理。