• Journal of KIISE
• /
• v.43 no.3
• /
• pp.304-313
• /
• 2016

We present a novel memory-efficient parallel ray casting algorithm for unstructured grid volume rendering on multi-core CPUs. Our method is based on the Bunyk ray casting algorithm. To solve the high memory overhead problem of the Bunyk algorithm, we allocate a fixed size local buffer for each thread and the local buffers contain information of recently visited faces. The stored information is used by other rays or replaced by other face's information. To improve the utilization of local buffers, we propose an image-plane based ray grouping algorithm that makes ray groups have high coherency. The ray groups are then distributed to computing threads and each thread processes the given groups independently. We also propose a novel hash function that uses the index of faces as keys for calculating the buffer index each face will use to store the information. To see the benefits of our method, we applied it to three unstructured grid datasets with different sizes and measured the performance. We found that our method requires just 6% of the memory space compared with the Bunyk algorithm for storing face information. Also it shows compatible performance with the Bunyk algorithm even though it uses less memory. In addition, our method achieves up to 22% higher performance for a large-scale unstructured grid dataset with less memory than Bunyk algorithm. These results show the robustness and efficiency of our method and it demonstrates that our method is suitable to volume rendering for a large-scale unstructured grid dataset.

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.5
• /
• pp.36-42
• /
• 2003

A CFD simulation technique has been developed to handle the unsteady body motion with large amplitude by use of overlapping multi-block grid system. The three-dimensional, viscous and incompressible flow around body is investigated by solving the Navier-Stokes equations, and the motion of body is represented by moving effect of the grid system. Composite grid system is employed in order to deal with both the body motion with large amplitude and the condition of numerical wave maker in convenience at the same time. The governing equations, Navier-Stokes (N-S) and continuity equations, are discretized by a finite volume method, in the framework of an O-H type boundary-fitted grid system (inner grid system including test model) and a rectangular grid system (outer grid system including simulation equipments for generation of wave environments). If this study, several flow configurations, such as an oscillating cylinder with large KC number, are studied in order to predict and evaluate the hydrodynamic forces. Furthermore, the motion simulation of a Series 60 model advancing in a uniform flow under the condition of enforced roll motion of angle 20$^{\circ}$ is performed in the developed numerical wave tank.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.483-487
• /
• 2011

In the present study, a CFD program based on a finite volume method was developed by using an unstructured polyhedral grid system for the accurate simulation with the complex geometry of computational domain. To simulate the transient flow induced by the moving solid object, the program used a fractional step method and a ALE (Algebric Lagrangian-Eulerian) method. The grid deformation for the moving of solid object were performed with a spring analogy based on the center coordinate of each computational grid. To verify the present program with these methodologies, the numerical results of the flow around the fixed and oscillating circular cylinder were compared with the previous numerical results.

• The Journal of the Korea Contents Association
• /
• v.9 no.12
• /
• pp.877-885
• /
• 2009

In this paper, we propose the SINDI-Grid which is a high-performance framework for scientific and technological knowledge discovery using the grid computing. By using the advantages of the grid computing providing data repository of large-volume and high-speed computing power, the SINDI-Grid framework provides a variety of grid services for distributed data analysis and scientific knowledge processing. And the SINDI-Workflow tool exploits these services so that performs the design and execution for scientific and technological knowledge discovery applications which integrate various information processing algorithms.

• Journal of Korean Port Research
• /
• v.15 no.1
• /
• pp.87-97
• /
• 2001

