• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.65-67
• /
• 2011

The moment-of-fluid (MOF) method is a new volume-tracking method that accurately treats evolving material interfaces. The MOF method uses moment data, namely the material volume fraction, as well as the centroid, for a more accurate representation of the material configuration, interfaces and concomitant volume advection. In this paper, unstructured mesh extension of the MOF method is to be presented. The MOF method is coupled with a stabilized finite element incompressible Navier-Stokes solver for two materials. The effectiveness of the MOF method is demonstrated with a free-surface dam-break problem.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.436-439
• /
• 2008

Computer simulation of multiphase flows has grown dramatically in the last two decades. In this work, we have studied the flow characteristics of immiscible two fluids in a 2-D micro channel driven by pressure gradient using multi-phase lattice Boltzmann method suggested by Shan and Chen(1993) considering the fluid-surface interaction. we tried to examine the effects of parameters related to the two phase flow characteristics and pressure drop in the micro channel like contact angle and channel configuration by changing their value. The results of current study could show the lattice Boltzmann method can simulate the behaviors of two phase flow in the region of micro fluidics well.

• Proceedings of the Korean Vacuum Society Conference
• /
• 1993.02a
• /
• pp.84-86
• /
• 1993

전자빔증착기를 이용하여 적외선영역에서 직접천이형 에너지갭을 갖는 III-V족 2원화합물반도체인 InSb박막을 제작하여 전기-자기적 및 광학적 특성을 조사하였다. 증착된 박막의 X-선회절법으로 분석한 결과, 열처리온도 5$25^{\circ}C$, 열처리시간 30분 열처리한 박막에서 In2O3피크가 없어지고, InSb피크만 나타났으며, 이때의 격자상수 a0=6.49$\AA$이었다. 기판온도가 증가할수록 InSb박막의 결정화가 일어나 전자이동도는 증가하고 비저항은 감소하였다. 온도 100~300K, 자계 500~9000 gauss범위에서 van der ppauw법에 의한 홀효과를 측정한 결과 증착된 박막의 전도형은 n형이었고, 상온에서 캐리어농도는 2.5$\times$10-16cm-3이었으며, 캐리어농도는 2.83$\times$104$\textrm{cm}^2$/V-sec이었다. 적외선 분광기로 측정한 InSb박막의 광학적 에너지갭은 기판온도가 증가할수록 InSb의 에너지갭에 해당하는 값으로 이동하였으며, 기판온도 30$0^{\circ}C$, 열처리온도 5$25^{\circ}C$일 때 측정한 값은 0.173eV였다.

• Journal of the Korean Vacuum Society
• /
• v.3 no.3
• /
• pp.346-354
• /
• 1994

삼원화합물 반도체 CuAlSe2 및 CuGaSe2 단결정의 solid solution 인 CuAl1-xGaxSe2 및 cobalt를 2.0 mol% 첨가한 CuAl1-xGaxSe2 : Co2+ 단결정을 iodine을 수송물질로 사용한 화학수송법으로 성장시켰 다. source material로는 CuAl1-xGaxSe2 의 화학조성비에서 Se를 3.0mol% 과잉으로 첨가하여 합성한 ingot를 사용하였으며 불순물이 첨가된 CuAl1-xGaxSe2:Co2+ 단결정성장시에는 source ma-terial 에 cobalt 분말을 2.0mol% 첨가하였다. X-선 회절무늬로부터 성장된 단결정들이 chalcopyrite 결정구조를 하고 있음을 확인하였으며 격자상수를 구하였다. 광흡수 spectra 측정으로부터 성장된 단결정에서 나타 는 cobalt 불순물에 의한 광흡수 peak가 Td 대칭점에위치한 Co2+ 이온의 에너지 준위들간 전자전이에 의해 나타남을 규명하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.10a
• /
• pp.95-102
• /
• 2002

Studies of the numerical characteristics of implicit upwind schemes, such as upwind ADI, Line Gauss-Seidel(LGS) and Point Gauss-Seidel(LU) algorithms, for preconditioned Navier-Stokes equations ate performed. All the algorithms are expressed in approximate factorization form and Von Neumann stability analysis and convergence studies are made. Preconditioning is applied for efficient convergence at low Mach numbers and low Reynolds numbers. For high aspect ratio computations, the ADI and LGS algorithms show efficient and uniform convergence up to moderate aspect ratio if we adopt viscous preconditioning based on min- CFL/max- VNN time-step definition. The LU algorithm, on the other hand, shows serious deterioration in convergence rate as the grid aspect ratio increases. Computations for practical applications also verify these results.

• The Journal of Engineering Geology
• /
• v.19 no.2
• /
• pp.205-215
• /
• 2009

