• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2005년도 춘계 학술대회논문집
• /
• pp.153-158
• /
• 2005

Numerical simulations of two-dimensional steady incompressible lid-driven flow in a square cavity are presented to verify the validity of a new solution code(PowerCFD) with unstructured grids. The code uses the non-staggered(collocated) grid approach which is very popular for incompressible flow analysis because of its numerical efficiency on the curvilinear or unstructured grids. Solutions are obtained for configurations with a Reynolds number as high as 10,000 with both rectangular and hybrid types of unstructured grid mesh. Interesting features of the flow are presented in detail and comparisons are made with benchmark solutions found in the literature. It is found that the code is capable of producing accurately the nature of the lid-driven cavity flow at high Reynolds numbers.

• Journal of Advanced Marine Engineering and Technology
• /
• 제20권5호
• /
• pp.22-29
• /
• 1996

Flow characteristics of two-dimensional closed square cavities near unsteady critical Reynolds numbers were studied numerically at four Reynolds numbers : $8{\times}10^3,\;8.5{\times}10^3,\;9{\times}10^3\;and\;9.5{\times}10^3.$ A convection conservative difference scheme based upon SOLA to maintain the nearly 2nd-order spatial accuracy is adopted on irregular grid formation. Irregular grid number is $80{\times}80$ and its minimum size is about 1/400 of the cavity height(H) and its maximum is about 1/53 H. The result shows that the critical Reynolds number indicating the emergence of flow wnsteadiness is ranging from Re=$8{\times}10^3\;to\;8.5{\times}10^3$ and their flow patterns reveal periodic fluctuation during transient and fully developed stages. But macroscopic flow behavior in terms of instantaneous and time-mean characteristics represent remarkable difference.

• 유변학
• /
• 제10권4호
• /
• pp.217-226
• /
• 1998

단섬유 강화 고분자의 사출성형시 고분자 수지의 유동에 의하여 섬유배향이 필연적으로 일어나며, 섬유배향에 의한 이방성 (anisotropy)은 최종성형물의 품질과 기계적인 특성 등에 많은 영향을 미친다. 사출공정 중에서 충전과정(filling stage)은 섬유배향에 지배적인 역할을 하므로, 충전과정의 유동장을 정확하게 해석하는 것은 매우 중요하다. 형상이 복잡한 캐비티(cavity)와 다중 캐비티에서는 먼저 충전이 완료되어 현탁액(suspension)이 압축을 받고 있는 영역들이 존재하게 된다. 기존의 방법처럼 사출성형의 충전과정을 비압축성 유동으로 가정하면 정확한 유동장을 계산할 수 없다. 본 연구에서는 충전과정에서 압축성을 고려한 혼합 유한요소법/유한차분법을 이용하여 유동장을 계산하였다. 충전이 완료되는 순간에서, 이차배향텐서에 대한 배향변화방정식을 4차 Runge-Kutta 방법을 이용하여 해석함으로써 3차원 섬유배향장을 예측하였다. 충전시간이 다른 4개의 캐비티를 갖는 금형을 제작하여 충전과정에서 압축성 효과를 고려한 수치해석 결과가 실험과 잘 일치함을 보였다. 또한, 압축성과 비압축성 유동장에서 이론적으로 계산된 섬유배향의 차이를 정성적 및 정량적으로 비교하였다.

• 유변학
• /
• 제8권2호
• /
• pp.78-91
• /
• 1996

단섬유 보강 플라스틱 복합재료에 사출성형에서 섬유배향은 금형 충전 공정 중의 유동장에 의해 결정되고 섬유의 배향 상태는 역으로 유동장에 영향을 미친다. 단섬유에 의 한 추가적인 응력을 포함하는 Dinh과 Armstrong의 이방성 구성방정식을 충전유동과 섬유 배향의 연계해석에 도입하였다. 충전유동의 해석은 새로운 압력 지배방정식과 에너지 방정 식을 유한요소법과 유한차분법을 이용하여 풀고 동시에 2차 배향 텐서의 변화방정식을 4차 Runge-kutta 방법으로 풀었다. 섬유의 배향상태를 구한 후에 일방향성 복합재료의 Halpin-Tsai 식과 배향 평균모델을 도입하여 사풀성형품의 이방성 기계적 성질이 예측되었 다. 직사각형 캐비티에서 수치해석결과를 실험결과와 비교하였다. 섬유배향과 유동과의 상호 연계작용을 특히 게이트 근처에서 섬유배향에 영향을 미치며 수치해석 결과는 벽면 근처에 서 유동방향으로 배향하는 shell층을 과대 예측함을 알수 있었는데 이는 배향 텐서 변화 방 정식의 최종근사에서 기인하는 오차로 판단된다. 수정된 복합최종 근사를 바탕으로 예측된 이방성 기계적 성질이 기존의 복합최종 근사에 기초한 예측보다 실험 결과에 정량적으로 보 다 잘 일치하였다.

• 태양에너지
• /
• 제17권3호
• /
• pp.59-66
• /
• 1997

