• 설비공학논문집
• /
• 제11권6호
• /
• pp.766-773
• /
• 1999

A numerical analysis has been performed solidification problem using the fixed grid-enthalpy method with enthalpy-porosity relation. A modified standard $k-\varepsilon$ model was applied to describe the influence of turbulent flow. Computational procedures are based on the finite volume method and the non-staggered grid system. Comparisons with the different three experimental results show that applying a modified standard $k-\varepsilon$model in mushyzone is better than the previous computation results. This paper includes another EMBR's influences such as change of velocity field, Increasement of temperature and dispersion of flow out of nozzle into the flow field. These EMBR's influences are compared to case without EMBR.

• 설비공학논문집
• /
• 제11권6호
• /
• pp.733-739
• /
• 1999

A modified low-Re $k-\varepsilon$ model is used for the calculation of drag-reducing turbulent flow by polymer injection in a pipe. With the viscoelastic model, molecular viscosity in the definition of turbulent viscosity is related to elongations viscosity of the solution to account for the effects of drag reduction. Finite volume method is used for the discretization, and power-law scheme is used as a numerical scheme. Computed dimensionless velocity profiles are in good agreements with the experimental data in case of low drag reductions. However, in case of high drag reductions, they deviate largely from the measurements in the central zone of the flow field.

• Ocean Systems Engineering
• /
• 제12권1호
• /
• pp.63-86
• /
• 2022

A refined projection-based purely Lagrangian meshfree method is presented towards reliable numerical analysis of fluid flow interactions with saturated/unsaturated porous media of uniform/spatially-varying porosities. The governing equations are reformulated on the basis of two-phase mixture theory with incorporation of volume fraction. These principal equations of mixture are discretized in the context of Incompressible SPH (Smoothed Particle Hydrodynamics) method. Associated with the consideration of governing equations of mixture, a new term arises in the source term of PPE (Poisson Pressure Equation), resulting in modified source term. The linear and nonlinear force terms are included in momentum equation to represent the resistance from porous media. Volume increase of fluid particles are taken into consideration on account of the presence of porous media, and hence multi-resolution ISPH framework is also incorporated. The stability and accuracy of the proposed method are thoroughly examined by reproducing several numerical examples including the interactions between fluid flow and saturated/unsaturated porous media of uniform/spatially-varying porosities. The method shows continuous pressure field, smooth variations of particle volumes and regular distributions of particles at the interface between fluid and porous media.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• 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.

• 대한기계학회논문집
• /
• 제19권8호
• /
• pp.1950-1963
• /
• 1995

A modification of the approximate-factorization method is made to accelerate the convergency rate and to take sufficiently large Courant number without loss of accuracy. And a stable implicit finite-difference scheme for solving the incompressible Navier-Stokes equations employed above modified method is developed. In the present implicit scheme, the volume fluxes with contravariant velocity components and the pressure formulation in curvilinear coordinates is adopted. In order to satisfy the continuity condition completely and to remove spurious errors for the pressure, the Navier-Stokes equations are solved by a modified SMAC scheme using a staggered gird. The upstream-difference scheme such as the QUICK scheme is also employed to the right hand side. The implicit scheme is unconditionally stable and satisfies a diagonally dominant condition for scalar diagonal linear systems of implicit operator on the left hand side. Numerical results for some test calculations of the two-dimensional flow in a square cavity and over a backward-facing step are obtained using both usual approximate-factorization method and the modified one, and compared with each other. It is shown that the present scheme allows a sufficiently large Courant number of O(10$^{2}$) and reduces the computing time.

• 대한기계학회논문집B
• /
• 제21권9호
• /
• pp.1149-1164
• /
• 1997

