• Journal of Ocean Engineering and Technology
• /
• v.11 no.4
• /
• pp.141-151
• /
• 1997

This paper deals with a numerical study of flow-guider applied to controlling current in a bay. Two dimensional numerical model for tidal currents based on the depth averaged equation is developed and standard k-.epsilon. model is adopted to determine the turbulence diffusion. Equations are described in a generalized coordinate system to be implemented by non-staggered grid system and discretized by using finite volume method. Unsteady flow is simulated by fully implicit scheme. Hybrid scheme and central differencing are used to compute the convective terms and source terms, respectively. The tidal current in a rectangular bay is simulated and it gives satisfactory results. The realistic and distinct models of a large structure placed in bay are also exemplified with or without flow-guiders. The simulation results show that the flow-guider gives the residual tidal current in the bay by the different flux with respect to the direction of tidal current.

• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.36-42
• /
• 2006

The existing numerical approximations of convection flux, especially the spatial higher-order difference schemes, in unstructured cell-centered finite volume methods are examined in detail with each other and evaluated with respect to the accuracy through their application to a 2-D benchmark problem. Six higher-order schemes are examined, which include two second-order upwind schemes, two central difference schemes and two hybrid schemes. It is found that the 2nd-order upwind scheme by Mathur and Murthy(1997) and the central difference scheme by Demirdzic and Muzaferija(1995) have more accurate prediction performance than the other higher-order schemes used in unstructured cell-centered finite volume methods.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.10a
• /
• pp.47-52
• /
• 2002

We present an improved Lagrangian vortex method in 2-D incompressible unsteady viscous flows, which is based on a mesh-free integral approach of the velocity-vorticity formulation. Vorticity fields are represented by discrete vortex blobs that are updated by the Lagrangian vorticity transport with the particle strength exchange scheme. Velocity fields are expressed in a form of the Helmholtz decomposition, which are calculated by a fast algorithm of the Biot-Savart integration with a smoothed kernel and by a well-established panel method. No-slip condition is enforced through viscous diffusion of vorticity from a solid body into field. The vorticity flux is determined in such a way that spurious slip velocity vanishes. Through the comparison with the existing finite volume scheme for the transient vortical flows around an impulsively started cylinder at Reynolds number Re=550, we would obtain a more accurate scheme for vortex methods in complicated flows.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2005.05a
• /
• pp.203-208
• /
• 2005

In this paper, an automated adaptive mesh generation scheme, based on an advancing-front-Delaunay method, is developed for finite element simulation of three dimensional bulk metal forming processes. During the simulation, the finite element mesh system is adaptively remeshed whenever the mesh is unacceptable. Several schemes are developed such as curvature compensation scheme to minimize volume loss, optimal smoothing scheme to improve element quality, etc. The presented approach is evaluated and applied to automatic forging simulation in order to demonstrate the effect of the developed schemes.

• Journal of Korea Water Resources Association
• /
• v.44 no.3
• /
• pp.249-262
• /
• 2011

In this study, in order to reduce the numerical oscillation due to the unbalance between source and flux terms as the HLLC scheme is applied to the flow analysis on the irregular bed topography, a unstructured finite volume model based on the well-balanced HLLC scheme and the shallow water equations is developed and applied to problems of dam-break waves. The well-balanced HLLC scheme considers directly the gradient of bed topography as the flux terms is calculated. This scheme provides the good numerical balance between the source and flux terms in the case of the application to the steady-state transcritical flow. To verify the numerical model developed in this study, it is applied to three cases of hydraulic model experiments and a field case study of Mapasset dam failure (France). As a result of the verification, the predicted numerical results agree relatively well with available laboratory and field measurements. The model provides slightly more accurate results compared with the existing models.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.8
• /
• pp.2681-2692
• /
• 1996

The finite-volume method for radiation in a three-dimensional non-orthogonal gas turbine combustion chamber with absorbing, emitting and anisotropically scattering medium is presented. The governing radiative transfer equation and its discretization equation using the step scheme are examined, while geometric relations which transform the Cartesian coordinate to a general body-fitted coordinate are provided to close the finite-volume formulation. The scattering phase function is modeled by a Legendre polynomial series. After a benchmark solution for three-dimensional rectangular combustor is obtained to validate the present formulation, a problem in three-dimensional non-orthogonal gas turbine combustor is investigated by changing such parameters as scattering albedo, scattering phase function and optical thickness. Heat flux in case of isotropic scattering is the same as that of non-scattering with specified heat generation in the medium. Forward scattering is found to produce higher radiative heat flux at hot and cold wall than backward scattering and optical thickness is also shown to play an important role in the problem. Results show that finite-volume method for radiation works well in orthogonal and non-orthogonal systems.

• 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 Ocean Engineering and Technology
• /
• v.14 no.1
• /
• pp.11-16
• /
• 2000

수중날개를 대상으로 임의로 변형하는 이층류의 수직계산을 수행하였다. 계산은 비구조격자와 접면포착법을 이용하여 자유표면의 형상을 결정하였다 본 연구에 의하면 낮은 물속 깊이와 높은 프루드수에서 쇄파현상이 발생하였다 수중날개에 대하여 쇄파형상의 효과를 연구하였고 실험값과 계산값을 상호비교함으로서 격자의 민감성을 수직적으로 확인하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.201-210
• /
• 2009

A volume capturing method using unstructured grid system for numerical analysis of multiphase flows is introduced in the present paper. This method uses an interface capturing method (CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The novelty of CICSAM lies in the adaptive combination of high resolution discretization scheme which ensures the preservation of the sharpness and shape of the interface while retaining boundedness of the field, and no explicit interface reconstruction which is perceived to be difficult to implement on unstructured grid system. Several typical test cases for multiphase flows are presented, which are simulated by an in-house solution code(PowerCFD). This code employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with CICSAM. It is found that the present method simulates efficiently and accurately complex free surface flows such as multiphase flows.

• Journal of Mechanical Science and Technology
• /
• v.14 no.9
• /
• pp.985-996
• /
• 2000

The issue of quench is related to safety operation of large-scale superconducting magnet system fabricated by cable-in-conduit conductor. A numerical method is presented to simulate the thermal hydraulic quench characteristics in the superconducting Tokamak magnet system, One-dimensional fluid dynamic equations for supercritical helium and the equation of heat conduction for the conduit are used to describe the thermal hydraulic characteristics in the cable-in-conduit conductor. The high heat transfer approximation between supercritical helium and superconducting strands is taken into account due to strong heating induced flow of supercritical helium. The fully implicit time integration of upwind scheme for finite volume method is utilized to discretize the equations on the staggered mesh. The scheme of a new adaptive mesh is proposed for the moving boundary problem and the time term is discretized by the-implicit scheme. It remarkably reduces the CPU time by local linearization of coefficient and the compressible storage of the large sparse matrix of discretized equations. The discretized equations are solved by the IMSL. The numerical implement is discussed in detail. The validation of this method is demonstrated by comparison of the numerical results with those of the SARUMAN and the QUENCHER and experimental measurements.