• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.3
• /
• pp.49-55
• /
• 2017

In this study, gas turbine engine inlet duct was designed to satisfy uniform flow at aerodynamic interface plane (AIP). Haack-series was selected as nose cone profile and duct outer radius($r_o$) was designed to satisfy to match with area change rate between the nose cone and outer duct wall by the 1-D sizing. The design object of the inlet duct wall profile which has the gradual area change rate was uniform Mach number in the core flow region and minimum boundary later thickness at the both inner nose wall and outer duct wall. The flow characteristics inside the inlet duct was evaluated using CFD. The static pressure distribution at the AIP showed uniform pattern within 0.16%. Based on Mach number profile, the boundary layer thickness was 2% of channel height. Kiel temperature rake location was decided less than 100 mm in front of nose cone where the Mach number is less than 0.1 in order to maximize the temperature probe recovery rate.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.7
• /
• pp.719-725
• /
• 2014

In the present study, a numerical analysis of an injection molding process was conducted for predicting the mold deformation considering non-Newtonian flow, heat transfer, and structural behavior. The accurate prediction of mold deformation during the filling stage is important to successfully design and manufacture a precision injection mold. While the local mold deformation can be caused by various factors, a pressure induced by the polymer melt is considered to be one of the most significant ones. In this regard, the numerical simulation considering both the melt filling and the mold deformation was carried out. A mold core for a 2D axisymmetric center-gated disk was used for the demonstration of the present study. The flow behavior inside the mold cavity and temperature distribution were analyzed along with the core displacement. Also, a Taguchi method was employed to investigate the influence of the relevant parameters including flow velocity, mold core temperature, and melt temperature.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.9
• /
• pp.799-805
• /
• 2013

An analytical study is conducted to assess the efficiency of a flat-plate solar collector using nanofluids. The nondimensionalized 2D heat diffusion equation is solved by assuming a wavelength-independent extinction coefficient and intensity to obtain the analytical solution of the temperature distribution in the flat-plate solar collector. The dimensionless temperature distribution is investigated as functions of the volume fraction of the nanofluids, magnitude of heat loss, and collector's depth based on the analytical solution when using water-based single-walled carbon nanohorn (SWCNH) nanofluids as a working fluid. Finally, the efficiency of the flat-plate solar collector using the nanofluids is predicted and compared with that of the conventional solar collector. The results indicate that the efficiency of the nanofluid solar collector is better than that of the conventional solar collector under specific geometrical conditions.

• Journal of Power System Engineering
• /
• v.11 no.4
• /
• pp.18-25
• /
• 2007

제 2세대 직접분사식 가솔린 기관에서 6공 연료분사기의 연료분무특성을 관찰하였다. 실험에 사용한 직접분사식 가솔린 기관은 2개의 흡입밸브와 2개의 배기밸브를 갖는 텀블형 Spray Guided 연소실과 Quartz로 제작된 실린더 라이너와 실린더 헤드 창으로 구성되어 있다. 선회유동을 유도하기 위하여 흡입매니폴드에 선회유동 제어밸브를 부착하였다. 2차원 Mie 스캐터링 기법을 이용하여 연료분사시기, 연료분사압력과 실린더 내 유동 및 냉각수 온도가 연료분무에 미치는 영향을 관찰하였다. 실험결과로는 흡기과정동안 흡기 선회유동은 분사된 연료의 공간적 분포에 크게 작용하였고, 압축과정동안에는 텀블 및 선회유동의 영향이 흡기과정에 비해 크지 않음을 확인하였다. 또한 성층연소를 위해서 압축과정에서 연료를 분사하는 경우 고압의 연료분사압은 분무도달거리의 성장을 촉진시키나 상승하는 피스톤과 이로 인한 실린더 압력의 상승으로 분무도달거리의 성장이 억제됨을 확인할 수 있었다.

• Journal of Korean Society of Environmental Engineers
• /
• v.31 no.10
• /
• pp.933-940
• /
• 2009

To evaluate characteristics in spray flame, laminar counterflow is investigated on the effects of equivalence ratio and fuel by a two-dimensional DNS (direct numerical simulation). For the gaseous phase, Eulerian mass, momentum, energy, and species conservation equations are solved. For the disperse phase, all individual droplets are calculated by the Lagrangian method without the parcel model. n-Decane ($C_{10}H_{22}$) and n-heptane ($C_7H_{16}$) is used as a liquid spray fuel, and a one-step global reaction is employed for the combustion reaction model. As equivalence ratio increases, the fuel ignites early and the high temperature region spreads wider. The peak value of temperature, however, tends to once increase and then decreases with increasing equivalence ratio. The decrease in the peak value of temperature for the higher equivalence ratio condition is caused by the cooling effect associated with droplet group combustion. Since the evaporation of n-heptane is early, the high temperature region spreads wider than ndecane, but the peak values of temperature for both n-heptane and n-decane is almost same.

