• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.6
• /
• pp.405-411
• /
• 2013

A sediment transport calculation considering cohesive force is proposed to deal with the transport phenomena of cohesive sediment. In the proposed calculation, each sand particle is assumed to be surrounded by a thin layer of mud. The critical Shields parameter and bed-load sediment transport rate are modified to include the cohesive force acting on the sand particle. The proposed calculation is incorporated into a two-way coupled fluid-structure-sediment interaction model, and applied to wave-induced topographic change of artificial shallows. Numerical results show that an increase in the content ratio of the mud, cohesive resistance force per unit surface area and water content cause increases in the critical Shields parameter and decreases in the bed-load sediment transport rate, reducing the topographic change of the shallow without changing its trend. This suggests that mixing mud in the pores of the sand particles can reduce the topographic change of shallows.

• KIEE International Transactions on Electrophysics and Applications
• /
• v.4C no.5
• /
• pp.220-229
• /
• 2004

We have examined a technique of one-dimensional (1D) fluid modeling for radio-frequency Ar capacitively coupled plasmas (CCP) between unequal-sized powered and grounded electrodes. In order to simulate a practical CCP reactor configuration with a grounded side wall by the 1D model, it has been assumed that the discharge space has a conic frustum shape; the grounded electrode is larger than the powered one and the discharge space expands with the distance from the powered electrode. In this paper, we focus on how much a 1D model can approximate a 2D model and evaluate their comparisons. The plasma density calculated by the 1D model has been compared with that by a two-dimensional (2D) fluid model, and a qualitative agreement between them has been obtained. In addition, 1D and 2D calculation results for another reactor configuration with equal-sized electrodes have also been presented together for comparison. In the discussion, four CCP models, which are 1D and 2D models with symmetric and asymmetric geometries, are compared with each other and the DC self-bias voltage has been focused on as a characteristic property that reflects the unequal electrode surface areas. Reactor configuration and experimental parameters, which the self-bias depends on, have been investigated to develop the ID modeling for reactor geometry with unequal-sized electrodes.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.292-298
• /
• 1991

Previous studies on the performance of the round type impactor were reexamined and extended to the case of high particle mass loading. It was pointed out that the previous numerical studies need to be supplemented in the numerical process. The impactor performance was calculated under the same conditions as previous studies by the exact calculation process and it was found out that a tail of the collection efficiency curve, which have not been found in the previous studies, appeared in the results of ours. Numerical results for high particle mass loading show that the value of the collection efficiency in the impactor decreases but better particle-cut characteristics can be obtained, as the amount of the particle mass loading increases.

• Proceedings of the KIEE Conference
• /
• 2009.07a
• /
• pp.847_848
• /
• 2009

This paper shows the temperature rise of the high voltage GIS bus bar. The temperature rise in GIS bus bar is due to Joule‘s losses in the conductor and the induced eddy current in the tank. The power losses of a bus bar calculated from the magnetic field analysis are used as the input data for the thermal analysis to predict the temperature. The required analysis is a couple-field Multiphysics that accounts for the interactions between three-dimensional AC harmonic magnetic and fluid fields. The heat transfer calculation using the fluid analysis is done by considering the natural convection and the radiation from the tank to the atmosphere. Consequently, because temperature distributions by couple-field Multiphysics (coupled magnetic-fluid) have good agreement with results of temperature rise test, the proposed couple-field Multiphysics technique is likely to be used in a conduction design of the single-pole and three pole-encapsulated bus bar in GIS..

• Wind and Structures
• /
• v.23 no.6
• /
• pp.523-542
• /
• 2016

Snowdrift formation on roofs should be considered in snowy and windy areas to ensure the safety of buildings. Presently, the prediction of snowdrifts on roofs relies heavily on field measurements, wind tunnel tests and numerical simulations. In this paper, a new snowdrift modeling method by using CFD (Computational Fluid Dynamics) coupled with DEM (Discrete Element Method) is presented, including material parameters and particle size, collision parameters, particle numbers and input modes, boundary conditions of CFD, simulation time and inlet velocity, and coupling calculation process. Not only is the two-way coupling between wind and snow particles which includes the transient changes in snow surface topography, but also the cohesion and collision between snow particles are taken into account. The numerical method is applied to simulate the snowdrift on a typical stepped flat roof. The feasibility of using coupled CFD with DEM to study snowdrift is verified by comparing the simulation results with field measurement results on the snow depth distribution of the lower roof.

