• 한국수자원학회논문집
• /
• 제41권12호
• /
• pp.1187-1198
• /
• 2008

국내 많은 저수지들과 그 하류하천은 높은 탁도를 가진 물의 장기적인 방류로 인해 수자원 이용과 수생태계 관리에 많은 어려움을 겪고 있다. 탁도($C_T$)는 물의 탁한 정도를 나타내는 척도로써, 수질과 수환경의 건강을 평가하는 매우 중요한 지표로써 광범위하게 사용되어 왔지만, 모형의 검증에 필요한 실험 자료의 부족으로 인해 탁도 모델링에 대한 연구는 지금까지 매우 부족하였다. 본 연구의 목적은 광범위한 현장 실측자료를 이용하여 성층화된 대청호로 유입한 탁수의 밀도류 거동 모의를 위한 3차원 수리-입자동력학 연동 수치모형인 ELCOM-CAEDYM의 적용성을 검증하는데 있다. 입자크기에 따라 구분된 3개 그룹의 부유물질 (SS) 농도가 모형의 모의변수로 사용되었으며, 모형 변수인 SS와 저수지내 실측값의 $C_T$의 변환을 위해 저수지 지점별로 측정한 SS-$C_T$ 상관관계를 사용하였다. 모의결과는 2004년에 대청호의 회남과 댐앞 지점에서 수심별로 실측한 수온과 탁도 자료와 비교함으로써 검증하였다. 모형은 저수지의 성층구조, 탁수를 포함한 하천 밀도류의 시간에 따른 진행과정을 잘 재현하였으며, 탁도의 수직분포와 크기도 실측값과 부합하였다. 본 연구에서 제시한 3차원 수치모형과 탁도 모델링 방법론은 유사한 탁도 문제를 가지고 있는 다른 저수지에서도 탁수의 최적관리를 위한 지원 도구로써 사용 가능하다고 사료된다.

• Ocean and Polar Research
• /
• 제25권1호
• /
• pp.63-74
• /
• 2003

The build-up of the heat field in shallow coastal water due to a point source has been investigated using an analytical solution of a time-integral form derived by extending the solutions by Holley(1969) and also presented in Harleman (1971). The uniform water depth is assumed with non-isotropic turbulent dispersion. The alongshore-flow is assumed to be uni-directional, spatially uniform and oscillatory. Due to the presence of the oscillatory alongshore-flow, the heat build-up occurs in an oscillatory manner, and the excess temperature thereby fluctuates in that course and even in the quasi-steady state. A series of calculations reveal that proper choices of the decay coefficient as well as dispersion coefficients are critical to the reliable prediction of the excess temperature field. The dispersion coefficients determine the absolute values of the excess temperature and characterize the shoreline profile, particularly within the tidal excursion distance, while the decay coefficient determines the absolute value of the excess temperature and the convergence rate to that of the quasi-steady state. Within the e-folding time scale $1/k_d$ (where $k_d$ is the heat decay coefficient), heat build-up occurs more than 90% of the quasi-steady state values in a region within a tidal excursion distance (L), while occurs increasingly less the farther we go to the downstream direction (about 80% at 1.25L, and 70% at 1.5L). Calculations with onshore and offshore discharges indicate that thermal spreading in the direction of the shoreline is reduced as the shoreline constraint which controls the lateral mixing is reduced. The importance of collecting long-term records of in situ meteorological conditions and clarifying the definition of the heat loss coefficient is addressed. Interactive use of analytical and numerical modeling is recommended as a desirable way to obtain a reliable estimate of the far-field excess temperature along with extensive field measurements.

• 대한토목학회논문집
• /
• 제33권3호
• /
• pp.965-974
• /
• 2013

