• 한국소음진동공학회논문집
• /
• 제24권3호
• /
• pp.205-212
• /
• 2014

이 연구에서 현대자동차의 단순실험모델(HSM)에 대한 썬루프 버페팅에 대한 수치해석이 수행되었다. 검증을 위하여 HSM 목부위의 경계층에 대한 속도분포 해석결과를 실험결과와 비교하였다. 썬루프 해석은 두 단계로 이루어졌다. 첫 번째로 난류 RANS 모델을 이용하여 정상상태 해석이 수행되었으며, 해석결과는 CAA++의 입력값으로 사용된다. 두 번째 단계는 유동속도에 대한 1차 최대 압력피크와 버페팅 주파수 해석을 위한 비정상상태 해석이 CAA++에서 이루어졌다. 주파수와 음향압력의 수치해석 결과는 타당한 물리적 현상을 보여주고 있으며, 현대 자동차의 실험결과와 잘 일치하는 결과를 보여주었다.

• 한국전산유체공학회지
• /
• 제14권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
• /
• 제9권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년도 PARALLEL CFD 2006
• /
• 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
• /
• 제20권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.

• 한국항공우주학회지
• /
• 제32권6호
• /
• pp.16-22
• /
• 2004

지면엔 근접한 사각주 후면에서 발생하는 비정상 와류 배출은 지상 운송체, 교량, 건물 등의 항력 증가뿐 아니라, 동안정성에도 큰 영향을 미친다. 비압축성 평균 Navier-Stokes 방정식에 수정된 ${\varepsilon}-SST$ 난류 모델을 적용하여 사각주 하부와 지면과의 간극 유동을 해석하였다. 사각주 후류에서 와류가 발생하는 경우에는 간극에서의 평균 최대 속도가 억제된 경우에 비하여 높으며, 또한 최대 속도의 위치 또한 사각주 하부에 근접한 것을 확인하였다. 본 연구에서는 사각주 하부에 수평, 수직의 펜스를 설치하는 수동 제어기법을 적용하여 사각 주 후류의 와류 배출용 억제할 수 있다.

• Nuclear Engineering and Technology
• /
• 제56권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
• /
• 제4권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
• /
• 제53권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
• /
• 제9권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.