• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.255-261
• /
• 2003

Tunnelling in poor and heterogeneous ground is a difficult task. Even with a good geological investigation, uncertainties with respect to the local rock mass structure will remain. Especially for such conditions, a reliable short-term prediction of the conditions ahead and outside the tunnel profile are of paramount importance for the choice of appropriate excavation and support methods. The information contained in the absolute displacement monitoring data allows a comprehensive evaluation of the displacements and the determination of the behaviour and influence of an anisotropic rock mass. Case histories and with numerical simulations show, that changes in the displacement vector orientation can indicate changing rock mass conditions ahead of the tunnel face (Schubert & Budil 1995, Steindorfer & Schubert 1997). Further research has been conducted to quantify the influence of weak zones on stresses and displacements (Grossauer 2001). Sellner (2000) developed software, which allows predicting displacements (GeoFit$\circledR$). The function parameters describe the time and advance dependent deformation of a tunnel. Routinely applying this method at each measuring section allows determining trends of those parameters. It shows, that the trends of parameter sets indicate changes in the stiffness of the rock mass outside the tunnel in a similar way, as the displacement vector orientation does. Three-dimensional Finite Element simulations of different weakness zone properties, thicknesses, and orientations relative to the tunnel axis were carried out and the function parameters evaluated from the results. The results are compared to monitoring results from alpine tunnels in heterogeneous rock. The good qualitative correlation between trends observed on site and numerical results gives hope that by a routine determination of the function parameters during excavation the prediction of rock mass conditions ahead of the tunnel face can be improved. Implementing the rules developed from experience and simulations into the monitoring data evaluation program allows to automatically issuing information on the expected rock mass quality ahead of the tunnel.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.9
• /
• pp.973-978
• /
• 2014

This study focused on the design of the reflecting layer of a light trapping system fora thin film solar cell using topology optimization based on the phase field method. Therefore, incident light was caused to propagate in the desired direction by reflecting it from this layer, which is the design domain. The same method was applied to the conceptual design of an infrared stealth structure in near infrared range. The results using the phase field method were compared with those using the density method. The design objective was to maximize the Poynting vector value representing the energy flux, which was measured in a measuring domain to control the reflected wave direction. A finite element analysis and optimization process were performed using the commercial package COMSOL combined with the MATLAB programming.

• Geophysics and Geophysical Exploration
• /
• v.7 no.2
• /
• pp.148-154
• /
• 2004

This article reviews the development of three-dimensional (3-D) magnetotelluric (MT) modeling. The 3-D modeling of electromagnetic fields is essential in understanding the physics of MT soundings, and in implementing an inversion method to reconstruct a 3-D resistivity image. Although various numerical schemes have been developed over the last two decades, practical methods have been quite limited. However, the recent rapid improvement in computer speed and memory, as well as the advance in iterative solution algorithms for a large system of equations, makes it possible to model the MT responses of complex 3-D structures, which have been very difficult to simulate before. The use of staggered grids in finite difference method has become popular, conserving a magnetic flux and an electric current and allowing for realistic discontinuous fields. The convergence of numerical solutions has been greatly accelerated by adopting Krylov subspace methods, proper preconditioning techniques, and static divergence corrections. The vector finite-element method using edge elements is also free from the discontinuity problem, and seems a natural choice for modeling complex structures including irregular topography because its flexibility allows one to capture full geometric complexity.

• Journal of the Korean GEO-environmental Society
• /
• v.4 no.3
• /
• pp.79-88
• /
• 2003

This paper presents a numerical study of the spatial behavior of a linear absorbing solute in a heterogeneous porous medium. The spatially correlated log-normal hydraulic conductivity field is generated in a given two-dimensional domain by using the geostatistical method (Turning Bands algorithm). The velocity vector field is calculated by applying the two-dimensional saturated groundwater flow equation to the Galerkin finite element method. The simulation of solute transport is carried out by using the random walk particle tracking model with CD(constant displacement) scheme in which the time interval is automatically adjusted. In this study, the spatial behavior of a solute is analyzed by the longitudinal center-of-mass displacement, longitudinal spatial spread moment and longitudinal plume skewness.

• Journal of the Korean Society for Railway
• /
• v.12 no.3
• /
• pp.437-442
• /
• 2009

In the case of a Linear Induction Motor (LIM), numerical analysis method like Finite Element Method (FEM) has been mainly used to analyze the travelling magnetic field problem which includes the velocity-induced electromotive force. If the problem including the velocity-induced electromotive force is analyzed by FEM using the Galerkin method, the solution can be oscillated according to the Peclet Number, which is determined by conductivity, permeability, moving velocity and size of mesh. Consequently, the accuracy of the solution can be low and the vortexes of flux can be occurred at the secondary back-iron. These vortexes of the flux occurred at the secondary back-iron does not exist physically, but it can be occurred in the analysis. In this case, the vortexes of the flux can be generally removed by using Up-Wind method which is impossible to apply a conventional S/W tool (Maxwell 2D). Therefore, in this paper, authors examined the vortexes of the flux occurred at the secondary back-iron of the LIM according to variations of the Peclet Number, and analyzed whether these vortexes of the flux affect on the dynamic force characteristics of the LIM or not.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.5
• /
• pp.345-352
• /
• 2014

