• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2012.05a
• /
• pp.304-308
• /
• 2012

This paper presents a numerical investigation on propagation of hydrocarbon (ethylene-air mixture) detonation in a deformable copper tube. In this study, we deal with interactions of multi-materials, gas and solid. In gas phase, the model consists of the reactive compressible Navier-Stokes equations and one step chemical reaction. Also we use Inviscid Euler equations in solid. In order to the interface tracking and the determination of boundary values, our model handle level-set and ghost fluid method. Through the numerical simulation results, we identify generations of expansion waves and interferences by the wall deformation. In addition, we predict the minimum copper tube thickness that ensures safety under an incident detonation.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.17 no.2
• /
• pp.114-121
• /
• 2013

In order to estimate the efficiency of an evaporative heat exchanger having mini channel, the equations to calculate heat exchanger properties, those are air temperatures and water temperatures etc, are derived from the governing equations based on the Navier-Stokes equation, even though there are several assumptions to make problem simplify. There are three heat transfer zones at the mini channel heat exchanger depending on the water condition. So, there are three governing equations and solutions to calculate the properties. As the results of this study, the equations to calculate a saturation point and a dry point are derived to evaluate an evaporative heat exchanger having micro channel. It is supposed to predict and evaluate the performance of a mini channel heat exchanger with evaporation of liquid.

• Journal of Korean Society of Forest Science
• /
• v.89 no.1
• /
• pp.116-125
• /
• 2000

This study was conducted to investigate the effects of CEPA(2-chloroethyl phosphonic acid) on bark physiology and anatomy of lactree(Rhus verniciflua Stokes). Sample trees of similar size and growth rate were selected from 7-year-old lactree plantation located in Hyengsung-kun, Kangwon-do. Lanolin pastes containing 0.1, 1, or 10% CEPA were put into the bark-removed hole made by corer(${\phi}1cm$) on the main stem at 1.2m above the ground on June 16, 1995. Five weeks after application of CEPA, bark thickness was markedly increased as a result of the increase in the amount of phloem and intercellular spaces, and correlated with the increased production of urushiol. By the application of 10% CEPA, bark thickness was increased approximately 2.5 times, and the urushiol content within bark was increased 2.8 times compared to that of untreated trees because CEPA stimulated the accumulation of urushiol within bark. Treatment of 10% CEPA also increased the size and the total number of secretory canals, and induced an increase in ray width. The phloem parenchyma cells of CEPA-treated trees were well-developed and closely packed with little intercellular space.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• 2005.10a
• /
• pp.153-159
• /
• 2005

Recently, various novel numerical models based on Navier-Stokes equation rave been developed for calculating wave motions in the sea with coastal or ocean structures. Among those models, Volume Of Fluid (VOF) method might be the most popular one, and it has been used for numerical simulations of wave motions including complicated phenomena of wave breakings. VOF method, however, needs enormous computation time and large computational storage memories in general, thus it is practically difficult to use VOF method for calculations in the case of random waves because long and stable computation ( e.g. for more than 100 significant wave periods) is required to obtain statistically meaningful results. On the other hand of the wave motion is potential motion, Boundary Element Method (BEM), which is a much faster and more accurate method than VOF method, am be effectively used. The aim of this study is to develop a new efficient model applicable to calculations of wave motion and/or wave-structure interactions under random waves. To achieve this, a strictly combined BEM-VOF model has been developed by making the best use of both methods' merits; VOF method is used in a restricted fluid domain around a structure where complicated phenomena of wave breakings may exist, and BEM is used in the other domains far from the disturbance where the wave motion may be assumed to be potential. The verification of the model was performed with numerical results for Stokes'5th order wave propagation and a random wave propagation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.1
• /
• pp.59-69
• /
• 2021

