• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.5 s.248
• /
• pp.430-437
• /
• 2006

The mixing effect is studied by comparing rotating and oscillating stirrers in the micro channel. The cases of Re=10 to 80 with various stirring speeds are considered to analysis the effect of Re and stirrer speed for the mixing. Under Re=20, the oscillating stirrer represents better mixing rate than the rotating stirrer up to the critical stirrer speed which has a maximum efficiency. Over Re=30, the results of oscillating and rotating stirrer show that the faster the stirrer speed, the higher the mixing effect within the concerned stirrer speed range and the oscillating stirrer keeps the higher mixing rate. It was found that the mixing effect is a function which has an optimum of the Reynolds number and the stirrer speed. The D2Q9 Lattice Boltzmann Method is used due to the merits of calculation for the unsteady flow with moving boundary.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.11
• /
• pp.1039-1048
• /
• 2008

new computational procedure for the Non-Linear Vortex Lattice Method (NLVLM) is suggested in this work. Conventional procedures suggested so far usually involves inner iteration loop to update free vortex shape and an under-relaxation based iteration loop to determine the free vortex shape. In this present work, we suggest a new formula based on quasi-steady concept to fix free vortex shape which eliminates the need for inner iteration loop. Further, the ensemble averaging of the induced velocities for a given free vortex segment evaluated at each iteration significantly improves the convergence property of the algorithm without resorting to the under-relaxation technique. Numerical experiments over several low aspect ratio wings are carried out to obtain optimal empirical parameters such as the length of the free vortex segment, the vortex core radius, and the rolled-up wake length.

• Nuclear Engineering and Technology
• /
• v.49 no.1
• /
• pp.1-5
• /
• 2017

The heavy water zero power reactor (HWZPR), which is a critical assembly with a maximum power of 100 W, can be used in different lattice pitches. The last change of core configuration was from a lattice pitch of 18-20 cm. Based on regulations, prior to the first operation of the reactor, a new core was simulated with MCNP (Monte Carlo N-Particle)-4C and WIMS (Winfrith Improved Multigroup Scheme)-CITATON codes. To investigate the criticality of this core, the effective multiplication factor ($K_{eff}$) versus heavy water level, and the critical water level were calculated. Then, for safety considerations, the reactivity worth of $D_2O$, the reactivity worth of safety and control rods, and temperature reactivity coefficients for the fuel and the moderator, were calculated. The results show that the relevant criteria in the safety analysis report were satisfied in the new core. Therefore, with the permission of the reactor safety committee, the first criticality operation was conducted, and important physical parameters were measured experimentally. The results were compared with the corresponding values in the original core.

• Nuclear Engineering and Technology
• /
• v.52 no.6
• /
• pp.1120-1130
• /
• 2020

The CREOLE experiment performed In the EOLE critical facility located In the Nuclear Center of CADARACHE - CEA have allowed us to get interesting and complete experimental information on the temperature effects in the light water reactor lattices. To analyze these experiments with accuracy an elaborate calculation scheme using the Monte Carlo method implemented in the MCNP6.1 code and the ENDF/B-VII.1 cross section library has been developed. We have used the ENDF/B-VII.1 data provided with the MCNP6.1.1 version in ACE format and the Makxsf utility to handle the data in the specific temperatures not available in the MCNP6.1.1 original library. The main purpose of this analysis is the qualification of the ENDF/B-VII.1 nuclear data for the prediction of the Reactivity Temperature Coefficient while ensuring the ability of the MCNP6.1 system to model such a complex experiment as CREOLE. We have analyzed the case of UO₂ lattice with 1166 ppm of boron in ordinary water moderator in specified temperatures. A detailed comparison of the calculated effective multiplication factors with the reference ones [1] in room temperature presented in this work shows a good agreement demonstrating the validation of our 3D calculation model. The discrepancies between calculations and the differential measurements of the Reactivity Temperature Coefficient for the analyzed configuration are relatively small: the maximum discrepancy doesn't exceed 1,1 pcm/℃. In addition to the analysis of direct differential measurements of the reactivity temperature coefficient performed in the poisoned UO₂ lattice configuration, we have also analyzed integral measurements in UO₂ clean lattice configuration using equivalency of the integral temperature reactivity worth with the driver core fuel reactivity worth and soluble boron reactivity worth. In this case both of the ENDF/B-VII.1 and JENDL.4 libraries were used in our analysis and the obtained results are very similar.

• Journal of the Korean Magnetics Society
• /
• v.23 no.1
• /
• pp.1-6
• /
• 2013

Spin-orbit coupling (SOC) effect is known to be the physical origin for various exotic magnetic phenomena in the low-dimensional systems. Recently, SOC also draws lots of attention in the study on magnetically doped thermoelectric alloys to determine their properties as the thermoelectric application as well as the topological insulator via the exact electronic structures determination near the Fermi level. In this research, aiming to investigate the spin-orbit coupling effect on the structural properties such as the lattice constants and the bulk modulus of the most widely investigated thermoelectric host material, $Bi_2Te_3$, we carried out the first-principles electronic structure calculation using the all-electron FLAPW (full-potential linearized augmented plane-wave) method. Employing both the local density approximation (LDA) and the generalized gradient approximation (GGA), the structural optimization is achieved by varying the in-plane lattice constant fixing the perpendicular lattice constant and vice versa, to find that the SOC effect increases the equilibrium lattices slightly in both directions while it markedly reduces the bulk modulus value implying the strong orientational dependence, which are attributed to the material's intrinsic structural anisotropy.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1240-1249
• /
• 2017

