• Transactions of the Korean hydrogen and new energy society
• /
• v.32 no.4
• /
• pp.245-255
• /
• 2021

Various studies about metallic bipolar plates have been conducted to improve fuel cell performance through flow field design optimization. These research works have been mainly focused on fuel cells for vehicle, but not fuel cells for building. In order to reduce the price and volume of fuel cell stacks for building, it is necessary to apply a metallic flow field, In this study, for a metallic flow field applied to a fuel cell for building, the effect of a change in the flow field shape on the performance of a polymer electrolyte membrane fuel cell was confirmed using a model and experiments with a down-sizing single cell. As a result, the flow field using a metal foam outperforms the channel type flow field because it has higher internal differential pressure and higher reactants velocity in gas diffusion layer, resulting in higher water removal and higher oxygen concentration in the catalyst layer than the channel type flow field. This study is expected to contribute to providing basic data for selecting the optimal flow field for the full stack of polymer electrolyte membrane fuel cells for buildings.

• Tunnel and Underground Space
• /
• v.25 no.5
• /
• pp.450-461
• /
• 2015

Extensive numerical experiments have been carried out to investigate effect of block size and fracture geometry on hydraulic characteristics of fractured rock masses based on connected pipe flow in DFN systems. Using two fracture sets, a total of 72 2-D fracture configurations were generated with different combinations of fracture size distribution and deterministic fracture density. The directional block conductivity including the theoretical block conductivity, principal conductivity tensor and average block conductivity for each generated fracture network system were calculated using the 2-D equivalent pipe network method. There exist significant effects of block size, orientation, density and size of fractures in a fractured rock mass on its hydraulic behavior. We have been further verified that it is more difficult to reach the REV size for the fluid flow network with decreasing intersection angle of two fracture sets, fracture plane density and fracture size distribution.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1584-1595
• /
• 1993

Many researches were carried out to estimate heat transfer rate on a circular cylinder in a uniform flow. Various empirical correlations were suggested in the past through expermental studies, however there are considerable discrepancies in the estimated values of heat transfer coefficient. Effects of variable properties of fluid on the heat transfer between a circular cylinder and the external uniform flow were numerically investigated in the present study. The flow and temperature fields were solved using a finite volume method for the uniform flow temperature of 200-900K and the wall temperature of 300-900K. The cold as well as the hot cylinders in the uniform flow of constant temperature were investigated. A unified correlation was obtained for the both cases.

• Journal of Korea Water Resources Association
• /
• v.37 no.6
• /
• pp.461-466
• /
• 2004

This study presents a combined experimental and numerical effort to investigate wave breaking of sinusoidal waves in a surf zone. Numerical predictions are verified by comparing to laboratory measurements. The model solves the Reynolds equations and$textsc{k}$-$\varepsilon$ models for the turbulence analysis. To track the free surface displacement, the volume of fluid method is employed. As the height of incident wave increases, the wave breaking occurs at a closer point of the slope in the numerical model and laboratory experiments with the same depth and period. When a wave breaking occurs, the ratio of wave height becomes larger, with the same wave height and depth, as the period increases.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.9 no.1
• /
• pp.35-45
• /
• 2005

Numerical study is carried out to investigate the effects of chemistry and thermal radiation on the rocket plume flow field at various altitudes. Navier-Stokes equations for compressible flows were solved by a fully-implicit TVD code based on the finite volume method. An infinitely fast chemistry module for hydrocarbon mixture with detailed thermo-chemical properties and a thermal radiation module for optically thick media were incorporated with the fluid dynamics code. The plume flow fields of a kerosene-fueled rocket flying at Mach number zero at sea-level, 1.16 at altitude of 5.06 km and 2.90 at 17.34 km were numerically analyzed. Results showed the plume structures at different altitude conditions with the effects of chemistry and radiation. It is understood that the excess temperature by the chemical reactions in the exhaust gas may not be ignored in the view point of propulsion performance and thermal protection of the rocket base, especially at higher altitude conditions.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.22 no.2
• /
• pp.96-107
• /
• 2018

Flow analysis using computational fluid dynamics was conducted to investigate the effect of the igniter flame caused by the gap between the igniter and the propellant grain in a solid rocket motor. Two propellant grain types were assumed; namely cylinder type (1 mm, 3 mm, and 5 mm gap) and the slot type. The slot type had two igniter hole locations. One was located at the small gap of the propellant grain, and the other one was located at the large gap. In the case of the cylinder type, the pressure in the igniter zone was higher with a thinner gap. Additionally, in the case of the cylinder type, the pressure difference between the igniter installed zone and the free volume was also higher as the gap became lower. The cylinder types were affected by the gap distance, but the slot types were not. Moreover, the results of the slot types were similar to the 5-mm gap case of the cylinder type.

