• Nuclear Engineering and Technology
• /
• v.54 no.8
• /
• pp.3059-3072
• /
• 2022

A Jacobian-Free Newton Krylov Two-Nodal Coarse Mesh Finite Difference algorithm based on Nodal Expansion Method (NEM_TNCMFD_JFNK) is successfully developed and proposed to solve the three-dimensional (3D) and multi-group reactor physics models. In the NEM_TNCMFD_JFNK method, the efficient JFNK method with the Modified Incomplete LU (MILU) preconditioner is integrated and applied into the discrete systems of the NEM-based two-node CMFD method by constructing the residual functions of only the nodal average fluxes and the eigenvalue. All the nonlinear corrective nodal coupling coefficients are updated on the basis of two-nodal NEM formulation including the discontinuity factor in every few newton steps. All the expansion coefficients and interface currents of the two-node NEM need not be chosen as the solution variables to evaluate the residual functions of the NEM_TNCMFD_JFNK method, therefore, the NEM_TNCMFD_JFNK method can greatly reduce the number of solution variables and the computational cost compared with the JFNK based on the conventional NEM. Finally the NEM_TNCMFD_JFNK code is developed and then analyzed by simulating the representative PWR MOX/UO₂ core benchmark, the popular NEACRP 3D core benchmark and the complicated full-core pin-by-pin homogenous core model. Numerical solutions show that the proposed NEM_TNCMFD_JFNK method with the MILU preconditioner has the good numerical accuracy and can obtain higher computational efficiency than the NEM-based two-node CMFD algorithm with the power method in the outer iteration and the Krylov method using the MILU preconditioner in the inner iteration, which indicates the NEM_TNCMFD_JFNK method can serve as a potential and efficient numerical tool for reactor neutron diffusion analysis module in the JFNK-based multiphysics coupling application.

• Advances in nano research
• /
• v.13 no.3
• /
• pp.207-216
• /
• 2022

In this paper, a novel model is developed for frequency behavior of single walled carbon nanotubes. The governing equation of motion is constructed method based on the Sander theory using Rayleigh-Ritz's method The frequencies enhances on increasing the power law index using simply supported, clamped and clamped free end conditions. The frequency curve for C-F is less than other conditions. It is due to the physical constraints which are applied on the edge of the CNT. It is observed that the C-F boundary condition have less frequencies from the other two conditions. The frequency phenomena for zigzag are insignificant throughout the aspect ratio. Moreover when index of power law is increased then frequencies increases for all boundary conditions. The natural frequency mechanism for the armchair (10, 10) for various values of power law index with different boundary conditions is investigated. Here frequencies decrease on increases the aspect ratio for all boundary conditions. The frequency curves of SS-SS edge condition is composed between the C-C and C-F conditions. The curves of frequency are less significant from small aspect ratio (L/d = 4.86 ~ 8.47) and decreases fast for greater ratios. It is found that the frequencies via aspect ratios, armchair (10, 10) have higher values from zigzag (10, 0). It is due to the material structure which is made by the carbon nanotubes. The power law index have momentous effect on the vibration of single walled carbon nanotubes. The present frequency result is also compared numerically experimentally with Raman Spectroscopy.

• Journal of the Korea Computer Graphics Society
• /
• v.27 no.3
• /
• pp.13-23
• /
• 2021

In this paper, we proposed a framework that generates the trajectory of a single rigid body based on its COM configuration and contact pose. Because we use a smaller input dimension than when we use a full body state, we can improve the learning time for reinforcement learning. Even with a 68% reduction in learning time (approximately two hours), the character trained by our network is more robust to external perturbations tolerating an external force of 1500 N which is about 7.5 times larger than the maximum magnitude from a previous approach. For this framework, we use centroidal dynamics to calculate the next configuration of the COM, and use reinforcement learning for obtaining a policy that gives us parameters for controlling the contact positions and forces.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.59 no.3
• /
• pp.189-195
• /
• 2023

