• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.130-130
• /
• 2023

제트류는 다양한 크기와 운동량의 에디가 복잡하게 혼합되어 이루어져 있으며, 이를 정확하게 모델링하고 이해하기 위해서는 제트류의 다양한 특성들을 잘 반영하여 연구를 수행해야 한다. 다양한 연구 수행 방법 중 수치해석 방법은 상대적으로 공간 및 시간적 비용이 적게 들어서 널리 사용되고 있다. 이러한 수치해석 방법에는 DNS(Direct Numerical Simulation), LES(Large Eddy Simulation), RANS(Reynolds Averaged Navier Stokes) 등이 있으며, 그중 LES는 난류 모델링을 사용하는 RANS 방법에 비해 더욱 정확한 흐름 모델링을 제공하는 장점이 있다. 이러한 LES는 대규모 에디는 직접 해석하면서, 일정 크기 이하의 에디는 모델링을 사용해 해석하는 것이 특징이다. 하지만, LES를 사용하기 위해서는 적절한 그리드 크기를 결정하는 것이 중요하며, 이는 모델의 정확성과 연산 비용에 큰 영향을 미친다. 하지만, 여전히 적절한 그리드 크기를 결정하는 것은 어려운 문제이다. 이러한 LES 모델링을 사용할 때 적절한 그리드 크기를 결정하기 위해서는 정확한 시간 평균 속도 변동을 연구하는 것이 앞서 선행되어야 한다. 따라서, 본 연구에서는 기계학습 기반 접근 방식을 사용하여 난류 제트 내 시간 평균 속도 변동을 예측하는 연구를 진행하였다. 즉, 난류 제트 역학을 이해하는 데 중요한 파라미터인 시간 평균 유속을 이용하여 시간 평균 속도 변동을 예측하는 데 초점을 맞추었다. 모델의 성능은 평균 제곱 오차와 R-제곱 등 다양한 지표를 사용하여 평가되었다.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.8
• /
• pp.829-837
• /
• 2012

Direct numerical simulations were carried out for turbulent channel flows with $Re_{\tau}$ = 180, 395 and 590 to investigate the turbulent flow structure related to the Reynolds shear stress. By examining the probability density function, the second quadrant (Q2) events with the largest contribution to the mean Reynolds shear stress were identified. The change in the inclination angle of Q2 events varies with wall units in $y^+<50$ and with the channel half height in y/h > 0.5. Conditionally averaged flow fields for the Q2 event show that the flow structures associated with Reynolds shear stress are a quasi-streamwise vortex in the buffer layer and a hairpin-shaped vortex in the outer layer. Three-dimensional visualization of the distribution of high Reynolds shear stress reveals that the organization of hairpin vortices in the outer layer having a size of 1.5~3 h is associated with large-scale motions with high Reynolds shear stress in the outer layer.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.11 no.4
• /
• pp.133-140
• /
• 2008

The purpose of this study is to provide response methods to minimize the damage from chemical terrorism in a naturally ventilated indoor system using several types of dispersion simulations. Three chemical warfare agents such as sarin(GB), phosgene and chlorine gas which have high potential to be used in terror or to be involved with accidents were selected in this simulation. Fire dynamic simulation based on Large Eddy Simulation which is effective because of less computational effort and detailed expression of the dispersion flow was adopted to describe the dispersion behavior of these agents. When the vent speed is 0.005m/s, the heights of 0.1 agent mass fraction are 0.9m for sarin, 1.0m for phosgene and 1.1m for chlorine gas, and the maximum mass fraction are 0.27 for all three agents. However, when the vent speed is increased to 0.05m/s, the heights of 0.1 agent mass fraction become 1.6m for all three agents and maximum mass fraction inside the room increase to 0.70 for sarin, 0.58 for phosgene and 0.53 for chlorine gas. It is shown that molecular weight of the agents has an important role for dispersion, and it is important to install ventilation system with height less than 1.6m to minimize the damage from chemical toxicity.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.43-58
• /
• 2022

