• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.1
• /
• pp.189-196
• /
• 1990

In the present study, a new method of generating boundary-fitted coordinates systems controlled by control function is introduced. Application of the method to a three-dimensional simply-connected region is the demonstrated. The numerical grid generation has following feat ures, (a) The generated boundary fitted coordinates is well concentrated in near wall region and satisfied orthogonality, (b) The grid control function is fully automatic and well controlled in sharp convex boundary.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.8 no.4 s.23
• /
• pp.136-145
• /
• 2005

특이응력해석을 위한 일반화된 가역상반일 경계적분식이 섬유강화복합재를 모형화한 직교 이방성 크랙평판의 수치해를 위하여 발전시켰다. 이 적분방정식은 평판경계에서의 탄성변위와 트랙션의 변수로 구성된 경계적분식의 형태로 하중이 없다는 두 크랙면의 경계조건과 유한의 탄성변형에너지의 개념에서 경계적분식에 필요한 특성해를 규정하고 대응되는 보조해를 계산하였다. 대칭모우드 I형의 중앙원공크랙평판 및 복합모우드형의 반원편측크랙 일단고정평판의 응력확대계수가 임의의 섬유방향각에 따라서 계산되었다.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.6 no.4
• /
• pp.279-290
• /
• 2004

The differential quadrature method (DQM) for a system of coupled differential equations governing the elastic stability of thin-walled curved members is presented, and is applied to computation of the eigenvalues of out-of-plane buckling of curved beams subjected to uniformly distributed radial loads including a warping contribution. Critical loads with warping, which were found to be significant, are calculated for a single-span wide-flange beam with various end conditions, opening angles, and stiffness parameters. The results are compared with the exact methods available. New results are given for the case of both ends clamped and clamped-simply supported ends without comparison since no data are available The differential quadrature method gives good accuracy and stability compared with previous theoretical results.

• Economic and Environmental Geology
• /
• v.44 no.2
• /
• pp.161-170
• /
• 2011

The multiphase flow simulator, CHEMPS, was developed based on the fractional flow approach reported in the petroleum engineering literature considering fully three phase flow in physically and chemically heterogeneous media. It is a extension of MPS developed by Suk and Yeh (2008) to include the effect of wettability on the migration of NAPL. The fractional flow approach employs water, total liquid saturation and total pressure as the primary variables. Most existing models are limited to two-phase flow and specific boundary conditions when considering physically heterogeneous media. In addition, these models focused mainly on the water-wet media. However, in a real system, variations in wettability between water-wet and oil-wet media often occur. Furthermore, the wetting of porous media by oil can be heterogeneous, or fractional, rather than uniform due to the heterogeneous nature of the subsurface media and the factors that affect the wettability. Therefore, in this study, the chemically heterogeneous media considering fractional wettability as well as physically heterogeneous media were simulated using CHEMPS. In addition, the general boundary conditions were considered to be a combination of two types of boundaries of individual phases, flux-type and Dirichlet type boundaries.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.4
• /
• pp.31-41
• /
• 1997

This paper describes a potential-based panel method formulated for the prediction of the steady performance of a waterjet propulsor. The method employs normal dipoles and sources distributed on the solid surfaces such as the impeller/stator blades, hub and duct, and normal dipoles in the shed wakes trailing the impeller and stator to represent the potential flow around the waterjet propulsor. To define a closed boundary surface, the inlet and outlet open boundary surfaces are introduced where the sources and dipoles are distributed. The kinematic boundary condition on the solid boundary surface is satisfied by requiring that the normal component of the total velocity should vanish. On the inlet surface, the total inflow flux into the duct is specified, and on the outlet surface the conservation of mass principle is applied to evaluate the source strength. The solid surfaces are discretized into a set of quadrilateral panel elements and the strengths of sources and dipoles are assumed constant at each panel. Applying this approximation to the boundary conditions leads to a set of simultaneous equations. Systematic numerical tests show that the present numerical method is fast and stable. In order to validate the present method, sample computations are carried out first for the case of a conventional axial flow fan which has a similar geometry as the waterjet propulsor, and then for the case of a waterjet propulsor on which experiments are carried out at KRISO(Korea Research Institute of Ships and Ocean Engineering).

• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.17 no.4
• /
• pp.418-425
• /
• 1988