• Journal of Korea Foundry Society
• /
• v.22 no.4
• /
• pp.192-199
• /
• 2002

Dimensional defects of castings are mainly due to the stresses and strains caused by a nonuniform temperature distribution and phase transformation during solidification and cooling, and by mechanical constraint between the mold and casting. It is, however, nearly impossible to trace movements of the casting and mold during solidification and cooling by experimental measurements for castings with complex shape. Two and three dimensional deformation analyses of the casting and the mold were performed using commercial finite element code, MARC. It was possible to calculate deformation and temperature distribution in the casting and mold simultaneously. Cooling curves of the casting obtained by calculation were close to that measured in the field since it was possible to treat latent heat evolution of the casting which could be divided into two parts, primary and eutectic parts. Mold bent inward just after pouring due to the temperature gradient across the mold thickness, and mold returned to its previous position with time. Plastic deformation occurred at the part of the casting where solidification was slow.

• Selected Papers of The Society of Naval Architects of Korea
• /
• v.2 no.1
• /
• pp.18-28
• /
• 1994

In the vibration analysis of submerged or floating bodies such as ships and offshore structures, the coupled system between fluid and structure should be considered using the compatibility conditions on the wetted surface. It is well known that the hydroelastic vibration analysis of structures in contact with fluid can be done by applying the finite element method(FEM) to structures and the boundary element method(BEM) to the fluid domain. However, such an approach is impractical due to the characteristics of the fully coupled added mass matrix of fluid on the entire wetted surface. To overcome this difficulty, an efficient approach based on a reanalysis scheme is proposed in this paper. The proposed method can be applied for cases of higher local modes and beam-like modes for which three-dimensional reduction factors are not known. The three dimensional reduction factors are not needled and thus the restrictions can be removed in the analyses of non-beam like modes or local vibration modes by considering fluid-structure interaction. The validity and calculation efficiency of the proposed method are proved through numerical examples.

• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.36D no.11
• /
• pp.33-44
• /
• 1999

The boundary element method was newly implemented into an optical lithography simulator so that it could evaluate rigorously the topological effects of 2dimensional phase-shifting masks. Both transparent and periodic boundary conditions were applied for the method, and the continuity conditions were used for treating interface nodes. The accuracy of the module developed for simulating aerial images was verified by comparison with analytic solutions and published results. In addition, it was found that our simulator would be more efficient than the conventional method based on the rigorous coupled wave analysis in views of the convergence and the calculation speed. Finally, the optimal design of two phase-shifting masks was performed.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.5
• /
• pp.755-761
• /
• 1987

Two-Dimensional incompressible laminar boundary layer with the reversed flow region is computed using the parially parabolized Navier-Stokes equations in primitive variables. The velocities and the pressure are explicity coupled in the difference equation and the resulting penta-diagonal matrix equations are solved by a streamwise marching technique. The test calculations for the trailing edge region of a finite flat plate and Howarth's linearly retarding flows demonstrate that the method is accurate, efficient and capable of predicting the reversed flow region.

• Journal of Magnetics
• /
• v.13 no.4
• /
• pp.144-148
• /
• 2008

The size and surface effects on the magnetism of a fcc Fe (001) surface was investigated by performing firstprinciples calculations on 3, 5, 7, and 9 monolayers fcc Fe (001) single slabs with two different two-dimensional lattice constants, ${\alpha}=3.44{\AA}$ (System I) and 3.65 ${\AA}$ (System II), using the all-electron full-potential linearized augmented plane wave method within a generalized gradient approximation. The surface layers were coupled ferromagnetically to the subsurface layer in both systems. However, the magnetism of the inner layers was quite different from each other. While all the inner layers of System II were ferromagnetically coupled in the same way as the surface layer, the inner layers of System I showed a peculiar magnetism, bilayer antiferromagnetism. The calculated spin magnetic moments per Fe atom were approximately 2.7 and 2.9 ${\mu}_B$ at the surface for Systems I and II, respectively, due to the almost occupied Fe d-state being in the majority spin state and band narrowing. The spin orientations of System I were out-of-plane regardless of its thickness, whereas the orientation of System II changed from out-of-plane to in-plane with increasing thickness.