An elliptic blending Reynolds stress transport equation model for Newtonian fluids has been extended to predict polymer-induced drag reduction FENE-P fluids. The conformation tensor equation which is related to the polymer stress is adopted from the model form of Resende et al., and the models of redistribution and dissipation rate terms for the Reynolds stress transport equation are considered by the elliptic blending equation. Also, the new model terms for viscoelastic turbulent transport and viscoelastic dissipation in the Reynolds stress transport equation are introduced to consider the polymer additives effect. The prediction results are directly compared to the DNS data to assess the performance of the present model predictions.

The characteristics of the turbulent MHD channel flows are analysed within the elliptic blending    model. The evaluation of additional terms representing the MDH interactions in the transport equation for the turbulent kinetic energy and dissipation rate is carried out firstly. And then, the improvement of model coefficients for the additional terms is achieved. Regardless of the magnetic field orientation, perpendicular and parrel to the main flow, the model coefficients for the additional terms are not changed to maintain the generality of the present model. The prediction results are directly compared to the DNS data to assess the performance of the present model.

An algebraic model for turbulent heat fluxes which is originally suggested by Suga & Abe is modified on the basis of the elliptic blending equation. In order to satisfy the heat transfer characteristics of near-wall region and the flow center region far away from the wall, the model coefficients of the algebraic heat flux model are modified by using the solution of the elliptic blending equation. The predictions of turbulent heat transfer in a plane channel flow have been carried out with constant wall heat flux and constant wall temperature difference boundary conditions respectively. Also, the predictions are performed at various Prandtl numbers to test the applicability of the model. The prediction results show that the distributions of the turbulent heat fluxes and mean temperature are well captured by the modified algebraic heat flux model

본 연구는 격자수, 첫 번째 격자까지의 거리(YP+), 난류모델 그리고 이산화 방법에 따른 해의 변화량을 조사하였다. 대상선박은 KVLCC이며, 격자구성과 유동해석은 상용코드인 Gridgen V15와 FLUENT를 사용하였다. 검토는 2가지 파트로 나누어서 수행하였다. 첫 번째 파트는 격자수, 난류모델 그리고 이산화 방법의 조합에 따른 해의 영향성을 평가하였다. 두 번째 파트는 적합한 YP+ 선정에 초점을 두었다. 격자수와 이산화 방법이 동일한 경우 마찰저항은 난류모델에 따라 약 1 % 내에서 차이를 보였으나, 압력저항은 약 9 %의 큰 차이를 보였다. YP+와 이산화 방법이 동일한 경우 YP+를 30과 50으로 설정하였을 때 마찰저항은 난류모델에 따라 약 1 % 내에서 차이를 보였으나, 100에서는 약 3 % 차이를 보였다. 반면, 압력저항은 YP+값에 무관하게 난류모델에 따라 약 10 % 차이를 보였다. 난류모델과 이산화 방법이 동일한 경우 격자 수 변화 따라 마찰저항, 압력저항 그리고 전 저항 모두 큰 차이를 보이지 않았다. 난류모델과 이산화 방법이 동일한 경우 YP+의 변화에 따라 마찰저항은 5~8 %의 큰 차이를 보였고, 압력저항은 큰 차이를 보이지 않았다.

This study describes the amendment of Durbin's k-∈-v2-f model and its application to turbulent channel flow to test the model’s performance. Modeling redistribution and dissipation rate terms for the scalar v2 transport equation is considered by the elliptic blending equation which is used in the second moment closure generally. The prediction results are directly compared to the DNS and Durbin's original k-∈-v2-f model to assess the performance of the new model predictions and to show their reasonable agreement with the DNS and Durbin's model for all the flow characteristics that are analyzed for the present study.

관망의 해석은 관망의 설치, 운영에서 매우 중요한 요소이다. 이를 위한 연구 방법은 1차원 모형과 3차원 모형이 있으나 각각은 정확도와 계산속도의 면에서 관망 해석에 장점과 한계점을 보인다. 2차원 모형의 경우 1차원 모형과 달리 난류 점성을 이용하여 보다 정확한 마찰을 모의할 수 있게 하는데, 이때 사용되는 난류 점성을 모의하는 모형은 Five-Region Turbulence Model을 이용하였다. 그러나 이들 모형에 사용되는 매개 변수들은 실험에 의한 값이기 때문에 이러한 매개변수에 따른 2차원 부정류 해석 방법의 응답 특성을 알아보았다. 이를 위해 실제 상수관망에서의 수격압 데이터와 부정류 해석 모형으로 부터의 수격압 데이터를 비교하고, 다양한 매개변수에 따른 속도 분포 모의 특성을 조사하였다. 이를 통해서, 난류 해석 모형의 마찰 거동과 매개변수와의 관계를 규명하였다.

