• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.6
• /
• pp.916-922
• /
• 2001

A mechanical face seal is a tribo-element intended to control leakage of working fluid at the interface between a rotating shaft and its housing. Leakage of working fluid decreases drastically as the clearance of the mating seal faces gets smaller. But the very small seal clearance results in an increased reduction of seal life because of high wear and heat generation. Therefore, in the design of mechanical face seals a compromise between low leakage and acceptable seal life is important, and it presents a difficult and practical design problem. A fluid film or sealing dam geometry of the seal clearance affects seal lubrication performance very much, and thereby it is one of the main design considerations. In this study the Reynolds equation for the sealing dam of mechanical face seals is numerically analyzed, using the Galerkin finite element method, which is readily applied to various seal geometries. Film pressures of the sealing dam are analyzed, including the effects of the seal face coning and tilt. Then, lubrication performances of the seals, such as opening forces, restoring moments, leakage, and dynamic coefficients, are calculated, and they are compared to the results obtained by the narrow seal approximation.

• Journal of electromagnetic engineering and science
• /
• v.19 no.1
• /
• pp.6-12
• /
• 2019

An IE-FFT algorithm is implemented and applied to the electromagnetic (EM) solution of perfect electric conducting (PEC) scattering problems. The solution of the method of moments (MoM), based on the magnetic field integral equation (MFIE), is obtained for PEC objects with closed surfaces. The IE-FFT algorithm uses a uniform Cartesian grid to apply a global fast Fourier transform (FFT), which leads to significantly reduce memory requirement and speed up CPU with an iterative solver. The IE-FFT algorithm utilizes two discretizations, one for the unknown induced surface current on the planar triangular patches of 3D arbitrary geometries and the other on a uniform Cartesian grid for interpolating the free-space Green's function. The uniform interpolation of the Green's functions allows for a global FFT for far-field interaction terms, and the near-field interaction terms should be adequately corrected. A 3D block-Toeplitz structure for the Lagrangian interpolation of the Green's function is proposed. The MFIE formulation with the IE-FFT algorithm, without the help of a preconditioner, is converged in certain iterations with a generalized minimal residual (GMRES) method. The complexity of the IE-FFT is found to be approximately $O(N^{1.5})$and $O(N^{1.5}logN)$ for memory requirements and CPU time, respectively.

• Journal of Korean Association for Spatial Structures
• /
• v.21 no.4
• /
• pp.65-72
• /
• 2021

The Precast concrete(PC) modular structures are a method of assembling pre-fabricated unit modules in the construction site. The essential aim of modular structures is to introduce a connection method that can ensure splicing performance and effectively resist shear strength. This study proposed PC module using a connecting plate that can replace splice sleeves and shear keys used in the conventional PC modular structures. To evaluate the splicing performance and shear capacity of the proposed method, the shear test was conducted by fabricating one monolithic reinforced concrete(RC) beam and two PC modular beams with a shear span-to-depth ratio as variables. The experimental results showed that the shear capacity of the PC modular beam was about 89% compared to that of the RC beam, and showed a failure of the RC beam according to the shear span-to-depth ratio. Therefore, it was considered that the connecting plate effectively transferred the stress between each PC module through the joint and ensure integrity. In addition, the applicability of shear strength equation of ACI 318-19 and Zsutty's equation to PC modular beams were evaluated. Results demonstrated that the improved shear strength equations are needed to consider reduction of shear strength in PC modules.

• Advances in Computational Design
• /
• v.9 no.1
• /
• pp.39-52
• /
• 2024

If the governing differential equation arising from engineering problems is treated as an analytic, continuous and derivable function, it can be expanded by one point as a series of finite numbers. For the function to be zero for each value of its domain, the coefficients of each term of the same power must be zero. This results in a recursive relationship which, after applying the natural conditions or the boundary conditions, makes it possible to obtain the values of the derivatives of the function with acceptable accuracy. The elastoplastic analysis of an inhomogeneous thick sphere of metallic materials with linear variation of the modulus of elasticity, yield stress and Poisson's ratio as a function of radius subjected to internal pressure is presented. The Beltrami-Michell equation is established by combining equilibrium, compatibility and constitutive equations. Assuming axisymmetric conditions, the spherical coordinate parameters can be used as principal stress axes. Since there is no analytical solution, the natural boundary conditions are applied and the governing equations are solved using a proposed new method. The maximum effective stress of the von Mises yield criterion occurs at the inner surface; therefore, the negative sign of the linear yield stress gradation parameter should be considered to calculate the optimal yield pressure. The numerical examples are performed and the plots of the numerical results are presented. The validation of the numerical results is observed by modeling the elastoplastic heterogeneous thick sphere as a pressurized multilayer composite reservoir in Abaqus software. The subroutine USDFLD was additionally written to model the continuous gradation of the material.

