• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.4
• /
• pp.58-65
• /
• 1994

It is very important to understand characteristics of solution due to the variation of computational grid sizes, especially when turbulence model not incorporating wall-function is used. The present paper performs numerical investigation on the grid dependency of numerical solution for three dimensional turbulent flow field around a ship. In the present study a finite volume method with a modified sub-grid scale turbulence model and a numerically constructed non-orthogonal curvilinear coordinate system capable of conforming complex ship geometries are used. Numerical studies are then performed for a mathematical Wigley hull and the Series 60, $C_B=0.8$ hull forms. The results for various grid sizes are compared with each other and with measured data to show grid dependencies of numerical solutions.

• Journal of radiological science and technology
• /
• v.36 no.2
• /
• pp.123-129
• /
• 2013

This study investigated the rate of setup variance by the rotating unbalance of gantry in image-guided radiation therapy. The equipments used linear accelerator(Elekta Synergy TM, UK) and a three-dimensional volume imaging mode(3D Volume View) in cone beam computed tomography(CBCT) system. 2D images obtained by rotating $360^{\circ}$and $180^{\circ}$ were reconstructed to 3D image. Catpan503 phantom and homogeneous phantom were used to measure the setup errors. Ball-bearing phantom was used to check the rotation axis of the CBCT. The volume image from CBCT using Catphan503 phantom and homogeneous phantom were analyzed and compared to images from conventional CT in the six dimensional view(X, Y, Z, Roll, Pitch, and Yaw). The variance ratio of setup error were difference in X 0.6 mm, Y 0.5 mm Z 0.5 mm when the gantry rotated $360^{\circ}$ in orthogonal coordinate. whereas rotated $180^{\circ}$, the error measured 0.9 mm, 0.2 mm, 0.3 mm in X, Y, Z respectively. In the rotating coordinates, the more increased the rotating unbalance, the more raised average ratio of setup errors. The resolution of CBCT images showed 2 level of difference in the table recommended. CBCT had a good agreement compared to each recommended values which is the mechanical safety, geometry accuracy and image quality. The rotating unbalance of gentry vary hardly in orthogonal coordinate. However, in rotating coordinate of gantry exceeded the ${\pm}1^{\circ}$ of recommended value. Therefore, when we do sophisticated radiation therapy six dimensional correction is needed.

• Journal of the Institute of Electronics Engineers of Korea SP
• /
• v.44 no.3
• /
• pp.48-60
• /
• 2007

In previous researches for iris feature extraction, they transform a original iris image into rectangular one by stretching and interpolation, which causes the distortion of iris patterns. Consequently, it reduce iris recognition accuracy. So we are propose the method that extracts iris feature by using polar coordinates without distortion of iris patterns. Our proposed method has three strengths compared with previous researches. First, we extract iris feature directly from polar coordinate circular iris image. Though it requires a little more processing time, there is no degradation of accuracy for iris recognition and we compares the recognition performance of polar coordinate to rectangular type using by Hamming Distance, Cosine Distance and Euclidean Distance. Second, in general, the center position of pupil is different from that of iris due to camera angle, head position and gaze direction of user. So, we propose the method of iris feature detection based on polar coordinate circular iris region, which uses pupil and iris position and radius at the same time. Third, we overcome override point from iris patterns by using polar coordinates circular method. each overlapped point would be extracted from the same position of iris region. To overcome such problem, we modify Gabor filter's size and frequency on first track in order to consider low frequency iris patterns caused by overlapped points. Experimental results showed that EER is 0.29%, d' is 5,9 and EER is 0.16%, d' is 6,4 in case of using conventional rectangular image and proposed method, respectively.

• Journal of Astronomy and Space Sciences
• /
• v.12 no.1
• /
• pp.78-89
• /
• 1995

We solved n-body problem about 9 planets, moon, and 4 minor planets with relativistic effect related to the basic equation of motion of the solar system. Perturbations including flgure potential of the earth and the moon and solid earth tidal effect were considered on this relativistic equation of motion. The orientations employed precession and nutation for the earth, and lunar libration model with Eckert's lunar libration model based on J2000.0 were used for the moon. Finally, we developed heliocentric ecliptic position and velocity of each planet using this software package named the SSEG (Solar System Ephemerides Generator) by long-term (more than 100 years) simulation on CRAY-2S super computer, through testing each subroutine on personal computer and short-time(within 800 dyas) running on SUN3/280 workstation. Epoch of input data JD2440400.5 were adopted in order to compare our results to the data archived from JPL's DE 200 by Standish and Newhall. Above equation of motion was integrated numerically having 1-day step-size interval through 40,000 days (about 110 years long) as total computing interval. We obtained high-precision ephemerides of the planets with maximum error, less $than\pm2\times10^{-8}AU(\approx\pm3km)$ compared with DE200 data (except for mars and moon).

