• 대한조선학회논문집
• /
• 제47권5호
• /
• pp.647-655
• /
• 2010

In this paper, a numerical study is carried out for super-pipe, flat plate and axisymmetric body flows to investigate a validity of using wall function and high $y_1^+$ in calculation of high Reynolds number flow. The velocity profiles in boundary layer agree well with the law of the wall. And it is found that the range of $y^+$��which validated the logarithmic law of the wall grows with increasing Reynolds number. From the result, an equation is suggested that can be used to estimate a maximum $y^+$ value of validity of the log law. And the slope(1/$\kappa$) of the log region of the numerical result is larger than that of experimental data. On the other hand, as $y_1^+$ is increasing, both the friction and the pressure resistances tend to increase finely. When using $y_1^+$ value beyond the range of log law, the surface shear stress shows a significant error and the pressure resistance increases rapidly. However, when using $y_1^+$ value in the range, the computational result is reasonable. From this study, the use of the wall function with high value of $y_1^+$ can be justified for a full scale Reynolds number ship flow.

• 한국태양에너지학회 논문집
• /
• 제35권1호
• /
• pp.35-43
• /
• 2015

A numerical simulation on the wake flow of a wind turbine which is a scaled version of a multi-megawatt wind turbine has been performed. Two different inlet conditions of averaged wind speed including one below and one above the rated wind speed were used in the simulation. Steady-state pitch angles of the blade associated with the two averaged wind speeds were imposed for the simulation. The steady state analysis based on the Reynolds averaged Navier-Stokes equations with the method of frame motion were used for the simulation to find the torque of the rotor and the wake field behind the wind turbine. The simulation results were compared with the results obtained from the wind tunnel testing. From comparisons, it was found that the simulation results on the turbine power are pretty close to the experimental values. Also, the wake results were relatively close to the experimental results but there existed some discrepancy in the shape of velocity deficit. The reason for the discrepancy is considered due to the steady state solution with the frame motion method used in the simulation. However, the method is considered useful for solutions with much reduced calculation time and reasonably good accuracy compared to the transient analysis.

• 한국철도학회:학술대회논문집
• /
• 한국철도학회 2008년도 추계학술대회 논문집
• /
• pp.66-72
• /
• 2008

본 연구에서는 대심도 지하역사에서 화재가 발생할 시 급/배기 팬의 동작 유무에 따른 승강장에서의 열 및 연기의 거시적인 거동을 화재시뮬레이션을 통하여 분석하였다. 시뮬레이션 분석결과를 토대로 현재 설치된 급/배기 팬에 대한 제연/배연능력에 대하여 고찰하였다. 본 연구의 대상은 숭실대 입구 역사(7호선, 도시철도공사운영)이며, 숭실대 역사의 승강장은 길이 165m, 폭 23.5m, 깊이 47m 이다. 본 연구에서 전산수치해석을 위한 모델은 선로부 지하터널를 감안하여 전후 각 100m를 추가하였다. 따라서 모델링의 크기는 길이 365m, 폭23.5m, 깊이 47m 이다. 격자는 육면체 정렬격자계를 사용하였으며, 격자의 수는 대략 10,000,000 개가 사용되었다. 빠른 수치전산처리를 위하여, 병렬처리기법을 적용하였다. 본 전산수치해석에 사용된 CPU자원에 Intel 3.0GHZ Dual CPU 6개(core 12개)가 사용되었다.

• Nuclear Engineering and Technology
• /
• 제18권2호
• /
• pp.78-84
• /
• 1986

일차섭동이론을 이용하여 물질밀도 감도 계수의 표현식을 유도하였다. Super-Phenix I 평형노심의 초기상태를 기준계로 택했으며 유효중배계수를 계의 응답함수로 정의했다. 볼츠만 연산자의 구성요 소인 물질밀도로 표현되는 핵연료의 농축도와 실효밀도를 입력변화로 선정했다. 위 계산을 수행하는데 전산코드시스템 (KAERI-26군 단면적 library/1DX/2DB/PERT-V)가 사용되었다. 핵연료 농축도의 유효증배계수에 대한 감도계수는 4.576로 계산되었으며, 핵연료 실효밀도의 감도 계수는 0.0756으로 계산되었다. 본 연구는 감도해석법이 대형전산코드를 이용한 직접반복계산법에 비해 계산시간의 단축과 아울러 많은 정보를 준다는 것을 보여준다.

