• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.3
• /
• pp.253-258
• /
• 2012

The Cartesian grid system has generally been used in casting simulations, even though it does not represent sloped and curved surfaces very well. These distorted boundaries cause several problems, and special treatment is necessary to resolve them. A cut cell method on a Cartesian grid has been developed for the simulation of threedimensional mold filling. Cut cells at a cast/mold interface are generated on Cartesian grids, and the governing equations are computed using the volume and areas of the cast at the cut cells. In this paper, we propose a new method based on the partial cell treatment (PCT) that can consider the cutting cells which are cut by the cast and the mold. This method provides a better representation of the surface geometry, and will be used in the computation of velocities that are defined on the cell boundaries in the Cartesian gird system. Various test examples for several casting process are computed and validated.

• Journal of Korea Foundry Society
• /
• v.43 no.6
• /
• pp.302-309
• /
• 2023

In general, castings often have complex shapes and significant variations in thickness within a single product, making grid generation for simulations challenging. Casting flows involve multiphase flows, requiring the tracking of the boundary between air and molten metal. Additionally, considerable time is spent calculating pressure fields due to density differences in a numerical analysis. For these reasons, the Cartesian grid system has traditionally been used in mold filling simulations. However, orthogonal grids fail to represent shapes accurately, leading to a momentum loss caused by the stair-like grid patterns on curved and sloped surfaces. This can alter the flow of molten metals and result in incorrect casting process designs. To address this issue, simulations in the Cartesian grid system involve creating a large number of grids to represent shapes more accurately. Alternatively, the Cut Cell method can be applied to address the problems arising from the Cartesian grid system. In this study, analysis results based on the number of grid in the Cartesian grid system for a casting flow analysis were compared with results obtained using the Cut Cell method. Casting flow simulations of actual products during various casting processes were also conducted, and these results were analyzed with and without applying the Cut Cell method.

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

A two-dimensional hybrid flaw solver has been developed for the accurate and efficient simulation of steady and unsteady flaw fields. The flow solver was cast to accommodate two different topologies of computational meshes. Triangular meshes are adopted in the near-body region such that complex geometric configurations can be easily modeled, while adaptive Cartesian meshes are, utilized in the off-body region to resolve the flaw more accurately with less numerical dissipation by adopting a spatially high-order accurate scheme and solution-adaptive mesh refinement technique. A chimera mesh technique has been employed to link the two flow regimes adopting each mesh topology. Validations were made for the unsteady inviscid vol1ex convection am the unsteady turbulent flaws over an NACA0012 airfoil, and the results were compared with experimental and other computational results.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.10
• /
• pp.17-25
• /
• 2005

A grid adaptation method within systems of chimera and patched grids is presented. Problem domains are divided into near-body and off-body fields. Near-body field is filled with curvilinear body-fitted grids that extend only a short distance from body surfaces and connected to other grid systems via chimera domain connectivity method. Off-body field is filled with patched uniform cartesian grids of varying levels of refinement. This method gives flexibility in grid generation and efficient adaptation capability. Several numerical experiments including 2D store separation were performed to show the performance of the proposed adaptation method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.5
• /
• pp.1158-1174
• /
• 1988

A numerical technique is developed for the solution of fully developed turbulent recirculating flow in the passage of variable area using the non-orthogonal coordinate transformation. In the numerical analysis, primitive pressure-velocity finite difference equations were solved by SIMPLER algorithm with 2-equation turbulence model and algebraic stress model (ASM). QUICK scheme on the differencing of convective terms which is free from the inaccuracies of numerical diffusion has been applied to the variable grids and the results compared with those from HYBRID scheme. In order to test the effect of streamline curvatures on turbulent diffusion Lee and Choi streamline curvature correction model which has been obtained by modifying the Leschziner and Rodi's model is testes. The ASM was also employed and the results are compared to those from another turbulence model. The results show that difference of convective differencing schemes and turbulence models give significant differences in the prediction of velocity fields in the expansion region and outlet region of the combustor, however show little differences in the parallel flow region.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.3
• /
• pp.561-565
• /
• 1988