2011년 7월 27일 우면산 유역에서 일련의 토석류가 발생하여 많은 인명과 재산 피해를 야기했다. 이러한 토석류 피해를 저감하기 위해서는 토석류가 어떻게 이동하고 퇴적되어 하류부에 피해를 주는지를 예측할 수 있는 기술이 필요하다. 이 연구에서는 우면산 토석류에 대해서 FLO-2D 수치모형을 적용하여, 모형의 성능을 평가하고, 이들 토석류의 유동 매개변수를 분석하였다. 우면산 북측사면에서 발생한 토석류의 현장 관측자료와 FLO-2D 모의 결과 분석을 통해서 역산된 우면산 토석류의 대표 항복응력은 1,022 Pa 그리고 대표 동점성은 652 $Pa{\cdot}s$ 인 것으로 나타났다. 선정된 매개변수를 이용하여 우면산 토석류의 퇴적영역을 산정한 결과, 래미안 유역의 경우 63 - 85% 정확도의 피해유역을 재현됨으로서 FLO-2D 모형은 피해유역 산정에 적합한 것으로 나타났다. 하지만, FLO-2D 모형이 침식물질과 우수의 연행작용을 고려할 수 없는 고유의 한계 때문에 연행작용이 현저한 래미안 및 신동아 유역 토석류의 유하속도는 과소평가하는 것으로 나타났다. 유하부에서 침식에 의한 연행작용이 활발하지 않은 형촌마을 유역 토석류는 수로와 저수지 상부 부근에서의 퇴적을 우수한 정확도로 모의되었다. 민감도 분석 결과 유사체적농도와 바닥조도계수 모두 토석류의 흐름두께와 유속 모두에 상당한 영향을 미치는 것으로 나났으므로, 이들 매개변수의 선정에 신중해야 함을 보여주었다.

• Ocean Systems Engineering
• /
• 제9권2호
• /
• pp.111-133
• /
• 2019

Flow through a scaled horizontal axis marine current turbine was numerically simulated after validation and the turbine design was optimized. The computational fluid dynamics (CFD) code Ansys-CFX 16.1 for numerical modeling, an in-house blade element momentum (BEM) code for analytical modeling and an in-house surrogate-based optimization (SBO) code were used to find an optimal turbine design. The blade-pitch angle (${\theta}$) and the number of rotor blades (NR) were taken as design variables. A single objective optimization approach was utilized in the present work. The defined objective function was the turbine's power coefficient ($C_P$). A $3{\times}3$ full-factorial sampling technique was used to define the sample space. This sampling technique gave different turbine designs, which were further evaluated for the objective function by solving the Reynolds-Averaged Navier-Stokes equations (RANS). Finally, the SBO technique with search algorithm produced an optimal design. It is found that the optimal design has improved the objective function by 26.5%. This article presents the solution approach, analysis of the turbine flow field and the predictability of various surrogate based techniques.

• 한국생산제조학회지
• /
• 제23권1호
• /
• pp.87-93
• /
• 2014

This paper provides a methodology for extracting design data from the three-dimensional design software SolidWorks$^{(R)}$, which is based on the existing trochoid pump design equations that are used by hydraulic field engineers. The data extracted from the SolidWorks$^{(R)}$ model are input to a hydraulic analysis software AMESim model to determine the design factors that can influence the properties of a trochoid pump. On the basis of the simulation results, this paper proposes a method to reduce the flow loss by adjusting the outlet angle of the trochoid pump. This proposal was verified by using actual experimental results, which confirmed that adjusting the outlet angle can increase the flow rate. Hence, the results presented in this paper can contribute to the prototyping of a trochoid pump by reducing the cost associated with a trial-and-error design.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2012년도 추계학술대회 논문집
• /
• pp.331-336
• /
• 2012

This paper proposes a new type of high-frequency directional valve controlled by the piezostack actuator associated with displacement amplifier. As a first step, a dynamic model of directional valve which can operate at 200 Hz with a flow rate of 12 l/min is derived by considering pressure drop and flow force. As a second step, an appropriate piezostack is selected by considering actuation force as well as field-dependent displacement. Subsequently, in order to control spool displacement and flow rate a proportional-derivative (PD) controller is designed based on the $3^{rd}$-order valve system. Control performances such as sinusoidal trajectory tracking of the spool displacement in time domain are evaluated. In addition, the field-dependent flow rate is also presented to verify the required performance of the valve system.

