• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.57-70
• /
• 2008

A high-order accurate Euler flow solver based on a discontinuous Galerkin finite-element method has been developed for the numerical simulations of blade-vortex interaction phenomena on unstructured meshes. A free vortex in freestream was investigated to assess the vortex-preserving property and the accuracy of the present flow solver. Blade-vortex interaction problems in subsonic and transonic freestreams were simulated by adopting a multi-level solution-adaptive dynamic mesh refinement/coarsening technique. The results were compared with those of other numerical and experimental methods. It was shown that the present discontinuous Galerkin flow solver can preserve the vortex structure for significantly longer vortex convection time and can accurately capture the complex unsteady blade-vortex interaction flows, including generation and propagation of acoustic waves.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2003.06a
• /
• pp.1605-1608
• /
• 2003

A unstructured tetrahedral mesh generation algorithm has been presented. To make better meshes in interior region using an advancing front technique, a connecting operator has been developed in addition to the existing operators. Before applying digging operators that generate new nodes inside of a meshing region, a connecting operator is employed that uses existing nodes which satisfy certain conditions for producing well-conditioned elements if possible. By introducing this new operator, tetrahedral meshing process becomes more robust and produces better quality of meshes.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.90-97
• /
• 1999

An implicit turbulent flow solver is developed for 2-D unstructured hybrid meshes. Spatial discretization is accomplished by a cell-centered finite volume formulation using an upwind flux differencing. Time is advanced by an implicit backward Euler time stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one equation model, which is coupled with wall function. The numerical method is applied for flows on a flat plate, the NACA 0012 airfoil, and the Douglas 3 element airfoil. The results are compared with experimental data.

• International Journal of CAD/CAM
• /
• v.1 no.1
• /
• pp.11-21
• /
• 2002

This paper presents a triangular-to-quadrilateral mesh conversion method that can control the directionality of the output quadrilateral mesh according to a user-specified vector field. Given a triangular mesh and a vector field, the method first scores all possible quadrilaterals that can be formed by pairs of adjacent triangles, according to their shape and directionality. It then converts the pairs into quadrilateral elements in order of the scores to form a quadrilateral mesh. Engineering analyses with finite element methods occasionally require a quadrilateral mesh well aligned along the boundary geometry or the directionality of some physical phenomena, such as in the directions of a streamline, shock boundary, or force propagation vectors. The mesh conversion method can control the mesh directionality according to any desired vector fields, and the method can be used with any existing triangular mesh generators.

• Journal of the Korean Society for Precision Engineering
• /
• v.21 no.3
• /
• pp.91-99
• /
• 2004

An unstructured tetrahedral mesh generation algorithm has been presented. In order to construct better meshes in interior region by using an advancing front technique, a connecting operator and a local finishing operator II have been developed in addition to the existing operators. Before applying digging operators that generate new nodes inside of a meshing region, a connecting operator is employed that uses existing nodes which satisfy certain conditions for producing well-conditioned elements. The local finishing operator II is introduced to terminate the meshing process more flexibly on remaining subregions. With these new operators, tetrahedral meshing process becomes more robust and good quality of meshes are constructed.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.27 no.2
• /
• pp.77-82
• /
• 2023

The structural analysis module is an essential part of any integrated structural system. Diverse integrated systems today require, from the analysis module, efficient real-time responses to real-time input such as earthquake signals, extreme weather-related forces, and man-made accidents. An integrated system may also be for the entire life span of a civil structure conceived during the initial conception, developed throughout various design stages, effectively used in construction, and utilized during usage and maintenance. All these integrated systems' essential part is the structural analysis module, which must be automated and computationally efficient so that responses may be almost immediate. The finite element method is often used for structural analysis, and for automation, many effective finite element meshes must be automatically generated for a given analysis. A computationally efficient finite element mesh generation scheme based on the r-h method of mesh refinement using strain deviations from the values at the Gauss points as error estimates from the previous mesh is described. Shape factors are used to sort out overly distorted elements. A standard cantilever beam analyzed by four-node plane stress elements is used as an example to show the effectiveness of the automated algorithm for a time-domain dynamic analysis. Although recent developments in computer hardware and software have made many new applications in integrated structural systems possible, structural analysis still needs to be executed efficiently in real-time. The algorithm applies to diverse integrated systems, including nonlinear analyses and general dynamic problems in earthquake engineering.

• Interaction and multiscale mechanics
• /
• v.3 no.3
• /
• pp.235-255
• /
• 2010

This paper investigates the heat equation for domains subjected to an internal source with a sharp spatial gradient. The solution is first approximated using linear finite elements, and sufficiently small time-step sizes to yield stable simulations. The main area of interest is then in the ability to approximate the solution using Generalized Finite Elements, and again explore the time-step limitations required for stable simulations. Both high order elements, as well as elements with special enrichments are used to generate solutions. When compared to linear finite elements, the high order elements deliver better accuracy at a given level of mesh refinement, but do not offer an increase in critical time-step size. When special enrichment functions are used, the solution can be approximated accurately on very coarse meshes, while yielding solutions which are both accurate and computationally efficient. The major conclusion of interest is that the significantly larger element size yields larger allowable time-step sizes while still maintaining stability of the time-stepping algorithm.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.7
• /
• pp.2078-2086
• /
• 1996

This paper introduces two three-dimensional eight-node hexagonal elements obtained by using multivariable variational mehtod. Both of them are based on the modified hellinger-reissner principle to employ incompatible displacements and assumed stresses of assumed strains. The internal functions of element are introduced to as element formulation through two different methods : the first one uses the functions determined directly from the element boundary condition of the incompatible displacements ; while the second, being a kind of B-bar mehtod, employs the modification technique of strain-displacement matrix to pass the patch test. The elements are evaluated on the selective problems of bending and material incompressibility with regular and distorted meshes. The results show that the new elements perform with good accuracy in both of deformation and stress calculation and they are insensitive to distorted geometry of element.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2002.04a
• /
• pp.105-112
• /
• 2002

A new grid-based approach is presented for automatic generation of hexahedral meshes for simulation of plastic deformation in metal forming. In this approach, special enveloping schemes are applied, to eradicate the sources of the degenerate elements that may appear in a generated mesh. The schemes are described in detail, along with a complete procedure for mesh generation. The capability of the approach to deal with an arbitrary, 3-D process geometry is demonstrated through application to a selected forming problem.

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

The Zienkiewicz-Zhu(Z/Z) error estimate is slightly modified for the hierarchical p-refinement, and is then applied to L-shaped plates subjected to bending to demonstrate its effectiveness. An adaptive procedure in finite element analysis is presented by p-refinement of meshes in conjunction with a posteriori error estimator that is based on the superconvergent patch recovery(SPR) technique. The modified Z/Z error estimate p-refinement is different from the conventional approach because the high order shape functions based on integrals of Legendre polynomials are used to interpolate displacements within an element, on the other hand, the same order of basis function based on Pascal's triangle tree is also used to interpolate recovered stresses. The least-square method is used to fit a polynomial to the stresses computed at the sampling points. The strategy of finding a nearly optimal distribution of polynomial degrees on a fixed finite element mesh is discussed such that a particular element has to be refined automatically to obtain an acceptable level of accuracy by increasing p-levels non-uniformly or selectively. It is noted that the error decreases rapidly with an increase in the number of degrees of freedom and the sequences of p-distributions obtained by the proposed error indicator closely follow the optimal trajectory.