• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2901-2906
• /
• 2007

A solver for icing and condensation of water has been developed. The phase change process was solved by the enthalpy method. For the code validation, the temperature and the phase change from water to ice of the driven cavity were calculated. Also, the melting process of the frost on the windshield glass of an automobile has been simulated. The calculation showed a good agreement with analytical solution and other numerical results. Using the present validated code, the condensation of water vapor has been first tried. The computed results provided some physical features of condensation phenomena even though experimental data and other numerical data were not available. For future work, it is recommended to throughly investigate the effects of boundary conditions on the solution.

• 한국연소학회:학술대회논문집
• /
• 2004.11a
• /
• pp.109-116
• /
• 2004

Through analysing the influence of steam flow direction on the liquid formation and motion behavior in the condenser shell side, the physical model for existing numerical simulation program of condenser is improved by introducing the correlations for flow resistance and condensation heat exchange coefficient in which the influences of steam flow direction are considered according to the available experimental data. Thus a more suitable and general condenser simulation approach is presented and a new condenser calculation program is developed. With the experimental data of a pannier- arrangement experimental condenser, the adaptability of the new condenser simulation approach is verified. General characteristics of this type of condenser are also revealed.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.10 no.5
• /
• pp.266-272
• /
• 1985

The analysis of a cylindrical waveguide which is partially filled with dilectric materials has been performed for many years. Hut most of the analyses were an approximation by the analytic method. In this paper a FORTRAN program for numerical analysis is introduced to calculate the propagation constant for TE. TM and Hybrid modes. The results of the numerical calculation can be used to determine the resonant frequencies for TM. TE and Hybrid modes in the dielectric resonator.

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

• Structural Engineering and Mechanics
• /
• v.66 no.5
• /
• pp.649-663
• /
• 2018

The scaled boundary finite element method is extended to evaluate the dynamic stress intensity factors and T-stress with a numerical procedure based on the improved continued-fraction. The improved continued-fraction approach for the dynamic stiffness matrix is introduced to represent the inertial effect at high frequencies, which leads to numerically better conditioned matrices. After separating the singular stress term from other high order terms, the internal displacements can be obtained by numerical integration and no mesh refinement is needed around the crack tip. The condition numbers of coefficient matrix of the improved method are much smaller than that of the original method, which shows that the improved algorithm can obtain well-conditioned coefficient matrices, and the efficiency of the solution process and its stability can be significantly improved. Several numerical examples are presented to demonstrate the increased robustness and efficiency of the proposed method in both homogeneous and bimaterial crack problems.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.3
• /
• pp.624-634
• /
• 1994

A numerical method which combines equal-order velocity-pressure formulation originated from SIMPLE algorithm and streamline upwinding method has been developed. To verify the proposed numerical method, we considered the lid-driven cavity flow and backward facing step flow. The trend of convergence history is stable up to the error criterion beyond which the maximum value of error is oscillatory due4 to the round-off error. In the present study, all results were obtained with the single precision calculation up to the given error criterion and it was found to be sufficient for our purpose. The present results were then compared with existing experimental results using laser doppler velocimetry and numerical results using finite difference method and mixed interpolation finite element method. It has been shown that the present method gives accurate results with less memories and execution time than the coventional finite element method.

• Journal of Advanced Marine Engineering and Technology
• /
• v.29 no.1
• /
• pp.78-85
• /
• 2005

A numerical calculation is performed to study the effects of buoyancy force on the heat transfer characteristics of laminar forced convection flow in inclined parallel plates with the upper part cooled and the lower Part heated uniformly. Numerical results are presented for the Reynolds number ranges from $4.0\times10^{-3}$ to $1.13\times10^{-1}$. the angle of inclination, $\theta$. from 0 to 90 degree and Pr of the high viscosity fluid is 909. It is found that the flow pattern of mixed convection in inclined parallel Plates can be classified into four patterns which affected by Reynolds number and the angle of inclination.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.10 s.241
• /
• pp.1083-1091
• /
• 2005

Turbulent flow around a rotating circular cylinder is investigated by Direct Numerical Simulation. The calculation is performed at three cases of low Reynolds number, Re=161, 348 and 623, based on the cylinder radius and friction velocity. Statistically strong similarities with fully developed channel flow are observed. Instantaneous flow visualization reveals that the turbulence length scale typically decreases as Reynolds number increases. Some insight into the spacial characteristics in conjunction with wave number is provided by wavelet analysis. The budget of dissipation rate as well as turbulent kinetic energy is computed and particular attention is given to the comparison with plane channel flow.

• Tunnel and Underground Space
• /
• v.7 no.4
• /
• pp.334-340
• /
• 1997

Alpha center theory which was first proposed by Stefanescu has been proved to be effective for the detection of the location of the conductive orebody. A numerical forward modeling was conducted to verify the effectiveness of this method. Field works were carried out along the three profiles in two different areas for the purpose of finding fractured zone which might be accompanied with the presence of the groundwater. And the results were modeled by alpha center method, which was later testified by wellproven 2-dimensional finite difference inversion scheme. Field data could be successfully modeled with this alpha center algorithm, especially for the smooth-varying resistivity models. For the abrupt change of the resistivity values, the alpha center coefficients have a tendency to be negative to simulate the steep resistivity gradients. This method is quite simple and easy for the future applications. The numerical calculation can be performed very quickly with the personal computers.

• Journal of Mechanical Science and Technology
• /
• v.16 no.10
• /
• pp.1327-1335
• /
• 2002

The shock wave process represents an abrupt change in fluid properties, in which finite variations in pressure, temperature, and density occur over the shock thickness which is comparable to the mean free path of the gas molecules involved. This shock wave fluid phenomenon is simulated by using the finite difference lattice Boltzmann method (FDLBM). In this paper, a new model is proposed using the lattice BGK compressible fluid model in FDLBM for the purpose of speeding up the calculation as well as stabilizing the numerical scheme. The numerical results of the proposed model show good agreement with the theoretical predictions.