The calculation of earthwork plays a major role in plan or design of many civil engineering projects, and thus it has become very important to advanced the accuracy of earthwork calculation. Current method used for estimating the volume of pit excavation assumes that the ground profile between the grid points is linear(trapezoidal rule), or nonlinear(simpson's formulas). In this paper the spot height method, least square method, and chamber formulas, Chen and Lin method are compared with the volumes of the pits in these examples. As a result of this study, algorithm of chen and Lin me쇙 by spline method should provide a better accuracy than the spot height method, least square method, chamber formulas. The Chen and Lin formulas can be used for estimating the excavation volume of a pit divide into a grid with unequal intervals. From the characteristics of the cubic spline polynomial, the modeling curve of the Chen and Lin method is smooth and matches the ground profile well. Generally speaking, the nonlinear profile formulas provide better accuracy than the linear profile formulas. The mathematical model mentioned make an offer maximum accuracy in estimating the volume of a pit excavation.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.1
• /
• pp.74-85
• /
• 2002

The numerical experiment has been conducted to investigate the unsteady shock wave reflecting phenomena. The cell-vertex finite-volume, Roe's upwind flux difference splitting method with unstructured grid is implemented to solve unsteady Euler equations. The $4^{th}$-order Runge-Kutta method is applied for time integration. A linear reconstruction of the flux vector using the least-square method is applied to obtain the $2^{nd}$-order accuracy for the spatial derivatives. For a better resolution of the shock wave and slipline, the dynamic grid adaptation technique is adopted. The new concept of grid adaptation technique, which is much simpler than that of conventional techniques, is introduced for the current study. Three error indicators (divergence and curl of velocity, and gradient of density) are used for the grid adaptation procedure. Considering the quality of the solution and the numerical efficiency, the grid adaptation procedure was updated up to $2^{nd}$ level at every 20 time steps. For the convenience of comparison with other experimental and analytical results, the case of interaction between the straight incoming shock wave and a sharp wedge is simulated for various flow conditions. The numerical results show good agreement with other experimental and analytical results, in the shock wave reflecting structure, slipline, and the trajectory of the triple points. Some critical cases show disagreement with the analytical results, but these cases also have been proven to show hysteresis phenomena.

• 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 KIISE:Computing Practices and Letters
• /
• v.15 no.8
• /
• pp.525-533
• /
• 2009

Staging technique is used to provide files for a job in the Grid. If a staged file has large volume, the start time of the job is delayed and the throughput of job in the Grid may decrease. Therefore, removal of staging overhead helps the Grid operate more efficiently. In this paper, we present two methods for efficient file provisioning to clear the overhead. First, we propose RA-RFT, which extends RFT of Globus Toolkit and enables it to utilize RLS with replica information. RA-RFT can reduce file transfer time by doing partial transfer for each replica in parallel. Second, we suggest Remote Link that uses remote I/O instead of file transfer. Remote link is able to save storage of computational nodes and enables fast file provisioning via prefetching. Through various experiments, we argue that our two methods have an advantage over existing staging techniques.

• Nuclear Engineering and Technology
• /
• v.53 no.4
• /
• pp.1119-1126
• /
• 2021

Spacer grid with mixing vane had been widely used in nuclear reactor core. One of the main feather of spacer grid with mixing vane was that strong swirl flow was formed after the spacer grid. The swirl flow not only changed the bubble generation in the near wall field, but also affected the bubble behaviors in the center region of the subchannel. The interaction between bubble and the swirl flow was one of the basic phenomena for the two phase flow modeling in fuel assembly. To obatin better understanding on the bubble behaviors in swirl flow, full three dimension numerical simulations were conducted in the present paper. The swirl flow was assumed in the cylindral calculation domain. The bubble interface was captured by Volume Of Fluid (VOF) method. The properties of saturated water and steam at different pressure were applied in the simulation. The bubble trajectory, motion, shape and force were obtained based on the bubble parameters captured by VOF. The simulation cases in the present study included single bubble with different size, at different angular velocity conditions and at different pressure conditions. The results indicated that bubble migrated to the center in swirl flow with spiral motion type. The lateral migration was mainly related to shear stress magnitude and bubble size. The bubble moved toward the center with high velocity when the swirl magnitude was high. The largest bubble had the highest lateral migration velocity in the present study range. The effect of pressure was small when bubble size was the same. The prelimenery simulation result would be beneficial for better understanding complex two phase flow phenomena in fuel assembly with spacer grid.

• Journal of the Society of Naval Architects of Korea
• /
• v.44 no.1 s.151
• /
• pp.16-25
• /
• 2007

Studies of unsteady-airfoil flows have been motivated mostly by efforts to avoid. or reduce such undesirable effects as flutter, noise and vibrations, dynamic stall. In this paper, we carry out a computational study of viscous flows around a two-dimensional oscillating airfoil to investigate unsteady effects in these important and challenging flows. A fully implicit incompressible RANS solver has been used for calculating unsteady viscous flows around an airfoil. The cell-centered End order finite volume method is utilized to discretize governing equations. in order to ease the flow computation for fluid region changing in time, improve the qualify of solution and simplify the grid generation for an oscillating airfoil flow, the computational method adopts a moving and deforming grid generation technique based on the multi-block grid topology. The numerical method is applied for calculating viscous flows of an oscillating NACA 0012 in uniform flow. The computational results are compared with available experimental data. Computed results are compared with experimental data and flow characteristics of the experiment are reproduced well In the computed results.