• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.4
• /
• pp.288-297
• /
• 2009

A serpentine channel geometry often used in a fuel cell has a strong pressure gradient between adjacent channels in specific regions. The pressure gradient helps some amount of reactant gas penetrate through a gas diffusion layer(GDL). As a result, the overall serpentine flow structure is slightly different from the intention of a designer. The purpose of this paper is to examine the effect of serpentine flow structure on current density distribution. By using a commercial code, STAR-CD, a numerical simulation is performed to analyze the fuel cell with high aspect ratio of active area. To increase the accuracy of the numerical simulation, GDL permeabilities are measured with various compressive forces. Three-dimensional flow field and current density distribution are calculated. For the verification of the numerical simulation results, water condensation process in the cathode channel is observed through a transparent bipolar plate. The result of this study shows that the region of relatively low current density corresponds that of dropwise condensation in cathode channels.

• Journal of Positioning, Navigation, and Timing
• /
• v.8 no.4
• /
• pp.201-208
• /
• 2019

Numerical integration is necessary for satellite orbit determination and its prediction. The numerical integration algorithm can be divided into single-step and multi-step method. There are lots of single-step and multi-step methods. However, the Runge-Kutta method in single-step and the Adams method in multi-step are generally used in global navigation satellite system (GNSS) satellite orbit. In this study, 4th and 8th order Runge-Kutta methods and various order of Adams-Bashforth-Moulton methods were used for GLObal NAvigation Satellite System (GLONASS) orbit integration using its broadcast ephemeris and these methods were compared with international GNSS service (IGS) final products for 7days. As a result, the RMSE of Runge-Kutta methods were 3.13m and 4th and 8th order Runge-Kutta results were very close and also 3rd to 9th order Adams-Bashforth-Moulton results. About result of computation time, this study showed that 4th order Runge-Kutta was the fastest. However, in case of 8th order Runge-Kutta, it was faster than 14th order Adams-Bashforth-Moulton but slower than 13th order Adams-Bashforth-Moulton in this study.

• 한국신재생에너지학회:학술대회논문집
• /
• 2007.11a
• /
• pp.147-151
• /
• 2007

A serpentine channel geometry often used in a polymer electrolyte membrane fuel cell has a strong pressure gradient between adjacent channels in specific regions. The pressure gradient helps some amount of reactant gas penetrate through a gas diffusion layer(GDL). As a result, the overall serpentine flow structure is slightly different from intention of a designer. The purpose of this paper is to examine the effect of serpentine flow structure on current density distribution. By using a commercial code, STAR-CD, a numerical simulation is performed to analyze the fuel cell with relatively high aspect ratio active area. To increase the accuracy of the numerical simulation, GDL permeabilities are measured with various compression conditions. Three-dimensional flow field and current density distribution are calculated. For the verification of the numerical simulation results, water condensation process in the cathode channel is observed through a transparent bipolar plate. The result of this study shows that the region of relatively low current density corresponds to that of dropwise condensation in cathode channels.

• Journal of computational fluids engineering
• /
• v.4 no.2
• /
• pp.57-66
• /
• 1999

Numerical Analysis has been done for the supersonic off-design jet flow due to the pressure difference between the jet and the ambient fluid. The difference of pressure generates an oblique shock or an expansion wave at the nozzle exit. The waves reflect repeatedly on the center axis and the sonic surface in the shear layer. The pressure difference is resolved across these reflected waves. In this paper, the axi-symmetric Navier-Stokes equation has been used with the κ-ε turbulence model. The second order TVD scheme with flux limiters, based on the flux vector split with the smooth eigenvalue split, has been used to capture internal shocks and other discontinuities. Numerical calculations have been done to analyze the off-design jet flow due to the pressure difference. The variation of pressure along the flow axis is compared with an experimental result and other numerical result. The characteristics of the interaction between the shock cell and the turbulence mixing layer have been analyzed.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.1631-1636
• /
• 2003

The purpose of this study is to analyze the phenomena of the thermally-induced vibration for the flexible space structure due to abrupt change of radiation heating circumstance using the numerical analyze and experiment test. In order to verify this structure, numerical approaches on the simplified flexible tube were compared with experimental test results at the ground experimental facility In this analyze, it was found that the thermal deformation occurs firstly due to fast radiation heating of flexible structure and then the thermally-induced vibration would be induced due to small periodic change of temperature. According to comparison of numerical and experimental result, in case of no tip mass, the first mode vibration by the numerical analyze was O.78Hz same as that of the experimental result However in case of increase tip-masses of 8g l6g, 50g and 100g, the first modes vibration theoretical analyze were 1.75Hz, 1.3Hz, 0.87Hz and O.73Hz, in decrease trend respectively and those by experimental test were 234Hz, 1.5Hz, O.78Hz and O.78Hz in decrease trend respectively Although using the simpled equation for the estimation, the estimation results were similar to experimental results.

• Journal of Power System Engineering
• /
• v.13 no.3
• /
• pp.27-32
• /
• 2009

