• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.24 no.3
• /
• pp.205-212
• /
• 2014

In this paper, sunroof buffeting analysis for Hyundai simple model(HSM) is studied computationally. For validation, the velocity profile of boundary layer around the opening of HSM was obtained and compared with experimental results. The analysis of sunroof buffeting is done in two parts. First a steady state solution is obtained using the Reynolds Averaged Navier Stokes (RANS) solver, and then the computed flow field information is used as input for CAA++. Second transient simulation by CAA++ is performed for the peak sound pressure levels and peak frequencies of buffeting noise over the ranges of flow velocities. The benchmark results of frequency and sound pressure levels showed the general phenomena and matched well with the experimental data obtained by Hyundai Motor Car.

• Journal of computational fluids engineering
• /
• v.14 no.4
• /
• pp.23-30
• /
• 2009

The existing code which solves two-dimensional RANS(Reynolds Averaged Navier-Stokes) equations and 2-equation turbulence model equations was modified to enable numerical simulation of various supersonic flows. For this, various boundary conditions have been implemented to the code. Bleed boundary condition was incorporated into the code for calculating wall mean flow quantities. Furthermore, boundary conditions for the turbulence quantities along rough surfaces as well as porous walls were applied to the code. The code was verified and validated by comparing the computational results against the experimental data for the supersonic flows over bleed region on a flat plate. Furthermore, numerical simulations for supersonic shock boundary layer interaction with a bleed region were performed and their results were compared with the existing computational results.

• Wind and Structures
• /
• v.9 no.1
• /
• pp.37-58
• /
• 2006

The paper describes a study about effects of upstream hills on design wind loads using two mathematical approaches: Computational Fluid Dynamics (CFD) and Artificial Neural Network (NN for short). For this purpose CFD and NN tools have been developed using an object-oriented approach and C++ programming language. The CFD tool consists of solving the Reynolds time-averaged Navier-Stokes equations and $k-{\varepsilon}$ turbulence model using body-fitted nearly-orthogonal coordinate system. Subsequently, design wind load parameters such as speed-up ratio values have been generated for a wide spectrum of two-dimensional hill geometries that includes isolated and multiple steep and shallow hills. Ground roughness effect has also been considered. Such CFD solutions, however, normally require among other things ample computational time, background knowledge and high-capacity hardware. To assist the enduser, an easier, faster and more inexpensive NN model trained with the CFD-generated data is proposed in this paper. Prior to using the CFD data for training purposes, extensive validation work has been carried out by comparing with boundary layer wind tunnel (BLWT) data. The CFD trained NN (CFD-NN) has produced speed-up ratio values for cases such as multiple hills that are not covered by wind design standards such as the Commentaries of the National Building Code of Canada (1995). The CFD-NN results compare well with BLWT data available in literature and the proposed approach requires fewer resources compared to running BLWT experiments.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.311-314
• /
• 2006

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split using van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge-Kutta method. The aspect ratio of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two-dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the fire cavity flow cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Journal of Mechanical Science and Technology
• /
• v.20 no.8
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split with van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge- Kutta method. The aspect ratios of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two- dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the first cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.6
• /
• pp.16-22
• /
• 2004

Unsteady vortex shedding behind a rectangular cylinder near a wall influences both increasing of drag and dynamic stability of heavy vehicle, bridge or building. Incompressible Averaged Navier-Stokes equation with modified ${\varepsilon}-SST$ turbulence model is adapted for investigating the flow field between the rectangular cylinder and the wall. In case the vortex shedding happens, not only the averaged maximum velocity is higher than other cases, but the position of the maximum velocity is closer to the lower surface of rectangular cylinder. On this study, it is confirmed that the vortex shedding behind a rectangular cylinder can be suppressed by the passive control method added by horizontal and vertical fences to the lower surface of rectangular cylinder.

• Nuclear Engineering and Technology
• /
• v.56 no.1
• /
• pp.9-18
• /
• 2024

In-vessel retention through external reactor vessel cooling (IVR-ERVC) is a severe accident management (SAM) strategy that has been adopted and used in many nuclear reactors such as AP1000, APR1400, and light water reactor etc. Some reactor accidents have raised concerns about nuclear reactors among residents, leading to a decrease in residents' acceptability and many studies on SAM are being conducted. Experiments on IVR-ERVC are almost impossible due to its specificity, so fluid characteristics are analyzed through BALI experiments with similar condition. In this study, computational fluid dynamics (CFD) via Reynolds-averaged Navier-Stokes (RANS) and large eddy simulation (LES) for BALI experiments were performed. Steady-state CFD analysis was performed on three turbulence models, and SST k-ω model was in good agreement with the experimental measurement temperature within the maximum error range of 1.9%. LES CFD analysis was performed based on the RANS analysis results and it was confirmed that the temperature and wall heat flux for depth was consistent within an error range of 1.0% with BALI experiment. The LES CFD analysis results were compared with those of the Lagrangian-based solver. LES matched the temperature distribution better than SOPHIA, but SOPHIA calculated the position of boundary between stratified layer and convective layer more accurately. On the other hand, Lagrangian-based solver predicted several small eddy behaviors of the convective layer and LES predicted large vortex behavior. The vibration characteristics near the cooling part of the BALI experimental device were confirmed through Fast Fourier Transform (FFT) investigation. It was found that the power spectral density for pressure at least 10 times higher near the side cooling than near the top cooling.