• Proceeding of EDISON Challenge
• /
• 2013.04a
• /
• pp.354-359
• /
• 2013

일반적으로 건축물의 설계시 풍동 실험을 통한 풍환경의 평가를 수행하고 있으며, 이는 환경 영향 평가법에서 정한 건축 사업 시행 시 수반되어야 할 자연환경, 생활환경 그리고 사회경제환경의 영향 평가의 일환으로 실시되고 있다. 그러나, 풍동 실험의 경우 여러 가지 현실적 제약조건으로 설계와 실험의 피드백 (Feedback)이 원활하지 못하며, 특히 대상 건축물이 공장과 같이 대기 오염원이 되는 경우 실험은 더욱 어려운 형편이다. 이에 대한 보완책으로 전산 유체 역학을 이용한 건축물의 풍압 해석에 의한 풍하중 추정이나 인접 지형-지물의 영향을 고려한 건축물 주위의 풍환경 평가가 있다. 전산 모사에 의해 풍동 실험의 미비점을 보완하고, 보다 상세한 정보를 확보함으로써 건축물의 구조적 안전성의 증대와 환경 피해 감소를 기할 수 있다. 그러나 복잡한 지형-지물이나 건축물 주위의 풍환경에 대한 전산 모사는 주로 두 가지의 기술적 어려움을 수반하게 되다. 그 중 하나는 고정 경계면을 이루는 형상의 복잡성으로 인해 기존에 많이 이용하고 있는 Body-fitted 격자계를 이용하는 경우, 격자 생성 과정이 매우 복잡하고 어려울 뿐 만 아니라 생성된 격자가 주로 비정렬 (unstructured) 특성을 갖게 되어 수치해석 과정의 효율을 저하시키는 요인이 되며, 격자의 형상도 수치해석의 수렴성을 저하시키는 예가 많다. 다른 어려움으로 풍환경은 전형적인 난류 유동장으로서 난류의 전산 해석은 아직도 해결하지 못한 부분이 많다는 점이다. 이에 본 논문에서는 복잡한 지형-지물이나 건축물의 풍하중과 풍환경의 전산 모사 기술 확보를 위하여 수행중인 연구의 일환으로 물체 형상의 기하학적 복잡성의 극복을 위한 가상경계법 (Immersed Boundary Method)과 난류 유동장의 물리적 엄밀성을 높이기 위한 다와동 모사 (Large Eddy Simulation)을 이용한 물체 형상과 무관한 유동장 해석 기술 개발에 대하여 다루고자 한다. 먼저 최근에 유동 해석에 이용되는 방법인 가상경계법(IBM)은 물체를 포함한 전체 전산 영역을 직교 좌표계에 의해 이산화하고, 유동장내 존재하는 물체의 표면에서의 점착 조건을 만족시키기 위하여 지배 방정식에 적절한 외력을 추가로 고려하는 방법이다. 본 연구에서는 가상경계법을 이용하여 경계층에 위치한 건물 형상의 각진 물체 주위 사이에 형성되는 공동 내부의 비정상 유속 및 압력에 대한 전산 해석을 수행하고, 풍상측 전면에 형성되는 경계층에 의한 영향을 분석하였다.

• Journal of the Korean GEO-environmental Society
• /
• v.25 no.4
• /
• pp.13-20
• /
• 2024

In this paper, a parallel analysis algorithm for Smoothed Particle Hydrodynamics (SPH), one of the numerical methods for fluidic materials, is introduced. SPH, which is a meshless method, can represent the behavior of a continuum using a particle-based approach, but it demands substantial computational resources. Therefore, parallel analysis algorithms are essential for SPH simulations. The domain decomposition algorithm, which divides the computational domain into partitions to be independently analyzed, is the most representative method among parallel analysis algorithms. In Discrete Element Method (DEM) and Molecular Dynamics (MD), the Cartesian coordinate-based domain decomposition method is popularly used because it offers advantages in quickly and conveniently accessing particle positions. However, in SPH, it is important to share particle information among partitioned domains because SPH particles are defined based on information from nearby particles within the smoothing length. Additionally, maintaining CPU load balance is crucial. In this study, a highly parallel efficient algorithm is proposed to dynamically minimize the size of orthogonal domain partitions to prevent excess CPU utilization. The efficiency of the proposed method was validated through numerical analysis models. The parallel efficiency of the proposed method is evaluated for up to 30 CPUs for fluidic models, achieving 90% parallel efficiency for up to 28 physical cores.

