• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.16 no.10 s.115
• /
• pp.1082-1089
• /
• 2006

In this study, advanced computational analysis system has been developed in order to investigate flow-induced vibration(FIV) phenomenon for general stator-rotor cascade configurations. Relative movement of the rotor with respect to stator is reflected by modeling Independent two computational domains. Fluid domains are modeled using the unstructured grid system with dynamic moving and local deforming methods. Unsteady, Reynolds-averaged Wavier-stokes equations with one equation Spalart-Allmaras and two-equation SST ${\kappa}-{\varepsilon}$ turbulence models are solved for unsteady flow problems and also relative moving and vibration effects of the rotor cascade are fully considered. A coupled implicit time marching scheme based on the Newmark integration method is used for computing the governing equations of fluid-structure interaction problems. Detailed vibration responses for different flow conditions are presented and then vibration characteristics are physically investigated in the time domain as computational virtual tests.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.4
• /
• pp.53-62
• /
• 2006

Two-dimensional steady flowfields generated by transverse injection jets into a supersonic mainstream are numerically simulated. Fine-scale turbulence effects are represented by a k-${\omega}$ SST two-equation closure model which includes $y^+$ effects on the turbulence model. Solution convergence is evaluated by using Grid Convergence Index(GCI), a measure of uncertainty of the grid convergence. Comparison is made with experimental data and other turbulence models in term of surface static pressure distributions, the length of the upstream separation region, and the penetration height. Results indicate that the k-${\omega}$ SST model correctly predicts the mean surface pressure distribution and the upstream separation length for low static pressure ratios. However, the numerical predictions become less consistent with experimental results as the static pressure ratio increases. All these results are taken within 1% error band of grid convergence.

• Journal of Institute of Convergence Technology
• /
• v.1 no.1
• /
• pp.5-8
• /
• 2011

Thermal analysis of a cooling module using Peltier elements are performed using a commercial software, CFD-ACE+. A standard k-e two-equation turbulent model is applied in order to represent the turbulent shear stress. Computed values are compared with the theoretical values for the validation. The effect of mass flow rates and transferred heat amounts on the temperature distributions inside the cooling system is analyzed. It was found that the increase in the mass flow rates causes the exit temperature rise. The increase in the absorbed heat amount diminished the overall temperature on the fin surfaces. In the present analysis, the material characteristics of the Peltier element itself are not considered. In the future, the effect of the turbulence models and material characteristics will be studied in detail.

• Nuclear Engineering and Technology
• /
• v.44 no.8
• /
• pp.871-888
• /
• 2012

This paper aims at presenting the state of the art, the recent progress, and the perspective for the future, in the modelling of two-phase flow in the horizontal legs of a PWR. All phenomena relevant for safety analysis are listed first. The selection of the modelling approach for system codes is then discussed, including the number of fluids or fields, the space and time resolution, and the use of flow regime maps. The classical two-fluid six-equation one-pressure model as it is implemented in the CATHARE code is then presented and its properties are described. It is shown that the axial effects of gravity forces may be correctly taken into account even in the case of change of the cross section area or of the pipe orientation. It is also shown that it can predict both fluvial and torrential flow with a possible hydraulic jump. Since phase stratification plays a dominant role, the Kelvin-Helmholtz instability and the stability of bubbly flow regime are discussed. A transition criterion based on a stability analysis of shallow water waves may be used to predict the Kelvin-Helmholtz instability. Recent experimental data obtained in the METERO test facility are analysed to model the transition from a bubbly to stratified flow regime. Finally, perspectives for further improvement of the modelling are drawn including dynamic modelling of turbulence and interfacial area and multi-field models.

• International Journal of Air-Conditioning and Refrigeration
• /
• v.10 no.3
• /
• pp.117-128
• /
• 2002

This paper describes the smoke filling process of a fire field model based on a self-deve-loped SMEP (Smoke Movement Estimating Program) code to the simulation of fire induced flows in the two types of atrium space containing a ceiling heat flux. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k-$\varepsilon$ turbulence model with buoyancy production term. Also it solves the radiation equation using the discrete ordinates method. Compressibility is assumed and the perfect gas law is used. Comparison of the calculated upper-layer average tempera-ture and smoke layer clear height with the zone models has shown reasonable agreement. The zone models used are the CFAST and the NBTC one-room. For atrium fires with ceiling glass the ceiling heat flux by solar heat causes a high smoke temperature near the ceiling. However, it has no effect on the smoke movement such as the smoke layer clear heights that are important in fire safety. In conclusion, the smoke layer clear heights that are important in evacuation activity except the early of a fire were not as sensitive as the smoke layer tem-perature to the nature of ceiling heat flux condition. Thus, a fire sensor in atrium with ceiling glass has to consider these phenomena.

• Journal of Korea Water Resources Association
• /
• v.41 no.2
• /
• pp.197-212
• /
• 2008

