• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.5
• /
• pp.604-611
• /
• 2004

Simulation of three-dimensional turbulent flow with LES and DNS lakes much time and expense with currently available computing resources and requires big computing resources especially for high Reynolds number. The emerging alternative to provide the required computing power and working environment is the Grid computing technology. We developed the CFD code which carries out the parallel computing under the Grid environment. We constructed the Grid environment by connecting different PC-cluster systems located at two different institutes of Pusan National University in Busan and KISTI in Daejeon. The specification of PC-cluster located at two different institutes is not uniform. We run our parallelized computer code under the Grid environment and compared its performance with that obtained using the homogeneous computing environment. When we run our code under the Grid environment, the communication time between different computer nodes takes much larger time than the real computation time. Thus the Grid computing requires the highly fast network speed.

• Wind and Structures
• /
• v.29 no.6
• /
• pp.405-416
• /
• 2019

Aerodynamic characteristic of a small scale wind turbine under the influence of an incoming uniform wind field is studied using k-ω Shear Stress Transport turbulence model. Firstly, the lift and drag characteristics of the blade section consisting of S826 airfoil is studied using 2D simulations at a Reynolds number of 1×10⁵. After that, the full turbine including the rotational effects of the blade is simulated using Multiple Reference Frames (MRF) and Sliding Mesh Interface (SMI) numerical techniques. The differences between the two techniques are quantified. It is then followed by a detailed comparison of the turbine's power/thrust output and the associated wake development at three tip speeds ratios (λ = 3, 6, 10). The phenomenon of blockage effect and spatial features of the flow are explained and linked to the turbines power output. Validation of wake profiles patterns at multiple locations downstream is also performed at each λ. The present work aims to evaluate the potential of the numerical methods in reproducing wind tunnel experimental results such that the method can be applied to full-scale turbines operating under realistic conditions in which observation data is scarce or lacking.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.05a
• /
• pp.97-103
• /
• 2002

Vibration isolation technology using an air spring and laminated rubber bearing is widely used because it has excellent vibration isolation characteristics. In the part of that, we usually make use of the self-damped air suing. It is occupied two chambers, restrictor, diaphragm and load plate. Two chambers contain compressed air and the volume of chambers and the area of load plate give a definition of stiffness and load. The restrictor and the volume ratio of two chambers give a definition of damping ratio. The conventional model of restrictor is made of one orifice and it causes turbulent flow in the orifice at the region of large deflection. The stillness of air suing is larger and the damping is lower in the region of large deflection. In the multi-orifice case, the stiffness is similar to air spring with one orifice but damping ratio is larger than conventional air spring. And damping ratio is smaller than conventional air suing in small deflection region. Deflection is small in the region of high frequency so small damping is better than large damping. As a result, we can reduce the storage stiffness of air suing in the wide region of deflection and increase the damping ratio in the region of large deflection. After this, we will try to and the relation of Reynolds Number and Flow Resistance then we are going to make another restrictor for air spring to improve damping ratio and stiffness.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.1
• /
• pp.8-18
• /
• 2015

A numerical simulation for a nonslender BWB UCAV configuration with a rounded leading edge and span of 1.0 m was performed to analyze its aerodynamic characteristics. Numerical results were compared with experimental data obtained at a free stream velocity of 50 m/s and at angles of attack from -4 to $26^{\circ}$. The Reynolds number, based on the mean chord length, is $1.25{\times}106$. 3D multi-block hexahedral grids are used to guarantee good grid quality and to efficiently resolve the boundary layer. Menter's shear stress transport model and two transition models (${\gamma}-Re_{\theta}$ model and ${\gamma}$ model) were used to assess the effect of the laminar/turbulent transition on the flow characteristics. Aerodynamic coefficients, such as drag, lift, and the pitching moment, were compared with experimental data. Drag and lift coefficients of the UCAV were predicted well while the pitching moment coefficient was underpredicted at high angles of attack and influenced strongly by the selected turbulent models. After assessing the pressure distribution, skin friction lines and velocity field around UCAV configuration, it was found that the transition effect should be considered in the prediction of aerodynamic characteristics of vortical flow fields.

