• Water Engineering Research
• /
• v.2 no.3
• /
• pp.179-185
• /
• 2001

The time consuming and expensive nature of experimental research on scouring processes caused by flowing water makes it attractive to develop numerical tools for the predication of the interaction of the fluid flow and the movable bed. In this paper the numerical simulation of scour by a wall jet is presented. The flow is assumed to be two-dimensional, and the alluvium is cohesionless. The solution process, repeated at each time step, involves simulation of a turbulent wall jet flow, solution of the convection-diffusion of sand concentration, and prediction of the bed deformation. For simulation of the jet flow, the governing equations for momentum, mass balance and turbulent parameters are solved by the finite volume method. The SIMPLE scheme with momentum interpolation is used for pressure correction. The convection-diffusion equation is solved for sediment concentration. A boundary condition for concentration at the bed, which takes into account the effect of bed-load, is implemented. The time rate of deposition and scour at the bed is obtained by solving the continuity equation for sediment. The shape and position of the scour hole and deposition of the bed material downstream of the hole appear realistic.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2005.10a
• /
• pp.269-272
• /
• 2005

For the three decades, the mold-filling of injection molding process was modeled as Hele-Shaw model. However, this model can not consider the 3D effect. In this paper, numerical simulations of three dimensional mold-filling during the filling phase were performed. The governing equations were discretized by segregated finite element method, which used equal order interpolation for pressure and velocity fields. The iterative linear equation solver (JCG, SOR) was employed for the solution of the momentum and pressure equations. Volume of Fluid (VOF) was employed for the melt front advancement. To check the validity of the numerical results, the results were compared with the experimental ones. The agreements between the experiment and the numerical results were found to be satisfactory.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05b
• /
• pp.173-178
• /
• 1996

A three-dimensional thermo-hydraulic computer code is developed for simulation of incompressible flows in complex geometries. The computer code employs a body-fitted, nonorthogonal grid system in order to efficiently handle the complex geometries encountered in many engineering applications. The finite volume method is used to discretize the governing equations and the convection term is treated by higher-order bounded schemes. The cell-centered, nonstaggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is avoided by use of momentum interpolation practice. The computer code employs the SIMPLE algorithm for pressure and velocity coupling and the k-$\varepsilon$ turbulence for turbulent calculation. The computer code has been tested through application to a variety of test problems and some results are presented in this paper

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.474-481
• /
• 2001

This paper addresses a numerical method for predicting transient temperature distributions in the wall of a curved pipe subjected to internally thermal stratification flow. A simple and convenient numerical method of treating the unsteady conjugate heat transfer in the non-orthogonal coordinate systems is presented. The proposed method is implemented in a finite volume thermal-hydraulic computer code based on a cell-centered, non-staggered grid arrangement, the SIMPLEC algorithm, a higher-order bounded convection scheme, and the modified version of momentum interpolation method. Calculations are performed for the transient evolution of thermal stratification in two curved pipes, where the one has thick wall and the other has so thin wall that its presence can be negligible in the heat transfer analysis. The predicted results show that the thermally stratified flow and transient conjugate heat transfer in a curved pipe with a finite wall thickness can be satisfactorily analyzed by the present numerical method, and that the neglect of wall thickness in the prediction of pipe wall temperature distributions can provide unacceptably distorted results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.9 s.240
• /
• pp.1049-1056
• /
• 2005

A conservative pressure-based finite-volume numerical method has been developed for computing flow and heat transfer by using an unstructured grid system. The method admits arbitrary convex polyhedra. Care is taken in the discretization and solution procedures to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are found by a novel second-order accurate spatial discretization. Momentum interpolation is used to prevent pressure checkerboarding and the SIMPLE algorithm is used for pressure-velocity coupling. The resulting set of coupled nonlinear algebraic equations is solved by employing a segregated approach, leading to a decoupled set of linear algebraic equations fer each dependent variable, with a sparse diagonally dominant coefficient matrix. These equations are solved by an iterative preconditioned conjugate gradient solver which retains the sparsity of the coefficient matrix, thus achieving a very efficient use of computer resources.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.900-907
• /
• 2000

The three-dimensional turbulent flow in curved pipes susceptible to flow-accelerated corrosion has been analyzed numerically to predict the pressure and shear stress distributions on the inner surface of the pipes. The analysis employs the body-fitted non-orthogonal curvilinear coordinate system and a standard $ {\kappa}-{\varepsilon}$ turbulence model with wall function method. The finite volume method is used to discretize the governing equations. The convection term is approximated by a high-resolution and bounded discretization scheme. The cell-centered, non-staggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is prevented by the application of a modified version of momentum interpolation scheme. The SIMPLE algorithm is employed for the pressure and velocity coupling. The numerical calculations have been performed for two curved pipes with different bend angles and curvature radii, and discussions have been made on the distributions of the primary and secondary flow velocities, pressure and shear stress on the inner surface of the pipe to examine applicability of the present analysis method. As the result it is seen that the method is effective to predict the susceptible systems or their local areas where the fluid velocity or local turbulence is so high that the structural integrity can be threatened by wall thinning degradation due to flow-accelerated corrosion.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.186-191
• /
• 2004