An essential component of the neutron transport solver is the resonance self-shielding calculation used to determine equivalence cross sections. The neutron transport code, MPACT, is currently using the subgroup self-shielding method, in which the method of characteristics (MOC) is used to solve purely absorbing fixed-source problems. Recent efforts incorporating multigroup kernels to the MOC solvers in MPACT have reduced runtime by roughly $2{\times}$. Applying the same concepts for self-shielding and developing a novel lumped parameter approach to MOC, substantial improvements have also been made to the self-shielding computational efficiency without sacrificing any accuracy. These new multigroup and lumped parameter capabilities have been demonstrated on two test cases: (1) a single lattice with quarter symmetry known as VERA (Virtual Environment for Reactor Applications) Progression Problem 2a and (2) a two-dimensional quarter-core slice known as Problem 5a-2D. From these cases, self-shielding computational time was reduced by roughly $3-4{\times}$, with a corresponding 15-20% increase in overall memory burden. An azimuthal angle sensitivity study also shows that only half as many angles are needed, yielding an additional speedup of $2{\times}$. In total, the improvements yield roughly a $7-8{\times}$ speedup. Given these performance benefits, these approaches have been adopted as the default in MPACT.

• Journal of the Korean Society of Mechanical Technology
• /
• v.13 no.2
• /
• pp.115-121
• /
• 2011

The size of hydrogen molecule is not so small as to invade into the lattice of material, and therefore, hydrogen invades into the material as atom. Hydrogen movement is done by diffusion or dislocation movement in the near crack tip or plastic deformation. Hydrogen appeared to have many effects on the mechanical properties of the Cr-Mo steel alloys. The materials for this study are 1.25Cr-0.5Mo and 2.25Cr-1Mo steels used at high temperature and pressure. The hydrogen amount obtained by theoretical calculation was almost same with the result solved by finite element analysis. The distribution of hydrogen concentration and average concentration was calculated for a flat specimen. Also, finite element analysis was employed to simulate the redistribution of hydrogen due to stress gradient. The calculation of hydrogen concentration diffused into the material by finite element method will provide the basis for the prediction of delayed fracture of notched specimen. The distribution of hydrogen concentration invaded into the smooth and notched specimen was obtained by finite element analysis. The hydrogen amount is much in smooth specimen and tends to concentrate in the vicinity of surface. Hydrogen embrittlement susceptibility of notched specimen after hydrogen charging is more remarkable than that of smooth specimen.

• Nuclear Engineering and Technology
• /
• v.10 no.4
• /
• pp.203-213
• /
• 1978

A combination of computer codes such as LATREP, HWR, AXAV and CITATION is utilized in an attempt to analyze the nuclear design characteristics of the CANDU-PHWR of the Wolsung Unit 1. The major nuclear properties to be computed are tile lattice properties of CANDU fuel channel and the core channel power distribution. The computed results are compared with the PSR documentation for the Wolsung reactor. The observed discrepancies between our computation and the PSB values are discussed in terms of incomplete information on the description of the core configuration in the PSR and the different calculation methods.

• Journal of the Korean Magnetic Resonance Society
• /
• v.10 no.1
• /
• pp.59-73
• /
• 2006

Spin-lattice relaxation pathway of $CH_2$ spin system by two dimensional NOESY sequence has been discussed. Two-dimensional spectra governed by dipolar relaxation mechanism were simulated in term of transfer mode, the generalization of conventionally used magnetization mode in one dimension. The transfer matrix directly related to the Redfield relaxation matrix can be constructed by the multiplet of transfer mode. The observable relaxation transfer modes causes to variation of the off-diagonal signal intensity of phase sensitive NOESY spectra from which variable spectral density can be extracted with simple group theoretical calculation. The variation of the J-coupling peak intensity as a function of the mixing time in 2-D spectra for $n-Undecane-5-^{13}C$ and Bromoacetic $2-^{13}C$ acid has been theoretically traced.

• Proceedings of the KIEE Conference
• /
• 1992.07b
• /
• pp.580-583
• /
• 1992

In calculating the hysteresis loops with the Preisach model for ferromagnetism, Everett function is used generally. Because the Everett function is usually given as a table of the lattice-shaped, it is very difficult to directly obtain the Everest table from the data. Therefore this makes some defects in the calculation processes or the accuracy of the results. In this study, using the data sufficiently obtained from the experiment by drawing up the Everett table in the triangle-shaped, and applying the generalized hysteresis model in which the magnetization is depend on the sum of the applied magnetic field and the molecular field, it is shown that our proposal is acceptible in calculating the hysteresis processes.