• Journal of Radiation Protection and Research
• /
• 제46권3호
• /
• pp.98-105
• /
• 2021

Background: In recent events of the coronavirus disease 2019 (COVID-19) pandemic, computed tomography (CT) scans are being globally used as a complement to the reverse-transcription polymerase chain reaction (RT-PCR) tests. It will be important to be aware of major organ dose levels, which are more relevant quantity to derive potential long-term adverse effect, for Korean pediatric and adult patients undergoing CT for COVID-19. Materials and Methods: We calculated organ dose conversion coefficients for Korean pediatric and adult CT patients directly from Korean pediatric and adult computational phantoms combined with Monte Carlo radiation transport techniques. We then estimated major organ doses delivered to the Korean child and adult patients undergoing CT for COVID-19 combining the dose conversion coefficients and the international survey data. We also compared our Korean dose conversion coefficients with those from Caucasian reference pediatric and adult phantoms. Results and Discussion: Based on the dose conversion coefficients we established in this study and the international survey data of COVID-19-related CT scans, we found that Korean 7-year-old child and adult males may receive about 4-32 mGy and 3-21 mGy of lung dose, respectively. We learned that the lung dose conversion coefficient for the Korean child phantom was up to 1.5-fold greater than that for the Korean adult phantom. We also found no substantial difference in dose conversion coefficients between Korean and Caucasian phantoms. Conclusion: We estimated radiation dose delivered to the Korean child and adult phantoms undergoing COVID-19-related CT examinations. The dose conversion coefficients derived for different CT scan types can be also used universally for other dosimetry studies concerning Korean CT scans. We also confirmed that the Caucasian-based CT organ dose calculation tools may be used for the Korean population with reasonable accuracy.

• Nuclear Engineering and Technology
• /
• 제53권6호
• /
• pp.1769-1785
• /
• 2021

The Integral Pressurized Water Reactors (IPWRs) as the innovative advanced and generation-III + reactors are under study and developments in a lot of countries. This paper is aimed at the thermal hydraulic study of the hot and average fuel sub-channel in a Generation III + IPWR by loose external coupling to the neutronic simulation. The power produced in fuel pins is calculated by the neutronic simulation via MCNPX2.6 then fuel and coolant temperature changes along fuel sub-channels evaluated by computational fluid dynamic thermal hydraulic calculation through an iterative coupling. The relative power densities along the fuel pin in hot and average fuel sub-channel are calculated in sixteen equal divisions. The highest centerline temperature of the hottest and the average fuel pin are calculated as 633 K (359.85 ℃) and 596 K (322.85 ℃), respectively. The coolant enters the sub-channel with a temperature of 557.15 K (284 ℃) and leaves the hot sub-channel and the average sub-channel with a temperature of 596 K (322.85 ℃) and 579 K (305.85 ℃), respectively. It is shown that the spacer grids result in the enhancement of turbulence kinetic energy, convection heat transfer coefficient along the fuel sub-channels so that there is an increase in heat transfer coefficient about 40%. The local fuel pin temperature reduction in the place and downstream the space grids due to heat transfer coefficient enhancement is depicted via a graph through six iterations of neutronic and thermal hydraulic coupling calculations. Working in a low fuel temperature and keeping a significant gap below the melting point of fuel, make the IPWR as a safe type of generation -III + nuclear reactor.

• Nuclear Engineering and Technology
• /
• 제53권12호
• /
• pp.3902-3909
• /
• 2021

SACOS series of subchannel analysis codes have been developed by XJTU-NuTheL for many years and are being used for the thermal-hydraulic safety analysis of various reactor cores. To achieve fine whole core pin-level analysis, the input preprocessing and parallel capabilities of the code have been developed in this study. Preprocessing is suitable for modeling rectangular and hexagonal assemblies with less error-prone input; parallelization is established based on the domain decomposition method with the hybrid of MPI and OpenMP. For domain decomposition, a more flexible method has been proposed which can determine the appropriate task division of the core domain according to the number of processors of the server. By performing the calculation time evaluation for the several PWR assembly problems, the code parallelization has been successfully verified with different number of processors. Subsequent analysis results for rectangular- and hexagonal-assembly core imply that the code can be used to model and perform pin-level core safety analysis with acceptable computational efficiency.