The reactive flowfield of the transverse injecting combustor has been studied using Euler-Lagrange method in order to develop an efficient solution procedure for the understanding of liquid spray combustion in the transverse injecting combustor which has been widely used in ramjets and turbojet afterburners. The unsteady two-dimensional gas-phase equations have been represented in Eulerian coordinates and the liquid-phase equations have been formulated in Lagrangian coordinates. The gas-phase equations based on the conservation of mass, momentum, and energy have been supplemented by combustion. The vaporization model takes into account the transient effects associated with the droplet heating and the liquid-phase internal circulation. The droplet trajectories have been determined by the integration of the Lagrangian equation in the flow field obtained from the separate calculation without considering the iterative effect between liquid and gas phases. The reported droplet trajectories had been found to deviate from the initial conical path toward the flow direction in the very end of its lifetime when the droplet size had become small due to evaporation. The integration scheme has been based on the TEACH algorithm for gas-phase equation, the second order Runge-Kutta method for liquid-phase equations and the linear interpolation between the two coordinate systems. The calculation results has shown that the characteristics of the droplet penetration and recirculation have been strongly influenced by the interaction between gas and liquid phases in such a way that most of the vaporization process has been confined to the wake region of the injector, thereby improving the flame stabilization properties of the flowfield.

• Journal of Mechanical Science and Technology
• /
• v.20 no.12
• /
• pp.2218-2229
• /
• 2006

A conservative finite-volume numerical method for unstructured grids with the cell-centered method has been developed for computing flow and heat transfer by combining the attractive features of the existing pressure-based procedures with the advances made in unstructured grid techniques. This method uses an integral form of governing equations for arbitrary convex polyhedra. Care is taken in the discretization and solution procedure to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. For both convective and diffusive fluxes the forms superior to both accuracy and stability are particularly adopted and formulated through a systematic study on the existing approximation ones. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are computed by using a linear reconstruction based on the divergence theorem. Momentum interpolation is used to prevent the pressure checkerboarding and a segregated solution strategy is adopted to minimize the storage requirements with the pressure-velocity coupling by the SIMPLE algorithm. An algebraic solver using iterative preconditioned conjugate gradient method is used for the solution of linearized equations. The flow analysis code (PowerCFD) developed by the present method is evaluated for its application to several 2-D structured-mesh benchmark problems using a variety of unstructured quadrilateral and triangular meshes. The present flow analysis code by using unstructured grids with the cell-centered method clearly demonstrate the same accuracy and robustness as that for a typical structured mesh.

• Korean Journal of Applied Biomechanics
• /
• v.15 no.3
• /
• pp.117-132
• /
• 2005

It was to study a following research of "A Kinematic Analysis of Air-rolling-breakfall in Judo". The purpose of this study was to analyze the Center of Gravity(COG) variables when performing Air-rolling-breakfall motion, while passing forward over(PFO) to the vertical-hurdles(2m height, take off board 1m height) in judo. Subjects were four males of Y. University squad, who were trainees of the demonstration exhibition team, representatives of national level judoists and were filmed by four 5-VHS 16mm video cameras(60field/sec.) through the three dimensional film analysis methods.COG variable were anterior-posterior directional COG and linear velocity of COG, vertical directional COG and linear velocity of COG. The data collections of this study were digitized by KWON3D program computed The data were standardized using cubic spline interpolation based by calculating the mean values and the standard deviation calculated for each variables. When performing the Air-rolling-breakfall, from the data analysis and discussions, the conclusions were as follows : 1. Anterior-posterior directional COG(APD-COG) when performing Air-rolling-breakfall motion, while PFO over to the vertical-hurdles(2m height) in judo. The range of APD-COG by forward was $0.31{\sim}0.41m$ in take-off position(event 1), $1.20{\sim}1.33m$ in the air-top position(event 2), $2.12{\sim}2.30m$ in the touch-down position(event 3), gradually and $2.14{\sim}2.32m$ in safety breakfall position(event 4), respectively. 2 The linear velocity of APD-COG was $1.03{\sim}2.14m/sec$. in take-off position(event 1), $1.97{\sim}2.22m/sec$. gradually in the air-top position(event 2), $1.05{\sim}1.32m/sec$. in the touch-down position (event 3), gradual decrease and $0.91{\sim}1.23m/sec$. in the safety breakfall position(event 4), respectively. 3. The vertical directional COG(VD-COG) when performing Air-rolling-breakfall motion, while PFO to the vertical-hurdles(2m height) in judo. The range of VD-COG toward upward from mat was $1.35{\sim}1.46m$ in take-off position(event 1), the highest $2.07{\sim}2.23m$ in the air-top position(event 2), and after rapid decrease $0.3{\sim}0.58m$ in the touch-down position(event 3), gradual decrease $0.22{\sim}0.50m$ in safety breakfall position(event 4), respectively. 4. The linear velocity of VlJ.COG was $1.60{\sim}1.87m/sec$. in take-off position(event 1), $0.03{\sim}0.08m/sec$. gradually in the air-top position(event 2), $-4.37{\sim}\;-4.76m/sec$. gradual decrease in the touch-down position(event 3), gradual decrease and -4.40${\sim}\;-4.77m/sec$. in safety breakfall position(event 4), respectively. When performing Air-rolling-breakfall showed parabolic movement from take-off position to air-top position, and after showed vertical fall movement from air-top position to safety breakfall. In conclusion, Ukemi(breakfall) is safety fall method Therefore, actions need for performing safety fall movement, that decrease and minimize shock and impact during Air-rolling-breakfall from take-off board action to air-top position must be maximize of angular momentum, and after must be minimize in touch-down position and safety breakfall position.