• Journal of the Korean Mathematical Society
• /
• v.49 no.4
• /
• pp.765-778
• /
• 2012

We construct a symmetric finite volume element (SFVE) scheme for a self-adjoint elliptic problem on tetrahedron grids and prove that our new scheme has optimal convergent order for the solution and has superconvergent order for the flux when grids are quasi-uniform and regular. The symmetry of our scheme is helpful to solve efficiently the corresponding discrete system. Numerical experiments are carried out to confirm the theoretical results.

• Journal of applied mathematics & informatics
• /
• v.29 no.5_6
• /
• pp.1205-1219
• /
• 2011

In this work, a stabilized characteristic finite volume method for the time-dependent Navier-Stokes equations is investigated based on the lowest equal-order finite element pair. The temporal differentiation and advection term are dealt with by characteristic scheme. Stability of the numerical solution is derived under some regularity assumptions. Optimal error estimates of the velocity and pressure are obtained by using the relationship between the finite volume and finite element methods.

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

• Proceedings of the Korean Nuclear Society Conference
• /
• 1995.10a
• /
• pp.480-486
• /
• 1995

Turbulent flow field in a subchannel of bare rod bundles has been numerically simulated using the control volume based finite element method. Launder & Ying model of Reynolds stress and Lam & Bremhorst low-Reynolds number model are implemented in k-$\varepsilon$ equations and momentum equations. Secondary flows are simulated using the stream function and vorticity approach. The control volume based finite element method enable to use the upwind scheme (donor cell scheme). Sensitivity of the constants in the models are studied, and proper values are found to get the close result to the measured flow distributions.

• 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 the Korean Mathematical Society
• /
• v.59 no.3
• /
• pp.519-548
• /
• 2022

In this paper, a stabilized-penalized collocated finite volume (SPCFV) scheme is developed and studied for the stationary generalized Navier-Stokes equations with mixed Dirichlet-traction boundary conditions modelling an incompressible biological fluid flow. This method is based on the lowest order approximation (piecewise constants) for both velocity and pressure unknowns. The stabilization-penalization is performed by adding discrete pressure terms to the approximate formulation. These simultaneously involve discrete jump pressures through the interior volume-boundaries and discrete pressures of volumes on the domain boundary. Stability, existence and uniqueness of discrete solutions are established. Moreover, a convergence analysis of the nonlinear solver is also provided. Numerical results from model tests are performed to demonstrate the stability, optimal convergence in the usual L² and discrete H¹ norms as well as robustness of the proposed scheme with respect to the choice of the given traction vector.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.51-56
• /
• 2000

The Finite Element Solutions Is reported on solid-liquid phase change in porous media with natural convection including freezing. The model is based on volume averaged transport equations, while phase change is assumed to occur over a small temperature range. The FEM (Finite Element Method) algorithm used in this study is 3-step time-splitting method which requires much less execution time and computer storage the velocity-pressure integrated method and the penalty method. And the explicit Lax-Wendroff scheme is applied to nonlinear convective term in the energy equation. For natural convection including melting and solidification the numerical results show reasonable agreement with FDM (Finite Difference Method) results.

• 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 the Korean Institute of Telematics and Electronics
• /
• v.26 no.1
• /
• pp.48-61
• /
• 1989

In this paper, two-dimensional numerical simulation of GaAs MESFFT with 0.7${\mu}m$ gate length is perfomed. Drift-diffusion model which consider that mobility is a function of local electric field, is used. As a discretization method, instead of FDM (finite difference method) and FEM (finite element method), the Control-Volume Formulation (CVF) is used and as a numerical scheme current hybrid scheme or upwind scheme is replaced by power-law scheme which is very approximate to exponential scheme. In the process of numerical analysis, Peclet number which represents the velocity ratio of drift and diffusion, is introduced. And using this concept a current equation which consider numerical scheme at the interface of control volume, is proposed. The I-V characteristics using the model and numerical method has a good agreement with that of previous paper by others. Therefore, it is confined that it may be useful as a simulator for GaAs MESFET. Besides I-V characteristics, the mechanism of both velocity saturation in drift-diffusion model is described from the view of velocity and electric field distribution at the bottom of the channel. In addition, the relationship between the mechanism and position of dipole and drain current, are described.

• Journal of the Korean Society for Precision Engineering
• /
• v.22 no.12 s.177
• /
• pp.100-107
• /
• 2005

The present work covers a new ray tracing scheme of an injection-molded plastic lens linked with finite element analysis fur injection molding processes. The traditional ray tracing schemes have been based on the assumption that optical property of the lens is homogeneous throughout the entire volume. However, this assumption is quite unrealistic for injection-molded plastic lenses since material properties vary at every point due to injection molding effects. In order to consider non-homogeneous property of a lens, a modified ray tracing method is proposed in connection with finite element analysis of injection molding. Through the analysis of the injection molding process, the distribution of refractive indices can be obtained. This information is then utilized in the proposed ray tracing scheme based on finite element meshes so as to take into account variation of the refractive indices. The effect of mold temperature is also investigated through finite element analysis, and the relevant optical quality is evaluated through the proposed ray tracing simulation.