• 대한조선학회논문집
• /
• 제55권6호
• /
• pp.474-481
• /
• 2018

This paper is to compare by numerical analysis the flow characteristics and propulsion performance of stern with the shape change of K-duct, a pre-swirl duct developed by Korea Research Institute of Ships & Ocean Engineering (KRISO). First, the characteristics of the propeller and the resistance and self-propulsion before and after the attachment of the K-duct to the ship were verified and the validity of the calculation method was confirmed by comparing this result with the model test results. After that, resistance and self-propulsion calculations were performed by the same numerical method when the K-duct was changed into five different shapes. The efficiency of the other five cases was compared using the delivery horsepower in the model scale and the flow characteristics of the stern were analyzed as the velocity and pressure distributions in the area between the duct end and the propeller plane. For the computation, STAR-CCM +, a general-purpose flow analysis program, was used and the Reynolds Averaged Navier-Stokes (RANS) equations were applied. Rigid Body Motion (RBM) method was used for the propeller rotating motion and SST $k-{\omega}$ turbulence model was applied for the turbulence model. As a result, the tangential velocity of the propeller inflow changed according to the position angle change of the stator, and the pressure of the propeller hub and the cap changes. This regulated the propeller hub vortex. It was confirmed that the vortex of the portion where the fixed blade and the duct meet was reduced by blunt change.

• 대한조선학회논문집
• /
• 제55권6호
• /
• pp.490-496
• /
• 2018

The consistency of the initially designed waves in the domain is essential for accurate calculation of the added resistance in waves through CFD. In particular, unwanted reflected waves at domain boundaries can cause incorrect numerical solutions due to the superposition with initially designed waves. Euler Overlay Method(EOM) is one of the methods for reducing wave reflections by adding an additional source term to momentum and phase conservation equations, respectively. In this study, we apply the Euler Overlay Method(EOM) to the open-source CFD library, OpenFOAM(R), to simulate the accurate free-surface waves in the domain and the parametric study is performed for efficient implementation of Euler Overlay Method(EOM). Considering that the damping efficiency depends on the selection of the overlay parameter in the added source terms, the size of overlay zone and the wave steepness, the influences of these factors are tested through the wave elevation measured at constant time intervals in the 2D numerical wave tank. Through this process, guidelines for selection of optimal overlay parameter and overlay zone size that can be applied according to the scaling law are finally presented.

• 융합신호처리학회논문지
• /
• 제19권4호
• /
• pp.199-205
• /
• 2018

기존의 전영역 탐색 방법의 많은 계산량으로 인하여 예측 화질 향상과 연산량 감축을 위한 연구가 진행되어 왔으며, 본 논문에서는 전영역 탐색기반의 방법과 비교하여 예측화질은 거의 유지하면서 효율적으로 계산량을 줄이는 방법을 제안한다. 제안하는 알고리즘은 다중 문턱치 필터를 이용하여 각 후보 지점에 대하여 부분 블록 에러 합을 계산하고, 이를 여러 문턱치 필터에 적용하여 각 후보들을 영역별로 분류 또는 제거하고, 이에 따라 다음 단계에서 진행할 후보들을 선별하고, 최소 에러지점의 최적후보에 대해 단계별 부동 회수를 비교 판단하여 그 다음 단계의 진행 여부를 결정함으로써 최적의 움직임 벡터를 고속으로 계산한다. 이를 통하여 전체의 최소블록매칭에러를 갖는 움직임 벡터를 조기에 발견하고, 불필요한 후보들을 더 빨리 제거함으로써 불필요한 계산량을 줄이고 계산속도 향상을 얻을 수 있다. 또한 제안하는 알고리즘은 단독으로 사용할 뿐만 아니라 기존의 고속 알고리즘들과 결합하여 사용해도 예측화질대비 우수한 연산량 감소를 얻을 수 있으며, 실험결과에서 이를 검증한다.