• Journal of Korean Society for Atmospheric Environment
• /
• v.9 no.E
• /
• pp.347-357
• /
• 1993

The Lagrangian grid-free numerical method, the discrete vortex method, was applied to solve the Navier-Stokes euqations. This method avoids the introduction of numerical viscosity swamping the real physical viscosity at high Reynolds number, unlike Eulerian method, e.g. finite difference and element methods. The boundary integral equation method for the potential flow solution was included to make the discrete vortex method more feasible for complex geometries. The fast adaptive multipole expansion method was incorporated to reduce the computational time from $O(N^2)$ to O(N) for the computations of vortex-vortex interactions. The test problems were air flow around one circular cylinder and two circular cylinders in tandem with various gaps. The numerical results were in excellent gareement with the experimental and other computational results. The applicabilty of the method was discussed with the indoor and the outdoor air pollution problems, especially the contaminant transport in the recirculation regions.

• The Journal of the Acoustical Society of Korea
• /
• v.39 no.1
• /
• pp.1-7
• /
• 2020

The acoustic parabolic equation method in the ocean is an efficient technique to calculate the acoustic field in the range-dependent environment, emanating from a point source. However, we often need to use the directional source with a main beam in the practical problem. In this paper, we present two methods to implement the directional source in the acoustic parabolic equation code easily. One is simply to filter the Delta function idealized as an omni-directional point source. Another method is based on the rational filtering of the self-starter solution. It has a limitation not to separate the up-going and the down-going wave for the depth, but would be useful in implementing the mode propagation. Numerical examples for validation are given in the Pekeris environment and the deep sea environment.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.1
• /
• pp.65-70
• /
• 2014

This paper introduces a preform design method for a ring rolling process with an outer step. Underfilling is one of the general defects of the profile ring rolling process. It occurs when the deformation amount is small or step depth of the profiled region is large. To prevent underfilling, increasing the deformation amount or using a preform of size similar to that of the final product are required. Furthermore, the filling limit equation is suggested based on the shape factor and deformation ratio for preventing defects in the products. The filling limit equation has been derived through finite element analyses and production tests for four different cases. For verifying the suggested method, realsized profile rolling tests were performed, and test results were compared with the predictions of the equation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.2
• /
• pp.194-200
• /
• 2002

A method for simulating flow of extremely low Reynolds number in a vortex-in-cell (VIC) was studied. The viscous diffusion was represented by the random walk method. The validity of this method was proven by applying it to the flows passing over a fence placed vertically at an extremely low Reynolds number. The effects of parameters of the VIC method on the random walk method, such as the number of point vortices, the mesh density and the time increment, were investigated by numerical analysis using a one-dimensional diffusion equation. Changes in the relative error of vorticity depending on those parameters were clarified.

• Journal of Korean Society of Steel Construction
• /
• v.23 no.3
• /
• pp.327-335
• /
• 2011

This study has two objectives: (1) to evaluate the degradation characteristics of symmetric unbraced steel frames by using analytical approach, and (2) to suggest equation which can approximately estimate the effect of degradation during the schematic design stage. For the analytical approach, the refined plastic hinge method with an arc length algorithm was adopted. The subject of analysis was one story one-bay, multistory one-bay, and multistory three-bay unbraced steel frames. The main parameters of the analytical approach include the stiffness ratio of column to beam and the axial force ratio. The study led to the following conclusions. The normalized stiffness of degradations is affected by both stiffness ratio of column to beam and the axial load ratio; however, the major influence on degradations is the axial force ratio. The equation, which can approximately estimate the effect of degradation, was suggested together with the research results.

• Journal of Korea Water Resources Association
• /
• v.55 no.10
• /
• pp.775-784
• /
• 2022

In this study, to determine the optimal order of the full-logged I-D-F polynomial equation, which is mainly used to calculate the probable rainfall over a temporal rainfall duration, the probable rainfall was calculated and the regression coefficients of the full-logged I-D-F polynomial equation was estimated. The optimal variable of the polynomial equation for each station was selected using a stepwise selection method, and statistical significance tests were performed through ANOVA. Using these results, the statistically appropriately calculated rainfall intensity equation for each station was presented. As a result of analyzing the variable selection outputs of the full-logged I-D-F polynomial equation at 9 stations in Gyeongbuk, the 1^st to 3^rd order equations at 6 stations and the incomplete 3^rd order at 1 station were determined as the optimal equations. Since the 1^st order equation is similar to the Sherman type equation and the 2^nd order one is similar to the general type equation, it was presented as a unified form of rainfall intensity equation for convenience of use by increasing the number of independent variables. Therefore, it is judged that there is no statistical problem in considering only the 3^rd order polynomial regression equation for the full-logged I-D-F.