• Journal of the Korean Geotechnical Society
• /
• v.32 no.9
• /
• pp.51-62
• /
• 2016

Parametric studies for various site conditions by using 3d numerical model were carried out in order to estimate dynamic behavior of soil-pile-structure system in dry soil deposits. Proposed model was analyzed in time domain using FLAC3D which is commercial finite difference code to properly simulate nonlinear response of soil under strong earthquake. Mohr-Coulomb criterion was adopted as soil constitutive model. Soil nonlinearity was considered by adopting the hysteretic damping model, and an interface model which can simulate separation and slip between soil and pile was adopted. Simplified continuum modeling was used as boundary condition to reduce analysis time. Also, initial shear modulus and yield depth were appropriately determined for accurate simulation of system's nonlinear behavior. Parametric study was performed by varying weight of superstructure, pile length, pile head fixity, soil relative density with proposed numerical model. From the results of parametric study, it is identified that inertial force induced by superstructure is dominant on dynamic behavior of soil-pile-structure system and effect of kinematic force induced by soil movement was relatively small. Difference in dynamic behavior according to the pile length and pile head fixity was also numerically investigated.

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

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.17 no.3
• /
• pp.341-351
• /
• 2015

When rainfall occurred on road tunnel, that is likely to have influence upon ventilation force in the tunnels but the tunnels ventilation system did not consider factors of rainfall. Thus, this study investigated effects of rainfall upon ventilation force in the tunnels at no rainfall and changing of rainfall by 3 dimensional numerical method. Flow rate into road tunnels decreased as many as 52.34% at rainfall of 150 mm/hr, and pressure drop of road tunnel between entrance and exit decreased as many as 22.22%, so that rainfall had influence upon ventilation force in the tunnel. The number of necessary jet fan in road tunnels is 12 at no rainfall but, when rainfall of 80 mm/hr on road tunnels, the number of necessary jet fan in road tunnels is 16, when rainfall of 150 mm/hr on road tunnels, the number of necessary jet fan in road tunnels is 17. So, factor of rainfall should be considered at estimation of ventilation system of road tunnel.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.23 no.4
• /
• pp.10-20
• /
• 2019

Recently, fluidic thrust vector control has become a core technique to control multifarious air vehicles, such as supersonic aircraft and modern rockets. Fluidic thrust vector control using the shock vector concept has many advantages for achieving great vectoring performance, such as fast vectoring response, simple structure, and low weight. In this paper, computational fluid dynamics methods are used to study a three-dimensional rectangular supersonic nozzle with a slot injector. To evaluate the reliability and stability of computational methodology, the numerical results were validated with experimental data. The pressure distributions along the upper and lower nozzle walls in the symmetry plane showed an excellent match with the test results. Several numerical simulations were performed based on the shear stress transport(SST) $k-{\omega}$ turbulence model. The effect of the momentum flux ratio was investigated thoroughly, and the performance variations have been clearly illustrated.

• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.241-249
• /
• 2006

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms commonly used, however, may not give sufficiently precise computational results of traveltime data particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. Considering the complex geology of Korea, we assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.3
• /
• pp.31-39
• /
• 1999

In this paper, numerical experiment is attempted to analyze and compare the combustion efficiency of the burning sprays due to OFO, FOF triplet / FOOF split doublet injectors. Preconditioned Wavier-Stokes equation system with low Reynolds number $\kappa$-$\varepsilon$ model for turbulence closure, is LU-SGS time-integrated. Spray processes are modeled by DSF analysis with experimentally determined injection characteristics. n-heptane/air global reaction model approximates the combustion for simplicity, and the influence of turbulence on the chemical reaction is included using eddy dissipation model. The results showed the FOF triplet injector of highest combustion efficiency, whereas the OFO type of poet performance. It was also observed that the droplet mean diameter and the average gas temperature due to the mixing efficiency, are the representative parameters for the performance design of combustion.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.1
• /
• pp.39-46
• /
• 1993

The need to develop method which can improve the shape design efficiency using high quality approximation is being brought up. In this study, to perform shape optimal design of three-dimensional continuum structures an efficient approximation method for stress constraints is proposed, based on expanding the nodal forces in Taylor series with respect to shape variables. Numerical examples are performed using the 3-D cantilever beam and fixed-fixed beam and compared with other method to demonstrate the efficiency and convergence rate of the Force Approximation method. It is shown that by taking advantage of this high quality approximation, the total number of finite element analysis required for shape optimization of 3-D continuum structures can be reduced significantly, resulting to the same level of efficiency achieved previously in sizing optimization problems. Also, shape representation by super curve technique applied to obtain optimal shape finds useful method.

• Tunnel and Underground Space
• /
• v.15 no.2 s.55
• /
• pp.145-156
• /
• 2005

MODFLOW, 3-D finite difference code, is widely used to model groundwater flow and has been used to assess the effect of excavations on the groundwater system due to construction of subways and mountain tunnels. The results of numerical analysis depend on boundary conditions, initial conditions, conceptual models and hydrogeological properties. Therefore, its accuracy can only be enhanced using more realistic and field oriented input parameters. In this study, SA(simulated annealing) was used to integrate hydraulic conductivities from a few of injection tests with geophysical reference images. The realizations of hydraulic conductivity random field are obtained and then groundwater flows in each geostatistically equivalent media are analyzed with a numerical simulation. This approach can give probabilistic results of groundwater flow modeling considering the uncertainty of hydrogeological medium. In other words, this approach makes it possible to quantify the propagation of uncertainty of hydraulic conductivities into groundwater flow.

• Journal of Advanced Marine Engineering and Technology
• /
• v.37 no.4
• /
• pp.344-350
• /
• 2013

The purpose of this study is to examine the structural stability of a low speed 200 W class gyromill type vertical axis wind turbine system. For the analysis, a commercial code is adopted. The pressure distribution on the rotor blade surface is examined in detail. In order to perform unidirectional FSI(Fluid-Structure Interaction) analysis, the pressure resulted from CFD analysis has been mapped on the surface of wind turbine as load condition. The rotational speed and gravitational force of wind turbine are also considered. The results of FSI analysis show that the wind turbine reveals an enough structural margin. The maximum structural displacement occurs at trailing edge of blade and the maximum stress occurs at the strut.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.3
• /
• pp.237-250
• /
• 2003

Since an orthotropic steel deck bridge has large number of design variables and shows complex structural behavior, it would be very difficult and impractical to directly use a Conventional Single Level (CSL) optimization algorithm for its optimum design. Thus, in this paper, an Improved Multi Level Design Synthesis (IMLDS) optimization algorithm is proposed to improve the computational efficiency. In the proposed IMLDS algorithm, a coordination method is introduced to divide the bridge into main girders and orthotropic steel deck with preserving the characteristics of the structural behavior. For an efficient optimization of the bridge, the IMLDS algorithm incorporates the various crucial approximation techniques such as constraints deletion, Automatic Differentiation (AD), stress reanalysis, and etc. In the case of orthotropic steel deck system, optimum design problems are characterized by mixed continuous discrete variables and discontinuous design space. Thus, a modified Genetic Algorithm (GA) is also applied to optimize discrete member design for orthotropic steel deck. From the numerical example, the efficiency and convergency of the IMLDS algorithm proposed in this paper is investigated. It may be positively stated that the IMLDS algorithm will lead to more effective and practical design compared with previous algorithms.