• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.320-325
• /
• 2010

A multiphase CFD analysis is performed to investigate the effect of near-wall grid for simulating a subcooled boiling flow in vertical tube. The multiphase flow model used in this CFD analysis is the two-fluid model in which liquid(water) and vapor(steam) are considered as continuous and dispersed fluids, respectively. A wall boiling model is also used to simulate the subcooled boiling heat transfer at the heated wall boundary. The diameter and heated length of tube are 0.0154 m and 2 m, respectively. The system pressure in tube is 4.5 MPa and the inlet subcooling is 60 K. The near-wall grid size in the non-dimensional wall unit ($y_{w}^{+}$) was examined from 64 to 172 at the outlet boundary. The CFD calculations predicted the void distributions as well as the liquid and wall temperatures in tube. The predicted axial variations of the void fraction and the wall temperature are compared with the measured ones. The CFD prediction of the wall temperature is shown to slightly depend on the near-wall grid size but the axial void prediction has somewhat large dependency. The CFD prediction was found to show a better agreement with the measured one for the large near-wall grid, e.g., $y_{w}^{+}$ > 100.

• Journal of computational fluids engineering
• /
• v.15 no.3
• /
• pp.24-31
• /
• 2010

boiling flow in vertical tube. The multiphase flow model used in this CFD analysis is the two-fluid model in which liquid(water) and gas(vapour) are considered as continuous and dispersed fluids, respectively. A wall boiling model is also used to simulate the subcooled boiling heat transfer at the heated wall boundary. The diameter and heated length of tube are 0.0154 m and 2 m, respectively. The system pressure in tube is 4.5 MPa and the inlet subcooling is 60 K. The near-wall grid size in the non-dimensional wall unit for lqiuid phase ($y^+_{w,l}$) was examined from 101 to 313 at the outlet boundary. The CFD calculations predicted the void distributions as well as the liquid and wall temperatures in tube. The predicted axial variations of the void fraction and the wall temperature are compared with the measured ones. The CFD prediction of the wall temperature is shown to slightly depend on the near-wall grid size but the axial void prediction has somewhat large dependency. The CFD prediction was found to show a better agreement with the measured one for the large near-wall grid, e.g., $y^+_{w,l}$ > 300 at the tube exit.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.7
• /
• pp.1030-1040
• /
• 1998

One of the main unresolved issues in large-eddy simulation(LES) of wall-bounded turbulent flows is the requirement of high spatial resolution in the near-wall region, especially in the spanwise direction. Such high resolution required in the near-wall region is generally used throughout the computational domain, making simulations of high Reynolds number, complex-geometry flows prohibitive. A grid-embedding strategy using a nonconforming spectral domain-decomposition method is proposed to address this limitation. This method provides an efficient way of clustering grid points in the near-wall region with spectral accuracy. LES of transitional and turbulent channel flow has been performed to evaluate the proposed grid-embedding technique. The computational domain is divided into three subdomains to resolve the near-wall regions in the spanwise direction. Spectral patching collocation methods are used for the grid-embedding and appropriate conditions are suggested for the interface matching. Results of LES using the grid-embedding strategy are promising compared to LES of global spectral method and direct numerical simulation. Overall, the results show that the spectral domain-decomposition grid-embedding technique provides an efficient method for resolving the near-wall region in LES of complex flows of engineering interest, allowing significant savings in the computational CPU and memory.

• Proceedings of the KSR Conference
• /
• 2011.10a
• /
• pp.1594-1600
• /
• 2011

Recently developed FDS5 CFD code has employed a near-wall flow treatment method which is Werner-Wengle wall law provided by NIST(National Institute of Standards and Technology). In this study, the wall law has been verified against DNS(Direct Numerical Simulation) data in the parallel plate. The $y^+$ was kept above 11 to fulfill the near-wall flow requirement in the grid generation. The total grid was $32{\times}32{\times}32$. The boundary condition for inlet and outlet was periodic condition and for both side, symmetric condition was used. The fully developed turbulent flow was generated and Re = 10,700. The simulated results were compared with DNS data. RANS results were also used for verification.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.2
• /
• pp.107-115
• /
• 2009

A Computational study was carried out in order to investigate the aerodynamic characteristics of circular cylinder moving near the wavy wall at a low Reynolds number of 50. Lattice Boltzmann method was used to simulate the flow field and immersed boundary method was combined to represent the moving cylinder and wavy wall regardless of the constructed grid in the domain. The aerodynamics characteristics of the cylinder moving near the wavy wall were represented by the comparing the lifting coefficients with various altitudes (H/D) and wave length and amplitudes of wavy wall. It indicated that the twice of increasing-decreasing variations of lifting coefficient are obtained while the cylinder moves near the wavy wall. The first variation is obtained where the cylinder locates near the peak of the wavy wall. Another variation occurs when the distance to the wavy wall becomes longer after passing the peak. It was also classified that three different patterns of relation between the lifting and drag coefficient of the cylinder. However, the classification is limited to the case of the same order of altitude, amplitude and wave length of the wavy wall.