The BARC flow is studied via Direct Numerical Simulation at a relatively low turbulent Reynolds number, with focus on the geometrical representation of the leading-edge (LE) corners. The study contributes to further our understanding of the discrepancies between existing numerical and experimental BARC data. In a first part, rounded LE corners with small curvature radii are considered. Results show that a small amount of rounding does not lead to abrupt changes of the mean fields, but that the effects increase with the curvature radius. The shear layer separates from the rounded LE at a lower angle, which reduces the size of the main recirculating region over the cylinder side. In contrast, the longitudinal size of the recirculating region behind the trailing edge (TE) increases, as the TE shear layer is accelerated. The effect of the curvature radii on the turbulent kinetic energy and on its production, dissipation and transport are addressed. The present results should be contrasted with the recent work of Rocchio et al. (2020), who found via implicit Large-Eddy Simulations at larger Reynolds numbers that even a small curvature radius leads to significant changes of the mean flow. In a second part, the LE corners are fully sharp and the exact analytical solution of the Stokes problem in the neighbourhood of the corners is used to locally restore the solution accuracy degraded by the singularity. Changes in the mean flow reveal that the analytical correction leads to streamlines that better follow the corners. The flow separates from the LE with a lower angle, resulting in a slightly smaller recirculating region. The corner-correction approach is valuable in general, and is expected to help developing high-quality numerical simulations at the high Reynolds numbers typical of the experiments with reasonable meshing requirements.

• Nuclear Engineering and Technology
• /
• v.56 no.4
• /
• pp.1310-1319
• /
• 2024

The reactor cavity cooling system (RCCS) is a passive reactor safety system commonly present in the designs of High-Temperature Gas-cooled Reactors (HTGR) that removes heat from the reactor pressure vessel by means of natural convection and radiation. It is one of the factors responsible for ensuring that the reactor does not melt down under any plausible accident scenario. For the simulation of accident scenarios, which are transient phenomena unfolding over a span of up to several days, intermediate fidelity methods and system codes must be employed to limit the models' execution time. These models can quantify radiation heat transfer well, but heat transfer caused by natural convection must be quantified with the use of correlations for the heat transfer coefficient. It is difficult to obtain reliable correlations for HTGR RCCS heat transfer coefficients experimentally due to such a system's size. They could, however, be obtained from high-fidelity steady-state simulations of RCCSs. The Rayleigh number in RCCSs is too high for using a Direct Numerical Simulation (DNS) technique; thus, a Reynolds-Averaged Navier-Stokes (RANS) approach must be employed. There are many RANS models, each performing best under different geometry and fluid flow conditions. To find the most suitable one for simulating an RCCS, the RANS models need to be validated. This work benchmarks various RANS models against three experiments performed on the HTTR RCCS Mockup by the Japanese Atomic Energy Agency (JAEA) in 1993. This facility is a 1/6 scale model of a vessel cooling system (VCS) for the High Temperature Engineering Test Reactor (HTTR), which is operated by JAEA. Multiple RANS models were evaluated on a simplified 2d-axisymmetric geometry. They were found to reproduce the experimental temperature profiles with errors of up to 22% for the lowest temperature benchmark and 15% for the higher temperature benchmarks. The results highlight that the pragmatic turbulence models need to be validated for high Rayleigh natural convection-driven flows and improved accordingly, more publicly available experimental data of RCCS resembling experiments is needed and indicate that a 2d-axisymmetric geometry approximation is likely insufficient to capture all the relevant phenomena in RCCS simulations.

• Journal of the Korean Society of Combustion
• /
• v.18 no.1
• /
• pp.13-20
• /
• 2013

Leading edge statistics are obtained by direct numerical simulation(DNS) of freely propagating incompressible and stagnating compressible turbulent premixed flames. Conditional averages of velocities in terms of reaction progress variable, c, and local flame surface density, ${\sum}^{\prime}_f$, are defined and compared through the flame brush. It holds asymptotically that $<u>_f=<S_d>_f$ and $<u>_u-<u>_b=D_t/L_w$ with the characteristic length scale of $\bar{c}$ variation, $L_w$. It also holds that $<u>_b=<u>_f$ for a freely propagating flame under no mean strain rate. The turbulent burning velocity, $S_T$, is determined by the conditional statistics at the leading edge under large activation energy.

• Journal of Mechanical Science and Technology
• /
• v.19 no.9
• /
• pp.1801-1810
• /
• 2005