Since conventional computer program is workable only with velocity boundary condition, in practical fluid passage such as clean room which usually have wide inlets and outlets, it is not easy to measure velocity itself because of its vector property. Furthermore a certain assumption of velocity at boundaries may lead to physically unreasonable results. From this motivation, we have developed a computer program to predict whole flow field imposed on pressure-based boundary condition which can be measured by relatively simple method. The only additional velocity boundary condition that should be imposed on to make the problem unique, are no slip condition at all walls and zero cross stream velocity at inlet. The result of present study was compared with that by Bernoulli equation being used practically. They were coincident well each other within 5%, therefore the validity of the present method is proved. In the present work, the flow field in a clean room subject to pressure-based boundary conditions at an inlet and two exits was predicted numerically. The pressure difference between the inlet and the left exit which keeps relatively low pressure among two exits is fixed as 150[Pa] and the pressure at the right exit is varied from zero to 150[Pa] by the increment of 25[Pa]. For each cases the flow characteristics in the clean room, the velocity profile at the inlet, and the flow rate through the two exits was predicted. The flow rate through the right exit imposed on relatively higher pressure than the left exit decreased linearly according to the increase of pressure of the right exit. When the pressure of the right exit is increased enough to cause back flow at the exit, the flow rate is rapidly decreased.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.17 no.6 s.123
• /
• pp.531-538
• /
• 2007

A new method for free nitration analysis using series functions is proposed to obtain the eigenvalues of arbitrarily shaped, polygonal plates with clamped edges. Since a general solution used in the method satisfies the equation of motion for the transverse vibration of a plate, the method offers very accurate eigenvalues, compared to FEM or BEM results. In addition, the method can minimize the amount of numerical calculation because it has the advantage of not needing to divide the plate of interest. Two case studies show that the proposed method is valid and accurate when the eigenvalues by the proposed method are compared to those by FEM (NASTRAN) or another analytical method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.7
• /
• pp.1255-1265
• /
• 1992

Under uniform heat flow, the thermal stress intensity factors for the interfacial crack on a rigid cusp-type inclusion are determined by Hilbert problem expressed with complex variable. The thermal stress intensity factors are expressed in terms of the periodic function of heat flow angle. When the tip of the interfacial crack meets that of the cusp crack, the thermal stress intensity factors have singularities. The thermal stress intensity factors at the interfacial crack tip located in the distance from the cusp crack tip vary with the location of the interfacial crack tip. From the results of the analysis, the complex potential functions and the thermal stress intensity factors for the cusp-type inclusion without the interfacial crack are derived under the cusp surface boundary conditions insulated or fixed to zero relative temperature.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.9 no.4
• /
• pp.159-168
• /
• 2005

This study carries out a finite difference nonlinear analysis of anisotropic advanced composite plate structures with various layer sequences. In the numerical analysis of various mechanical problems involving complex partial differential equations, the finite difference method (FDM) developed in this study has an advantage over the finite element method in its ability to avoid mesh generation and numerical integration. Many studies in FDM have been made on clamped or simple boundary conditions using merely an energy approach. These approaches cannot be satisfied, however, with pivotal points along the free boundary. Therefore, this study addresses the nonlinear problem of anisotropic plates by adopting a finite difference modeling elimination of pivotal difference points in the case of a free boundary condition. Complex nonlinear behaviors of composite plate structures for various parameters, especially for layer sequences, are analyzed using the proposed approach.

• Tunnel and Underground Space
• /
• v.12 no.3
• /
• pp.142-151
• /
• 2002

Claesson(2001)'s analytical solution, and two numerical models with Dirichlet and Neuman interior boundary condition respectively were investigated to estimate the transient temperature distribution with distances from the Taejon underground food cold storage pilot cavern. Claesson's solution, which is based on constant temperature boundary condition at the rock wall during a temperature decline step, showed relatively good agreement with temperature measurements in the rock mass in order of average error difference, 0.89$\^{C}$ without any adjustments on laboratory thermal properties to represent the rock mass. For the numerical model with heat flux through the rock wall, a boundary condition setting technique was newly proposed to overcome the difficulty of prescribing variable convective heat tranfer coefficient and far-field air temperature inside the cavern as they may be certainly changed according to the cooling-down time. The results showed also good agreement with measurements in order of average error difference, 1.58$\^{C}$, and were compared to those of the numerical model with fixed temperature at the rock wall. Finally, the most proper procedure to precisely predict the temperature profile around a cavern was proposed as a series of analysis steps including an analytical exact solution and numerical models.