In this paper, the CAD data for the optimal shape design obtained by isogeometric shape optimization is directly used to fabricate the specimen by using 3D printer for the experimental validation. In a conventional finite element method, the geometric approximation inherent in the mesh leads to the accuracy issue in response analysis and design sensitivity analysis. Furthermore, in the finite element based shape optimization, subsequent communication with CAD description is required in the design optimization process, which results in the loss of optimal design information during the communication. Isogeometric analysis method employs the same NURBS basis functions and control points used in CAD systems, which enables to use exact geometrical properties like normal vector and curvature information in the response analysis and design sensitivity analysis procedure. Also, it vastly simplify the design modification of complex geometries without communicating with the CAD description of geometry during design optimization process. Therefore, the information of optimal design and material volume is exactly reflected to fabricate the specimen for experimental validation. Through the design optimization examples of elasticity problem, it is experimentally shown that the optimal design has higher stiffness than the initial design. Also, the experimental results match very well with the numerical results. Using a non-contact optical 3D deformation measuring system for strain distribution, it is shown that the stress concentration is significantly alleviated in the optimal design compared with the initial design.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.22 no.1
• /
• pp.47-57
• /
• 2010

Subgrid scale stabilization method is applied to Woo and Liu(2004)'s extended Boussinesq FEM numerical model to eliminate the 2dx wiggles. In order to optimize the computational efficiency, Hessian operator is introduced and the matrix of velocity vector is combined to one matrix for solving matrix equations. The mass lumping technique is also applied to the matrix equations of auxiliary variables. The newly developed code is applied to simulate Vincent and Briggs(1989)' wave transformation experiments and the results show that the numerical solution is almost wiggle-free and it matches very well with experimental data. Due to improvement of computational efficiency and wiggle reduction, it is plausible to apply this model to a realistic problem such as harbor oscillation problems.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.155-162
• /
• 2014

Using an isogeometric approach, a continuum-based shape design optimization method is developed for steady state power flow problems at high frequencies. In case the isogeometric method is employed to the shape design optimization, the NURBS basis functions used in CAD geometric modeling are directly utilized to embed the exact geometry into the computational framework so that the design parameterization for shape optimization is much easier than that in the finite element method and consequently provides the enhanced smoothness of design perturbations. Thus, exact geometric models can be used in both the response and the shape sensitivity analyses, where normal vector and curvature are continuous over the whole design space so that enhanced shape sensitivity can be expected. Through numerical examples, the developed isogeometric sensitivity is compared with finite difference one to provide excellent agreement. Also, it turns out that the proposed method works very well in the shape optimization problems.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.63 no.10
• /
• pp.1377-1383
• /
• 2014

This paper deals with eddy current loss of magnetic coupling with radial permanent magnet (PM) using analytical method such as a space harmonic method. Superposition of two kinds analysis model is used to analyze eddy current loss induced in inner PM and outer PM of magnetic coupling. When the eddy current is induced, the environmental temperature increases, and the permanent magnet(PM) characteristics are degraded because the performance of PM is greatly influenced by temperature rise. Hence, the calculation of eddy current loss becomes an important factor in the magnetic coupling. In order to analyze eddy current loss, first, on the basis of the magnetic vector potential and two-dimensional(2-D) polar-coordinate system, the magnetic field solutions of the radial magnetized PM are obtained. And we obtain the analytical solutions for the eddy current density produced by permanent magnet. Lastly, analytical solutions for eddy current loss are derived by using equivalent, electrical resistance calculated from magnet volume and analytical solution for eddy current density. This analytical results are validated by comparing with the 2-D finite element analysis (FEA).

• The Transactions of the Korean Institute of Electrical Engineers P
• /
• v.55 no.3
• /
• pp.141-145
• /
• 2006

The line current problem(2-dimensional space : point source) is not easy to analyze the magnetic field using the standard finite element method(FEM), such as overhead trolley line or transmission line. To supplement such a defect this paper is proposed the coupling scheme of analytical solution and FEM. In analysis of the magnetic field using the standard FEM. If the current region is a relatively small compared to the whole region. Therefore the current region must be finely divided using a large number of elements. And the large number of elements increase the number of unknown variables and the use of computer memories. In this paper, an analytical solution is suggested to supplement this weak points. When source is line current and the part of interest is far from line current, the analytical solution can be coupling with FEM at the boundary. Analytical solution can be described by the multiplication of two functions. One is power function of radius, the other is a trigonometric function of angle in the cylindrical coordinate system. There are integral constants of two types which can be established by fourier series expansion. Also fourier series is represented as the factor to apply the continuity of the magnetic vector potential and magnetic field intensity with tangential component at the boundary. To verify the proposed algorithm, we chose simplified model existing magnetic material in FE region. The results are compared with standard FE solution. And it is good agreed by increasing harmonic order.