• Fire Science and Engineering
• /
• v.32 no.4
• /
• pp.7-16
• /
• 2018

In this study, the prediction performance of design fire curves (DF) was evaluated for gas fires in a compartment by using CFAST. The CFAST simulations adopted the 2-stage DF suggested by the previous study and the Quadratic and Exponential DF suggested by Ingason. It was found by comparing the simulation and experimental results that the overall prediction performance of the design fire cures for the spatially-averaged temperature and concentrations of $O_2$ and $CO_2$ was, from the most reasonable to the most inaccurate, 2-stage DF > Quadratic DF > Exponential DF. The CFAST simulation could not predict for the difference in the spatially-averaged temperature and concentrations of $O_2$ and $CO_2$ at door and inner side locations in a compartment. The CFAST simulations also showed a limitation in the prediction of the spatially-averaged temperature at lower layer and the concentration of CO.

• Wind and Structures
• /
• v.18 no.5
• /
• pp.549-566
• /
• 2014

In spite of progress in the numerical simulation of typhoon wind field in atmospheric boundary layer (ABL), using typhoon wind field model in conjunction with Monte Carlo simulation method can only accurately evaluate typhoon wind field over a general terrain. This method is not enough for a reliable evaluation of typhoon wind field over the actual complex terrain with surface roughness and topography variations. To predict typhoon wind field over the actual complex terrain in ABL, a hybrid numerical simulation method combined typhoon simulation used the typhoon wind field model proposed by Meng et al. (1995) and CFD simulation in which the Reynolds averaged Navier-Stokes (RANS) equations and k-${\varepsilon}$ turbulence model are used. Typhoon wind filed during typhoon Dujuan and Imbudo are simulated using the hybrid numerical simulation method, and compared with the results predicted by the typhoon wind field model and the wind field measurement data collected by Fugro Geotechnical Services (FGS) in Hong Kong at the bridge site from the field monitoring system of wind turbulence parameters (FMS-WTP) to validate the feasibility and accuracy of the hybrid numerical simulation method. The comparison demonstrates that the hybrid numerical simulation method gives more accurate prediction to typhoon wind speed and direction, because the effect of topography is taken into account in the hybrid numerical simulation method.

• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.4
• /
• pp.361-372
• /
• 2019

The flow pattern of air layers and skin-friction drag reduction by air injection are investigated to find the suitable multiphase flow model using unstructured finite-volume CFD solver for the Reynolds-averaged Navier-Stokes equations. In the present computations, two different multiphase flow modeling approaches, such as the Volume of Fluid (VOF) and the Eulerian Multi-Phase (EMP), are adopted to investigate their performances in resolving the two-phase flow pattern and in estimating the frictional drag reduction. First of all, the formation pattern of air layers generated by air injection through a circular opening on the bottom of a flat plate are investigated. These results are then compared with those of MMkiharju's experimental results. Subsequently, the quantitative ratios of skin-friction drag reduction including the behavior of air layers, within turbulent boundary layers in large scale and at high Reynolds number conditions, are investigated under the same conditions as the model test that has been conducted in the US Navy's William B. Morgan Large Cavitation Channel (LCC). From these results, it is found that both VOF and EMP models have similar capability and accuracy in capturing the topology of ventilated air cavities so called'air pockets and branches'. However, EMP model is more favorable in predicting quantitatively the percentage of frictional drag reduction by air injection.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.6 no.2
• /
• pp.12-20
• /
• 1998

Exter ballistics of a typical high-speed projectile is studied through a flow-visualization experiment and an unstructured grid Navier-Srokes computation. Experiment produced a schlieren photograph that adequately shows the characteristic features of this complex flow, namely two kinds of oblique cone shocks and turbulent wake developing into the downstream. A hybrid scheme of finite volume-element method is used to simulate the compressible Reynolds-Averaged Navier-Stok- es solution on unstructured grids. Osher's approximate Riemann solver is used to discretize the cinvection term. Higher-order spatial accuracy is obtained by MUSCL extension and van Albada ty- pe flux limiter is used to stabilize the numerical oscillation near the solution discontinuity. Accurate Gakerkin method is used to discretize the viscous term. Explict fourth-order Runge-Kutta method is used for the time-stepping, which simplifies the application of MUSCL extension. A two-layer k-$\varepsilon$ turbulence model is used to simulate the turbulent wakes accurately. Axisymmetric folw and two-dimensional flow with an angle of attack have been computed. Grid-dependency is also checked by carrying out the computation with doubled meshes. 2-D calculation shows that effect of angle of attack on the flow field is negligible. Axi-symmetric results of the computation agrees well with the flow visualization. Primary oblique shock is represented within 2-3 meshes in numerical results, and the varicose mode of the vortex shedding is clearly captured in the turbulent wake region.

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.7-13
• /
• 2010

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 5.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}10^6$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental datum in the low aspect ratio cavity (L/D = ~4.5). In the high aspect ratio cavity, however, there are other low dominant frequencies of the leading edge shear layer with the dominant frequencies of the feedback mechanism.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.13-18
• /
• 2009

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought about by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 2.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}106$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental data in the low aspect ratio cavity (L/D = ~ 4.5). In the large aspect ratio cavity, however, there are other low dominant frequencies due to the leading edge shear layer with the dominant frequencies of the feedback mechanism. The characteristics of the acoustic wave propagation are analyzed using the Correlation of Pressure Distribution (CPD).

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.5
• /
• pp.323-334
• /
• 2020

