• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.78-85
• /
• 2012

Cavitation causes a great deal of noise, damage to components, vibrations, and a loss of efficiency in devices, such as propellers, pump impellers, nozzles, injectors, torpedoes, etc., Thus, cavitating flow simulation is of practical importance for many engineering systems. In this study, a two-phase flow solver based on the homogeneous mixture model has been developed. The flow characteristics around an axisymmetric cylinder were calculated and then validated by comparing with the experimental results in the cavitation water tunnel at the Korea Ocean Research & Development Institute. The results show that this solver is highly suitable for simulating the cavitating flows. After the code validation, the cavity length with changes of water depth, angle of attack and velocity were obtained.. Cavitation inception was also calculated for various operational conditions.

• Journal of Mechanical Science and Technology
• /
• v.16 no.3
• /
• pp.386-394
• /
• 2002

The present work is aimed at investigating an unusual variation in flow and performance characteristics of a small propeller fan at low flow rates. A performance test of the fan showed dual performance characteristics, i.e., radial type characteristics at low flow rates and axial type at high flow rates. Dual performance characteristics of the fan are numerically investigated using viscous flow calculations. The Finite Volume Method is used to solve the continuity and Navier-Stokes equations in the flow domain around a fan. The performance parameters and the circumferentially averaged velocity components obtained from the calculations are compared with the experimental results. Numerical values of the performance parameters show good agreement with the measured values. The calculation simulates the steep variations of performance parameters at low flow rates and shows the difference in the flow structure between high and low flow rates. At a low flow coefficient of $\Phi$=0.2, the flow enters the fan in an axial direction and is discharged radially outward at its tip, which is much like the flow characteristics of a centrifugal fan. The centrifugal effect at low flow rates makes a significant difference in performance characteristics of the fan. As the inlet flow rate increases, flow around the fan changes into the mixed type at $\Phi$=0.24 and the axial discharge at $\Phi$=0.4.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.8
• /
• pp.1950-1963
• /
• 1995

A modification of the approximate-factorization method is made to accelerate the convergency rate and to take sufficiently large Courant number without loss of accuracy. And a stable implicit finite-difference scheme for solving the incompressible Navier-Stokes equations employed above modified method is developed. In the present implicit scheme, the volume fluxes with contravariant velocity components and the pressure formulation in curvilinear coordinates is adopted. In order to satisfy the continuity condition completely and to remove spurious errors for the pressure, the Navier-Stokes equations are solved by a modified SMAC scheme using a staggered gird. The upstream-difference scheme such as the QUICK scheme is also employed to the right hand side. The implicit scheme is unconditionally stable and satisfies a diagonally dominant condition for scalar diagonal linear systems of implicit operator on the left hand side. Numerical results for some test calculations of the two-dimensional flow in a square cavity and over a backward-facing step are obtained using both usual approximate-factorization method and the modified one, and compared with each other. It is shown that the present scheme allows a sufficiently large Courant number of O(10$^{2}$) and reduces the computing time.

• Journal of Mechanical Science and Technology
• /
• v.16 no.1
• /
• pp.39-45
• /
• 2002

Initiation and propagation of interfacial crack along bimaterial interface are considered in this study. A biaxial loading device for a single specimen is used for obtaining a wide range of mode-mixities. The specimen is an edge-cracked bimaterial strip of glass and epoxy; the biaxial loading device, being capable of controlling displacements in two perpendicular directions, is developed. A series of interfacial crack initiation and Propagation experiments are conducted using the biaxial loading device for various mixed modes. Normal crack opening displacement (NCOD) is measured near crack front by a crack opening interferometry and used for extracting fracture parameters. From mixed mode interfacial crack initiation experiments, large increase in toughness with shear components is observed. The behavior of interfacial crack propagation analyzed as a function of mode-mix shows that initial crack propagation is delayed with increase of mode-mixity, and its velocity is increased with positive mode-mixity but decreased with negative case. However, it is found that crack propagation is less accelerated with positive mode-mixity than the negative mode-mixity, which may be caused by contact and/or effects of friction between far field and near-tip Held along the interfacial crack.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.6
• /
• pp.819-828
• /
• 2001

The operating point of a small-sized axial fan is strongly dependent upon the system resistance. Therefore, the turbulent flow characteristics around a small-sized axial fan may change significantly according to the operating point. This study represents three-dimensional turbulent flow characteristics around a small-sized axial fan measured at the ideal design point $\phi$=0.25, which is equivalent to the maximum total efficiency point, by using three dimensional fiber-optic type LDA system. This LDA system is composed of a 5 W Argon-ion laser, two optics in back-scatter mode, three BSAs, a PC, and a three-dimensional automatic traversing system. A kind of paraffin fluid is used to supply particles by means of fog generator. Mean velocity profiles downstream of a small-sized axial fan along the radial distance show that the streamwise and the tangential components exist in a predominant manner, while the radial component has a small scale distribution and shows the inflection which its flow direction is inward or outward. Moreover, the turbulent intensity profiles show that the radial component exists the most greatly among turbulent energies.