• 유체기계공업학회:학술대회논문집
• /
• 2003.12a
• /
• pp.385-390
• /
• 2003

Shape of a multi-blades centrifugal fan is optimized by response surface method based on three-dimensional Navier-Stokes analysis. For numerical analysis, Reynolds-averaged Wavier-Stokes equations with standard $k-{\varepsilon}$ turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in this centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. Optimizations with and without constraints are carried out. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. The correlation of efficiency with relative size of inactive zone at the exit of impeller is discussed as well as with average momentum fluxes in the scroll.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.327-330
• /
• 2006

This paper represents the numerical analysis on effects of radius ratio in a concentric annulus with a rotating inner cylinder. The numerical model consisted of two cylinder which inner cylinder is rotating and outer cylinder is fix, and the axial direction is used the cyclic condition because of the length for axial direction is assumed infinite. The diameter of inner cylinder is assumed 86.8mm, the numerical parameters are angular velocity and radius ratio. Also, the whole walls of numerical model have no-slip and the working fluid is used water at $20^{\circ}C$. The numerical analysis is assumed the transient state to observe the flow variations by time and the 3-D cylindrical coordinate system. The calculation grid adopted a non-constant grid for dense arrangement near the wall side of cylinder, the standard $k-{\omega}$ high Reynolds number model to consider the effect of turbulence flow and wall, the fully implicit method for time term and the quick scheme for momentum equation. The numerical method is compared with the experimental results by Wereley and Lueptow, and the results are very good agreement. As the results, TVF isn't appeared when Re is small because of the initial flow instability is disappear by effect of the centrifugal force and viscosity. The vortex size is from 0.8 to 1.1 for TVF at various $\eta$, and the traveling distance for wavy vortex have the critical traveling distance for each case.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.5 s.248
• /
• pp.446-456
• /
• 2006

The effect of relative position of the stationary turbine blade for the fixed vane has been investigated on blade tip and shroud heat transfer. The local mass transfer coefficients were measured on the tip and shroud fur the blade fixed at six different positions within a pitch. A low speed stationary annular cascade with a single turbine stage was used. The chord length of the tested blade is 150 mm and the mean tip clearance of the blade having flat tip is 2.5% of the blade chord. A naphthalene sublimation technique was used for the detailed mass transfer measurements on the tip and the shroud. The inlet flow Reynolds number based on chord length and incoming flow velocity is fixed to $1.5{\times}10^5$. The results show that the incoming flow condition and heat transfer characteristics significantly change when the relative position of the blade changes. On the tip, the size of high heat/mass transfer region along the pressure side varies in the axial direction and the difference of heat transfer coefficient is up to 40% in the upstream region of the tip because the position of flow reattachment changes. On shroud, the effect of tip leakage vortex on the shroud as well as tip gap entering flow changes as the blade position changes. Thus, significantly different heat transfer patterns are observed with various blade positions and the periodic variation of heat transfer is expected with the blade rotation.

• Wind and Structures
• /
• v.29 no.1
• /
• pp.75-84
• /
• 2019

A steady slot suction near the free-end leading edge of a finite-length square cylinder was used to control its aerodynamic forces and vortex-induced vibration (VIV). The freestream oncoming flow velocity ($U_{\infty}$) was from 3.8 m/s to 12.8 m/s. The width of the tested cylinder d = 40 mm and aspect ratio H/d = 5, where H was the height of the cylinder. The corresponding Reynolds number was from 10,400 to 35,000. The tested suction ratio Q, defined as the ratio of suction velocity ($U_s$) at the slot over the oncoming flow velocity at which the strongest VIV occurs ($U_{\nu}$), ranged from 0 to 3. It was found that the free-end slot suction can effectively attenuate the VIV of a cantilevered square cylinder. In the experiments, the RMS value of the VIV amplitude reduced quickly with Q increasing from 0 to 1, then kept approximately constant for $Q{\geq}1$. The maximum reduction of the VIV occurs at Q = 1, with the vibration amplitude reduced by 92%, relative to the uncontrolled case. Moreover, the overall fluctuation lift of the finite-length square cylinder was also suppressed with the maximum reduction of 87%, which occurred at Q = 1. It was interesting to discover that the free-end shear flow was sensitive to the slot suction near the leading edge. The turbulent kinetic energy (TKE) of the flow over the free end was the highest at Q = 1, which may result in the strongest mixing between the high momentum free-end shear flow and the near wake.