Since boundary layer transition has a significant impact on the aero-thermodynamic performance of hypersonic flight vehicles, capability of accurate prediction of transition location is essential for design and performance analysis. In this study, γ-Re_θt model is improved to predict transition of hypersonic boundary layers and validated. A coefficient in the production term of the intermittency transport equation that affects the transition onset location is constructed and applied as a function of Mach number, wall temperature, and freestream stagnation temperature based on the similarity numerical solution of compressible boundary layer. To take into account a Mach number dependency of transition onset momentum thickness Reynolds number and transition length, additional correlation equations are determined as function of Mach number and applied to Re_θc and F_length correlations of the baseline model. The suggested model is implemented to a commercial CFD code in consideration of practical use. Analysis of hypersonic flat plate and circular cone boundary layers is carried out by using the model for validation purpose. An improvement of prediction capability with respect to variation of Mach number and unit Reynolds number is identified from the comparison with experimental data.

• 한국해양학회지
• /
• v.29 no.2
• /
• pp.119-131
• /
• 1994

Vertical structures of wind-driven and tidal currents in a rectangular shaped uniform-depth basin of homogeneous water have been investigated using a mode-splitted, multi-level grid-box, hydrodynamic numerical model. The model was verified using analytical solutions for various vertical eddy viscosity profiles such as: a constant eddy viscosity, a linearly decreasing or increasing variation with depth, a quadratic variation with depth and an exponential variation with depth. Particular attention has been paid on the effects of "near-surface wall layer" on vertical shear of velocity. In numerical calculations, the whole water depth was divided into 13 levels with an unequal grid spacing. the model satisfactorily reproduces the velocity profile, but in case the eddy viscosity decreases rapidly with depth as in quadratical or exponential variation with depth, the vertical gradient of velocity near the bottom became very steep, and analytical solutions and numerical results showed some discrepancy. The vertical structures of horizontal velocity vary with both the depth-averaged value of eddy viscosity and its profiles. the velocity near the sea surface and near the bottom responded sensitively to the eddy viscosity of wall layer. For wind-driven current, the strong velocity shear was generated near the sea surface as eddy viscosity near the surface became small. For tidal current, the velocity above the sea bottom layer was almost constant regardless of the profiles of vertical eddy viscosity, but velocity in the sea bottom layer showed strong shear as eddy viscosity became small.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.297-306
• /
• 2014

The paper presents the numerical analysis of wave run-up over rubble-mound breakwaters covered by antifer units using a technique integrating Computer-Aided Design (CAD) and Computational Fluid Dynamics (CFD) software. Direct application of Navier-Stokes equations within armour blocks, is used to provide a more reliable approach to simulate wave run-up over breakwaters. A well-tested Reynolds-averaged Navier-Stokes (RANS) Volume of Fluid (VOF) code (Flow-3D) was adopted for CFD computations. The computed results were compared with experimental data to check the validity of the model. Numerical results showed that the direct three dimensional (3D) simulation method can deliver accurate results for wave run-up over rubble mound breakwaters. The results showed that the placement pattern of antifer units had a great impact on values of wave run-up so that by changing the placement pattern from regular to double pyramid can reduce the wave run-up by approximately 30%. Analysis was done to investigate the influences of surface roughness, energy dissipation in the pores of the armour layer and reduced wave run-up due to inflow into the armour and stone layer.

• Journal of computational fluids engineering
• /
• v.10 no.4 s.31
• /
• pp.51-58
• /
• 2005

The unsteady supersonic flow over two- and three-dimensional cavities has been analyzed by the integration of unsteady Reynolds-Averaged Navier-Stokes(RANS) with the k-$\omega$ turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in the cavity. An explicit 4th order Runge-Kutta scheme and an upwind TVD scheme based on the flux vector split with the van Leer limiters are used for time and space discritizations, respectively. The cavity has a L/D ratio of 3 for two-dimensional case, and same L/D and W/D ratio of I for three-dimensional case. The Mach and Reynolds numbers are 1.5 and 450000 respectively. In the three-dimensional flow, the field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follows Rossiter's formula. In the two-dimensional simulation, the self-sustained oscillating flow has more violent fluctuation inside the cavity. The primary fluctuating frequencies of two- and three- dimensional flow agree very well with the 2nd mode of Rossiter's frequency. In the three-dimensional flow, the 1st mode of frequency could be seen.