This study quantitatively evaluated size selectivity for three netting shapes (T0; regular, T45, T90) and hanging ratio (35%, 70%) of T0 netting used for trawl codend. The size selectivity experiment was performed in a tank using a cube experimental model with a length of 50 cm on one side and 389 experimental individuals, jack mackerel (Trachurus japonicus). In the selectivity analysis, a selectivity curve was created based on the selection ratio using a logistic function, and the 25%, 50%, and 70% selection length and selection range (SR) were obtained. The T0 netting was 19.54 cm when the 50% selective length, which is a selectivity evaluation index, had a hanging ratio of 35%, a selection range of 0.51 cm, and 22.70 cm and 3.08 cm for the hanging ratio of 70%. The T45 netting was 24.34 cm and 2.13 cm, and the T90 netting was 23.51 cm and 2.84 cm. The results of the T45 netting and the T90 netting are similar, and the 50% selection length and selection range were relatively larger than the T0 netting. There was a significant difference in the correlation between the circumference of the inner circle of the mesh by the shape of the netting and the body girth of the experimental individual (Pearson test, r = 0.86, p < 0.05). There was no significant difference in the correlation between the selection ratio by the T0 netting, T45 netting, and T90 netting with a 70% hanging ratio (one-way ANOVA, p > 0.05). The results of this study showed that selectivity such as T45 netting and T90 netting appeared when the hanging ratio, which maximizes the area of T0 netting, was maintained at 70%.

• Nuclear Engineering and Technology
• /
• v.54 no.12
• /
• pp.4698-4707
• /
• 2022

As part of the ICRP Task Group 103 project, we developed ten thyroid models for the pediatric mesh-type reference computational phantoms (MRCPs). The thyroid is not only a radiosensitive target organ needed for effective dose calculation but an important source region particularly for radioactive iodines. The thyroid models for the pediatric MRCPs were constructed by converting those of the pediatric voxel-type reference computational phantoms (VRCPs) in ICRP Publication 143 to a high-quality mesh format, faithfully maintaining their original topology. At the same time, we improved several anatomical parameters of the thyroid models for the pediatric MRCPs, including the mass, overlying tissue thickness, location, and isthmus dimensions. Absorbed doses to the thyroid for the pediatric MRCPs for photon external exposures were calculated and compared with those of the pediatric VRCPs, finding that the differences between the MRCPs and VRCPs were not significant except for very low energies (<0.03 MeV). Specific absorbed fractions (target ⟵ thyroid) for photon internal exposures were also compared, where significant differences were frequently observed especially for the target organs/tissues close to the thyroid (e.g., a factor of ~1.2-~327 for the thymus as a target) due mainly to anatomical improvement of the MRCP thyroid models.

• Nuclear Engineering and Technology
• /
• v.55 no.2
• /
• pp.452-459
• /
• 2023

The core power control is an important issue for the study of dynamic characteristics in China initiative accelerator driven subcritical system (CiADS), which has direct impact on the control strategy and safety analysis process. The CiADS is an experimental facility that is only controlled by the proton beam intensity without considering the control rods in the current engineering design stage. In order to get the optimized operation scheme with the stable and reliable features, the variation of beam intensity using the continuous and periodic control approaches has been adopted, and the change of collimator and the adjusting of duty ratio have been proposed in the power control process. Considering the neutronics and the thermal-hydraulics characteristics in CiADS, the physical model for the core power control has been established by means of the point reactor kinetics method and the lumped parameter method. Moreover, the multi-inputs single-output (MISO) logical structure for the power control process has been constructed using proportional integral derivative (PID) controller, and the meta-heuristic algorithm has been employed to obtain the global optimized parameters for the stable running mode without producing large perturbations. Finally, the verification and validation of the control method have been tested based on the reference scenarios in considering the disturbances of spallation neutron source and inlet temperature respectively, where all the numerical results reveal that the optimization method has satisfactory performance in the CiADS core power control scenarios.

• Nuclear Engineering and Technology
• /
• v.55 no.4
• /
• pp.1439-1447
• /
• 2023

The fuel cycle characteristics of the IVG.1M reactor were studied within the framework of the research reactor conversion program to modernize the IVG.1M reactor. Optimum use of the nuclear fuel and reactor was achieved through routine methods which included partial fuel reloading combined with scheduled maintenance operations. Since, the additional problem in planning the fuel cycle of the IVG.1M reactor was the poisoning of the beryllium parts of the core, reflector, and control system. An assessment of the residual power and composition of spent fuel is necessary for the selection and justification of the technology for its subsequent management. Computational studies were performed using the MCNP6.1 program and the neutronics model of the IVG.1M reactor. The proposed scheme of annual partial fuel reloading allows for maintaining a high reactor reactivity margin, stabilizing it within 2-4 β_eff for 20 years, and achieving a burnup of 9.9-10.8 MW × day/kg U in the steady state mode of fuel reloading. Spent fuel immediately after unloading from the reactor can be placed in a transport packaging cask for shipping or safely stored in dry storage at the research reactor site.