• Coupled systems mechanics
• /
• v.7 no.4
• /
• pp.407-420
• /
• 2018

In solid-liquid two phase flow, the knowledge of how descending solid particles affected by the presence of downstream wall is important. This work studies at what interstitial distance the velocity of a vertically descending sphere is affected by a downstream wall as a consequence of wall-modified hydrodynamic forces through a validated dynamic model. This interstitial distance-the hydrodynamic coupling distance ${\delta}_c-is$ found to decay monotonically with the approach Stokes number St which compares the particle inertia to viscous drag characterized by the quasi-steady Stokes' drag. The scaling relation ${\delta}_c-St-1$ decays monotonically as literature below the value of St equal to 10. However, the faster diminishing rate is found above the threshold value from St=10-40. Furthermore, an empirical relation of ${\delta}_c-St$ shows dependence on the drop height which clearly indicates the non-negligible effect of unsteady hydrodynamic force components, namely the added mass force and the history force. Finally, we attempt a fitting relation which embedded the particle acceleration effect in the dependence of fitting constants on the diameter-scaled drop height.

• Journal of Institute of Control, Robotics and Systems
• /
• v.2 no.3
• /
• pp.174-180
• /
• 1996

The friction is one of the nonlinearities to be considered in the precise position control of a system which has electromechanical components. The friction has complicated nonlinear characteristics and depends on the velocity, the position and the time. The conventional fixed friction compensator and the controller based on linear control theory may cause the steady state position error or oscillation. The plant to be controlled in this study is a positioning system with a linear brushless DC motor(LBLDCM). The system behaves like a 4th-order model including the compliance and the friction. In this study, the plant model is simplified to a 2nd-order model to reduce the computation in on- line estimation. Also, to reduce the computation time, only the friction is estimated on-line while the mass and the viscous damping coefficient are fixed to the values obtained from off-line estimation. The validity of the proposed scheme is illustrated with the computer simulation and the experiment where the friction is compensated by using the estimation.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2001.09a
• /
• pp.36-43
• /
• 2001

In Korea, we are absolutely short of earthquake data in good quality from moderate and large earthquakes, which are needed fur the study of strong ground motion characteristics. This means that the best use of the available data is needed far the time being. In this respect, several methods are suggested in this paper, which can be applied in the process of data selection and analysis. First, it is shown that the calibration status of seismic stations can be easily checked by comparing the spectra from accelerometer and velocity sensor both of which are located at the same location. Secondly, it is recommended that S/N ratio in the frequency domain should be checked before discarding the data by only look of the data in time domain. Thirdly, the saturated earthquake data caused by ground motion level exceeding the detection limit of a seismograph are considered to see if such data can be used for spectrum analysis by performing numerical simulation. The result reveals that the saturated data can still be used within the dominant frequency range according to the levels of saturation. Finally, a technique to minimize the window effect that distorts the low frequency spectrum is suggested. This technique involves detrending in displacement domain once the displacement data are obtained by integration of low frequency components of the original data in time domain. Especially, the low frequency component can be separated by using discrete wavelet transform among many alternatives. All of these methods mentioned above may increase the available earthquake data and frequency range.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.11 no.2 s.54
• /
• pp.1-9
• /
• 2007

The method of generating the artificial acceleration time histories for seismic analysis based on genetic algorithms is presented. For applying to the genetic algorithms, the frequencies are selected as the decision variables eventually to be genes. An arithmetic average crossover operator and an arithmetic ratio mutation operator are suggested in this study. These operators as well as the typical simple crossover operator are utilized in generating the artificial acceleration time histories corresponding to the specified design response spectrum. Also these generated artificial time histories are checked whether their outward features are to be coincident with the recorded earthquake motion or not. The features include envelope shape, correlation condition between 2 horizontal components of motion, and the relationship of max. acceleration, max. velocity and max. displacement of ground.

• Proceedings of the KSR Conference
• /
• 2007.05a
• /
• pp.1827-1832
• /
• 2007

The purpose of a railway track is to provide a smooth surface for safe and economical train transportation. The performance of the track results from a complex interaction of the track and subgrade components in response to train loading and environmental actions. In the past, the role of subgrade as the track foundation were not recognized adequately. There are insufficient information and inadequate methods for subgrade design, assessment and improvement. This situation has survived for a long time largely because a subgrade defect can often be adjusted by adding more ballast under the ties or applying more frequent track maintenance. Therefore, the application of reinforced roadbed technology will be expected to increase in the future. The reinforced roadbed thickness is set depending on subgrade reaction modulus$(K_{30})$ in the condition of upper subgrade through PBT in both conventional railroad and KTX railroads. As train velocity (V), train passing tonnage (N), and train axial load (P) are not considered in design, the roadbed thickness could be overestimated (or underestimated). Therefore, In this study, the computer model, GEOTRACK, was analyzed the influence of reinforced roadbed thickness factors on track modulus and the characteristics of stress pulses in track and subgrade generated by repeated axle loading.