• Wind and Structures
• /
• 제14권5호
• /
• pp.435-447
• /
• 2011

A Computational Fluid Dynamics model is presented in this study for the simulation of the complex fluid flows with free surfaces inside the Tuned Liquid Column Dampers in horizontal motion. The characteristics of the fluid model of the TLCD in horizontal motion include the free surface of the multiphase flow and the horizontal moving frame. In this study, the time depend unsteady Standard ${\kappa}-{\varepsilon}$ turbulent model based on Navier-Stokes equations is chosen. The volume of fluid (VOF) method and sliding mesh technique are adopted to track the free surface of water inside the vertical columns of TLCD and treat the moving boundary of the walls of TLCD in horizontal motion. Several model solution parameters comprising different time steps, mesh sizes, convergence criteria and discretization schemes are examined to establish model parametric independency results. The simulation results are compared with the experimental data in the dimensionless amplitude of the water column in four different configured groups of TLCDs with four different orifice areas. The predicted natural frequencies and the head loss coefficient of TLCDs from CFD model are also compared with the experimental data. The predicted numerical results agree well with the available experimental data.

• 산업기술연구
• /
• 제22권A호
• /
• pp.59-65
• /
• 2002

The objective of this study is to provide the quantitative and qualitative computational data about the aerodynamic performance of Gurney flap on NACA 00XX airfoils and to show the optimum Gurney flap height for each airfoil. The test was performed on 7 different airfoils from NACA 0006 to NACA0024, which have a 3% chord(=c) thickness interval. For every NACA 00XX airfoil, Gurney flap heights were changed by 0.5% or 0.25% chord interval from 0 to 2.0%c to study their effects. The aerodynamic characteristics of clean and Gurney flap airfoil were compared, and the influences of Gurney flap on each airfoil were compared. As a CFD (Computational Fluid Dynamics) solver, FLUENT, based on Navier-Stokes code, was used to calculate the flow field around the airfoil. The fully-turbulent results were obtained using the standard $k-{\varepsilon}$ two-equation turbulence model. The test results showed that Gurney flap increased the lift coefficient much more than the drag coefficient over a certain range of the lift coefficient, so the lift-to-drag ratio, which is the important index of airfoil performance, was increased. Based on the test results, the relationship between the airfoil thickness and the optimum Gurney flap heights was suggested.

• 한국자동차공학회논문집
• /
• 제11권3호
• /
• pp.115-123
• /
• 2003

With the advent of high performance computers and more efficient numerical algorithms, computational fluid dynamics(CFD) has come out as a modem alternative for reducing the use of wind tunnels test in automotive engineering. However, in spite of the fact that many competent researchers have made all their talents in developing turbulence model over since the past dozen or more years, it has been an important impediment in using the CFD effectively to design machinery and to diagnose or to improve engineering problems in the industry since the turbulence model has been acting as the Achilles' tendon in aspect of the reliability even to this time. In this study, Reynolds-averaged Wavier-Stokes equations were solved to simulate an incompressible turbulent flow around a bus-like bluff body near ground plane. In order to investigate the effect of the discretisation schemes and turbulence model on the aerodynamic forces several turbulence models with five convective difference schemes are adopted. From the results of this study, it is clear that choice of turbulence model and discretisation scheme profoundly affects the computational outcome. The results also show that the adoption of RNG $k-\varepsilon$ turbulence model and nonlinear quadratic turbulence model with the second order accurate discretisation scheme predicts fairly well the aerodynamic coefficients.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1995년도 추계 학술대회논문집
• /
• pp.160-168
• /
• 1995

This paper describes computational efforts on the various energy and environmental problems using Patankar's SIMPLE method. The specific problems included in this study are : pollutant and flammable material dispersions in open and confined areas, aerator-induced flow in a lake for DO(dissolved oxygen) concentration, primary clarifier for water and waste water treatment, hood ventilation in workplace, cyclone and LNG combustors and Dow chlorination reactor. A control-volume based finite-difference method is employed together with the power-law scheme. The pressure-velocity coupling is resolved by the use of the revised version of SIMPLE, says SIMPLER and SIMPLEC. The Reynolds stresses are closed using the standard or the RNG $k-{\varepsilon}$ models. Turbulent reaction is modeled using two fast chemistry methods such as eddy breakup and conserved scalar models. Further, a nonequilibrium model is developed for the application of the chlorination process in the Dow reactor. Other important empirical models and physical insights appeared in this study are presented and discussed in a brief note. The computational method developed in this study is considered, in general, as a viable tool for the design and determination of the optimal condition of various engineering system of interest.

