• Journal of the Korean Geosynthetics Society
• /
• v.15 no.2
• /
• pp.45-54
• /
• 2016

Double deck tunnels are considered to have a large demand in a near future for solving traffic congestion problems and overcoming the limitations in constructing new tunnels. This study presents a numerical investigation using finite element (FE) analysis on the behaviors of the tunnels and the stability of pillars in a divergence section where single tunnel is diverged from a main line double deck tunnel. The effects of the separation distance between the diverged and the main tunnels and the ground condition were examined through the FE analysis by varying the separation distance from 0.1D to 2.0D (D: diameter of main tunnel) and the rock class from class I to V, respectively, and the analysis results were compared with those using empirical methods, strength-stress ratio, and the volume of interference. The FE analysis results indicated that the separation distance has a larger effect on tunnel behaviors, compared with the rock strength, and a single tunnel with a large cross section is more favorable than two separated tunnels for tunnel stability when the separation distance is below 0.7D.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.35 no.4
• /
• pp.353-358
• /
• 2011

In this paper, we describe a process for optimally designing a ring-type permanent magnet thrust bearing. The bearing consists of two sets of permanent magnet rings. One set is located inside the other set. An axial offset between the two sets creates axial force, which results in a thrust bearing function. In order to realize an optimal design of the bearing where the required load capacity of the bearing is achieved with the least magnet volume, we derived analytical design equations by adopting the equivalent current sheet (ECS) method. We considered the following two types of magnet arrays: axial arrays and Halbach arrays. These two types of arrays are optimized using the analytical design equations. The results of the optimization are verified using three dimensional (3D) finite element analyses (FEA). The results show that the Halbach array can achieve the required load capacity with less amount of permanent magnet than the axial array does. The efficacy of the ECS method is also verified by using 3D FEA. It is found that the accuracy of ECS method is more sensitive to the underlying assumptions for the Halbach array than for the axial array.

• Journal of Biomedical Engineering Research
• /
• v.23 no.6
• /
• pp.467-475
• /
• 2002

Hyperthermia using transrectal thermal probes has been used for a noninvasive treatment of prostate diseases. However it is known that heating the rectal wall at excessively high temperature can lead to destruction of the rectal mucous membrane. and it is difficult to maintain an optimum temperature over the entire prostate. Thus, a more accurate understanding of the heat transfer mechanism between prostate and hyperthermia system is needed Numerical analysis was performed to investigate how the cold/warm stimulations on the prostate surface affect the temperature distribution in the prostate model. The general purpose software "FLUENT" was used for obtaining a finite volume solution to the unsteady conduction equation and to calculate the time-varying temperature in the prostate. Effects of the warm/cold stimulations and the stimulation frequency on the temperature distribution were simulated. and we visualized how hyperthermia affected the inside of the prostate. It was found that the effect of hyperthermia by using a typical heating method is limited due to the low thermal conductivity of the prostate. Consecutive repetitions of warm and cold stimulations were considered to provide the thermal irritations inside a prostate. The effects of temperature difference and duration of warm/cold stimulations were investigated, and basic data for the optimum period and effective patterns of stimulations were obtained. A simplified bioheat equation was also solved to describe effects of the blood flow on the blood-tissue heat transfer. The effect of blood flow was not dominant compared to that of warm/cold stimulations. These results might be used as data for design of prostate treating probe, prostatic therapy and thermal stimulation effects on the prostate.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.36 no.4
• /
• pp.329-336
• /
• 2000

Welding is a important technological method in mechanical engineering. $CO_2$MAG(metal active gas) welding means that metal part in double capstan drum for the inshore and costal vessels are joined by melting(with or without a filler material) or that new material is added to a metal part by melting. The thermal stresses appear due to a non-uniform temperature field, inhomogeneous material properties, external restraint and volume changes during phase transformations. In this study analysis the elasto-plastic thermal stresses distribution of welding part in double capstan drum for the inshore and costal vessels using finite element method (FBM). Therefore it calculates the numerical value that can be applied to the optimum design of welding parts and the shapes. The significant results obtained in this study are summarized as fellows. At early stage of the cooling after welding process, the abrupt thermal stresses gradient has been shown in the vicinity of welding part. In the thermal stresses analysis due to temperature gradient and heat shocking maximum stress was occurred of welding part and stresses were distributed from 54MPa~48MPa.