In the present study, a new method of generating a nearly orthogonal boundary-fitted coordinate systems with automatic grid spacing control is introduced. Applications of the method to a two dimensional simply-connected region is then demonstrated. The nearly orthogonal boundary-fitted method has the following features, (a) Strong grid control in the .eta.-direction can be made, (b) The generated boundary-fitted coordinates are nearly orthoronal, (c) Both the .xi.-and .eta.-direction control function are mathematically derived. Especially the .eta.-direction control function is derived under the assumption that the .eta.-direction grid spacing is by far smaller than the .xi.-direction grid spacing when the .eta.-direction grid line is strongly clustered. (d) The grid control functions are dynamically adjusted by the metric scale factors imposed on the boundary. The control function is fully automatic and eliminates the need of user manipulation of the control function.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.402-407
• /
• 2006

본 연구는 곡선좌표계에서 유한차분기법(finite difference method)을 이용하여 2차원 흐름이 모의가능한 수치모형을 개발하는 것이다. 기존의 연구는 대부분 직교좌표계(cartesian coordinate system)에서의 격자망을 대상으로 개발되고 적용되었기 때문에 불규칙한 흐름의 경계 및 형상을 올바로 표현하기 어려웠다. 유한요소법이나 유한체적법같은 수치모의기법들이 개발되어 비구조격자체계를 구성하고 자연현상에 가까운 경계 표현할 수 있도록 개발되었다. 하지만 위의 기법들은 질량과 운동량과 같은 물리량을 보존하기 위해서 매우 조밀한 격자체계를 가져야만 한다. 이에 본 연구에서는 기존의 문제점들을 해결하기 위하여 곡선좌표계(curvilinear coordinate system)를 이용하여 지배방정식을 표현하고 2차원 흐름을 모의할 수 있는 모형을 구축한다. 수치모형은 leap-frog기법과 1차 정확도의 풍상차분기법(upwind scheme)을 사용하여 구성하였다. 본 연구에서 개발된 모형을 사각수조 및 만곡수로흐름에 적용하여 모의결과를 해석해 및 실험관측값과 비교하였다. 이로부터 본 수치모형이 해석해 및 실측치와 잘 일치하고 있음을 알 수 있었다.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.6
• /
• pp.561-567
• /
• 2011

The inverse radiation analysis of a two-dimensional irregular configuration using unstructured triangular meshes is presented. In this study, an enclosure filled with an absorbing, emitting and scattering medium with diffusely emitting and reflecting opaque boundaries is considered. The finite volume method is applied to solve the radiative transfer equation in order to simulate the measured incident radiation values which are used as input data for the inverse analysis. The conjugate gradient method is adopted for the estimation of wall emissivities by minimizing the objective function at each iteration step. To verify the performance of the unstructured grid system, we compare the results with those using a structured grid system for the two-dimensional lopsided shape. The effect of measurement errors on the estimation accuracy is also investigated.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1368-1372
• /
• 2007

바닥이 언덕과 저면으로 이루어진 연속적인 횡방향 바닥형상을 갖는 개수로 흐름의 부유사 농도분포를 수치모의 하였다. 비선형 ${\kappa}-{\varepsilon}$ 모형을 이용하여 흐름의 지배방정식인 곡선 직교좌표계에 대한 Navier-Stokes 방정식을 해석하였으며, 와점성계수 개념을 이용하여 부유사 수송 방정식을 해석하였다. 기존의 실험결과와 비교하여 모형이 격자형 이차흐름을 비교적 잘 예측하는 것을 확인하였으며, 동일한 흐름에 대하여 부유사 농도 분포를 계산하였다. 부유사량의 계산 결과 언덕 위의 부유사 농도분포가 저면 위의 농도분포보다 크게 나타났다.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.419-428
• /
• 2014

This paper presents a numerical simulation method for the flow around advancing ships in regular waves by using a rectilinear grid system. Because the grid lines do not consist with body surface in the rectilinear grid system, the body geometries are defined by the interaction points of those grid lines and the body surface. For the satisfaction of body boundary conditions, no-slip and divergence free conditions are imposed on the body surface and body boundary cells, respectively. Meanwhile, free surface is defined with the modified marker density method. The pressure on the free surface is determined to make the pressure gradient terms of the governing equations continuous, and the velocity around the free surface is calculated with the pressure on the free surface. To validate the present numerical method, a vortex induced vibration (VIV) phenomenon and flows around an advancing Wigley III ship model in various regular waves are simulated, and the results are compared with existing and corresponding research data. Also, to check the applicability to practical ship model, flows around KRISO Container Ship (KCS) model advancing in calm water are numerically simulated. On the simulations, the trim and the sinkage are set free to compare the running attitude with some other experimental data. Moreover, flows around the KCS model in regular waves are also simulated.