Recently, stream flow analysis has been accomplished by one or two dimensional equations and was applied by simple momentum equations and fixed energy conservations which contain many condition limits. In this study, FLOW-3D using CFD (Computational Fluid Dynamics) was applied to stream flow analysis which can solve three dimensional RANS (Reynolds Averaged Navier-Stokes Equation) control equation to find out physical behaviors and the effect of hydraulic structures. Numerical simulation accomplished those results was compared by using turbulence models such as ${\kappa}-{\varepsilon}$, RNG (Renormalized Group) ${\kappa}-{\varepsilon}$ and LES (Large Eddy Simulation). Numerical analysis results have been illustrated by the turbulence energy effects, velocity of flow, water level pressure and eddy flows around the piers and transverse overflow type structures. These results will be able to used by basis data that catch hold of effects on long-term bed elevation changes, sediment accumulations, scours and water aggravations by removal of obsolete transverse over flow type structures in urban stream.

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

In this paper, a new computational code was developed using Chorin's artificial compressibility method to solve the two-dimensional incompressible Navier-Stokes equations. In spatial derivatives, Roe's flux difference splitting was used for the inviscid flux, while central differencing was used for the viscous flux. Furthermore, AF-ADI with dual time stepping method was implemented for accurate unsteady computations. Two-equation turbulence models, Menter's $k-{\omega}$ SST model and Coakley's $q-{\omega}$ model, hae been adopted to solve high-Reynolds number flows. A number of numerical simulations were carried out for steady laminar and turbulent flow problems as well as unsteady flow problem. The code was verified and validated by comparing the results with other computational results and experimental results. The results of numerical simulations showed that the present developed code with the artificial compressibility method can be applied to slve steady and unsteady incompressible flows.

• Journal of the Korean Solar Energy Society
• /
• v.27 no.4
• /
• pp.1-9
• /
• 2007

Sensitivity analysis of wind resource micrositing has been performed through the application case at the Antarctic King Sejong station with the most representative micrositing softwares: WAsP, WindSim and Meteodyn WT. The wind data obtained from two met-masts separated 625m were applied as a climatology input condition of micro-scale wind mapping. A tower shading effect on the met-mast installed 20m apart from the warehouse has been assessed by the CFD software Fluent and confirmed a negligible influence on wind speed measurement. Theoretically, micro-scale wind maps generated by the two met-data located within the same wind system and strongly correlated meteor-statistically should be identical if nothing influenced on wind prediction but orography. They, however, show discrepancies due to nonlinear effects induced by surrounding complex terrain. From the comparison of sensitivity analysis, Meteodyn WT employing 1-equation turbulence model showed 68% higher RMSE error of wind speed prediction than that of WindSim using the ${\kappa}-{\epsilon}$ turbulence model, while a linear-theoretical model WAsP showed 21% higher error. Consequently, the CFD model WindSim would predict wind field over complex terrain more reliable and less sensitive to climatology input data than other micrositing models. The auto-validation method proposed in this paper and the evaluation result of the micrositing softwares would be anticipated a good reference of wind resource assessments in complex terrain.

• International Journal of Aerospace System Engineering
• /
• v.2 no.2
• /
• pp.67-72
• /
• 2015

The main concern of this study is to integrate the EARSM into an industrial RANS solver in conjunction with the $k-{\omega}$ model, as proposed by Hellsten (EARSMKO2005). In order to improve the robustness, particular limiters are introduced to turbulent conservative variables, and a suitable full-approximation storage (FAS) multi-grid (MG) strategy is designed to incorporate turbulence model equations. The present limiters and MG strategy improve both robustness and efficiency significantly but without degenerating accuracy. Two discretization approachs for velocity gradient on cell interfaces are implemented and compared with each other. Numerical results of a three-dimensional supersonic square duct flow show that the proper discretization of velocity gradient improves the accuracy essentially. To assess the capability of the resulting EARSM $k-{\omega}$ model to predict complex engineering flow, the case of Common Research Model (CRM, Wing-Body) is performed. All the numerical results demonstrate that the resulting model performs well and is comparable to the standard two-equation models such as SST $k-{\omega}$ model in terms of computational effort, thus it is suitable for industrial applications.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.3
• /
• pp.444-453
• /
• 1987

In the present study, a new computational closure model is proposed in order to contain physical models in the k- and .epsilon.- equations. The time scale of the third-order diffusive transport of turbulent kinetic energy in a curved streamline flow field is assumed as a function of a velocity time scale and a curvature time scale, the latter being derived from the analogy between buoyancy and streamline curvature effects on turbulence. The curvature time scale is represented by a combination of Brunt-Vaisala frequency of the curvature instability and the velocity time scale. Besides the modification of diffusive transport time scale, the destruction term in the dissipation rate equation is modeled to incorporate the streamline curvature effect on the dissipation rate of turbulent kinetic energy as a function of the ratio between velocity time scale and curvature time scale. The new curvature dependent 2-equation model is found to yield very good prediction accuracy for the various turbulent recirculating flows. Particurarly, the recovery of the mean velocity profile in the redeveloping region after the reattachment is correctly simulated by the present model.