• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.5
• /
• pp.415-423
• /
• 2015

In this study, we present a grid refinement model in the lattice Boltzmann method (LBM) for two-dimensional incompressible fluid flow. That is, the model combines the desirable features of the lattice Boltzmann method and stream function-vorticity formulations. In order to obtain an accurate result, very fine grid (or lattice) is required near the solid boundary. Therefore, the grid refinement model is used in the lattice Boltzmann method for stream function-vorticity formulation. This approach is more efficient in that it can obtain the same accurate solution as that in single-block approach even if few lattices are used for computation. In order to validate the grid refinement approach for the stream function-vorticity formulation, the numerical simulations of lid-driven cavity flows were performed and good results were obtained.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.3
• /
• pp.541-550
• /
• 1988

Natural convection in a square enclosure was investigated numerically for low prandtl number fluids using Finite Element Method. In case of Ra=10$^{4}$, 10$^{5}$ and 10$^{6}$ the temperature gradient decreases gradually at the lower end of the hot wall(or at the upper end of the cold wall) as prandtl number decreases in the range of 0.01 .leq. Pr .leq. 10. Maximum heat transfer occurs at a somewhat higher point from the lower end of hot wall(or at somewhat lower point from the upper end of the cold wall) and it draws near to the lower end of the hot wall(or draws near to the upper end of the cold wall) with increasing prandtl number. The flow in the enclosure appears as an Unicell Pattern for Ra .leq. 10$^{4}$ and secondarily flows(or tertiary flows) appears in the core region for Ra .geq. 10$^{5}$ . The line joining the center of secondary cells skewes in a clockwise direction as the Prandtl number decreases.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.105-105
• /
• 2012

최근 우리나라는 집중호우로 인한 토석류의 발생이 현저하게 증가하고 있으며, 2002년 태풍 루사와 2003년 매미 그리고 2006년 7월 집중호우 등으로 인한 피해의 규모는 여러 조사와 문헌에서도 확인할 수 있다. 2011년 7월, 서울 우면산 일대에서는 집중호우로 인한 토석류가 발생하였으며, 16명의 인명손실을 포함한 큰 피해가 있었다. 우면산 토석류는 십여 개 지역에서 동시다발적으로 발생하였으며, 토석류 발생 유역 특성과 토석류의 유동특성을 분석하기 위한 현장 조사가 진행되었다. 조사한 자료에 따르면 우면산 일대의 일일 강우량은 서초 기상측정소를 기준으로 최대 24시간 누적 324mm 그리고 시간당 최대 68.5mm/hr를 기록하였다. 상업용 소프트웨어인 FLO-2D는 유사농도의 함수로서 점성(viscous)응력, 항복(yield)응력, 난류 및 분산(dispersive) 응력항을 포함하는 2차 유변학(quadratic rheology) 모델을 기본으로 사용하여 이류(mudflow)와 토석류(debris flow)를 모의할 수 있다. FLO-2D는 흐름의 운동량 및 에너지 보존을 고려하여 격자와 시간에 관계없이 유동심도, 속도, 압력을 예측할 수 있으며, 격자 기반의 모델로서 GIS 및 기타 응용 프로그램들과 연동이 쉽다는 장점이 있다. 그러나 하상침식에 의해 유발된 토석류의 체적 증가는 고려 할 수 없으므로 토석류의 전파 및 퇴적영역에서의 토석류 모의에만 사용할 수 있는 단점이 있다. 이 연구의 목적은 FLO-2D 소프트웨어를 이용하여 우면산에서의 토석류 현상을 재현하는 것이다. 우면산 일대에서 발생한 토석류 중 서초구 방배동의 래미안 아파트 부근에서 발생한 토석류에 대하여 수치지도(DEM)와 현장조사를 통해 얻은 지형자료, 해당 지역의 강우량 및 지질 특성 자료 등을 토대로 FLO-2D 모델을 적용하여 토석류의 흐름특성을 검토한다. 토석류 유동 및 퇴적에 대한 가용한 현장관측 자료와의 비교 분석을 통하여 토석류 특성 값을 산정하고, 모델의 적용성을 검증한다.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.7
• /
• pp.625-638
• /
• 2014

In this study, the uncertainties in the critical flow functions (CFFs) calculated by the AGA8-dc equation of state were estimated. To this end, the formulas for enthalpy, entropy, and speed of sound, which are used in calculating the CFF, were expressed in the form of dimensionless Helmholtz free energy and its derivatives, and the uncertainty in Helmholtz free energy was inferred. To consider the variations in the compressibility-dependent variables induced by the variation (i.e., uncertainty) in compressibility, the form of the AGA8-dc equation was modified to have a deviation equal to the uncertainty under each flow condition. For each independent uncertainty component of the CFF, a model for uncertainty contribution was developed. All these changes were applied to GASSOLVER, which is KOGAS's thermodynamic database. As a result, the uncertainties in the CFF were estimated to be 0.025, 0.055, and 0.112 % at 10, 50, and 100 bar, respectively, and are seen to increase with the increase in pressure. Furthermore, these results could explain the deviations in the CFFs across the different labs in which the CFF international comparison test was conducted under the ISO management in 1999.

• Journal of Environmental Science International
• /
• v.6 no.4
• /
• pp.323-333
• /
• 1997

An objective method for the generation of velocity streamfunction is presented for dealing with discretely sampled oceauc data. The method treats a Poisson equation (forced by vorticity) derived from Helmholtz theorem In which streamfunction is obtained by isolating the non-divergent part of the two-dimensional flow field. With a mixed boundary condition and vorticity field estimated from observed field, the method Is Implemented over the Texas-Louisiana show based on the current meter data of the Texas-Louisiana Shelf Circulation and Transport Processes Study (LATEX) measured at 31 moorings for 32 months (April 1992 - November 1994). The resulting streamfunction pattern is quote consistent with observations. The streamfunction field by this method presents an opportunity to initiauze and to verier computer models for local forecasts of enoronmental flow conditions for ell spill, nutrient and plankton transports as well as opportuuty to understand shelf-wlde low-frequency currents.