• Environmental Sciences Bulletin of The Korean Environmental Sciences Society
• /
• v.4 no.1
• /
• pp.31-42
• /
• 2000

The dry deposition velocities and fluxes of air pollutants such as SO2(g), O3(g), HNO3(g), sub-micron particulates, NO3(s), and SO42-(s) were estimated according to local meteorological elements in the atmospheric boundary layer. The model used for these calculations was the multiple layer resistance model developed by Hicks et al.1). The meteorological data were recorded on an hourly basis from July, 1990 to June, 1991 at the Austin Cary forest site, near Gainesville FL. Weekly integrated samples of ambient dry deposition species were collected at the site using triple-fiter packs. For the study period, the annual average dry deposition velocities at this site were estimated as 0.87$\pm$0.07 cm/s for SO2(g), 0.65$\pm$0.11 cm/s for O3(g), 1.20$\pm$0.14cm/s for HNO3(g), 0.0045$\pm$0.0006 cm/s for sub-micron particulates, and 0.089$\pm$0.014 cm/s for NO3-(s) and SO42-(s). The trends observed in the daily mean deposition velocities were largely seasonal, indicated by larger deposition velocities for the summer season and smaller deposition velocities for the winter season. The monthly and weekly averaged values for the deposition velocities did not show large differences over the year yet did show a tendency of increased deposition velocities in the summer and decreased values in the winter. The annual mean concentrations of the air pollutants obtained by the triple filter pack every 7 days were 3.63$\pm$1.92 $\mu\textrm{g}$/m3 for SO42-, 2.00$\pm$1.22 $\mu\textrm{g}$/m-3 for SO2, 1.30$\pm$0.59 $\mu\textrm{g}$/m-3 for HNO3, and 0.704$\pm$0.419 $\mu\textrm{g}$/m3 for NO3-, respectively. The air pollutant with the largest deposition flux was SO2 followed by HNO3, SO42-(S), and NO3-(S) in order of their magnitude. The sulfur dioxide and NO3- deposition fluxes were higher in the winter than in the summer, and the nitric acid and sulfate deposition fluxes were high during the spring and summer.

• Nuclear Engineering and Technology
• /
• v.53 no.5
• /
• pp.1676-1685
• /
• 2021

We develop a new high-fidelity multiphysics model to simulate boron chemistry in the porous Chalk River Unidentified Deposit (CRUD) deposits. Heat transfer, capillary flow, solute transport, and chemical reactions are fully coupled. The evaporation of coolant in the deposits is included in governing equations modified by the volume-averaged assumption of wick boiling. The axial offset anomaly (AOA) of the Seabrook nuclear power plant is simulated. The new model reasonably predicts the distributions of temperature, pressure, velocity, volumetric boiling heat density, and chemical concentrations. In the thicker CRUD regions, 60% of the total heat is removed by evaporative heat transfer, causing boron species accumulation. The new model successfully shows the quantitative effect of coolant evaporation on the local distributions of boron. The total amount of boron in the CRUD layer increases by a factor of 1.21 when an evaporation-driven increase of soluble and precipitated boron concentrations is reflected. In addition, the concentrations of B(OH)₃ and LiBO₂ are estimated according to various conditions such as different CRUD thickness and porosity. At the end of the cycle in the AOA case, the total mass of boron incorporated in CRUD deposits of a reference single fuel rod is estimated to be about 0.5 mg.

• Advances in aircraft and spacecraft science
• /
• v.9 no.5
• /
• pp.415-431
• /
• 2022

Secondary flows have a huge impact on losses generation in modern low pressure gas turbines (LPTs). At design point, the interaction of the blade profile with the end-wall boundary layer is responsible for up to 40% of total losses. Therefore, predicting accurately the end-wall flow field in a LPT is extremely important in the industrial design phase. Since the inlet boundary layer profile is one of the factors which most affects the evolution of secondary flows, the first main objective of the present work is to investigate the impact of two different inlet conditions on the end-wall flow field of the T106A, a well known LPT cascade. The first condition, labeled in the paper as C1, is represented by uniform conditions at the inlet plane and the second, C2, by a flow characterized by a defined inlet boundary layer profile. The code used for the simulations is based on the Discontinuous Galerkin (DG) formulation and solves the Reynolds-averaged Navier-Stokes (RANS) equations coupled with the Spalart Allmaras turbulence model. Secondly, this work aims at estimating the influence of viscosity and turbulence on the T106A end-wall flow field. In order to do so, RANS results are compared with those obtained from an inviscid simulation with a prescribed inlet total pressure profile, which mimics a boundary layer. A comparison between C1 and C2 results highlights an influence of secondary flows on the flow field up to a significant distance from the end-wall. In particular, the C2 end-wall flow field appears to be characterized by greater over turning and under turning angles and higher total pressure losses. Furthermore, the C2 simulated flow field shows good agreement with experimental and numerical data available in literature. The C2 and inviscid Euler computed flow fields, although globally comparable, present evident differences. The cascade passage simulated with inviscid flow is mainly dominated by a single large and homogeneous vortex structure, less stretched in the spanwise direction and closer to the end-wall than vortical structures computed by compressible flow simulation. It is reasonable, then, asserting that for the chosen test case a great part of the secondary flows details is strongly dependent on viscous phenomena and turbulence.