• Journal of the Korean Geotechnical Society
• /
• v.29 no.12
• /
• pp.35-44
• /
• 2013

In this study, a particle method without using grid was applied for analysing large deformation problems in soil flows instead of using ordinary finite element or finite difference methods. In the particle method, a continuum equation was discretized by various particle interaction models corresponding to differential operators such as gradient, divergence, and Laplacian. Soil behavior changes from solid to liquid state with increasing water content or external load. The Mohr-Coulomb failure criterion was incorporated into the particle method to analyze such three-dimensional soil behavior. The yielding and hardening behavior of soil before failure was analyzed by treating soil as a viscous liquid. First of all, a sand column test without confining pressure and strength was carried out and then a self-standing clay column test with cohesion was carried out. Large deformation from such column tests due to soil yielding or failure was used for verifying the developed particle method. The developed particle method was able to simulate the three-dimensional plastic deformation of soils due to yielding before failure and calculate the variation of normal and shear stresses both in sand and clay columns.

• Journal of the Korea Computer Graphics Society
• /
• v.2 no.2
• /
• pp.1-9
• /
• 1996

It is important issue to increase the efficiency of visibility operations in rendering since the most of rendering time is spent for performing visibility operations. In this paper, we present light buffers of area light sources, which provide efficient visibility operations from the sources. The light buffers are extended to support visibility operations not originated from the sources. We also explore their applications such as shadow generation, ray tracing and radiosity. Our experimental results show that the light buffers support more efficient visibility operations than the previously developed data structures including BSP trees, linear grids, and octrees.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.8
• /
• pp.1122-1133
• /
• 2003

Numerical characteristics of implicit upwind schemes, such as upwind ADI, line Gauss-Seidel (LGS) and point Gauss-Seidel (LU) algorithms, for Navier-Stokes equations have been investigated. Time-derivative preconditioning method was applied for efficient convergence at low Mach/Reynolds number regime as well as at large grid aspect ratios. All the algorithms were expressed in approximate factorization form and von Neumann stability analysis was performed to identify stability characteristics of the above algorithms in the presence of high grid aspect ratios. Stability analysis showed that for high aspect ratio computations, the ADI and LGS algorithms showed efficient damping effect up to moderate aspect ratio if we adopt viscous preconditioning based on min-CFL/max-VNN time-step definition. The LU algorithm, on the other hand, showed serious deterioration in stability characteristics as the grid aspect ratio increases. Computations for several practical applications also verified these results.

• Journal of computational fluids engineering
• /
• v.19 no.4
• /
• pp.61-67
• /
• 2014

Since the most of the existing non-Newtonian models are not adequate to apply to the lattmce Boltzmann method, it is a challenging task from both the theoretical and the numerical points of view. In this research the hydro-kinetic model was modified and applied to the 3-D moving sphere in the circular channel flow and the characteristics of the shear thinning effect by the HK-model was evaluated and the condition of ${\Gamma}$ in the model was suggested for the stable simulation to generate non-trivial prediction in three dimension strong shear flows. On the wall boundaries of circular channel the curved wall surface treatment with constant velocity condition was applied and the bounceback condition was applied on the sphere wall to simulate the relative motion of the sphere. The condition is adequate at the less blockage than 0.7 but It may need to apply a multi-scale concept of grid refinement at the narrow flow region. to obtain the stable numerical results.

• Journal of computational fluids engineering
• /
• v.19 no.4
• /
• pp.20-28
• /
• 2014

An efficient solution algorithm for simulating free surface problem is presented. Navier-Stokes equations for variable density incompressible flow are employed as the governing equation on Cartesian meshes. In order to describe the free surface motion efficiently, VOF(Volume Of Fluid) method utilizing THINC(Tangent of Hyperbola for Interface Capturing) scheme is employed. The most time-consuming part of the current free surface flow simulations is the solution step of the linear system, derived by the pressure Poisson equation. To solve a pressure Poisson equation efficiently, the PCG(Preconditioned Conjugate Gradient) method is utilized. This study showed that the proper application of the preconditioner is the key for the efficient solution of the free surface flow when its pressure Poisson equation is solved by the CG method. To demonstrate the efficiency of the current approach, we compared the convergence histories of different algorithms for solving the pressure Poisson equation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.4
• /
• pp.395-405
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among four filling patterns at each triangular control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. Using the proposed numerical technique, the collapse of a water dam has been analyzed to predict flow phenomenon of fluid and the predicted front positions with respect to time have been compared with the reported experimental results.

• Journal of the Korean Society for Aviation and Aeronautics
• /
• v.18 no.1
• /
• pp.33-38
• /
• 2010

Developments of numerical methods are very important to design and analysis for a high subsonic turbine blade. In general, Analysis by experimental investigation has needed a lot of human resources and required time, indispensably, and equipments still have a limit to measure in conditions of high temperature. Rapid technological developments of CPU and integration level of memory make it possible to advance computation with almost exactly simulation so, recent developments of numerical methods are in spotlight. In the present study, the panel method, which is well-known as relatively simplified numerical method, and 2-dimensional ordinary differential Falkner-Skan equation were computed in order to analyze the outer flow, and FVM-based solid heat transfer equation, was also computed to forecast the temperature distribution of the airfoil and the turbine blade. Unstructured grid was constructed in the turbine blade, which has double cooling holes, in order to analyze the internal heat transfer. Cooling fluid was assumed as fully-developed turbulent flow and that circulated in cooling holes.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.126-131
• /
• 2004

The interaction between the hypersonic free stream and the side jet flow at high altitudes is investigated by direct simulation Monte Carlo(DSMC) method. Since there is a great difference in density between the free stream and the side jet flow, the weighting factor technique which could control the number of simulation particles, is applied to calculate these two flows simultaneously. Chemical reactions are not considered in the calculation. For validation, the corner flow passing between a pair of plates that are perpendicularly attached is solved. The side jet flow is then injected into this comer flow and solution is found for the merged flow. Results are compared with the experiments. For a more realistic rocket model, the flow past a blunted cone cylinder shape is solved. The leeward or windward jet injection is merged with this flow. The effect on the rocket surface is observed at various flow angles. The lambda effect and the wake structure are found like low attitudes. High interaction between the free stream and the side jet flow is observed when the side jet is injected in the windward direction.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.2 no.1
• /
• pp.26-37
• /
• 1994

In this study, flow of a model intake port/valve system is analyzed by using low Reynolds number $k-{\varepsilon}$ model. Discharge coefficient was obtained from computational results for the various cases of valve lifts. Discharge coefficient becomes maximum when the valve lift is 20mm, and does not increase or decrease in proportional to valve lift. Most of pressure drop and production of turbulent kinetic energy occur at the edge points of the valve and the valve seat Thus, in order to improve discharge coefficient, rounding of edge points in valve and valve seat is recommended. As valve lift is increased, the velocity of the intake jet in the valve passage decreases, and the direction of the jet is more inclined toward the valve seat.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.11a
• /
• pp.78-84
• /
• 2011

CFD simulation has been conducted on a small scaled 4 nozzle clustered engine operating with air. In the present paper, the effects of grid size, turbulence models, flux difference methods have been compared. The results show that the base flows are somewhat different as the turbulence models, while Roe and AUSM flux differences produced almost the same results. Spalart-Allmaras turbulence model produces more accurate results rather than famous SST k-w model. The calculated Mach number and pressure profile in the engine base reveal the complex base flow structure, which is somewhat different from the generally estimated flow fields.