• Wind and Structures
• /
• v.13 no.1
• /
• pp.57-82
• /
• 2010

The 486m-long roof of Shenzhen Citizens Centre is one of the world's longest spatial lattice roof structures. A comprehensive numerical study of wind effects on the long-span structure is presented in this paper. The discretizing and synthesizing of random flow generation technique (DSRFG) recently proposed by two of the authors (Huang and Li 2008) was adopted to produce a spatially correlated turbulent inflow field for the simulation study. The distributions and characteristics of wind loads on the roof were numerically evaluated by Computational Fluid Dynamics (CFD) methods, in which Large Eddy Simulation (LES) and Reynolds Averaged Navier-Stokes Equations (RANS) Model were employed. The main objective of this study is to explore a useful approach for estimations of wind effects on complex curved roof by CFD techniques. In parallel with the numerical investigation, simultaneous pressure measurements on the entire roof were made in a boundary layer wind tunnel to determine mean, fluctuating and peak pressure coefficient distributions, and spectra, spatial correlation coefficients and probability characteristics of pressure fluctuations. Numerical results were then compared with these experimentally determined data for validating the numerical methods. The comparative study demonstrated that the LES integrated with the DSRFG technique could provide satisfactory prediction of wind effects on the long-span roof with complex shape, especially on separation zones along leading eaves where the worst negative wind-induced pressures commonly occur. The recommended LES and inflow turbulence generation technique as well as associated numerical treatments are useful for structural engineers to assess wind effects on a long-span roof at its design stage.

• Journal of computational fluids engineering
• /
• v.20 no.2
• /
• pp.23-36
• /
• 2015

The $2^{nd}$ generation ice accretion analysis program has been developed and validated for various icing conditions. The essential feature of the $2^{nd}$ generation code lies in its capability of handling general 3-D geometry and improved accuracy. The entire velocity fields are obtained based on Navier-Stokes equations in order to take the massively separated flow field into account. Unlike $1^{st}$ generation code, the droplet trajectories are calculated using Eulerian approach, which is adopted to yield appropriate collection efficiency even in the shadow region. For improved thermodynamic analysis on the surfaces, water film model and modified Messinger model are newly included in the present analysis. The ice shape for a given time step is obtained by considering the exact amount of ice accreted on the surface. Each module of the icing analysis code has been seamlessly integrated on the OpenFOAM platform. The developed code was validated against available experimental data for 2D airfoils and 3D DLR-F4. Due to the lack of experimental data, the computed results of DLR-F4 were compared with those obtained from FENSAP-ICE, which is state-of-the-art 3D icing analysis code. It was clearly shown that the present code produces comparable results to those of FENSAP-ICE, in terms of prediction accuracy and the capability of handling general 3-D geometries.

• Transactions of the Korean hydrogen and new energy society
• /
• v.31 no.4
• /
• pp.387-394
• /
• 2020

A parametric study of a 2.5 kW class propeller type micro hydraulic turbine was performed. In order to analyze the internal flow characteristics in the hydraulic turbine, three dimensional Reynolds-averaged Navier-Stokes equations with shear stress transport turbulence model were used and the hexahedral grid system was used to construct computational domain. To secure the reliability of the numerical analysis, the grid dependency test was performed using the grid convergence index method based on the Richardson extrapolation, and the grid dependency was removed when about 1.7 million nodes were used. For the parametric study, the axial distance at shroud span (L) between the inlet guide vane and the runner, and the inlet and outlet blade angles (β₁, β₂) of the runner were selected as the geometric parameters. The inlet and outlet angles of the runner were defined in the 3 spans from the hub to tip, and a total of 7 geometric parameters were investigated. It was confirmed that the outlet angles of the runner had the most sensitive effect on the power and efficiency of the micro hydraulic turbine.

• Proceedings of the Korean Society of Computational and Applied Mathematics Conference
• /
• 2003.09a
• /
• pp.1-1
• /
• 2003

