• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.9
• /
• pp.605-612
• /
• 2016

In this study, numerical simulations were conducted for conjugate heat transfer around ice balls in an encapsulated ice thermal storage system. Four shapes of ice balls were modeled; the default one was a sphere, and the other three shapes were designed to enhance convective heat transfer through the ball surface. The flow around the ball was laminar, for which the Reynolds number was 300, and both forced and natural convections inside and outside the balls were considered. The simulations revealed that the magnitude of convective heat transfer for the different shapes decreased in the following order: bone, dimple, hole, and sphere. For the entire simulation, the maximum difference in the average temperatures of water inside the capsules was found to be $0.9^{\circ}C$. Therefore, it can be said that the effect of ice-ball shape on the performance of the ice thermal storage system is significant, considering that more than 0.3 million balls are used in this system.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.18 no.2
• /
• pp.235-242
• /
• 2017

A driving controller for active destruction protections can be applied to machinery, aerospace and military fields. In particular, this controller can be used to track and attack enemy flying objects through the active control. It is important to ensure reliability of the driving controller since its operation should be kept with precision to the target point. The temperature of the environment where the driving controller is used is about -32 C ~ 50 C (241~323 ). Heat generated in the driving controller should be maintained below a certain threshold (85 C (358 )) to ensure reliability; therefore, the study and analysis of the heat flow characteristics in the driving controller are required. In this research, commercial software Solid-Works Flow Simulation was used for the numerical simulation assuming a low Reynolds number turbulence model and an incompressible viscous flow. The goal of this paper is to design the driving controller safely by analyzing the characteristics of the heat flow inside of the controller composed of chips or boards. Our analysis shows temperature distributions for boards and chips below a certain threshold.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.10
• /
• pp.1139-1146
• /
• 2012

An alternative approach for topology optimization of steady incompressible Navier-Stokes flow problems is presented by using P1 nonconforming finite elements. This study is the extended research of the earlier application of P1 nonconforming elements to topology optimization of Stokes problems. The advantages of the P1 nonconforming elements for topology optimization of incompressible materials based on locking-free property and linear shape functions are investigated if they are also valid in fluid equations with the inertia term. Compared with a mixed finite element formulation, the number of degrees of freedom of P1 nonconforming elements is reduced by using the discrete divergence-free property; the continuity equation of incompressible flow can be imposed by using the penalty method into the momentum equation. The effect of penalty parameters on the solution accuracy and proper bounds will be investigated. While nodes of most quadrilateral nonconforming elements are located at the midpoints of element edges and higher order shape functions are used, the present P1 nonconforming elements have P1, {1, x, y}, shape functions and vertex-wisely defined degrees of freedom. So its implentation is as simple as in the standard bilinear conforming elements. The effectiveness of the proposed formulation is verified by showing examples with various Reynolds numbers.

• Wind and Structures
• /
• v.7 no.5
• /
• pp.317-332
• /
• 2004

Wind tunnel experiments were conducted under highly turbulent and disturbed flow conditions over a solid/perforated plate with a long splitter plate in its plane of symmetry. The effect of varied level of perforation of the normal plate on fluctuating velocities and fluctuating pressures measured across and along the separation bubble was studied. The different perforation levels of the normal plate; that is 0%, 10%, 20%, 30%, 40% and 50% are studied. The Reynolds number based on step height was varied from $4{\times}10^3$ to $1.2{\times}10^4$. The shape and size of the bubble vary with different perforation level of the normal plate that is to say the bubble is reduced both in height and length up to 30% perforation level. For higher perforation of the normal plate, bubble is completely swept out. The peak turbulence value occurs around 0.7 to 0.8 times the reattachment length. The turbulence intensity values are highest for the case of solid normal plate (bleed air is absent) and are lowest for the case of 50% perforation of the normal plate (bleed air is maximum in the present study). From the analysis of data it is observed that $\sqrt{\overline{u^{{\prime}2}}}/(\sqrt{\overline{u^{{\prime}2}}})_{max}$, (the ratio of RMS velocity fluctuation to maximum RMS velocity fluctuation), is uniquely related with dimensionless distance y/Y', (the ratio of distance normal to splitter plate to the distance where RMS velocity fluctuation is half its maximum value) for all the perforated normal plates. It is interesting to note that for 50% perforation of the normal plate, the RMS pressure fluctuation in the flow field gets reduced to around 60% as compared to that for solid normal plate. Analysis of the results show that the ratio [$C^{\prime}_p$ max/$-C_{pb}(1-{\eta})$], where $C^{\prime}_p$ max is the maximum coefficient of fluctuating pressure, $C_{pb}$ is the coefficient of base pressure and ${\eta}$ is the perforation level (ratio of open to total area), for surface RMS pressure fluctuation levels seems to be constant and has value of about 0.22. Similar analysis show that the ratio $[C^{\prime}_p$ max/$-C_{pb}(1-{\eta})]$ for flow field RMS pressure fluctuation levels seems to be constant and has a value of about 0.32.

• Fire Science and Engineering
• /
• v.33 no.4
• /
• pp.35-40
• /
• 2019

In this study, the mass flow rate of the heat release rate equation, which is the major factor of the oxygen consumption method, was analyzed for the fundamental investigation of the cone-calorimeter (5 m length and 0.3 m diameter). The shapes of a completely empty inside, 3 mm pore diameter mesh and pore diameter 10 mm honeycomb with 0.76 porosity were constructed using the cone-calorimeter. To calculate the mass flow rate, four bi-directional probes and thermocouples were installed in a uniform position in the vertical direction of flow. The velocity gradient and flow perturbation were measured from the increase in Reynolds number. As the flow capacity increased, the speed gradient increased in all three shapes relative to the turbulence intensity. In addition, the deviation of extended uncertainty to the mass flow was completely low in the order of empty space, mesh (d_p = 3 mm) and honeycomb (d_p = 10 mm and 𝜖 = 0.76) at the 95% confidence level. The results can be used in designs to improve the flow stability of the cone calorimeter.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.76-85
• /
• 1993

