• 한국군사과학기술학회지
• /
• 제2권2호
• /
• pp.279-286
• /
• 1999

Wall-bounded 유동과 달리 자유 전단 유동은 Hyperbolic Tangent Profile을 가지고 비점성 불안정에 의해 지배된다. 따라서 자유 전단 유동에서 난류로의 천이과정은 비점성 불안정성 이론에 의해 해석되어 진다. 본 연구는 분리판(Splitter plate)에 의해 분리된 속도가 다른 자유 유속의 혼합에 의해 형성되는 혼합층에서 와류병합과정에 대한 속도비의 영향에 대하여 연구한다. 속도비는, R, $ \frac{U_1-U_2}{U_1+U_2}$ 로 정의되며, 여기서 $U_1$은 분리판 위에서의 자유 유속을 그리고 $U_2$ 는 분리판 아래에서의 자유 유속을 나타낸다. 본 연구에서 와류구조의 병합작용을 분석하기 위하여 2차원 비정상 Large-Eddy Simulation 방법을 적용하였다. 속도비의 변화에 따라 혼합층에서 불안정 Wave가 성장하게 되고, 유체는 2차원 와류구조에서 Roll-up한다. 이러한 2차원 와류구조는 주위의 다른 와류구조와 상호작용을 하게 되고 하나의 커다란 와류구조를 형성하는 것을 볼 수 있다. 혼합층에서 와류병합과정은 반복적으로 일어나는 것을 알 수 있었고, 이 결과를 이용하여 혼합층의 성장을 제어할 수 있다.

• 유변학
• /
• 제3권1호
• /
• pp.56-67
• /
• 1991

압출공정에서 발생하는 전달현상을 이해하기 위하여 싱글-스크류 압출기에서의 비 뉴톤성 유체의 비등온유동에 대한 수치해석을 위한 유한요서법(FEM)을 개발하였다. FEM 은 준-3차원 물리적 모델을 기본으로 하였고 온도분포해석을 위하여 전진방법을 도입하였 다. 본 연구에서 도입한 수치해석적 모델과 더 간단한 준-2차원 모델을 온도분포의 예측에 관하여 초점을 맞추어서 비교 하였다. 압출기의 내부채넬에서 순환유동이 온도분포에 미치 는 영향의 중요성이 강조되었다. 순환유동의 효과를 파악하기 위하여 유한 요소식에서 순환 과 관련된 대류항을 빼도록 유한요소 프로그램을 개조 하였고 개발된 원래의 프로그램과 개 조된 프로그램의 해석을 수행한후 해석결과 중에서 특히 온도분포에 있어서의 차이점을 비 교하였다. 이로부터 간단한 준-2차원 모델은 순환효과를 무시했기 때문에 속도분포의 예측 에서는 별문제가 없으나 온도분포의 예측에 있어서는 크게 잘못될 수 있음을 알수 있다.

• 대한조선학회논문집
• /
• 제32권4호
• /
• pp.38-47
• /
• 1995

연안해역에서의 해수의 유동을 추정하기 위하여, $\sigma$ 좌표계를 이용한 3차원 해수유동 수치계산 기법을 개발하여, 이를 다양한 해저지형을 가진 정방형의 유체영역 및 인천항부근의 실해역에 적용한 바 있다. 수치해의 수렴성 및 안정성에 있어서 몇가지 제약이 뒤따랐던, 기존의 수치계산기법을 다음과 같은 몇가지 개선을 하여, 수치해의 수렴성 및 안정성을 도모하였다. (1) 개방경계조건으로서 무반사 경계조건을 도입하였다. (2) 시간의 전개에 있어서, 속도장과 수면변위의 추정시각을 교대로 취하였다. (3) 운동방정식중의 이류항의 공간차분을 중앙차분법에서 상류차분법으로 치환하였다. 그 결과, 유체영역내부의 질량 및 운동량보존의 향상이 얻어졌으며, 시간영역에서의 진동이 대부분 억제되었고, 수치해의 수렴성 및 계산정도 등에 있어서 괄목할 만한 향상이 얻어졌다. 또한, 개선된 3차원 해수유동 수치계산기법과 동경대학의 3차원 해수유동 수치계산기법 상호간의 유효성을 검증하기 위하여 공동연구를 수행하였다. 물리좌표계를 이용한 동경대학의 해수유동 추정결과와, $\sigma$ 좌표계를 이용한 울산대학교의 해수유동 추정결과를 비교한 결과, 일부 다소간의 차가 생기는 점을 제외하고는 대부분 상당한 일치를 보이는 유용한 결과를 얻었다.

