• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.7
• /
• pp.849-861
• /
• 1997

In this study, the effect of hole geometry of the cooling system on the flow and temperature field was numerically calculated. The finite volume method was employed to discretize the governing equation based on the non-orthogonal coordinate with non-staggered variable arrangement. The standard k-.epsilon. turbulence model was used and also the predicted results were compared with the experimental data to validate numerical modeling. The predicted results showed good agreement in all cases. To analyze the effect of the discharge coefficient for slots of different length to width, the inlet chamfering and radiusing holes were considered. The discharge coefficient was increased with increment of the chamfering ratio, radiusing ratio and slot length to width and also the effect of radiusing showed better result than chamfering in all cases. In order to analyze the difference between the predicted results with plenum region and without plenum region, the velocity profiles of jet exit region for a various flow conditions were calculated. The normal velocity components of jet exit showed big difference for the low slot length to width and high blowing rate cases. To analyze the flow phenomena injected from a row of inclined holes in a real turbine blade, three dimensional flow and temperature distribution of the region including plenum, hole and cross stream with flow conditions were numerically calculated. The results have shown three-dimensional flow characteristics, such as the development of counter rotating vortices, jetting effect and low momentum region within the hole in addition to counter rotating vortex structure in the cross stream.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.12
• /
• pp.4013-4026
• /
• 1996

A numerical simulation has been carried out for the jet impinging on a flat plate and a semi-circular concave surface. In this computation finite volume method was employed to solve the full Navier-Stokes equation based on a non-orthogonal coordinate with non staggered variable arrangement. The standard k-.epsilon. turbulent model and low Reynolds number k-.epsilon. model(Launder-Sharmar model) with Yap's correction were adapted. The accuracy of the numerical calculations were compared with various experimental data reported in the literature and showed good predictions of centerline velocity decay, wall pressure distribution and skin friction. For the jet impingement on a semi-circular concave surface, potential core length was calculated for two different nozzle(round edged nozzle and rectangular edged nozzle) to consider effects of the nozzle shape. The result showed that round edged nozzle had longer potential core length than rectangular edged nozzle for the same condition. Heat transfer rate along the concave surface with constant heat flux was calculated for various nozzle exit to surface distance(H/B) in the condition of same jet velocity. The maximum local Nusselt number at the stagnation point occurred at H/B = 8 where the centerline turbulent intensity had maximum value. The predicted Nusselt number showed good agreement with the experimental data at the stagnation point. However heat transfer predictions along the downstream were underestimated. This results suggest that the improved turbulence modeling is required.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.9
• /
• pp.1183-1192
• /
• 2001

The numerical simulation has been conducted for the investigation of flow and heat transfer characteristics of an oblique impingement jet injected to a flat plate. The finite volume method was used to discretize the governing equations based on the non-orthogonal coordinate with non-staggered variable arrangement. The $textsc{k}$-$\varepsilon$-ν(sup)'2 turbulence model was employed to consider the consider the anisotropic flow characteristics generated by the impingement jet flow. The predicted results were compared with the experimental data and those of the standard $textsc{k}$-$\varepsilon$ turbulence model. The results of the $textsc{k}$-$\varepsilon$-ν(sup)'2 model showed better agreement with the experimental data than those of the standard $textsc{k}$-$\varepsilon$ model. In order to get the optimum condition, the flow and temperature fields were calculated with a variation of inclined angle($\alpha$=30$^{\circ}$~90$^{\circ}$) and the distance between the jet exit and impingement plate-to-diameter (L/D=4~10) at a fixed Reynolds number(Re=20,000). For a small L/D, the near-peak Nusselt numbers were not significantly effected by the inclined angle. The near-peak Nusselt numbers were not significantly affected by the L/D in the case of a large $\alpha$. The overall shape of the local Nusselt numbers was influenced by both the jet orifice-to-plate spacing and the jet angle.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.2
• /
• pp.70-76
• /
• 2000

A numerical study on the viscous flow around a 2-dimensional circulation control foil is carried out for application on the field of naval architecture and ocean engineering. The governing equations are the RANS and the continuity equations. The equations are discretized by finite difference method and MAC method and the pressure poisson equation is calculate by a SOR method and an O-type non-staggered boundary fitted coordinate system which is overlapped near the slot is used to improve the numerical accuracy. Turbulence is approximated by a modified Baldwin-Lomax turbulence model. In the present paper, the Coanda effect on a 2-dimensional foil of a 20% thickness ellipse with modified rounded trailing edge has been numerically studied. The change in drag and lift of the foil with various jet momentums are calculated and compared to the experimental results to show good agreements.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.7
• /
• pp.703-714
• /
• 2010