• 대한기계학회논문집
• /
• 제18권7호
• /
• pp.1851-1865
• /
• 1994

The objective of the present study is to predict the film cooling effectiveness by a row of holes at various injection ratios and injection angles. Numerical calculations have been performed to investigate the characteristics of flow and temperature distributions in a region near the down-stream of injection hole including the region of adverse pressure gradient. The elliptic turbulent 3-dimensional governing equations with variable thermal properties using the low-Reynolds number k-$\bar{varepsilon}$ model was solved by SIMPLE algorithm. The results showed that the presence of adverse pressure gradient and secondary vortex in the region near the downstream of injection hole induces large temperature gradent. The $45^{\circ}$ injection has higher averaged film cooling effectiveness than $60^{\circ}$ injection. But neverthless the $90^{\circ}$ injection has greater deviation from a flat plate than $45^{\circ}$ and $60^{\circ}$ injection, the $90^{\circ}$ injection has higher averaged film cooling effectiveness than $45^{\circ}$ and $60^{\circ}$ injection in the region near the downstream of injection hole.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 추계학술대회논문집B
• /
• pp.568-572
• /
• 2000

A numerical investigation was performed to determine the effect of airfoil thickness on the optimum Gurney flap height using NACA 00XX series airfoils. Seven airfoils which have 3% chord thickness difference were used. These were NACA 0006, 0009, 0012, 0015, 0018, 0021, and 0024. A Navier-Stokes code, FLUENT, was used to calculate the flow field about airfoil. The fully turbulent results were obtained using the standard $k-{\varepsilon}$ two-equation turbulence model. To provide a check case fur our computational method, numerical studies for NACA 4412 airfoil were made and compared with already existing experimental data for this airfoil by Wadcock. For every NACA 00XX airfoil, Gurney flap heights ranging from 0.5% to 2.0% chord were changed by 0.5% chord interval and their effects were studied. With the numerical solutions, the relationship between $(L/D)_{max}$ and airfoil thickness as a function of flap height and the relationship between $(L/D)_{max}$ and flap height as a function of airfoil thickness were investigated. The same relationship for $(C_l)_{max}$ also were shown. From these results, the optimum flap size for each airfoil thickness can be determined and vice versa.

• 한국산업융합학회 논문집
• /
• 제25권1호
• /
• pp.61-70
• /
• 2022

Fire damage time in high-rise buildings and wildland fire increasing every year. The use of high-pressure fire pumps is required to effectively extinguish fires. Reflecting the curvature effect of the fire hose occurring at the actual fire fighting site, this study provides a database of pressure drop, discharge velocity and maximum discharge height through C FD numerical analysis and it can provide using standards for fire extinguishing. Two Reynolds numbers of 200000 and 400000 were numerically analyzed at 0° -180° bending with water of 25℃ as a working fluid in hoses with a diameter of 65mm, a length of 15m, and a radius of curvature of 130mm. Realizable k-ε turbulence model was used and standard wall function was used. The pressure drop increases as the bending angle increases, and the maximum value at 90° and then decreases. The increasing rate is greater than the decrease. The velocity of the secondary flow also decreases after having the maximum value at 90°. The decreasing rate is greater than the increase. The turbulent kinetic energy increases to 120° and decreases with the maximum value. Pressure drop, velocity of the secondary flow, and turbulence kinetic energy are measured larger in the second bending region than in the first bending region.

• 한국화재소방학회논문지
• /
• 제11권1호
• /
• pp.21-35
• /
• 1997