In order to predict the free surface signature of turbulent ship wakes two things are essential; a basic understanding of the mechanism of turbulent vortical flow/free surface interactions and a mathematical model to accurately predict the signature. The goal of the study described here is both to supplement experimental work to obtain basic understanding, as well as to condense this understanding in a model(or models) that captures the essential phenomena and thus allows predictions. To do so we followed two main paths guided by experimental observations. One is full simulations of the flow using the clavier-Stokes equations. The other is a vortex modeling, where the vortical structures of the flows are approximated by idealized structures, an the interaction assumed to be essentially inviscid. These approaches complement each other. Full simulations are only applicable to small scale phenomena, where the system is simple, and the Reynolds number is low. The vortex modeling, on the other hand, cannot represent essentially viscous aspects of the problem such as the effect of contamination gradient. Obviously, the modeling is what may eventually lead to a prediction method; the full simulations-too limited to mimic all but the simplest circumstances-are to aid and support the construction of realistic models. We address two-dimensional aspects of the vortex/free surface interaction first. Secondly we obtain some basic understanding of the interaction process through an experiment and then talk about several three-dimensional problems hoping to develop a successful prediction model.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.48-56
• /
• 2002

The flow characteristics around KRISO 3600TEU container ship model have been experimentally investigated in a circulating water channel. The instantaneous velocity vectors were measured using 2-frame PIV measurement system. The mean velocity fields and turbulent statistics including turbulent kinetic energy and vorticity were obtained by ensemble-averaging 400 instantaneous velocity fields. The free stream velocity was fixed at 0.6m/s and the corresponding Reynolds number was $9{\times}10^5$. The test sections were divided into two regions, three transverse sections of the wake region(Station -0.5767, -1, -3) and five longitudinal sections of the wake((Z/(B/2)=0, 0.1, 0.2, 0.4, 0.6). In the wake region, large-scale longitudinal vortices of nearly same strength are symmetric with respect to the wake centerline and a relatively weak secondary vortex is formed near the waterline. With going downstream, the strength of longitudinal vortex is decreased and the wake region expands.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.7 no.2
• /
• pp.21-27
• /
• 2003

In general, simple fluid added mass method is used for the seismic and vibration analysis of the immersed structure to consider the fluid-structure interaction effect. Actually, the structural response of the immersed structure can be affected by both the fluid added mass and damping caused by the fluid viscosity. These variables appeared as a consistent matrix form with the coupling terms. In this paper, finite element formula for the inviscid fluid case and viscous fluid case are derived from the linearized Navier Stoke's equations. Using the finite element program developed in this paper, the analyses of fluid added mass and damping for the hexagon core structure of the liquid metal reactor are carried out to investigate the effect of fluid viscosity with variation of the fluid gap and Reynolds number. From the analysis results, it is verified that the viscosity significantly affects the fluid added mass and damping as the fluid gap size decrease. From the analysis results of eccentricity effect on the fluid added mass and damping of the concentric cylinders, the fluid added mass increase as the eccentricity increases, however the fluid damping increases only when the eccentricity is very severe.

• Journal of Korea Water Resources Association
• /
• v.52 no.10
• /
• pp.697-706
• /
• 2019

A multiphase flow modeling approach equipped with a hybrid turbulence modeling method is applied to compute the gravity currents in a rectangular channel. The present multiphase solver considers the dense fluid, the less-dense ambient fluid and the air above free surface as three phases with separate flow equations for each phase. The turbulent effect is simulated by the IDDES (improved delayed detach eddy simulation), a hybrid RANS/LES, approach which resolves the turbulent flow away from the wall in the LES mode and models the near wall flow in RANS mode on moderately fine computational meshes. The numerical results show that the present model can successfully reproduce the gravity currents in terms of the propagation speed of the current heads and the emergence of large-scale Kelvin-Helmholtz type interfacial billows and their three dimensional break down into smaller turbulent structures, even on the relatively coarse mesh for wall-modeled RANS computation with low-Reynolds number turbulence model. The present solutions reveal that the modeling approach can capture the large-scale three dimensional behaviors of gravity current head accompanied by the lobe-and-cleft instability at affordable computational resources, which is comparable to the LES results obtained on much fine meshes. It demonstrates that the multiphase modeling method using the hybrid turbulence model can be a promising engineering solver for predicting the physical behaviors of gravity currents in natural environmental configurations.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2018.05a
• /
• pp.71-71
• /
• 2018

Sea pollution accidents have been occurred due to the increase of marine ship traffic. Oil spill from the hull hole induced by tanker collision results in the huge sea pollution. Proper and prompt reaction on such oil spill disaster is needed to minimize the damage. Thru-hull emergency wood plug is typically used to manually close small holes, while it is required to develop some mechanical devices for closing large holes in the hull due to huge fluid pressure. Accurate estimation of oil discharge and velocity from such holes are important to develop proper device to control hull hole damage. High resolution CFD modeling investigation on the configurations of hull hole of 7.5 m initial depth and 30 cm diameter, which was observed in the oil spill accident of the Hebei Sprit off the west coast of Korea in 2007, has been carried out to compute the oil spill velocity distribution in terms of flow depth. Friction loss due to the viscous flow and the discharge coefficient of crude oil with specific gravity SG = 0.85 and viscosity of $4-12cP(mPa{\cdot}s)$ at the temperature of $20^{\circ}C-100^{\circ}C$ are presented in terms of Reynolds number based on the results of high-resolution CFD modeling.