• The Bulletin of The Korean Astronomical Society
• /
• v.46 no.2
• /
• pp.47.3-48
• /
• 2021

We studied the large scale dynamo process in a system forced by helical magnetic field. The dynamo process is basically nonlinear, but can be linearized with 𝛼&𝛽 coefficients and large scale magnetic field $\bar{B}$. This is very useful to the investigation of solar (stellar) dynamo. A coupled semi-analytic equations based on statistical mechanics are used to investigate the exact evolution of 𝛼&𝛽. This equation set needs only magnetic helicity ${\bar{H}}_M({\equiv}{\langle}{\bar{A}}{\cdot}{\bar{B}}{\rangle},\;{\bar{B}}={\nabla}{\times}{\bar{A}})$ and magnetic energy ${\bar{E}}_M({\equiv}{\langle}{\bar{B}}^2{\rangle}/2)$. They are fundamental physics quantities that can be obtained from the dynamo simulation or observation without any artificial modification or assumption. 𝛼 effect is thought to be related to magnetic field amplification. However, in reality the averaged 𝛼 effect decreases very quickly without a significant contribution to ${\bar{B}}$ field amplification. Conversely, 𝛽 effect contributing to the magnetic diffusion maintains a negative value, which plays a key role in the amplification with Laplacian ∇²(= - k²) for the large scale regime. In addition, negative magnetic diffusion accounts for the attenuation of plasma kinetic energy E_V(= 〈 U² 〉/2) (U: plasma velocity) when the system is saturated. The negative magnetic diffusion is from the interaction of advective term - U • ∇ B from magnetic induction equation and the helical velocity field. In more detail, when 'U' is divided into the poloidal component U_pol and toroidal one U_tor in the absence of reflection symmetry, they interact with - B • ∇ U and - U • ∇ B from ∇ × 〈 U × B 〉 leading to 𝛼 effect and (negative) 𝛽 effect, respectively. We discussed this process using the theoretical method and intuitive field structure model supported by the simulation result.

• Nuclear Engineering and Technology
• /
• v.55 no.6
• /
• pp.2112-2124
• /
• 2023

With the rise of economic and safety standards for nuclear reactors, new concepts of Gen-IV reactors and modular reactors showed more complex designs that challenge current tools for reactor physics analysis. A Monte Carlo (MC) two-step method was proposed in this work. This calculation scheme uses the continuous-energy MC method to generate multi-group cross-sections from heterogeneous models. The multi-group MC method, which can adapt locally-heterogeneous models, is used in the core calculation step. This calculation scheme is verified using a Gen-IV modular lead-based fast reactor (LFR) benchmark case. The influence of homogenized patterns, scatter approximations, flux separable approximation, and local heterogeneity in core calculation on simulation results are investigated. Results showed that the cross-sections generated using the 3D assembly model with a locally heterogeneous representation of control rods lead to an accurate estimation with less than 270 pcm bias in core reactivity, 0.5% bias in control rod worth, and 1.5% bias on power distribution. The study verified the applicability of multi-group cross-sections generated with the MC method for LFR analysis. The study also proved the feasibility of multi-group MC in core calculation with local heterogeneity, which saves 85% time compared to the continuous-energy MC.

• Proceedings of the Korean Institute of Information and Commucation Sciences Conference
• /
• 2021.05a
• /
• pp.330-331
• /
• 2021

Newly, the fourth industrial revolution is a way of describing the blurring of boundaries between the physical, digital, and biological worlds. It's a fusion of advances in AI (artificial intelligence), robotics, the IoT (internet of things), 3d printing, genetic engineering, quantum computing, and other technologies. At the world economic forum in Davos, switzerland, in january 2016, chairman professor (klaus schwab) proposed the fourth industrial revolution for the first time. In order to apply the AHP (analytic hierarchy process) analysis method, the first stage factors were designed as Digital Technology, Physics Technology and Biological Technology. In addition, the second stage factors were organized into 8 detailed services presented in the conceptual model. Thus, we present the theoretical and practical implications of these results.