The natural convection and combined heat transfer induced by fire in a rectangular enclosure is numerically studied. The model for this numerical analysis is partially opened right wall. The solution procedure includes the standard k-$\varepsilon$ model for turbulent flow and the discrete ordinates method (DOM) is used for the calculation of radiative heat transfer equation. In numerical study, SIMPLE algorithm is applied for fluid flow analysis, and the investigations of combustion gas induced by fire is performed by FAST model of HAZARD I program. In this study, numerical simulation on the combined naturnal convection and radiation is carried out in a partial enclosure filled with absorbed-emitted gray media, but is not considered scattering problem. The streamlines, isothermal lines, average radiation intensity and kinetic energy are compared the results of pure convection with those of the combined convection-radiation, the combined heat transfer. Comparing the results of pure convection with those of the combined convection-radiation, the combined heat transfer analysis shows the stronger circulation than those of the pure convection. Three different locations of heat source are considered to observe the effect of heat source location on the heat transfer phenomena. As the results, the circulation and the heat transfer in the left region from heating block are much more influenced than those in the right region. It is also founded that the radiation effect cannot be neglected in analyzing the building in fire. And as the results of combustion gas analysis from FAST model, it is found that O2 concentration is decreased according to time. While CO and CO2 concentration are rapidly increased in the beginning(about 100sec), but slowly decreased from that time on.

• Ocean Systems Engineering
• /
• 제6권1호
• /
• pp.23-60
• /
• 2016

The major objective of this study was to develop further understanding of 3D nearshore hydrodynamics under a variety of wave and tidal forcing conditions. The main tool used was a comprehensive 3D numerical model - combining the flow module of Delft3D with the WAVE solver of XBeach - of nearshore hydro- and morphodynamics that can simulate flow, sediment transport, and morphological evolution. Surf-swash zone hydrodynamics were modeled using the 3D Navier-Stokes equations, combined with various turbulence models (${\kappa}-{\varepsilon}$, ${\kappa}-L$, ATM and H-LES). Sediment transport and resulting foreshore profile changes were approximated using different sediment transport relations that consider both bed- and suspended-load transport of non-cohesive sediments. The numerical set-up was tested against field data, with good agreement found. Different numerical experiments under a range of bed characteristics and incident wave and tidal conditions were run to test the model's capability to reproduce 3D flow, wave propagation, sediment transport and morphodynamics in the nearshore at the field scale. The results were interpreted according to existing understanding of surf and swash zone processes. Our numerical experiments confirm that the angle between the crest line of the approaching wave and the shoreline defines the direction and strength of the longshore current, while the longshore current velocity varies across the nearshore zone. The model simulates the undertow, hydraulic cell and rip-current patterns generated by radiation stresses and longshore variability in wave heights. Numerical results show that a non-uniform seabed is crucial for generation of rip currents in the nearshore (when bed slope is uniform, rips are not generated). Increasing the wave height increases the peaks of eddy viscosity and TKE (turbulent kinetic energy), while increasing the tidal amplitude reduces these peaks. Wave and tide interaction has most striking effects on the foreshore profile with the formation of the intertidal bar. High values of eddy viscosity, TKE and wave set-up are spread offshore for coarser grain sizes. Beach profile steepness modifies the nearshore circulation pattern, significantly enhancing the vertical component of the flow. The local recirculation within the longshore current in the inshore region causes a transient offshore shift and strengthening of the longshore current. Overall, the analysis shows that, with reasonable hypotheses, it is possible to simulate the nearshore hydrodynamics subjected to oceanic forcing, consistent with existing understanding of this area. Part II of this work presents 3D nearshore morphodynamics induced by the tides and waves.

• 한국안전학회지
• /
• 제15권1호
• /
• pp.72-79
• /
• 2000

In recent, occupational diseases in harmful working places become a social issue. It is the well-known fact that a respiration in polluted working places exert a serious effect on health of workers. Accordingly, the cutting off contaminants air originally is the best way to improve working environments. In these cases, ventilation systems should be essentially installed to dilute or exhaust the contaminated indoor air. In this study, we investigated the characteristics of ventilation system of the noxious gas in working indoor places with natural ventilation by using COMET. The numerical simulations were carried out the natural ventilation with two phase(air, dust). For turbulent flow, Reynolds stresses were closed by the standard $\kappa$-$\varepsilon$ model. The results are as follows ; 1) In the natural exhaust in the working place, the flows of the central region have a more rapid velocity vector than the right and left one. 2) Numerical results show that the distribution of contaminants concentration have greater influence on convection than the case of diffusion by government of velocity vectors. 3) To observe the velocity variation with distance, three location of distance are considered. As results, it shows that the velocity are 0.075(m/s) at y=5(m), 10(m) and mean concentration are raised 10.6% at y=5(m), 10(m). 4) We have presented the useful data for the adequate counterplan in the harmful working places by carrying out the various investigation of the natural ventilation.