• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.28 no.6
• /
• pp.332-349
• /
• 2016

OLAFOAM is the powerful CFD code and is an expanded version of $OpenFOAM^{(R)}$, for wave mechanics simulation. The $OpenFOAM^{(R)}$ does provide many solvers to correspond to each object of the numerical calculation in a variety of fields. OLAFOAM's governing equation bases on VARANS (Volume-Averaged Reynolds-Averaged Navier-Stokes) equation, and the finite volume method is applied to numerical techniques. The program is coded in C++ and run on the Linux operating system. First of all, in this study, OLAFOAM was validated for 1) wave transformation inside porous structure under bore and regular wave conditions, 2) wave transformation by submerged breakwater under regular wave condition, and 3) regular wave transformation and resultant vertical velocity distribution under current by comparison with existing laboratory measurements. Hereafter, this study, which is almost no examination carried out until now, analyzed closely variation characteristics of water surface level, wave height, frequency spectrum, breaking waves, averaged velocity and turbulent kinetic energy around porous submerged breakwater in the wave and current coexisting field for the case of permeable or impermeable rear beach. It was revealed that the wave height fluctuation according to current direction(following or opposing) was closely related to the turbulent kinetic energy, and others.

• 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.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.25 no.4
• /
• pp.7-12
• /
• 1988

The Resistance Committee of Korea Towing Tank Conference extended the Cooperative Experimental Study Program(1985)[1] to perform the geosim tests by exchanging the five different scaled Series 60, $C_B=0.60$ models between the participating organizations and 13 sets of resistance data have been obtained. The test results are compared among the participating towing tanks and also with the results given in the report of the Resistance and Flow Com-identify the mittee of 18th ITTC. The form factor of each model is derived by Prohaska's method to investigate its dependency on $R_n$. On the other hand, at each $F_n$, form factors are also derived by Telfer's method to relation between $F_n$ and form factor. For this hull form, form factors show relatively weak dependency on $R_n$ and strong dependency on $F_n$. And it is also found that dependencies on both have a cross relation. It seems that further study should be continued to understand more clearly the physical phenomena involved in this problem.

• Solar Energy
• /
• v.13 no.1
• /
• pp.11-21
• /
• 1993

The purpose of this study was to investigate the enhancement of heat transfer without additional external power in the case of rectangular impingement air jet vertically on the flat heating surface. The technique used in the present study was placement of square rod bundles as a turbulence promoter in front of the heat transfer surface. The results obtained through this study were summerized as follws. High heat transfer enhancement was achieved by inserting rods in front of the heating flat plate. According to visulaization, it was examined because of flow acceleration and separation and disturbance of boundary layer. The smaller clerance between rod and heating plate was, the larger heat transfer effect became at each H/B. Arverage Nusselt number reached maximum at H/B=10 and the local augmentation rate of heat transfer became maximum at H/B=2. The maximum average heat transfer enhancement rate increase about 43% for the case of X/B=2 and C=1mm, compared to a flat plate without rods. The correlating equation of average Nusselt number and Reynolds number was obatined. As follws : ${\overline{Nu}}_0=1.249Re^{0.465}(C/A)^{-0.033}(H/B)^{0.013}$.

• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.7
• /
• pp.441-448
• /
• 2016

The thrust force created by a propeller depends on the incoming flow velocity and the rotational velocity of the propeller. The performance of the propeller can be described by dimensionless variables, advanced ratio, thrust coefficient, and power coefficient. This study included the application of the immersed boundary lattice Boltzmann method (IBLBM) with the stereo lithography (STL) file of the rotating object for performance analysis. The immersed boundary method included the addition of the external force term to the LB equation defined by the velocity difference between the lattice points of the propeller and the grid points in the domain. The flow by rotating a 4-blade propeller was simulated with various Reynolds numbers (Re) (including 100, 500 and 1000), with advanced ratios in the range of 0.2~1.4 to verify the suggested method. The typical tendency of the thrust efficiency of the propeller was obtained from the simulation results of different advanced ratios. It was also necessary to keep the maximum mesh size ratio of the propeller surface to a grid size below 3. Additionally, a sufficient length of the downstream region in the domain was maintained to ensure the numerical stability of the higher Re and advanced ratio flow.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.1
• /
• pp.11-23
• /
• 1997

This paper contains a new approach to blade section design method for marine propellers. The hydrodynamic characteristics of 2-D section are highly influenced by its geometrical parameters i.e., thickness and camber distributions and leading edge radius etc. To consider fully turbulent flow field near 2-D section. the finite volume method with k-${\varepsilon}$ turbulent model which solve Reynolds time averaged Navier-Stokes(RANS) equation is applied. In this study, O-type grid system that can provide many calculation points on blade surface is used. The results were compared with those of the experiment of NACA0012 to confirm the accuracy of the developed codes. The goal of this study is the development of a blade section with high efficiency and low drag. To achieve this, we carried out the tests of lift, drag and cavitation characteristics in cavitation tunnel. The results of experiment were compared with numerical results in order to validate the proposed blades design method. By comparing the numerical results with the experiments, we found that the new blade section, KH28 allows superior performance in efficiency and cavitation avoidance characteristics. We further investigated the blade section design method and an application study of this section, KH28 to apply to the marine propeller. In order to improve the accuracy of numerical results on prediction of lift and drag, we conclude here that the 2-layer boundary model must be used.