The conventional numerical models to analyze flow in subsurface porous media under the transient state usually generate numerical oscillation and unstability due to local flux domain for critical cases such as infiltration into initially dry soil during rainfall period. In this case, it is required refined mesh and small time step, but it decrease efficiency of computation. In this study, numerical unstability in discontinuity domain is removed by applying particle tracking algorithm to simulate unsteady subsurface flow with inflow boundary condition. Finally the hybrid LE FEM improving numerical stability is proposed. The hypothetical domains with unsteady uniform and nonuniform flow field were used to demonstrated algorithm verification. In comparison with analytic solution, we obtained reasonable results and conducted simulation of hypothetical 3-D recharge/pumping area. The proposed algorithm can simulate saturated/unsaturated porous media with more practical problems and will greatly contribute to accuracy and stability of numerical computation.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.7 no.3
• /
• pp.69-78
• /
• 2007

In this study, a time-dependent aspect of an embankment failure is considered to simulate a flood inundation map and calculate overtopping discharge induced by an embankment failure. A numerical model has been developed by solving the two dimensional nonlinear shallow water equations with a finite volume method on unstructured grids. To analyze a Riemann problem, the HLLC approximate Riemann solver and the Weighted Averaged Flux method are employed by using a TVD limiter and the source term treatment is also employed by using the operator splitting method. Firstly, the numerical model is applied to a dam break problem and a sloping seawall. Obtained numerical results show good agreements with experimental data. Secondly, the model is applied to a flow induced by an embankment failure by assuming that the width and elevation of embankment are varied with time-dependent functions. As a result of the comparison with each numerical overtopping discharge, established flood inundation discharges in the previous studies are overestimated than the result of the present numerical model.

• Journal of Advanced Marine Engineering and Technology
• /
• v.35 no.6
• /
• pp.809-813
• /
• 2011

Turbulent flows are numerically simulated in the three dimensional inlet duct for heat recovery steam generator. The present study is aimed to analyze the effect of a variation in turbulent flow pattern by the change of roof angle in the transition duct. The finite volume based Navier-Stokes equations with unstructured grids are solved to make clear the flow dynamic phenomena. Reviews are made on with the data of path lines, velocity vectors, dynamic pressure, residuals for numerical convergence and so on. The k-epsilon, k-omega, Reynolds stress and RNG k-epsilon are used for generation of turbulence. Two types of roof angle are applied with and without the swirl in the duct. Turbulent flow patterns could be investigated for the optimum duct design based on the computational results.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.7
• /
• pp.265-274
• /
• 2018

The increasing use of curved beams in buildings, vehicles, ships, and aircraft has prompted studies directed toward the development of an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic curved beams has been the subject of a large number of investigations. Solutions of the relevant differential equations have been obtained traditionally using standard finite difference or finite element methods. These techniques require a great deal of computer time as the number of discrete nodes becomes relatively large under the conditions of complex geometry and loading. One of the efficient procedures for the solution of partial differential equations is the method of differential quadrature. The differential quadrature method (DQM) has been applied to a large number of cases to overcome the difficulties of the complex algorithms of programming for the computer, as well as the excessive use of storage due to the conditions of complex geometry and loading. The in-plane buckling of curved beams considering the extensibility of the arch axis was analyzed under uniformly distributed radial loads using the DQM. The critical loads were calculated for the member with various parameter ratios, boundary conditions, and opening angles. The results were compared with the precise results by other methods for cases, in which they were available. The DQM, using only a limited number of grid points, provided results that agreed very well (less than 0.3%) with the exact ones. New results according to diverse variations were obtained, showing the important roles in the buckling behavior of curved beams, and can be used in comparisons with other numerical solutions or with experimental test data.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.10
• /
• pp.300-309
• /
• 2016

Curved beams are increasingly used in buildings, vehicles, ships, and aircraft, which has resulted in considerable effort towards developing an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic curved beams has been the subject of many investigations. Solutions to the relevant differential equations have traditionally been obtained by the standard finite difference or finite element methods. However, these techniques require a great deal of computer time for a large number of discrete nodes with conditions of complex geometry and loading. One efficient procedure for the solution of partial differential equations is the differential quadrature method (DQM). This method has been applied to many cases to overcome the difficulties of complex algorithms and high storage requirements for complex geometry and loading conditions. Out-of-plane buckling of curved beams with rotatory inertia were analyzed using DQM under uniformly distributed radial loads. Critical loads were calculated for the member with various parameter ratios, boundary conditions, and opening angles. The results were compared with exact results from other methods for available cases. The DQM used only a limited number of grid points and shows very good agreement with the exact results (less than 0.3% error). New results according to diverse variation are also suggested, which show important roles in the buckling behavior of curved beams and can be used for comparisons with other numerical solutions or experimental test data.