• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1999.10a
• /
• pp.18-18
• /
• 1999

IRR형태의 액체 램제트 추진기관의 공기 흡입구 유동과 내부 연소 유동을 파악하기 위한 수치적 해석을 수행하였다. 해석은 다원 혼합기체에 대한 압축성 Navier-Stoke 방정식과 공기/Kerosene에 대한 화학 반응을 고려하였으며, 결합된 형태의 k-$\omega$/k-$\varepsilon$ 2 방정식 난류모델을 이용하였다. 기본 유동 해법으로는 고차의 시간 및 공간 정확도를 가지는 근사 Riemann 해법과 LU-SGS 방법을 이용하였다.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.2 no.1
• /
• pp.9-25
• /
• 1994

In the present study, the Reynolds-averaged Navier-Stokes equations, together with the equations of the $k-{\varepsilon}$ model of turbulence, were solved numerically in a general body-fitted coordinate system for three-dimensional turbulent flows around the six basic shapes of the magnetically levitated train(MAGLEV). The numerical computations were conducted on the MAGLEV model configurations to provide information on shapes of this type very near the elevated track at a constant Reynolds number of $1.48{\times}10^{6}$ based on the body length. The coordinate system was generated by numerically solving a set of Poisson equations. The convective transport equations were discretized using the finite-analytic scheme which employed analytic solutions of the locally-linearized equations. A time marching algorithm was employed to enable future extensions to be made to handle unsteady and fully-elliptic problems. The pressure-velocity coupling was treated with the SIMPLER-algorithm. Of particular interests were wall effect by the elevated track on the aerodynamic forces and flow characteristics of the six models calculated. The results indicated that the half-circle configuration with extended sides and with smooth curvature of sides was desirable because of the low aerodynamic forces and pitching moment. And it was found that the separation bubble was occured at wake region in near the elevated track.

• Journal of Biosystems Engineering
• /
• v.25 no.4
• /
• pp.279-286
• /
• 2000

Most of the domestic cold storage rooms are inefficient for agricultural products because of temperature gradients inside the storage rooms. Temperature gradients are developed mainly by improper airflow pattern inside the storage room, which is a main cause of the spoilage of the agricultural products. There proper airflow pattern is essential to minimize these temperature gradients and the spoilage. The performance and characteristics of a cold storage room were determined as a function of airflow pattern and temperature distribution in forced circulation cold storage room. A commercial CFD(computational fluid dynamics) code was used to simulate 3-D airflow in the cold storage room. Solving the flow equations for the storage room, a standard k-$\varepsilon$ turbulent model was implemented to calculate steady state turbulent velocity distribution. The CFD prediction results were compared with temperature measurements inside the cold storage room. In case of pallet storage, Temperature gradients inside pallet storage was reduced because the contact area of cold air expanded through an alley of airflow in storage. But is case of bulk storage, the last temperature of storage considerably rose more than the initial temperature of storage. The reason was that bulk storage didn't include any alley of airflow in storage.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.11
• /
• pp.925-933
• /
• 2014

Today, computational fluid dynamics (CFD) is widely used in industry because of the availability of high-performance computers. However, full-scale analysis poses problems owing to the limited resources and time. In this study, the performance and optimal size of a heat exchanger were calculated using the effectiveness-number of transfer units (${\varepsilon}-NTU$) method and a database of characteristics heat exchanger. Information about the geometry and performance of various heat exchangers is collected, and the performance of the heat exchanger is calculated under the given operating conditions. To determine the optimal size of the heat exchanger, a Genetic Algorithm (GA) is used, and MATLAB and REFPROP are used for the calculation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.10
• /
• pp.1457-1463
• /
• 2003

In this paper, the response surface method using three-dimensional Navier-Stokes analysis to optimize the shape of a multi-blade centrifugal fan, is described. For numerical analysis, Reynolds-averaged Navier-Stokes equations with standard k - c turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in this centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. Linear Upwind Differencing Scheme(LUDS) is used to approximate the convection terms in the governing equations. SIMPLEC algorithm is used as a velocity-pressure correction procedure. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. It was found that the optimization process provides reliable design of this kind of fans with reasonable computing time.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.9
• /
• pp.2294-2303
• /
• 1993