• Journal of Korea Water Resources Association
• /
• v.50 no.2
• /
• pp.89-97
• /
• 2017

Turbulent flow structure in the high amplitude meandering channel is complex due to secondary recirculation with helicoidal motions and shear layers formed by flow separation from the curved sidewall. In this work, the secondary flow and the superelevation of the water surface produced in the high-amplitude Kinoshita channel are reproduced by the unsteady Reynolds-averaged Navier-Stokes (RANS) computations using the VOF technique for resolving the variation of water surface elevation and three statistical turbulence models ($k-{\varepsilon}$, RNG $k-{\varepsilon}$, $k-{\omega}$ SST). The numerical results computed by a second-order accurate finite volume method are compared with an existing experimental measurement. Among applied turbulence models, $k-{\omega}$ SST model relatively well predicts overall distribution of the secondary recirculation in the Kinoshita channel, while all three models yield similar prediction of water superelevation transverse slope. The secondary recirculation driven by the radial acceleration in the upstream bend affects the flow structure in the downstream bend, which yields a pair of counter-rotating vortices at the bend apex. This complex flow pattern is reasonably well reproduced by the $k-{\omega}$ SST model. Both $k-{\varepsilon}$ based models fail to predict the clockwise-rotating vortex between a pair of counter-rotating vortices which was observed in the experiment. Regardless of applied turbulence models, the present computations using the VOF method appear to well reproduce the superelevation of water surface through the meandering channel.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.2
• /
• pp.63-70
• /
• 2014

Design sensitivity analysis and topology design optimization for a piezoelectric crystal resonator are developed. The piezoelectric crystal resonator is deformed mechanically when subjected to electric charge on the electrodes, or vice versa. The Mindlin plate theory with higher-order interpolations along thickness direction is employed for analyzing the thickness-shear vibrations of the crystal resonator. Thin electrode plates are masked on the top and bottom layers of the crystal plate in order to enforce to vibrate it or detect electric signals. Although the electrode is very thin, its weight and shape could change the performance of the resonators. Thus, the design variables are the bulk material densities corresponding to the mass of masking electrode plates. An optimization problem is formulated to find the optimal topology of electrodes, maximizing the thickness-shear contribution of strain energy at the desired motion and restricting the allowable volume and area of masking plates. The necessary design gradients for the thickness-shear frequency(eigenvalue) and the corresponding mode shape(eigenvector) are computed very efficiently and accurately using the analytical design sensitivity analysis method using the eigenvector expansion concept. Through some demonstrative numerical examples, the design sensitivity analysis method is verified to be very efficient and accurate by comparing with the finite difference method. It is also observed that the optimal electrode design yields an improved mode shape and thickness-shear energy.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.19 no.2
• /
• pp.265-281
• /
• 2017

For shallow tunnel constructions, settlement of the ground surface is a main issue. Recent technical developments in shield TBM tunneling technique have enabled a decrease in such settlements based on tunneling with ground deformation controls. For this objective, the tail void grouting is a common practice. Generally surface settlements in a soil of low permeability occur during a tunnel construction but also during a long period after completion of the tunnel. The long-term settlements occur mainly due to consolidation around the tunnel. The consolidation process is caused and determined by the tail void grouting which leads to an excess pore water pressure in the vicinity of the tunnel. Because of this, the grouting pressure has a strong effect on the long-term settlements in the shield tunneling. In order to investigate this effect, a series of coupled hydro-mechanical 3D finite element simulations have been performed. The results show that an increase in grouting pressure reduces the short-term settlements, but in many cases, it doesn't lead to a reduction of the final settlements after the completion of consolidation. Thereby, the existence of a critical grouting pressure is identified, at which the minimal settlements are expected.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.5
• /
• pp.441-449
• /
• 2015

