• 생물환경조절학회지
• /
• 제3권1호
• /
• pp.10-19
• /
• 1994

Wind load is known to be one of major forces to influence the stability of agricultural structures. General flow fields were calculated to determine flow characteristics over the envelop of the following three types of greenhouses with arched roof : single span, twin span greenhouses, and two single span greenhouses apart 3m inbetween. Pressure coefficients along the envelop of greenhouse were numerically calculated by the k-$\varepsilon$ turbulence model, which lead to determine wind forces on it. Curvilinear coordinate for an arched roof and the upwind scheme were adopted for the study. The calculated pressure coefficients were validated with the avaliable data of Japanese Standard and NGAM Standard. The Magnitude of calculated forces over the envelop was not in good accordance with data except the windward wall. Even tile data of Japanese and NGAM Standard for validation deviated a lot from each other in quantity and quality. Such discrepancy may be attributed to different geometric and/or flow configuration conditions for experiments, or the insenstivity of the k-$\varepsilon$ turbulence model to recirculation flow.

• 설비공학논문집
• /
• 제14권11호
• /
• pp.914-922
• /
• 2002

Thermo-acoustic waves can be thermally generated in a compressible flow field by rapid heating and cooling, and chemical reaction near the boundary walls. This mechanism is very important in the space environment in which natural convection does not exist. Also this may be a significant factor for heat transfer when the fluids are close to the thermodynamic critical point. In this study, the generation and transmission characteristics of thermo-acoustic waves in an air-filled confined domain with two-step pulsed heating are studied numerically. The governing equations are discretized using control volume method, and are solved using PISO algorithm and second-order upwind scheme. For the purpose of stable solution, time step was set to the order of $1\times10_-9s,\;and\;grids\;are\;50\times2000$. Results show that temperature and pressure distributions of fluid near the boundary wall subjected to a rapid heating are increased abruptly, and the induced thermo-acoustic wave propagates through the fluid until it decays due to viscous and heat dissipation. Pressure waves have sharp front shape and decay with a long tail in the case of step heating, but these waves have sharp pin shape in the case of pulsed heating.

• 대한설비공학회:학술대회논문집
• /
• 대한설비공학회 2008년도 하계학술발표대회 논문집
• /
• pp.1366-1371
• /
• 2008

Four simulation models of plastic plate heat exchangers are designed and simulated. The flat plate type heat exchanger is designed as the reference model in order to evaluate how much thermal performance increases. The turbulence promoter type heat exchanger is fabricated with cylindrical-type vortex generators and rib-type turbulence promoters. The corrugate type is obtained from the conventional stainless steel compact heat exchangers, which are called the herringbone-type compact heat exchangers. The dimple type heat exchanger has a number of dimples on its surface. In this study, the flow and heat transfer characteristics of the plastic plate heat exchanger are investigated using numerical simulation and compared with experimental results. The flows are assumed as a three-dimensional, incompressible and turbulent model. The standard k-$\varepsilon$ model is used as the turbulent flow modeling, the SIMPLE algorithm is used to treat the coupling between pressure and velocity, and first order upwind scheme is used for discretization of momentum, turbulent and energy. The computational analysis and experimental results both show that the friction coefficient and Nu number is highest in the corrugate type.

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

The three-dimensional turbulent cascade flows with and without endwall fences are numerically investigated by solving the incompressible Navier-Stokes equations with a high-Reynolds number $k-{\varepsilon}$ turbulence closure model. A projection method based algorithm is used in the finite-volume formulation, with the second order upwind-differencing scheme for the convective terms. First, assessments on accuracy of the present method are made by comparing the static pressure distributions at the mid-span of the cascade with measured data, and also by confirming the experimental observations on the choice of an optimal fence height for the secondary flow control. In understanding the three-dimensional nature of the secondary flow in turbine cascade, the limiting streamline patterns and the static pressure contours at the suction surface of the blade as well as on the cascade endwall are employed to visualize the effectiveness of the endwall fence for the secondary flow control. Analysis on the streamwise vorticity contour maps along the cascade with the three-dimensional representation of their iso-surfaces reveals the strucuture of the complicated vortical flow in the turbine cascade with endwall fence, and also leads to an understanding on formation of the counter-rotating streamwise vortex over the endwall fence, in explaining the mechanisms of controlling the secondary flow and also for the proper selection of an optimal fence height.

• International Journal of Aeronautical and Space Sciences
• /
• 제3권1호
• /
• pp.74-85
• /
• 2002

The numerical experiment has been conducted to investigate the unsteady shock wave reflecting phenomena. The cell-vertex finite-volume, Roe's upwind flux difference splitting method with unstructured grid is implemented to solve unsteady Euler equations. The $4^{th}$-order Runge-Kutta method is applied for time integration. A linear reconstruction of the flux vector using the least-square method is applied to obtain the $2^{nd}$-order accuracy for the spatial derivatives. For a better resolution of the shock wave and slipline, the dynamic grid adaptation technique is adopted. The new concept of grid adaptation technique, which is much simpler than that of conventional techniques, is introduced for the current study. Three error indicators (divergence and curl of velocity, and gradient of density) are used for the grid adaptation procedure. Considering the quality of the solution and the numerical efficiency, the grid adaptation procedure was updated up to $2^{nd}$ level at every 20 time steps. For the convenience of comparison with other experimental and analytical results, the case of interaction between the straight incoming shock wave and a sharp wedge is simulated for various flow conditions. The numerical results show good agreement with other experimental and analytical results, in the shock wave reflecting structure, slipline, and the trajectory of the triple points. Some critical cases show disagreement with the analytical results, but these cases also have been proven to show hysteresis phenomena.

