• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2010.11a
• /
• pp.631-632
• /
• 2010

• 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.

• 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.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 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.