In this study, the applicability of the RNG k-.epsilon. model to the analysis of the complex flows is studied. The governing equations based on a non-orthogonal coordinate formulation with Cartesian velocity components are used and discretized by the finite volume method with non-staggered variable arrangements. The predicted results using the RNG k-.epsilon. model of three complex flows, i.e., the flow over a backward-facing step and a blunt flat plate, the flow around a 2D model car are compared to these from the standard k-.epsilon. model and experimental data. That of the unsteady axisymmetric turbulent flow within a cylinder of reciprocating model engine including port/valve assembly and the spray characteristics within a chamber of direct injection model engine are compared to these from the standard k-.epsilon. model and experimental data. The results of reattachment length, separated eddy size, average surface pressure distribution using the RNG k-.epsilon. model show more reasonable trends comparing with the experimental data than those using the modified k-.epsilon. model. Although the predicted rms velocity using the modified k-.epsilon. model is lower considerably than the experimental data in incylinder flow with poppet valve, predicted axial and radial velocity distributions at the valve exit and in-cylinder region show good agreements with the experimental data. The spray tip penetration predicted using the RNG k-.epsilon. model is more close to the experimental data than that using the modified k-.epsilon. model. The application of the RNG k-.epsilon. model seems to have some potential for the simulations of the unsteady turbulent flow within a port/valve-cylinder assembly and the spray characteristics over the modified k-.epsilon. model.

• Journal of Hydrospace Technology
• /
• 제2권2호
• /
• pp.57-67
• /
• 1996

The formulation for hybrid-QUICK scheme of convective transport terms in finite-volume calculation procedure is presented. Source terms are modified to apply the hybrid-QUICK scheme. Test calculations are performed for wall-driven cavity flow at Re=$10_2$, $10_3$, and $10_4$. These include the evaluation of boundary conditions approximated by third-order finite difference scheme. The stable and converged solutions are obtained without unsteady terms in the momentum equations. The results using hybrid-QUICK scheme show no difference with those using hybrid scheme at low Re ($=10_2$) and are better at higher Re ($10_3$, and $10_4$).

• 대한기계학회논문집
• /
• 제15권2호
• /
• pp.668-676
• /
• 1991

Turbulent flows around tube banks and in the diffuser were studied using a non-orthogonal boundary fitted coordinate system and the modified K-.epsilon. turbulence model. In these cases, many problems emerge which stem from the geometrical complexity of the flow domain and the physical complexity of turbulent flow itself. To treat the complex geometry, governing equations were reformulated in a non-orthogonal coordinate system with Cartesian velocity components and discretised by the finite volume method with a non-staggered variable arrangement. The modified K-.epsilon. model of Hanjalic and Launer was applied to solve above two cases under the condition of strong and mild pressure gradient. The results using the modified K-.epsilon. model results in both test cases.

• Journal of Advanced Marine Engineering and Technology
• /
• 제23권5호
• /
• pp.662-670
• /
• 1999

A numerical prediction was performed to clarify the air motion in the cylinder of an axisymmet-ric four-stroke reciprocating engine at its intake and compression stage. A scheme of finite volume method is used for the calculation. Modified $k-{\varepsilon}$ turbulence model is adopted and wall function is applied to the grids near the wall. The predicted mean velocity and rms velocity profiles showed a reasonable agreement with an available experimental data at its intake and compression stage. The predicted in-cylinder flow fields show that a strong turbulent twin vortex structure is pro-duced during induction but it commences to decay rapidly around inlet valve closure. The mean velocity continues to fall to a low level during compression but the turbulence intensity attains an approximate constant level.

• 대한조선학회지
• /
• 제24권3호
• /
• pp.9-16
• /
• 1987

In this paper, the fluid flow in the two-dimensional tank is analyzed by the Finite Difference Method. The Navier-Stokes equation is modified for the tank fixed coordinate system. For the treatment of the free surface, the Volume of Fluid Method by Hirt and Nichols is adopted. The continuity equation and the Poisson equation which is derived from the Navier-Stokes equation to find the pressure are solved by the Successive-Line-Overrelaxation Method. The comparison of the calculated results with experimental data show a favorable agreement. The fluid flow in the two-dimensional tank can be predicted reasonably before the free surface reaches breaking by this numerical method.