Numerical analysis information of a complex discharge-passage will be very useful to improve hydrogen compression system. General information about an internal gas flow is presented by numerical analysis approach. Relating with hydrogen compressing system, which have an important role in hydrogen energy utilization, this should be a useful tool to observe the flow quickly and clearly. Flow characteristic analysis, including pressure and turbulence kinetic energy distribution of hydrogen gas from cylinder going to the chamber of a reciprocating compressor are presented in this paper. Discharge-passage model is designed based on real model of hydrogen compressor. Pressure boundary conditions are applied considering the real condition of operating system. The result shows pressure and turbulence kinetic energy are not distributed uniformly along the passage of the hydrogen compressing system. Path line or particles tracks help to demonstrate flow characteristics inside the passage. The existence of vortices and flow direction can be precisely predicted. Based on this result, the design improvement should be done. Consequently, development of the better hydrogen compressing system will be achieved.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.21 no.3
• /
• pp.502-506
• /
• 2012

Air bearing is characterized by its extremely low friction and cleanliness such that it is widely used especially for spindles with ultra-high rotational speed at several tens of thousands rpm. This paper contributes to design of a static radial air bearing suggesting numerical analysis to anticipate its performances. The numerical analysis is an iteration method based on finite difference formulation of the Reynolds equation. A prototype air bearing has been designed and manufactured. Its load capacity has been measured and compared with the numerical solutions. The result shows good consistency between the experiment and theory, which informs that the numerical analysis can be used as an useful tool to anticipate the performances. Effects of design variables on the bearing performance have been examined by Taguchi's experimental methods using orthogonal array. Number of holes for supplying pressurized air, clearance between shaft and bearing, the hole diameter and bearing length are chosen for the design variables. The result shows that the clearance and the bearing length are the most influential variables while the others can be considered as almost negligible.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.2 no.2
• /
• pp.112-118
• /
• 2010

The numerical study about the vortex-induced vibration and vortex shedding in the wake has been presented. Prior to the numerical simulation of flexible riser systems concerning engineering conditions, efficiency validating of the proposed FSI solution method have been performed. The comparison between numerical simulation and published experimental data shows that the CFD method designed for FSI solution could give acceptable result for the VIV prediction of flexible riser/pipe system. As meaningful study on VIV and vortex shedding mode with the focus on flexible riser model systems, two kinds of typical simulation cases have been carried out. One was related to the simulation of vortex visualization in the wake for a riser model subject to forced oscillation, and another was related to the simulation of fluid-structure interaction between the pipes of coupled multi-assembled riser system. The result from forced oscillation simulation shows that the vortex-induced vibration with high response frequency but small instantaneous vibration amplitude contributes to vortex conformation as much as the forced oscillation with large normalized amplitude does, when the frequency of forced oscillation was relatively high. In the multi-assembled riser systems, it has been found that the external current velocity and the distance between two pipes are the critical factors to determine the vibration state and the steady vibration state emerging in quad-pipe system may be destroyed more easily than dual-pipe system.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2009.09a
• /
• pp.482-491
• /
• 2009

In this study, the mechanical behavior of very soft ground that is reinforced on the surface has been investigated with the aid of a series of numerical analyses. Key material properties of each dredged soft ground, reinforcement and backfill sand mat have been parametrically estimated in the numerical analysis. Along with the result of the study previously performed, a series of in-situ loading conditions and settlement exerted by surface reinforcing operation by construction vehicles has been numerically simulated. These result have been used to evaluate the limit bearing capacity for the unreinforced and reinforced soft ground. Also, the results of the numerical analysis obtained in this research were compared with Yamanouchi's empirical correlation for the limit bearing capacity. Engineering charts listed in this paper for estimating the limit bearing capacity provide field engineers with preliminary design tool for surface reinforcement of very soft ground.

• Structural Engineering and Mechanics
• /
• v.62 no.2
• /
• pp.151-162
• /
• 2017

Two $Chang-{\alpha}$ dissipative family methods and two $KR-{\alpha}$ family methods were developed for time integration recently. Although the four family methods are in the category of the dissipative structure-dependent integration methods, their performances may be drastically different due to the detrimental property of weak instability or overshoot for the two $KR-{\alpha}$ family methods. This weak instability or overshoot will result in an adverse overshooting behavior or even numerical instability. In general, the four family methods can possess very similar numerical properties, such as unconditional stability, second-order accuracy, explicit formulation and controllable numerical damping. However, the two $KR-{\alpha}$ family methods are found to possess a weak instability property or overshoot in the high frequency responses to any nonzero initial conditions and thus this property will hinder them from practical applications. Whereas, the two $Chang-{\alpha}$ dissipative family methods have no such an adverse property. As a result, the performances of the two $Chang-{\alpha}$ dissipative family methods are much better than for the two $KR-{\alpha}$ family methods. Analytical assessments of all the four family methods are conducted in this work and numerical examples are used to confirm the analytical predictions.