• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.5-9
• /
• 2010

A numerical wind tunnel simulation is performed in order to predict wind loads acting on a building. The aim of the present study is to suggest a guideline for the numerical wind tunnel analysis, which could provide more detail wind load distributions compared to the wind code and expensive wind tunnel experiments. To validate the present numerical simulation, wind-induced loads on a 6 m cube model is predicted. Atmospheric boundary layer is used as a inlet boundary condition. Various effect of numerical methods are investigated such as size of computational domain, grid density, turbulence model and discretization scheme. The appropriate procedure for the numerical wind tunnel analysis is suggested through the present study.

• Structural Engineering and Mechanics
• /
• v.21 no.3
• /
• pp.315-332
• /
• 2005

The Element-free Galerkin Method has become a very popular tool for the simulation of mechanical problems with moving boundaries. The internally applied Moving Least Squares interpolation uses in general Gaussian or cubic weighting functions and has compact support. Due to the approximative character of this interpolation the obtained shape functions do not fulfill the interpolation conditions, which causes additional numerical effort for the application of the boundary conditions. In this paper a new weighting function is presented, which was designed for meshless shape functions to fulfill these essential conditions with very high accuracy without any additional effort. Furthermore this interpolation gives much more stable results for varying size of the influence radius and for strongly distorted nodal arrangements than existing weighting function types.

• Proceedings of the Korean Society of Machine Tool Engineers Conference
• /
• 2002.10a
• /
• pp.3-10
• /
• 2002

This study proposed a technique for inverse problem, linear approximation of contact position and loading in single and double meshing of spur gear, using 2-dimension model considered near the tooth by root stress. Determine root stress is carried out for the gear tooth by finite element method and boundary element method. Boundary element discretization near contact point is carefully performed to keep high computational accuracy. And from those estimated results, the comparing estimate value with boundary element method value was discussed.

• Korean Journal of Mathematics
• /
• v.7 no.2
• /
• pp.245-269
• /
• 1999

Multigrid method and acceleration by conjugate gradient method for first-order system least squares (FOSLS) using bilinear finite elements are developed for various boundary value problems of planar linear elasticity. They are two-stage algorithms that first solve for the displacement flux variable, then for the displacement itself. This paper focuses on solving for the displacement flux variable only. Numerical results show that the convergence is uniform even as the material becomes nearly incompressible. Computations for convergence factors and discretization errors are included. Heuristic arguments to improve the convergences are discussed as well.

• Structural Engineering and Mechanics
• /
• v.25 no.2
• /
• pp.161-180
• /
• 2007

A first endeavor is made to exploit the differential quadrature method (DQM) as a simple, accurate, and computationally efficient numerical tool for the large deformation analysis of thin laminated composite skew plates, which has very strong singularity at the obtuse vertex. The geometrical nonlinearity is modeled by using Green's strain and von Karman assumption. A recently developed DQ methodology is used to exactly implement the multiple boundary conditions at the edges of skew plates, which is a major draw back of conventional DQM. Using oblique coordinate system and the DQ methodology, a mapping-DQ discretization rule is developed to simultaneously transform and discretize the equilibrium equations and the related boundary conditions. The effects of skew angle, aspect ratio and different types of boundary conditions on the convergence and accuracy of the presented method are studied. Comparing the results with the available results from other numerical or analytical methods, it is shown that accurate results are obtained even when using only small number of grid points. Finally, numerical results for large deflection behavior of antisymmetric cross ply skew plates with different geometrical parameters and boundary conditions are presented.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.1
• /
• pp.10-17
• /
• 1985

A new method is suggested for the solution of plane elasticity problems. With use of the dislocation model in the crack problems, the basic scheme of this method is to find equilibrium Burgers vectors of dislocations which are distributed along the boundary of the first fundamental boundary value problems. The stress distribution in the region can be found by superposition of the contributions of each dislocation. The method is applied to three cases with known analytical solutions, and to a V-notched specimen under uniaxial tension. The numerical results are compared with other available solutions. This method is effective and simple in its use, compared with other numerical methods. The method also provides very accurate solutions in the region except near the boundary where the discretization error is significant. The extrapolation method is suggested for the stresses in the boundary region. Extensive application are also suggested for a general estimate of the computational efficiency of the method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.5
• /
• pp.787-797
• /
• 1986

An implicit approach is employed to obtain a general boundary integral formulation of axisymmetric elastic problems in terms of a pair of singular integral equations. The corresponding kernel functions from the solutions of Navier's equation are derived by applying a three dimensional integral and a direct axisymmetrical approach. A numerical discretization schem including the evaluation of Cauchy principal values of the singular integral is described. Finally the typical axisymmetric elastic models are analyzed, i.e. the hollow sphere, the constant thickness and the V-notched round bar.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.3 s.258
• /
• pp.300-305
• /
• 2007

Many practical flow problems are defined with the circular boundary. Fluid flows within a circular boundary are however susceptible to a singularity problem when the cylindrical coordinates are employed. To remove this singularity a method has been developed in this study which uses the local harmonic functions in discretization of derivatives as well as interpolation. This paper describes the basic reason for introducing the harmonic functions and the overall numerical methods. The numerical methods are evaluated in terms of the accuracy and the stability. The Lamb-dipole flow is selected as a test flow. We will see that the harmonic-function method indeed gives more accurate solutions than the conventional methods in which the polynomial functions are utilized.

• Journal of Hydrospace Technology
• /
• v.2 no.2
• /
• pp.36-43
• /
• 1996

In this work, a panel method is described, which cart solve the flow field round a surface-piercing body that experiences lift and wave resistance. As the body boundary condition, a Dirichlet type is employed, and as the free surface boundary condition the Poisson type is implemented, while in its discretization Dawson's 4-point upwind difference scheme is utilized, and as the Kutta condition a Morino-Kuo type is chosen. As to the type of singularity, source panels are distributed on the free surface, and source and dipole panels on the body surface, and dipole panels on the wake surface. For a sample run, a catamaran of the parabolic Wigley hull is chosen, for which experimental data are available, and the predictions by the numerical means and by the experiment are compared for a wide range of parameters.

• Nuclear Engineering and Technology
• /
• v.51 no.1
• /
• pp.13-30
• /
• 2019

This paper presents the recent improvements in the DeCART code for HTGR analysis. A new 190-group DeCART cross-section library based on ENDF/B-VII.0 was generated using the KAERI library processing system for HTGR. Two methods for the eigen-mode adjoint flux calculation were implemented. An azimuthal angle discretization method based on the Gaussian quadrature was implemented to reduce the error from the azimuthal angle discretization. A two-level parallelization using MPI and OpenMP was adopted for massive parallel computations. A quadratic depletion solver was implemented to reduce the error involved in the Gd depletion. A module to generate equivalent group constants was implemented for the nodal codes. The capabilities of the DeCART code were improved for geometry handling including an approximate treatment of a cylindrical outer boundary, an explicit border model, the R-G-B checker-board model, and a super-cell model for a hexagonal geometry. The newly improved and implemented functionalities were verified against various numerical benchmarks such as OECD/MHTGR-350 benchmark phase III problems, two-dimensional high temperature gas cooled reactor benchmark problems derived from the MHTGR-350 reference design, and numerical benchmark problems based on the compact nuclear power source experiment by comparing the DeCART solutions with the Monte-Carlo reference solutions obtained using the McCARD code.