본 연구에서는 중층 밀도류를 모의할 수 있는 k-ε 난류모형의 지배방정식을 제시하고 수치모의를 수행하였다. 깊은 수체에 모형을 적용하여 중층 밀도류를 모의하고 게산된 유속과 초과밀도 분포를 분석하였다. 밀도류의 주 흐름방향을 따라 물 연행으로 인해 유속이 감소되는 것과 Richardson 수의 증가로 인해 유속 변화율이 감소되는 것을 관찰하였다. 유속과 초과밀도의 유사성을 확인하였으나, 난류운동에너지와 소산율의 유사성에서는 보이지 않았다. k-ε 모형의 모의 결과를 이용하여 중층 밀도류의 층적분 모형에서 사용될 수 있는 형상계수를 계산하였다. 또한, 층적분 모형을 이용하여 k-ε 모형에서 사용되는 부력관련 모형상수 (c₃ɛ)와 부피팽창계수 (β0)를 계산하였다.

이군집 응집현상은 수자원환경에서 점착성 유사가 결합-해체의 과정을 통해 응집핵-응집체의 이군집 입자크기분포(Biomodal Floc Size Distribution)를 형성하는 일련의 과정을 의미한다. 본 연구는 저난류 및 고난류 두 가지 조건에서 수행한 응집-침전관 실험결과를 바탕으로 이군집 응집모형(TCPBE: Two Class Population Balance Equation)의 적용성을 단일군집 응집모형(SCPBE: Single Class Population Balance Equation) 및 다군집 응집모형(MCPBE: Multi Class Population Balance Equation)과 비교·평가하였다. 기존 SCPBE에 비하여, TCPBE는 응집핵-응집체의 상호작용 및 침강속도차에 따른 응집 기작을 모의할 수 있었다. 또한, 3개의 연립미분방정식을 가진 TCPBE는 30개 미분방정식을 가진 다군집 응집모형(MCPBE: Multi Class Population Balance Equation)과 대등한 모의 결과를 나타내었다. 따라서 TCPBE는 이군집 응집현상을 모의 할 수 있는 가장 단순한 모델로 검증되었고, 향후 수자원환경이나 수처리 공정에 다양하게 적용할 수 있으리라 판단된다.

본 연구는 난류현상의 모형화를 위해 널리 이용되는 k-과 k- 난류모형을 비교하는 것이 목적으로, 횡방향 흐름이 무시될 수 있는 U-튜브 모양의 터널형 수로 내 높은 레이놀즈수를 가진 진동 경계층 흐름에 두 난류해석방법을 적용하였다. 난류모형의 적용은 1차원 연직 모형을 통해 이루어지며, 수치 모의 결과, 유속의 분포와 난류운동에너지 (turbulent kinetic energy) 모두에서 두 모형이 매우 유사한 결과를 나타낸다. 이를 통하여, 횡방향

복잡한 연료실 내의 유동현상을 난류모형을 적용하여 해석하였다. 적용한 모형은 k-ε, k-Ω, spalart-allmaras, reynolds stress이고, 연로실내의 격자는 혼합격자(hybrid grid) 이다. 속도벡터, 압력분포, 반복계산(iteration)에 의한 잔류치(residual), 동력학적 양정(dynamic head) 등을 모사하였다. 4개의 난류모형을 연료실 유동에 적용하였다. 3차원 수치실험에 앞서, 수치검증을 위하여 2차원 물체 주위의 점성유동을 k-ε 난류모형을 적용하여 모사하였고 항력계수를 비교하였다.

장애물이 있는 배관속의 점성유동을 다양한 난류모형을 적용하여 해석하였다. 적용한 난류모형은 k-ε, k-Ω, Spalart-Allmaras, Reynolds stress 이고, 배관내의 격자는 구조격자(structured grid) 이다. 속도벡터, 압력분포 반복계산(iteration)에 의한 잔류치(residual), 양정(dynamic head) 등을 모사하였다. 4개의 난류모형을 배관유동에 적용하였고 상용 프로그램을 사용하여 해석을 수행하였다.

The flow of non-rotation atmosphere with uniform stratification and wind past an isolated three dimensional topography obstacle is investigated with three-dimensional hydrostatic and non- hydrostatic numerical model. The characteristic of turbulence created the back of topography obstacle is usually defined by Froude number which is the function of upstream wind speed, the height of topography obstacle, and atmospheric stability. Turbulence tends to be formed more easily at the non-hydrostatic model than hydrostatic model. Especially, the difference between flow patterns of two models generated by isolated obstacle is more clear under low Froude number. The difference of flow patterns can be only seen at relatively low altitude, but at high altitude the patterns of two models are almost same. In this research, wind velocity in the parameters related with Froude number have great sensitivity at responsibility of numerical models. and slop of obstacle is also important factor at the flow pattern regardless of the species of numerical model

가로흐름이 존재하는 수역으로 방출되는 온배수 해석에 3가지 난류모형의 사용성 평가를 위해서 근역 2차원 수치모형을 개발하였다. 적용한 난류모형은 2-방정식 난류모형인 {{{{ { k}_{ } }}}}- 및 {{{{ { k}_{ } }}}}-ι 난류모형과 {{{{ { k}_{ } }}}}- 난류모형에 부력생성 항 및 난류 열 플럭스 항 결정을 위한 변동온도 평균자승항 및 이의 감쇠율에 대한 전달 방정식을 추가한 4-방정식 난류모형이다. 개발된 모형은 간단