This paper investigated the characteristics of the turbulent incompressible flow past the orifice ring in an axi-symmetric pipe. The flow field was the turbulent pulsatile flow for Reynolds number of $2{\times}10^{5}$ which was defined based on the maximum velocity and the pipe diameter at the inlet, with oscillating frequence $(f_{os})=1/4{\pi}$ which was considered as quasi-steady state frequence. In the present investigation, finite analytic method was used to solve the governing equations in Navier Stokes and turbulent transport formulations. Particularly at high Reynolds number and low oscillation frequency, the effects of orifice ring on the flow were numerically investigated. The separation zone behind the orifice ring during the acceleration phase was found to be decreased. However, during the deceleration phase, the separation behind the orifice ring for pulsatile flow continuously grow to a size even larger than that in steady flow. The pressure drop in steady flow was found to be constant and always positive while for pulsatile flow the pressure drop change with time. And large turbulent kinetic energy, dissipation rate were found to be located in the region where the flow passes through the orifics ring. The maximum turbulent kinetic energy, generally occurs along the shear layer where the velocity gradient is large.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.9
• /
• pp.2315-2328
• /
• 1993

The turbulence effect of particle deposition on a horizontal free-standing wafer in a vertical flow has been studied numerically by using the low-Reynolds-number k-.epsilon. turbulence model. For both the upper and lower surfaces of the wafer, predictions are made of the averaged particle deposition velocity and its radial distribution. Thus, it is now possible to obtain local information about the particle deposition on a free-standing wafer. The present result indicates that the particle deposition velocity on the lower surface of wafer is comparable to that on the upper one in the diffusion controlled deposition region in which the particle sizes are smaller than $0.1{\mu}m$. And it is found in this region that, compared to the laminar flow case, the averaged deposition velocity under the turbulent flow is about two times higher, and also that the local deposition velocity at the center of wafer is high equivalent to that the wafer edge.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.8 no.3
• /
• pp.53-60
• /
• 2004

In order to investigate the nozzle - rotor interactions and the effect of partial admission, the flows in supersonic turbine rotor cascades with a nozzle have been computed. Extensive computations of partial admission supersonic turbines provide the shock structures and flow patterns in the nozzle and rotor. The governing equations were discretized with Euler implicit method in time and 2nd-order upwind scheme of FVM in space. The $\kappa$-$\varepsilon$ turbulence model was utilized to describe the turbulent flow field. It is clearly shown that the nozzle flow is highly affected by the shocks or expansion waves propagated from the rotor leading edge. And the rotor flow is also affected by the shocks or wakes originated from the nozzle.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.10a
• /
• pp.254-259
• /
• 2003

The Ejector is used to get low pressure, and it has been applied to a lot of industry field like the heat engine, the fluid instrument power plant, the food industry, environment industry etc... because there are not any problem even it is mixed with a any kind of liquid, gas, and solid. The flow characteristics in ejector are investigated by PIV and CFD. The experiment using PIV measurement for mixing pipe’s flow characteristics acquired velocity distribution, .Condition : when mixing pipe’s diameter ratio is 1:1.9, and the flux is $Q_{1}=1.136\;l/s$, $Q_{2}=1.706\;l/s$, $Q_{3}=2.276\;l/s$. Based on the PIV and the CFD results, the flow characteristics in ejector are discussed, and it shows the validity of this study.

• Journal of computational fluids engineering
• /
• v.6 no.2
• /
• pp.22-31
• /
• 2001

The two-dimensional unsteady incompressible viscous flow behind rectangular bluff bodies between two parallel plates was analyzed numerically. The steady state flow and the vortex flow behind rectangular bluff bodies in the channel were investigated for two regimes i.e., the laminar(Re = 100, 300, 500) and the turbulent flows(Re = 10⁴∼10/sup 6/). The vortex shedding was generated by a physical disturbance(6%) numerically imposed at the rear of the bluff bodies for a short time. It was observed that the perturbed flow became periodic after a transient period. And in the case of unsteady inflow, the sinusoidal pulsatile flow was applied as the inlet condition in the turbulent flow of Reynolds number of 1.0×10/sup 5/. FLUENT code was employed to solve the problems. The power-law scheme was used to get stable linearized equations and the PISO algorithm was applied to finding the solution of them.