• Proceedings of the Korea Institute of Fire Science and Engineering Conference
• /
• 2008.04a
• /
• pp.234-237
• /
• 2008

A direct numerical simulation (DNS) code suitable for the prediction of buoyant jet diffusion flames was developed in this study. The thermodynamic and transport properties were evaluated using CHEMKIN package to enhance the prediction performance of the developed DNS code. A two dimensional simulations were performed for the jet diffusion flames in normal and zero-gravity conditions where the Froude numbers are 5 and infinity, respectively. The simulated buoyant jet diffusion flame in normal gravity showed that the unsteady and dynamic motion although the reynolds number is low (400). It was identified that the flame in normal gravity flickered periodically. The periodic motion of the flame disappeared in zero-gravity condition. The dynamic motion of the buoyant jet diffusion flame could be well understood by comparing the flame structures obtained by the simulations of normal and zero-gravity conditions.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.23 no.2
• /
• pp.60-67
• /
• 2019

To understand the low frequency instability(LFI) characteristics in hybrid rockets combustion, effects of equivalence ratio variations on the phase shift between pressure and heat release oscillations were investigated by using the direct numerical simulation. The change in the equivalence ratio of the main chamber was simulated by the temperature and composition variation of the combustion gas introduced into the post-combustion chamber. In the results, additional combustion appeared along with vortex generation at the backward step, and combustion pressure and heat release oscillations were observed as the vortex moved. In addition, the results confirmed that the phase difference between the pressure and heat release oscillation shifts because of the changes in the propagation velocity of pressure wave as the temperature of combustion gas changes.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.18 no.4
• /
• pp.383-389
• /
• 2006

In three dimensional wave field, a direct numerical simulation model, which is able to handle free-fall and slope-fall of a waterdrop due to overtopping, is proposed to evaluate the overtopping rate on the coastal structures like an absorbing revetment. A comparison between the numerical model and existing experimental results for overtopping rate was made to validate the proposed numerical model's accuracy, and showed fairly good agreement between them. It is confirmed in numerical and hydraulic(existing) model test that the overtopping quantity on a absorbing revetment becomes larger with an increase in Ursell number. Also, the overtopping rate estimated by 3-D numerical model is compared with it obtained by 2-D numerical model.

• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.12
• /
• pp.977-982
• /
• 2009

Direct numerical simulation (DNS) of a turbulent boundary layer with moderate Reynolds number was performed to scrutinize streamwise-coherence of hairpin packet motions. The Reynolds number based on the momentum thickness (${\theta}_{in}$) and free-stream velocity (U${\infty}$) was varied in the range $Re_{\theta}$=1410${\sim}$2540 which was higher than the previous numerical simulations in the turbulent boundary layer. In order to include the groups of hairpin packets existing in the outer layer, large computational domain was used (more than 50${\delta}_o$, where ${\theta}_o$ is the boundary layer thickness at the inlet in the streamwise domain). Characteristics of packet motions were investigated by using instantaneous flow fields, two-point correlation and conditional average flow fields in xy-plane. The present results showed that a train of hairpin packet motions was propagating coherently along the downstream and these structures induced the very large-scale motions in the turbulent boundary layer.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.6
• /
• pp.860-869
• /
• 2000

Numerical and experimental investigations are performed for the molecular transition and slip flows in pumping channels of a disk-type drag pump. The flow occurring in the pumping channel develops from the molecular transition to the slip flow traveling downstream. Two different numerical methods are used in this analysis: the first one is a continuum approach in solving the Navier-Stokes equations with slip boundary conditions, and the second one is a stochastic approach through the use of the direct simulation Monte Carlo method. In the experimental study, the inlet pressures are measured for various outlet pressures in the range of 0.1{\sim}4Torr. From the present study, the numerical results of predicting the performance, obtained by both methods, agree well with the experimental data for the range of Knudsen number $Kn{\leq}0.1$ (i.e., the slip flow regime). But the results from the second method only agree with the experimental data for Kn>0.1(i.e., the molecular transition regime)

• Proceedings of the KSR Conference
• /
• 2011.10a
• /
• pp.1594-1600
• /
• 2011

Recently developed FDS5 CFD code has employed a near-wall flow treatment method which is Werner-Wengle wall law provided by NIST(National Institute of Standards and Technology). In this study, the wall law has been verified against DNS(Direct Numerical Simulation) data in the parallel plate. The $y^+$ was kept above 11 to fulfill the near-wall flow requirement in the grid generation. The total grid was $32{\times}32{\times}32$. The boundary condition for inlet and outlet was periodic condition and for both side, symmetric condition was used. The fully developed turbulent flow was generated and Re = 10,700. The simulated results were compared with DNS data. RANS results were also used for verification.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.10
• /
• pp.1401-1411
• /
• 2003

A two-dimensional direct numerical simulation was performed to investigate the flame structure of CH$_4$$N_2$-air counterflow nonpremixed flame interacting with a single vortex. The detailed transport properties and a modified 16-step augmented reduced mechanism based on Miller and Bowman's detailed chemistry were adopted in this simulation. The characteristic vortex and chemical time scales were introduced to quantify and investigate the extinction phenomenon during a flame-vortex interaction. The results showed that fuel- and air-side vortex cause an unsteady extinction. In this case, the flame interacting with a vortex was extinguished at much larger scalar dissipation rate than steady flame. It was also found that the air-side vortex extinguished a flame more rapidly than the fuel-side vortex. Furthermore, it was noted that the degree of unsteady effect experienced by a flame can be investigated by comparing the above two characteristic time scales, and this analysis could give an appropriate reason for the results of the previously reported experiment.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.8
• /
• pp.1068-1076
• /
• 2001

The relations between wall pressure fluctuations and near-wall streamwise vortices are investigated in a spatially-developing turbulent boundary layer using the direct numerical simulation. The power spectra and two-point correlations of wall pressure fluctuations are presented to validate the present simulation. Emphasis is placed on the identification of the correlation between wall pressure fluctuations and streamwise vorticities. It is shown that wall pressure fluctuations are directly linked with the upstream streamwise vortices in the buffer region of the turbulent boundary layer. The maximum correlation occurs with the spanwise displacement from the location of wall pressure fluctuations. The conditionally-averaged vorticity field and the quadrant analysis of Reynolds shear stress indicate that the sweep events due to streamwise vortices generate positive wall pressure fluctuations, while negative wall pressure fluctuations are created beneath the ejection events and vortex cores. The instantaneous flow field and time records reveal that the rise of high wall pressure fluctuations coincide with the passages of the upstream streamwise vortices.

• Nuclear Engineering and Technology
• /
• v.42 no.6
• /
• pp.608-619
• /
• 2010

A multi-scale analysis of water-cooled reactor thermalhydraulics can be used to take advantage of increased computer power and improved simulation tools, including Direct Numerical Simulation (DNS), Computational Fluid Dynamics (CFD) (in both open and porous mediums), and system thermalhydraulic codes. This paper presents a general strategy for this procedure for various thermalhydraulic scales. A short state of the art is given for each scale, and the role of the scale in the overall multi-scale analysis process is defined. System thermalhydraulic codes will remain a privileged tool for many investigations related to safety. CFD in porous medium is already being frequently used for core thermalhydraulics, either in 3D modules of system codes or in component codes. CFD in open medium allows zooming on some reactor components in specific situations, and may be coupled to the system and component scales. Various modeling approaches exist in the domain from DNS to CFD which may be used to improve the understanding of flow processes, and as a basis for developing more physically based models for macroscopic tools. A few examples are given to illustrate the multi-scale approach. Perspectives for the future are drawn from the present state of the art and directions for future research and development are given.