• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.15 no.1
• /
• pp.57-65
• /
• 2011

In this paper, a good quality mesh generation for the finite element method is investigated for artificial hip joint simulations. In general, bad meshes with a large aspect ratio or mixed elements can give rise to excessively long computational running times and extremely high errors. Typically, hexahedral elements outperform tetrahedral elements during three-dimensional contact analysis using the finite element method. Therefore, it is essential to mesh biologic structures with hexahedral elements. Four meshing schemes for the finite element analysis of an artificial hip joint are presented and compared: (1) tetrahedral elements, (2) wedge and hexahedral elements, (3) open cubic box hexahedral elements, and (4) proposed hexahedral elements. The proposed meshing scheme is to partition a part before seeding so that we have a high quality three-dimensional mesh which consists of only hexahedral elements. The von Mises stress distributions were obtained and analyzed. We also performed mesh refinement convergence tests for all four cases.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.4
• /
• pp.395-405
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among four filling patterns at each triangular control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. Using the proposed numerical technique, the collapse of a water dam has been analyzed to predict flow phenomenon of fluid and the predicted front positions with respect to time have been compared with the reported experimental results.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.55 no.6
• /
• pp.306-313
• /
• 2006

A three dimensional adaptive finite element refinement algorithm is developed for non-linear magnetostatic field problems. In the method, the edge elements are used for finite element formulation, and the local error in each element is estimated from the fact that the tangential components of magnetic field intensity and the normal components of magnetic flux density should be continuous at the interface of the two adjacent elements. Based on the estimated error, the elements which have big error are divided into several elements using bisection method. The effectiveness of the developed algorithm is proved through numerical examples.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.4 s.175
• /
• pp.1016-1023
• /
• 2000

The present work is concerned with three-dimensional finite element analysis of the hollow section extrusion process using a porthole die. The effects of related design parameters are discussed through the finite element simulation for extrusion of a triply-connected rectangular tubular section. For economic computation, mismatching refinement, an efficient domain decomposition method with different mesh density for each subdomain, is implemented. In order to obtain the uniform flow at the outlet, design parameters such as the hole size and the hole position are investigated and compared through the numerical analysis. Comparing the velocity distribution with that of the original design, it is concluded that the design modification enables more uniform flow characteristics. The analysis results are then successfully reflected on the industrial porthole die design.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.146-152
• /
• 2000

Two and Three dimensional dynamic adaptation code is developed for transient computations. This code involves mesh refinement and coarsening to either add points in high gradient regions of flow or remove points where they are not needed, for high spatial accuracy. Temporary cell algorithm is used to maintain the original grid quality. To show the assessment of the accuracy and efficiency, two dimensional study and unsteady flows are computed. Also, three dimensional steady computations are made to assess the refinement using temporary cell algorithm. The result shows the high spatial accuracy primarily in discontinuity regions in steady and unsteady computation.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1996.10a
• /
• pp.45-51
• /
• 1996

The combined shear and compression behaviors of seismic isolation rubber bearings are analyzed using the hyperelasticity material option of the ABAQUS computer program. The purpose of the analysis is to predict the behavior of laminated rubber bearing before the several tests. Some kinds of strain energy density functions are used as constitutive law for rubber itself having the hyperelasticity. The results are compared with test data peformed in Italy The analysis results show a little different with experimental results depending on the constitutive model and the refinement of finite element. The high order form of strain energy density functions results in good agreements and the mesh refinement above two for one rubber layer is enough to get good results.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.127-132
• /
• 1998

The present study is aimed to investigate flow characteristics of Two dimensional backward-facing step by numerical approach. A convection conservative difference scheme based upon SOLA algorithm is used for the solution of the two-dimensional incompressible Navier-Stokes equations to simulate the laminar, transitional and turbulent flow conditions at which the experimental data can be available for the backward-facing step. The twenty kinds of Reynolds number are used for the calculations. In an effort to demonstrate that the reported solutions are dependent on the mesh refinement, computations are performed on seven different meshes of uniformly increasing refinement. Also to investigate the result of inflow dependence, two kinds of the inflow profile are chosen for the laminar flow. As criterion of benchmarking the result of numerical simulation, reattachment length is used for the selected Reynolds numbers.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.271-274
• /
• 2007

The finite element linear buckling analysis of folded plate structures using adaptive h-refinement methods is presented in this paper. The variable-node flat shell element used in this study possesses the drilling D.O.F. which, in addition to improvement of the element behavior, permits an easy connection to other elements with six degrees of freedom per node. The Box-typed structures can be analyzed using these developed flat shell elements. By introducing the variable node elements some difficulties associated with connecting the different layer patterns, which are common in the adaptive h-refinement on quadrilateral mesh, can be overcome. To obtain better stress field for the error estimation, the super-convergent patch recovery is used. The convergent buckling modes and the critical loads associated with these modes can be obtained.

• The Bulletin of The Korean Astronomical Society
• /
• v.38 no.2
• /
• pp.101.2-101.2
• /
• 2013

We present a new computational fluid dynamic (CFD) simulation code. The code employs the moving and polyhedral unstructured mesh scheme, which is known as a superior approach to the conventional SPH (smoothed particle hydrodynamics) and AMR (adaptive mesh refinement) schemes. The code first generates unstructured meshes by the Voronoi tessellation at every time step, and then solves the Riemann problem for surfaces of every Voronoi cell to update the hydrodynamic states as well as to move former generated meshes. For the second-order accuracy, the MUSCL-Hancock scheme is implemented. To increase efficiency for generating Voronoi tessellation we also develop the incremental expanding method, by which the CPU time is turned out to be just proportional to the number of particles, i.e., O(N). We will discuss the applications of our code in the context of cosmological simulations as well as numerical experiments for galaxy formation.

• International Journal of Aeronautical and Space Sciences
• /
• v.8 no.1
• /
• pp.66-78
• /
• 2007

A multi-point aerodynamic shape optimization technique has been developed for helicopter rotor blades in hover based on a continuous adjoint method on unstructured meshes. The Euler flow solver and the continuous adjoint sensitivity analysis were formulated on the rotating frame of reference. The 'objective function and the sensitivity were obtained as a weighted sum of the values at each design point. The blade section contour was modified by using the Hicks-Henne shape functions. The mesh movement due to the blade geometry change was achieved by using a spring analogy. In order to handle the repeated evaluation of the design cycle efficiently, the flow and adjoint solvers were parallelized based on a domain decomposition strategy. A solution-adaptive mesh refinement technique was adopted for the accurate capturing of the wake. Applications were made to the aerodynamic shape optimization of the Caradonna-Tung rotor blades and the UH-60 rotor blades in hover.