• 한국정밀공학회지
• /
• 제11권1호
• /
• pp.157-165
• /
• 1994

A general-purpose program for automatic two-dimensional finite-element modelling with quadrilateral elements was developed in this research. The conventional looping method employed in the program was introduced with emphasis on a new splitting criterion and a splitting scheme developed for improving the method. Some application examples were given, which show versatility and applicability of the developed program.

• 한국전산유체공학회지
• /
• 제4권2호
• /
• pp.39-46
• /
• 1999

The FVM(Finite Volume Method) have been used mainly for the flow analyses in the piston-cylinder. The objective of the present study is to analyze numerically the piston-driven intake flows using the FEM(Finite Element Method). The FEM algorithm used in this study is 4-step time-splitting method which requires much less execution time and computer storage than the velocity-pressure integrated method and the penalty method. And the explicit Lax-Wendroff scheme is applied to nonlinear convective term in the momentum equations to prevent checkerboard pressure oscillations. Also, the ALE(arbitrary Lagrangian Eulerian) method is adopted for the moving grids. The calculated results show good agreement in comparison with those by the FVM and the experimental results by the LDA.

• Computers and Concrete
• /
• 제33권3호
• /
• pp.265-273
• /
• 2024

Predicting fracture properties requires an understanding of structural failure behaviour in relation to specimen type, dimension, and notch length. Facture properties are evaluated using various testing methods, wedge splitting test being one of them. The wedge splitting test was numerically modelled three dimensionally using the finite element method on self compacting concrete specimens with varied specimen and notch depths in the current work. The load - Crack mouth opening displacement curves and the angle of rotation with respect to notch opening till failure are used to assess the fracture properties. Furthermore, based on the simulation results, failure curve was built to forecast the fracture behaviour of self-compacting concrete. The fracture failure curve revealed that the failure was quasi-brittle in character, conforming to non-linear elastic properties for all specimen depth and notch depth combinations.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제6권1호
• /
• pp.1-15
• /
• 2002

Advection-dominated transport problems possess difficulties in the design of numerical methods for solving them. Because of the hyperbolic nature of advective transport, many characteristic numerical methods have been developed such as the classical characteristic method, the Eulerian-Lagrangian method, the transport diffusion method, the modified method of characteristics, the operator splitting method, the Eulerian-Lagrangian localized adjoint method, the characteristic mixed method, and the Eulerian-Lagrangian mixed discontinuous method. In this paper relationships among these characteristic methods are examined. In particular, we show that these sometimes diverse methods can be given a unified formulation. This paper focuses on characteristic finite element methods. Similar examination can be presented for characteristic finite difference methods.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.84-85
• /
• 2003

A parallelized FEM code based on domain decomposition method has been recently developed for a large scale computational fluid dynamics. A 4-step splitting finite element algorithm is adopted for unsteady computation of the incompressible Navier-Stokes equation, and Smagorinsky LES(Large Eddy Simulation) model is chosen for turbulent flow computation. Both METIS and MPI library are used for domain partitioning and data communication between processors respectively. Tiburon of Hyundai-motor is chosen as the computational model at $Re=7.5{\times}10^{5}$, which is based on the car height. It is confirmed that the drag under road condition is smaller than that of wind tunnel condition.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 1997년도 추계학술대회논문집
• /
• pp.17-23
• /
• 1997

In the two-dimensional Finite Element Method for forming simulation, mesh generation and remeshing process are very significant. In this paper, using the modified splitting mesh generation algorithm, we can overcome the limitation of existing techniques and acquire mesh, which has optimal mesh density. A modified splitting algorithm for automatically generating quadrilateral mesh within a complex domain is described. Unnecessary meshing process for density representation is removed. Especially, during the mesh generation with high gradient density like as shear band representation, the modified mesh density scheme, which will generate quadrilateral mesh with the minimized error, which takes effect on FEM solver, is introduced.

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

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2005년도 추계 학술대회논문집
• /
• pp.117-124
• /
• 2005

Splitting algorithms of the incompressible Navier-Stokes equations using P1P1/P2P1 finite element formulation are newly proposed. P1P1 formulation allocates velocity and pressure at the same nodes, while P2P1 formulation allocates pressure only at the vertex nodes and velocity at both the vertex and mid nodes. For comparison of the elapsed time and accuracy of the two methods, they have been applied to the well-known benchmark problems. The three cases chosen are the two-dimensional steady and unsteady flows around a fixed cylinder, decaying vortex, and impinging slot jet. It is shown that the proposed P2P1 semi-splitting method performs better than the conventional P1P1 splitting method in terms of both accuracy and computation time.

• 한국전산유체공학회지
• /
• 제5권1호
• /
• pp.1-7
• /
• 2000

In general, FDM(finite difference method) and FVM(finite volume method) are used for analyzing the fluid flow numerically. However it is difficult to apply them to problems involving complex geometries, multi-connected domains, and complex boundary conditions. On the contrary, FEM(finite element method) with coordinates transformation for the unstructured grid is effective for the complex geometries. Most of previous studies have used commercial codes such as KIVA or STAR-CD for the flow analyses in the engine cylinder, and these codes are mostly based on the FVM. In the present study, using the FEM for three-dimensional, unsteady, and incompressible Navier-Stokes equation, the velocity and pressure fields in the engine cylinder have been numerically analyzed. As a numerical algorithm, 4-step time-splitting method is used and ALE(arbitrary Lagrangian Eulerian) method is adopted for moving grids. In the Piston-Cylinder, the calculated results show good agreement in comparison with those by the FVM and the experimental results by the LDA.

• Korean Journal of Mathematics
• /
• 제22권3호
• /
• pp.429-441
• /
• 2014

In this work we study two a posteriori error estimators of hierarchical type for lowest-order mixed finite element methods. One estimator is computed by solving a global defect problem based on the splitting of the lowest-order Brezzi-Douglas-Marini space, and the other estimator is locally computable by applying the standard localization to the first estimator. We establish the reliability and efficiency of both estimators by comparing them with the standard residual estimator. In addition, it is shown that the error estimator based on the global defect problem is asymptotically exact under suitable conditions.