• Journal of the Institute of Convergence Signal Processing
• /
• v.22 no.3
• /
• pp.122-127
• /
• 2021

Flow features near the ship bow such as wave breaking, small scale phenomena have been studied using numerical methods. In this study, the bow shaped wedge was adopted which is from previous paper [1, 2] and the conditions of simulation were Re = 1.64 × 10⁵) and Fr = 2.93. Star CCM+, one of the commercial CFD programs has been used for the simulations. Simulation results such as wave profiles near the ship bow, shape of plunging jet, air entrainment, and wave breaking process have been compared with previous experimental and numerical studies. Overall results showed good agreements with previous studies. Profiles of bow waves showed that overturning jet has been created and broken along the wedge. Plunging wave breaking has been observed along the wedge and four components of plunging wave breaking process were shown. It is confirmed that velocity near the overturing jet significantly increased during plunging wave breaking.

• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.1
• /
• pp.30-45
• /
• 2013

The present work focuses on the unsteady aerodynamics and aeroelastic properties of a small-medium sized wind-turbine blade operating under ideal conditions. A tapered/twisted blade representative of commercial blades used in an experiment setup at the National Renewable Energy Laboratory is considered. The aerodynamic loads are computed using Computational Fluid Dynamics (CFD) techniques. For this purpose, FLUENT$^{(R)}$, a commercial finite-volume code that solves the Navier-Stokes and the Reynolds-Averaged Navier-Stokes (RANS) equations, is used. Turbulence effects in the 2D simulations are modeled using the Wilcox k-w model for validation of the CFD approach. For the 3D aerodynamic simulations, in a first approximation, and considering that the intent is to present a methodology and workflow philosophy more than highly accurate turbulent simulations, the unsteady laminar Navier-Stokes equations were used to determine the unsteady loads acting on the blades. Five different blade pitch angles were considered and their aerodynamic performance compared. The structural dynamics of the flexible wind-turbine blade undergoing significant elastic displacements has been described by a nonlinear flap-lag-torsion slender-beam differential model. The aerodynamic quasi-steady forcing terms needed for the aeroelastic governing equations have been predicted through a strip-theory based on a simple 2D model, and the pertinent aerodynamic coefficients and the distribution over the blade span of the induced velocity derived using CFD. The resulting unsteady hub loads are achieved by a first space integration of the aeroelastic equations by applying the Galerkin's approach and by a time integration using a harmonic balance scheme. Comparison among two- and three- dimensional computations for the unsteady aerodynamic load, the flap, lag and torsional deflections, forces and moments are presented in the paper. Results, discussions and pertinent conclusions are outlined.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.2
• /
• pp.204-213
• /
• 2000

This study deals with jet impingement, which is extensively used in the process industries to achieve intense heating, cooling or drying rates and also widely employed as a test flow for turbulent models due to its complex flow configuration, on a flat plate by numerical methods. In this calculation, the finite volume method was employed to solve the Navier-stokes equation based on the non-orthogonal coordinate with non-staggered variable arrangement. To get a better understanding for the fluid flow and heat transfer characteristics of the turbulent jet impingements, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model was adapted and compared with the experimental data and the result of standard $k-{\varepsilon}$ turbulent model. Numerical calculations were carried out with various flow rates, nozzle to plate distances. In the case of the axisymmetric jet impingement on a flat plate, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model showed better agreement with the experimental data than the standard $k-{\varepsilon}$ turbulent model in the prediction of the mean velocity profiles, the turbulent velocity profiles. the turbulent shear stress and the heat transfer rate. The highest heat transfer rate can be obtained when the impingement occurs within the potential core..

• Journal of Advanced Marine Engineering and Technology
• /
• v.36 no.8
• /
• pp.1091-1096
• /
• 2012

The aerodynamic characteristics of circular cylinder in a channel are studied to make clear the flow feature by solving the Navier-Stokes equation based on the finite volume method with unstructured grids. Reviews are made on with the vorticity, velocity, dynamic pressure, residual and drag, where the Reynolds numbers are 50 and 100. The flows for $Re{\succeq}50$ shows the vortex shedding in the wake, and the result is the same as the case of moving cylinder. The ground effect of flat bottom results in the growth of vortex, being generated in the upper side of the cylinder and elongated in the rear. As the cylinder approaches to wall, for example 0.6, the cylinder plays as a role of blockage to obstruct the flow between the cylinder and wall. The drag coefficients are compared with others' results to confirm the validity of the present numerical simulation.