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

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.91-98
• /
• 2004

In two dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a stabilized finite element formulation. The stabilization has been peformed by a modified residual method proposed by Illinca et. al.[3]. The stabilization which is necessary to escape from the LBB constraint renders an equal order formulation. In this paper, we increased the order of interpolation whithin an element up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.57-65
• /
• 2004

The incompressible Navier-Stokes equations in two dimensions are stabilized by a modified residual method, and then discretized by hierarchical elements. The stabilization is necessary to escape from the Ladyzhenskaya-Babuska-Brezzi(LBB) constraint and hence to achieve an equal order formulation. To expedite a standard iterative method such as the conjugate gradient squared(CGS) method, a preconditioning technique called the Hierarchical Iterative Procedure(HIP) has been applied. In this paper, we increased the order of interpolation within an element up to cubic. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far The numerical results by the present HIP for the lid driven cavity flow and others showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.67-76
• /
• 1992

A hierarchical formulation based on p-version of the finite element method for linear elastic axisymmetric stress analysis is presented. This is accomplished by introducing additional nodal variables in the element displacement approximation on the basis of integrals of Legendre polynomials. Since the displacement approximation is hierarchical, the resulting element stiffness matrix and equivalent nodal load vectors are hierarchical also. The merits of the propoosed element are as follow: i) improved conditioning, ii) ease of joining finite elements of different polynomial order, and iii) utilizing previous solutions and computation when attempting a refinement. Numerical examples are presented to demonstrate the accuracy, efficiency, modeling convenience, robustness and overall superiority of the present formulation. The results obtained from the present formulation are also compared with those available in the literature as well as with the analytical solutions.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.05a
• /
• pp.216-221
• /
• 2004

In two dimensional incompressible flaws, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint and hence achieve an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flaw analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flaw showed the present procedure to be stable, very efficient and useful in flaw analyses in conjunction with hierarchical elements.

• Journal of Ocean Engineering and Technology
• /
• v.17 no.5 s.54
• /
• pp.11-18
• /
• 2003

In two-dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure (HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint hence, achieving an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared (CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper and efficient iterative procedure for higher order finite element formulation has not been available, thus far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient, and useful in flow analyses, in conjunction with hierarchical elements.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.9
• /
• pp.1209-1217
• /
• 2002

The use of a two dimensional hierarchical elements are investigated for the incompressible flow computation. The construction of hierarchical elements are explained by both a geometric configuration and a determination of degrees of freedom. Also 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 shows that the computation with hierarchical higher order elements can increase the convergence rate and accuracy of finite element solutions in more efficient manner than the use of standard first order element. for Stokes and Cavity flow cases, a mixed version of penalty function approach has been introduced in connection with the hierarchical elements. Solutions from hierarchical elements showed better resolutions with consistent trends in both mesh shapes and the order of elements.

• Journal of Power System Engineering
• /
• v.4 no.1
• /
• pp.68-73
• /
• 2000

A mixed degree finite element solutions using hierarchical elements are investigated for convergences on a 2-D simple cases. Elements are generated block by block and each block is assigned an arbitrary solution degree. The numerical study showed that a well constructed blocks can increase the convergence and accuracy of finite element solutions. Also, it has been found that for higher order elements, the convergence trends can be deteriorated for smaller mesh sizes. A procedure for a variable fixed boundary condition has been included.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.1
• /
• pp.117-124
• /
• 2009

The p-convergent boundary element method has been proposed to analyze two-dimensional potential problem on the basis of high order Legendre shape functions that have different property comparing with the shape functions in conventional boundary element method. The location of nodes corresponding to high order shape function are not defined along the boundary, called by nodeless node, similar to the p-convergent finite element method. As the order of shape function increases, the collocation point method is used to solve linear simultaneous equations. The collocation patterns of p-convergent boundary element method consist of non-symmetric hierarchial or symmetric non-hierarchical. As the order of shape function increases, the number of collocation point increases. The singular integral that appears in p-convergent boundary element has been calculated by special numeric quadrature technique and semi-analytical integration technique. The L-shape domain problem including singularity in the vicinity of reentrant comer is analyzed and the numerical results show that the relative error is smaller than $10^{-2}%$ range as compared with other results in literatures. In case of same condition, the symmetric p-collocation point pattern shows high accuracy of solution.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.8
• /
• pp.1007-1013
• /
• 2010

A new local refinement scheme for spline finite element method has been proposed; this scheme involves the use of hierarchical B-spline. NURBS has been widely used in CAD; however, the local refinement of NURBS is difficult due to its tensor-product property. In this study, we attempted to use hierarchical B-splines as local refinement strategy in spline FEM. The regions of high gradients are overlapped by hierarchically-created local meshes. Knot vectors and control points in local meshes are extracted from global meshes, and they are refined using specific schemes. Proper compatibility conditions are imposed between global and local meshes. The effectiveness of the proposed method is verified on the basis of numerical results. Further, it is shown that by using a proposed local refinement scheme, the accuracy of the solution can be improved and it could be higher than that of the solution of a conventional spline FEM with relatively lower degrees of freedom.