• Journal of Electrochemical Science and Technology
• /
• v.14 no.1
• /
• pp.38-50
• /
• 2023

Proton exchange membrane fuel cell (PEMFC) is a portable and clean power generation device. The structural arrangement of the flow field has a significant influence on the delivery efficiency of PEMFC. In this article, a new bionic flow channel is designed based on the inspiration of a spider shape. The branch channel width and branch corner are studied as the focus, and its simulation is carried out by the method of computational fluid dynamics (CFD). The results show that when channel width/rib width and corner of the branch are 1.5 and 130° , respectively, it is the best numerical combination and the cell comprehensive performance is excellent. The final model using this numerical combination is compared with the traditional flow channel model to verify the advancement of this scheme.

• Ocean Systems Engineering
• /
• v.13 no.4
• /
• pp.325-347
• /
• 2023

The present work investigates the wave generation and propagation in a 2-D wave flume to assess the effect of wave reflection for varying beach slopes by using a numerical tool based on computational fluid dynamics. At first, a numerical wave flume (NWF) is created with different mesh sizes to select the optimum mesh size for time efficient simulation. In addition, different beach slope conditions are introduced such as 1:3, 1:5 and numerical beach at the far end of the NWF to optimize the wave reflection solutions. In addition, several parameters are analysed in order to optimize the solutions. The developed numerical model and its key findings are compared with analytical and experimental surface elevation results and it reveals a good correlation. Finally, the recommended numerical solutions are validated with the experimental findings.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.30 no.1
• /
• pp.108-117
• /
• 2024

In this study, simulations based on computational fluid dynamics were performed for self-propulsion performance prediction of a catamaran that has asymmetrical inside and outside hull form and numerous knuckle lines. In the simulations, the Moving Reference Frame (MRF) or Sliding Mesh (SDM) techniques were used, and the rotation angle of the propeller per time step was different to identify the difference using the analysis technique and condition. The propeller rotation angle used in the MRF technique was 1˚ and those used in the SDM technique were 1˚, 5˚, or 10˚. The torque of the propeller was similar in both the techniques; however, the thrust and resistance of the hull were computed lower when the SDM technique was applied than when the MRF technique was applied, and higher as the rotation angle of the propeller per time step in the SDM technique was smaller in the simulations for several revolutions of the propeller to estimate the self-propulsion condition. The revolutions, thrust, and torque of the propeller in the self-propulsion condition obtained using linear interpolation and the delivered power, wake fraction, thrust deduction factor, and revolutions of the propeller obtained using the full-scale prediction method showed the same trend for both the techniques; however, most of the self-propulsion efficiency showed the opposite trend for these techniques. The accuracy of the propeller wake was low in the simulations when the MRF technique was applied, and slight difference existed in the expression of the wake according to the rotation angle of the propeller per time step when the SDM technique was applied.

• Journal of Advanced Research in Ocean Engineering
• /
• v.3 no.4
• /
• pp.165-173
• /
• 2017

This study aims at validating simulations of the forced and freely vibrating cylinders at Reynolds number of approximately 500 in order to identify the capability of the CFD code, and to establish the analysis process of the vortex-induced vibration (VIV). The direct numerical and large eddy simulations were employed to resolve the various length scales of the vortices, and the morphing technique was used to consider a motion of the circular cylinder. For the forced vibration case, both in- and anti-phase VIV processes were observed regarding the frequency ratio. Namely, when the frequency ratio approaches to unity, the synchronization/lock-in process occurs, leading to substantial increases in drag and lift coefficients. This is strongly linked with the switch in timing of the vortex formation, and this physical tendency is consistent with that of Blackburn and Henderson (J. Fluid Mech., 1999, 385, 255-286) as well as force coefficients. For the free oscillation case, the mass and damping ratio of 50.8 and 0.0024 were considered based on the study of Blackburn et al. (J. Fluid Struct., 2000, 15, 481-488) to allow the direct comparison of simulation results. The simulation results for a peak amplitude of the cylinder and a shedding mode are reasonably comparable to that of Blackburn et al. (2000). Consequently, based on aforementioned results, it can be concluded that numerical methods were successfully validated and the calculation procedure was well established for VIV analysis with reasonable results.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.7
• /
• pp.731-738
• /
• 2012

In this study, the evaluation of the aerodynamic characteristics of the NREL Phase VI Rotor System has been performed, for the 7 m/s upwind case using commercial FEA and CFD tools which are ANSYS Mechanical 12.1 and CFX 12.1. The initial operating conditions of the rotor blade include a $3^{\circ}$ tip pitch angle. A numerical simulation was carried out on only the rotor parts, excluding the tower structure based on the equivalent stiffness model, to consider the aeroelastic effect for the numerical simulation using the loosely coupled 2-way fluid-structure interaction method. The blade root bending moment was monitored in real time to obtain reasonable results. To verify the analysis results, the numerical simulation results were compared with the measurements in the form of the root bending moment and the pressure distributions of the NREL/NASA Ames wind tunnel test.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.5
• /
• pp.282-290
• /
• 2013

