• Bulletin of the Society of Naval Architects of Korea
• /
• v.19 no.2
• /
• pp.33-39
• /
• 1982

It is one of the basic problems of naval architecture and ocean engineering how to describe the wave kinematics normally under the assumption of an ideal fluid. At present, there are many wave theories available for design purposes. These can be classified into two groups: One is the analytic theory and the other is the numerical theory. This paper briefly introduces the stream function method of R.G. Dean which belongs to the latter group and shows its numerical evaluations exemplified for two cases: One is applied to observed waves and the other is for design waves. In the former case, the wave profiles are calculated by the stream function method and compared with those of the observed waves and also with the results of R.G. Dean. They show good agreement. In the latter case, the wave kinematics and wave loads on a column of diameter 1m are calculated by the stream function method and these are compared with those resulted from the 5th-order gravity wave theory. As a result of comparison the values by the stream function method are slightly larger than those by the 5th-order gravity wave theory but the difference are negligible. From this it is concluded that the stream function method is very useful. And as characteristics of the numerical theories, the stream function method of R.G. Dean can be easily extended to the higher order terms and can include easily the current velocity and the pressure distribution on the free surface. In addition, when the data of observed wave profile are given, this method can reproduced the observed wave profile as closely as possible so that this method seems to describe the ocean wave more realistically. And from standpoint of a mathematical principle the stream function method exactly satisfies the kinematic free-surface boundary condition.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.3 no.4
• /
• pp.21-31
• /
• 1983

Circulation phenamena in open channel with abrupt variation in flow section are studied by solving numerically two dimensional Navier-Stokes equations integrated over depth. Galerkin type finite element method is used as numerical scheme. Numerical results by both implicit and explicit schemes tested in one-demensional rectangular channel agree closely with the known solution. The numerical experiments carded out in the open channel with a pool indicate the expected flow pattern and the center of the circulation coincides with the geometrical center, but the vectors of velocity appear father small, and it remains to be further investigated. Numerically simulated flow profiles along the channel with constrictions such as bridge piers and abutments are shown to be close to hydraulic experimental results. Thus further refined numerical technique is expected to be able to serve as a tool to evaluate the effect of bridge backwater.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.3
• /
• pp.301-307
• /
• 2012

The characteristics of pressure loss in an asymmetric bifurcating tube were investigated numerically for steady inspiratory conditions. The loss coefficient K calculated for various asymmetry and flow-distribution ratios found in human lung airways showed a power-law dependence on the Reynolds number (Re) and length-to-diameter ratio (L/d), with different exponents for Re $\geq$ 100 and Re < 100. The fundamental characteristics of the asymmetric bifurcation are similar to the case of symmetric bifurcation. In addition, the effect of skewed inlet velocity profiles on the pressure loss was weak, and decreased with an increasing number of bifurcations.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1940-1954
• /
• 1992

Fine grid computations were attempted to analyze the turbulent flows in the near wall low Reynolds number region and the numerical analyses were incorporated by a finite-volume discretization with full find grid system and low Reynolds number k-.epsilon. model was employed in this region. For the improvement of low Reynolds number k-.epsilon. model, modification coefficient of eddy viscosity $f_{\mu}$ was derived as a function of turbulent Reynolds number $R_{+}$ and nondimensional length $y^{+}$ from the concept of two length scales of dissipation rate of turbulent kinetic energy. The modification coefficient $f_{\epsilon}$ in .epsilon. transport equation was also derived theoretically. In the turbulent kinetic energy equation, pressure diffusion term was added in order to consider low Reynolds number region effect. The main characteristics of this low Reynolds number k-.epsilon. model were founded as : (1) In high Reynolds number region, the present model has limiting behavior which approaches to the high Reynolds number model. (2) Present low Reynolds number k-.epsilon. model dose not need additional empirical constants for the transport equations of turbulent kinetic energy and dissipation of turbulent kinetic energy in order to consider wall effect. Present low Reynolds number turbulence model was tested in the pipe flow and obtained improved results in velocity profiles and Reynolds stress distributions compared with those from other k-.epsilon. models.s.s.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.8
• /
• pp.20-28
• /
• 2004

For the prediction on the onset of oblique shock wave-induced vortex breakdown, computational studies on the Oblique Shock wave/Vortex Interaction (OSVI) are conducted and compared with both experimental results and analytic mode1. A Shock-stable numerical scheme, the Roe scheme with Mach number-based function (RoeM), and a two-equation eddy viscosity-transport approach arc used for three-dimensional turbulent flow computations. The computational configuration is identical to available experiment, and we attempt to ascertain the effect of parameters such as a vortex strength, streamwise velocity deficit, and shock strength at a freestream Mach number of 2.49. Numerical simulations using the k-w SST turbulence model and suitably modeled vortex profiles are able to accurately reproduce many fine features through a direct comparison with experimental observations. The present computational approach to determine the criterion on the onset of oblique shock wave-induced vortex breakdown is found to be in good agreement with both the experimental result and the analytic prediction.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.8
• /
• pp.1577-1584
• /
• 2002

