• Bulletin of the Korean Mathematical Society
• /
• v.51 no.6
• /
• pp.1655-1668
• /
• 2014

In this article, we propose and analyze a new nonconforming primal mixed finite element method for the stationary Stokes equations. The approximation is based on the pseudostress-velocity formulation. The incompressibility condition is used to eliminate the pressure variable in terms of trace-free pseudostress. The pressure is then computed from a simple post-processing technique. Unique solvability and optimal convergence are proved. Numerical examples are given to illustrate the performance 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 the Korean Society for Industrial and Applied Mathematics
• /
• v.16 no.1
• /
• pp.1-13
• /
• 2012

In this paper we present an a posteriori error estimator for the stabilized P1 nonconforming finite element method of the linear elasticity problem based on a nonsymmetric H(div)-conforming approximation of the stress tensor in the first-order Raviart-Thomas space. By combining the equilibrated residual method and the hypercircle method, it is shown that the error estimator gives a fully computable upper bound on the actual error. Numerical results are provided to confirm the theory and illustrate the effectiveness of our error estimator.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.12 no.2
• /
• pp.69-79
• /
• 2008

Multigrid methods finite element/finite volume methods and their convergence properties are reviewed in a general setting. Some early theoretical results in simple finite element methods in variational setting method are given and extension to nonnested-noninherited forms are presented. Finally, the parallel theory for nonconforming element[13] and for cell centered finite difference methods [15, 23] are discussed.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.1 no.1
• /
• pp.105-125
• /
• 1997

In this paper we consider multigrid algorithms for solving elasticity problems by using Wilson's nonconforming finite element method. We consider two types of intergrid transfer operators which is needed to define the multigrid algorithm and prove convergence of $\mathcal{W}$-cycle mutigrid algorithm and uniform condition number estimates for the variable $\mathcal{V}$-cycle multigrid preconditioner.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.5 no.1
• /
• pp.55-78
• /
• 2001

Given the anisotropic Poisson equation $-{\nabla}{\cdot}{\mathcal{K}}{\nabla}p=f$, one can convert it into a system of two first order PDEs: the Darcy law for the flux $u=-{\mathcal{K}{\nabla}p$ and conservation of mass ${\nabla}{\cdot}u=f$. A very natural mixed finite volume method for this system is to seek the pressure in the nonconforming P1 space and the Darcy velocity in the lowest order Raviart-Thomas space. The equations for these variables are obtained by integrating the two first order systems over the triangular volumes. In this paper we show that such a method is really a standard finite element method with local recovery of the flux in disguise. As a consequence, we compare two approaches in analyzing finite volume methods (FVM) and shed light on the proper way of analyzing non co-volume type of FVM. Numerical results for Dirichlet and Neumann problems are included.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.10
• /
• pp.1139-1146
• /
• 2012

An alternative approach for topology optimization of steady incompressible Navier-Stokes flow problems is presented by using P1 nonconforming finite elements. This study is the extended research of the earlier application of P1 nonconforming elements to topology optimization of Stokes problems. The advantages of the P1 nonconforming elements for topology optimization of incompressible materials based on locking-free property and linear shape functions are investigated if they are also valid in fluid equations with the inertia term. Compared with a mixed finite element formulation, the number of degrees of freedom of P1 nonconforming elements is reduced by using the discrete divergence-free property; the continuity equation of incompressible flow can be imposed by using the penalty method into the momentum equation. The effect of penalty parameters on the solution accuracy and proper bounds will be investigated. While nodes of most quadrilateral nonconforming elements are located at the midpoints of element edges and higher order shape functions are used, the present P1 nonconforming elements have P1, {1, x, y}, shape functions and vertex-wisely defined degrees of freedom. So its implentation is as simple as in the standard bilinear conforming elements. The effectiveness of the proposed formulation is verified by showing examples with various Reynolds numbers.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2001.10a
• /
• pp.191-194
• /
• 2001

A partially coupled thermo-piezoelectric-mechanical triangular finite element model of composite laminates with surface bonded piezoelectric actuators, subjected to externally applied mechanical load, temperature change load, electric field load is developed. The governing differential equations are obtained by applying the principle of free energy and variational techniques. A higher order zigzag theory displacement field is employed to accurately capture the transverse shear and normal effects in laminated composite plates of arbitrary thickness. Nonconforming shape functions by Specht are employed in the transverse displacement variables. Numerical examples demonstrate the accuracy and efficiency of the proposed triangular plate element.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.4
• /
• pp.323-329
• /
• 2013

This study attempts to simulate the structure of a woody leaf netted venation system by using topology optimization techniques. Based on finite element method (FEM) analysis of an incompressible fluid, a topology optimal design is applied to those woody leaf netted venation models. To solve the transverse shear locking problem of a thin plate caused by the Mindlin-Reissner plate model where a leaf netted venation is assumed to be a thin plate, a P1-nonconforming element and selective reduced integration are employed. Topology optimal design is applied to multiple physical domains. Combined with the Darcy-Stokes flow problems and extended to the optimal design of fluid channels, the multiple physical models of the flow system are analyzed and venation patterns of leafs are simulated. The calculated optimal shapes are compared with the natural shapes of woody leaf venation patterns. This interdisciplinary approach may improve our understanding of the leaf venation system.

• Computational Structural Engineering
• /
• v.7 no.2
• /
• pp.101-109
• /
• 1994

A degenerated cylindrical shell element for modeling glass fiber reinforced plastic pipes is developed and its performance for static structural analysis under internal uniform pressure is evaluated. The element is a nine node degenerated solid shell element with reduced integration technique, addition of nonconforming displacement modes, and assumed strain method to improve convergence of analysis. Several numerical examples are solved and compared with analytical solutions and other F.E.M programs, The results show that the increment of fiber orientation in the GFRP pipes with reference to the longitudinal axis cause less radial displacements and much stiffness in the pipes. This is reasonable since the internal pressure will primarily cause hoop stresses in the ring and 90-angle ply GFRP ring carry these efficiently in pure tension.