• Journal of Korea Water Resources Association
• /
• v.31 no.3
• /
• pp.279-290
• /
• 1998

The height and speed of the shock wave are critical data in flood-control operations or in the design of channel walls and bridges along rivers with high flow velocities. Therefore, a numerical model is needed for simulating flow discontinuity over a wide range of conditions. In this study, a governing equation. As a Riemann solver Roe(1981)'s one is used. The model employs the modified MUSCL for handling the unstructured grids in this research. this model that adopts the explicit tradditional twl dimmensional dam break problems, two hydraulic dam break model is simulations, and a steady state simulation in a curved channel. Conclusions of this research are as follows : 1) the finite volume method can be combined with the Godonov-type method that is useful for modeling shocks. Hence, the finite volume method is suitable for modeling shocks. 2) The finite volume model combined with the modified MUSCL is successful in modeling shock. Therefore, modified MUSCL is proved to be valid.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1230-1243
• /
• 1997

A numerical technique for simulating incompressible viscous flow with free surface is presented. The flow field is obtained by penalty finite element formulation. In this work, a modified volume of fluid (VOF) method which is compatible with 4-node element is proposed to track the moving free surface. This scheme can be applied to irregular mesh system, and can be easily extended to three dimensional geometries. Numerical analyses were done for two benchmark examples, namely the broken dam problem and the solitary wave propagation problem. The numerical results were in close agreement with the existing data. Illustrative examples were studied to show the effectiveness of the proposed numerical scheme.

• Journal of computational fluids engineering
• /
• v.14 no.1
• /
• pp.35-44
• /
• 2009

To adequately analyze flows in pipe or duct network system, traditional node-based junction coupling methods require the junction loss which is specified by empirical or analytic correlations. In this paper, a new finite volume junction coupling method using a ghost junction cell is developed by considering the interchange of linear momentum as well as the important wall-effect at junction without requiring any correlation on the junction loss. Also, boundary treatment is modified to preserve the stagnation enthalpy across boundaries, such as pipe-end and the interface between junction and branch. Also, the computational accuracy and efficiency of the Godunov-type finite volume schemes are investigated by tracing the total mechanical energy of rapid transients due to sudden closure of valve at downstream end.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1996.04a
• /
• pp.787-791
• /
• 1996

In injection molding, location of gates have great influence on the quality of plastic parts. Usually, they are located by releated trial and errors of experienced mold designers. In this topic we will present the numerical algorithm for finding the optimal gate locations. Optimization algorithm is devided into two stages. In the first stage, candidated optimal gate locations can be found by geometry of part only; whereas in the next step, more acculate gate locations are selected byiterative computation with optimization part and analysis part. So from the following study, we suggested the modified flow-volume method, which will define the optimal gate locations in injection mold design.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.10
• /
• pp.2762-2772
• /
• 1994

Branch flows for Newtonian and non-Newtonian fluids are simulated by the finite volume method. The modified power-law model is employed as a constitutive equation of the non-Newtonian fluids. Numerical analyses are focused on understanding of flow patterns for different values of branch angles, diameter ratios and Reynolds numbers. The numerical results are compared with the existing experimental data. The calculated velocity profiles and pressure variations are in good agreement with available experimental results.

• Korea-Australia Rheology Journal
• /
• v.13 no.2
• /
• pp.97-106
• /
• 2001

Three-dimensional flow analysis was performed by using the control volume finite-element method for design of a profile extrusion die. Because polymer melt behavior is complicated and cross-sectional shape of the profile extrusion die is changing continuously, the fluid flow within the die must be analyzed three-dimensionally. A commercially available polypropylene is used for theoretical and experimental investigations. Material properties are assumed to be constant except for the viscosity. The 5-constant modified Cross model is used for the numerical analysis. A test problem is examined in order to verify the accuracy of the numerical method. Simulations are performed for conditions of three different screw speeds and three different die temperatures. Predicted pressure distribution is compared with the experimental measurements and the results of the previous two-dimensional study. The computational results obtained by using three dimensional CVFEM agree with the experimental measurements and are more accurate than those obtained by using the two-dimensional cross-sectional method. The velocity profiles and the temperature distributions within several cross-sections of the die are given as contour plots.

• Journal of Advanced Marine Engineering and Technology
• /
• v.22 no.4
• /
• pp.436-444
• /
• 1998

Applicability of the RNG k-$\varepsilon$ model to the analysis of unsteady axisymmetric turbulent flow of a reciprocating engine including port/valve assembly is studied numerically. The governing equations based on non-orthogonal including port/valve assembly is studied numerically. The governing equations based on a non-orthogonal coordinate formulation with Cartesian velocity components are used and discretised by the finite volume method with non-staggered variable arrangements. The predicted results using the RNG k-$\varepsilon$ model of the unsteady axisymmetric turbulent flow within a cylinder of reciprocating model engine including port/valve assembly are compared to these from the modified k-$\varepsilon$ model and experimental data. Using the RNG k-$\varepsilon$ model seems the have some potential for the simulations of the unsteady turbulent flow within a port/valve-cylinder assembly over the modified k-$\varepsilon$model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.9
• /
• pp.1148-1156
• /
• 2000

The QUICK scheme for convection terms is modified for unstructured finite volume method by using linear reconstruction technique and validated through the computation of two well defined laminar flows. It uses two upstream grid points and one downstream grid point in approximating the convection terms. The most upstream grid point is generated by considering both the direction of flow and local grid line. Its value is calculated from surrounding grid points by using a linear construction method. Numerical error by the modified QUICK scheme is shown to decrease about 2.5 times faster than first order upwind scheme as grid size decreases. Computations are also carried out to see effects of the skewness and irregularity of grid on numerical solution. All numerical solutions show that the modified QUICK scheme is insensitive to both the skewness and irregularity of grid in terms of the accuracy of solution.

• 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.

• 한국연소학회:학술대회논문집
• /
• 1998.10a
• /
• pp.139-154
• /
• 1998

Under certain circumstance, premixed turbulent flame can be treated as wrinkled thin laminar flame and its motion in a hydrodynamic flow field has been investigated by employing G-equation. Past studies on G-equation successfully described certain aspects of laminar flame propagation such as effects of stretch on flame speed. In those studies, flames were regarded as a passive interface that does not influence the flow field. The experimental evidences, however, indicate that flow field can be significantly modified by the propagation of flames through the volume expansion of burned gas. In the present study, a new method to be used with G -equation is described to include the effect of volume expansion in the flame dynamics. The effect of volume expansion on the flow field is approximated by Biot-Savart law. The newly developed model is validated by comparison with existing analytical solutions of G -equation to predict flames propagating in hydrodynamic flow field without volume expansion. To further investigate the influence of volume expansion, present method was applied to initially wrinkled or planar flame propagating in an imposed velocity field and the average flame speed was evaluated from the ratio of flame surface area and projected area of unburned stream channel. It was observed that the initial wrinkling of flame cannot sustain itself without velocity disturbance and wrinkled structure decays into planar flame as the flame propagates. The rate of decay of the structure increased with volume expansion. The asymptotic change in the average burning speed occurs only with disturbed velocity field. Because volume expansion acts directly on the velocity field, the average burning speed is affected at all time when its effect is included. With relatively small temperature ratio of 3, the average flame speed increased 10%. The combined effect of volume expansion and flame stretch is also considered and the result implied that the effect of stretch is independent of volume release.