• 한국소음진동공학회논문집
• /
• 제22권10호
• /
• pp.1020-1026
• /
• 2012

This paper proposes a new type of high-frequency directional valve controlled by the piezostack actuator associated with displacement amplifier. As a first step, a dynamic model of directional valve which can operate at 200 Hz with a flow rate of 12 litter/min is derived by considering pressure drop and flow force. As a second step, an appropriate piezostack is selected by considering actuation force as well as field-dependent displacement. Subsequently, in order to control spool displacement and flow rate a proportional-derivative(PD) controller is designed based on the 3rd-order valve system. Control performances such as sinusoidal trajectory tracking of the spool displacement in time domain are evaluated. In addition, the field-dependent flow rate is also presented to verify the required performance of the valve system.

• Structural Engineering and Mechanics
• /
• 제81권4호
• /
• pp.523-527
• /
• 2022

The objective of this study is to simulate the flow in pipes with various boundary conditions. Free-pressure fluid model, is used in the pipe based on Navier-Stokes equation. The models are solved by using the numerical method. A problem called "stability of pipes" is used in order to compare frequency and critical fluid velocity. When the initial conditions of problem satisfied the instability conditions, the free-pressure model could accurately predict discontinuities in the solution field. Employing nonlinear strains-displacements, stress-strain energy method the governing equations were derived using Hamilton's principal. Differential quadrature method (DQM) is used for obtaining the frequency and critical fluid velocity. The results of this paper are analyzed by hyperbolic numerical method. Results show that the level of numerical diffusion in the solution field and the range of well-posedness are two important criteria for selecting the two-fluid models. The solutions for predicting the flow variables is approximately equal to the two-pressure model 2. Therefore, the predicted pressure changes profile in the two-pressure model is more consistent with actual physics. Therefore, in numerical modeling of gas-liquid two-phase flows in the vertical pipe, the present model can be applied.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2005년도 제25회 추계학술대회논문집
• /
• pp.450-454
• /
• 2005

축대칭 동축형 램제트 연소기에 대한 반응 유동장 해석을 수행하였다. 2차원 축대칭 Wavier-Stokes 방정식과 낮은 레이놀즈 수 $k-\varepsilon$ 난류 모델을 이용하였고, 유한반응률 화학반응 모델을 적용하였다. 난류 연소 모델인 EDM (Eddy-Dissipation Model)과 층류 반응 모델을 적용한 경우를 서로 비교하였다. 급확대 연소기와 wedge형 보염기를 장착한 동축형 램제트 연소기에 대한 반응 유동장 수치해석을 통해 두 가지 결과를 화염안정 측면에서 서로 비교하였다.

• Korea-Australia Rheology Journal
• /
• 제21권1호
• /
• pp.47-58
• /
• 2009

We have investigated the transient behavior of 1D fully developed Poiseuille viscoelastic flow under finite pressure gradient described by the Oldroyd-B and Leonov constitutive equations. For analysis we employ a simple $2^{nd}$ order discretization scheme such as central difference for space and the Crank-Nicolson for time approximation. For the analysis of the Oldroyd-B model, we also apply the analytical solution, which is obtained again in this work in terms of elementary solution procedure simpler than the previous one (Waters and King, 1970). Both models demonstrate qualitatively similar solutions, but their eventual steady flowrate exhibits noticeable difference due to the absence or presence of shear thinning behavior. In the inertialess flow, the flowrate instantaneously attains a large value corresponding to the Newtonian creeping flow and then decreases to its steady value when the applied pressure gradient is low. However with finite liquid density the flow field shows severe fluctuation even accompanying reversals of flow directions. As the assigned pressure gradient increases, the flowrate achieves its steady value significantly higher than its value during oscillations after quite long period of time. We have also illustrated comparison between 1D and 2D results and possible mechanism of complex 2D flow rearrangement employing a previous solution of [mite element computation. In addition, we discuss some mathematical points regarding missing boundary conditions in 2D modeling due to the change of the type of differential equations when varying from inertialess to inertial flow.