• Proceedings of the KIEE Conference
• /
• 2004.07c
• /
• pp.1933-1935
• /
• 2004

본 연구에서는 변압기형 플라즈마 전류 모델을 기초로 한 평판형 유도 결합 플라즈마 장치에 대한 회로를 분석하여 임피던스 매칭 특성을 조사하였다. 장치 임피던스는 collisional surface impedance를 기반으로 계산된 플라즈마 임피던스와 안테나 임피던스로 결정된다. 매칭 network에 사용된 회로는 Altcmatc-typc의 회로이고, 매칭 소자인 $C_T$와 $C_L$은 임피던스 매칭 조건을 이용하여 계산하였다. 완전 매칭의 경우에는 $C_T$와 $C_L$을 플라즈마 변수들의 함수로 표현하여 의존성을 분석하고, 불완전 매칭의 경우에는 반사파에 대한 반사계수, 반사율을 계산하였다.

• Journal of the Korean Geotechnical Society
• /
• v.27 no.6
• /
• pp.49-61
• /
• 2011

This paper presents the applicability of rheological models for describing fine-laden debris flows and analyzes the flow characteristics as a function of grain size. Two types of soil samples were used: (1) clayey soils - Mediterranean Sea clays and (2) silty soils - iron ore tailings from Newfoundland, Canada. Clayey soil samples show a typical shear thinning behavior but silty soil samples exhibit the transition from shear thinning to the Bingham fluid as shear rate is increased. It may be due to the fact that the determination of yield stress and plastic viscosity is strongly dependent upon interstructrual interaction and strength evolution between soil particles. So grain size effect produces different flow curves. For modeling debris flows that are mainly composed of fine-grained sediments (<0.075 mm), we need the yield stress and plastic viscosity to mimic the flow patterns like shape of deposition, thickness, length of debris flow, and so on. These values correlate with the liquidity index. Thus one can estimate the debris flow mobility if one can measure the physical properties.

• Journal of the Korean Society of Groundwater Environment
• /
• v.2 no.2
• /
• pp.49-57
• /
• 1995

Areas of Chungyang in Chungcheongnamdo, Cheju island and Georgia in U.S.A. were selected to calculate their groundwater ages and dispersion parameters using tritium. Piston flow model, Completely mixing model and Dispersion binomial model which calculate men residence times, and Dispersion normal model and Dispersion model(C$\_$FF) which calculate mean residence times and dispersion parameters simultaneously were applied. Since the input data, tritium concentrations of rainwaters, lack in part, tritium input function was prepared using the correlation of tritium concentrations of rainwaters of Pohang, Korea and Ottawa, Canada. Similar results of PFM and DBM reflect the intrinsical similarity of two models, assumption of low dispersion. The ages of sites of showing relatively higher tritium concentrations than other sites in the sam, region were not calculated by CMM. Since the calculations of DNM and DM(C$\_$FF/) provide the combination of wide ranges of parameters and groundwater ages, the ranges of dispersion parameters were narrowed down under the assumption that ages calculated by PFM and DBM are correct. Since large variation of tritium concentrations of outflows in a same region may reflect the different characteristics of each groundwater flow regime, using only one specific model on a whole region is not recommended.

• Journal of computational fluids engineering
• /
• v.20 no.3
• /
• pp.94-103
• /
• 2015

The turbulent flow characteristics in the channel flow are investigated using large eddy simulation(LES) of FDS code, built in NIST(USA), in which the near-wall flow is solved by Werner-Wengle wall function. The periodic flow condition is applied in streamwise direction to get the fully developed turbulent flow and symmetric condition is applied in lateral direction. The height of the channel is H=1m, and the length of the channel is 6H, and the lateral length is H. The total grid is $32{\times}32{\times}32$ and $y^+$ is kept above 11 to fulfill the near-wall flow requirement. The Smagorinsky model is used to solve the sub-grid scale stress. Smagorinsky constant $C_s$ is 0.2(default in FDS). Three cases of Reynolds number(10,700, 26,000, 49,000.), based on the channel height, are analyzed. The simulated results are compared with direct numerical simulation(DNS) and particle image velocimetry(PIV) experimental data. The linear low-Re eddy viscosity model of Launder & Sharma and non-linear low-Re eddy viscosity model of Abe-Jang-Leschziner are utilized to compare the results with LES of FDS. Reynolds normal stresses, Reynolds shear stresses, turbulent kinetic energys and mean velocity flows are well compared with DNS and PIV data.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.5
• /
• pp.323-334
• /
• 2020

Since boundary layer transition has a significant impact on the aero-thermodynamic performance of hypersonic flight vehicles, capability of accurate prediction of transition location is essential for design and performance analysis. In this study, γ-Re_θt model is improved to predict transition of hypersonic boundary layers and validated. A coefficient in the production term of the intermittency transport equation that affects the transition onset location is constructed and applied as a function of Mach number, wall temperature, and freestream stagnation temperature based on the similarity numerical solution of compressible boundary layer. To take into account a Mach number dependency of transition onset momentum thickness Reynolds number and transition length, additional correlation equations are determined as function of Mach number and applied to Re_θc and F_length correlations of the baseline model. The suggested model is implemented to a commercial CFD code in consideration of practical use. Analysis of hypersonic flat plate and circular cone boundary layers is carried out by using the model for validation purpose. An improvement of prediction capability with respect to variation of Mach number and unit Reynolds number is identified from the comparison with experimental data.