• Nuclear Engineering and Technology
• /
• v.56 no.4
• /
• pp.1310-1319
• /
• 2024

The reactor cavity cooling system (RCCS) is a passive reactor safety system commonly present in the designs of High-Temperature Gas-cooled Reactors (HTGR) that removes heat from the reactor pressure vessel by means of natural convection and radiation. It is one of the factors responsible for ensuring that the reactor does not melt down under any plausible accident scenario. For the simulation of accident scenarios, which are transient phenomena unfolding over a span of up to several days, intermediate fidelity methods and system codes must be employed to limit the models' execution time. These models can quantify radiation heat transfer well, but heat transfer caused by natural convection must be quantified with the use of correlations for the heat transfer coefficient. It is difficult to obtain reliable correlations for HTGR RCCS heat transfer coefficients experimentally due to such a system's size. They could, however, be obtained from high-fidelity steady-state simulations of RCCSs. The Rayleigh number in RCCSs is too high for using a Direct Numerical Simulation (DNS) technique; thus, a Reynolds-Averaged Navier-Stokes (RANS) approach must be employed. There are many RANS models, each performing best under different geometry and fluid flow conditions. To find the most suitable one for simulating an RCCS, the RANS models need to be validated. This work benchmarks various RANS models against three experiments performed on the HTTR RCCS Mockup by the Japanese Atomic Energy Agency (JAEA) in 1993. This facility is a 1/6 scale model of a vessel cooling system (VCS) for the High Temperature Engineering Test Reactor (HTTR), which is operated by JAEA. Multiple RANS models were evaluated on a simplified 2d-axisymmetric geometry. They were found to reproduce the experimental temperature profiles with errors of up to 22% for the lowest temperature benchmark and 15% for the higher temperature benchmarks. The results highlight that the pragmatic turbulence models need to be validated for high Rayleigh natural convection-driven flows and improved accordingly, more publicly available experimental data of RCCS resembling experiments is needed and indicate that a 2d-axisymmetric geometry approximation is likely insufficient to capture all the relevant phenomena in RCCS simulations.

• Journal of Korea Water Resources Association
• /
• v.53 no.1
• /
• pp.1-13
• /
• 2020

The most common approach for computing engineering flow problems at high Reynolds number is still the Reynolds-averaged Navier-Stokes (RANS) computations based on turbulence models with wall functions. The recently developed generalized wall functions blending between the wall-limiting viscous and the outer logarithmic relations ensure a smooth transition of flow quantities across two regions. The performances and convergence properties of widely used turbulence models with wall functions that are applicable for turbulence kinetic energy (TKE), turbulent and specific dissipation rates, and eddy viscosity are presented through a series of near wall flow simulations. The present results show that RNG k-𝜖 model should be carefully applied with small tolerance to get the stable solution when the first grid lies in the buffer layer. The standard k-𝜖 and RNG k-𝜖 models are not sensitive to the selection of wall functions for both TKE and eddy viscosity, while the k-ω SST model should be applied together with kL-wall function for TKE and nutUB-wall functions for eddy viscosity to ensure accurate and stable boundary conditions. The applications to a backward-facing step flow at Re=155,000 reveal that the reattachment length is reasonably well predicted on appropriately refined mesh by all turbulence models, except the standard k-𝜖 model which about 13% underestimates the reattachment length regardless of the grid refinement.