본 연구에서는 PTV기법에 의해 속도를 계측하였으며 계측시스템은 레이저를 이용한 시이트라이트와 이미지보오드를 내장한 퍼스널컴퓨터 그리고 관련 소프트웨어로 구성하였다. 속도벡터는 PTV 시스템에 의해 구하였으며 그 결과를 압력에 관한 프와송방정식에 적용하여 전유동장의 압력을 구하였다. 이때의 경계조건으로서는 벽면에서는 부착조건을 부여하였고 구동류가 흐르는 캐비티의 상부에서는 내부의 속도 값을 외삽하였다. 압력계산에 있어서 격자의 해상도는 $40{\times}40$이며 등간격의 엇갈림격자를 사용하였다. 계측결과는 유동장의 속도 및 압력분포를 잘 나타내었다.

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

Measurements of the leakage flow in the shrouded cavity were performed in axial compressor cascades at $Re=2.6{\times}10^5$. This paper describes the effects of the leakage flow tangential velocity on kinematics of the leakage flow in the shrouded cavity and consequent overall loss and exit flow turning at stator blade row downstream. Flow data and flow visualization images consistently indicate that leakage flow circumferentially migrates 2, 4 and 5 blade passages in the direction of rotation for ${\upsilon}_y/c=0.09$, 0.35 and 0.45, respectively where ${\upsilon}_y$ is the leakage tangential velocity and c is the mainstream velocity. The leakage flow contracts to a jet across the seal-tooth resulting in an increase in the leakage axial velocity-doubling the leakage axial velocity in upstream cavity compared to that in the downstream cavity. Consequently, two flow regions are distinguished before and after the seal-tooth. As increasing the leakage tangential velocity, the overall loss downstream of stator blade row decreases and the exit flow turning in the range of span. from the hub endwall to 15% increases while the decreases in the flow turning from 15% to 30% span is observed.

• 설비공학논문집
• /
• 제7권4호
• /
• pp.622-632
• /
• 1995

The present numerical study is aimed to investigate time-dependent characteristics of a two-dimensional lid-driven square cavity flow of three high Reynolds numbers, $7.5{\times}10^3$, $10^4$ and $3{\times}10^4$. A conservative convection term on irregular grids was adopted by renewing the MAC type difference schemes on regular grids. Relaxation of velocity and pressure is implemented by SOLA algorithm. In case of $Re=7.5{\times}10^3$, flow behavior converges to steady state after a transient period. But for $Re=10^4$, periodic unsteady sinusoidal fluctuation of local velocity and kinetic energy is found and continuous movements of small eddies in the secondary flow regions are also discovered. Random generation of eddies and their active migrating behavior are detected for $Re=3{\times}10^4$, resulting in complete unsteady and non-linear flow characteristics. And, an organized structure similar to a Moffat vortex is also observed from the time-mean flow patterns. Furthermore, a typoon-like vortex(TLV) appears intemittently and rotates along the separation regions and boundary layers.

• 한국해양공학회지
• /
• 제12권1호
• /
• pp.113-119
• /
• 1998

An experimental study was carried out in a three-dimensional cubic cavity driven by 2-dimensional plane Poiseuille flow for three kinds of Reynolds number, $10^4$, 3 $\times$ $10^4$ and 5 $\times$ $10^4$ based on the cavity width and cavity inlet mean flow velcoity. Instant simultaneous velocity vectors at whole field were measured by 2-D PIV system. Laser based illumination and two-frame grey-level cross correlation algorithm are adopted. Severe unsteady flow fluctuation within the cavity are remarkable at above Re = 3 $\times$ $10^4$ Reynolds numbers and sheared mixing layer phenomena are also found at the region where inlet driving Poiseuille flow is collided with the clock-wise rotating main primary vortex at upper center area. Instant velocity profiles reveal that deformed forced vortex formation is observed throughout the separate two areas.

• Journal of Advanced Marine Engineering and Technology
• /
• 제22권6호
• /
• pp.935-941
• /
• 1998

Experiments were carried out for a cubic cavity flow. Contrinuous shear stress is supplied by driven flow for high Reynolds number and three kinds of aspect ratios. Velocity vectors are obtained by PIV and they are used as velocity components for Poisson equation for pressure, Related boundary conditions and no-slip condition at solid wall and the linear velocity extrapolation on the upper side of cavity are well examined for the present study. For calculation of pressure resolution of grid is basically $40{\times}40$ and 2-dimensional uniform mesh using MSC staggered grid is adopted. The flow field within the cavity maintains a forced-vortex formation and almost of the shear stress from the driving inflow is transformed into rotating flow energy and the size of the distorted forced-vortex increases with increment of Reynolds number

• 동력기계공학회지
• /
• 제11권4호
• /
• pp.56-64
• /
• 2007

A numerical simulation is carried out mixed convection in horizontal channel with a heat source from below of rectangular cavity. Finite volume method was employed for the discretization and PISO algorithm was used for calculating pressure term. The parameters governing the problem are the Reynolds number ($10^{-2}{\leq}Re{\leq}50$), the Rayleigh number ($10^3{\leq}Ra{\leq}2.06{\times}10^5$), the Prandtl number ($0.72{\leq}Pr{\leq}909$), the aspect ratio ($0.5{\leq}AR=W/H{\leq}2$) and the angle of inclination ($0^{\circ}{\theta}60^{\circ}$). Mean Nusselt number distributions were obtained and effect of Reynolds number, Rayleigh number and Prandtl number on mixed convection in the horizontal channel with rectangular cavity were investigated.