• 한국전산유체공학회:학술대회논문집
• /
• 2000.10a
• /
• pp.135-140
• /
• 2000

Numerical calculations are carried out in order to evaluate the performance of low-Re Reynolds stress model based on SSG model for a swirling turbulent flow in a pipe. The results are compared with those of $\kappa-\epsilon$ model and GL model, and the experimental data. The finite volume method is used for the discretization, and the power-law scheme is employed as a numerical scheme. The SIMPLE algorithm is used for velocity-Pressure correction in the governing equations.

• Journal of the Semiconductor & Display Technology
• /
• v.3 no.3
• /
• pp.19-21
• /
• 2004

This paper presents a 3-dimensional diffusion simulation with adaptive solution strategy. The developed diffusion simulator VLSIDIF-3 was designed to re-refine areas. Refine scheme was calculated by the difference of doping concentration between any of two nodes. Each element is greater than tolerance and redo diffusion process until error is tolerable. Numerical experiment in low doping diffusion problem showed that this adaptive solution strategy is very efficient in both memory and time, and expected this scheme would be more powerful in complex diffusion model.

• Journal of Korea Water Resources Association
• /
• v.36 no.1
• /
• pp.1-11
• /
• 2003

This paper describes a finite volume, two-dimensional model. It adopts a recently developed essentially non-oscillatory(ENO) schemes based on the Lax-Friedrichs solver, which was modified for a finite volume grid, and employs a modified MUSCL(Monotonic Upstream centered Scheme for Conservation Law) for second-order accuracy in space. To demonstrate the applications of the model, it is applied to solve the 1-D and 2-D dam-break problems. The model in conjunction with the modified MUSCL showed a better agreement with analytical solutions than the minmod function in 1-D dam-break problems and is satisfactorily validated with documented published data in 2-D dam-break problems. The model was applied to tidal wane entering channel at one end, and the results showed a good agreement with analytical solutions. In the channel with reflective boundary conditions specified at the extremities, the model was capable of accurately simulating the wave propagation.

• Journal of electromagnetic engineering and science
• /
• v.8 no.3
• /
• pp.100-109
• /
• 2008

The finite volume time domain(FVTD) technique faces serious limitations in simulating electromagnetic scattering at high frequencies due to requirements related to discretization. A modified FVTD method is proposed for electrically large, perfectly conducting scatterers by partially incorporating a time-domain physical optics(PO) approximation for the surface current. Dominant specular returns in the modified FVTD method are modeled using a PO approximation of the surface current allowing for a much coarser discretization at high electrical sizes compared to the original FVTD scheme. This coarse discretization can be based on the minimum surface resolution required for a satisfactory numerical evaluation of the PO integral for the scattered far-field. Non-uniform discretization and spatial accuracy can also be used in the context of the modified FVTD method. The modified FVTD method is aimed at simulating electromagnetic scattering from geometries containing long smooth illuminated sections with respect to the incident wave. The computational efficiency of the modified FVTD method for higher electrical sizes are shown by solving two-dimensional test cases involving electromagnetic scattering from a circular cylinder and a symmetric airfoil.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.9 s.180
• /
• pp.2397-2405
• /
• 2000

The remeshing is a method to replace a distorted mesh by a new mesh without interrupting the finite element calculation. The remeshing procedure in this paper refers to the rezoning, for which a sm oothing process is developed to alleviate the distortions of mesh. In the paper, an automatic finite element rezoning system with tetrahedral elements for large deformation analysis has been developed. Our smoothing process is composed of two steps, a surface smoothing and a volume smoothing. In the surface smoothing, checking the dihedral angle and projection on surface patch reduced the change of shape and nodes penetrating die. The constrained Laplacian smoothing has been employed for the volume smoothing process. The state variables are mapped from old mesh to new mesh by using volume coordinates within a tetrahedral element. All these procedures have been linked to the NIKE3D program As illustrated in the examples the overall strategy ensures a robust and efficient rezoning scheme for finite element simulation of metal-forming processes.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.4
• /
• pp.19-30
• /
• 2000

A finite-volume method is developed to solve turbulent flows around modern commercial hull forms with bow and stern bulbs. The RANS equations are solved. The cell-centered finite-volume method employs QUICK and central difference scheme for convective and diffusive flux discretization, respectively. The SIMPLEC method is adopted for the velocity-pressure coupling. The developed numerical methods are applied to calculate turbulent flow around KRISO 3600TEU container ship. Surface meshes are generated into five blocks: bow and stern bulbs, overhang, fore and afterbody. 3-D field grid system with O-H topology is generated using elliptic grid generation method. Surface friction lines and wake distribution at propeller plane is compared with experiment. The calculated results show that the present method can be used to predict flow around a modern commercial hull forms with bulbs.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.12
• /
• pp.1555-1569
• /
• 2000

Numerical simulation of fluid flow with moving free surface has been carried out. For the free surface flow, a VOF(Volume of Fluid)-based algorithm utilizing a fixed grid system has been investigated. In order to reduce numerical smearing at the free surface represented on a fixed grid system, a new free surface tracking algorithm based on the donor-acceptor scheme has been presented. Novel features of the proposed algorithm are characterized as two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed algorithm can be easily implemented in any irregular non-uniform grid systems that are usual in finite element method (FEM). Moreover, the proposed algorithm can be extended and applied to the 3-D free surface flow problem without additional efforts. For computation of unsteady incompressible flow, a finite element approximation based on the explicit fractional step method has been adopted. In addition, the SUPG(streamline upwind/Petrov-Galerkin) method has been implemented to deal with convection dominated flows. Combination of the proposed free surface tracking scheme and explicit fractional step formulation resulted in an efficient solution algorithm. Validity of the present solution algorithm was demonstrated from its application to the broken dam and the solitary wave propagation problems.

• Water for future
• /
• v.29 no.5
• /
• pp.173-184
• /
• 1996

In this study, a two-dimensional shallow-water equation was used to develop the mathematical model for computing water levels and flow distribution. In the discretization equations, based on the finite volume method (FVM), the third order Runge-Kutta method and the third order upwind scheme were introduced to handle the unsteady and vconvective terms in the governing equations. To determine the accuracy of the developed model, it was applied to the rectangular horizontal channel in a frictionless flow. The water depth and velocity obtained by the numerical model were found to agree closely with the exact solution. The model was also applied to the rectangular channel with both the symmetric and the non symmetric constriction. The velocity distribution of the flow and the propagation of the flood wave were simulated and the results well described the flow characteristics.

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

• Journal of the Korean Society of Hazard Mitigation
• /
• v.9 no.6
• /
• pp.89-97
• /
• 2009

In this study, a two-dimensional unstructured finite volume model based on the shallow-water equations and well-balanced HLLC scheme is developed. The model is verified by applying to various one- and two-dimensional problems related to the analyses of dam-break wave. The predicted numerical results agree very well with available analytical solutions and laboratory measurements. The model provides slightly more accurate results compared with the existing models.