Recently, large-scale offshore and coastal structures have been constructed owing to the increasing interest in eco-friendly energy development. To achieve this, precise simulations of waves are necessary to ensure the safe operations of marine structures. Several experiments are required in the field to understand the offshore wave; however, in terms of scale, it is difficult to control variables, and the cost is significant. In this study, numerical waves under various wave conditions are produced using a piston-type wavemaker, and the produced wave profiles are verified by comparing with the results from a numerical wave tank (NWT) modeled using the smoothed particle hydrodynamics (SPH) method and theoretical equations. To minimize the effect by the reflected wave, a mass-weighted damping zone is set at the right end of the NWT, and therefore, stable and uniform waves are simulated. The waves are generated using the linear and Stokes wave theories, and it is observed that the numerical wave profiles calculated by the Stokes wave theory yield high accuracy. When the relative depth is smaller than two, the results show good agreement irrespective of the wave steepness. However, when the relative depth and wave steepness are larger than 2 and 0.04, respectively, the errors are negligible if the measurement position is close to the excitation plate. However, the error is 10% or larger if the measurement position is away from the excitation location. Applicable target wave ranges are confirmed through various case studies.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.10
• /
• pp.939-946
• /
• 2008

Euler codes or Navier-Stokes codes for compressible flows suffer severe degradation in convergence as Mach number approaches zero. The convergence problem arose from the wide disparity in characteristic speeds can be solved using preconditioning methods without large modifications. In this paper, a preconditioned RANS(Reynolds Averaged Navier-Stokes) solver is developed for analysis of low Mach number flows. In order to validate the method, computational examples are chosen and the results are compared with the experimental data and the existing computed results showing a good accuracy and convergence characteristics for steady inviscid, laminar and turbulent flows at low Mach number.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.6
• /
• pp.8-16
• /
• 2003

Aerodynamic characteristics over Micro Aerial Vehicle(MAV) in low Reynolds number regime are numerically studied using 3-D unsteady, incompressible Navier-Stokes flow solver with single partitioning method for multi-block grid. For more efficient computation of unsteady flows, this flow solver is parallel-implemented with MPl(Message Passing Interface) programming method. Firstly, MAV wing with not complex geometry is considered and then, we analyze aerodynamic characteristics over full MAV configuration varying the angle of attack. Present computational results show a better agreement with the experimental data by MACDL(Micro Aerodynamic Control and Design Lab.), Seoul National University. We can also find the conceptually designed MAV by MACDL has the static stability.

• Aerospace Engineering and Technology
• /
• v.11 no.2
• /
• pp.65-72
• /
• 2012

Numerical simulations were performed to study the stall characteristics of turboprop aircraft. Stall characteristics were qualitatively investigated using the computational results of various configurations based on the combinations of propeller and high lift device. For the analysis of stall characteristics, three-dimensional Navier-Stokes solver with Spalart-Allmaras turbulence model was used and the relative motion between propeller and wing was simulated using sliding mesh technique. For the cruise configurations, major flow separation was occurred at the fuselage/wing fairing and the separation was reduced under propeller slipstream condition. For the high lift device configuration without propeller, major flow separation was occurred at the outboard side of nacelle. With rotating propeller, early stall onset due to low relative velocity and high effective angle of attack was observed on the outboard wing section. Regarding rotating direction of propeller, inboard-down direction was preferred due to the stall delay effect of propeller slipstream.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.45 no.9
• /
• pp.766-774
• /
• 2017

When a plume flow exhausted from a lunar lander descent engine impinges on the lunar surface, regolith particles on the lunar surface will be dispersed due to a plume-surface interaction. If the dispersed particles collide with the lunar lander, some adverse effects such as a performance degradation can be caused. Thus, this study tried to predict the plume flow behaviors using the CFD methods. A nozzle inside region was analyzed by a continuum flow model based on the Navier-Stokes equations while the plume behaviors of the outside nozzle was performed by comparing and analyzing the individual results using the continuum flow model and the DSMC method. As a result, it was possible to establish an optimum procedure of the plume analysis for the lunar lander descent engine in the vacuum condition. In the future, it is expected to utilize the present results for the development of the Korean lunar lander.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.3
• /
• pp.8-14
• /
• 2003

The particle dispersion in the turbulent mixing layer has been numerically investigated to clarify the effect of the velocity ratio in the large-scale vortical structures. In this study the LES with subgrid-scale model is employed. The Lagrangian method to predict the particle motion is applied. The particles of 10, 50, 150, 200${\mu}m$ in mean diameter were loaded into the origin of the mixing layer. It is shown that the characteristics of flow and growth rate are strongly dependent on the variation of the velocity ratio. It is also shown the relationship between the Stokes number and the particle dispersion. As a result, in the case of St~1 the particle dispersion is faster than the diffustion of the flow field while in the cases of both St<<1 and St>>1 it is shown that the particle dispersion in lower than the diffusion of the flow filed.