This paper develops a finite element model for studying the occupant behavior and injury cofficients of a large-sized cab-over type truck. Since it does not have a room to absorb collision energy and deformation in front of the passenger compartment the deformation is directly transmitted to the passenger compartment. Moreover, since its steering column is attached on the frame, severe deformation of the frame directly affects on the steering wheel's movement. Therefore, if the occupant behavior and injury coefficients analysis is performed using a finite element model developed based on a sled test, it is very difficult to expect acquiring satisfactory results. Thus, the finite element model developing in this paper is based on the frontal crash test in order to overcome the inherent problems of the sled test based model commonly used in the passenger car. The occupant behavior and injury coefficients analysis is performed using PAM-CRASH installed in super-computer SP2. In order to validate the reliability of the developed finite element model, a frontal crash test is carried out according to a test method used fur developing truck occupant's secondary safety system in european community and japan. That is, test vehicle's collision direction is vertical to the rigid barrier and collision velocity is 45kph. Thus, measured vehicle pulses at the lower parts of the left and right B-pilla., dummy chest and head deceleration profiles, HIC(head injury criterial) and CA(chest acceleration) values, and dummy behavior from the frontal crash test are compared to the analysis results to validate reliability of the developed model.

• Journal of the Korean Society of Visualization
• /
• v.1 no.1
• /
• pp.98-106
• /
• 2003

A PLLIF (Planar Liquid Laser Induced Fluorescence) technique has been known to be a useful tool for the measurement of the spray patterns for various spray injectors because it can obtain two-dimensional images with high spatial resolutions without any intrusion on the spray field. In case of dense spray, however, the secondary emission as well as the extinction of an incident laser beam or a fluorescence signal can cause errors in quantifying a mass distribution. Unfortunately, a like-doublet injector which has a dense spray zone at the center may not be a good example or the application of the PLLIF technique. Therefore, we took PLLIF data for the like-doublet injector with a 12 bit color CCD camera by varying laser power, and then assessed their accuracy by comparing with the data obtained with a mechanical patternator and a PDPA (Phase Doppler Particle Analyzer). The experimental results showed that the gray level of fluorescence signal increases nonlinearly due to a secondary emission at the dense spray zone but this nonlinearity can be avoided by reducing the incident laser beam power. In addition, the mass flux distribution of the spray could be obtained by using the mass concentration data from PLLIF technique and the velocity profiles of liquid drops, and this distribution showed good agreement with that of mechanical pattemator. Therefore, it is possible that the PLLIF technique can be successfully applied to finding the mass distributions of like-doublet injectors.

• Korean Journal of Applied Biomechanics
• /
• v.21 no.3
• /
• pp.289-296
• /
• 2011

The purpose of this study was to quantify the biomechanical characteristics of the ground reaction force(GRF) during the Taekwondo's Apkubi, one of the basic movement in Taekwondo and the walking. The GRF profiles under the stance foot of Apkubi movement and walking were directly measured in sample of 20 healthy older persons. In the anterior-posterior and vertical direction, the GRF of the Apkubi movement reached to the peak braking force at 10% of the normalized stance time percent and the peak driving force at 90% of stance time, but that of the walking reached to the peak braking force at 20% of stance time and the peak driving force at 80% of stance time. In vertical force, the GRF of the walking showed two peak values, but that of the Apkubi movement seemed three peak values. Moreover the first peak vertical force was significantly(t=6.085, p<.001) greater in the walking(about 1.8 times of body weight) than the Apkubi(about 1.4 times of body weight). The walking velocity was affected significantly(over p<.05) by the braking impulse, the peak braking force and the first peak vertical force. Futhermore the peak braking force in the Apkubi showed a significant effect on the Apkubi's stride length(p<.01). So, we concluded that the braking force after the right touch down, the stance foot on the ground contributed to move the leg forward.

• Journal of the Korean Institute of Gas
• /
• v.7 no.3 s.20
• /
• pp.44-50
• /
• 2003

The axisymmetric methane-air counterflow flame was simulated to investigate changes in the flame structure due to the fuel concentration and to evaluate the numerical method. The global strain rates $a_g=20,\;60,\;90\;s^{-1}$ and the mole fractions of methane $x_m=20,\;50,\;80\%$ in the fuel stream were taken to be numerical parameters. The axisymmetric simulation was conducted by using the Fire Dynamics Simulator (FDS) which employed a mixture fraction combustion model, and the results were compared with those of OPPDIF, which is an one-dimensional flamelet code and includes detail chemical reactions. In all the cases tested, there was good agreement in the temperature and axial velocity profiles between the axisymmetric and one-dimensional simulations. It was shown that the flame thickness and peak flame temperature increase and the flame radius decreases as the fuel concentration increases.

• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.1-10
• /
• 2016

The analysis of characteristics of turbulent flow and thermal boundary layer for natural convection caused by fire along vertical wall is performed. The 4m-high vertical copper plate is heated and kept at a uniform surface temperature of $60^{\circ}C$ and the surrounding fluid (air) is kept at $16.5^{\circ}C$. The flow and temperature is solved by large eddy simulation(LES) of FDS code(Ver.6), in which the viscous-sublayer flow is calculated by Werner-Wengle wall function. The whole analyzed domain is assumed as turbulent region to apply wall function even through the laminar flow is transient to the turbulent flow between $10^9$<$Gr_z$<$10^{10}$ in experiments. The various grids from $7{\times}7{\times}128$ to $18{\times}18{\times}128$ are applied to investigate the sensitivity of wall function to $x^+$ value in LES simulation. The mean velocity and temperature profiles in the turbulent boundary layer are compared with experimental data by Tsuji & Nagano and the results from other LES simulation in which the viscous-sublayer flow is directly solved with many grids. The relationship between heat transfer rate($Nu_z$) and $Gr_zPr$ is investigated and calculated heat transfer rates are compared with theoretical equation and experimental data.