• EDISON SW 활용 경진대회 논문집
• /
• 제6회(2016년)
• /
• pp.49-54
• /
• 2016

본 연구에서는 유/무인 우주발사체의 재진입모듈에 적용되는 Heat-Shield(열차단막)의 형상에 따른 유동특성을 수치해석으로 분석했다. 재진입모듈이 지구의 대기권을 다시 진입하는 환경(고도 70km)을 기준으로 해석을 진행했다. 열차단막의 형상은 평판, 곡률이 다른 타원으로 나누었으며, 각 형상별 유동특성을 확인했다. 결론적으로 재진입 모듈의 형상에 따라 압력과 속도분포를 기준으로 열적인 분포를 예상했으며, 계산된 항력계수를 비교했다. 단순한 유동으로 열적 분포를 예상한 것에 한계가 있지만 대기권 재진입 모듈의 2차원 설계에 도움이 될 형상기준을 제시했다. 수치해석은 모두 Edison_CFD에서 제공하는 툴을 이용해 수행했으며, 전처리에는 e-mega (structured)을 후처리에는 e-dava를 이용했다. 해석 solver는 '정렬격자기반 2차원 압축성 유동 범용해석 소프트웨어'를 이용했다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2007년도 추계 학술대회논문집
• /
• pp.194-199
• /
• 2007

IB (immersed boundary) method is one of the prominent tool in computational fluid dynamics for the analysis of flows over complex geometries. The IB technique simplyfies the solution procedure by eliminating the requirement of complex body fitted grids and it is also superior in terms of memory requirement. In this study we have developed numerical code (FOTRAN) for the analysis of 2D flow over a cylinder using IB technique. The code is validated by comparing the wake lengths and separation angles given by Guo et. al. We employed fractional-step procedure for solving the Navier-Stokes equations governing the flow and discrete forcing IB technique for imposing boundary conditions. Also we have developed a 3D code for the backward-facing-step flow and flow over a sphere. The reattachment length in backward-facing-step flow was compared with the one given by Nie and Armaly, which has proven the validity of our code.

• 융합신호처리학회논문지
• /
• 제24권2호
• /
• pp.97-102
• /
• 2023

방향타는 조선 분야에서 중요한 역할을 하기 때문에 방향타의 유체역학적 특성을 조사하기 위해 수치 시뮬레이션이 수행되었다. 유체역학적 힘과 같은 일부 값은 예인 탱크에서 쉽게 측정할 수 있지만, 실험을 통해 압력 분포, 속도 분포, 와류 발생과 같은 유동장에 대한 자세한 정보를 얻기는 어렵다. 본 연구에서는 전산유체역학(CFD)을 이용하여 방향타에 작용하는 유체역학 계수와 레이놀즈수가 미치는 영향을 연구하였다. 상용 전산유체역학 프로그램 Ansys Fluent를 이용하여 방향타 주위의 유동 특성을 연구하였고, 이때 k-epsilon 난류 모델이 사용되었다. 먼저 CFD 상용코드를 이용하여 KCS 방향타의 받음각에 따른 항력계수와 양력계수를 구하였다. 둘째, 동일한 조건에서 2차원 양력계수와 항력계수를 3차원 계수와 비교되었다. 셋째, 레이놀즈수가 유체역학적 힘에 미치는 영향이 연구되었다.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 1998년도 유체기계 연구개발 발표회 논문집
• /
• pp.27-32
• /
• 1998

