• Nuclear Engineering and Technology
• /
• v.53 no.7
• /
• pp.2377-2386
• /
• 2021

To model the atmospheric dispersion of radionuclides released from nuclear accident is very important for nuclear emergency. But the uncertainty of model parameters, such as source term and meteorological data, may significantly affect the prediction accuracy. Data assimilation (DA) is usually used to improve the model prediction with the measurements. The paper proposed a parameter bias transformation method combined with Lagrangian puff model to perform DA. The method uses the transformation of coordinates to approximate the effect of parameters bias. The uncertainty of four model parameters is considered in the paper: release rate, wind speed, wind direction and plume height. And particle swarm optimization is used for searching the optimal parameters. Twin experiment and Kincaid experiment are used to evaluate the performance of the proposed method. The results show that the proposed method can effectively increase the reliability of model prediction and estimate the parameters. It has the advantage of clear concept and simple calculation. It will be useful for improving the result of atmospheric dispersion model at the early stage of nuclear emergency.

• Nuclear Engineering and Technology
• /
• v.52 no.9
• /
• pp.1990-1997
• /
• 2020

To support the safe operation of the Miniature Neutron Source Reactor (MNSR), a thermo-hydraulic transient model using the RELAP5/Mod3.2 code was simulated. The model was verified by comparing the results with the measured and the previously calculated data. The comparisons consisted of comparing the MNSR parameters under normal constant power operation and reactivity insertion transients. Reactivity Insertion Accident (RIA) for three different initial reactivity values of 3.6, 6.0, and 6.53 mk have been simulated. The calculated peaks of the reactor power, fuel, clad and coolant temperatures in hot channel were calculated in this model. The reactor power peaks were: 103 kW at 240 s, 174 kW at 160 s and 195 kW at 140 s, respectively. The fuel temperature reached its maximum value of 116 ℃ at 240 s, 124 ℃ at 160 s and 126 ℃ at 140 s respectively. These calculation results ensured the high inherently safety features of the MNSR under all phases of the RIAs.

• Nuclear Engineering and Technology
• /
• v.55 no.4
• /
• pp.1280-1286
• /
• 2023

MgF₂ as a moderator material has been extensively used in the beam shaping assembly (BSA) that plays an important role in the boron neutron capture therapy (BNCT). Regarded as important for applications, the thermal neutron scattering data of MgF₂ were calculated, based on the phonon expansion model. The structural properties of MgF₂ were researched by the VASP code based on the ab-initio methods. The PHONOPY code was employed to calculate the phonon density of states. Furthermore, the NJOY code was used to calculate the thermal neutron scattering data of MgF₂. The calculated inelastic cross sections plus absorption cross sections are in agreement with the available experimental data. The neutron transport in the BSA has been simulated by using a hybrid Monte-Carlo-Deterministic code NECP-MCX. The results indicated that compared with the calculation of the free gas model, the thermal neutron flux and epithermal neutron flux at the BSA exit port calculated by using the thermal neutron scattering data of MgF₂ were reduced by 27.7% and 8.2%, respectively.

• Nuclear Engineering and Technology
• /
• v.51 no.5
• /
• pp.1261-1271
• /
• 2019

Droplet size and distribution are important parameters determining venturi scrubber performance. In this paper, we proposed physical models for a maximum stable droplet size prediction and upper limit log-normal (ULLN) distribution parameters. For the proposed maximum stable droplet size prediction model, a Eulerian-Lagrangian framework and a Reitz-Diwakar breakup model are solved simultaneously using CFD calculations to reflect the effect of multistage breakup and droplet acceleration. Then, two ULLN distribution parameters are suggested through best fitting the previously published experimental data. Results show that the proposed approach provides better predictions of maximum stable droplet diameter and Sauter mean diameter compared to existing simple empirical correlations including Boll, Nukiyama and Tanasawa. For more practical purpose, we developed the simple, one dimensional (1-D) calculation of Sauter mean diameter.

• Nuclear Engineering and Technology
• /
• v.42 no.2
• /
• pp.211-218
• /
• 2010

An analysis model on debris transport in the containment floor of pressurized water reactors is developed in which the flow field is calculated by Eulerian conservation equations of mass and momentum and the debris particles are traced by Lagrange equations of motion using the pre-determined flow field data. For the flow field calculation, two-dimensional Shallow Water Equations derived from Navier Stokes equations are solved using the Finite Volume Method, and the Harten-Lax-van Leer scheme is used for accuracy to capture the dry-to-wet interface. For the debris tracing, a simplified two-dimensional Lagrangian particle tracking model including drag force is developed. Advanced schemes to find the positions of particles over the containment floor and to determine the position of particles reflected from the solid wall are implemented. The present model is applied to calculate the transport fraction to the Hold-up Volume Tank in Advanced Power Reactors 1400. By the present model, the debris transport fraction is predicted, and the effect of particle density and particle size on transport is investigated.

• Nuclear Engineering and Technology
• /
• v.54 no.3
• /
• pp.834-841
• /
• 2022

System thermal-hydraulic (STH) code is adopted for nuclear safety analysis. The critical flow model (CFM) is significant for the accuracy of STH simulation. To overcome the defects of current CFMs (low precision or long calculation time), a CFM based on a genetic neural network (GNN) has been developed in this work. To build a powerful model, besides the critical mass flux, the critical pressure and critical quality were also considered in this model, which was seldom considered before. Comparing with the traditional homogeneous equilibrium model (HEM) and the Moody model, the GNN model can predict the critical mass flux with a higher accuracy (approximately 80% of results are within the ±20% error limit); comparing with the Leung model and the Shannak model for critical pressure prediction, the GNN model achieved the best results (more than 80% prediction results within the ±20% error limit). For the critical quality, similar precision is achieved. The GNN-based CFM in this work is meaningful for the STH code CFM development.

