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

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.1
• /
• pp.49-56
• /
• 2024

An automated dynamic structural analysis module stands as a crucial element within a structural integrated mitigation system. This module must deliver prompt real-time responses to enable timely actions, such as evacuation or warnings, in response to the severity posed by the structural system. The finite element method, a widely adopted approximate structural analysis approach globally, owes its popularity in part to its user-friendly nature. However, the computational efficiency and accuracy of results depend on the user-provided finite element mesh, with the number of elements and their quality playing pivotal roles. This paper introduces a computationally efficient adaptive mesh generation scheme that optimally combines the h-method of node movement and the r-method of element division for mesh refinement. Adaptive mesh generation schemes automatically create finite element meshes, and in this case, representative strain values for a given mesh are employed for error estimates. When applied to dynamic problems analyzed in the time domain, meshes need to be modified at each time step, considering a few hundred or thousand steps. The algorithm's specifics are demonstrated through a standard cantilever beam example subjected to a concentrated load at the free end. Additionally, a portal frame example showcases the generation of various robust meshes. These examples illustrate the adaptive algorithm's capability to produce robust meshes, ensuring reasonable accuracy and efficient computing time. Moreover, the study highlights the potential for the scheme's effective application in complex structural dynamic problems, such as those subjected to seismic or erratic wind loads. It also emphasizes its suitability for general nonlinear analysis problems, establishing the versatility and reliability of the proposed adaptive mesh generation scheme.

• Journal of electromagnetic engineering and science
• /
• v.1 no.1
• /
• pp.78-82
• /
• 2001

For fast capacitance computation, a simple mesh refinement technique on MLFMA(Multi-Level Fast Multipole Algorithm) is proposed The triangular meshes are refined mainly in the area which has heavy charge density. The technique is applied to the capacitance extraction of three dimensional conductors. The results show good convergence with comparable accuracy. An adaptive technique concerned with MLFMA is useful to reduce computation time and the number of elements without additional computational efforts in large three dimensional problems.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.4
• /
• pp.295-306
• /
• 2013

In this paper, we consider the adaptive multigrid method for solving the Black-Scholes equation to improve the efficiency of the option pricing. Adaptive meshing is generally regarded as an indispensable tool because of reduction of the computational costs. The Black-Scholes equation is discretized using a Crank-Nicolson scheme on block-structured adaptively refined rectangular meshes. And the resulting discrete equations are solved by a fast solver such as a multigrid method. Numerical simulations are performed to confirm the efficiency of the adaptive multigrid technique. In particular, through the comparison of computational results on adaptively refined mesh and uniform mesh, we show that adaptively refined mesh solver is superior to a standard method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.11
• /
• pp.2932-2943
• /
• 1994

An efficient approach to the automatic construction of effective quadrilateral finite element meshes for two-dimensional analysis is presented. The procedure is composed of, firstly, an initial mesh generation and, secondly, an h-version of adaptive refinement based on error analysis. As for an initial mesh generation scheme, a modified looping algorithm has been employed. For the adaptive refinement process, an error indicator obtained by computing the residual error of the equilibrium equations in the energy norm with a relaxation factor has been employed. Examples of mesh generation and self-adaptive mesh improvements are given. These example solutions demonstrate that an effective mesh for a given error tolerance can be obtained in a few steps of the analysis processes.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.23 no.2
• /
• pp.27-32
• /
• 1986

This paper describes an application of the adaptive finite element computations to a free surface flow problem in a canal. A-posteriori error estimates for the adaptive finite element computations are based on the dual extremum principles. Previously the dual extremum principles were applied to compute the upper and lower bounds of the added mass of two-dimensional cylinders in a canal[1,2]. However, the present method improves the convergence of the computed results by utilizing the local error estimates and by applying the adaptive meshes in the finite element computations. In a test result using triangular elements it is shown that the numerical error in the adaptive finite elements reduces quadratically compared with that in a uniform mesh subdivision.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.8
• /
• pp.37-45
• /
• 2002

A two-dimensional unsteady, inviscid flow solver has been developed for the simulation of airfoil-vortex interactions on unstructured dynamically adapted meshes. The Euler solver is based on a second-order accurate implicit time integration using a point Gauss-Seidel relaxation scheme and a dual time-step subiteration. A vertex-centered, finite-volume discretization is used in conjunction with the Roe's flux-difference splitting. An unsteady solution-adaptive dynamic mesh scheme is used by adding and deleting mesh points to take account of both spatial and temporal variations of the flow field. The effect of vortex interaction on the aeroelastic displacement of an airfoil attached to the idealized two degree-of-freedom spring system is investigated.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.3
• /
• pp.139-170
• /
• 2013

We extend a p-hierarchical decomposition of the second degree finite element space of N$\acute{e}$d$\acute{e}$lec for tetrahedral meshes in three dimensions given in [1] to meshes with hexahedral elements, and derive p-hierarchical decompositions of the second degree finite element space of Raviart-Thomas in two dimensions for triangular and quadrilateral meshes. After having proved stability of these subspace decompositions and requiring certain saturation assumptions to hold, we construct a local a posteriori error estimator for fem and bem coupling of a time-harmonic electromagnetic eddy current problem in $\mathbb{R}^3$. We perform some numerical tests to underline reliability and efficiency of the estimator and test its usefulness in an adaptive refinement scheme.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1994.10a
• /
• pp.64-73
• /
• 1994

An updated Lagrangian finite element analysis with automatic remeshing scheme is applied to the three-dimensional hot extrusion through landless square dies. In the remeshing procedure, it is very difficult that the meshes are generated automatically with consideration of physical characteristics. In the presented study, the mesh generation is accomplished by modular concept. The generated meshes by modular concept have advantages, especially for three-dimensional problems, such as economized computational time and consideration of physical characteristic. In the problem, orifice shapes of square die are divided into two for the extrusion of solid sections. The orifice adaptive modules are developed for each type and the numerical examples are carried out for each type.