• Journal of Korean Association for Spatial Structures
• /
• v.15 no.1
• /
• pp.75-84
• /
• 2015

This paper presents the dimensionless wall distance, y+ effect on SST turbulent model for wind turbine blade. The National Renewable Energy Laboratory (NREL) Phase VI wind turbine was used for the study, which the wind tunnel and structural test data has publicly available. The near wall treatment and turbulent characteristics have important role for proper CFD simulation. Most of the CFD development in this area is focused on advanced turbulence model closures including second moment closure models, and so called Low-Reynolds (low-Re) number and two-layer turbulence models. However, in many cases CFD aerodynamic predictions based on these standard models still show a large degree of uncertainty, which can be attributed to the use of the $\epsilon$-equation as the turbulence scale equation and the associated limitations of the near wall treatment. The present paper demonstrates the y+ definition effect on SST (Shear Stress Transport) turbulent model with advanced automatic near wall treatment model and Gamma theta transitional model for transition from lamina to turbulent flow using commercial ANSYS-CFX. In all cases the SST model shows to be superior, as it gives more accurate predictions and is less sensitive to grid variations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.3
• /
• pp.157-164
• /
• 2016

The performance of convective heat transfer in rod bundle flow was experimentally evaluated using a twist-vane spacer grid. A $4{\times}4$ square-arrayed rod bundle was prepared as the test section, with a pitch-to-diameter ratio(P/D) of ~1.35. To check the convective heat transfer performance, the circumferential and longitudinal variations in rod-wall temperatures were measured downstream of the twist-vane spacer grid. In the circumferential measurements, the rod-wall temperature toward the twist-vane tip showed the lowest value, which might be due to the deflected water flow caused by the twist-vane. On the other hand, the wall temperature of the longitudinal measurements near the twist-vane spacer grid decreased dramatically, which implies that the convective heat transfer performance was enhanced. A heat transfer enhancement of ~35 % was achieved near downstream of the twist-vane spacer grid, as compared with the upstream value. Based on the present experimental data, a correlation for predicting the heat transfer performance of a twist-vane spacer grid was proposed.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.1
• /
• pp.37-50
• /
• 2016

A large eddy simulation (LES) of a turbulent channel flow is performed by using the third order low-storage Runge-Kutta method in time and second order finite difference formulation in space with staggered grid at a Reynolds number, $Re_{\tau}=590$ based on the channel half width, ${\delta}$ and wall shear velocity, $u_{\tau}$. To reduce the calculation cost of LES, algebraic wall model (AWM) is applied to approximate the near-wall region. The computation is performed in a domain of $2{\pi}{\delta}{\times}2{\delta}{\times}{\pi}{\delta}$ with $32{\times}20{\times}32$ grid points. Standard Smagorinsky model is used for subgrid-scale (SGS) modeling. Essential turbulence statistics of the flow field are computed and compared with Direct Numerical Simulation (DNS) data and LES data using no wall model. Agreements as well as discrepancies are discussed. The flow structures in the computed flow field have also been discussed and compared with LES data using no wall model.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.15 no.1
• /
• pp.67-82
• /
• 2011

Calculation of accurate wall heat flux for high enthalpy flows requires a dense grid system, which leads to significantly large computational time. A high-order scheme can improve the efficiency of calculation because wall heat flux can be obtained accurately even with a relatively coarse grid system. However, conventional high order schemes have some drawbacks such as oscillations near a discontinuity and instability in multi-dimensional problem. To resolve these problems, enhanced Multi-dimensional Limiting Process(e-MLP) was applied as a high-order scheme. It could provide robust and accurate solutions with high order accuracy in calculation of high enthalpy flows within a short time. We could confirm the efficiency of the high order e-MLP scheme through grid convergence tests with different grid densities in a hypersonic blunt nose problem.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.17-23
• /
• 2011

The turbulent flow and vortex shedding phenomena around pantograph panhead of high speed train were investigated and compared with available experimental data and other simulations. The pantograph head was simplified to be a square-cross-section pillar and assumed to be no interference with other bodies. The Reynolds number (Re) was 22,000. The LES(large eddy simulation) of FDS code was applied to solve the momentum equations and the Wener-Wengle wall model was employed to solve the near wall turbulent flow. Smagorinsky model($C_s$=0.2) was used as SGS(subgrid scale) model. The total grid numbers were about 9 millions and the analyzed domain was divided into 12 multi blocks which were communicated with each other by MPI. The time-averaged mainstream flows were calculated and well compared with experimental data. The phased-averaged quantities had also a good agreement with experimental data. The near-wall turbulence should be carefully treated by wall function or direct resolution to get successful application of LES methods.