• Journal of the Korean Association of Geographic Information Studies
• /
• v.12 no.3
• /
• pp.128-142
• /
• 2009

A lot of efforts have been made to resolve balancing problems between flux and source term and most finite volume models under development have been applied to simple channels such as imaginary and experimental channels. However, a number of numerical problems which can not be found in simple channels occur in the application to natural rivers and the problems should be resolved to apply finite volume models to natural rivers. In this study, 2D finite volume model which is applicable to natural rivers was developed and the accuracy of the developed model was validated through the application of partial dam break In addition, a simple and efficient 2D mesh generation method was suggested and the method can be accurately reflected to 2D mesh converted from surveyed cross sections in Han-river using GIS. The accuracy and applicability of the developed model on natural rivers were verified by performing simulation on Han-river using the generated mesh and comparing computed water elevation with measured water elevation.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.75-80
• /
• 2006

Navier-Stokes simulation of the flow in a micro viscous pump is carried out. The micro viscous pump consists of a rotating circular rotor placed in a two dimensional channel. All simulation is carried out by using a finite volume approach, at the Reynolds number of 0.5, to study the performance of the micro viscous pump. Length of channel of the pump is varied to simulate the effects of the pumping load. Numerical solutions show that the net flow of the pump is realized by two counter rotating vortices formed on both sides of the rotor. The volume flow rate of the pump is decreased as length of the channel is increased, while the static pressure difference across the rotor is increased. The static pressure difference across the rotor is observed to be inversely proportional to the volume flow rate as inertia effects are negligibly small. The efficiency of the pump is found to reach a maximum when two counter rotating vortices on both sides of the rotor becomes to merge forming an outer enveloping vortex.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1999.04a
• /
• pp.195-202
• /
• 1999

In this Paper, a so-called volume-control method for nonlinear failure analysis of reinforced concrete cylindrical structures is proposed. The pressure node which defines uniform change of pressure on finite element is added into layered shell element utilizing in-plane constitutive models of reinforced concrete and layered formulation. With the pressure node formulation, one can control the change in volume enclosed by the cylindrical structures and determine the required change in pressure. An algorith of volume-control method is employed and failure analyses for RC cylindrical structures are carried out using proposed method.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.9 no.2
• /
• pp.286-290
• /
• 2008

This paper summarizes the components of an explicit aeroelastic solver developed especially for the simulation of dynamic fracture events occurring during the flight of solid propellant rockets. The numerical method combines an explicit Arbitrary Lagrangian Eulerian (ALE) version of the Cohesive Volumetric Finite Element (CVFE) scheme, used to simulate the spontaneous motion of one or more cracks propagating dynamically through a domain with regressing boundaries, and an explicit unstructured finite volume Euler code to follow the flow field during the failure event. A key feature of the algorithm is the ability to adaptively repair and expand the fluid mesh to handle the large geometrical changes associated with grain deformation and crack motion.

• 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 B
• /
• v.28 no.4
• /
• pp.395-405
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among four filling patterns at each triangular control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. Using the proposed numerical technique, the collapse of a water dam has been analyzed to predict flow phenomenon of fluid and the predicted front positions with respect to time have been compared with the reported experimental results.

• Journal of Mechanical Science and Technology
• /
• v.18 no.6
• /
• pp.1026-1035
• /
• 2004

An immersed boundary method for solving the Navier-Stokes and thermal energy equations is developed to compute the heat transfer over or inside the complex geometries in the Cartesian or cylindrical coordinates by introducing the momentum forcing, mass source/sink, and heat source/sink. The present method is based on the finite volume approach on a staggered mesh together with a fractional step method. The method of applying the momentum forcing and mass source/sink to satisfy the no-slip condition on the body surface is explained in detail in Kim, Kim and Choi (2001, Journal of Computational Physics). In this paper, the heat source/sink is introduced on the body surface or inside the body to satisfy the iso-thermal or iso-heat-flux condition on the immersed boundary. The present method is applied to three different problems : forced convection around a circular cylinder, mixed convection around a pair of circular cylinders, and forced convection around a main cylinder with a secondary small cylinder. The results show good agreements with those obtained by previous experiments and numerical simulations, verifying the accuracy of the present method.

• Journal of the Korean Society for Advanced Composite Structures
• /
• v.4 no.4
• /
• pp.15-21
• /
• 2013

This study carried out finite element deflection analysis of cylindrical shell structures made of composite materials, which is based on the micro-mechanical approach for different fiber-volume fractions. The finite element (FE) models for composite structures using multi-scale approaches described in this paper is attractive not only because it shows excellent accuracy in analysis but also it shows the effect of the material combination. New results reported in this paper are focused on the significant effects of the fiber-volume fraction for various parameters, such as fiber angles, layup sequences, and length-thickness ratios. It may be concluded from this study that the combination effect of fiber and matrix, largely governing the dynamic characteristics of composite shell structures, should not be neglected and thus the optimal combination could be used to design such civil structures for better dynamic performance.

• Journal of Mechanical Science and Technology
• /
• v.19 no.9
• /
• pp.1773-1780
• /
• 2005

The process of mould design in the foundry industry has been based on the intuition and experience of foundry engineers and designers. To bring the industry to a more scientific basis the design process should be integrated with scientific analysis such as heat transfer. The production by foundry techniques is influenced by the geometry configuration, which affects the solidification conditions and subsequent cooling. Numerical simulation and/or experiments make possible the selection of adequate materials, reducing cycle times and minimizing production costs. The main propose of this work is to study the heat transfer phenomena in the mould considering the phase change of the cast-part. Due to complex geometry of the mould, a block unstructured grid and a generalized curvilinear formulation engaged with the finite volume method is described and applied. Two types of boundary conditions, diffusive and Newtonian, are used and compared. The developed numerical code is tested in real case and the main results are compared with experimental data. The results showed that the solidification time is about 6 seconds for diffusive boundary conditions and 14.8 seconds for Newtonian boundary conditions. The use of the block unstructured grid in combination with a generalized curvilinear formulation works well with the finite volume method and allows the development of more efficient algorithms with better capacity to describe the part contours through a lesser number of elements.

• Journal of ILASS-Korea
• /
• v.5 no.1
• /
• pp.13-22
• /
• 2000

The present study has numerically analyzed the vaporization characteristics of fuel droplets in the high temperature convective flow field. The axisymmetric governing equations for mass, momentum, energy, and species are solved by an iterative and implicite unstructured finite-volume method. The moving boundary due to vaporization is handled by the deformable unstructured grid technique. The pressure-velocity coupling in the density-variable flows is treated by the SIMPLEC algorithm. In terms of the matrix solver, Bi-CGSTAB is employed for the numerically efficient and stable convergence. The n-decane is used as a liquid fuel and the initial droplet temperature is 300K. Computations are performed for the nonevaporating and evaporating droplets with the relative interphase velocity(25m/s). The unsteady vaporization process has been simulated up to the nondimensional time, 25. Numerical results indicate that the mathematical model developed in this study succesfully simulates the main features of the droplet vaporization process in the convective environment.