Park, Il Heum;Cho, Young Jun;Kim, Tae Yun;Lee, Moon Ock;Hwang, Sung Su
해양환경안전학회지
/
제20권3호
/
pp.323-333
/
2014
In this paper, the wakes behind a square cylinder were simulated using two kinds of different turbulence models for the eddy viscosity concept such as the zero- and the one-equation model in which the former is the mixing length model and the latter is the k-equation model. For comparison between numerical and analytical solutions, we employed three skill assessments: the correlation coefficient(r) for the similarity of the wake shape, the error of maximum velocity difference(EMVD) for the accuracy of wake velocity and the ratio of drag coefficient(RDC) for the pressure distribution around the structure. On the basis of the numerical results, the feasibility of each model for wake simulation was discussed and a suitable value for the empirical constant was suggested in these turbulence models. The zero-equation model, known as the simplest turbulence model, overestimated the EMVD and its absolute mean error(AME) for r, EMVD and RDC was ranging from 20.3 % to 56.3 % for all test. But the AME by the one-equation model was ranging from 3.4 % to 19.9 %. The predicted values of the one-equation model substantially agreed with the analytical solutions at the empirical mixing length scale $L=0.6b_{1/2}$ with the AME of 3.4 %. Therefore it was concluded that the one-equation model was suitable for the wake simulation behind a square cylinder when the empirical constant for eddy viscosity would be properly chosen.
The wakes behind a square cylinder were simulated using two-equation turbulence models, $k-{\varepsilon}$ and RNG $k-{\varepsilon}$ models. For comparisons between the model predictions and analytical solutions, we employed three skill assessments:, the correlation coefficient for the similarity of the wake shape, the error of maximum velocity difference (EMVD) of the accuracy of wake velocity, and the ratio of drag coefficient (RDC) for the flow patterns as in the authors' previous study. On the basis of the calculated results, we discussed the feasibility of each model for wake simulation and suggested a suitable value for an eddy viscosity related constant in each turbulence model. The $k-{\varepsilon}$ model underestimated the drag coefficient by over 40 %, and its performance was worse than that in the previous study with one-equation and mixing length models, resulting from the empirical constants in the ${\varepsilon}-equation$. In the RNG $k-{\varepsilon}$ model experiments, when an eddy viscosity related constant was six times higher than the suggested value, the model results were yielded good predictions compared with the analytical solutions. Then, the values of EMVD and RDC were 3.8 % and 3.2 %, respectively. The results of the turbulence model simulations indicated that the RNG $k-{\varepsilon}$ model results successfully represented wakes behind the square cylinder, and the mean error for all skill assessments was less than 4 %.
A modified k-$\epsilon$ model is proposed for calculation of transitional boundary layer flows. In order to develop the eddy viscosity model for the problem, the flow is divided into three regions; namely, pre-transition region, transition region and fully turbulent region. The pre-transition eddy-viscosity is formulated by extending the mixing Length concept. In the transition region, the eddy-viscosity model employs two length scales, i.e., pre-transition length scale and turbulent length scale pertaining to the regions upstream and the downstream, respectively, and a university model of stream-wise intermittency variation is used as a function bridging the pre-transition region and the fully turbulent region. The proposed model is applied to calculate three benchmark cases of the transitional boundary layer flows with different free-stream turbulent intensity ( $1\%{\~}6\%$ ) under zero-pressure gradient. It was found that the profiles of mom velocity and turbulent intensity, local maximum of velocity fluctuations, their locations as well as the stream-wise variation of integral properties such as skin friction, shape factor and maximum velocity fluctuations are very satisfactorily Predicted throughout the flow regions.
A new eddy viscosity equation was formulated from assumption of turbulence length scale equation and specific dissipation ratio equation. Then, a set of turbulence model equations for the turbulent kinetic energy ${\kappa}$, the viscosity ${\nu}_t$, and the intermittency factor ${\gamma}$ is proposed by considering the entrainment effect. Closure coefficients are determined by experimental data and resorting to numerical optimization. Present model has been applied to compute four representative cases of free shear flows and successfully compared with experimental data. In particular, the spreading rate, the centreline mean velocity and the profiles of intermittency are calculated with improved accuracy. Also, the proposed ${\nu}_t-{\kappa}-{\gamma}$ model was applied to channel flow by considering the wall effect and the results show good agreements with the Direct Numerical Simulation data.
Recent observations reveal that the Tsushima Current has a double-cored structure downstream of the Tsushima Island. To explain this, a simple analytical model is proposed based on the assumption of small lateral eddy viscosity. This model suggests that an otherwise uniform current becomes to have a stream core immediately after it enters a channel due to the action of lateral friction. The core is initially broad but becomes sharper downstream. The speed at which the core develops depends on the intensity of lateral eddy viscosity. Likewise, a single-cored stream changes rapidly to a double-cored stream when it passes through an island located in the center of the channel. When the stream leaves the island behind, the reverse process from the double-to single-cored structures takes place. In this case, however, the double-cored structure is retained for a significant distance from the island. Overall, this model suggests that the double-cored structure of the Tsushima Current observed downstream of the Tsushima Island Is created by the lateral friction exerted by the Tsushima Island.
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.
난류모형을 도입하여 풍성류에 대한 3차원 수치모형을 수립하고, 풍성류 계산에 있어서 연직 와점성계수 산정방법들에 대한 비교평가를 수행하였다. 검토된 방법은 함수형(0-방정식), 1-방정식, 2개의 2-방정식 난류모형($ extsc{k}$-$\varepsilon$과 $textsc{k}$-ι)이며, 난류모형은 난류확산 특성길이의 연직분포에 대한 비교를 통하여 각 방법의 특징이 검토되었다. 양단이 막힌 수로와 순환수로에서 유속의 연직분포 계산결과를 수리실험자료와 비교검증하였다. 분석결과에 의하면 1-방정식 난류모형에 적합한 난류확산 특성길이 산정식은 포물선형이었으며, 2-방정식 난류모형($textsc{k}$-ι 모형과 $textsc{k}$-$\varepsilon$모형)은 대체로 실험치와 일치하는 경향을 보였다. 바람에 의한 영향이 전수심에 미치지 않는 경우에 가정된 연직 와점성계수의 분포는 바람이 영향을 미치는 수심까지만 적용되며, 영향이 없는 수심에서는 연직확산이 거의 일어나지 않아 적정성계수의 크기가 0에 가까웠다.
Present study investigated the influence of time step size, turbulent eddy viscosity, and the number of layer on rapid and unsteady propagation of dam break flow. When the time step size had a value such that it resulted in Cr of 0.89, a significant numerical oscillation was observed in the vicinity of the wave front. Higher turbulent viscosity ensured smooth and mild slope of velocity and water stage compared with the flow behavior by no viscosity. The vertical velocity at the lower layer positioned near the bottom showed lower velocity compared with other layers.
The turbulent flows in a tunnel mock-up($10L{\times}0.5W{\times}0.25H$ m3 : scale reduction 1/20) with rectangular cross section were investigated. The instantaneous velocity fields of Re = 49,029, 89,571 were measured by the 2-D PIV system which is consisted of double pulsed Nd:Yag laser and the tracer particles in the straight-duct mock-up where the flows were fully developed. The mean velocity profiles were taken from the ensemble averages of 1,000 instantaneous velocity fields. Simultaneously, numerical simulations(RANS) were performed to compare with experimental data using STREAM code. Non-linear eddy viscosity model (NLEVM : Abe-Jang-Leschziner Eddy Viscosity Model) was employed to resolve the turbulent flows in the duct. The calculated mean velocity profiles were well compared with PIV results. In the log-law profiles, the experimental data were in good agreement with numerical simulations all the way to the wake region except the viscous sub-layer (near wall region).
VIC (Vortex-In-Cell) method for viscous incompressible flow is presented to simulate the wake behind a modified NACA16 foil. With uniform rectangular grid, the velocity in field is calculated using streamfunction from vorticity field by solving the Poisson equation in which FFT(Fast Fourier Transform) is combined with 2nd order finite difference scheme. Here, LES(Large Eddy Simulation) with Smagorinsky model is applied for turbulence calculation. Effective viscosity is formulated using magnitude of strain tensor(or vorticity). Then the turbulent diffusion as well as viscous diffusion becomes particle strength exchange(PSE) with averaged eddy viscosity. The well-established panel method is combined to obtain the irrotational velocity and to apply the no-penetration boundary condition on the body panel. And wall diffusion is used for no-slip condition numerical results of turbulent stresses are compared with experimental results (Bourgoyne, 2003). Before comparing process, LES(Large Eddy Simulation) SGS(Subgrid scale) stress is transformed Reynolds averaged stress (Winckelmans, 2001).
본 웹사이트에 게시된 이메일 주소가 전자우편 수집 프로그램이나
그 밖의 기술적 장치를 이용하여 무단으로 수집되는 것을 거부하며,
이를 위반시 정보통신망법에 의해 형사 처벌됨을 유념하시기 바랍니다.
[게시일 2004년 10월 1일]
이용약관
제 1 장 총칙
제 1 조 (목적)
이 이용약관은 KoreaScience 홈페이지(이하 “당 사이트”)에서 제공하는 인터넷 서비스(이하 '서비스')의 가입조건 및 이용에 관한 제반 사항과 기타 필요한 사항을 구체적으로 규정함을 목적으로 합니다.
제 2 조 (용어의 정의)
① "이용자"라 함은 당 사이트에 접속하여 이 약관에 따라 당 사이트가 제공하는 서비스를 받는 회원 및 비회원을
말합니다.
② "회원"이라 함은 서비스를 이용하기 위하여 당 사이트에 개인정보를 제공하여 아이디(ID)와 비밀번호를 부여
받은 자를 말합니다.
③ "회원 아이디(ID)"라 함은 회원의 식별 및 서비스 이용을 위하여 자신이 선정한 문자 및 숫자의 조합을
말합니다.
④ "비밀번호(패스워드)"라 함은 회원이 자신의 비밀보호를 위하여 선정한 문자 및 숫자의 조합을 말합니다.
제 3 조 (이용약관의 효력 및 변경)
① 이 약관은 당 사이트에 게시하거나 기타의 방법으로 회원에게 공지함으로써 효력이 발생합니다.
② 당 사이트는 이 약관을 개정할 경우에 적용일자 및 개정사유를 명시하여 현행 약관과 함께 당 사이트의
초기화면에 그 적용일자 7일 이전부터 적용일자 전일까지 공지합니다. 다만, 회원에게 불리하게 약관내용을
변경하는 경우에는 최소한 30일 이상의 사전 유예기간을 두고 공지합니다. 이 경우 당 사이트는 개정 전
내용과 개정 후 내용을 명확하게 비교하여 이용자가 알기 쉽도록 표시합니다.
제 4 조(약관 외 준칙)
① 이 약관은 당 사이트가 제공하는 서비스에 관한 이용안내와 함께 적용됩니다.
② 이 약관에 명시되지 아니한 사항은 관계법령의 규정이 적용됩니다.
제 2 장 이용계약의 체결
제 5 조 (이용계약의 성립 등)
① 이용계약은 이용고객이 당 사이트가 정한 약관에 「동의합니다」를 선택하고, 당 사이트가 정한
온라인신청양식을 작성하여 서비스 이용을 신청한 후, 당 사이트가 이를 승낙함으로써 성립합니다.
② 제1항의 승낙은 당 사이트가 제공하는 과학기술정보검색, 맞춤정보, 서지정보 등 다른 서비스의 이용승낙을
포함합니다.
제 6 조 (회원가입)
서비스를 이용하고자 하는 고객은 당 사이트에서 정한 회원가입양식에 개인정보를 기재하여 가입을 하여야 합니다.
제 7 조 (개인정보의 보호 및 사용)
당 사이트는 관계법령이 정하는 바에 따라 회원 등록정보를 포함한 회원의 개인정보를 보호하기 위해 노력합니다. 회원 개인정보의 보호 및 사용에 대해서는 관련법령 및 당 사이트의 개인정보 보호정책이 적용됩니다.
제 8 조 (이용 신청의 승낙과 제한)
① 당 사이트는 제6조의 규정에 의한 이용신청고객에 대하여 서비스 이용을 승낙합니다.
② 당 사이트는 아래사항에 해당하는 경우에 대해서 승낙하지 아니 합니다.
- 이용계약 신청서의 내용을 허위로 기재한 경우
- 기타 규정한 제반사항을 위반하며 신청하는 경우
제 9 조 (회원 ID 부여 및 변경 등)
① 당 사이트는 이용고객에 대하여 약관에 정하는 바에 따라 자신이 선정한 회원 ID를 부여합니다.
② 회원 ID는 원칙적으로 변경이 불가하며 부득이한 사유로 인하여 변경 하고자 하는 경우에는 해당 ID를
해지하고 재가입해야 합니다.
③ 기타 회원 개인정보 관리 및 변경 등에 관한 사항은 서비스별 안내에 정하는 바에 의합니다.
제 3 장 계약 당사자의 의무
제 10 조 (KISTI의 의무)
① 당 사이트는 이용고객이 희망한 서비스 제공 개시일에 특별한 사정이 없는 한 서비스를 이용할 수 있도록
하여야 합니다.
② 당 사이트는 개인정보 보호를 위해 보안시스템을 구축하며 개인정보 보호정책을 공시하고 준수합니다.
③ 당 사이트는 회원으로부터 제기되는 의견이나 불만이 정당하다고 객관적으로 인정될 경우에는 적절한 절차를
거쳐 즉시 처리하여야 합니다. 다만, 즉시 처리가 곤란한 경우는 회원에게 그 사유와 처리일정을 통보하여야
합니다.
제 11 조 (회원의 의무)
① 이용자는 회원가입 신청 또는 회원정보 변경 시 실명으로 모든 사항을 사실에 근거하여 작성하여야 하며,
허위 또는 타인의 정보를 등록할 경우 일체의 권리를 주장할 수 없습니다.
② 당 사이트가 관계법령 및 개인정보 보호정책에 의거하여 그 책임을 지는 경우를 제외하고 회원에게 부여된
ID의 비밀번호 관리소홀, 부정사용에 의하여 발생하는 모든 결과에 대한 책임은 회원에게 있습니다.
③ 회원은 당 사이트 및 제 3자의 지적 재산권을 침해해서는 안 됩니다.
제 4 장 서비스의 이용
제 12 조 (서비스 이용 시간)
① 서비스 이용은 당 사이트의 업무상 또는 기술상 특별한 지장이 없는 한 연중무휴, 1일 24시간 운영을
원칙으로 합니다. 단, 당 사이트는 시스템 정기점검, 증설 및 교체를 위해 당 사이트가 정한 날이나 시간에
서비스를 일시 중단할 수 있으며, 예정되어 있는 작업으로 인한 서비스 일시중단은 당 사이트 홈페이지를
통해 사전에 공지합니다.
② 당 사이트는 서비스를 특정범위로 분할하여 각 범위별로 이용가능시간을 별도로 지정할 수 있습니다. 다만
이 경우 그 내용을 공지합니다.
제 13 조 (홈페이지 저작권)
① NDSL에서 제공하는 모든 저작물의 저작권은 원저작자에게 있으며, KISTI는 복제/배포/전송권을 확보하고
있습니다.
② NDSL에서 제공하는 콘텐츠를 상업적 및 기타 영리목적으로 복제/배포/전송할 경우 사전에 KISTI의 허락을
받아야 합니다.
③ NDSL에서 제공하는 콘텐츠를 보도, 비평, 교육, 연구 등을 위하여 정당한 범위 안에서 공정한 관행에
합치되게 인용할 수 있습니다.
④ NDSL에서 제공하는 콘텐츠를 무단 복제, 전송, 배포 기타 저작권법에 위반되는 방법으로 이용할 경우
저작권법 제136조에 따라 5년 이하의 징역 또는 5천만 원 이하의 벌금에 처해질 수 있습니다.
제 14 조 (유료서비스)
① 당 사이트 및 협력기관이 정한 유료서비스(원문복사 등)는 별도로 정해진 바에 따르며, 변경사항은 시행 전에
당 사이트 홈페이지를 통하여 회원에게 공지합니다.
② 유료서비스를 이용하려는 회원은 정해진 요금체계에 따라 요금을 납부해야 합니다.
제 5 장 계약 해지 및 이용 제한
제 15 조 (계약 해지)
회원이 이용계약을 해지하고자 하는 때에는 [가입해지] 메뉴를 이용해 직접 해지해야 합니다.
제 16 조 (서비스 이용제한)
① 당 사이트는 회원이 서비스 이용내용에 있어서 본 약관 제 11조 내용을 위반하거나, 다음 각 호에 해당하는
경우 서비스 이용을 제한할 수 있습니다.
- 2년 이상 서비스를 이용한 적이 없는 경우
- 기타 정상적인 서비스 운영에 방해가 될 경우
② 상기 이용제한 규정에 따라 서비스를 이용하는 회원에게 서비스 이용에 대하여 별도 공지 없이 서비스 이용의
일시정지, 이용계약 해지 할 수 있습니다.
제 17 조 (전자우편주소 수집 금지)
회원은 전자우편주소 추출기 등을 이용하여 전자우편주소를 수집 또는 제3자에게 제공할 수 없습니다.
제 6 장 손해배상 및 기타사항
제 18 조 (손해배상)
당 사이트는 무료로 제공되는 서비스와 관련하여 회원에게 어떠한 손해가 발생하더라도 당 사이트가 고의 또는 과실로 인한 손해발생을 제외하고는 이에 대하여 책임을 부담하지 아니합니다.
제 19 조 (관할 법원)
서비스 이용으로 발생한 분쟁에 대해 소송이 제기되는 경우 민사 소송법상의 관할 법원에 제기합니다.
[부 칙]
1. (시행일) 이 약관은 2016년 9월 5일부터 적용되며, 종전 약관은 본 약관으로 대체되며, 개정된 약관의 적용일 이전 가입자도 개정된 약관의 적용을 받습니다.