In this talk, a reduced-order modeling methodology based on centroidal Voronoi tessellations (CVT's)is introduced. CVT's are special Voronoi tessellations for which the generators of the Voronoi diagram are also the centers of mass (means) of the corresponding Voronoi cells. The discrete data sets, CVT's are closely related to the h-means clustering techniques. Even with the use of good mesh generators, discretization schemes, and solution algorithms, the computational simulation of complex, turbulent, or chaotic systems still remains a formidable endeavor. For example, typical finite element codes may require many thousands of degrees of freedom for the accurate simulation of fluid flows. The situation is even worse for optimization problems for which multiple solutions of the complex state system are usually required or in feedback control problems for which real-time solutions of the complex state system are needed. There hava been many studies devoted to the development, testing, and use of reduced-order models for complex systems such as unsteady fluid flows. The types of reduced-ordered models that we study are those attempt to determine accurate approximate solutions of a complex system using very few degrees of freedom. To do so, such models have to use basis functions that are in some way intimately connected to the problem being approximated. Once a very low-dimensional reduced basis has been determined, one can employ it to solve the complex system by applying, e.g., a Galerkin method. In general, reduced bases are globally supported so that the discrete systems are dense; however, if the reduced basis is of very low dimension, one does not care about the lack of sparsity in the discrete system. A discussion of reduced-ordering modeling for complex systems such as fluid flows is given to provide a context for the application of reduced-order bases. Then, detailed descriptions of CVT-based reduced-order bases and how they can be constructed of complex systems are given. Subsequently, some concrete incompressible flow examples are used to illustrate the construction and use of CVT-based reduced-order bases. The CVT-based reduced-order modeling methodology is shown to be effective for these examples and is also shown to be inexpensive to apply compared to other reduced-order methods.

• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.61-70
• /
• 2016

Dong-tan Station, shared by high-speed railway and urban express railway, is a very complicated underground station having 6 tracks together with barrier and shafts between them, therefore it seems very hard to investigate the aerodynamic effects including the pressure variation and train gust in the station when a high-speed train runs through it. In this study, the aerodynamic effects on the structures and platform passengers when a high-speed train runs through an underground station have been studied using Computational Fluid Dynamics. STAR-CCM+ has been employed for numerical simulation based on Navier-Stokes equation and 2-equation turbulence model and moving mesh scheme supported by STAR-CCM+ has also been used to represent the relative motion between a train and station. Based on the simulation results, the unsteady flow fields in the underground station induced by the high-speed train have been analyzed and the pressures on the PSDs and pressure variation at the platform have quantitatively assessed.

• 한국신재생에너지학회:학술대회논문집
• /
• 2009.06a
• /
• pp.373-373
• /
• 2009

Relative humidity, membrane conductivity and water activity are critical parameters of polymer electrolyte membrane fuel cells (PEMFC) for high performance and reliability. These parameters are closely related with temperature. Moreover, the ideal values of these parameters are not always identical along the channels. Therefore, the cooling channel design and its operating condition should be well optimized along the all location of the channels. In the present study, we have performed a numerical investigation on the effects of cooling channels on performance of a PEMFC. Three-dimensional Navier-Stokes equations are solved with the energy equation including heat generated by the electrochemical reactions in the fuel cell. The present numerical model includes the gas diffusion layers (GDL) and serpentine channels for both anode and cathode gas flows, as well as cooling channels. To accurately predict the water transport across the membrane, the distribution of water content in the membrane is calculated by solving a nonlinear differential equation with a nonlinear coefficient, i.e., the water diffusivity which is a function of water content as well as temperature. Main emphasis is placed on the heat transfer between the solid bipolar plate and coolant flow. The present results show that local current density is affected by cooling channels due to the change of the oxygen concentration and the membrane conductivity as well as the water content. It is also found that the relative humidity is influenced by the generated water and the gas temperature and thus it affects the distribution of fuel concentration and the conductivity of the membrane, ultimately fuel cell performance. Unit-cell experiments are also carried out to validate the numerical models. The performance curves between the models and experiments show reasonable results.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.6
• /
• pp.470-475
• /
• 2017

Since the intake air of gas turbine engine of marine purpose contains water particles, inertial separator for separating the air and water particles are provided. Saw type and wave type separator are now used to separate inflow water particle from the gas. In this paper, the design parameters of saw type separator are studied by numerical simulations. Using the commercial CFD program, Star-CCM+, Lagrangian-Eulerian method was used to perform the analysis of two phase flow of the mist in the air. This method solves Reynolds-Averaged Navier-Stokes equations in Eulerian framework for the continuous phase, while solves equation of motion for individual particles in Lagrangian framework. Lagrangian multiphase method was applied to monitor the particles of different sizes and shapes and to verify collision between particles by chasing particles. Water particles were injected through injectors located at the inlet of the separator and escape mode was used which assumes that the particles attached on the surface of inertial separator were removed from the simulation, effectively escaping the solution domain. Through the numerical computations with the inlet condition of constant water particle size in the wetness fraction of 85%, efficiency of eliminating the water particle and the pressure drop between the inlet and outlet were examined.

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

This paper presents numerical analysis and design optimization of various turbine blade cooling techniques with three-dimensional Reynolds-averaged Navier-Stokes(RANS) analysis. The fluid flow and heat transfer have been performed using ANSYS-CFX 11.0. A fan-shaped hole for film-cooling has been carried out to improve film-cooling effectiveness with the radial basis neural network method. The injection angle of hole, lateral expansion angle of hole and ratio of length-to-diameter of the hole are chosen as design variables and spatially averaged film-cooling effectiveness is considered as an objective function which is to be maximized. The impingement jet cooling has been performed to investigate heat transfer characteristic with geometry variables. Distance between jet nozzle exit and impingement plate, inclination of nozzle and aspect ratio of nozzle hole are considered as geometry variables. The area averaged Nusselt number is evaluated each geometry variables. A rotating rectangular channel with staggered array pin-fins has been investigated to increase heat transfer performance ad to decrease friction loss using KRG modeling. Two non-dimensional variables, the ratio of the eight diameter of the pin-fins and ratio of the spacing between the pin-fins to diameter of the pin-fins selected as design variables. A rotating rectangular channel with staggered dimples on opposite walls are formulated numerically to enhance heat transfer performance. The ratio of the dimple depth and dimple diameter are selected as geometry variables.

• Journal of the Korea Computer Graphics Society
• /
• v.24 no.4
• /
• pp.21-27
• /
• 2018

In this paper, we propose a real-time interaction method using user's mouth wind in mobile device. In mobile and virtual reality, user interaction technology is important, but various user interface methods is still lacking. Most of the interaction technologies are hand touch screen touch or motion recognition. In this study, we propose an interface technology that can interact with real time using user's mouth wind. The direction of the wind is determined by using the angle and the position between the user and the mobile device, and the size of the wind is calculated by using the magnitude of user's mouth wind. To show the superiority of the proposed technique, we show the result of visualizing the flow of the vector field in real time by integrating the mouth-wind interface into the Navier-Stokes equations. We show the results of the paper on mobile devices, but can be applied in the Agumented reality(AR) and Virtual reality(VR) fields requiring interface technology.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.4
• /
• pp.244-252
• /
• 2018

This paper provides quantification of the effects of the turbulence model and grid refinement on the analysis of tip vortex flows by using the RANS(Reynolds averaged Navier-Stokes) method. Numerical simulations of the tip vortex flows of the NACA $66_2$-415 elliptic hydrofoil were conducted, and two turbulence models for RANS closure were tested, i.e., the Realizable $k-{\varepsilon}$ model and the Reynolds stress transport model. Numerical results were compared with available experimental data, and it was shown that the data for the Reynolds stress transport model that were computed on the finest grid system had better agreement in reproducing the development and propagation of the tip vortex. The Realizable $k-{\varepsilon}$ model overestimated the turbulence level in the vortex core and showed a diffusive behavior of the tip vortex. The tip vortex cavitation on the hydrofoil and its trajectory also showed good agreement between the current numerical results that were obtained using the Reynolds stress transport model and the results observed in the experiment.