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

• Journal of the Korean Geotechnical Society
• /
• v.28 no.4
• /
• pp.91-100
• /
• 2012

Soil-water characteristic curves (SWCCs), which represent a physical property in partially saturated soils, show the relation between volumetric water content and matric suction. The SWCCs exhibit hysteresis during wetting and drying, however experimental expressions used to describe SWCCs have generally ignored the hysteresis. In addition, the shape of SWCC may depend on the void ratio which is changed by soil skeleton deformations or hysteretic behavior under various loading conditions. Hence, it is necessary to understand, both empirically and analytically, the relationship between soil skeleton deformations and the SWCCs of various soils. The typical SWCCs experimentally have drying, wetting, and the second drying curve. The measurement of a complete set of hysteretic curves is severely time-consuming and difficult works, then the first drying curve of SWCC is generally determined to estimate the hydraulic conductivity and shear strength function of partially saturated soils. This paper presents the hydraulic-mechanical behavior of partially saturated soils (weathered soil and silty soil) for volume changes and hysteresis in SWCCs regarding the difference between the first drying and wetting curve.

• Polymer(Korea)
• /
• v.29 no.6
• /
• pp.549-556
• /
• 2005

Effects of the dimensionless variables on the heat transport phenomena in the extrusion process of a single screw extruder have been studied numerically. Based on the understanding of the solids conveying related to the geometrical structure and characteristics of the screw, the heat balance equation for the solids conveying zone was established and normalized. The finite volume method and power-law scheme were applied to derive a discretized equation and the equation was solved using the alternating direction iterative method with relaxation. Effects of the dimensionless parameters, Biot and Peclet numbers, that define the heat transfer characteristics of the solids conveying zone have been investigated with respect to the temperature of the feeding zone and the length of the solids conveying zone. As the Biot number is increased, the heat loss by cooling dominates to decrease the temperature of the barrel but it has little effects on the temperature of the solids bed and the length of the solids conveying zone. On the other hand, if the Peclet number is increased, the convection term dominates to decrease the temperature of the solids bed and it results in an increase in the length of the solids conveying zone.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.4
• /
• pp.771-781
• /
• 1994

The heat generation of cement causes the internal temperature rise and volume change at early age, particulary in massive concrete structures. As the results of the temperature rise and restraint conditions, the thermal stress may induce cracks in concrete. Therefore, the prediction of the thermal stress is very important in the design and construction in order to control the cracks developed in mass concrete. In case of young concrete, creep effect by the temperature load is larger than that of old concrete. Thus, the effect of creep must be considered for checking the cracks, serviceability, durability and leakage. This paper is concentrated on the development of a finite element program which is capable of simulating the temperature history and the thermal stress considering creep and the modified elastic modulus due to inner temperature change and maturity. The analytical results in the inner parts highest important to control cracks are in good agreement with experimental data. Therefore this study may provide available method to control the cracks.

• Magazine of the Korea Concrete Institute
• /
• v.7 no.6
• /
• pp.176-185
• /
• 1995

The heat generation of hydration of cement causes the internal temperature rise and volum& change at early age, particularly in massive concrete structures. As the results of the temperature rise and external restraint conditions, the thermal stress may induce cracks in concrete. Therefore, various techniques of the thermal stress control of the mass concrete have been wid'dy used. One of them is pipecooling which reduces the temperature of concrete with flowing water. The objective of this paper' is to develop a finite element program which is capable of simulating the temperature history considering pipe-cooling effect. The numerical results in this study are in good agreement with experimental data measured in the footing(l1 x22m). Therefore, this study may provide available method to predict the hydration temperature of concrete with pip:-cooling.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1149-1164
• /
• 1997

In this study, the applicability of the RNG k-.epsilon. model to the analysis of the complex flows is studied. The governing equations based on a non-orthogonal coordinate formulation with Cartesian velocity components are used and discretized by the finite volume method with non-staggered variable arrangements. The predicted results using the RNG k-.epsilon. model of three complex flows, i.e., the flow over a backward-facing step and a blunt flat plate, the flow around a 2D model car are compared to these from the standard k-.epsilon. model and experimental data. That of the unsteady axisymmetric turbulent flow within a cylinder of reciprocating model engine including port/valve assembly and the spray characteristics within a chamber of direct injection model engine are compared to these from the standard k-.epsilon. model and experimental data. The results of reattachment length, separated eddy size, average surface pressure distribution using the RNG k-.epsilon. model show more reasonable trends comparing with the experimental data than those using the modified k-.epsilon. model. Although the predicted rms velocity using the modified k-.epsilon. model is lower considerably than the experimental data in incylinder flow with poppet valve, predicted axial and radial velocity distributions at the valve exit and in-cylinder region show good agreements with the experimental data. The spray tip penetration predicted using the RNG k-.epsilon. model is more close to the experimental data than that using the modified k-.epsilon. model. The application of the RNG k-.epsilon. model seems to have some potential for the simulations of the unsteady turbulent flow within a port/valve-cylinder assembly and the spray characteristics over the modified k-.epsilon. model.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.1
• /
• pp.29-41
• /
• 2008

The present paper provides a CFD analysis of diffraction problem for a ship with forward speed using an unsteady RANS simulation method, a WAVIS code. The WAVIS viscous solver adopting a finite volume method has second order accuracy in time and field discretizaions for the RANS equations. A two phase level-set method and a realizable ${\kappa}-{\varepsilon}$ turbulence model are adopted to compute the free surface and to meet the turbulence closure, respectively. To validate the capability of the present numerical methods for the simulation of an unsteady progressive regular wave, computations are performed for three grid sets with refinement ratio of ${\sqrt{2}}$. The main simulation is performed for a DTMB5512 model with a forward speed in a regular head sea condition. Validation of the present numerical method is carried out by comparing the present CFD results with available unsteady experimental data published in the 2005 Tokyo CFD Workshop: resistance, heave force, pitch moment, unsteady free surface elevations and velocity fields.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.7
• /
• pp.629-636
• /
• 2013

Pneumatic levitation is based upon Bernoulli's principle. However, this method is known to require a large gas flow rate that can lead to an increase in the cost of products. In this case, the gas flow rate should be increased, and the compressible effects of the gas may be of practical importance. In the present study, a computational fluid dynamics method has been used to obtain insights into Bernoulli levitation flows. Three-dimensional compressible Navier-Stokes equations in combination with the SST k-${\omega}$ turbulence model were solved using a fully implicit finite volume scheme. The gas flow rate, workpiece diameter,and clearance gap between the workpiece and the circular cylinder were varied to investigate the flow characteristics inside. It is known that there is an optimal clearance gap for the lifting force and that increasing the supply gas flow rate results in a larger lifting force.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.