• Journal of the Korean Mathematical Society
• /
• v.44 no.5
• /
• pp.1103-1119
• /
• 2007

We consider multiscale mortar mixed finite element discretizations for slightly compressible Darcy flows in porous media. This paper is an extension of the formulation introduced by Arbogast et al. for the incompressible problem [2]. In this method, flux continuity is imposed via a mortar finite element space on a coarse grid scale, while the equations in the coarse elements (or subdomains) are discretized on a fine grid scale. Optimal fine scale convergence is obtained by an appropriate choice of mortar grid and polynomial degree of approximation. Parallel numerical simulations on some multiscale benchmark problems are given to show the efficiency and effectiveness of the method.

• Journal of the Korean Society of Industry Convergence
• /
• v.7 no.4
• /
• pp.355-361
• /
• 2004

The finite element method were used to determine the stress intensity factor of cracked plate. The stress method, displacement method and J Integral are most popular finte element method. ANSYS proposed another a kind of displacement method. In this paper, it was examined that the accuracy and utility of the ANSYS method could believable to determine the stress intensity factors of centered inclined crack. Generally, inclined crack has two portion of stress intensity factors, tensile mode F1 and shear mode F2. For the purpose of increasing the accuracy of stress intensity factors, examined the effect of the numbers of nodes and elements, crack tip element size and number of partition of the crack tip vicinity. It was found that the method proposed by ANSYS is useful and has high accuracy. Accuracy of calculated stress intensity factors was increased by increase of the number of nodes and elements, and at the small size of crack tip elements can get more highly accuracy.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.6 no.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.

• Structural Engineering and Mechanics
• /
• v.73 no.1
• /
• pp.1-16
• /
• 2020

Finite elements based on the partition of unity (PU) approximation have powerful capabilities for p-adaptivity and solutions with high smoothness without remeshing of the domain. Recently, the PU approximation was successfully applied to the three-node shell finite element, properly eliminating transverse shear locking and showing excellent convergence properties and solution accuracy. However, the enrichment with the PU approximation results in a significant increase in the number of degrees of freedom; therefore, it requires greater computational cost, thus making it less suitable for practical engineering. To circumvent this disadvantage, we propose a new strategy to decrease the total number of degrees of freedom in the existing PU-based shell element, without loss of optimal convergence and accuracy. To alleviate the locking phenomenon, we use the method of mixed interpolation of tensorial components and perform convergence studies to show the accuracy and capability of the proposed shell element. The excellent performances of the new shell elements are illustrated in three benchmark problems.

• Coupled systems mechanics
• /
• v.7 no.6
• /
• pp.649-668
• /
• 2018

In this paper, we present a numerical model for fluid-structure interaction between structure built of porous media and acoustic fluid, which provides both pore pressure inside porous media and hydrodynamic pressures and hydrodynamic forces exerted on the upstream face of the structure in an unified manner and simplifies fluid-structure interaction problems. The first original feature of the proposed model concerns the structure built of saturated porous medium whose response is obtained with coupled discrete beam lattice model, which is based on Voronoi cell representation with cohesive links as linear elastic Timoshenko beam finite elements. The motion of the pore fluid is governed by Darcy's law, and the coupling between the solid phase and the pore fluid is introduced in the model through Biot's porous media theory. The pore pressure field is discretized with CST (Constant Strain Triangle) finite elements, which coincide with Delaunay triangles. By exploiting Hammer quadrature rule for numerical integration on CST elements, and duality property between Voronoi diagram and Delaunay triangulation, the numerical implementation of the coupling results with an additional pore pressure degree of freedom placed at each node of a Timoshenko beam finite element. The second original point of the model concerns the motion of the outside fluid which is modeled with mixed displacement/pressure based formulation. The chosen finite element representations of the structure response and the outside fluid motion ensures for the structure and fluid finite elements to be connected directly at the common nodes at the fluid-structure interface, because they share both the displacement and the pressure degrees of freedom. Numerical simulations presented in this paper show an excellent agreement between the numerically obtained results and the analytical solutions.

• Journal of the Korean GEO-environmental Society
• /
• v.18 no.2
• /
• pp.53-58
• /
• 2017