• Journal of the Korea Society of Computer and Information
• /
• v.13 no.7
• /
• pp.117-126
• /
• 2008

Recent advances in medical imaging technologies have enabled the high-resolution data acquisition. Therefore visualization of such large data set on standard graphics hardware became a popular research theme. Among many visualization techniques, we focused on bricking method which divided the entire volume into smaller bricks and rendered them in order. Since it switches bet\W8n bricks on main memory and bricks on GPU memory on the fly, to achieve better performance, the number of these memory swapping conditions has to be minimized. And, because the original bricking algorithm was designed for regular volume data such as CT and MR, when applying the algorithm to ultrasound volume data which is based on the toroidal coordinate space, it revealed some performance degradation. In some areas near bricks' boundaries, an orthogonal viewing ray intersects the single brick twice, and it consequently makes a single brick memory to be uploaded onto GPU twice in a single frame. To avoid this redundancy, we divided the volume into bricks allowing overlapping between the bricks. In this paper, we suggest the formula to determine an appropriate size of these shared area between the bricks. Using our formula, we could minimize the memory bandwidth. and, at the same time, we could achieve better rendering performance.

• Journal of Korea Water Resources Association
• /
• v.45 no.4
• /
• pp.349-360
• /
• 2012

A horizontally curvilinear non-hydrostatic free surface model that was applicable to three-dimensional viscous flows was developed. The proposed model employed a top-layer equation to close kinematic free-surface boundary condition, and an isotropic k-${\varepsilon}$ model to close turbulence viscosity in the Reynolds averaged Navier-Stokes equation. The model solved the governing equations with a fractional step method, which solved intermediate velocities in the advection-diffusion step, and corrects these provisional velocities by accounting for source terms including pressure gradient and gravity acceleration. Numerical applications were implemented to the wind-driven currents in a two-dimensional closed basin, the flow in a steep-sided trench, and the flow in a strongly-curved channel accounting for secondary current by the centrifugal force. Through the numerical simulations, the model showed its capability that were in good agreement with experimental data with respect to free surface elevation, velocity, and turbulence characteristics.

• The Journal of the Korea Contents Association
• /
• v.5 no.4
• /
• pp.94-102
• /
• 2005

In this paper, three axial accelerometer was used to develop a small sensor module, which was attached to human body to calculate the acceleration in gravity direction by human motion, when it was positioned in any direction. To measure its wearer's walking or running motion using the sensor module, the acquired sensor data was pre-processed to enable its quantitative analysis. The acquired digital data was transformed to orthogonal coordinate value in three dimension and calculated to be single scalar acceleration data in gravity direction and normalized to be physical unit value. The normalized sensor data was used to detect walking pattern and calculate their step counts. Developed algorithm was implemented in the form of PDA application. The accuracy of the developed sensor to detect step count was about 97% in laboratory experiment.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.2
• /
• pp.204-213
• /
• 2000

This study deals with jet impingement, which is extensively used in the process industries to achieve intense heating, cooling or drying rates and also widely employed as a test flow for turbulent models due to its complex flow configuration, on a flat plate by numerical methods. In this calculation, the finite volume method was employed to solve the Navier-stokes equation based on the non-orthogonal coordinate with non-staggered variable arrangement. To get a better understanding for the fluid flow and heat transfer characteristics of the turbulent jet impingements, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model was adapted and compared with the experimental data and the result of standard $k-{\varepsilon}$ turbulent model. Numerical calculations were carried out with various flow rates, nozzle to plate distances. In the case of the axisymmetric jet impingement on a flat plate, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model showed better agreement with the experimental data than the standard $k-{\varepsilon}$ turbulent model in the prediction of the mean velocity profiles, the turbulent velocity profiles. the turbulent shear stress and the heat transfer rate. The highest heat transfer rate can be obtained when the impingement occurs within the potential core..