• Nuclear Engineering and Technology
• /
• v.47 no.6
• /
• pp.662-668
• /
• 2015

A suspension of titanium nanoparticles (Ti NPs) in liquid sodium (Na) has been proposed as a method to mitigate the violent sodium-water reaction (SWR). The interparticle potential between Ti NPs in liquid Na may play a significant role in the agglomeration of NPs on the reaction surface and in the bulk liquid Na, since the potential contributes to a reduction in the long-term dispersion stability. For the effective control of the SWR with NPs, a physical understanding of the molecular dynamics of NPs in liquid Na is key. Therefore in this study, the nonretarded Van der Waals model and the solvation potential model are employed to analyze the interparticle potential. The ab initio calculation reveals that a strong repulsive force driven by the solvation potential exceeds the interparticle attraction and predicts the agglomeration energy required for two 10-nm Ti NPs to be $4{\times}10^{-17}J$. The collision theory suggests that Ti NPs can be effective suppressors of the SWR due to the high energy barrier that prevents significant agglomeration of Ti NPs in quiescent liquid Na.

• Nuclear Engineering and Technology
• /
• v.49 no.7
• /
• pp.1495-1504
• /
• 2017

In a previous study, a set of polygon-mesh (PM)-based skin models including a $50-{\mu}m-thick$ radiosensitive target layer were constructed and used to calculate skin dose coefficients (DCs) for idealized external beams of electrons. The results showed that the calculated skin DCs were significantly different from the International Commission on Radiological Protection (ICRP) Publication 116 skin DCs calculated using voxel-type ICRP reference phantoms that do not include the thin target layer. The difference was as large as 7,700 times for electron energies less than 1 MeV, which raises a significant issue that should be addressed subsequently. In the present study, therefore, as an extension of the initial, previous study, skin DCs for three other particles (photons, protons, and helium ions) were calculated by using the PM-based skin models and the calculated values were compared with the ICRP-116 skin DCs. The analysis of our results showed that for the photon exposures, the calculated values were generally in good agreement with the ICRP-116 values. For the charged particles, by contrast, there was a significant difference between the PM-model-calculated skin DCs and the ICRP-116 values. Specifically, the ICRP-116 skin DCs were smaller than those calculated by the PM models-which is to say that they were under-estimated-by up to ~16 times for both protons and helium ions. These differences in skin dose also significantly affected the calculation of the effective dose (E) values, which is reasonable, considering that the skin dose is the major factor determining effective dose calculation for charged particles. The results of the current study generally show that the ICRP-116 DCs for skin dose and effective dose are not reliable for charged particles.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3450-3463
• /
• 2023

In the two-step analysis of Pebble-Bed type High-Temperature Gas-Cooled Reactor (PB-HTGR), the lattice physics calculation for the generation of homogenized cross-sections is based on the fuel pebble. However, the randomly-dispersed fuel particles in the fuel pebble introduce double heterogeneity and randomness. Compared to the deterministic method, the Monte Carlo method which is flexible in geometry modeling provides a high-fidelity treatment. Therefore, the Monte Carlo code NECP-MCX is extended in this study to perform the lattice physics calculation of the PB-HTGR. Firstly, the capability for the simulation of randomly-dispersed media, using the explicit modeling approach, is developed in NECP-MCX. Secondly, the capability for the generation of the homogenized cross-section is also developed in NECP-MCX. Finally, simplified PB-HTGR problems are calculated by a two-step neutronics analysis tool based on Monte Carlo homogenization. For the pebble beds mixed by fuel pebble and graphite pebble, the bias is less than 100 pcm when compared to the high-fidelity model, and the bias is increased to 269 pcm for pebble bed mixed by depleted fuel pebble. Numerical results show that the Monte Carlo lattice physics calculation for the two-step analysis of PB-HTGR is feasible.

• Nuclear Engineering and Technology
• /
• v.52 no.4
• /
• pp.721-733
• /
• 2020

In nuclear engineering, the occurrence of critical heat flux (CHF) is complicated for rod bundle, and it is much more difficult to predict the CHF when it is in natural circulation under motion condition. In this paper, the dryout-type CHF is investigated for the rod bundle in a natural circulation loop under rolling motion condition based on the coupled analysis of subchannel method, a one-dimensional system analysis method and a CHF mechanism model, namely the three-fluid model for annular flow. In order to consider the rolling effect of the natural circulation loop, the subchannel model is connected to the one-dimensional system code at the inlet and outlet of the rod bundle. The subchannel analysis provides the local thermal hydraulic parameters as input for the CHF mechanism model to calculate the occurrence of CHF. The rolling motion is modeled by additional motion forces in the momentum equation. First, the calculation methods of the natural circulation and CHF are validated by a published natural circulation experiment data and a CHF empirical correlation, respectively. Then, the CHF of the rod bundle in a natural circulation loop under both the stationary and rolling motion condition is predicted and analyzed. According to the calculation results, CHF under stationary condition is smaller than that under rolling motion condition. Besides, the CHF decreases with the increase of the rolling period and angular acceleration amplitude within the range of inlet subcooling and mass flux adopted in the current research. This paper can provide useful information for the prediction of CHF in natural circulation under motion condition, which is important for the nuclear reactor design improvement and safety analysis.