In this paper, we propose a new method to debond between mixed finite elements for porous media in ABAQUS (2014). ABAQUS just provides debonding algorithm for the u-p model using cohesive elements in standard version. However, this approach has a drawback that it is hard to simulate complex debonding problems like element separation, rigid body motion, and contact between separated elements in standard version. ABAQUS-explicit can resolve these complex problems, but cohesive elements for the u-p model cannot be applied. We introduce a new algorithm for debonding for porous media instead of using cohesive elements. In this method, subroutines VUMAT to apply constitutive models and VDISP to separate elements in ABAQUS are used to simulate debonding problems. In addition, a simple 2-D example is demonstrated in the ABAQUS-explicit solver.

• Steel and Composite Structures
• /
• v.28 no.6
• /
• pp.655-670
• /
• 2018

A coupled method, that combines the Ritz method and the finite element (FE) method, is proposed to solve the vibration problem of rectangular thin and thick plates with general boundary conditions. The eigenvalue partial differential equation(s) of the plate is (are) first reduced to a set of eigenvalue ordinary differential equations by the application of the Ritz method. The resulting eigenvalue differential equations are then reduced to an eigenvalue algebraic equation system using the finite element method. The natural boundary conditions of the plate problem including the free edge and free corner boundary conditions are also implemented in a simple and accurate manner. Various boundary conditions including simply supported, clamped and free boundary conditions are considered. Comparisons with existing numerical and analytical solutions show that the proposed mixed method can produce highly accurate results for the problems considered using a small number of Ritz terms and finite elements. The proposed mixed Ritz-FE formulation is also compared with the mixed FE-Ritz formulation which has been recently proposed by the present author and his co-author. It is found that the proposed mixed Ritz-FE formulation is more efficient than the mixed FE-Ritz formulation for free vibration analysis of rectangular plates with Levy-type boundary conditions.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.425-429
• /
• 2006

The understanding and prediction of the behavior of flow in open channels are important to the solution of a wide variety of practical flow problems in water resources engineering. Recently, frequent drought has increased the necessity of an effective water resources control and management of river flows for reserving instream flow. The objective of this study is to develop an efficient and accurate finite element model based on Streamline Upwind/Petrov-Galerkin(SU/PG) scheme for analyzing and predicting two dimensional flow features in complex natural rivers. Several tests were performed in developed all elements(4-Node, 6-Node, 8-Node elements) for the purpose of validation and verification of the developed model. The U-shaped channel of flow and natural river of flow were performed for tests. The results were compared with these of laboratory experiments and RMA-2 model. Such results showed that solutions of high order elements were better accurate and improved than those of linear elements. Also, the suggested model displayed reasonable velocity distribution compare to RMA-2 model in meandering domain for application of natural river flow. Accordingly, the developed finite element model is feasible and produces reliable results for simulation of two dimensional natural river flow. Also, One contribution of this study is to present that results can lead to significant gain in analyzing the accurate flow behavior associated with hydraulic structure such as weir and water intake station and flow of chute and pool.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.422-429
• /
• 2001

A two dimensional hierarchical elements are investigated for a use on the incompressible flow computation. The construction of hierarchical elements are explained through the tensor product of 1-D hierarchical functions, and a systematic treatment of essential boundary values has been developed for the degrees of freedom corresponding to higher order terms. The numerical study for the poisson problem showed that the present scheme can increase the convergence and accuracy of finite element solutions, and can be more efficient than the standard first order with many elements. Also, for Stokes and cavity flow cases, solutions from hierarchical elements showed better resolutions and future promises for higher order solutions.

• Journal of Korea Water Resources Association
• /
• v.42 no.12
• /
• pp.1053-1067
• /
• 2009

The objective of this study is to develop an efficient and accurate quadratic finite element model based on Streamline Upwind/Petrov Galerkin (SU/PG) scheme for analyzing and predicting two dimensional flow features in complex natural rivers. For a development of model, quadratic tin, quadrilateral and mixed elements as well as linear tin, quadrilateral and mixed elements were used in the model. Also, this model was developed through reinforcement of Gauss Quadrature which was necessary to integral of governing equation. Several tests for bottom-rising channel and U-type channel were performed for the purpose of validation and verification of the developed model. Such results showed that solutions of second order elements are better accurate and improved than those of linear elements. Results obtained by the developed model and RMA-2 model are compared, and the results for the developed model were better accurate than those of RMA-2 model. In the future if the developed model is applied in natural rivers, it can provide better accurate results than those of existing model.