• Nuclear Engineering and Technology
• /
• v.25 no.1
• /
• pp.166-177
• /
• 1993

A simplified thermal analysis method to evaluate the maximum temperature of the CANDU 37-element fuel bundle within a fuel basket in a given spent fuel dry storage canister has been presented along with the results of sample analyses performed to examine the parametric effects of the ambient conditions on the maximum fuel temperature within a canister. To solve the multi-dimensional heat transfer problem of the complex geometry of rod bundles within a canister where three modes of heat transfer are superimposed, the CANDU spent fuel bundles stored in the dry storage canister are first replaced by equivalent concentric fuel cylinders. The simplified axi-symmetric two-dimensional multi-mode heat transfer problem of the equivalent fuel cylinders is then analyzed with an existing computer code, HEATING5, using additional input data and heat transfer correlations. A comparison between the predicted temperature profile and the mock-up test results shows that the agreement is quite satisfactory.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.5
• /
• pp.952-958
• /
• 1988

To analyze the nonflatness and residual stress in accelerated cooled plate, a numerical analysis model has been developed. Two factors, i.e. temperature and phase transformation, are considered in calculating the stress distribution that develops during cooling. The plastic strain and plate-buckling, which are often shown in accelerated cooled plate, were determined from this stress. Mean temperature in through thickness direction and temperature difference in width direction are considered in the model to simplify the calculation. The temperature and stress distribution changes caused by phase transformation are involved in terms of the effective specific heat and the effective thermal expansion coefficient. With the model, accelerated cooling of 10mm(t) $^{*}$3000mm(w) plate was simulated. The condition of accelerated cooling was .deg. C/sec from just after hot rolling to 500.deg. C. The initial temperature-difference ratio, .DELTA.Tr, in width direction is an important factor in evaluating the stress distribution. When .DELTA.Tr is 0.08, buckling occurs during cooling and 7kgf/m $m^{2}$ of residual stress develops at the edge of plate. To secure the flatness, .DELTA.Tr should be less than 0.07. Small scaled cooling test was conducted to verify the exactness of the model and the results proved the usefulness of this numerical analysis model.l.

• Nuclear Engineering and Technology
• /
• v.14 no.2
• /
• pp.51-63
• /
• 1982

Stability of a PWR core against xenon-induced axial power oscillation is studied using one-dimensional xenon trausient analysis code, DD1D, that has been developed and verified at KAERI. Analyzed by DD1D utilizing the Kori Unit 1 design and operating data is the sensitivity of axial stability in a PWR core to the changes in core physical parameters including core power level, moderator temperature coefficient, core inlet temperature, doppler power coefficient and core average turnup. Through the sensitivity study the Kori Unit 1 core is found to be stable against axial xenon oscillation at the beginning of cycle 1. But, it becomes less stable as turnup progresses, and unstable at the end of the cycle. Such a decrease in stability is mainly due to combined effect of changes in axial power distribution, moderator temperature coefficient and doppler power coefficient as core turnup progresses. It is concluded from the stability analysis of the Kori Unit 1 core that design of a large PWR with high power density and increased dimension can not avoid xenon-induced axial power instabilities to some extents, especially at the end of cycle.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.2 no.2
• /
• pp.59-66
• /
• 1990

Various types of ice distribution under low temperature greatly influence the environment of the Arctic and Antarctic Oceans. To understand fundamentals of ice properties such as Polar glaciers, icebergs and sea ice, this study focuses on the material behaviors and failure mechanisms of polycrystalline ice. Utilizing the continuum damage theory, a three-dimensional constitutive model to describe creep deformation characteristics in the glacial flow is developed in consideration of micro-cracking as the major physical process of ice deformation. The numerical model is compared with the published experimental data especially in uniaxial constant stress creep tests. The model can simulate primary and secondary creeps as well as tertiary creep characteristics due to the microcrack accumulation.

• Fire Science and Engineering
• /
• v.15 no.1
• /
• pp.7-15
• /
• 2001

This paper describes the smoke filling process of a fire field model based on a self-developed SMEP (Smoke Movement Estimating Program) codo to the simulation of fire induced flows in the atrium space (SIVANS atrium at Japan) containing smoke radiation effect. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k-$\varepsilon$ turbulence model with buoyancy term. Also it solves the radiation equation using the discrete ordinates method. The result of the calculated smoke temperature containing radiation effect has shown a better prediction than the result calculated by only convection effect in comparison with the experimental data. This seems to come from the radiation effect of $H_2$O and $CO_2$ gas under smoke productions. Thus, the consideration of the radiation effect under smoke in fire should be necessary in order to get more realistic result. Also the numerical results indicated that the smoke layer is developing at a rate of about 0.1 m/s. It would take about 450 seconds after starting the ultra fast fire of 560 kW that the smoke layer move down to 1.5m above the escape level.