The effect of longitudinal thin rectangular riblets aligned with the flow direction on turbulent channel flow has been investigated using direct numerical simulation. The thin riblets have been modeled using the immersed boundary method (IBM) where the velocities at only one set of vertical nodes at the riblets positions are enforced to be zeros. Different spacings, ranging between 11 and 43 wall units, have been simulated aiming at getting the optimum spacing corresponding to the maximum drag reduction while keeping the height/spacing ratio at 0.5. Reynolds number based on the friction velocity ${\mu}_\tau$ and the channel half depth $\delta$ is set to 150. The flow is driven by adjusted pressure gradient so that the mass flow rate is kept constant in all the simulations. This study shows similar trend of the drag ratio to that of the experiments at the different spacings. Also, this research provides an optimum spacing of around 17 wall units leading to maximum drag reduction as experimental data. Explanation of drag increasing/decreasing mechanism is highlighted.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2004.04a
• /
• pp.67-70
• /
• 2004

An investigation was performed to study the impact damage of the laminated composite plates caused by a low- velocity foreign object with multi-scale modeling based on the concepts of Direct Numerical Simulation (DNS)[4]. In the micro-scale part, we discretize the composite plates through separate modeling of fiber and matrix for the local microscopic analysis. A micro-scalemodel was developed for predicting the initiation of the damage and the extent of the final damage as a function of material properties, laminate configuration and the impactor's mass, etc. Anda macro-scale model was developed for description of global dynamic behavior. The connection betweenmicroscopic and macroscopic is implemented by the tied interface constraints of LS-DYNA contact card. A transient dynamic finite element analysis was adopted for calculating the contact force history and the stresses and strains inside the composites during impact resulting from a point-nose impactor. The low-velocity impact events such as contact force, deformation, etc. are simulated in the macroscopic sense and the impact damages, fiber-breakage, matrix cracking and delamination etc. are examined in the microscopic sense.

• Wind and Structures
• /
• v.30 no.5
• /
• pp.465-474
• /
• 2020

Vortex-induced vibrations of a yawed flexible cylinder near a plane boundary are numerically investigated at a Reynolds number Re_n= 500 based on normal component of freestream velocity. Free to oscillate in the in-line and cross-flow directions, the cylinder with an aspect ratio of 25 is pinned-pinned at both ends at a fixed wall-cylinder gap ratio G/D = 0.8, where D is the cylinder diameter. The cylinder yaw angle (α) is varied from 0° to 60° with an increment of 15°. The main focus is given on the influence of α on structural vibrations, flow patterns, hydrodynamic forces, and IP (Independence Principle) validity. The vortex shedding pattern, contingent on α, is parallel at α=0°, negatively-yawed at α ≤ 15° and positively-yawed at α ≥ 30°. In the negatively- and positively-yawed vortex shedding patterns, the inclination direction of the spanwise vortex rows is in the opposite and same directions of α, respectively. Both in-line and cross-flow vibration amplitudes are symmetric to the midspan, regardless of α. The RMS lift coefficient C_L,rms exhibits asymmetry along the span when α ≠ 0°, maximum C_L,rms occurring on the lower and upper halves of the cylinder for negatively- and positively-yawed vortex shedding patterns, respectively. The IP is well followed in predicting the vibration amplitudes and drag forces for α ≤ 45° while invalid in predicting lift forces for α ≥ 30°. The vortex-shedding frequency and the vibration frequency are well predicted for α = 0° - 60° examined.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.4
• /
• pp.294-306
• /
• 2005

Most immersed tunnels investigated have been investigated based on the engineer's experience with design and construction. From engineering point of view, it is very important to understand the wave interaction with the seabed and immersed tunnel, since the stability of an immersed tunnel depends largely on the behavior of the seabed foundation. In this study, for the first stage research to find out the mechanism of the wave interaction with the seabed and immersed tunnel, the benchmarking method called as direct numerical simulation (DNS) was employed to analyze comprehensively the wave-induced pore water pressures, vorticity and flows in seabed or inside rubble stone around the immersed tunnel. The immersed tunnel is modeled based on Busan-Geoje fixed link project in Korea, which is now on the stage of planning. Moreover, the nonlinear water wave interaction with an immersed tunnel/its seabed foundation was thoroughly examined with regard to the stabilities of the immersed tunnel subjected to various water wave conditions, median grain size and so forth.