• Journal of applied mathematics & informatics
• /
• 제15권1_2호
• /
• pp.1-28
• /
• 2004

We first apply a first order splitting to a semilinear reaction-diffusion equation and then discretize the resulting system by an $H^1$-Galerkin mixed finite element method in space. This semidiscrete method yields a system of differential algebraic equations (DAEs) of index one. A priori error estimates for semidiscrete scheme are derived for both differ-ential as well as algebraic components. For fully discretization, an implicit Runge-Kutta (IRK) methods is applied to the temporal direction and the error estimates are discussed for both components. Finally, we conclude the paper with a numerical example.

• Interaction and multiscale mechanics
• /
• 제6권2호
• /
• pp.83-105
• /
• 2013

In this paper, a multi-scale meshfree-enriched finite element formulation is presented for the analysis of acoustic wave propagation problem. The scale splitting in this formulation is based on the Variational Multi-scale (VMS) method. While the standard finite element polynomials are used to represent the coarse scales, the approximation of fine-scale solution is defined globally using the meshfree enrichments generated from the Generalized Meshfree (GMF) approximation. The resultant fine-scale approximations satisfy the homogenous Dirichlet boundary conditions and behave as the "global residual-free" bubbles for the enrichments in the oscillatory type of Helmholtz solutions. Numerical examples in one dimension and two dimensional cases are analyzed to demonstrate the accuracy of the present formulation and comparison is made to the analytical and two finite element solutions.

• 대한기계학회논문집A
• /
• 제32권6호
• /
• pp.496-503
• /
• 2008

With the improvement of computer power and technology, fracture modelling by finite element methods has become a topic of extensive studies. However, fracture simulation much limited to an academic study of crack propagation with a fine mesh. Element deletion method is a useful tool for estimating damage due to accidental or extreme loads on structures, provided that an effective and realistic criterion is established for simulating the material failure and subsequent element deletion. In this study, ABAQUS/Explicit is used to simulate the material failure on the basis of experimental results by X. Huang et al. Through numerical experiments, we suggest a formulation to determine the failure strain associated with the size and thickness of removed elements.

• 대한기계학회논문집
• /
• 제17권12호
• /
• pp.2995-3006
• /
• 1993

An automatic mesh generation scheme has been developed for finite element analysis with two-dimensional, quadrilateral elements. The basic strategies of the method are to transform the analysis domain into loops with key nodes and the loops are recursively subdivided into subloops with the use of best split lines. Finally by using the basic loop operators, the meshes are completed. In this algorithm an eight-node loop operator is proposed, which is useful in the area where the change of element size is large and the splitting criteria for subdividing the loops have also been modified to the existing algorithms. Lines, arcs, and cubic spline curves are used to define the boundaries of analysis domain. Sample meshes for several geometries are presented to demonstrate the robustness of the algorithm.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.223-227
• /
• 2009

In the present study, a least square weighted residual method and Taylor-Galerkin method were formulated and tested for the discretization of the two hyperbolic type equations of level set method; advection and reinitialization equations. The two approaches were compared by solving a time reversed vortex flow and three-dimensional broken dam flow by employing a four-step splitting finite element method for the solution of the incompressible Navier-Stokes equations. From the numerical experiments, it was shown that the least square method is more accurate and conservative than Taylor-Galerkin method and both methods are approximately first order accurate when both advection and reinitialization phase are involved in the evolution of free surface.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.55-58
• /
• 2002

Parallel finite element code has been recently developed for the analysis of the incompressible Wavier-Stokes equations using domain decomposition method. Metis and MPI libraries are used for the domain partitioning of an unstructured mesh and the data communication between sub-domains, respectively. For unsteady computation of the incompressible Navier-Stokes equations, 4-step splitting method is combined with P1P1 finite element formulation. Smagorinsky and dynamic model are implemented for the simulation of turbulent flows. For the validation performance-estimation of the developed parallel code, three-dimensional Laplace equation has been solved. It has been found that the speed-up of 40 has been obtained from the present parallel code fir the bench mark problem. Lastly, the turbulent flows around the MIRA model and Tiburon model have been solved using 32 processors on IBM SMP cluster and unstructured mesh. The computed drag coefficient agrees better with the existing experiment as the mesh resolution of the region increases, where the variation of pressure is severe.

• Structural Engineering and Mechanics
• /
• 제18권5호
• /
• pp.691-707
• /
• 2004

This paper presents the micro-mechanical modeling for predicting concrete behavior under compressive loading. The model is able to represent the heterogeneities in the microstructure up to three phases, i.e., aggregate particles, matrix and interfaces. The smeared crack concept based on non-linear fracture mechanics is implemented in order to formulate the constitutive relation for each component. The splitting tensile strength is considered as a fracture criterion for cracking in micro-level. The finite element method is employed to simulate the model based on plane stress condition by using quadratic triangular elements. The validation of the model is verified by comparing with the experimental results. The influence of tensile strength from both aggregate and matrix phases on the concrete compressive strength is demonstrated. In addition, a guideline on selecting appropriate tensile strength for each phase to obtain specified concrete compressive strength is also presented.

• 대한전기학회논문지
• /
• 제41권7호
• /
• pp.735-741
• /
• 1992

This paper presents a parallel algorithm for the finite element analysis using relatively inexpensive transputer parallel system. The substructure method, which is highly parallel in nature, is used to improve the parallel computing efficiency by splitting up the whole structure into substructures. The proposed algorithm is applied to a simple two-dimensional magnetostatic problem. It is found that the more the number of transputer is increased, the more the total computation time is reduced. And the computational efficiency becomes better as the number of internal boundary nodes becomes smaller.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2002년도 춘계 학술대회논문집
• /
• pp.8-15
• /
• 2002

For the large scale computation of turbulent flows around an arbitrarily shaped body, a parallel LES (large eddy simulation) code has been recently developed in which domain decomposition method is adopted. METIS and MPI (message Passing interface) libraries are used for domain partitioning and data communication between processors, respectively. For unsteady computation of the incompressible Wavier-Stokes equation, 4-step splitting finite element algorithm [1] is adopted and Smagorinsky or dynamic LES model can be chosen fur the modeling of small eddies in turbulent flows. For the validation and performance-estimation of the parallel code, a three-dimensional laminar flow generated by natural convection inside a cube has been solved. Then, we have solved the turbulent flow around MIRA (Motor Industry Research Association) model at $Re = 2.6\times10^6$, which is based on the model height and inlet free stream velocity, using 32 processors on IBM SMP cluster and compared with the existing experiment.

• 한국해안해양공학회지
• /
• 제10권1호
• /
• pp.27-36
• /
• 1998

물의 유동을 해석하기 위한 3차원 수치모형이 개발되었다 모형은 균질류에 대한 $\sigma$-좌표에서 방정식들 을 유한요소법을 사용하여 해석한다. 모형의 정확성을 정토하기 위하여 1차원 수로에서 취송류 분포, 정사각형 호수에서 취송류 분포를 해석하고 해석해와 비교 검증하였으며, 마산-진해만에서 조류분포를 계산하고 현장 관측자료와 비교 검증하였다. 계산결과가 비교된 해석해 및 관측치와 대체로 일치하는 양호한 결과를 보였다. 따라서, 개발된 모형은 복잡한 육지경계를 갖는 자연 수괴의 3차원적 순환현상을 해석하는 데 널리 활용될 수 있을 것이다