• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.5
• /
• pp.904-913
• /
• 2002

In this paper, the dynamic stress intensity factors of fracture mechanics are numerically computed in time domain using the FEM. For which the finite element formulations are derived applying the Galerkin method to the equations of motion in convolution integral as has been presented in the previous paper. To assure the strain fields of r$^{-1}$ 2/ singularity near the crack tip, the triangular quarter-point singular elements are imbedded in the finite element mesh discretized by the isoparametric quadratic quadrilateral elements. Two-dimensional problems of the elastodynamic fracture mechanics under the impact load are solved and compared with the existing numerical and analytical solutions, being shown that numerical results of good accuracy are obtained by the presented method.

• Journal of applied mathematics & informatics
• /
• v.31 no.1_2
• /
• pp.55-63
• /
• 2013

In this paper, a general technique is proposed for solving a system of fourth-order boundary value problems. The solution is given in the form of series and its approximate solution is obtained by truncating the series. Advantages of the method are that the representation of exact solution is obtained in a new reproducing kernel Hilbert space and accuracy of numerical computation is higher. Numerical results show that the method employed in the paper is valid. Numerical evidence is presented to show the applicability and superiority of the new method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.9
• /
• pp.1148-1156
• /
• 2000

The QUICK scheme for convection terms is modified for unstructured finite volume method by using linear reconstruction technique and validated through the computation of two well defined laminar flows. It uses two upstream grid points and one downstream grid point in approximating the convection terms. The most upstream grid point is generated by considering both the direction of flow and local grid line. Its value is calculated from surrounding grid points by using a linear construction method. Numerical error by the modified QUICK scheme is shown to decrease about 2.5 times faster than first order upwind scheme as grid size decreases. Computations are also carried out to see effects of the skewness and irregularity of grid on numerical solution. All numerical solutions show that the modified QUICK scheme is insensitive to both the skewness and irregularity of grid in terms of the accuracy of solution.

• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.9-15
• /
• 2007

A three dimensional numerical model to predict the flow and transport of mixtures and also the electrochemical reactions in polymer electrolyte membrane (PEM) fuel cells is developed. The numerical computation is base on vorticity- velocity method. Governing equations for the flow and transport of mixtures are coupled with the equations for electrochemical reactions and are solved simultaneously including production and condensation of vapor. Fuel cell performance predicted by this calculation is compared with the experimental results and resonable agreements are achieved.

• Computational Structural Engineering
• /
• v.6 no.2
• /
• pp.103-114
• /
• 1993

An advanced time-domain finite element formulation is presented for the displacement and stress analysis of isotropic, linear viscoelastic problems of a hereditary-type constitutive law. The semidiscrete finite element method with linear time-stepping scheme and an elastic-viscoelastic correspondence principle are used in the theoretical development. An efficient treatment of hereditary integral is introduced to improve the numerical accuracy, to reduce the computation time, and to avoid the use of large memory storage. Two-dimensional numerical examples of plane strain and plane stress are solved under the assumption on the material property of being elastic in dilatation and like three-element Voigt model in distorsion, and compared with the analytical solutions and the past numerical results to show the versatility and efficiency of the proposed method.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.04a
• /
• pp.1-8
• /
• 2004

For the linearized differential algebraic equation of the nonlinear constrained system, exact initial values of the acceleration are needed to solve itself. It may be very troublesome to perform the inverse operation for obtaining the incremental quantities since the mass matrix contains the zero element in the diagonal. This fact makes the mass matrix impossible to be positive definite. To overcome this singularity phenomenon the mass matrix needs to be modified to allow the feasible application of predictor and corrector in the iterative computation. In this paper the proposed numerical algorithm based on the modified mass matrix combines the conventional implicit algorithm, Newton-Raphson method and Newmark method. The numerical example presents reliabilities for the proposed algorithm via comparisons of the 4th order Runge-kutta method. The proposed algorithm seems to be satisfactory even though the acceleration, Lagrange multiplier, and energy show unstable behaviour. Correspondingly, it provides one important clue to another algorithm for the enhancement of the numerical results.

• Wind and Structures
• /
• v.9 no.6
• /
• pp.473-492
• /
• 2006

This paper presents numerical results of pedestrian-level wind environment around the base of a row of tall buildings by CFD. Four configurations of building arrangement are computed including a single square tall building. Computed results of pedestrian-level wind flow patterns and wind speeds are compared to previous wind tunnel measurement data to enable an assessment of CFD predictions. The CFD model uses the finite-volume method with RNG $k-{\varepsilon}$ model for turbulence closure. It is found that the numerical results can reproduce key features of pedestrian-level wind environment such as corner streams around corners of upwind building, sheltered zones behind buildings and channeled high-speed flow through a building gap. However, there are some differences between CFD results and wind tunnel data in the wind speed distribution and locations of highest wind speeds inside the corner streams. In locations of high ground-level wind speeds, CFD values match wind tunnel data within ${\pm}10%$.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.148-153
• /
• 2001

A computational investigation of the effect of the electromagnetic force(or Lorentz force) on the flow behavior around a circular cylinder, a typical model of bluff bodies, is conducted. Two-dimensional unsteady flow computation for $Re=10^2$ is carried out using a numerical method of finite difference approximation in a curvilinear body-fitted coordinate system by solving the momentum equations including the Lorentz force as a body force. The effect of the spatial variations of the Lorentz forcing region and forcing direction along the cylinder circumference is investigated. The numerical results show that the Lorentz force can effectively suppress the flow separation and oscillation of the lift force of the circular cylinder cross-flow, leading to the reduction of the drag.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.3
• /
• pp.359-368
• /
• 2010

This paper describes the generation and propagation of internal solitary wave in a two-layer fluid system by numerical analysis. Characteristics of interfacial soliton such as wave type, wave height, wave celerity are investigated numerically with respect to an extent of initial disturbance, fluid thicknesses of the two fluids and etc. The difference between the internal wave propagation on sloping beach and flat bottom was also examined. Laboratory experiments were conducted in the wave flume and compared with the results of numerical computation for verification.

• Proceedings of the KSME Conference
• /
• 2001.06c
• /
• pp.604-609
• /
• 2001

In the sheet metal forming process, the drawbead is used to control the flow of material during the forming process. The drawbead provides proper restraining force to the material and prevents defects such as wrinkling or breakage. For these reasons, many studies for designing the effective drawbead have been conducted. For the analysis, the numerical method called the static-explicit finite element method was used. The finite element analysis code for this method has been developed and applied to the drawbead process problems. In result, convergence problem and computation time due to large non-linearity in the existing numerical analysis methods were no longer a critical problem. Futhermore, this approach could treat the contact friction problem easily by applying very small time intervals. It is expected that various results from the numerical analysis will give very useful information for the design of tools in sheet metal forming process.