• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.13 no.2
• /
• pp.10-15
• /
• 1984

The aims of this study are to obtain a suitable method and a proper mesh for investigation of the temperature distribution and heat transfer. The relative errors of temperature distribution and heat transfer for each mesh are acquired in accordance with linear finite element (FEM 3), square finite element (FEM 6), cubic finite element (FEM 10), and finite difference method (FDM). It has been found that FEM 10 is the most accurate measure to obtain the temperature distribution and heat transfer. However, no significant results have been obtained successfully, because when higher order finite element methods are used the more computational efforts are necessary due to the distribution of elements. The results of this study are as follows ; 1 . In case of a=b=L, meshes for less than $1\%$ relative errors (temperature distribution) acquired in various methods to exact solution are $2\times2,\;4\times4,\;8\times8\;and\;8\tiems8$ for each FEM 10, FEM 6, FEM 3 and FDM and a=L, b=1/2L are $10\times5$ for each FEM 3 and FDM. And the tendency of results acquired of heat transfer is similar to those above. 2 . In computational efforts (a=b=L), FEM 6 has taken 21 times. and FEM 10 154times FEM 3 and FDM and FEM 3 is the sane as FDM.

• Journal of Korean Society of Steel Construction
• /
• v.11 no.4 s.41
• /
• pp.373-384
• /
• 1999

The shell structures with composite materials have the advantages in strength, corrosion resistance, and weight reduction. The objective of this study is to analyze anisotropic composite circular cylindrical shells with shear deformation theory. In applying numerical methods to solve differential equations of anisotropic shells, this paper use finite difference method. The accuracy of the numerical method can be improved by taking higher order of interval ${\Delta}$ to reduce error. This study compares the results of finite difference method with the results of ANSYS based on finite element method. Several numerical examples show the advantages of the stiffness increasement when the composite materials aroused. Therefore, it is expected that results of this study give various guides for change of the subtended angles, load cases, boundary conditions, and side-to-thickness ratio.

• Proceedings of the KSME Conference
• /
• 2001.06a
• /
• pp.686-691
• /
• 2001

Structural optimization often require the evaluation of design sensitivities. The Semi Analytic method(SAM) is popular for shape optimization because this method has several advantages. But when relatively large rigid body motions are identified for individual elements, the SA method shows severe inaccuracy. In this paper, the improvement of design sensitivities corresponding to the rigid body mode is evaluated by exact differentiation of the rigid body modes. Moreover, the error of the SA method caused by numerical difference scheme is alleviated by using a series approximation for the sensitivity derivatives and considering the higher order terms.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.29 no.3
• /
• pp.225-232
• /
• 2018

This paper proposed the optimization method of the extremely low numerical-dispersion finite-difference time-domain (ELND-FDTD) method based on the H(2,4) scheme for wideband and extremely accurate electromagnetic properties of lossy material, which has a constant conductivity and relative permittivity. The optimized values of three variables are calculated for the minimum numerical dispersion errors of the proposed FDTD method. The excellent accuracy of the proposed method is verified by comparing the calculated results of three different FDTD methods and the analytical results of the two-dimensional dielectric cylinder scattering problem.

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.26-34
• /
• 2001

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Navier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to 25% of the cylinder diameter and in the case of the lock-in region it is 60%. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.181-188
• /
• 2001

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Wavier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to $25\%$ of the cylinder diameter and in the case of the lock-in region it is $60\%$. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.255-262
• /
• 2013

In this paper, a shape design sensitivity analysis(DSA) method is presented for Mindlin plates using an isogeometric approach. The isogeometric method possesses desirable advantages; the representation of exact geometry and the higher order inter-element continuity, which lead to the fast convergence of solution as well as accurate sensitivity results. Unlike the finite element methods using linear shape functions, the isogeometric method considers the exact normal vector and curvature of the CAD geometry, taking advantages of higher order NURBS basis functions. A selective reduced integration(SRI) technique is incorporated to overcome the difficulty of 'shear locking' phenomenon. This simple technique is surprisingly helpful for the accuracy of the isogeometric shape sensitivity without complicated formulation. Through the numerical examples of plate bending problems, the accuracy of the proposed isogeometric analysis method is compared with that of finite element one. Also, the isogeometric shape sensitivity turns out to be very accurate when compared with finite difference sensitivity.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.1
• /
• pp.8-15
• /
• 2002

A two-dimensional higher order panel method using B-splines has been developed to overcome the disadvantages of the low order panel method and to obtain more accurate solution. The sources and the normal dipoles are distributed on both the body and the free surface. Instead of applying the upwind finite difference schemes to satisfy the linearized free surface and the radiation condition, the derivatives of the basis functions of the B-splines are directly applied to the linearized free surface condition. Numerical damping in the Dawson's method are avoided in the Present computations. In order to validate the present method, numerical computations are carried out for a submerged cylinder and a two-dimensional hydrofoil steadily moving beneath a free surface. The numerical results show that fast convergence and better accuracies have been achieved by the present method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.11
• /
• pp.3031-3040
• /
• 1995

Two-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with the unstructured triangular meshes. The pressure-velocity coupling is handled by the artificial compressibility algorithm due to its computational efficiency associated with the hyperbolic nature of the resulting equations. The convective fluxes are obtained by the Roe's flux difference splitting scheme using edge-based connectivities and higher-order differences are achieved by a reconstruction procedure. The time integration is based on an explicit four-stage Runge-Kutta scheme. Numerical procedures with local time stepping and implicit residual smoothing have been implemented to accelerate the convergence for the steady-state solutions. Comparisons with experimental data and other numerical results have proven accuracy and efficiency of the present unstructured approach.

• Journal of Advanced Marine Engineering and Technology
• /
• v.40 no.4
• /
• pp.264-269
• /
• 2016

This paper reports a fundamental study of temperature characteristics of a solar pond with an insulation layer. Further, these characteristics were compared with those of a solar pond without the insulation layer. The governing equation was discretized via finite difference method. The governing equations are two-dimensional unsteady-state second-order partial differential equations. The conclusions of the study are as follows: 1) If the depth of the solar pond was increased, the desired effect of increase in temperature was not produced because the amount of solar insolation received by the bottom of the solar pond decreased. 2) As the temperature of the soil during winter is higher than the temperature of the water in a solar pond, heat was transferred from the soil to the solar pond. 3) For the case of the solar pond with insulation layer, it was estimated that the dependence rate of solar energy was 83.3% and that of the boiler was 16.7%.