This paper presents the development of a two-dimensional model for investigating the fluid flow in water jet and calculating the velocity and pressure distributions. The mathematical formulation as a standard k-$\epsilon$ model was solved employing a general thermofluid-mechanics computer program, PHOENICS code, which is based on the Semi-Implicit Method Pressure Linked Equations(SIMPLE) algorithm. The developed code was applied to water jet design to determine the nozzle size, and investigated the effect of the change of nozzle location. Calculated results showed that the flow pattern is not changed as the change of nozzle location.

• 기계저널
• /
• 제29권4호
• /
• pp.360-375
• /
• 1989

이 논평에서는 수치해석에 관련된 국내의 난류모델의 연구현황을 살펴보았다. 국내에서 연구가 많이 되었던 혼합거리모델, R, .epsilon. 2차방정식모델, 대수응력모델의 연구현황을 주로 살펴 보았다. 국내에서는 아직 연구인력이 적기 때문에 꼭 필요한 부분이지만 연구가 되지 않은 부 분이 많았으며 현재 연구가 되고 있는 부분도 그 연구량이 적은 형편이다. 앞으로 전산유체역 학에 의해 열유체기계내의 복잡한 난류유동을 해석하여 기계설계와 운전을 개선하는데 활용하기 위해서는 더욱 많은 연구가 요구된다. 그 중에서 특히 더 관심을 갖고 연구될 분야를 열거 해 보면 (1) 벽면부근의 저 레이놀즈 R, .epsilon. 방정식모델과 대수응력모델에 관한 연구 (2) 복잡한 3차원유동을 정확히 해석할 난류모델 개발에 관한 연구 (3) 복잡한 난류유동에서 열전 달모델에 관한 연구 (4) 레이놀즈응력모델에 관한 연구 (5) Large Eddy Simulation에 관한 연구를 들 수 있다. 이들 연구에서는 난류모델의 물리적 의미에 관한 이론적 연구와 이 모델들 을 여러 경우의 실제유동에 적용하여 검증하는 연구가 병행되어야 할 것이다.

• 한국전산유체공학회지
• /
• 제3권2호
• /
• pp.9-16
• /
• 1998

Predictive behaviors by the extended k-${\varepsilon}$ turbulence model and the standard k-${\varepsilon}$ turbulence model are compared. Grid dependency is tested with the H-type grid as well as the O-type grid. Computations have been performed on a circular-to-rectangular transition duct. The Reynolds number is 390,000 based on the bulk velocity at the inlet. The computed axial velocity contours, transverse velocity profiles, static pressure contours, peripheral skin friction coefficient, peripheral wall static pressure distributions and turbulence kinetic energy have been compared with experimental results. The computed results than those obtained with the standard k-${\varepsilon}$ turbulence model. Comparing to the computed results obtained with the H-type grid and O-type grid, those with H-type grid seem to agree well with experimental results.

• 한국전산유체공학회지
• /
• 제7권3호
• /
• pp.53-59
• /
• 2002

A two-dimensional laminar flow through a channel with a pair of symmetric vertical fins is investigated. At far up- and down-stream from the fins, the plane Poiseuille flow exists in the channel. The Stokes flow for this channel is first investigated analytically and then the other laminar flows by numerical method. For analytic method, the method of eigen function expansion and collocation method are employed. In numerical solution for laminar flows, finite difference method(FDM) is used to obtain vorticity and stream function. From the results, the streamline patterns are shown and the additional pressure drop due to the attached fins and the force exerted on the fin are calculated. It is clear that the force depends on the length of fins and Reynolds number. When the Reynolds number exceeds a critical value, the flow becomes asymmetric. This critical Reynolds number Re/sub c/ depends on the length of the fins.