The effects of the combination of blade number for forward and after propeller on the propeller shaft forces of a contra-rotating propeller (CRP) system are presented in the paper. The research is performed through the numerical simulations based on the Reynolds-Averaged Navier-Stokes equations (RANS). The simulation results of the present method in open water condition are validated comparing with the experimental data as well as the other numerical simulation results based on the potential method for 4-0-4 CRP (3686+3687A) and 4-0-5 CRP (3686+3849) of DTNSRDC. Two sets of CRP are designed and simulated to study the effect of the combination of blade number in behind-hull condition. One set consists of 3-blade and 4-blade, while the other is 4-blade and 4-blade. A full hull body submerged under the free surface is modeled in the computational domain to simulate directly the wake field of the ship at the propeller plane. From the simulation results, the fluctuations of axial force and moment are dominant in the case of same blade numbers for forward and after propellers, whereas the fluctuations of horizontal and vertical forces and moments are very large in the case of different blade numbers.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.11 no.4
• /
• pp.192-205
• /
• 2009

The unexpected wind over the Mt. Hwawang on 9 February 2009 was deadly when many spectators were watching a traditional event to burn dried grasses and the fire went out of control due to the wind. We analyzed the fatal wind based on wind flow simulations over a digitized complex terrain of the mountain with a localized heating area using a three dimensional computational fluid dynamics model, CFD_NIMR_SNU (Computational Fluid Dynamics_National Institute of Meteorological Research_Seoul National University). Three levels of fire intensity were simulated: no fire, $300^{\circ}C$ and $600^{\circ}C$ of surface temperature at the site on fire. The surface heat accelerated vertical wind speed by as much as $0.7\;m\;s^{-1}$ (for $300^{\circ}C$) and $1.1\;m\;s^{-1}$ (for $600^{\circ}C$) at the center of the fire. Turbulent kinetic energy was increased by the heat itself and by the increased mechanical force, which in turn was generated by the thermal convection. The heating together with the complex terrain and strong boundary wind induced the unexpected high wind conditions with turbulence at the mountain. The CFD_NIMR_SNU model provided valuable analysis data to understand the consequences of the fatal mountain fire. It is suggested that the place of fire was calm at the time of the fire setting due to the elevated terrain of the windward side. The suppression of wind was easily reversed when there was fire, which caused updraft of hot air by the fire and the strong boundary wind. The strong boundary wind in conjunction with the fire event caused the strong turbulence, resulting in many fire casualties. The model can be utilized in turbulence forecasting over a small area due to surface fire in conjunction with a mesoscale weather model to help fire prevention at the field.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.167-173
• /
• 1999

An Euler-based CFD tool has been developed for the performance evaluation of a thrust reverser mounted on a high bypass ratio turbofan engine. The computational domain surrounded by the ground and non-reflection boundary includes the whole nacelle configuration with a deployed thrust reverser. The numerical algorithm is based on the modified Godunovs scheme to allow the second order accuracy in both space and time. The grid system is generated by using eleven multi-blocks, of which the total cell number is 148,400. The thrust reverser is modeled as if it locates at the nacelle simply in all circumferential direction. The existence of a fan and an OGV(Outlet Guide Vane) is simulated by adopting the actuator disk concept, in which predetermined radial distributions of stagnation pressure ratio and adiabatic efficiency coefficient are used for the rotor type disk, and stagnation pressure losses and flow outlet angles for the stator type disk. All boundary conditions including the fan and OGV simulation are treated by Riemann solver. The developed solver is applied to a turbofan engine with a bypass ratio of about 5.7 and the diameter of the fan cowl of 83 inch. The computational results show that the Euler-based inviscid method is very useful and economical to evaluate the performance of thrust reversers.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.43 no.4
• /
• pp.304-310
• /
• 2015

Pintle thrusters use pintle stroke to change nozzle throat area, and this controls thrust. Using MATLAB, one-dimensional simulation has been investigated and the results are compared to those of cold flow tests and computational fluid dynamics for the pintle thruster of Chungnam National University. The prediction based on one-dimensional flow theory shows good agreement with measurements for chamber pressure, but deviates for thrust, partly because of nozzle wall separation. Computational results show that nozzle wall separation occurs at an early stage of nozzle expansion, near the design nozzle throat, for the course of pintle strokes. Empirical thrust prediction incorporates nozzle wall separation, and thus 1-D simulation using empirical thrust prediction showed good results for an early stage of pintle stroke.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.5
• /
• pp.381-389
• /
• 2007

6-DOF simulation program is developed in order to increase the efficiency of the store separation analysis. This S/W is much faster than a method based on CFD(Computational Fluid Dynamics) technology, and allows the simulation of stores with fixed shape as well as with extensible wings, because it uses aerodynamic databases which are prepared beforehand. In this paper, aerodynamic databases of stores are obtained with MSAP(Multi-body Separation Analysis Program), and unsteady damping coefficients are modeled with Missile Datcom. These databases and the 6-DOF simulation program are used to predict the trajectory of an external store, while its wings are being deployed. The analysis results indicate that the safe separations of the store can be achieved not only with the wing fixed but with the wings being deployed.