This paper describes the design optimization of a rotating rectangular channel with staggered arrays of pin-fins by Kriging metamodeling technique. Two non-dimensional variables, the ratio of the height to the diameter of the pin-fins and the ratio of the spacing between the pin-fins to the diameter of the pin-fins are chosen as the design variables. The objective function that is a linear combination of heat transfer and friction loss related terms with a weighting factor is selected for the optimization. To construct the Kriging model, objective function values at 20 training points generated by Latin hypercube sampling are evaluated by a three-dimensional Reynolds-averaged Navier-Stokes (RANS) analysis method with the SST turbulence model. The Kriging model predicts the objective function value that agrees well with the value calculated by the RANS analysis at the optimum point. The objective function is reduced by 11% by the optimization of the channel.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.2
• /
• pp.79-87
• /
• 2013

In this paper, an Euler equation solver based on a Cartesian-grid method and non-uniform staggered grid system is applied to predict the ship motion response and added resistance in waves. Water, air, and solid domains are identified by a volume-fraction function for each phase and in each cell. For capturing the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with a weighed line interface calculation (WLIC) method. The volume fraction of solid body embedded in a Cartesian-grid system is calculated by a level-set based algorithm, and the body boundary condition is imposed by volume weighted formula. Added resistance is calculated by direct pressure integration on the ship surface. Numerical simulations for a Wigley III hull and an S175 containership in regular waves have been carried out to validate the newly developed code, and the ship motion responses and added resistances are compared with experimental data. For S175 containership, grid convergence test has been conducted to investigate the sensitivity of grid spacing on the motion responses and added resistances.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.181-186
• /
• 2006

An imaging method of seismic sources using time-reversal wave propagation is presented. The method is based on the time-reversal invariance and the spatial reciprocity of the wave equation. Time-reversal wave propagation has been used to image anomalous features of a midium in medical imaging, non destructive testing and waveform tomography. Seismogram is the record whose energy is propagated from the seismic source. If time-reversed seismogram propagates back into the medium, seismic energy is concentrated at the origin time of the event and at the source location. In this work, a staggered-grid finite-difference method of the elastic wave equation is parallelized for 3-D wave propagation simulation. With numerical experiments, we show that the time-reversal imaging will enable us to explore the spatio-temporal history of complex earthquake.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.395-406
• /
• 2011

In this research a numerical simulation method is developed for moving body in free surface flows using fixed staggered rectangular grid system. The non-linear free surface near the body is defined by marker-density method. The body boundary is defined by line segment connecting the points where the body surface and grid line meet. Continuity equation and Navier-Stokes equations are used as governing equations and the equations are coupled with two-step projection method. The velocities and pressures of body boundary and free surface cells are calculated with simultaneous iterative method. To treat a body movement in a fixed grid system, the volume displaced by moving body is added to the divergence of the body boundary cell. For the verification of the present numerical method. vortex shedding period of advancing cylinder is calculated and the period is compared with existing experiment results. Moreover, added mass and damping coefficients of a vertically excited box are calculated and the computed results are compared with published experiment results. Impulsive pressure and water level variation due to sloshing phenomenon are simulated and the results are compared with published experiment results. Varying the plunger shape, the waves generated by plunging type wave maker are compared with the 2nd order Stokes wave theory The plunger shape generating the wave that shows the best agreement with the theory is represented.

• Journal of Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• Journal of Information Display
• /
• v.4 no.2
• /
• pp.1-6
• /
• 2003

In this paper, we demonstrate that the organic electrophosphorescent device is driven by the organic thin film transistor with spin-coated photoacryl gate insulator. It was found that electrical output characteristics in our organic thin film transistors using the staggered-inverted top-contact structure showed the non-saturated slope in the saturation region and the sub-threshold nonlinearity in the triode region, where we obtained the maximum power luminance that was about 90 $cd/m^2$. Field effect mobility, threshold voltage, and on-off current ratio in 0.45 ${\mu}m$ thick gate dielectric layer were 0.17 $cm^2/Vs$, -7 V, and $10^6$ , respectively. In order to form polyimide as a gate insulator, vapor deposition polymerization process was also introduced instead of spin-coating process, where polyimide film was co-deposited by high-vacuum thermal evaporation from 4,4'-oxydiphthalic anhydride (ODPA) and 4,4'-oxydianiline (ODA) and cured at 150${\sqsubset}$for 1hr. It was also found that field effect mobility, threshold voltage, on-off current ratio, and sub-threshold slope with 0.45 ${\mu}m$ thick gate dielectric films were 0.134 $cm^2/Vs$, -7 V, and $10^6$ A/A, and 1 V/decade, respectively.