In this paper we perform a linearized buckling analysis using the Kirchhoff plate theory and the von Karman nonlinear strain-displacement relation. Design sensitivity analysis(DSA) expressions for plane elasticity and buckling problems are derived with respect to Young's modulus and thickness. Using the design sensitivity, we can formulate the topology optimization method for minimizing the compliance and maximizing eigenvalues. We develop a topology optimization method applicable to plate buckling problems using the prestress for buckling analysis. Since the prestress is needed to assemble the stress matrix for buckling problem using the von Karman nonlinear strain, we introduced out-of-plane motion. The design variables are parameterized into normalized bulk material densities. The objective functions are the minimum compliance and the maximum eigenvalues and the constraint is the allowable volume. Through several numerical examples, the developed DSA method is verified to yield very accurate sensitivity results compared with the finite difference ones and the topology optimization yields physically meaningful results.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.1
• /
• pp.11-23
• /
• 1997

This paper contains a new approach to blade section design method for marine propellers. The hydrodynamic characteristics of 2-D section are highly influenced by its geometrical parameters i.e., thickness and camber distributions and leading edge radius etc. To consider fully turbulent flow field near 2-D section. the finite volume method with k-${\varepsilon}$ turbulent model which solve Reynolds time averaged Navier-Stokes(RANS) equation is applied. In this study, O-type grid system that can provide many calculation points on blade surface is used. The results were compared with those of the experiment of NACA0012 to confirm the accuracy of the developed codes. The goal of this study is the development of a blade section with high efficiency and low drag. To achieve this, we carried out the tests of lift, drag and cavitation characteristics in cavitation tunnel. The results of experiment were compared with numerical results in order to validate the proposed blades design method. By comparing the numerical results with the experiments, we found that the new blade section, KH28 allows superior performance in efficiency and cavitation avoidance characteristics. We further investigated the blade section design method and an application study of this section, KH28 to apply to the marine propeller. In order to improve the accuracy of numerical results on prediction of lift and drag, we conclude here that the 2-layer boundary model must be used.

• Journal of the Society of Naval Architects of Korea
• /
• v.55 no.3
• /
• pp.252-264
• /
• 2018

The calculation of steady-state cavitating flows around Supercavitating Underwater Bodies (SUB's), which consist of a circular disk head (cavitator), a conical fore-body, a cylindrical middle-body and either a boat-tail or a flare-tail, are carried out. To calculate the axisymmetric cavitating flow, used is a commercial computational fluid dynamics code based on the finite volume method, Fluent. From the analysis of numerical results, the cavity and drag, affected by the fore-body and tail of the SUB's, are investigated. Firstly, the effect of the fore-body shape is investigated with the same disk cavitator and a cylindrical rear-body of fixed diameter. Then with the same cavitator and a fixed fore-body, the effect of the rear-body shape is investigated. Before the cavity generated by the cavitator covers the slant of fore-bodies sufficiently, the larger the cone angle of the fore-body(i.e., the shorter the slant length), the larger the drag and the slower the development of cavity. After the cavity covers the fore-body completely so that the pressure drag component of the body is vanished, the characteristics of drag-velocity curves are identical. Also, as the tail angle is bigger, the cavity generated by the cavitator is suppressed further and the drag becomes larger. The peak of the drag appears for the flare-tail, i.e., when the tail angle is positive(+). On the contrary, the trough of the drag appears for the boat-tail, i.e., when the tail angle is negative(-). When the tail angle is 5 degrees, the peak of the drag appears at the body speed of 80m/s and the value of the drag is 43% larger than that at the design speed of 100m/s. When the tail angle is -5 degrees, the trough of the total drag appears at 75m/s and that drag is 30% smaller than that of the cavitator, which means the rest of the body has a negative drag.