• 대한기계학회논문집
• /
• 제16권3호
• /
• pp.527-539
• /
• 1992

본 연구에서는 터빈익렬의 입구유동면에 주어지는 끝벽 경계층유동에 의하여 익렬 내의 유동에서 발생하는 여러 와류들에 의한 2차 유동과 이와 연관된 여러가지 3차원 점성유동 현상 그리고 이에 따른 유동손실을 보다 정확히 예측하기 위한 수치해 석적 연구를 수행하였으며, 이에 필요한 수치해석적 연구를 수행하였으며, 이에 필요 한 수치해석코드를 작성하였다.유동특성에 대하여 상세한 연구결과가 보고되어 있 는 UTRC(United Technologies Research Center) 평면 터빈익렬을 연구대상으로 채택하 여 익렬 내의 3차원 유동특성을 연구하고 계산한 결과를 기존의 결과와 비교 검토하였 다. 강한 2차유동이 존재하는 경우에 발생하는 수치확산을 감소시키기 위하여 대류 항에 대하여 2차 정확도(second-order accuracy)의 선형상류도식(linear upwind sche- me)을 사용하여 일반적으로 널리 사용되는 하이브리드도식(hybrid scheme)에 의한 해 석결과와 비교하였다. 터빈익렬 내의 난류 유동은 익렬의 회전과 유선의 만곡 등에 의한 영향으로 복잡한 유동현상을 나타내지만, 터빈익렬 내의 난류유동 특성에 대한 실험결과가 아직까지는 부족하고 또한 본 연구에서는 평균유동값의 정확한 해석에 중 점을 두었으므로 표준 k-.epsilon. 모델을 사용하였다.

• 한국전산유체공학회지
• /
• 제19권3호
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

• 한국수자원학회논문집
• /
• 제36권6호
• /
• pp.937-947
• /
• 2003

지진해일 또는 조석과 같은 파의 거동을 수치모의할 때 해안선의 위치는 파랑의 움직임에 따라 끊임없이 이동하게 된다. 따라서, 이를 수치모형에 효과적으로 반영하기 위해서는 수치모형의 경계선을 파랑의 움직임에 따라 이동시켜야 한다. 본 연구에서는 경사지형을 계단지형으로 단순화한 이동경계조건에 대해 해안선의 이동을 보다 정확히 모의하기 위하여 천수방정식의 비선형항을 2차 정확도의 풍상차분기법으로 차분화하여 해석하였다. 개발된 수치모형을 검증하기 위하여 바닥마찰이 없는 것으로 가정된 타원형 수조에서의 유체의 주기적인 거동을 수치해석하였으며, 원형섬에서의 처오름높이를 계산하여 수리실험결과와 비교하였다. 수치해석의 결과는 수리실험 결과와 양호하게 일치하였다.

• 한국항공우주학회지
• /
• 제45권3호
• /
• pp.233-240
• /
• 2017

본 연구에서는 전산유체역학의 특징에 대한 이해를 위해 천음속 날개-동체 주위의 유동장을 In-house 전산유체 코드로 해석하여 시험 결과와 비교하였다. 날개는 RAE 101 익형 단면을 가진 RAE Wing 'A'이며 동체는 축대칭 형상이다. In-house 코드는 비정렬 격자 기반의 압축성 Euler/Navier-Stokes 해석 코드이다. 격자에 대한 의존도, 난류 모형, 공간차분 기법, 점성/비점성의 영향을 시험 결과와 비교하여 살펴보았다. 난류 모형은 $k-{\omega}$ 모형, Spalart-Allmaras 모형, $k-{\omega}$ SST을 적용하였고, 공간차분 기법은 Jameson의 인공 점성를 도입한 중앙 차분 기법과 Roe의 풍상 차분 기법을 적용하였다. 대체적으로 시험 결과를 잘 예측하였으나, 압력분포 및 충격파의 위치가 난류 모형 및 공간 차분 기법에 따라 조금씩 다르게 예측되었으며, 정확한 충격파 위치를 예측하기 위해서는 난류 점성 효과가 고려되어야 함을 알 수 있다.

• 에너지공학
• /
• 제5권1호
• /
• pp.50-55
• /
• 1996

상용 code인 Fluent를 이용하여 도자기의 1차 소성로에 대해 온도장과 속도장을 계산하고 열효율을 산출하였다. 수치해석의 변수로는 배출구 및 도자기의 위치를 변화시켰으며, 수치해석 방법은 검사체적에 기초한 유한차분방법 및 Upwind scheme과 SIMPLEC Algorithm을 사용하였고 난류모델로는 표준 k-$\varepsilon$ 모델을 사용하였다. 계산결과 출구 위치가 전체 소성로내 벡터유동의 양상을 크게 좌우하는 것으로 나타났으며, 전체 온도장에 대해서는 복사의 영향으로 큰 차이를 보이지 않았으나 예열대 상부에 출구가 있는 경우와 비교할 때 예열대 또는 냉각대의 측면에 출구가 있는 경우에 그 영역의 온도가 다소 높게 나타났다. 소성품의 위치는 로내 유동장 및 온도 분포에 크게 영향을 끼치지는 않으나 소성품 내 온도는 그 위치하는 영역의 온도 분포에 크게 영향을 받는 것으로 나타났다. 예열대 상부에 출구가 있는 경우 열효율은 17%로 매우 저조하였고 출구에서의 가스온도는 약 1000 K로 매우 높았다.