• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.693-702
• /
• 2010

In this paper, FETI(Finite Element Tearing and Interconnecting) method is introduced in order to improve numerical efficiency of Staggered method. The porous media theory, the Staggered method and the FETI method are briefly introduced in this paper. In addition, we account for the MPI(Message Passing Interface) library for parallel analysis, and the proposed combined Staggered method with FETI method. Finally Lagrange multipliers and CG(Conjugate Gradient) algorithm to solve decomposed domain are proposed, and then the proposed method is verified to be numerically efficient by MPI library.

• 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 for Industrial and Applied Mathematics
• /
• v.12 no.4
• /
• pp.249-260
• /
• 2008

In this work we propose a mixed finite volume method for the Signorini problem which are based on the idea of Keller's finite volume box method. The triangulation may consist of both triangles and quadrilaterals. We choose the first-order nonconforming space for the scalar approximation and the lowest-order Raviart-Thomas vector space for the vector approximation. It will be shown that our mixed finite volume method is equivalent to the standard nonconforming finite element method for the scalar variable with a slightly modified right-hand side, which are crucially used in a priori error analysis.

• Journal of applied mathematics & informatics
• /
• v.15 no.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.

• Structural Engineering and Mechanics
• /
• v.63 no.2
• /
• pp.137-146
• /
• 2017

The size-dependent behavior of single walled carbon nanotubes (SWCNT) embedded in the elastic medium and subjected to the initial axial force is investigated using the mixed finite element method. The SWCNT is assumed to be Euler-Bernoulli beam incorporating nonlocal theory developed by Eringen. The mixed finite element model shows its great advantage of dealing with nonlocal behavior of SWCNT subjected to a concentrated load owing to the existence of two coefficients ${\alpha}_1$ and ${\alpha}_2$. This is the first numerical approach to deal with a puzzling fact of nonlocal theory with concentrated load. Numerical examples are performed to show the accuracy and efficiency of the present method. In addition, parametric study is carefully carried out to point out the influences of nonlocal effect, the elastic medium, and the initial axial force on the behavior of the carbon nanotubes.

• Korean Journal of Mathematics
• /
• v.13 no.2
• /
• pp.127-137
• /
• 2005

The objective of numerical analysis is to devise and analyze efficient algorithms or numerical methods for equations arising in mathematical modeling for science and engineering. In this article, we present some recent topics in computational mathematics, specially in the finite element method and overview the development of the mixed finite element method in the context of second order elliptic and parabolic problems. Multiscale methods such as MsFEM, HMM, and VMsM are included.

• Structural Engineering and Mechanics
• /
• v.20 no.4
• /
• pp.405-420
• /
• 2005

In this paper, a hybrid/mixed nonlinear shell element is developed in polar coordinate system based on Hellinger/Reissner variational principle and the large-deflection theory of plate. A numerical solution scheme is formulated using the hybrid/mixed finite element method (HMFEM), in which the nodal values of bending moments and the deflection are the unknown discrete parameters. Stability of the present element is studied. The large-deflection analyses are performed for simple supported and clamped circular plates under uniformly distributed and concentrated loads using HMFEM and the traditional displacement finite element method. A parametric study is also conducted in the research. The accuracy of the shell element is investigated using numerical computations. Comparisons of numerical solutions are made with theoretical results, finite element analysis and the available numerical results. Excellent agreements are shown.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.4
• /
• pp.369-378
• /
• 2010

The mixed finite element analysis is the most widely used method for saturated porous media. Generally, in this method, direct method and iterative method are proposed to obtain unknown variable, however, the iterative method is recommended because the method provide numerical stability and accuracy under the material properties for solid and fluid are different. In this paper, we introduce staggered method which has strong numerical stability, and FETI(Finite Element Tearing and Interconnecting) which is one of decomposition methods are applied into the method in order to obtain numerical efficiency. In which, Lagrange Multipliers and conjugated gradient method to solve decomposed domain are proposed, and then, the proposed method is verified numerical efficiency by point to point MPI(Message Passing Interface) library.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.20-25
• /
• 2008

The finite element analysis for porous media is severe job because constituents have different physical peoperties, and element's continuity and stability should be considered. Thus, we propose the new mixed finite element method in order to overcome the problems. In this method, multi time step, remeshing step, and sub iteration step are introduced. The multi time step and remeshing step make it possible to satisfy a stability and an accuracy during sub iteration in which global time is determined. Finally, the proposed method is compared with the ABAQUS(2007) software and exact solution(Schiffman 1967) through two dimensional consolidation model.

• Structural Engineering and Mechanics
• /
• v.13 no.3
• /
• pp.277-298
• /
• 2002

In this study, free vibration analysis of Reissner plates on Pasternak foundation is carried out by mixed finite element method based on the G$\hat{a}$teaux differential. New boundary conditions are established for plates on Pasternak foundation. This method is developed and applied to numerous problems by Ak$\ddot{o}$z and his co-workers. In dynamic analysis, the problem reduces to the solution of a standard eigenvalue problem and the mixed element is based upon a consistent mass matrix formulation. The element has four nodes and bending and torsional moments, transverse shear forces, rotations and displacements are the basic unknowns. The element performance is assessed by comparison with numerical examples known from literature. Validity limits of Kirchhoff plate theory is tested by dynamic analysis. Shear locking effects are tested as far as $h/2a=10^{-6}$ and it is observed that REC32 is free from shear locking.